Next Article in Journal
Design, Synthesis, and Evaluation of Acetylcholinesterase and Butyrylcholinesterase Dual-Target Inhibitors against Alzheimer’s Diseases
Previous Article in Journal
Huntington’s Disease: A Review of the Known PET Imaging Biomarkers and Targeting Radiotracers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Predicting the Blood-Brain Barrier Permeability of New Drug-Like Compounds via HPLC with Various Stationary Phases

by
Małgorzata Janicka
1,
Małgorzata Sztanke
2,* and
Krzysztof Sztanke
3,*
1
Department of Physical Chemistry, Faculty of Chemistry, Institute of Chemical Science, Maria Curie-Skłodowska University, Maria Curie-Skłodowska Sq. 3, 20-031 Lublin, Poland
2
Chair and Department of Medical Chemistry, Medical University, 4A Chodźki Street, 20-093 Lublin, Poland
3
Laboratory of Bioorganic Synthesis and Analysis, Chair and Department of Medical Chemistry, Medical University, 4A Chodźki Street, 20-093 Lublin, Poland
*
Authors to whom correspondence should be addressed.
Molecules 2020, 25(3), 487; https://doi.org/10.3390/molecules25030487
Submission received: 23 November 2019 / Revised: 15 January 2020 / Accepted: 21 January 2020 / Published: 23 January 2020
(This article belongs to the Section Analytical Chemistry)

Abstract

:
The permeation of the blood-brain barrier is a very important consideration for new drug candidate molecules. In this research, the reversed-phase liquid chromatography with different columns (Purosphere RP-18e, IAM.PC.DD2 and Cosmosil Cholester) was used to predict the penetration of the blood-brain barrier by 65 newly-synthesized drug-like compounds. The linear free energy relationships (LFERs) model (log BB = c + eE + sS + aA + bB + vV) was established for a training set of 23 congeneric biologically active azole compounds with known experimental log BB (BB = Cblood/Cbrain) values (R2 = 0.9039). The reliability and predictive potency of the model were confirmed by leave-one-out cross validation as well as leave-50%-out cross validation. Multiple linear regression (MLR) was used to develop the quantitative structure-activity relationships (QSARs) to predict the log BB values of compounds that were tested, taking into account the chromatographic lipophilicity (log kw), polarizability and topological polar surface area. The excellent statistics of the developed MLR equations (R2 > 0.8 for all columns) showed that it is possible to use the HPLC technique and retention data to produce reliable blood-brain barrier permeability models and to predict the log BB values of our pharmaceutically important molecules.

1. Introduction

The biological activity of drugs depends primarily on their pharmacokinetics. The expected pharmacological effect of a given drug can be observed if the pharmacokinetic processes provide its high concentration within the range of the receptor. The amount of drug in tissue and the time that remains an effective concentration depend on the fundamental processes that make up the pharmacokinetic phase of the drug’s action, i.e., liberation (L), absorption (A), distribution (D), metabolism (M) and excretion (E). It is extremely difficult to predict the processes mentioned above because all of them are concentration-dependent and connected with the chemical structure of the agent. Since most drugs must pass through at least one cell membrane to provide the desired effect, for the rational design of drugs, it is vitally important to understand and to be able to predict the solute partitioning in the biomembranes. Drugs can cross membranes by passive or active transport [1,2,3]. While active transport is determined by compounds’ affinities for specific transporters and it uses energy, the most common mode for the passage of xenobiotics is passive transport, which depends on physicochemical properties of the compound, such as lipophilicity, size of the molecule, ionization state, as well as the diffusion coefficient through the membrane and the concentration gradient of the compound [4,5].
One of the most important properties of a potential drug is the ability of its molecule to penetrate the blood-brain barrier (BBB). Potential effective agents that are intended to interact with the central nervous system must be able to cross the BBB and satisfactory transport through the blood-brain barrier is an essential prerequisite for a potential drug to affect the central nervous system. However, the agents that act peripherally should not cross the BBB in order to avoid side effects. In both cases, the permeability of the BBB must be known and it should be evaluated at the earliest possible stage of testing. Doing so allows scientists to choose drug candidates that have more selective pharmacologic properties with fewer side effects and lower toxicities [6].
The common measure of the extent permeation of the blood-brain barrier is the ratio of the concentration of the drug in the brain (Cbrain) to the concentration of the drug in the blood (Cblood) or in the plasma (Cplasma), which is expressed as the log BB (BB = Cblood/Cbrain) [7,8,9]. Although measurement of the blood-brain barrier penetration in vivo is essential, the procedure is time-consuming, expensive and difficult. In addition, in recent years, an emphasis has been placed on modelling the log BB permeation to avoid unethical animal testing. Over the past three decades, various models have been proposed for predicting the BBB permeation and they have suggested different descriptors of the physicochemical properties of substances. Research shows that the penetration of the blood-brain barrier of the compound depends on its hydrogen bonding potential, lipophilicity and size. The BBB penetration is promoted by a weak potential for hydrogen bonding and high lipophilicity [4,5,7,8,9].
As noted, lipophilicity is one of the most important features affecting the biological activity. Octanol-water partition coefficients (logs P) are the most extensively used measure of lipophilicity in modelling the biological partition/distribution. This value can be determined by the classical shake-flask method, which is a time-consuming and tedious procedure. Contrarily, the liquid chromatography is a convenient, reliable and efficient method for assessing the partition parameters that describe the lipophilic properties of organic compounds. The background to this is that the same basic molecular interactions determine the behaviour of the solute in both biological and chromatographic systems. Moreover, there is an increasing evidence of the convenience for modelling pharmacokinetic processes chromatographically, especially by reversed-phase liquid chromatography using biomimetic stationary phases. The octadecylsilyl (ODS) stationary phase provides a fast approach, but immobilized artificial membranes (IAMs) are more similar to the membranes of eukaryotic cells and therefore better mimic biological systems [10,11,12]. Artificial membranes are more similar to biological systems because they anchor synthetic phosphatidylcholine analogues to silica [13,14,15,16,17,18]. Cholesterol is one of the major components of many eukaryotic membranes and it seems highly likely that cholesterol immobilized on silica would offer similar possibilities. Currently, the stationary phases with immobilized cholesterol are becoming more and more popular and therefore they are increasingly used to study biological properties of different organic compounds [19,20].
All pharmaceutically relevant compounds [19,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35] were resynthesized in our laboratory for the current research purposes and their structures, belonging to particular groups, are listed in Table 1. They have been shown to possess mainly promising anticancer [19,21,22,23,24,25,26,27,28,29,30,32], analgesic [19,21,31,33], antiviral and antihaemolytic [27] activities. Small molecules 16 (group I) and 6165 (group VII) are of particular importance as possible anticancer drug candidates for the treatment of human tumours of lung, cervix, breast and ovary [22,23,24]. In addition, the most promising structures 16 (i.e., showing the minimum embryotoxic concentration higher or comparable to that of aciclovir as well as protective effects on oxidatively-stressed erythrocytes) revealed significant anticancer activities in human tumour cells of pharynx and tongue. The majority of them proved to be more selective than a clinically useful anticancer agent—hydroxycarbamide. Besides, the compound 6 has been shown to possess the remarkable concentration-dependent potency against Herpes simplex virus type 1, while revealing a low toxicity to normal Vero cells and inhibiting the oxidatively-induced haemolysis of erythrocytes [27]. In turn, the confirmed remarkable antiproliferative effects of compounds 1217 (group III) may be of benefit in the treatment of human multiple myeloma cells that are susceptible and resistant to thalidomide as well as in human tumour cells of cervix and breast [25,26]. Molecules 1, 2, 4, 5, 6, 15, 19, 21, 22, 24, 28, 39, 63 and 64 have been reported as promising anticancer drug candidates, due to not only their proven significant antiproliferative activities in some human cancer cells but also less toxic effects for normal cells [22,23,24,25,26,27]. Furthermore, test compounds proved to be in vivo active when investigated in the central nervous system. Among analgesic active and relatively low toxic molecules (811, 32, 34, 39, 42, 48, 51 and 53), the structures 8, 42 and 51 have been shown to produce the strongest antinociceptive effect in mice [19,21,31,33].

2. Results

2.1. Chromatographic Results

Retention parameters reported as the log k values were calculated by the expression:
log k = log ( t R t 0 ) t 0
where tR and t0 are the retention times of the solute and a non-retained compound (citric acid), respectively. They were used to calculate the log kw values, i.e., logarithms of retention parameter in the buffer as the mobile phase. For this purpose the Soczewiński-Wachtmeister’s equation was used [36]:
log k = log kwsφ
where φ is the volume fraction of organic modifier in the mobile phase; k and kw are retention parameters corresponding to mixed effluent and buffer as the mobile phase, respectively. The slope s is characteristic of a given solute and chromatographic system. Strong linear relationships between log k and φ values were found for all the compounds in the range of effluent composition examined (R2 > 0.9). The log kw and s values obtained from particular chromatographic systems are presented in Table 2. The log kw values determined for the ODS (log kw, ODS), IAM (log kw, IAM) and Cholester (log kw, Cholester) columns were intercorrelated and the following linear relationships with a very good statistical quality were obtained:
log kw, ODS = 0.691(± 0.094) + 0.982(± 0.044) log kw, IAM
SD = 0.257; R2 = 0.8880; N = 65
log kw, ODS = 0.911(± 0.062) + 0.765(± 0.025) log kw, Cholester
SD = 0.192; R2 = 0.9371; N = 65
log kw, IAM = 0.356(± 0.073) + 0.721(± 0.029) log kw, Cholester
SD = 0.224; R2 = 0.9068; N = 65
Moreover highly significant linear relationships were obtained between log kw and s values (intercepts and slopes of Equation (2)):
for the ODS column:
log kw, ODS = −1.214(± 0.062) + 0.823(± 0.013) sODS
SD = 0.094; R2 = 0.9850; N = 65
for the IAM column:
log kw, IAM = −0.779(± 0.026) + 0.638(± 0.006) sIAM
SD = 0.052; R2 = 0.9950; N = 65
and for the Cholester column:
log kw, Cholester = −0.892(± 0.040) + 0. 798(± 0.010) sCholester
SD = 0.092; R2 = 0.9910; N = 65

2.2. Establishment of the LFER Model

The property of the substance can be predicted on the basis of the linear free-energy relationships (LFER), but to do so, the relationship between the chemical structure and the desired property should be identified [37]. A symbolic representation of LFERs model is the equation originally employed by Abraham et al. [38,39,40,41,42]:
SP = c + eE + sS + aA + bB + vV
Here SP is a set of solute properties in a given system, e.g., log BB values. The independent values are solute descriptors: E is an excess molar refraction, S is the dipolarity/polarizability, A and B are the hydrogen bond acidity (donating ability) and basicity (accepting ability), respectively and V is the solute McGowan volume (cm3∙mol−1/100). Coefficients c, e, s, a, b and v are characteristic for a given biphasic system and solute.
In this study, this equation was established for a training set of 23 azole compounds that were congeneric with those tested in our investigations. For these compounds, we obtained from the literature [43] the experimental log BB (BB = Cbrain/Cplasma) values for rats (Table 3). The actual values of log BB ranged from −0.82 to 0.58. The following MLR equation was obtained:
log BBexp. = 0.934(± 0.166) + 0.191(± 0.107) E − 0.605(± 0.134) S − 0.743(± 0.137) A
0.768(± 0.177) B + 0.545(± 0.104) V
N = 23; SD = 0.134; R2 = 0.9039; F = 32; p < 10−5
Irrelevant cross correlations between the descriptors were observed; the values of R2 between pairs of descriptors were E/S 0.66, E/A 0.01, E/B 0.25, E/V 0.60, S/A 0.06, S/B 0.69, S/V 0.72, A/B 0.20, A/V 0.05 and B/V 0.49. This is an important information indicating that there are no significant inter-correlations of structural descriptors, i.e., E, S, A, B and V. On this basis, we can conclude that non-physical factors do not affect the parameters of Equation (10). The reliability and predictive potency of the model expressed by Equation (10) were estimated by leave-one-out (LOO) cross validation and the parameters that were obtained are presented in Table 4. Figure 1 shows the standardized coefficients for particular descriptors and it confirms the well-known qualitative relationships: compound polarity, i.e., dipolarity/polarizability, hydrogen bond acidity and hydrogen bond basicity expressed by S, A and B values, respectively, decreases the BBB permeation, while the compound size measured by the McGovan volume (V) as well as the excess molar refraction E (in a minor degree) contribute to increase of log BB values [41]. The PLS response plot is presented in Figure 2, which shows the linear regression between the predicted (calculated response) and experimental log BB values of the 23 compounds from Table 3. The residual versus leverage plot in Figure 3 proves that the model that was obtained is valid for the domain in which it was developed [44]. The warning leverage limit (h*) was calculated according to:
h * = 3 m n
where m is the number of descriptors and n is the number of observations (compounds) in the dataset.
Equation (10) was used to calculate the log BB values for our 65 pharmaceutically relevant compounds (Table 5).

2.3. Establishment of QSARs Models

Efforts to predict the biological activity (including the BBB permeation) based on the properties of substances led to the development of Quantitative Structure Activity Relationships (QSARs) and Quantitative Retention Activity Relationships (QRARs) methodology [37,44,45,46,47,48]. Various models and approaches have been developed to predict the penetration of the blood-brain barrier based on various physicochemical properties of molecules, including the lipophilicity, molecular size, polarizability, polar surface area and the number of groups that can establish potential hydrogen bonds [49,50,51,52,53,54,55,56,57,58,59,60,61,62,63]. It is reasonable to assume that the combination of theoretical and experimental data increases the reliability of the anticipated transport of the drug across the blood-brain barrier [64,65,66]. The chromatographic retention parameter is one of the most popular experimental values used to characterize the properties (lipophilicity/hydrophobicity) of compounds used in QRAR and QSAR studies. The solute retention depends on the changes in free energy that are associated with the distribution of the solute between the mobile and stationary phases in a given chromatographic system. Thus, it is possible to use the values obtained on HPLC columns that imitate biomembranes for modelling the blood-brain barrier permeation.
In our investigations in which we modelled the blood-brain permeation of 65 biologically active molecules, the chromatographic lipophilicity (log kw) (Table 2), polarizability (α) and topological polar surface area (TPSA) (Table 5) were considered. The following MLR equations corresponding to various stationary phases were obtained:
for ODS column:
log BB = 0.587(± 0.127) + 0.011(± 0.020) log kw. ODS − 0.013(± 0.001) TPSA + 0.008(± 0.004) α
N = 65; SD = 0.069; R2 = 0.8474; F = 113; p < 10−6; VIF < 3.1
for IAM column:
log BB = 0.578(± 0.127) + 0.007(± 0.019) log kw, IAM − 0.013(± 0.001) TPSA + 0.009(± 0.004) α
N = 65; SD = 0.069; R2 = 0.8469; F = 113; p < 10−6; VIF < 2.7
and for Cholester column:
log BB = 0.595(± 0.139) + 0.008(± 0.017) log kw, Cholester − 0.013(± 0.001) TPSA +
0.009(± 0.004) α
N = 65; SD = 0.069; R2 = 0.8471; F = 113; p < 10−6; VIF < 3.7
The statistics of Equations (12)–(14) were very good, i.e., R2 > 0.8 for all chromatographic systems and all variance inflation factors (VIF < 5) indicated that the variables were correlated moderately. The reliability of the MLR models, as expressed by Equations (12)–(14), were estimated by leave-one-out as well as leave-50%-out cross validation (Table 4). In each case, the log kw and α values provided positive inputs to the log BB, while TPSA decreased the permeation of the blood-brain barrier (Figure 4).
The relationships showed that molecular polarizability and lipophilicity promote increases in the log BB, while polar surface area decreased the ability of a substance to cross the BBB. The correlations between the log BB values calculated according to the LFER model (Equation (10)) and the optimized QSARs models, i.e., by Equations (12)–(14) are presented in Figure 5, Figure 6 and Figure 7, which shows the response plots obtained by PLS for particular stationary phases.
To assess the significance of chromatographic parameters, the lipophilicity descriptors, log kw were compared with partition coefficients in the n-octanol-water system, i.e., the log P values. The relationships between log P values calculated from molecular structures of the tested compounds and obtained by use of Alog Ps algorithm [67,68] and their chromatographic factors (i.e., log kw, ODS, log kw, IAM and log kw, Cholester, respectively) were established. Good linear correlations between these descriptors (R > 0.8), that are observed in Figure 8, confirmed that chromatographic parameters can be used as lipophilicity descriptors in case of the studied compounds (Figure 8).
Moreover the following QSAR equation including Alog Ps, TPSA and α descriptors was established:
log BB = 0.627(± 0.115) + 0.077(± 0.028) Alog Ps − 0.012(± 0.001) TPSA + 0.001(± 0.004) α
N = 65; SD = 0.065; R2 = 0.8639; F = 130; p < 10−6; VIF < 4.3; Q2 = 0.8572; MSE =
0.004287; MSEcv* = 0.004218; MSEcv** = 0.004218; PRESS* = 0.2900; PRESS** = 0.2698.
The statistics of Equation (15) proved to be very good and similar to those obtained for Equations (12)–(14), which confirms their ability in predicting the blood brain barrier permeation.

3. Materials and Methods

3.1. Reagents

Acetonitrile (HPLC grade) was purchased from Merck (Lublin, Poland). Citric acid and Na2HPO4 (both pure) were supplied from POCh (Lublin, Poland). Distilled water was obtained from Direct-Q3 UV apparatus (Millipore, Warsaw, Poland).

3.2. Instrumental

Shimadzu Vp (Shimadzu, Izabelin, Poland) liquid chromatographic system equipped with LC 10AT pump, SPD 10A UV-Vis detector, SCL 10A system controller, CTO-10 AS chromatographic oven and Rheodyne injector valve with a 20 µL loop was applied in HPLC measurements. Three stationary phases were employed:
  • Purosphere RP-18e (ODS), 125 × 4 mm i.d., 5 µm (Merck);
  • IAM.PC.DD2 100 × 4.6 mm i.d., 10 µm (Regis Chemicals Company, Morton Grove, IL, USA);
  • Cosmosil Cholester, 75 × 2 mm i.d., 2.5 µm (Genore, Warsaw, Poland).

3.3. Chromatographic Conditions and Test Substances

As mobile phases buffer-acetonitrile mixtures were used. The buffer was prepared from 0.01 mol L−1 solutions of Na2HPO4 and citric acid and the pH 7.4 value was fixed before mixing with organic modifier. With the ODS column acetonitrile concentration in the effluent, expressed as a volume fraction (φ, v/v), was changed in the range 0.3–0.6, with the constant step of 0.1. The flow rate was 1 mL min−1. With the IAM column acetonitrile concentration was changed in the range 0.2–0.5, also with the constant step of 0.1 and the flow rate was 1.3 mL min−1. With the Cosmosil Cholester column acetonitrile concentration was changed in the range 0.4–0.6, with the constant step of 0.05 and the flow rate was 0.4 mL min−1. Samples of test compounds were dissolved in acetonitrile—c.a. 0.005 mg mL−1. All the compounds proved to be in the neutral form in solution under experimental conditions and had values of peak asymmetry factor in the acceptable range. They were detected under UV light at 210 and 254 nm. All measurements were carried out at 25 °C. The dead time values were measured from non-retained compound (citric acid) peaks. All reported log kw values are the average of at least three independent measurements. The extrapolated retention coefficients (logs kw) achieved by HPLC on ODS, IAM and Cholester stationary phases were used for modelling the log BB permeation of 65 drug-like compounds employed as a whole test set (Table 1).

3.4. In silico Calculations

Molecular descriptors (E, S, A, B and V), molecular weight (MW), topological polar surface area (TPSA) and polarizability (α) of compounds were evaluated by ACD/Percepta software (Łodź, Poland).

3.5. Statistical Analysis

Linear regression, multiple linear regression (MLR), partial last square (PLS), leave-one-out (LOO) and leave-50%-out cross validation were performed using Minitab 16 software (Minitab Inc., State College, PA, USA).

4. Discussion

The usefulness of HPLC with three different reverse-phase columns (including two imitating biosystems) for predicting the blood-brain barrier (BBB) permeability of 65 structurally related drug-like molecules was highlighted in our investigations. QSAR models predicting the BBB permeation were built on the basis of experimentally accessible log kw values of the test molecules together with their important in silico molecular descriptors. The obtained results confirmed that all three stationary phases, i.e., octadecylsilane, immobilized artificial membrane and cholesterol immobilized on silica gel analogously described the lipophilic properties of the studied solutes (Equations (3)–(5)).
In our studies the extrapolated retention parameters (logs kw) were used as they are preferred in QSARs instead of the isocratic log k values and usually employed in correlation studies with in silico molecular descriptors and drug-likeness properties in case of pharmaceutics and drug-like molecules. It should be noted that chromatographically derived retention parameters are very useful descriptors in QSAR modelling as the partitioning process between the stationary and mobile phase of a solute investigated mimics a membrane penetration process of a pharmaceutic or potential drug candidate [19,21,33,37].
Highly significant linear relationships were obtained between intercepts and slopes of Equation (2), i.e., the log kw and s values, for particular reversed-phase stationary phases. These correlations confirmed not only the congenereity between compounds that were investigated, but also suggested that the log kw and s values may be considered as alternative lipophilicity descriptors [19,21,37] in case of structurally related small molecules that were tested. In our studies compounds bearing more hydrophobic substituents revealed the greater s values. This observation is consistent with other research findings showing that more hydrophobic molecules reveal greater slopes [21]. According to the background retention theory the s values are related to the solute/mobile phase and the solvent/stationary phase net interactions [21,36,37].
The obtained results showed that the chromatographic retention parameters obtained using three stationary phases recruited as well as important in silico molecular descriptors can be recommended to derive the reliable QSAR models for predicting the blood-brain barrier permeability in case of our structurally related small molecules, considered as a test set of potential drug candidates.
The calculated log P values were obtained by using Alog Ps algorithm and compared with the experimental log kw, ODS, log kw, IAM and log kw, Cholester values, respectively. This estimation was essential to check the validity of the obtained results through correlation with the log BB values.

5. Conclusions

Experimental literature data of the log BB (for rats) for a training set of 23 biologically active compounds (including drugs) were correlated against five solute descriptors (E, S, A, B and V) and the established equation (Equation (10)) was validated. The standardized coefficients obtained for particular descriptors confirmed that the compound polarity, i.e., dipolarity/polarizability (S), hydrogen bond acidity (A) and hydrogen bond basicity (B) decreases the blood-brain barrier (BBB) permeation while the compound size measured by the McGovan volume (V) and the excess molar refraction (E) contribute to increase in the log BB values. The Equation (10) was used to calculate the log BB values for 65 newly synthesized compounds being considered as potential drugs.
The blood-brain barrier permeability of the compounds that were tested was modelled by three descriptors, i.e., the chromatographic lipophilicity (log kw), polarizability (α) and topological polar surface area (TPSA). Using a simple statistical model (i.e., MLR), three structure-activity relationships were obtained for each chromatographic system on endcapped octadecylsilyl, immobilized artificial membrane and cholesteryl stationary phases (Equations (12)–(14)). The log BB values calculated according to Equation (10) and predicted from Equations (12)–(14) were compared and highly significant relationships were obtained between them. The relationships were confirmed by leave-one-out as well as leave-50%-out cross validation, implying that the models are robust and reliable. The results that were obtained showed that, in the case of the compounds that were studied (65 weak organic bases), each of the stationary phases used in the chromatographic measurements was equally useful. From practical and economic perspectives, the Cholester microcolumn is recommended because it allows the acquisition of more data in a shorter time with lower costs.
We showed that it is possible to use the HPLC technique to build a reliable model for predicting which our organic compounds (congeneric in structure to the training set of azole molecules) can penetrate the blood-brain barrier in rats. The investigations highlighted the key role of chromatographic techniques and QSARs methods in reducing unethical animal testing.
The presented results will be particularly useful in further more extensive in vivo research, which is planned to be carried out on our small molecules considered as potential drugs.

Author Contributions

Conceptualization, M.J.; Methodology, M.J.; Investigation, M.J.; Formal analysis, M.J.; Funding acquisition, M.S., K.S.; Resources, M.J., M.S., K.S.; Writing—Original Draft, M.J., M.S., K.S.; Writing—Review & Editing, M.J., M.S., K.S. 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. Sugano, K.; Kansy, M.; Artursson, P.; Avdeef, A.; Bendels, S.; Di, L.; Ecker, G.F.; Faller, B.; Fischer, H.; Gerebtzoff, G.; et al. Coexistence of passive and carrier-mediated processes in drug transport. Nat. Rev. Drug Discov. 2010, 9, 597–614. [Google Scholar] [CrossRef]
  2. Sugiyama, Y.; Kusuhara, H.; Suzuki, H. Kinetic and biochemical analysis of carrier-mediated efflux of drugs through the blood-brain and blood-cerebrospinal fluid barriers: Importance in the drug delivery to the brain. J. Cont. Rel. 1999, 62, 179–186. [Google Scholar] [CrossRef]
  3. Seddon, A.M.; Casey, D.; Law, R.V.; Gee, A.; Templera, R.H.; Cesab, O. Drug interactions with lipid membranes. Chem. Soc. Rev. 2009, 38, 2509–2519. [Google Scholar] [CrossRef] [PubMed]
  4. Wolak, D.J.; Thorne, R.G. Diffusion of macromolecules in the brain: Implications for drug delivery. Mol. Pharm. 2013, 10, 1492–1504. [Google Scholar] [CrossRef] [PubMed]
  5. Banks, W.A. Characteristics of compounds that cross the blood-brain barrier. BMC Neurol. 2009, 9 (Suppl. 1), S3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Konovalov, D.A.; Coomans, D.; Deconinck, E.; Heyden, Y.V. Benchmarking of QSAR Models for Blood-Brain Barrier Permeation. J. Chem. Inf. Model. 2007, 47, 1648–1656. [Google Scholar] [CrossRef] [PubMed]
  7. Goodwin, J.T.; Clark, D.E. In Silico Predictions of Blood-Brain Barrier Penetration: Considerations to “Keep in Mind”. J. Pharm. Exp. Therap. 2005, 315, 477–483. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Abbott, N.J. ADME-Tox Approaches. In Comprehensive Medicinal Chemistry II, 2nd ed.; Taylor, J.B., Triggle, D.J., Eds.; Elsevier LTD: Amsterdam, The Netherlands, 2007; pp. 301–320. [Google Scholar]
  9. Bickel, U. How to Measure Drug Transport across the Blood-Brain Barrier. NeuroRx 2005, 2, 15–26. [Google Scholar] [CrossRef]
  10. Liu, X.; Testa, B.; Fahr, A. Lipophilicity and its relationship with passive drug permeation. Pharm. Res. 2011, 28, 962–977. [Google Scholar] [CrossRef]
  11. Grumetto, L.; Russo, G.; Barbato, F. Immobilized Artificial Membrane HPLC Derived Parameters vs PAMPA-BBB Data in Estimating in Situ Measured Blood-Brain Barrier Permeation of Drugs. Mol. Pharm. 2016, 13, 2808–2816. [Google Scholar] [CrossRef]
  12. Russo, G.; Grumetto, L.; Szucs, R.; Barbato, F.; Lynen, F. Screening therapeutics according to their uptake across the blood-brain barrier: A high throughput method based on immobilized artificial membrane liquid chromatography-diode-array-detection coupled to electrospray-time-of-flight mass spectrometry. Eur. J. Pharm. Biopharm. 2018, 127, 72–84. [Google Scholar] [CrossRef] [PubMed]
  13. 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–128. [Google Scholar] [CrossRef]
  14. Yang, C.Y.; Cai, S.J.; Liu, H.; Pidgeon, C. Immobilized Artificial Membranes – screens for drug membrane interactions. Adv. Drug Del. Rev. 1996, 23, 229–256. [Google Scholar] [CrossRef]
  15. Taillardat-Bertschinger, A.; Galland, A.; Carrupt, P.-A.; Testa, B. Immobilized artificial membrane liquid chromatography: Proposed guidelines for technical optimization of retention measurements. J. Chromatogr. A 2002, 953, 39–53. [Google Scholar] [CrossRef]
  16. Verzele, D.; Lynen, F.; De Vrieze, M.; Wright, A.G.; Hanna-Brown, M.; Sandra, P. Development of the first sphingomyelin biomimetic stationary phase for immobilized artificial membrane (IAM) chromatography. Chem. Commun. 2012, 48, 1162–1164. [Google Scholar] [CrossRef]
  17. De Vrieze, M.; Lynen, F.; Chen, K.; Szucs, R.; Sandra, P. Predicting drug penetration across the blood-brain barrier: Comparison of micellar liquid chromatography and immobilized artificial membrane liquid chromatography. Anal. Bioanal. Chem. 2013, 405, 6029–6041. [Google Scholar] [CrossRef]
  18. De Vrieze, M.; Verzele, D.; Szucs, R.; Sandra, P.; Lynen, F. Evaluation of sphingomyelin, cholester and phosphatidylcholine-based immobilized artificial membrane liquid chromatography to predict drug penetration across the blood-brain barrier. Anal. Bioanal. Chem. 2014, 406, 6179–6188. [Google Scholar] [CrossRef]
  19. Janicka, M.; Sztanke, M.; Sztanke, K. Reversed-phase liquid chromatography with octadecylsilyl, immobilized artificial membrane and cholesterol columns in correlation studies with in silico biological descriptors of newly synthesized antiproliferative and analgesic active compounds. J. Chromatogr. A 2013, 1318, 92–101. [Google Scholar] [CrossRef]
  20. Janicka, M.; Pachuta-Stec, A. Retention-property relationships of 1,2,4-triazoles by micellar and reversed-phase liquid chromatography. J. Sep. Sci. 2014, 31, 1419–1428. [Google Scholar] [CrossRef]
  21. Sztanke, M.; Tuzimski, T.; Janicka, M.; Sztanke, K. Structure-retention behaviour of biologically active fused 1,2,4-triazinones – Comparison with in silico molecular properties. Eur. J. Pharm. Sci. 2015, 68, 114–126. [Google Scholar] [CrossRef]
  22. Sztanke, M.; Rzymowska, J.; Janicka, M.; Sztanke, K. Synthesis, structure confirmation, identification of in vitro antiproliferative activities and correlation of determined lipophilicity parameters with in silico bioactivity descriptors of two novel classes of fused azaisocytosine-like congeners. Arab. J. Chem. 2016. [Google Scholar] [CrossRef]
  23. Sztanke, M.; Sztanke, K.; Medical University of Lublin. 3-(2-Phenylethyl)-8-aryl-7,8-dihydroimidazo[2,1-c][1,2,4]triazin-4(6H)-ones, Method for Obtaining them and Medical Application. Polish Patent 224678, 31 January 2017. [Google Scholar]
  24. Sztanke, M.; Sztanke, K.; Medical University of Lublin. 3-Ethyl-8-aryl-7,8-dihydroimidazo[2,1-c][1,2,4]triazin-4(6H)-ones, Method for Obtaining them and Medical Application. Polish Patent 224679, 31 January 2017. [Google Scholar]
  25. Sztanke, K.; Sztanke, M.; Medical University of Lublin. Ethyl 2-(4-oxo-4,6,7,8-tetrahydroimidazo[2,1-c][1,2,4]triazin-3-yl)acetates, Method for Obtaining them and Medical Application. Polish Patent 219424, 30 April 2015. [Google Scholar]
  26. Sztanke, M.; Rzymowska, J.; Sztanke, K. Synthesis, structure elucidation and in vitro anticancer activities of novel derivatives of diethyl (2E)-2-[(2E)-(1-arylimidazolidin-2-ylidene)hydrazono]succinate and ethyl (4-oxo-8-aryl-4,6,7,8-tetrahydroimidazo[2,1-c][1,2,4]triazin-3-yl)acetate. Bioorg. Med. Chem. 2013, 21, 7465–7480. [Google Scholar] [CrossRef] [PubMed]
  27. Sztanke, M.; Sztanke, K.; Rajtar, B.; Świątek, Ł.; Boguszewska, A.; Polz-Dacewicz, M. The influence of some promising fused azaisocytosine-containing congeners on zebrafish (Danio rerio) embryos/larvae and their antihaemolytic, antitumour and antiviral activities. Eur. J. Pharm. Sci. 2019, 132, 34–43. [Google Scholar] [CrossRef] [PubMed]
  28. Sztanke, K.; Tuzimski, T.; Sztanke, M.; Rzymowska, J.; Pasternak, K. Synthesis, structure elucidation, determination of the lipophilicity and identification of antitumour activities in vitro of novel 3-(2-furanyl)-8-aryl-7,8-dihydroimidazo[2,1-c][1,2,4]triazin-4(6H)-ones with a low cytotoxicity towards normal human skin fibroblast cells. Bioorg. Med. Chem. 2011, 19, 5103–5116. [Google Scholar]
  29. Sztanke, K.; Sztanke, M.; Pasternak, K.; Medical University of Lublin. 3-(2-Furanyl)-7,8- dihydroimidazo[2,1-c][1,2,4]triazin-4(6H)-ones Substituted with Mono- or Dichlorophenyl and Process for the Preparation Thereof. Polish Patent PL 212442, 31 October 2012. [Google Scholar]
  30. Sztanke, K.; Sztanke, M.; Pasternak, K.; Medical University of Lublin. Derivatives of 3-(2-furanyl)-7,8-dihydroimidazo[2,1-c][1,2,4]triazin-4(6H)-one Substituted with phenyl, Alkylphenyl, Alkoxyphenyl and Process for the Preparation Thereof. Polish Patent PL 212447, 31 October 2012. [Google Scholar]
  31. Sztanke, K.; Medical University of Lublin. New 8-aryl-3-phenyl-6,7-dihydro-4H-imidazo[2,1-c][1,2,4]triazine-4-ones and Methods for their Manufacture. Polish Patent PL 199750, 31 October 2008. [Google Scholar]
  32. Sztanke, K.; Pasternak, K.; Sztanke, M.; Kandefer-Szerszeń, M.; Kozioł, A.E.; Dybała, I. Crystal structure, antitumour and antimetastatic activities of disubstituted fused 1,2,4-triazinones. Bioorg. Med. Chem. Lett. 2009, 19, 5095–5100. [Google Scholar] [CrossRef]
  33. Tuzimski, T.; Sztanke, K. Retention data for some carbonyl derivatives of imidazo[2,1-c][1,2,4]triazine in reversed-phase systems in TLC and HPLC and their use for determination of lipophilicity. Part 1. Lipophilicity of 8-aryl-3-phenyl-6,7-dihydro-4H-imidazo[2,1-c][1,2,4]triazin-4-ones. J. Planar Chromatogr. 2005, 18, 274–281. [Google Scholar]
  34. Łyszczek, R.; Bartyzel, A.; Głuchowska, H.; Mazur, L.; Sztanke, M.; Sztanke, K. Thermal investigations of biologically important fused azaisocytosine-containing congeners and the crystal structure of one representative. J. Anal. Appl. Pyrolysis 2018, 135, 141–151. [Google Scholar] [CrossRef]
  35. Bartyzel, A.; Sztanke, M.; Sztanke, K. Thermal behaviour of antiproliferative active 3-(2-furanyl)-8-aryl-7,8-dihydroimidazo[2,1-c][1,2,4]triazin-4(6H)-ones. J. Therm. Anal. Calorim. 2017, 130, 1541–1551. [Google Scholar] [CrossRef] [Green Version]
  36. Soczewiński, E.; Wachtmeister, C.A. The relation between the composition of certain ternary two-phase solvent systems and RM values. J. Chromatogr. 1962, 7, 311–320. [Google Scholar] [CrossRef]
  37. Kaliszan, R. Quantitative Structure-Chromatographic Retention Relationships; John Wiley & Sons: New York, NY, USA, 1987. [Google Scholar]
  38. Abraham, M.H.; Chadha, H.S.; Mitchell, R.C. Hydrogen bonding. 33. Factors that influence the distribution of solutes between blood and brain. J. Pharm. Sci. 1994, 83, 1257–1268. [Google Scholar] [CrossRef]
  39. Abraham, M.H. Scales of Solute Hydrogen-bonding: Their Construction and Application to Physicochemical and Biochemical Processes. Chem. Soc. Rev. 1993, 22, 73–83. [Google Scholar] [CrossRef]
  40. Abraham, M.H.; Ibrahim, A.; Zissimos, A.M. Determination of sets of solute descriptors from chromatographic measurements. J. Chromatogr. A 2004, 1037, 29–47. [Google Scholar] [CrossRef] [PubMed]
  41. Abraham, M.H. The factors that influence permeation across the blood-brain barrier. Eur. J. Med. Chem. 2004, 39, 235–240. [Google Scholar] [CrossRef] [PubMed]
  42. Vitha, M.; Carr, P.W. The chemical interpretation and practice of linear solvation energy relationships in chromatography. J. Chromatogr. A 2006, 1126, 143–194. [Google Scholar] [CrossRef] [PubMed]
  43. Mente, S.R.; Lombardo, F. A recursive-partitioning model for blood-brain barrier permeation. J. Comp. Aided Mol. Design. 2005, 19, 465–481. [Google Scholar] [CrossRef]
  44. Bujak, R.; Struck-Lewicka, W.; Kaliszan, M.; Kaliszan, R.; Markuszewski, M.J. Blood-brain barrier permeability mechanisms in view of quantitative structure-activity relationships (QSAR). J. Pharm. Biomed. Anal. 2015, 108, 29–37. [Google Scholar] [CrossRef]
  45. Carpenter, T.S.; Kirshner, E.D.A.; Lau, 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–641. [Google Scholar] [CrossRef] [Green Version]
  46. Chen, Y.; Zhu, Q.-J.; Pan, J.; Wu, X.-P. A prediction model for blood-brain barrier permeation and analysis on its parameter biologically. Comp. Meth. Prog. Biomed. 2009, 95, 280–287. [Google Scholar] [CrossRef] [PubMed]
  47. Ooms, F.; Weber, T.; Carrupt, P.A.; Testa, B. A simple model to predict blood-brain barrier permeation from 3D molecular fields. Biochim. Biophys. Acta 2002, 1587, 118–125. [Google Scholar] [CrossRef]
  48. Sagandykova, G.N.; Pomastowski, P.P.; Kaliszan, R.; Buszewski, B. Modern analytical methods for consideration of natural biological activity. Trends Anal. Chem. 2018, 109, 198–213. [Google Scholar] [CrossRef]
  49. Kaliszan, R.; Markuszewski, M. Brain/blood distribution described by a combination of partition coefficient and molecular mass. Int. J. Pharm. 1996, 145, 9–16. [Google Scholar] [CrossRef]
  50. Österberg, T.; Norinder, U. Prediction of polar surface area and drug transport processes using simple parameters and PLS statistics. J. Chem. Inf. Comput. Sci. 2000, 40, 1408–1411. [Google Scholar] [CrossRef] [PubMed]
  51. Clark, D.E. Rapid calculation of polar molecular surface area and its application to the prediction of transport phenomena. 2. Prediction of blood-brain barrier penetration. J. Pharm. Sci. 1999, 88, 815–821. [Google Scholar] [CrossRef] [PubMed]
  52. Norinder, U.; Haeberlein, M. Computational approaches to the prediction of the blood-brain distribution. Adv. Drug Deliv. Rev. 2002, 54, 291–313. [Google Scholar] [CrossRef]
  53. Luco, J.M. Prediction of the brain-blood distribution of a large set of drugs from structurally derived descriptors using partial least-squares (PLS) modeling. J. Chem. Inf. Comput. Sci. 1999, 39, 396–404. [Google Scholar] [CrossRef]
  54. Feher, M.; Sourial, E.; Schmidt, J.M. A simple model for the prediction of blood-brain partitioning. Int. J. Pharm. 2000, 201, 239–247. [Google Scholar] [CrossRef]
  55. Fu, X.-C.; Song, Z.-F.; Fu, C.-Y.; Liang, W.-Q. A simple predictive model for blood-brain barrier penetration. Pharmazie 2005, 60, 354–358. [Google Scholar]
  56. Fu, X.-C.; Wang, G.-P.; Shan, H.-L.; Liang, W.-Q.; Gao, J.-Q. Predicting blood-brain barrier penetration from molecular weight and number of polar atoms. Eur. J. Pharm. Biopharm. 2008, 70, 462–466. [Google Scholar] [CrossRef]
  57. Young, R.C.; Mitchell, R.C.; Brown, T.H.; Ganellin, C.R.; Griffiths, R.; Jones, M.; Rana, K.K.; Saunders, D.; Smith, I.R.; Sore, N.E.; et al. Development of a new physicochemical model for brain penetration and its application to the design of centrally acting H2 receptor histamine antagonists. J. Med. Chem. 1988, 31, 656–671. [Google Scholar] [CrossRef]
  58. van de Waterbeemd, H.; Kansy, M. Hydrogen-bonding capacity and brain penetration. Chimia 1992, 46, 299–303. [Google Scholar]
  59. Lombardo, F.; Blacke, J.F.; Curatolo, W.J. Computation of brain-blood partitioning of organic solutes via free energy calculations. J. Med. Chem. 1996, 39, 4750–4755. [Google Scholar] [CrossRef] [PubMed]
  60. Kelder, J.; Grootenhuis, P.D.J.; Bayada, D.M.; Delbressine, L.P.C.; Ploemen, J.P. Polar molecular surface as a dominating determinant for oral absorption and brain penetration of drugs. Pharm. Res. 1999, 16, 1514–1519. [Google Scholar] [CrossRef] [PubMed]
  61. Ertl, P.; Rohde, B.; Selzer, P. Fast calculation of molecular polar surface area as a sum of fragment-based contributions and its application to the prediction of drug transport properties. J. Med. Chem. 2000, 43, 3714–3717. [Google Scholar] [CrossRef] [PubMed]
  62. Vilar, S.; Chakrabarti, M.; Constanzi, S. Prediction of passive blood-brain partitioning: Straightforward and effective classification models based on in silico derived physicochemical descriptors. J. Mol. Graph. Model. 2010, 28, 899–903. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Raevsky, O.A.; Solodova, S.L.; Raevskaya, O.E.; Mannhold, R. Molecular biology problems of drug design and mechanism of drug action. Quantitative interactions between the structures of organic compounds and their abilities to penetrate the blood-brain barrier. Pharm. Chem. J. 2012, 46, 133–138. [Google Scholar] [CrossRef]
  64. Ruiz-Ángel, M.J.; Garcia-Álvarez-Coque, M.C.; Berthod, A. New insights and recent developments in micellar liquid chromatography (review). Sep. Purif. Rev. 2009, 38, 45–96. [Google Scholar] [CrossRef]
  65. Mensch, J.; Oyarzabal, J.; Mackie, C.; Augustijns, P. In vitro, in vitro and in silico methods for small molecule transfer across the BBB. J. Pharm. Sci. 2009, 98, 4429–4468. [Google Scholar] [CrossRef]
  66. Naik, P.; Cucullo, L. In vitro blood-brain barrier models: Current and perspective technologies. J. Pharm. Sci. 2012, 101, 1337–1354. [Google Scholar] [CrossRef] [Green Version]
  67. Tetko, I.V.; Tanchuk, V.Y. Application of associative neural networks for prediction of lipophilicity in ALOGPS 2.1 program. J. Chem. Inf. Comput. Sci. 2002, 42, 1136–1145. [Google Scholar] [CrossRef]
  68. Tetko, I.V.; Tanchuk, V.Y.; Villa, A.E. Prediction of n-octanol/water partition coefficients from PHYSPROP database using artificial neural networks and E-state indices. J. Chem. Inf. Comput. Sci. 2001, 41, 1407–1421. [Google Scholar] [CrossRef]
Sample Availability: Samples of the compounds 165 are available from the authors.
Figure 1. PLS standardized coefficient plot (Equation (10)).
Figure 1. PLS standardized coefficient plot (Equation (10)).
Molecules 25 00487 g001
Figure 2. PLS response plot (Equation (10)).
Figure 2. PLS response plot (Equation (10)).
Molecules 25 00487 g002
Figure 3. PLS residual vs. leverage (Equation (10)).
Figure 3. PLS residual vs. leverage (Equation (10)).
Molecules 25 00487 g003
Figure 4. PLS standardized coefficient plots obtained for Equation (12) (A), Equation (13) (B) and Equation (14) (C).
Figure 4. PLS standardized coefficient plots obtained for Equation (12) (A), Equation (13) (B) and Equation (14) (C).
Molecules 25 00487 g004aMolecules 25 00487 g004b
Figure 5. PLS response plots obtained for Equation (12).
Figure 5. PLS response plots obtained for Equation (12).
Molecules 25 00487 g005aMolecules 25 00487 g005b
Figure 6. PLS response plots obtained for Equation (13).
Figure 6. PLS response plots obtained for Equation (13).
Molecules 25 00487 g006aMolecules 25 00487 g006b
Figure 7. PLS response plots obtained for Equation (14).
Figure 7. PLS response plots obtained for Equation (14).
Molecules 25 00487 g007aMolecules 25 00487 g007b
Figure 8. Alog Ps vs. log kw relationships obtained for ODS, IAM and Cholester columns.
Figure 8. Alog Ps vs. log kw relationships obtained for ODS, IAM and Cholester columns.
Molecules 25 00487 g008
Table 1. Compounds tested.
Table 1. Compounds tested.
GroupCompound No.; R, R’Chemical NameReferences
Group I
Molecules 25 00487 i001
1: R = H
2: R = 4-CH3
3: R = 2-Cl
4: R = 3-Cl
5: R = 4-Cl
6: R = 3,4-Cl2
8-(R-phenyl)-3-ethyl-7,8-dihydroimidazo[2,1-c][1,2,4]triazin-4(6H)-ones[22,24,27]
Group II
Molecules 25 00487 i002
7: R = H
8: R = 4-CH3
9: R = 4-OCH3
10: R = 4-OC2H5
11: R = 4-Cl
Methyl [4-oxo-8-(R-phenyl)-4,6,7,8-tetrahydroimidazo[2,1-c][1,2,4]triazin-3-yl]acetates[21]
Group III
Molecules 25 00487 i003
12: R = H
13: R = 4-CH3
14: R = 4-OCH3
15: R = 3-Cl
16: R = 4-Cl
17: R = 3,4-Cl2
Ethyl [4-oxo-8-(R-phenyl)-4,6,7,8-tetrahydroimidazo[2,1-c][1,2,4]triazin-3-yl]acetates[25,26]
Group IV
Molecules 25 00487 i004
18: R = H
19: R = 2-CH3
20: R = 4-CH3
21: R = 2,3(-CH3)2
22: R = 2-OCH3
23: R = 4-OCH3
24: R = 2-Cl
25: R = 3-Cl
26: R = 4-Cl
27: R = 3,4-Cl2
28: R = 2,6-Cl2
8-(R-phenyl)-3-(2-furanyl)-7,8-dihydroimidazo[2,1-c][1,2,4]triazin-4(6H)-ones[19,28,29,30]
Group V
Molecules 25 00487 i005
29: R = H
30: R = 2-CH3
31: R = 3-CH3
32: R = 4-CH3
33: R = 2-OCH3
34: R = 4-OCH3
35: R = 4-OC2H5
36: R = 2,3(-CH3)2
37: R = 2-Cl
38: R = 3-Cl
39: R = 4-Cl
40: R = 3,4-Cl2
8-(R-phenyl)-3-phenyl-7,8-dihydroimidazo[2,1-c][1,2,4]triazin-4(6H)-ones[19,31,32,33]
Group VI
Molecules 25 00487 i006
41: R = H; R’ = H
42: R = H; R’ = 2-Cl
43: R = H; R’ = 3-Cl
44: R = H; R’ = 4-Cl
45: R = 4-CH3; R’ = H
46: R = 4-CH3; R’ = 4-CH3
47: R = 4-CH3; R’ = 3-CH3
48: R = 4-CH3; R’ = 2-Cl
49: R = 4-CH3; R’ = 3-Cl
50: R = 4-CH3; R’ = 4-Cl
51: R = 4-OC2H5; R’ = H
52: R = 4-OC2H5; R’ = 4-CH3
53: R = 4-OC2H5; R’ = 2-Cl
54: R = 4-OC2H5; R’ = 3-Cl
55: R = 4-OC2H5; R’ = 4-Cl
56: R = 2-CH3; R’ = 2-Cl
57: R = 4-Cl; R’ = H
58: R = 4-Cl; R’ = 2-Cl
59: R = 4-Cl; R’ = 3-Cl
60: R = 4-Cl; R’ = 4-Cl
8-(R-phenyl)-3-benzyl/3-(R’-benzyl)-7,8-dihydroimidazo[2,1-c][1,2,4]triazin-4(6H)-ones[19,21]
Group VII
Molecules 25 00487 i007
61: R = H
62: R = 4-CH3
63: R = 2-Cl
64: R = 4-Cl
65: R = 3,4-Cl2
8-(R-phenyl)-3-(2-phenylethyl)-7,8-dihydroimidazo[2,1-c][1,2,4]triazin-4(6H)-ones[22,23]
Table 2. Parameters of the Soczewiński-Wachtmeister’s equation obtained for various columns.
Table 2. Parameters of the Soczewiński-Wachtmeister’s equation obtained for various columns.
Compound TestedODSIAMCholester
log kw, ODSsODSlog kw, IAMsIAMlog kw, CholestersCholester
10.972.700.552.100.201.25
21.303.121.002.880.571.84
31.223.020.802.510.461.70
41.333.150.842.600.621.82
51.783.561.243.251.182.61
62.534.491.854.112.053.61
71.283.090.812.380.581.92
81.903.821.333.331.292.75
91.303.020.812.450.571.82
101.953.911.423.391.362.81
112.364.441.824.391.823.55
121.753.611.213.111.112.55
132.374.361.703.851.923.58
141.373.150.912.580.651.92
152.854.902.054.452.414.15
162.804.882.114.512.374.19
173.165.332.425.012.774.61
182.124.061.293.231.433.12
191.743.480.962.721.202.78
202.434.391.693.902.033.68
212.264.201.423.511.663.33
221.673.481.113.021.222.68
231.913.821.193.131.533.05
242.254.221.273.251.342.95
252.754.852.084.552.364.05
262.654.592.024.442.294.15
273.465.683.236.302.924.85
282.534.451.703.902.053.66
292.494.511.914.222.344.02
302.124.011.533.521.883.50
312.935.022.264.682.774.42
322.935.062.224.812.864.65
331.983.991.473.621.692.94
342.594.621.713.912.273.80
353.095.332.144.482.804.61
362.744.851.854.052.343.60
372.434.601.733.841.963.65
383.295.702.645.413.205.20
393.415.202.575.323.195.30
403.485.903.266.403.866.02
412.905.023.266.452.233.75
423.285.102.344.782.834.71
433.595.852.485.182.904.80
443.605.622.394.813.014.70
453.355.552.154.692.744.59
463.045.302.515.223.155.05
473.095.242.595.143.095.06
483.125.112.655.213.244.90
493.315.412.785.503.405.33
503.225.552.815.623.395.34
513.285.612.114.512.684.55
523.145.362.435.023.075.02
532.935.022.515.173.215.15
543.355.552.735.613.315.22
553.325.662.725.443.365.35
563.014.952.024.312.344.05
573.115.052.525.062.984.85
583.555.653.226.323.745.62
593.615.853.156.013.675.75
604.026.453.176.183.685.82
612.804.922.074.452.334.02
623.235.442.415.012.844.70
632.464.441.814.081.943.66
643.565.812.685.443.225.12
654.296.663.406.554.126.20
Table 3. Experimental log BB values [43] and molecular descriptors (A, B, S, E, V) calculated using ACD Percepta for a training set of compounds.
Table 3. Experimental log BB values [43] and molecular descriptors (A, B, S, E, V) calculated using ACD Percepta for a training set of compounds.
No.CAS #ABSEVlog BBexp.
123830-88-80.450.901.221.5601.53170.16
221571-08-40.450.861.301.6901.65410.47
338941-33-20.450.861.542.2401.81190.58
44205-93-00.390.901.361.9201.63690.33
540065-09-60.450.861.381.8701.70670.41
64205-90-70.551.161.341.6001.53170.19
773590-58-60.352.053.182.6702.5161−0.82
828981-97-701.552.502.8962.2041−0.04
984379-13-501.552.842.5202.7008−0.09
1078755-81-401.502.631.9102.0884−0.29
1159467-70-801.382.012.5702.26280.32
1299632-94-701.482.521.9202.2773−0.25
132507-81-50.750.941.521.9061.6051−0.18
14112598-30-80.401.692.162.0702.0043−0.66
157120-01-60.750.801.001.3051.1382−0.04
16104076-38-20.401.382.642.6892.99460.14
17104076-32-60.401.402.692.6942.88980.22
18133099-04-40.491.582.822.8002.2978−0.62
19142494-12-00.001.731.831.4902.65770.16
20486-56-60.001.381.491.0491.38670.04
2154-11-50.001.080.920.8651.37100.56
22494-97-30.130.851.020.9901.23010.32
2358-08-20.051.281.721.5001.36320.01
Table 4. Statistical parameters of cross-validation of MLR models described by Equations (10), (12), (13) and (14); MSE—mean square error, MSEcv*—mean square error of leave-one-out cross validation, MSEcv**—mean square error of leave-50%-out cross validation, PRESS*—predicted residual sum of squares of leave-one-out cross validation, PRESS**—predicted residual sum of squares of leave-50%-out cross validation.
Table 4. Statistical parameters of cross-validation of MLR models described by Equations (10), (12), (13) and (14); MSE—mean square error, MSEcv*—mean square error of leave-one-out cross validation, MSEcv**—mean square error of leave-50%-out cross validation, PRESS*—predicted residual sum of squares of leave-one-out cross validation, PRESS**—predicted residual sum of squares of leave-50%-out cross validation.
MLR ModelStatistical ParametersValues
Equation (10)N23
R20.9039
Q20.8756
MSE0.01784
PRESS*0.55198
MSEcv*0.01784
Equation (12)N65
R20.8474
Q20.8398
MSE0.00481
PRESS*0.3272
PRESS**0.3076
MSEcv*0.004799
MSEcv**0.004799
Equation (13)N65
R20.8469
Q20.8394
MSE0.00482
PRESS*0.3293
PRESS**0.3100
MSEcv*0.004841
MSEcv**0.004841
Equation (14)N65
R20.8471
Q20.8396
MSE0.00482
PRESS*0.3270
PRESS**0.3067
MSEcv*0.004817
MSEcv**0.004817
Table 5. Values of A, B, S, E, V, topological polar surface area (TPSA), polarizability (α), molecular weight (MW) and log BB of tested compounds.
Table 5. Values of A, B, S, E, V, topological polar surface area (TPSA), polarizability (α), molecular weight (MW) and log BB of tested compounds.
Comp. TestedABSEVTPSA
[A2]
α
[A3]
MW
[g/mol]
log BB calculated
Equation (10)Equation (12)Equation (13)Equation (14)
101.341.731.881.814448.2727.55242.280.210.190.200.22
201.341.671.901.955348.2729.30256.300.320.210.220.24
301.331.812.011.936848.2729.37276.720.260.210.220.24
401.281.802.011.936848.2729.37276.720.300.210.220.24
501.331.842.031.936848.2729.37276.720.240.210.220.24
601.261.912.142.059248.2731.20311.170.340.240.240.26
701.672.131.942.029774.5730.27286.29−0.16−0.13−0.11−0.10
801.672.071.972.170674.5732.02300.31−0.04−0.11−0.09−0.08
901.882.262.012.229383.8032.57316.31−0.28−0.23−0.21−0.20
1001.882.262.012.370283.8034.40330.34−0.20−0.21−0.19−0.17
1101.652.242.092.152174.5732.09320.73−0.12−0.10−0.09−0.07
1201.682.141.942.170674.5732.09300.31−0.10−0.11−0.09−0.07
1301.682.081.962.311574.5733.85314.380.02−0.09−0.07−0.05
1401.882.241.982.370283.8034.40330.34−0.20−0.21−0.20−0.18
1501.622.202.072.293074.5733.92334.760.00−0.08−0.07−0.05
1601.662.252.092.293074.5733.92334.76−0.05−0.08−0.07−0.05
1701.592.322.202.415474.5735.74369.200.05−0.06−0.05−0.03
1801.462.072.241.960361.4130.82280.310.060.060.070.09
1901.462.012.262.101261.4132.57294.340.170.070.080.10
2001.462.012.262.101261.4132.57294.340.170.080.080.11
2101.461.952.292.242161.4134.33308.370.290.090.100.12
2201.662.172.282.159970.6433.12310.34−0.04−0.05−0.03−0.02
2301.662.192.312.159970.6433.12310.34−0.05−0.05−0.03−0.01
2401.452.152.362.082761.4132.64314.750.110.070.080.10
2501.402.132.362.082761.4132.64314.750.160.080.090.11
2601.452.182.392.082761.4132.64314.750.090.080.090.11
2701.382.252.502.205161.4134.47349.200.190.100.110.13
2801.382.222.482.205161.4134.47349.200.210.090.100.12
2901.422.192.462.140448.2733.92290.320.150.260.270.29
3001.432.132.482.281348.2735.68304.350.260.270.280.30
3101.432.132.482.281348.2735.68304.350.260.280.290.31
3201.432.132.482.281348.2735.68304.350.260.280.290.31
3301.632.292.502.340057.5036.23320.350.050.150.170.19
3401.632.322.522.340057.5036.23320.350.040.160.170.19
3501.632.322.522.480957.5038.05334.410.110.180.190.21
3601.432.072.502.422248.2737.43318.370.380.290.300.32
3701.412.272.582.262848.2735.75324.760.200.270.280.30
3801.372.252.582.262848.2735.75324.760.250.280.290.31
3901.412.302.612.262848.2735.75324.760.190.280.290.31
4001.342.372.712.385248.2737.57359.210.290.300.310.34
4101.432.192.452.281348.2735.75304.350.220.280.300.31
4201.422.272.602.403748.2737.57338.820.280.300.310.33
4301.422.272.602.403748.2737.57338.820.280.300.310.33
4401.422.272.602.403748.2737.57338.820.280.300.310.33
4501.432.142.482.422248.2737.50318.410.330.300.300.33
4601.432.082.502.563148.2739.26332.440.450.310.320.35
4701.432.082.502.563148.2739.26332.440.450.310.320.35
4801.432.222.632.544648.2739.26352.850.380.310.320.35
4901.432.222.632.544648.2739.33352.850.380.310.320.35
5001.432.222.632.544648.2739.33352.850.380.310.320.35
5101.632.322.522.621857.5039.88348.440.190.190.200.23
5201.642.272.542.903657.5041.63362.470.370.210.220.25
5301.632.412.672.744257.5041.70382.880.230.210.220.25
5401.632.412.672.744257.5041.70382.880.230.210.220.25
5501.632.412.612.744257.5041.70382.880.220.210.220.25
5601.432.222.632.544648.2739.33352.850.380.310.320.34
5701.412.302.602.403748.2737.57338.820.270.290.310.33
5801.412.382.752.526148.2739.40373.270.310.310.330.35
5901.412.382.752.526148.2739.40373.270.310.310.330.35
6001.412.382.752.526148.2739.40373.270.310.320.330.35
6101.432.202.452.422248.2737.58318.370.290.290.300.32
6201.432.142.482.563148.2739.33332.400.410.310.320.34
6301.422.282.582.544648.2739.40352.820.340.300.320.34
6401.422.312.602.544648.2739.40352.820.330.310.320.35
6501.352.382.712.667048.2741.22387.260.430.340.350.37

Share and Cite

MDPI and ACS Style

Janicka, M.; Sztanke, M.; Sztanke, K. Predicting the Blood-Brain Barrier Permeability of New Drug-Like Compounds via HPLC with Various Stationary Phases. Molecules 2020, 25, 487. https://doi.org/10.3390/molecules25030487

AMA Style

Janicka M, Sztanke M, Sztanke K. Predicting the Blood-Brain Barrier Permeability of New Drug-Like Compounds via HPLC with Various Stationary Phases. Molecules. 2020; 25(3):487. https://doi.org/10.3390/molecules25030487

Chicago/Turabian Style

Janicka, Małgorzata, Małgorzata Sztanke, and Krzysztof Sztanke. 2020. "Predicting the Blood-Brain Barrier Permeability of New Drug-Like Compounds via HPLC with Various Stationary Phases" Molecules 25, no. 3: 487. https://doi.org/10.3390/molecules25030487

APA Style

Janicka, M., Sztanke, M., & Sztanke, K. (2020). Predicting the Blood-Brain Barrier Permeability of New Drug-Like Compounds via HPLC with Various Stationary Phases. Molecules, 25(3), 487. https://doi.org/10.3390/molecules25030487

Article Metrics

Back to TopTop