Next Article in Journal
Heat Source Models in Numerical Simulations of Laser Welding
Previous Article in Journal
The Photocatalytic Activity of Titania Coatings Produced by Electrochemical and Chemical Oxidation of Ti6Al4V Substrate, Estimated According to ISO 10678:2010
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Teicoplanin-Modified HPLC Column as a Source of Experimental Parameters for Prediction of the Anticonvulsant Activity of 1,2,4-Triazole-3-Thiones by the Regression Models

1
Department of Analytical Chemistry, Medical University of Lublin, Chodźki 4a, 20-093 Lublin, Poland
2
Faculty of Medicine, Medical University of Lublin, Aleje Racławickie 1, 20-059 Lublin, Poland
3
Department of Biochemical Diagnostics, Medical University of Lublin, Staszica 16, 20-081 Lublin, Poland
4
Chair and Department of Biophysics, Medical University of Lublin, Jaczewskiego 4, 20-090 Lublin, Poland
5
Department of Applied Mathematics, University of Technology, Nadbystrzycka 38D, 20-618 Lublin, Poland
6
Department of Pharmacology, Faculty of Nursing and Health Sciences, Medical University of Lublin, Chodźki 4a, 20-093 Lublin, Poland
7
Jerzy Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, Niezapominajek 8, 30-239 Krakow, Poland
*
Author to whom correspondence should be addressed.
Materials 2020, 13(11), 2650; https://doi.org/10.3390/ma13112650
Submission received: 2 May 2020 / Revised: 27 May 2020 / Accepted: 29 May 2020 / Published: 10 June 2020
(This article belongs to the Section Materials Simulation and Design)

Abstract

:
The cell membrane is a complex system that consists of lipids, proteins, polysaccharides, and amphiphilic phospholipids. It plays an important role in ADME processes that are responsible for the final pharmaceutical effects of xenobiotics (bioavailability, activity). To study drug-membrane interaction at the molecular level, several high-performance liquid chromatography (HPLC) membrane model systems have been proposed which are mimicking mainly its lipid character. The aim of this work was to study interactions of new synthesized antiepileptic compounds of 4-alkyl-5-(3-chlorophenyl)-2,4-dihydro-3H-1,2,4-triazole-3-thione derivatives with Chirobiotic column containing glycoprotein ligand attached to the silica matrix. The affinity of the analytes to immobilized glycoprotein ligand was examined chromatographically in reversed-phase mode. The thermodynamics of interactions between bioactive compounds and teicoplanin was studied in terms of the van’t Hoff linear relationship ln k vs. 1/T in the range of 5–45 °C. Change in enthalpy (ΔH°), change in entropy (ΔS°) and change in Gibbs free energy (ΔG°) were estimated utilizing graphical extrapolation and interpolation methods. The density functional theory (DFT) approach and docking simulations were used to get the molecular interpretation and prove the obtained experimental results. Cross-correlations of chromatographic and thermodynamic parameters with non-empirical topological and quantum chemical indices suggest that the polarizability of analytes appears to be responsible for the interactions of the tested molecules with teicoplanin and, ultimately, their retention on the column. Experimental and theoretical parameters were subjected to statistical analysis using regression models. Partial least squares (PLS) regression model showed the usefulness of the experimentally measured parameter φ0 (MeOH) to discriminate between anticonvulsant active and inactive 1,2,4-triazole-3-thione derivatives. Obtained results point out the usefulness of interaction of potential anticonvulsants with glycoprotein class of compounds to anticipate their activity.

Graphical Abstract

1. Introduction

In a living organism, xenobiotics undergo many biochemical and biophysical processes that are abbreviated as ADME, i.e., adsorption, distribution, metabolism, and elimination. Each of these stages depends in different ways on the lipophilicity of the xenobiotic. Therefore, many empirical and non-empirical lipophilicity scales are commonly used to predict the pharmacological activity of potential drugs. The process of penetration of exogenous substances through semipermeable membranes can be reflected by a partition coefficient. The oldest extra-thermodynamic extraction system, n-octanol-water, is considered to be the first model to simulate biological permeability. This two-phase system is a source of a parameter reflecting the relative lipophilic character of the compound, which is defined by log P, where P is the octanol-water partition coefficient [1,2]. However, in fact, the cell membrane is a complex system that consists of lipids, proteins and polysaccharides, and amphiphilic phospholipids. The drug’s ability to cross the membrane barrier is responsible for the observed range of pharmaceutical effects (bioavailability, activity). To study drug-membrane interaction at the molecular level, several membrane models have been created that mimic natural cell membranes, such as stationary phases in chromatography. Owing to this, it was possible to perform experiments under strictly defined and controlled conditions. Examples of surface-confined membranes are lipid vesicles, liposomes, Langmuir monolayers, and solid-supported lipid monolayers/bilayers (Scheme 1).
In the beginning, lipophilicity tests were performed by HPLC on octadecyl-dimethyl silane (ODS) phases with silica as solid support [3,4]. Chromatographic retention parameters (k, tr, logk, logkw) have been recognized as independent parameters of lipophilicity and are still successfully used to model biological lipophilicity [5,6,7,8]. In the 1990s, in order to model drug transport, solid phases known as artificial membranes (IAM) were synthesized [9,10,11,12,13]. Monolayers of phospholipid analogs were covalently bound to the silica surface or immobilized on conventional ODS columns. IAM-type stationary phases containing a phosphatidylcholine ligand or lipid-like, such as a cholesteric moiety covalently bonded to silica through a propylamine chain, are currently commercially available [14,15]. The retention mechanism on IAM columns involves significant intermolecular forces such as ionic bonds, electrostatic interactions between drug molecules, and the membrane, which is significant progress over historical octanol/water partition systems, or RP-HPLC exposing mainly hydrophobic interactions. That is why the relationships between the retention parameters (log k) obtained in the IAM column and log Po/w coefficients were rather weakly correlated [11,14,16].
Recently, immobilized liposome chromatography (ILC) raises great interest and hopes in simulating membrane-drug interactions. Many researchers emphasize its usefulness for the screening of drug partition in lipid bilayers [17,18,19,20,21]. The potential of this method is associated with the possibility of introducing additional ligands to the liposome surface. However, the widespread use of the ILC method limits the poor stability of liposomes and the need for self-preparation of the column.
To mimic the process of transport of xenobiotics across biological membranes, many solid carriers imitating biological membranes have been used to predict the interaction between the membrane and the drug, such as human skin, intestinal epithelium, and microvascular vascular endothelium [22,23,24].
The current work concerns the bioactivity prediction of newly synthesized substances with activity in the central nervous system (CNS). The therapeutic effect in CNS diseases is limited by a poor response to drug treatment. The main causes of therapeutic failure are mainly changes in drug uptake into target brain cells as a result of limited blood-brain barrier (BBB) ​​permeability, which impairs drug penetration through the brain [25,26,27]. Despite the fact that ordinary antiepileptic drugs (phenytoin, carbamazepine, oxcarbazepine, lamotrigine, phenobarbital, felbamate, valproic acid, topiramate) penetrate well into the brain, they have low efficiency. As a result, the therapy appears to be ineffective in approximately 30–40% of patients with epilepsy [25,28]. Therefore, reliable tests are required to predict whether newly synthesized compounds potentially exhibit anti-epileptic activity or only the ability to cross the blood-brain barrier (BBB). The most common experimental measures of BBB permeability are log BB (drug concentration in the brain divided by blood concentration) and log PS (product of the permeability surface). Both methods are time-consuming and expensive [29]. Other valuable alternatives to predicting BBB permeability can be provided by in silico calculation methods.
In this work, a Chirobiotic T column and a water-organic mobile phase operating in reversed-phase high-performance liquid chromatography (RP-HPLC) were used as a model to study a series of new candidates for anticonvulsants. The proposed RPLC is a dynamic two-phase system capable of mimicking biological membranes. Until now, various sorbents and artificial membranes have been used to simulate biological processes occurring at the interface. They partly reflected the structure of biological barriers, such as cholesterol in phosphatidylcholine columns, or focused only on demonstrating lipophilic properties, as in the case of classic ODS columns. The Chirobiotic T column was chosen for its structure representing the glycoprotein class of chemical compounds, anisotropic properties, reverse phase capability, and both types of chemical interaction, i.e., polar and hydrophobic, which makes it similar to complex and heterogeneous biological membranes.
Chirobiotic columns are obtained by attaching high purity macrocyclic glycoproteins to a spherical silica gel. The Chirobiotic T column has an attached amphoteric glycopeptide-eicoplanin. It has 23 chiral centers surrounding four cavities. Sites that can be both hydrogen bond donors and acceptors are located near seven aromatic rings. In 1994, Armstrong et al. described the first use of vancomycin, thiostreptone and rifamycin B macrocyclic antibiotic as chiral stationary phases (CSP) [30,31]. The chemical structure of these stationary phases enables various interactions with separated analytes, such as hydrophobic, π-π, dipole-dipole, hydrogen, electrostatic, ionic, and van der Waals interactions [32], as well as complex formation [33]. In turn, chiral recognition is possible due to a large number of stereogenic centers [34,35]. Teicoplanin contains amide linkages and sugar residues. It possesses an element of tertiary structure similar to proteins. In comparison to very large molecules like proteins, the teicoplanin-based stationary phase shows the dynamic mass transfer allowing separations for a wide range of compound classes.
The titled compounds belong to quite a new group of voltage-gated sodium channel blockers, the 4-alkyl-5-substituted-1,2,4-triazole-3-thione derivatives with a wide log P range (1.96–6.46). Docking simulations confirmed that these compounds bind to voltage-sensor domain IV (VSD4) located in mammalian NaV 1.7 proteins [36]. Importantly, the 1,24-triazole-3-thione derivatives turned out to possess strong anticonvulsant activity in various animal models of epilepsy, including MES and 6 Hz tests [37,38,39,40,41]. Simultaneously, they demonstrated many beneficial ADME-Tox properties, i.e., low neurotoxicity, lack of genotoxic effects observed in HepG2 cells, and the lack of negative impact on the viability of human cells when administered in concentrations necessary to achieve antiepileptic effects. Both pharmacological and physicochemical properties of the 4-alkyl-5-substituted-1,2,4-triazole-3-thione derivatives were highly affected by the structure of the alkyl substituents. The compounds containing unbranched chains produced stronger anticonvulsant effects than respective branched analogs. Their affinity towards respective molecular targets (i.e., VGSCc) was also stronger. However, the ability of the 1,2,4-triazole-based molecules to pass through the biological barriers (e.g., intestinal and blood-brain barriers) was much higher when the unbranched alkyl chain was replaced by the branched one [36,37,38,39,40,41]. Such a characteristic is of crucial importance in the development of new antiepileptics since the treatment of epilepsy is based mainly on oral medication. Therefore, the investigated 1,2,4-triazole derivatives (named TPR) seem to be an interesting group of antiepileptic drug candidates. In our previous papers [15,42], we examined these compounds on the cholesterol-based stationary phase, and the commercially available phosphatidylcholine biomimetic stationary phase [42]. Molecular dynamics simulations revealed the selective preferences of the phosphatidylcholine with respect to the compounds with either ether or sulfide moieties.
The aim of the current work concerns the retention mechanism of 1,2,4-triazole-3-thione derivatives on teicoplanin used as the stationary phase modifier. The thermodynamic van’t Hoff relationships and docking simulation were used for the above purpose. Retention behavior of studied molecules in this new dynamic two-phase system revealed some new factors influencing the observed chromatographic behavior as well as pharmacological activity. Furthermore, presented experimental (HPLC, DCS) and theoretical (docking simulation, DFT calculations, classification analysis) results are helpful for an understanding of the relationship between the structure vs. the properties of studied 1,2,4-triazole-3-thiones, and the function of teicoplanin as a new membrane mimicking the biological barrier.

2. Methods and Materials

2.1. Solutes and Reagents

The synthesis procedure of the derivatives of the 4-alkyl-5-substituted-1,2,4-triazole-3-thione derivatives (Table 1) was described previously [38]. Their structures were confirmed by gas chromatography-mass spectrometry (GC-MS), proton nuclear magnetic resonance (1H-NMR), carbon-13 nuclear magnetic resonance (13C-NMR), infrared spectroscopy (IR), and their purity was assessed by elemental analysis. The teicoplanin reference standard for DSC measurements was purchased from Sigma-Aldrich, Poznan, Poland). HPLC grade methanol (MeOH), and acetonitrile (ACN) were obtained from Merck (Darmstadt, Germany). Water purified by ULTRAPURE Millipore Direct-Q 3UV-R (Merck, Darmstadt, Germany) of the resistivity 18.2 MΩ cm was used to prepare all the aqueous solutions.

2.2. HPLC Measurement

The retention factors were measured with a liquid chromatograph, an Elite LaChrom HPLC Merck-Hitachi (Merck, Darmstadt, Germany) with DAD (Diode Array Detector L-2455) detector, pomp L-2130 and a manual sample injection valve equipped with a 20 µL loop and EZChrom Elite software (Merck) system manager as a data processor. The column was an Astec CHIROBIOTICTM (SUPELCO, Bellefonte, PA, USA); (150 mm × 4.6 mm I.D., 5 µm), pore size: 100Å. The injection of blank mobile phase volumes produced visible detector fluctuations that were used as the hold-up volume. The column was thermostated at 20 °C ± 0.1 using column thermostat Jetstream 2 Plus (100375, Knauer, Sigma-Aldrich, Poznan, Poland). The elution was carried out in the isocratic mode by the mobile phase consisting of methanol or acetonitrile in water. The mobile phase was filtered through a Nylon 66 membrane filter (0.45 μm) Whatman (Maidstone, England) by the use of a filtration apparatus.
The retention data were recorded at a flow-rate of 1 mL min−1 with online degassing using L-7612 solvent degasser at a wavelength chosen accordingly with the recorded spectra. The solutions (10−4 M to 10−3 M) were prepared by dissolving the solutes in the mobile phase. Typical injection volumes were 3 µL.

2.3. Differential Scanning Calorimetry (DSC) Experiments

Differential scanning calorimetry (DSC Q200 calorimeter, TA Instruments, Mettler-Toledo, Warsaw, Poland) was applied to study thermal transformations of the dry stationary phase. 5–8 mg samples were closed in aluminum pans (Tzero Pan) and temperature scanned from 0–210 °C at heating 5 °C/min in a nitrogen atmosphere. An empty pan was used as a reference.

2.4. The Methodology of Statistical Analysis

All calculations and visualizations were conducted in R environment [43] with a few packages which extend the possibilities of the base program, like caret [44], DALEX [45], ggpubr [46], random Forest [47], glmnet [48,49], and pls [50]. All linear regression analyses were performed using Excel (version from Office 97).
Two types of models were used to predict activity: classification and regression. For classification purposes, Regularized Logistic Regression, Random Forest, and Partial Least Squares Discriminant Analysis were used. In the regression approach, Random Forest, Partial Least Squares Regression, and Penalized Regression were applied to predict the dependent variable. The intuition standing behind all applied methods is to reduce the dimensionality and the collinearity in the pre-processing stage of analysis. To avoid estimation problems resulting from the high dimensionality of the predictive space and the lack of a sufficient number of observations, two methods of feature selection were used: Recursive Feature Elimination [51] and Simulated Annealing Feature Selection [52]. The recursive process and regularization were used to manage collinearity, while the high dimensionality was reduced using the partial least squares technique.
Data were randomly divided into train and test samples (they have 9 and 5 observations, respectively). The obtained models were verified on the test sample. During the model tuning procedure, a repeated 3-fold cross-validation method was used [53,54]. Models fit were assessed by Root Mean Square Error (RMSE), R2, and Mean Absolute Error (MAE), but the best fit was chosen based on RMSE, considering ED50 to predict response, while classification model’s accuracy was used to assess the model performance.

2.5. The Methodology of Molecular Modeling

The ligand molecules were prepared by using the Avogadro 1.1.1 software [55] and optimized within the UFF force field [56] (steepest descent algorithm). The flexible, optimized ligands were docked to the teicoplanin molecule according to the 1:1 binding stoichiometry. The teicoplanin structures found in the four following PDB entries were applied: 4mfl, 4pk0, 5awv, and 4pzj (X-ray resolutions varying from 0.16 to 0.23 nm). The AutoDock Vina software [57] was applied for docking simulations. The procedure of docking was carried out within the cuboid region covering the whole teicoplanin molecule. All the default procedures and algorithms implemented in AutoDock Vina were applied during docking simulations. The predicted binding energies were averaged over all the four teicoplanin structures. The more favorable binding mode is associated with the lower binding energy value; only the lowest energy value corresponding to the given ligand/teicoplanin system was accounted for during subsequent analyses. In accordance with the validated methodology of Vina-based docking, no explicit solvent was considered. However, the solvent effects are partially accounted for in an implicit manner by using the potentials calibrated against the data measured for solvent-containing systems.
The docking procedure was followed by the density functional theory (DFT)-based calculations. The geometry optimization was performed with the B3LYP nonlocal exchange correlation functional [58] and the 6-311 + G basis set [59]. The presence of the implicit solvent (methanol) was accounted for by applying the COSMO approach. The possibility of introducing the explicit solvent is limited due to both computational efficiency reasons and the large configurational complexity of the solvent-containing systems. The default convergence criteria implemented in Gaussian09 [60] were accepted. Calculations involved the teicoplanin molecule, molecules of all ligands, and the ligand-teicoplanin complexes. The initial structures were taken from the docking results; wherever necessary, the part of the system was removed to obtain the derided composition. The binding energies were calculated from the energy difference between (i) energy of the complex; (ii) the sum of the energies of isolated teicoplanin and ligand molecules.

2.6. Evaluation of Anticonvulsant Activity

The details of the anticonvulsant activity evaluation procedure and subsequent results are partially described in our previous works [15,37,42]. The anticonvulsant activity of the investigated compounds was measured in mice maximal electroshock (MES)-induced seizure model of epilepsy. The MES model is considered to mimic generalized tonic-clonic seizures in humans. The experiments were conducted on adult Swiss albino mice (males) weighing 20–25 g, which were kept in colony cages with free access to food and tap water. After one week of acclimatization, the animals were randomly divided into groups consisting of 8 mice. The investigated compounds were suspended in a 1% solution of Tween 80 in distilled water and administered intraperitoneally in a volume of 5 mL/kg. Electroconvulsions were produced by an alternating current (0.2 s stimulus duration; 500 V, 50 Hz, fixed current intensity of 25 mA) delivered via ear-clip electrodes by a Type 221 Rodent Shocker generator (Hugo Sachs Elektronik, Freiburg, Germany). The abolition of a hind-limb tonic extension (HLTE) was taken as a criterion of the anticonvulsant activity. The tested compounds were administered in increasing doses in order to calculate ED50 values using the log-probit method by Litchfield and Wilcoxon. Each ED50 reflects the dose of a compound required to protect 50% of the animals tested against maximal electroshock (MES)-induced seizures.

3. Results and Discussion

Two classical models usually used to interpret the dependence of the retention factor (k) on the mobile phase composition and the temperature are described by the Soczewiński-Wachtmeister equation (Section 3.1) known as the linear solvent strength model (LSSM), and the second one concerning thermodynamic studies represented by the Van’t Hoff plot (Section 3.2).

3.1. Retention Behavior on Teicoplanin Column

The main purpose of the initial experiments was to examine the retention behavior of new 1,2,4-triazole-3-thion derivatives on the teicoplanin column. Teicoplanin stationary phase may be used under reversed-phase (RP) conditions with low proportions of polar organic solvent to an aqueous component. Retention and selectivity on the teicoplanin column depend primarily on the concentration and nature of the organic modifier (methanol or acetonitrile), similarly to the most reversed-phase systems. Higher concentrations of the organic modifier resulted in lower retention, typically for reversed-phase (RP) systems (Figure 1).
The retention behavior of the analyzed substances was tested in polycratic conditions. Methanol concentration was ranging from 0.15 to 0.40 of the volume fraction of the organic modifier in the eluent (φ), whereas acetonitrile concentration was from 0.05 to 0.3 φ.
Obtained linear dependences of log k on φ, are described by the Soczewiński-Wachtmeister Equation [61]:
log k = − Sφ + log kw
Based on these relationships, log k values extrapolated to 100% water (log kw) and the S values expressing the slope of the following equation were determined. The chromatographic lipophilicity parameters, for two selected eluent systems (A and B), are presented in Table 2. The coefficients of determination (R2) of the above correlations were higher than 0.9 for the acetonitrile-water mobile phase (A) as well as in the case of methanol as an eluent modifier (B). The high values of F and small standard correlation errors (se < 0.05) also testify to the perfect matching of the dependencies to the straight line.
The obtained retention parameters appear to favor the acetonitrile-type mobile phase. The average slope values of the relationship log k vs. φ are 25% higher for acetonitrile- (S = −3.816) containing eluent system in comparison to methanol (S = −2.846).
The log kw parameter is the theoretical retention of the analyte in pure water or buffer. It is very often used in quantitative structure–activity relationships (QSAR) as a lipophilicity parameter. Theoretically, this parameter eliminates the influence of the organic modifier, however, it should be remembered that it refers rather to interphase, which is the stationary phase changed by the presence of mobile phase molecules.
In 1993, Valko and Sleger [62] introduced another parameter that is used to compare lipophilicity, i.e., φ0. This parameter illustrates the composition of the eluent for which k is 1, i.e., the amount of substance in the mobile and stationary phases equals each other. The φ values are determined by the interpolation method of the relationship log k = f(φ). It turns out that this parameter is less dependent on the type of organic modifier used in the eluent, and it is better reproducible and repeatable compared to log kw [63,64].
The log kw vs. S correlation graphs are used to study the structural similarity of analytes. For related substances that do not differ in donor-acceptor properties, statistically, significant relationship S vs. log kw should be observed. For both elution systems (methanol-water, acetonitrile-water) satisfactory correlations described by the following linear regression equations were obtained:
MeOH/Water: log kw = −0.4205(±0.03939)S + 0.4540(±0.1139)
n = 14, R2 = 0.8976, F = 113.9818, se = 0.062873
ACN/Water: log kw = −0.3359(0.009338)S + 0.3145(0.36408)
n = 14, R2 = 0.9901, F = 1294.178, se = 0.019135
where numbers in parentheses correspond to 95% confidence limits; n is the number of the compounds; R2 is the squared correlation coefficient; se is the standard deviation and F is the Fisher’s test.
The obtained results prove the usefulness of the Chirobiotic T column for the assessment of structural similarities of the investigated 1,2,4-triazole-3-thiones. Furthermore, significant correlations, presented above, additionally imply that the retention of the investigated analytes is controlled by the same intermolecular forces.
Regardless of the type of organic modifier in the mobile phase, the log kw values and the S parameters are very similar to each other. These similarities are reflected by the relationship of log kw (MeOH/Water) vs. log kw (ACN/Water), and S (MeOH/Water) vs. S (ACN/Water), (Figure 2A,C) which has a slope close to one, a low value of the intercept, close to one correlation coefficient, and a small standard error. Only the φ0 values are not correlated with each other in the tested elution systems.
log kw (MeOH/Water) = 0.9952(±0.06705) log kw (ACN/Water) + 0.0541(±0.1087)
n = 14, R2=0.9443; se = 0.04638; F = 220.3161
S (MeOH/Water) = 0.7682(±0.0777) S (ACN/Water) + 0.08605(±0.2988)
n = 14, R2 = 0.8907; se = 0.1384; F = 97.8039
φ0 (MeOH/Water) = 1.1149(±0.6355) φ0 (ACN/Water) − 0.1033(±0.2716)
n = 14, R2 = 0.2041; se = 0.02442; F = 3.0779
This result is related to the linearity of the Soczewiński-Wachtmeister relationship. Linear relationships are obtained in limited regions of the organic modifier concentration. Usually, deviations from linearity occur at higher modifier concentrations. This effect is due to excessive adsorption of the organic solvent in the stationary phase. The concept of the adsorbed layer was first suggested by Knox and Pryde [65] and confirmed by Kazakevich et al. [66]. As with other RP systems, the retention process on the teicoplanin column with elution systems containing methanol/acetonitrile can also be modeled using the Soczewiński-Wachtmeister relationship. However, the obtained relationships were better correlated in the case of methanol playing the role of solvent (with a correlation coefficient approaching 0.99 in the entire range of tested concentrations). Thus, in this case, the results of extrapolation and interpolation cannot be questioned.
In turn, correlations between log k and acetonitrile content were slightly worse with an average R2 value approaching 0.98, suggesting that acetonitrile may form a multilayer which is adsorbed on top of the bound phase. In this case, the Soczewiński-Wachtmeister relationship may not be true in a wide range of organic content. Unlike methanol, excessive adsorption of acetonitrile at the surface can cause a reduction in analyte-surface interaction and, consequently, the corresponding retention. This is the main reason for the disturbances from the linear shape of the Soczewiński-Wachtmeister graphs appearing in the form of a concave shape. Since deviations from linearity usually occur at higher acetonitrile concentrations, the most disturbing is the interpolation results obtained in acetonitrile (φ0).
The lipophilicity parameters can also be determined through quantum-chemical calculations. The theoretical method comprised calculations of AC log P using the in silico model, based on Hansch and Leo’s approach [67], from the VCCLAB online calculator [68].
As can be seen in Figure 3, the theoretical hydrophobicity values expressed by AC log P are much higher than the chromatographic scale of lipophilicity represented by log kw values, either in the methanol or acetonitrile eluent. Thus, in a two-phase system in an aqueous-organic environment, the tested compounds showed more uniform hydrophobicity, being at least two times lower in comparison to the calculated values. This may indicate the induction of the dipole moment of the studied molecules by the polar environment, which results in the weakening of hydrophobic interactions.

3.2. Thermodynamic Studies

Temperature effects also provide very useful information to explain retention mechanisms and better understand analyte behavior in the chromatographic systems. Moreover, elevated temperatures offer many advantages in liquid chromatography. An increase in temperature can improve the sample throughput because van Deemter minima are transferred to higher flow rates. Besides, as the temperature rises, the water behaves more like an organic solvent, therefore, it is possible to reduce the amount of organic modifier at elevated temperatures. The thermodynamic parameters were obtained based on the temperature dependence of the retention factor (k) known as the van‘t Hoff relationship:
ln k = −ΔH°/RT + ΔS°/R + ln Φ
where Φ is the phase ratio; (Vs/VM) is the ratio of the stationary phase volume to the mobile phase volume, the standard enthalpy change (ΔH°), and the standard entropy change (ΔS°); R is the universal gas constant (8.31 J mol−1 K−1); T is the absolute temperature and k is the retention factor in chromatography.
The individual van‘t Hoff plots representing log k vs. 1/T were generated under isocratic conditions involving the following elution systems: 15% ACN and 15% MeOH (Figure 4A,B). The obtained graphs were almost linear with the sufficient values of the determination coefficient (R2 > 0.9) only for methanol-containing eluent systems. Sufficient linearity of the van’t Hoff graphs proves that the retention mechanism is constant over the temperature range under consideration, and makes them useful in measuring thermodynamic parameters such as the ΔH°, ΔS°, and free energy of transfer (ΔG°) of 1,2,4-triazoles from mobile phases to the investigated sorbent. In the further static analysis, only parameters obtained in the methanol-water system were taken into account as the most reliable. From ln k vs. 1/T plots, the slope and intercept are −ΔH°/R and ΔS°/R + ln Φ, respectively. The Gibbs free energy change (ΔG°) was calculated by the Gibbs-Helmholtz equation:
ΔG° = ΔH° − TΔS°
In contrast to methanol, acetonitrile-containing mobile phases formed non-linear van’t Hoff plots. According to Chester and Coym’s study [69], the disturbances in a van’t Hoff plot are the consequences of a phase transition in the stationary phase or changes in the phase ratio. In this case, the thermodynamic parameters derived from the van’t Hoff plots cannot be trusted. Determination of the phase ration (Φ) in chromatography is the subject of debate. There are some methods of its determination proposed among others by Dorsey [70,71]. However, most thermodynamic studies assume that the phase ratio values are constant. Although, there is some evidence that the volume of the stationary phase, which is accessible to analyzed molecules, changes with the mobile phase composition and temperature. There is evidence, especially for the ODS bonded column, of structure collapse changing its conformation, in the case of high water content in the eluent system [66] or temperature increase [72]. Moreover, there is evidence that the phase ratio is also solute-dependent [69]. As interphase changes its conformation and finally, its accessible volume with temperature, it will further affect interactions with different analytes by changing steric conditions and efficiency of solute diffusion into the stationary phase. We also cannot ignore the accessibility of the mobile phase to solutes, because some of them are permanently excluded from the liquid phase owing to unfavorable charges.
The impact of various factors, i.e., phase transition, type of analyte, and especially their charges, temperature change and, even rarely considered in HPLC, pressure changes that affect the volume of sorbent available for interaction through the compression effect, on thermodynamic parameters including the enthalpy value, is very complex and not completely clear. In addition, according to Chester and Coym [69], these factors often offset each other. Efforts to obtain reliable values for the corresponding parameters are only relevant when we compare thermodynamic parameters measured on different columns. As in our work only one stationary phase was used, no phase transition for teicoplanin was observed in the studied temperature range (Section 3.3), so then the thermodynamic parameters were calculated directly from the van’t Hoff curves if they were sufficiently linear. For this reason, the further static analysis included only the parameters obtained in the methanol-water system as the most reliable.
All parameters calculated based on the simple linear regression analysis have been collected in Table 3. The obtained ΔS° values (ΔS°*) are biased by the constant error, namely R ln Φ which could be completely canceled for ΔΔS calculations.
The obtained values of the enthalpy change were all negative, but the more negative ΔH° values were obtained for methanol-containing eluent systems in comparison to acetonitrile ones, indicating the existence of more favorable, exothermic interactions between the 1,2,4-triazole derivatives and the teicoplanin stationary phase. The changes in enthalpy are 1000 times higher as compared to changes in entropy. Therefore, the free energy is favorable (ΔG° < 0), mainly due to the enthalpic contribution for both eluents.
The ΔH° values for a series of the compounds containing unbranched chains TPR 14, 16, 28, 30 became the most negative (exothermic). It was proved previously that these compounds possess stronger anticonvulsant effects, the highest affinity towards respective molecular targets, and the smaller ability to pass through the biological barriers than respective branched analogs [37]. The obtained chromatographic results show that the incorporation of an unbranched, long hydrocarbon chain into the entire structure significantly enhanced interactions with the model membrane. The ΔH° values for increasingly hydrophobic analytes containing more carbon atoms should become more negative leading to larger values of k. This effect is shown in Figure 1. The more hydrophobic analytes containing 7 carbon atoms (TPR 34) in substituent have larger retention than those with shorter chain lengths containing 4 carbon atoms (TPR 2, TPR 20). To study the interactions’ pattern between examined compounds and teicoplanin, docking simulations have been performed.

3.3. The Phase Transition of Teicoplanin

To examine a phase transition, differential scanning calorimetry (DSC) was applied. Obtained results showed, over a temperature range of 30–150 °C, that the melting temperature of neat teicoplanin equals 124.47 °C (Figure 5). Dry teicoplanin did not undergo any phase transition under studied conditions in the range from 5 to 40 °C. As it is above the examined range, it indicates that rather the phase ratio or structure of the individual solutes accessible for the different stationary phase volume cause observed non-linearities in the van’t Hoff plots.

3.4. Selectivity of Teicpolanin Column towards Studied Triazoles

Chester and Coym elaborated a procedure [69] useful for estimation of transfer enthalpy which is completely independent of phase ratio. The molecular difference method enables to measure transfer enthalpy for the molecular difference between two solutes. The method can be applied for similar molecules possessing access to the same stationary phase volume. In this case, we will be able to assume that they encounter the same thermodynamic phase ratio. Chester and Coym’s method is valid at constant pressure at which compression effects of the stationary phase cannot be observed. The slope of the relationship between ln α versus 1/T equals:
s l o p e = ( ln α ( 1 / T ) ) = Δ H d 0 R
where ΔH°d is the partial molar enthalpy of transfer of the molecular difference for each of two selected series members; α is the quotient of the retention coefficients (k) for selected pairs of compounds showing structural differences (α = k2/k1, where k2 > k1).
For most distinguished pairs, some degree of the phase-ratio dependence causes the not linear shape of ln α = f(1/T) relationships (Figure 6). It suggests that they access to different stationary phase volume. Therefore, we selected only the linear dependencies of ln α versus 1/T for both eluent systems studied, characterized by the determination coefficients higher than 0.97 (Figure 7).
Figure 7 presents the relationship between ln α for analytes differing with the heteroatoms at R1 substituents (TPR 14/16 and 28/30) in the imidazole ring and selectivity between TPR 4/22 and TPR 10/34 (since the molecular difference between them is a methylene group), against the inverse temperature (1/T). When analyzing Figure 6, it can be concluded that chosen pairs access the same stationary phase volume. This suggests also that the transfer enthalpy for heteroatoms: sulfur and oxygen being their molecular difference, differs from each other. Since the selectivity plot is not affected by the phase ratio, the presence of the heteroatoms is responsible for access of these compounds to the stationary phase.
The transfer enthalpy from the mobile phase containing methanol to the teicoplanin stationary phase for the heteroatoms calculated from the slope of ln α against 1/T is about −3.77 kJ/mol for series without methylene linker and −5.46 kJ/mol for series containing methylene linker. For pairs differing in methylene linker: TPR 4/22 and TPR 10/34, enthalpy of transfer varies from −1.71 to −3.28 kJ/mol. Corresponding values of the transfer enthalpy from mobile phase containing ACN are 30–40% less negative, and they are expressed by the following for heteroatom selectivity: −2.80 kJ/mol (TPR14/16), −3.29 kJ/mol (TPR 28/30), and for methylene selectivity: −0.93 kJ/mol (TPR 4/22), −2.26 kJ/mol (TPR 10/34).

3.5. Cross-Correlation of Chromatographic and Thermodynamic Parameters with Theoretical Structural Descriptors

In the aim to find a descriptor which is responsible for the interactions of the tested molecules with teicoplanin and ultimately contributes their retention on the column, cross-correlations of chromatographic and thermodynamic parameters with non-empirical topological and quantum chemical indices have been performed. In Table 4, cross-correlation results represented by the determination coefficient (R2) have been collected for the relationship between experimental retention parameter log kw and the Gibbs free energy change ΔG° estimated in each eluent system studied and theoretical structural descriptors, such as topological and quantum chemical indicators.
The statistically significant correlation was obtained for polarizability values and ΔG° estimated for the methanol-containing mobile phase. It could be concluded that the strength of the interactions between teicoplanin and examined molecules is related to their polarizability, in that case. It means that the induced dipole moment of examined molecules formed in the methanol/water solvent system is associated with an appropriate and proportional strength of interactions with teicoplanin. This conclusion could be explained by solvent polarity expressed by the empirical ENT scale. In this scale, methanol polarity is 0.762 whereas acetonitrile is 0.046. The polarity of water was assumed to be 1. So, methanol in water will have higher solvation properties in comparison to acetonitrile/water.
These correlations are in line with previous suggestions regarding dipole moment induction by the polar environment. What is more, the theoretically calculated polarization value seems to be a good tool in predicting the retention properties of tested molecules on a teicoplanin column with a methanol/water solvent, and thus, to assess the thermodynamics of absorption of these new therapeutic substances. It should also be noted that the low correlation with the BBB parameter representing the possibility of crossing the blood-brain barrier confirms the observation of other researchers that the passage of BBB is not always associated with anticonvulsant activity. To examine the mechanism of interaction between teicoplanin and the derivatives tested, a docking simulation was performed.

3.6. Study of Interactions of Teicoplanin with Investigated Triazoles by Docking Simulation

The density functional theory (DFT) approach and docking simulations were used to get the molecular interpretation and prove the obtained experimental results. The binding energies found during docking simulations and DFT calculations are given in Table 5 and graphically illustrated in Figure 7. In general, the obtained values exhibit very similar magnitude and vary from each other by less than 1 kJ/mol. The latter value is well below the accuracy threshold characteristic of the majority of standard molecular modeling methods, including those applied in our study. Thus, we refrain from a detailed, quantitative analysis of the data collected in Table 5. However, it is worth noting that both the magnitude of the calculated binding energies as well as their relatively minor scatter are in line with the corresponding experimental data (see Table 3 and Figure 3). Let us mention that the docking procedure included the whole teicoplanin molecule and resulted in nearly uniform binding modes, depending solely on the molecular topology of the studied ligands (see discussion below); this additionally indicates the meaningfulness of the obtained results.
The energies extracted from the DFT calculations exhibit a much higher correlation with experimentally-determined energy values (see Figure 8). However, the determined values are systematically higher, which can be ascribed to the lack of both the explicit solvent and any additional protocol to compensate the effects following this approximation. Such compensating approaches are inherent to the Vina-based docking, described above, which is probably the main reason for the discrepancies between DFT and docking-based energies.
The results of both the docking and DFT studies have been analyzed with respect to the mechanistic interaction pattern that may be significant in the context of the interpretation of the obtained binding energy values. The graphical illustration of the docking results is given in Figure 9A,B whereas the exemplary visualization of the DFT-based geometry optimization is given in Figure 9C,D. In general, the two distinct binding patterns were identified during the docking procedure which were maintained after subsequent DFT calculations. The first of them is exhibited by compounds with direct (i.e., via single covalent bond) linkage between triazole and phenyl rings (TRP-2-TRP-16) whereas the second, by compounds that contain the additional methylene moiety between those two rings (TRP-20-TRP-34). Within each of these groups, the binding mode is very similar and the main differences concern only the region of the triazole ring substituent (−R in Figure 9), the type of which differs from one compound to another. These small structural divergences may be responsible for the resulting minor differences in binding energies. Moreover, the mentioned ring substituent does not exhibit any direct contact with the teicoplanin molecule which, additionally, can be associated with its limited influence on the total binding energy. Finally, the –R substituents are usually flexible, which implies a high entropic component of the associated binding free energies; this factor is hardly reproduced by the static docking simulations and DFT calculations.
As identified by the docking study, the binding pattern exhibited by the first group of compounds is characterized by the close contact of the ligand chlorophenyl group with the dihydroxy phenyl moiety of teicoplanin. The parallel arrangement of these two rings indicates the presence of π-π interactions. At the same time, the triazole ring is rotated by about 60 deg with respect to the phenyl ring and maintains close contact with aliphatic and aromatic hydrogen atoms that belong to the hydroxyl methylene and chlorophenyl groups of teicoplanin, respectively. The arrangement of those contacts speaks for the presence of the attractive H-π stacking interactions. The remaining part of the ligand molecule (i.e., the −R substituent) does not exhibit any close contacts with teicoplanin but rather tends to create intramolecular interactions with the chlorophenyl group.
Most of the above-mentioned patterns are maintained after DFT-based geometry optimization. The most important differences include: (i) The disruption of the initially parallel arrangement of the ligand chlorophenyl group with respect to the dihydroxy phenyl moiety of teicoplanin. The twisted alignment speaks for the contribution of both π-π interactions and H-π stacking; (ii) creating the hydrogen bonding between the triazole ring and the carboxyl group of teicoplanin. At the same time, the −R substituent is moved even further from the teicoplanin molecule but this type of interaction does not seem to be relevant for the main binding driving force.
Ligands in the second group contain additional methylene moiety between two ring structures in their molecule. This has significant implications for the binding pattern in comparison to the first group, discussed above. Namely, now the triazole ring interacts by the π-π interactions with the chlorophenyl group of teicoplanin. The attractiveness of this interaction pattern is additionally enhanced by the proximity of the hydroxyl group of the teicoplanin with electronegative nitrogen atoms that create the triazole ring. The second ring moiety creates the H-π stacking interactions that involve aliphatic and aromatic hydrogen atoms belonging to the hydroxy methylene and chlorophenyl groups of teicoplanin, respectively. Again, the –R substituent does not maintain any close contacts with teicoplanin molecule but, in the case of larger substituents, may interact with the chlorophenyl group of its own molecule.
The DFT calculations reveal a series of slight changes in the described above patterns of interactions. Namely, the triazole ring arrangement with respect to the chlorophenyl group of teicoplanin is twisted in relation to the initially parallel arrangement. This speaks for the larger contribution of the H-π stacking interactions involving those two groups. Secondly, the hydroxyphenyl group of the teicoplanin forms a fully developed hydrogen bond with the triazole ring. At the same time, the −R substituent and the group responsible for the H-π stacking interactions identified during docking simulations are moved further away from the teicoplanin molecule, which results in reducing the magnitude of the corresponding interactions.
Summarizing, the interactions of the studied compounds with teicoplanin exhibit well-defined structural patterns, common for all compounds following the same molecular topology. The main driving force for binding seems to be the π-π and H-π interactions, supported by hydrogen bonding which involves the triazole ring. The limited contacts between the −R substituent and teicoplanin explain the small scatter of binding free energies determined across the studied group of compounds.
As can be seen in Table 3; Table 5, Gibbs free energy values obtained by docking simulation differ from those obtained experimentally in a chromatographic system. The docking-derived binding energies exhibit a very weak correlation with those obtained experimentally, however, some analogy between them can be noticed (Figure 8). The ΔΔG° values representing teicoplanin selectivity in relation to the heteroatom type (oxygen, sulfur), i.e., ΔΔG° TPR 14/16 and ΔΔG° TPR 28/30 are 0.23 and 0.25 kJ/mol in docking, 1.34 and 1.37 kJ/mol in the methanol system, and in the system with acetonitrile, 1.02 and 1.03 kJ/mol. As for the selectivity towards the methylene group, the values obtained by docking simulation for pairs TPR 4/22 and TPR 10/34 are 0.35 and 0.62 kJ/mol, respectively, while obtained in the methanol system is 0.65 kJ/mol for both pairs of compounds, and in the system with acetonitrile is 0.4 kJ/mol. The ΔΔG° values obtained experimentally are a little higher than those obtained by the docking method, which once again confirms the contribution of the solvent involved in the interaction process.
In contrast, correlations of the DFT results versus chromatographically-obtained ΔG° in both eluent systems are surprisingly consistent. The R2 values are the following: 0.8595 for MeOH- containing eluent and a little bit worse for ACN, namely 0.6814 (Figure 8).

3.7. Statistical Analysis of Thermodynamic and Chromatographic Parameters for Activity Evaluation

All collected data on chromatographic retention, thermodynamic parameters, and anticonvulsant activity expressed as IC50 were subjected to statistical analysis to find the method and parameter most useful for screening candidates for anticonvulsants.

3.7.1. Activity Division

Figure 10 illustrates the distribution of the activity parameter. The clear bi-modal distribution suggests that there is a mixture of two populations of observation, one with the lower activity described by higher ED50 parameter, and the second of higher activity individuals with lower ED50. Based on the distribution, the optimal discrimination boundary for ED50 is 350.

3.7.2. Comparison of Classification Models

During Regularized Logistic Regression (PLR) model tuning process, best accuracy (0.83) was obtained for λ = 0.1 and α ≈ 0.016 on the training dataset. Validation procedure conducted on test samples shows 0.8 model prediction, so quite similar to the accuracy of train samples. Therefore, we can say our model is not overfitted (Table 6).
Two methods of feature selection (Recursive Feature Elimination (RFE) and Simulated Annealing (SA) Feature Selection) were applied to extract the best set of predictors. Considering both selection methods, the following set of predictors was selected as optimal: RFE: (MeOH), S (MeOH), (ACN), log k (25% ACN); and SA: G (15% MeOH), log k (20%MeOH), log k (10% ACN), log k (15% ACN), log k 25% ACN), log kw (MeOH). Better precision of the RF model on train samples was obtained for the RFE algorithm of the initial selection of variables. So, the final RF model was built with RFE selection.
Calibration of the PLS-DA model shows that the best accuracy was obtained for two components (Figure 11A), which gathered information from all explanatory variables. The overall model precision was 0.84 on the training dataset and 0.8 on the test samples (Table 6).
A comparison of the model’s quality was shown on the plot representing the distribution of residuals (Figure 11A). Figure 11A shows that the PLS-DA model slightly outperforms the rest of them in a sense of residuals, but also has the highest variance.
The two classification models (PLR, RF) are black boxes, which means it is hard to trace the link between input variables and model outcomes [53]. Only regularized logistic regression (PR) is the kind of model in which the relation between explanatory and dependent variables is evident. To unify the assessment of variable importance, universal method was used [44]. All applied classification models (PLS-DA, PR, RF) confirm that the ϕ0 MeOH value is the most important predictor among all. The second most important is sMeOH (Figure 11B).

3.7.3. Regression Model Comparisons

Similarly, to the classification approach, RF models were estimated for predictors from two types of feature selection. One could be noticed that feature algorithms were applied to the same data, but in the regression task, so the sets of variables chosen by them are different than in classification. Again, the RF model with initial RFE selection gives the best prediction. RMSE for the final model on the training dataset was equal to 145.72 but on test samples, an even better RMSE = 106.17.
The calibration of PLS regression shows that the best model was obtained for two components, where RMSE = 154.72 on the training dataset. Overall performance for test samples was better, because RMSE = 115.74 with R2 = 0.48, and MAE = 99.53.
Finally, estimation of Penalized Regression parameters was conducted and the best combination of hyperparameters was α = 10−5 and λ = 0.001. RMSE for the final model was 142.06 on the training dataset. RMSE on test samples was even lower at 112.20 with R2 = 0.44 and MAE = 92.75.
Based on the distribution of residuals, evaluated for each model (RF, PR, PLS), one could notice that RF is the best regression model because it outperforms the rest of the models in the sense of distribution of residuals (Figure 12A). Assessment of the variable importance in the regression models again pointed out that the most important variable is ϕ0 MeOH and the second most important is S MeOH (Figure 12B). Key variables were common to both types of models.

4. Conclusions

The drug, after administration, contacts many different biological membranes, i.e., macrophage cells, vascular endothelium, and complex absorption barriers. The interaction of the biomolecule with these membranes affects its further biodistribution in the body. To date, due to the complexity of the structure and function of various biological barriers, many simpler models have been proposed to study the in vitro interaction between bioactive compounds and biological membrane models. Due to the analogy between interactions occurring in the biological environment and chromatographic systems, the potential of the teicoplanin-modified HPLC column for predicting physicochemical parameters significant in ADMET processes has been studied in the current work.
In this paper, the chromatographic system with the teicoplanin modified column (Chirobiotic T) and organic-aqueous mobile phase was used to model the biological membrane. Although teicoplanin-bonded silica was first introduced as an enantioselective stationary phase, it has shown that it may also be useful to the analysis of achiral molecules.
Thermodynamic parameters obtained experimentally and calculated based on molecular modeling consistently confirm weak interactions of 1,2,4-triazolo-3-thiones with teicoplanin moiety. This is evidenced especially by small ΔG° values: theoretically estimated by the docking simulation in the range of −5.38 to −4.65 kJ/mol, experimentally on teicoplanin column and methanol/water eluent system in the range of −7.64 to −5.23 kJ/mol, and acetonitrile/water eluent system in the range of −6.02 to −4.16 kJ/mol. Despite this, the DFT-derived binding energies are systematically higher, and their correlation with experimental data is better (R2 = 0.8595 for MeOH, and 0.6814 for ACN). The weak interaction behavior of investigated analytes reflects the results of earlier pharmacological studies that showed poor permeability of these compounds through the blood-brain barrier despite the confirmed anticonvulsant activity.
The experimentally measured values of the so-called chromatographic lipophilicity parameter (log kw) are also significantly smaller than theoretically predicted, which further confirms the ability of teicoplanin to predict the behavior of the studied xenobiotics at the interface of biological barriers.
The usefulness of this model should be emphasized, especially with a methanol/water eluent system which is able to very precisely reflect the polarizability of molecules. This parameter appears to be important for the future selection of anticonvulsant drug candidates.
Teicoplanin generates both polar and hydrophobic chemical interactions, ensuring not only a partitioning phenomenon but also a binding pathway. Only on the basis of measurements of retention in the teicoplanin/methanol-water system, it is possible to initially anticipate both pharmacokinetic aspects of the behavior of the potential drug in vivo. It was proved by PLS-DA, that there is a potential to discriminate between active and inactive 1,2,4-triazole-3-thiones due to the chromatographically measured variables (φ0 MeOH, S MeOH). Summarizing, it can be concluded that the teicoplanin HPLC column possesses the potential to estimate 1,2,4-triazole-3-thiones’ antiepileptic properties.

Author Contributions

Conceptualization, J.F. and A.O.; methodology, J.F., A.O.; software, W.P., D.M.; validation, A.K.-K., M.P. and H.T.; formal analysis, A.K.-K., M.P. and H.T.; investigation, A.K.-K. and M.P.; resources, T.P.; data curation, D.M.; writing—original draft preparation, J.F.; writing—review and editing, J.F. and W.P.; visualization, W.P.; supervision, J.F.; project administration, J.F.; funding acquisition, J.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kirchmair, J. Drug Metabolism Prediction, 1st ed.; Wiley-VCH Verlag GmbH & Co. KGaA: Weinheim, Germany, 2014. [Google Scholar]
  2. Pliska, V.; Testa, B.; van de Waterbeemd, H. Lipophilicity in Drug Action and Toxicology; John Wiley & Sons: Weinheim, Germany, 1996. [Google Scholar]
  3. Kaliszan, R. Structure and Retention in Chromatography. A Chemometric Approach; Harwood Academic Publishers: Amsterdam, The Netherlands, 1997. [Google Scholar]
  4. Van de Waterbeemd, H.; Carter, R.E.; Grassy, G.; Kubinyi, H.; Martin, Y.C.; Tute, M.S.; Willett, P. Glossary of terms used in computational drug design (IUPAC Recommendations 1997). Pure Appl. Chem. 1997, 69, 1137. [Google Scholar] [CrossRef]
  5. Lambert, W.J. Modeling Oil-Water Partitioning and Membrane Permeation Using Reversed-Phase Chromatography. J. Chromatogr. A 1993, 656, 469. [Google Scholar] [CrossRef]
  6. Dorsey, J.G. Implications for Biological Partitioning Processes. J. Chromatogr. A 1993, 656, 485. [Google Scholar] [CrossRef]
  7. Kaliszan, R. Effect of Separation Conditions on Chromatographic Determination of Hydrophobicity of Acidic Xenobiotics. J. Chromatogr. B 1998, 717, 125. [Google Scholar] [CrossRef]
  8. Markuszewski, M.; Kaliszan, R. Quantitative Structure-Retention Relationships in Affinity High-Performance Liquid Chromatography. J. Chromatogr. B 2002, 768, 55. [Google Scholar] [CrossRef]
  9. Flieger, J.; Tatarczak-Michalewska, M.; Wujec, M.; Pitucha, M.; Świeboda, R. RP-HPLC analysis and in vitro identification of antimycobacterial activity of novel thiosemicarbazides and 1,2,4-triazole derivatives. J. Pharm. Biomed. Anal. 2015, 107, 501. [Google Scholar] [CrossRef]
  10. Ong, S.; Liu, H.; Pidgeon, C. Immobilized-Artificial Membrane Chromatography: Measurements of Membrane Partition Coefficient and Predicting Drug Membrane Permeability. J. Chromatogr. A 1996, 728, 113. [Google Scholar] [CrossRef]
  11. Yang, C.Y.; Cai, S.J.; Liu, H.; Pidgeon, C. Immobilized Artificial Membranes-Screens for Drug Membrane Interactions. Adv. Drug Deliv. Rev. 1997, 23, 229. [Google Scholar] [CrossRef]
  12. Lepont, C.; Poole, C.F. Retention Characteristics of an Immobilized Artificial Membrane Column in Reversed-Phase Liquid Chromatography. J. Chromatogr. A 2002, 946, 107. [Google Scholar] [CrossRef]
  13. Turnhofer, H.; Schnabel, J.; Betz, M.; Lipka, G.; Pidgeon, C.; Hauser, H. Cholesterol transfer protein located in the intestinal brush-border membrane. Partial purification and characterization. Biochim. Biophys. Acta 1991, 1064, 275. [Google Scholar] [CrossRef]
  14. Taillardat-Bertschinger, A.; Carrupt, P.A.; Barbato, F.; Testa, B. Immobilized Artificial Membrane HPLC in Drug Research. J. Med. Chem. 2003, 46, 1. [Google Scholar] [CrossRef] [PubMed]
  15. Flieger, J.; Trębacz, H.; Pizoń, M.; Kowalska, A.; Szczęsna, A.; Plech, T. High-performance liquid chromatography thermodynamic study of new potential antiepileptic compounds on a cholesterol column using isocratic elution with methanol/water and acetonitrile/water eluent systems. J. Sep. Sci. 2017, 40, 4176. [Google Scholar] [CrossRef] [PubMed]
  16. Flieger, J.; Pizon, M.; Plech., T. Chromatographic behavior of new antiepileptic active compounds on different reversed-phase materials. J. Chromatogr. A 2014, 1338, 188. [Google Scholar] [CrossRef] [PubMed]
  17. Hewen Li, H.; Zhao, T.; Sun, Z. Analytical techniques and methods for study of drug-lipid membrane interactions. Rev. Anal. Chem. 2018, 37, 20170012. [Google Scholar]
  18. Pereiraleite, C.; Nunes, C.; Reis, S. Interaction of Nonsteroidal Anti-Inflammatory Drugs with Membranes: In Vitro Assessment and Relevance for Their Biological Actions. Prog. Lipid Res. 2013, 52, 571–584. [Google Scholar] [CrossRef] [PubMed]
  19. Moravcova, D.; Planeta, J.; Wiedmer, S.K. Silica-Based Monolithic Capillary Columns Modified by Liposomes for Characterization of Analyte-Liposome Interactions by Capillary Liquid Chromatography. J. Chromatogr. A 2013, 1317, 159–166. [Google Scholar] [CrossRef] [PubMed]
  20. Zhang, C.; Li, J.; Xu, L.; Shi, Z. Fast Immobilized Liposome Chromatography Based on Penetrable Silica Microspheres for Screening and Analysis of Permeable Compounds. J. Chromatogr. A 2012, 1233, 78–84. [Google Scholar] [CrossRef] [PubMed]
  21. Liu, X.; Fan, P.; Chen, M.; Hefesha, H.; Scriba, G.K.E.; Gabel, D.; Fahr, A. Drug-Membrane Interaction on Immobilized Liposome Chromatography Compared to Immobilized Artificial Membrane (IAM), Liposome/Water, and Octan-1-ol/Water Systems. Helv. Chim. Acta 2010, 93, 203–211. [Google Scholar] [CrossRef]
  22. Baldwin, T.D.; Kelly, J.W. Application of Enzyme Biotechnology; Plenum Press: New York, NY, USA, 1992. [Google Scholar]
  23. Pidgeon, C.; Ong, S.; Liu, H.; Qiu, X.; Pidgeon, M.; Dantzig, A.H.; Munroe, J.; Hornback, W.J.; Kasher, J.S.; Glunz, L.; et al. IAM Chromatography: An in Vitro Screen for Predicting Drug Membrane Permeability. J. Med. Chem. 1995, 38, 590. [Google Scholar] [CrossRef] [PubMed]
  24. Stewart, B.H.; Chan, H. Use of Immobilized Artificial Membrane Chromatography for Drug Transport Applications. J. Pharm. Sci. 1998, 87, 1471. [Google Scholar] [CrossRef] [PubMed]
  25. Schäfer-Korting, M. Drug Delivery; Springer: Berlin/Heidelberg, Germany, 2010; p. 197. [Google Scholar]
  26. Thuerauf, N.; Fromm, M.F. The role of the transporter P-glycoprotein for disposition and effects of centrally acting drugs and for the pathogenesis of CNS diseases. Eur. Arch. Psychiatry Clin. Neurosci. 2006, 256, 281. [Google Scholar] [CrossRef] [PubMed]
  27. Persidsky, Y.; Ramirez, S.H.; Haorah, J.; Kanmogne, G.D. Blood-brain barrier; structural components and function under physiologic and pathologic conditions. J. Neuroimmune Pharmacol. 2006, 1, 223. [Google Scholar] [CrossRef] [PubMed]
  28. Marchi, N.; Betto, G.; Fazio, V.; Fan, Q.; Ghosh, C.; Machado, A.; Janigro, D. Blood-brain barrier damage and brain penetration of antiepileptic drugs: Role of serum proteins and brain edema. Epilepsia 2009, 50, 664. [Google Scholar] [CrossRef] [PubMed]
  29. Carpenter, T.S.; Kirshner, D.A.; Lau, E.Y.; Wong, S.E.; Nilmeier, J.P.; Lightstone, F.C. A Method to Predict Blood-Brain Barrier Permeability of Drug-Like Compounds Using Molecular Dynamics Simulations. Biophys. J. 2014, 107, 630. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Pignatello, R.; Musumeci, T.; Basile, L.; Carbone, C.; Puglisi, G. Biomembrane models and drug-biomembrane interaction studies: Involvement in drug design and development. J. Pharm. Bioall. Sci. 2011, 3, 4. [Google Scholar] [CrossRef] [PubMed]
  31. Armstrong, D.W.; Tang, Y.; Chen, S.; Zhou, Y.; Bagwill, C.; Chen, J. Macrocyclic antibiotics as a new class of chiral selectors for liquid chromatography. Anal. Chem. 1994, 66, 1473. [Google Scholar] [CrossRef]
  32. Ward, T.J.; Farris, A.B. Chiral separations using the macrocyclic antibiotics: A review. J. Chromatogr. A 2001, 906, 73. [Google Scholar] [CrossRef]
  33. Ilisz, I.; Berkecz, R.; Péter, A. Retention mechanism of high-performance liquid chromatographic enantioseparation on macrocyclic glycopeptide-based chiral stationary phases. J. Chromatogr. A 2009, 1216, 1845. [Google Scholar] [CrossRef] [PubMed]
  34. Aburachid Cardoso, P.; Costa César, I. Chiral method development strategies for HPLC using macrocyclic glycopeptide-based stationary phases. Chromatographia 2018, 81, 841. [Google Scholar] [CrossRef]
  35. Moldoveanu, S.C.; David, V. Stationary Phases and Columns for Chiral Chromatography. HPLC Method in Chemical Analysis; Elsevier: Amsterdam, The Netherlands, 2017; pp. 363–376. [Google Scholar]
  36. Tang, M.; Zhang, J.; Zhuang, S.; Liu, W. Development of chiral stationary phases for high-performance liquid chromatographic separation. Trends Anal. Chem. 2012, 39, 180. [Google Scholar] [CrossRef]
  37. Kaproń, B.; Łuszczki, J.J.; Płazińska, A.; Siwek, A.; Karcz, T.; Gryboś, A.; Nowak, G.; Makuch-Kocka, A.; Walczak, K.; Langner, E.; et al. Development of the 1,2,4-triazole-based anticonvulsant drug candidates acting on the voltage-gated sodium channels. Insights from in-vivo, in-vitro, and in-silico studies. Eur. J. Pharm. Sci. 2019, 129, 42. [Google Scholar] [PubMed]
  38. Plech, T.; Łuszczki, J.J.; Wujec, M.; Flieger, J.; Pizoń, M. Synthesis, characterization and preliminary anticonvulsant evaluation of some 4-alkyl-1,2,4-triazoles. Eur. J. Med. Chem. 2013, 60, 208. [Google Scholar] [CrossRef] [PubMed]
  39. Plech, T.; Kaproń, B.; Łuszczki, J.J.; Paneth, A.; Siwek, A.; Kołaczkowski, M.; Żołnierek, M.; Nowak, G. Studies on the anticonvulsant activity of 4-alkyl-1,2,4-triazole-3-thiones and their effect on gabaergic system. Eur. J. Med. Chem. 2014, 86, 690. [Google Scholar] [CrossRef] [PubMed]
  40. Kaproń, B.; Łuszczki, J.J.; Siwek, A.; Karcz, T.; Nowak, G.; Zagaja, M.; Andres-Mach, M.; Stasiłowicz, M.; Cielecka-Piontek, J.; Kocki, J. Preclinical evaluation of 1,2,4-triazole-based compounds targeting voltage-gated sodium channels (VGSCs) as promising anticonvulsant drug candidates. Bioorg. Chem. 2020, 94, 103355. [Google Scholar] [CrossRef] [PubMed]
  41. Kaproń, B.; Łuszczki, J.; Paneth, A.; Wujec, M.; Siwek, A.; Karcz, T.; Mordyl, B.; Głuch-Lutwin, M.; Gryboś, A.; Nowak, G. Molecular mechanism of action and safety of 5-(3-chlorophenyl)-4-hexyl-2,4-dihydro-3H-1,2,4-triazole-3-thione—A novel anticonvulsant drug candidate. Int. J. Med. Sci. 2017, 14, 741. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Flieger, J.; Trębacz, H.; Pizoń, M.; Płazińska, A.; Płaziński, W.; Kowalska, A.; Szczęsna, A.; Plech, T. Thermodynamic study of new antiepileptic compounds by combining chromatography on the phosphatidylcholine biomimetic stationary phase and differential scanning calorimetry. J. Sep. Sci. 2019, 42, 2628. [Google Scholar] [CrossRef] [PubMed]
  43. R Core Team. A Language and Environment for Statistical Computing. 2020. Available online: https://www.r-project.org/ (accessed on 2 May 2020).
  44. Kuhn, M. Caret: Classification and Regression Training. 2020. Available online: https://CRAN.R-project.org/package=caret (accessed on 2 May 2020).
  45. Biecek, P. DALEX: Explainers for Complex Predictive Models in R. J. Mach. Learn. Res. 2018, 19, 1. [Google Scholar]
  46. Kassambara, A. Ggpubr: ’Ggplot2’ Based Publication Ready Plots. Available online: https://CRAN.R-project.org/package=ggpubr (accessed on 2 May 2020).
  47. Liaw, A.; Wiener, M. Classification and Regression by random Forest. R. News 2002, 2, 18. [Google Scholar]
  48. Friedman, J.; Hastie, T.; Tibshirani, R. Regularization Paths for Generalized Linear Models via Coordinate Descent. J. Stat. Softw. 2010, 33, 1. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Noah, S.; Friedman, J.; Hastie, T.; Tibshirani, R. Regularization Paths for Cox’s Proportional Hazards Model via Coordinate Descent. J. Stat. Softw. 2011, 39, 1. [Google Scholar]
  50. Mevik, B.-H.; Wehrens, R.; Hovde Liland, K. Pls: Partial Least Squares and Principal Component Regression. Available online: https://CRAN.R-project.org/package=pls (accessed on 2 May 2020).
  51. Granitto, P.M.; Furlanello, C.; Biasioli, F.; Gasperi, F. Recursive Feature Elimination with Random Forest for PTR-MS Analysis of Agroindustrial Products. Chemom. Intell. Lab. 2006, 83, 83. [Google Scholar] [CrossRef]
  52. Kirkpatrick, S.; Gelatt, C.D.; Vecchi, M.P. Optimization by Simulated Annealing. Science 1983, 220, 671. [Google Scholar] [CrossRef] [PubMed]
  53. Browne, M.W. Cross-Validation Methods. J. Math. Psychol. 2000, 44, 108. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Tibshirani, R.; Bien, J.; Friedman, J.; Hastie, T.; Simon, N.; Taylor, J.; Tibshirani, R.J. Strong Rules for Discarding Predictors in Lasso-Type Problems: Strong Rules for Discarding Predictors. J. R. Stat. Soc. 2012, 74, 245. [Google Scholar] [CrossRef] [PubMed]
  55. Hanwell, M.D.; Curtis, D.E.; Lonie, D.C.; Vandermeersch, T.; Zurek, E.; Hutchison, G.R. Avogadro: An advanced semantic chemical editor, visualization, and analysis platform. J. Cheminform. 2012, 4, 17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Rappe, A.K.; Casewit, C.J.; Colwell, K.S.; Goddard, W.A.; Skiff, W.M. UFF, a full periodic table force field for molecular mechanics and molecular dynamics simulations. J. Am. Chem. Soc. 1992, 114, 10024. [Google Scholar] [CrossRef]
  57. Trott, O.; Olson, A.J. AutoDock Vina: Improving the Speed and Accuracy of Docking with a New Scoring Function, Efficient Optimization, and Multithreading. J. Comput. Chem. 2010, 31, 455. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Lee, C.; Yang, W.; Parr, R.G. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B 1988, 37, 785. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Hehre, W.J.; Radom, L.; Schleyer, P.R.; Pope, J. Ab Initio Molecular Orbital Theory; Wiley: New York, NY, USA, 1986; pp. 86–87. [Google Scholar]
  60. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09, Revision E.01; Gaussian, Inc.: Wallingford, UK, 2009. [Google Scholar]
  61. Soczewinski, E.; Wachtmeister, C.A. The relation between the composition of certain ternary two-phase solvent systems and RM values. J. Chromatogr. 1962, 7, 311. [Google Scholar] [CrossRef]
  62. Valko, K.; Sleger, P. New chromatographic hydrophobicity index (&) based on the slope and the intercept of the logk’ versus organic phase concentration plot. J. Chromatogr. A 1993, 631, 49. [Google Scholar]
  63. Barbaro, G.L.; Sapone, A.M.; Borea, A.; Variani, P.A.; Recanatini, M. Study of lipophilic character of serotonic ligand. J. Chromatogr. A 1996, 723, 135. [Google Scholar]
  64. Biagi, L.; Barbaro, A.M.; Sapone, A.; Recanatini, M. Determination lipophilicity by means of reversed phase thin-layer chromatography. Influence of the organic modifier on the slope of the thin-layer chromatographic equation. J. Chromatogr. A 1994, 669, 246. [Google Scholar] [CrossRef]
  65. Knox, H.; Pryde, A. Performance and selected applications of a new range of chemically bonded packing materials in high-performance liquid chromatography. J. Chromatogr. 1975, 112, 171. [Google Scholar] [CrossRef]
  66. Kazakevich, Y.V.; LoBrutto, R.; Chan, F.; Patel, T. Interpretation of the excess adsorption isotherms of organic eluent components on the surface of reversed-phase adsorbents: Effect on the analyte retention. J.Chromatogr. A 2001, 913, 75. [Google Scholar] [CrossRef]
  67. Hansch, C.; Leo, A.; Hoekman, D. Exploring QSAR—Hydrophobic, Electronic, and Steric Constants; American Chemical Society: Washington, DC, USA, 1995. [Google Scholar]
  68. Wildman, S.A.; Crippen, G.M. Prediction of physicochemical parameters by atomic contributions. J. Chem. Inf. Comput. Sci. 1999, 39, 868. [Google Scholar] [CrossRef]
  69. Chester, T.L.; Coym, J.W. Effect of phase ratio on van’t Hoff analysis in reversed-phase liquid chromatography, and phase-ratio independent estimation of transfer enthalpy. J. Chromatogr. A 2003, 1003, 101. [Google Scholar] [CrossRef]
  70. Dorsey, J.G.; Dill, K.A. The molecular mechanism of retention in reversed-phase liquid chromatography. Chem. Rev. 1989, 89, 331. [Google Scholar] [CrossRef]
  71. Sentell, K.B.; Dorsey, J.G. On the Calculation of the Stationary Phase Volume in Reversed Phase Chromatography. J. Liq. Chromatogr. 1988, 11, 1875. [Google Scholar] [CrossRef]
  72. Doyle, C.A.; Vickers, T.J.; Mann, C.K.; Dorsey, J.G. Characterization of C18-bonded liquid chromatographic stationary phases by Raman spectroscopy: The effect of temperature. J. Chromatogr. A 2000, 877, 41. [Google Scholar] [CrossRef]
Scheme 1. Schematic representations of stationary phases (A) simulating biological membrane (B). Surface-confined membranes at the top were obtained by covalently bonding phosphatidylcholine octadecylated silica (ODS), liposomes, cholesterol, proteins, and teicoplanin (the aim of current work) to porous silica.
Scheme 1. Schematic representations of stationary phases (A) simulating biological membrane (B). Surface-confined membranes at the top were obtained by covalently bonding phosphatidylcholine octadecylated silica (ODS), liposomes, cholesterol, proteins, and teicoplanin (the aim of current work) to porous silica.
Materials 13 02650 sch001
Figure 1. Plots of log k against φ (v/v) methanol (MeOH), acetonitrile (ACN) in the mobile phase obtained on CHIROBIOTICTM column.
Figure 1. Plots of log k against φ (v/v) methanol (MeOH), acetonitrile (ACN) in the mobile phase obtained on CHIROBIOTICTM column.
Materials 13 02650 g001
Figure 2. The cross relationship between the chromatographic lipophilicity parameters: S (A), φ0 (B), and log kw (C), measured on chirobiotic column by the use of different eluent systems: ACN/water and MeOH/water.
Figure 2. The cross relationship between the chromatographic lipophilicity parameters: S (A), φ0 (B), and log kw (C), measured on chirobiotic column by the use of different eluent systems: ACN/water and MeOH/water.
Materials 13 02650 g002
Figure 3. Comparison of different lipophilicity scales: chromatographic (log kw MeOH, log kw ACN, and calculated theoretically (AC log P). Values were taken from the experimental database VCCLAB (ALIGPS 2.1).
Figure 3. Comparison of different lipophilicity scales: chromatographic (log kw MeOH, log kw ACN, and calculated theoretically (AC log P). Values were taken from the experimental database VCCLAB (ALIGPS 2.1).
Materials 13 02650 g003
Figure 4. Van’t Hoff plots between the natural logarithms of the retention factors vs. 1000/T for the examined compounds. (A) 15% methanol as the mobile phase; (B) 15% acetonitrile as the mobile phase.
Figure 4. Van’t Hoff plots between the natural logarithms of the retention factors vs. 1000/T for the examined compounds. (A) 15% methanol as the mobile phase; (B) 15% acetonitrile as the mobile phase.
Materials 13 02650 g004
Figure 5. Phase transitions for the dry stationary phase checked by DSC (differential scanning calorimetry).
Figure 5. Phase transitions for the dry stationary phase checked by DSC (differential scanning calorimetry).
Materials 13 02650 g005
Figure 6. The selectivity plot for distinguishing pairs that differ with methylene or methyl moieties, heteroatoms, or isomeric form examined on Chirobiotic column by the use of 15% ACN (A) and 15% MeOH (B).
Figure 6. The selectivity plot for distinguishing pairs that differ with methylene or methyl moieties, heteroatoms, or isomeric form examined on Chirobiotic column by the use of 15% ACN (A) and 15% MeOH (B).
Materials 13 02650 g006
Figure 7. The phase ratio-independent selectivity curves for analytes differing with the heteroatoms at R1 substituents (TPR 14/16 and 28/30) and the methylene group (TPR 4/22 and TPR 10/34) obtained on Chirobiotic column. Experimental conditions: teicoplanin column (CHIROBIOTICTM), 15% ACN (red lines) and 15% MeOH (black lines) mobile phases, at a flow rate of 1 mL/min.
Figure 7. The phase ratio-independent selectivity curves for analytes differing with the heteroatoms at R1 substituents (TPR 14/16 and 28/30) and the methylene group (TPR 4/22 and TPR 10/34) obtained on Chirobiotic column. Experimental conditions: teicoplanin column (CHIROBIOTICTM), 15% ACN (red lines) and 15% MeOH (black lines) mobile phases, at a flow rate of 1 mL/min.
Materials 13 02650 g007
Figure 8. The comparison of the Gibbs free energy change (ΔG°) estimated experimentally on Chirobiotic column with methanol (MeOH)- or acetonitrile (ACN)-containing eluent system and obtained by the use of docking simulation (insert) and DFT calculations.
Figure 8. The comparison of the Gibbs free energy change (ΔG°) estimated experimentally on Chirobiotic column with methanol (MeOH)- or acetonitrile (ACN)-containing eluent system and obtained by the use of docking simulation (insert) and DFT calculations.
Materials 13 02650 g008
Figure 9. (A,B) Graphical illustration of the docking results. The two groups of ligands are considered separately, i.e., those that either do (A) or do not (B) contain a methylene group between triazole and phenyl rings. The pictures show the identified locations of the ligand molecules (stick representation) interacting with teicoplanin (stick plus surface representation). (C,D) Graphical illustration of the DFT-optimized structures of complexes, shown on the example of the TRP-10 (C) and TRP-26 (D) compounds. The ligands are shown in the ball-and-stick representation.
Figure 9. (A,B) Graphical illustration of the docking results. The two groups of ligands are considered separately, i.e., those that either do (A) or do not (B) contain a methylene group between triazole and phenyl rings. The pictures show the identified locations of the ligand molecules (stick representation) interacting with teicoplanin (stick plus surface representation). (C,D) Graphical illustration of the DFT-optimized structures of complexes, shown on the example of the TRP-10 (C) and TRP-26 (D) compounds. The ligands are shown in the ball-and-stick representation.
Materials 13 02650 g009
Figure 10. The distribution of activity parameter.
Figure 10. The distribution of activity parameter.
Materials 13 02650 g010
Figure 11. Distribution of residuals (A) and variable importance for the set of optimal RFE (Recursive Feature Elimination) predictors (B) of three estimated models.
Figure 11. Distribution of residuals (A) and variable importance for the set of optimal RFE (Recursive Feature Elimination) predictors (B) of three estimated models.
Materials 13 02650 g011
Figure 12. Distribution of residuals (A), variable importance for set of optimal RFE predictors (B) of three estimated models.
Figure 12. Distribution of residuals (A), variable importance for set of optimal RFE predictors (B) of three estimated models.
Materials 13 02650 g012
Table 1. Investigated compounds.
Table 1. Investigated compounds.
General Formula (1)General Formula (2)
Materials 13 02650 i001 Materials 13 02650 i002
No.RNo.R
TPR-2
AC log P = 3.93 *
310.3 mg/kg bw **
Materials 13 02650 i003TPR-20
AC log P = 3.83 *
(−) **
Materials 13 02650 i004
TPR-4
AC log P = 4.33 *
221.7 mg/kg bw **
Materials 13 02650 i005TPR-22
AC log P = 4.23 *
130.4 mg/kg bw **
Materials 13 02650 i006
TPR-6
AC log P = 4.08 *
295.1 mg/kg bw **
Materials 13 02650 i007TPR-24
AC log P = 3.98 *
(−) **
Materials 13 02650 i008
TPR-8
AC log P = 4.79 *
(−) **
Materials 13 02650 i009TPR-26
AC log P = 4.69 *
(−) **
Materials 13 02650 i010
TPR-10
AC log P = 5.32 *
(−) **
Materials 13 02650 i011TPR-34
AC log P = 5.22 *
(−) **
Materials 13 02650 i012
TPR-14
AC log Ps = 3.49 *
208.9 mg/kg bw **
Materials 13 02650 i013TPR-28
AC log P = 3.39 *
132.4 mg/kg bw **
Materials 13 02650 i014
TPR-16 AC
log P = 4.16 *
201.0 mg/kg bw **
Materials 13 02650 i015TPR-30 AC
log P = 4.06 *
320.0 mg/kg bw **
Materials 13 02650 i016
* AC log P—The logarithm of the theoretical partition coefficients of the examined solutes were calculated using ALOGPS 2.1–VCCLab, available free of charge on the following website: www.vcclab.org. ** The anticonvulsant activity of the compounds was expressed as ED50 values. Each ED50 reflects the dose of a compound required to protect 50% of the animals tested against maximal electroshock (MES)-induced seizures. The data were obtained from our previous paper [39]. The symbol: (−) means that the compound is inactive.
Table 2. Statistical estimates for the linear relationships (log k vs. the volume fraction of the organic modifier: methanol and acetonitrile) obtained on the Chirobiotic T column.
Table 2. Statistical estimates for the linear relationships (log k vs. the volume fraction of the organic modifier: methanol and acetonitrile) obtained on the Chirobiotic T column.
SymbolMeOHACN
log kwSφ0R2Fselog kwSφ0R2Fse
TPR 21.3574−2.37280.57210.99024040.02471.3794−3.14280.43890.97841820.0488
TPR 41.4772−2.57120.57450.99255260.02341.4850−3.46470.42860.97961920.0523
TPR 61.4200−2.19860.64590.98042000.03251.4859−3.32700.44660.98432500.0440
TPR 81.6778−2.94970.56880.99376300.02461.6600−3.94920.42030.98632890.0486
TPR 101.7921−3.27120.54780.99569130.02261.8393−4.48320.41030.98562730.0567
TPR 141.4238−2.40810.59120.992025190.01261.4227−3.31810.42880.97323750.0419
TPR 161.7177−2.84690.60340.99609850.01901.6785−3.98310.42140.98182160.0567
TPR 201.4968−2.66510.56160.997415180.01431.4711−3.33630.44090.97651660.0541
TPR 221.6746−2.93380.57080.998322850.01281.5954−3.66360.43550.98152130.0526
TPR 241.6998−3.01680.56350.998425190.01261.6752−3.88150.43160.98953750.0419
TPR 261.8619−3.24600.57360.997918840.01561.8080−4.32760.41780.98582780.0543
TPR 281.5547−2.63480.59010.997918840.01561.5212−3.55730.42760.98582780.0543
TPR 301.8527−3.05140.60720.997214090.01701.7957−4.26440.42110.98512640.0549
TPR 342.0070−3.67460.54620.996611610.02261.9553−4.72970.41340.99154660.0458
Descriptive Statistics
Mean1.644−2.8460.5800.9951345.60.0191.627−3.8160.4270.983272.760.050
Median1.676−2.8900.5730.9961285.30.0181.628−3.7730.4280.985268.810.052
Standard deviation0.1880.3880.0250.005755.20.0060.1690.4760.0100.00581.190.005
Skewness0.211−0.3161.195−2.1600.1850.7090.367−0.4220.183−0.3990.972−0.557
Kurtosis-0.873−0.0881.9315.447−1.2400.022−0.916−0.9680.714−0.1320.723−1.171
Coefficient of variation11.85414.1404.5360.48758.24430.85110.78012.9532.4940.51130.8910.487
Table 3. Thermodynamic parameters obtained for investigated compounds by the use of the teicoplanin stationary phase and 15% organic-aqueous eluents with methanol or acetonitrile.
Table 3. Thermodynamic parameters obtained for investigated compounds by the use of the teicoplanin stationary phase and 15% organic-aqueous eluents with methanol or acetonitrile.
MeOHACN
No
TPR
∆H°
kJ/mol
∆S°^
kJ/mol*K
∆G°
kJ/mol
R2FSe∆H°
kJ/mol
∆S°
kJ/mol*K
∆G°
kJ/mol
R2FSe
2−20.6090−0.0521−5.23030.997928820.017−8.3367−0.0142−4.16760.8175270.072
4−20.5715−0.0502−5.76690.9999504840.004−7.7317−0.0111−4.48720.8404320.062
6−16.4897−0.0355−5.98240.98243340.040−4.11040.0022−4.74080.396040.093
8−23.5131−0.0569−6.74580.9999461200.005−8.8404−0.0129−5.06770.8828450.059
10−22.4397−0.0524−6.98650.9996169370.008−10.0121−0.0151−5.59320.8824450.067
14−23.4591−0.0598−5.82280.9995215500.007−9.8373−0.0194−4.16420.96701310.045
16−27.5706−0.0693−7.16380.9995129070.011−13.4231−0.0281−5.18330.96291560.048
20−21.4193−0.0534−5.65400.9997194020.007−10.2859−0.0198−4.47650.96681750.035
22−22.7078−0.0552−6.42150.9999444850.005−10.2560−0.0183−4.88150.96241540.037
24−22.0695−0.0528−6.48150.9997215500.007−11.9776−0.0231−5.19550.95621310.047
26−24.5931−0.0586−7.31470.9994108950.011−11.0007−0.0189−5.46150.95971430.041
28−25.2271−0.0643−6.27510.9995120460.010−11.2507−0.0232−4.46120.97792660.031
30−30.1099−0.0763−7.64800.9996134130.012−14.9553−0.0323−5.49330.98143170.038
34−25.7792−0.0615−7.64010.9999483300.005−13.1747−0.0244−6.02260.96181510.048
Table 4. Cross-correlations of ΔG° and log kw estimated in methanol and acetonitrile eluent system with the theoretical descriptors.*
Table 4. Cross-correlations of ΔG° and log kw estimated in methanol and acetonitrile eluent system with the theoretical descriptors.*
Log PhyperHEPolarizabilitySAVolµHFHOMOLUMOALOGPsAC logPALOG PMLOG PXLOG P2XLOG P3SVM BBBAdaBoost BBB
Cross-correlation with ΔG° estimated in methanol/water eluent system
0.0050.0000.9210.0370.6700.0430.0010.4060.0110.1430.3610.2820.3850.3020.3390.3970.119
Cross-correlation with ΔG° estimated in acetonitrile/water eluent system
0.1130.1050.6840.0060.3330.0000.1390.1670.1290.3940.6260.5840.6870.6850.6180.2150.007
Cross-correlation with log kw estimated in methanol/water eluent system
0.0360.0150.8860.0070.6070.0150.0030.4200.0030.2450.4270.3770.4580.4340.4170.3280.043
Cross-correlation with log kw estimated in acetonitrile/water eluent system
0.1100.0940.6890.0160.3350.0020.0680.2020.0650.4060.6420.5980.6810.6870.6300.2340.009
* The theoretical descriptors were calculated with software available free of charge available at https://preadmet.bmdrc.kr/adme/, https://www.cbligand.org/BBB/predictor.php.
Table 5. The docking and DFT (density functional theory) results obtained with respect to the studied ligands. The docking calculations were averaged over four teicoplanin structures available in the PDB (Protein Data Bank) records; the corresponding standard deviations are given.
Table 5. The docking and DFT (density functional theory) results obtained with respect to the studied ligands. The docking calculations were averaged over four teicoplanin structures available in the PDB (Protein Data Bank) records; the corresponding standard deviations are given.
LigandBinding Energy, Docking (kJ/mol)Binding Energy, DFT (kJ/mol)LigandBinding Energy, Docking (kJ/mol)Binding Energy, DFT (kJ/mol)
TRP-2−5.15 ± 0.53−10.22TRP-20−5.05 ± 0.29−11.47
TRP-4−4.95 ± 0.33−10.88TRP-22−5.30 ± 0.49−13.57
TRP-6−4.98 ± 0.11−12.62TRP-24−5.38 ± 0.29−12.10
TRP-8−4.90 ± 0.41−13.47TRP-26−5.05 ± 0.42−13.96
TPR-10−4.68 ± 0.23−14.93TPR-28−5.05 ± 0.38−12.56
TRP-14−4.88 ± 0.31−12.98TRP-30−4.80 ± 0.37−16.46
TRP-16−4.65 ± 0.21−14.59TRP-34−5.30 ± 0.44−15.96
Table 6. Metrics for classification of PR, RF and PLS-DA on test data.
Table 6. Metrics for classification of PR, RF and PLS-DA on test data.
Prediction: Active/NonactivePR—Regularized Logistic RegressionRF—Random ForestPLS-DA
Accuracy0.80.80.8
95% CI(0.2836, 0.9949)(0.2836, 0.9949)(0.2836, 0.9949)
No Information Rate0.60.60.6
p-Value (Acc > NIR)0.3370.3370.337
Kappa0.54550.54550.5455
McNemar’s Test p-Value1.0001.0001.000
Sensitivity1.001.001.00
Specificity0.500.500.50
Pos Pred Value0.750.750.75
Neg Pred Value1.001.001.00
Prevalence0.600.600.60
Detection Rate0.600.600.60
Detection Prevalence0.800.800.80
Balanced Accuracy0.750.750.75
Positive classactiveactiveactive

Share and Cite

MDPI and ACS Style

Flieger, J.; Orzeł, A.; Kowalska-Kępczyńska, A.; Pizoń, M.; Trębacz, H.; Majerek, D.; Plech, T.; Płaziński, W. Teicoplanin-Modified HPLC Column as a Source of Experimental Parameters for Prediction of the Anticonvulsant Activity of 1,2,4-Triazole-3-Thiones by the Regression Models. Materials 2020, 13, 2650. https://doi.org/10.3390/ma13112650

AMA Style

Flieger J, Orzeł A, Kowalska-Kępczyńska A, Pizoń M, Trębacz H, Majerek D, Plech T, Płaziński W. Teicoplanin-Modified HPLC Column as a Source of Experimental Parameters for Prediction of the Anticonvulsant Activity of 1,2,4-Triazole-3-Thiones by the Regression Models. Materials. 2020; 13(11):2650. https://doi.org/10.3390/ma13112650

Chicago/Turabian Style

Flieger, Jolanta, Anna Orzeł, Anna Kowalska-Kępczyńska, Magdalena Pizoń, Hanna Trębacz, Dariusz Majerek, Tomasz Plech, and Wojciech Płaziński. 2020. "Teicoplanin-Modified HPLC Column as a Source of Experimental Parameters for Prediction of the Anticonvulsant Activity of 1,2,4-Triazole-3-Thiones by the Regression Models" Materials 13, no. 11: 2650. https://doi.org/10.3390/ma13112650

APA Style

Flieger, J., Orzeł, A., Kowalska-Kępczyńska, A., Pizoń, M., Trębacz, H., Majerek, D., Plech, T., & Płaziński, W. (2020). Teicoplanin-Modified HPLC Column as a Source of Experimental Parameters for Prediction of the Anticonvulsant Activity of 1,2,4-Triazole-3-Thiones by the Regression Models. Materials, 13(11), 2650. https://doi.org/10.3390/ma13112650

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop