Next Article in Journal
[68Ga]FAPI PET for Imaging and Treatment Monitoring in a Preclinical Model of Pulmonary Fibrosis: Comparison to [18F]FDG PET and CT
Next Article in Special Issue
Discovery of 1H-benzo[d]imidazole-(halogenated) Benzylidenebenzohydrazide Hybrids as Potential Multi-Kinase Inhibitors
Previous Article in Journal
Chemical and Biological Characterization of the Ethyl Acetate Fraction from the Red Sea Marine Sponge Hymedesmia sp.
Previous Article in Special Issue
Unveiling the Multifaceted Capabilities of Endophytic Aspergillus flavus Isolated from Annona squamosa Fruit Peels against Staphylococcus Isolates and HCoV 229E—In Vitro and In Silico Investigations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study of the Lipophilicity and ADMET Parameters of New Anticancer Diquinothiazines with Pharmacophore Substituents

by
Daria Klimoszek
1,
Małgorzata Jeleń
2,*,
Małgorzata Dołowy
1 and
Beata Morak-Młodawska
2
1
Department of Analytical Chemistry, Faculty of Pharmaceutical Sciences in Sosnowiec, Medical University of Silesia in Katowice, Jagiellońska Street 4, 41-200 Sosnowiec, Poland
2
Department of Organic Chemistry, Faculty of Pharmaceutical Sciences in Sosnowiec, Medical University of Silesia in Katowice, Jagiellońska Street 4, 41-200 Sosnowiec, Poland
*
Author to whom correspondence should be addressed.
Pharmaceuticals 2024, 17(6), 725; https://doi.org/10.3390/ph17060725
Submission received: 16 May 2024 / Revised: 28 May 2024 / Accepted: 30 May 2024 / Published: 3 June 2024
(This article belongs to the Special Issue Heterocyclic Compounds in Medicinal Chemistry)

Abstract

:
Lipophilicity is one of the principal parameters that describe the pharmacokinetic behavior of a drug, including its absorption, distribution, metabolism, elimination, and toxicity. In this study, the lipophilicity and other physicochemical, pharmacokinetic, and toxicity properties that affect the bioavailability of newly synthesized dialkylaminoalkyldiquinothiazine hybrids as potential drug candidates are presented. The lipophilicity, as RM0, was determined experimentally by the RP-TLC method using RP18 plates and acetone–TRIS buffer (pH 7.4) as the mobile phase. The chromatographic parameters of lipophilicity were compared to computationally calculated partition coefficients obtained by various types of programs such as iLOGP, XLOGP3, WLOGP, MLOGP, SILCOS-IT, LogP, logP, and milogP. In addition, the selected ADMET parameters were determined in silico using the SwissADME and pkCSM platforms and correlated with the experimental lipophilicity descriptors. The results of the lipophilicity study confirm that the applied algorithms can be useful for the rapid prediction of logP values during the first stage of study of the examined drug candidates. Of all the algorithms used, the biggest similarity to the chromatographic value (RM0) for certain compounds was seen with iLogP. It was found that both the SwissADME and pkCSM web tools are good sources of a wide range of ADMET parameters that describe the pharmacokinetic profiles of the studied compounds and can be fast and low-cost tools in the evaluation of examined drug candidates during the early stages of the development process.

Graphical Abstract

1. Introduction

Heterocyclic compounds are some of the best-known and most important structural components of drugs. Of these, nitrogen-containing heterocycles are particularly important. As the FDA (Food and Drug Administration) data show, 59% of all unique small-molecule drugs contain a nitrogen atom, and it should also be noted that 4 of the 10 most commonly used nitrogen heterocycles also contain a sulfur atom [1]. Phenothiazine is considered to be the third most commonly used six-membered nonaromatic nitrogen heterocycle and is present in 16 unique small-molecule drugs with various effects ranging from antihistaminic, sedative, and antipsychotic effects to anti-neurodegenerative (i.e., Parkinson’s and Alzheimer’s diseases) effects [1,2].
Classical phenothiazines, mainly used as neuroleptics, are substituted at position 10 with dialkylaminoalkyl groups and additionally at position 2 with small groups. These substances have significant neuroleptic, antiemetic, antihistaminic, antipruritic, analgesic, and anthelmintic effects. Continuing research in new directions on the activity of neuroleptic phenothiazines and on the modification of their structures provides information on their anticancer, antiviral (including anti-SARS-CoV-2), antibacterial, and anti-inflammatory activities, and the reversal of multidrug resistance. These substances have antioxidant and antihyperlipidemic effects [3,4,5,6,7,8,9,10].
The lipophilicity of medicinal substances has a significant impact on their ADMET parameters (absorption, distribution, metabolism, excretion, toxicity), which refer to absorption, distribution, metabolism, elimination, and toxicity.
The assessment and consideration of the lipophilicity of medicinal substances are important during drug design, as it can have a significant impact on their pharmacokinetic properties and toxicity [11]. Medicinal substances with moderate lipophilicity tend to be better absorbed through cell membranes, which can affect their rate and absorption efficiency from the gastrointestinal tract or through the skin. Lipophilic substances can more easily penetrate cell membranes and migrate to lipid-rich tissues, which can affect their distribution in the body. Substances with increased lipophilicity may be more susceptible to metabolism in the liver through oxidation, reduction, and conjugation reactions. The impact of lipophilicity on metabolism can have consequences for pharmacological activity and toxicity. The lipophilicity of drugs can influence their excretion, as substances with increased lipophilicity can be stored in fatty tissues and exhibit a prolonged presence in the body. The lipophilicity of drugs may be related to their toxicity, as this can affect their accumulation in tissues and interactions with receptors and proteins in the body [12,13,14,15,16,17]. Therefore, to better understand the behavior of biologically active compounds, including new drug candidates, their lipophilic properties should be assessed. Theoretical and experimental methods are commonly used to describe the lipophilicity of compounds. Calculation methods are used to estimate the lipophilicity parameter quantified as P (partition coefficient) or its decimal logarithm (logP). The extensive development of chemoinformatics has an influence on the number of programs available for the online prediction (in silico) of this important parameter and other ADMET properties, and such platforms include ADMETlab, pkCSM platform, SwissADME, and MetaTox [18]. Calculation approaches are useful to rapidly predict logP values, especially during the early stages of drug development; thus, next, these values should be complemented with experimental data.
Among the experimental techniques, the classic shake-flask method and liquid chromatography play an important role in determining lipophilicity. The lipophilicity chromatographic parameters (RM0) that are obtained by RP-TLC (reversed-phase thin-layer chromatography) and logk0 and assessed by reversed-phase high-performance liquid chromatography (RP-HPLC) are commonly used to assess the lipophilic nature of compounds [19]. The standard shake-flask procedure recommended by the Organization for Economic Co-operation and Development involves the direct measurement of the partition coefficient [20]. It allows for the accurate measurement of the logP values in the range of −2 to 4 but requires relatively large amounts of pure compounds compared to other methods. The main disadvantage of this method is that it is time-consuming and requires the control of many parameters affecting the equilibrium state of the tested system, usually lasting from 1 h to 24 h [21]. Therefore, currently, most lipophilicity tests are conducted by means of chromatographic techniques. Chromatographic approaches in reversed-phase systems (RP-TLC and RP-HPLC) are the most widely used indirect methods to experimentally determine lipophilicity. Both of these chromatographic methods need a smaller amount of sample and a relatively shorter time for analysis compared to the classical shake-flask method. The obtained results are repeatable, and the accuracy of the partition coefficient values can be within ±1 unit in relation to the shake-flask value [22].
A comprehensive review of chromatographic procedures dedicated to the determination of the lipophilicity parameters of different drug substances as an essential tool in medicinal chemistry was performed by Soares and co-workers [23]. Taking into account the importance of lipophilicity parameters as key factors in drug chemistry, namely in the design of new drugs, the aim of this study was to assess the lipophilicity of a newly synthesized group of diquinothiazines by means of both the RP-TLC and calculation methods.
All diquinothiazines that were the subject of this study were tested early for their antiproliferative activity using cultured glioblastoma SNB-19, colorectal carcinoma Caco-2, breast cancer MDA-MB-231, and lung cancer A549 cell lines and NHDF normal fibroblasts [24]. They can therefore be considered as bioactive compounds. Most of the compounds were very active against at least one cancer cell line, with an IC50 value < 3 μM being more active than cisplatin. The most tested diquinothiazines showed higher activity against the A549 lung cancer cell line. Compounds 13, 58, and 1114 were very active against all cancer cells. As was stated, the most active were the dimethylaminopropyldiquinothiazine 2 against the A549 cell line and the pyrrolidinylethyldiquinothiazine 8 against the SNB-19 cell line, with an IC50 value of 0.3 μM. The mechanism of the antiproliferative effect was examined using the RT-QPCR method. It caused a significant reduction in CDKN1A expression in the MDA-MB-231, A549, and SNB-19 tumor lines. Compound 8 markedly reduced the expression of BCL-2 in A549 and SNB-19 and the expression of BAX in cancer cell lines [24].
In another study, the diquinothiazine 2 demonstrated significant in vitro anticancer activity against the human lung carcinoma A549 and non-small lung carcinoma H1299 lines and protective potential for the healthy cell lines BEAS-2B and NHDF. Using the 72 h MTT, strong cytotoxic activity was observed in the viability test (Promega). The weak lethal effect observed in NHDF or BEAS-2B cells at IC50 doses against A549 or H1299 cells confirms promising cancer selectivity. The cell cycle revealed that substance 2 activated the necrosis phase [25].
Continuing from our previous studies, the purpose of this work was to determine the lipophilicity parameters of fifteen newly developed anticancer, angularly condensed diquinothiazines, 115, with pharmacophore dialkylaminoalkyl substituents using combined computational and chromatographic approaches as logPcalcd, RM0, and logPTLC (Figure 1). The full spectral characteristics of these compounds and their synthesis have been well described previously [24]. The current work aimed to discuss the influence of the nature of the substituents and the method of condensation of rings in a five-ring molecule system on the value of lipophilicity indices, determined by both calculation and RP-TLC methods, as well as on other drug-likeness and ADME properties predicted by in silico studies that are key for describing the pharmacokinetic behavior of these drugs.
The usefulness of the RP-TLC technique as well as logP predictions using different computational software for the design of promising drug candidates belonging to the studied diquinothiazines with pharmacophore dialkylaminoalkyl substituents was evaluated.

2. Results

2.1. Lipophilicity Studies

Both the computational and chromatographic values of the lipophilicity parameters logPcalcd, RM0, and logPTLC were determined for fifteen 7- and 14-substituted angularly condensed diquinothiazines with pharmacophore dialkylaminoalkyl substituents on the thiazine nitrogen atom. The tested diquinothiazines were divided into three groups, differing in the way the quinoline rings are connected to the 1,4-thiazine ring: 7-substituted diquino[3,2-b;3′,4′-e]thiazines, 15; 7-substituted diquino[3,2-b;6′,5′-e]thiazines, 610; and 14-substituted diquino[3,2-b;8′,7′-e]thiazines, 1115 (Figure 1).
Firstly, the computational method was chosen to determine the lipophilicity parameters of the studied compounds. For this purpose, popular computational programs were used based on various mathematical algorithms available on the following platforms: SwissADME [26], pkCSM [27], Molinspiration [28], and the ChemDraw Ultra program [29].
The calculated logP values (logPcalcd) for the angularly fused diquinothiazines 1–15 are shown in Table 1 and differed depending on the substituents on the thiazine nitrogen atom, the shape of the five-ring diquinothiazine system, and on the calculation program.
The next step of our study focused on determining the more reliable lipophilicity parameters of the tested compounds. In order to obtain the relative lipophilicity, expressed as RM0, the chromatographic behavior of the fifteen tested diquinothiazines, 115, was investigated under proper RP-TLC conditions. RP-18 plates were used as the stationary phase, while organic modifiers containing acetone were used as the mobile phase. The linear relationship between RM0 and acetone concentration was determined on the basis of Equation (2) (Table 2). In addition to this, thanks to the relationship observed between RM0 and the slope of these linear plots (b), the lipophilicity parameter C0 was also determined (Table 2).
Then, the relative lipophilicity parameter RM0 was converted to an absolute value lipophilicity parameter, logPTLC, using a calibration curve determined under the same measurement conditions for a set of standards, IV, with the literature values of logPlit in the range of 1.21–6.38 (Table 3). The obtained values of the RM0 coefficient of the tested compounds were in the range of 3.01–3.83. The correlation between the logPlit values and the experimental RM0 values for standards IV gave the following calibration equation:
log PTLC = 1.2838 RM0 + 0.2138
(r = 0.9967; s = 0.1920; F = 459.32; p < 0.001)
The standard curve equation was used to obtain the logPTLC parameter for compounds 115, and the results are presented in Table 4.

2.2. Molecular Descriptors

For all tested compounds, 115, selected molecular descriptors such as molar mass (M), molar volume (VM), molar refraction (RefM), and surface area were calculated to check how they correlate with the experimentally determined lipophilicity parameter RM0 (Table 5).

2.3. In Silico ADME Prediction

Physicochemical parameters used to predict drug-like properties, i.e., Lipinski’s rule of five and the Ghose, Veber, Egan, and Muegge rules, were also calculated (Table 6).
ADME parameters were obtained using the pkCSM and PreADMET servers (Table 7, Table 8 and Table 9). As shown in Table 7, Table 8 and Table 9, the compounds studied showed significant differences in the molecular descriptors as well as in their ADME parameters.

2.3.1. Absorption

The potential absorption of 15 tested angularly condensed diquinothiazines, 1–15, was evaluated using the parameters of water solubility, Caco-2 cell permeability (human colon adenoma cells), absorption in the human intestine, and skin permeability, obtained using the pkCSM platform. The obtained results are summarized in Table 7. Water solubility was measured using the logS parameter (S—solubility, expressed in mol/L). When it comes to the absorption of administered drugs, the most frequently used method to determine this parameter is to test the permeability of potential drugs through a monolayer of Caco-2 cells. This is due to the similarity in structure and function of Caco-2 cells to the human intestinal epithelium. Since the main site of absorption of an orally administered drug is usually the intestine, it is important to determine the amount of the compound that is absorbed here. This method predicts the percentage of the compound that was absorbed. The skin permeability parameter is also significant for the effectiveness of some products. This parameter is expressed as logKp, and a compound is considered to have relatively low skin permeability if logKp > −2.5 [33,34].

2.3.2. Distribution

The potential distribution of the tested diquinothiazines 115 was assessed using the parameters of VDss, fraction unbound, BBB permeability, and CNS permeability, obtained using the pkCSM platform. The results obtained are summarized in Table 8. The volume of distribution (VDss) provides an indication of the distribution of the drug in the body and is a pharmacokinetic parameter representing the volume into which the dose of drug would have to be distributed to give rise to the same concentration observed in the blood plasma. A low VDss indicates high water solubility or high plasma protein binding because more of the drug remains in the plasma; a high VDss suggests significant concentration in tissues, for example, due to tissue binding or high lipid solubility [35,36]. According to the model used in the pkCMS software, VDss is considered low if log VDss < −0.15 and high if log VDss > 0.45.
Fu (unbound fraction) is also an important pharmacokinetic parameter because it affects various factors of drug effectiveness and side effects (including glomerular filtration in the kidneys, total clearance, and hepatic metabolism). Therefore, it is important to accurately predict Fu during drug development. Unbound drugs in plasma may exhibit pharmacological activity by interacting with targets such as proteins, enzymes, receptors, and channels; hence, the plasma unbound fraction (Fu) of a drug is an important factor in determining drug efficacy [37].
When considering a substance as a drug candidate, it is also important to determine the extent to which it will cross the blood–brain barrier. This parameter is measured in vivo as logBB, the logarithmic ratio of brain to plasma drug concentrations. According to the computational model used in pkCSM, if the value of the logBB parameter is greater than 0.3, the substance crosses the blood–brain barrier, while if logBB is lower than −1, the substance is distributed to the brain to a small extent [38].

2.3.3. Excretion and Toxicity

The potential excretion and toxicity of the tested diquinothiazines 115 were evaluated using the parameters of total clearance, maximum tolerated dose, acute oral toxicity of rats, Tetrahymena pyriformis toxicity, and minnow toxicity, obtained using the pkCSM platform. The results obtained are summarized in Table 9. The expected total clearance, expressed as log (ml/min/kg), is the volume of plasma completely cleared of the drug per unit of time by the organ eliminating the drug from the body. Knowledge of this parameter is necessary to determine the maintenance dose of the drug [39]. The maximum recommended tolerated dose (MRTD) allows us to initially determine the toxic dose of a chemical substance to humans. The model used in pkCSM was developed using 1222 experimental data points from human clinical trials. The obtained results of the MRTD parameter values are given as log (mg/kg/day). Calculating this parameter is helpful in determining the maximum recommended starting dose for drugs under investigation. If the calculated value is less than or equal to 0.477 log (mg/kg/day), it is considered low, and if it is higher, it is considered high. A model based on Tetrahymena pyriformis, a protozoal bacterium whose toxicity is often considered a toxic endpoint, was used to determine in silico toxicity in the pkCSM program. This method was built using the concentrations of 1571 compounds required to inhibit 50% of growth. This parameter is designated as plGC50 (negative logarithm of the concentration required to inhibit 50% growth in log µg/L), and a compound can be assessed as toxic if the value of this parameter is calculated to be greater than −0.5 log µg/L. The minnow toxicity parameter is based on measurements of the LC50 parameter, i.e., the concentration of the substance necessary to cause the death of 50% of flathead minnows, which were used as an animal model in this study. The model was based on tests of the LC50 parameter for over 550 substances. According to this model, an LC50 value below 0.5 mM (logLC50 < −0.3) is considered to indicate high acute toxicity.

3. Discussion

It is known that lipophilicity is one of the most frequently studied physicochemical parameters of new substances that are drug candidates [19,40]. It is most often defined to support quantitative structure–activity relationships (QSARs), including absorption, distribution in tissues, and drug transport across the biological barrier [41,42,43]. This parameter is characterized by the distribution of a dissolved substance in two-phase liquid–liquid or solid–liquid systems. In silico methods have also been proposed to assess lipophilicity. Several programs have been developed to calculate the value of the logP parameter [26,27,28,29]. It is still important to use experimental methods because most programs determining the logP parameter calculate lipophilicity using atomic methods that do not take into account some structural parameters. There may be discrepancies between the calculated data and those determined experimentally, which can be explained by the fact that newly synthesized drug candidates may contain substructures or heterocycle systems that are not covered in the software development training set [44,45].
The present research on lipophilicity began with in silico analyses using various computer programs available on the SwissADME [26], pkCSM [27], and Molinspiration [28] platforms and the ChemDraw program [29]. These platforms use various mathematical modules described on the website of the above-mentioned servers. The obtained results of the calculated lipophilicity are within a wide range of values. This is most likely due to differences in the calculation models. The comparison of the chromatographic parameters (logP) of the tested angularly condensed N-dialkylaminoalkyldiquinothiazines 115 is shown in Figure 2. As observed (Table 1), the lowest value of the logPcalc parameter was obtained for the 14-(3′-dimethylaminopropyl)diquino[3,2-b;8′,7′-e]thiazine 12 according to calculations with the iLOGP program (SwissADME) and the highest one was obtained for all three isomeric diquinothiazines 5, 10, and 15 with a piperidinylethyl substituent based on calculations with the logP program (pkCSM). This program also calculated the same values of the logP parameter for the remaining isomers (i.e., for diquinothiazines 1, 6, and 11, logP = 5.17; for 2, 7, and 12, logP = 4.60; for 3, 8, and 13, logP = 4.81; and for 4, 9, and 14, logP = 5.22).
In order to compare all obtained theoretical values of the partition coefficient of the examined drugs expressed in the form of logP, and to then estimate the lipophilic character of the studied N-dialkyloaminoalkyldiquinothiazines, a chemometric approach, i.e., cluster analysis (CA) with Euclidean distance, was conducted. The results are presented in Figure 3 and Figure 4.
Figure 3 shows the dendrogram of all the theoretical lipophilicity parameters of the examined angularly condensed N-dialkyloaminoalkyldiquinothiazines 115. The first cluster contains the partition coefficients LogPa, logPaverage, logPc, WLOGP, milogP, and XLOGP3. In the second cluster, the rest of the computed partition coefficients are placed: SILICOS-IT, MLOGP, and iLOGP. The reason for these differences is the prediction power of each software. This grouping of partition coefficient values shows the biggest similarity (the smallest distance on the dendrogram) of logPc with the logPaverage value calculated on the basis of all theoretical logP values. Therefore, the parameter logPc, may be a good alternative for estimation of the lipophilicity of the examined group of N–dialkyloaminoalkyldiquinothiazines.
Figure 4 shows the similarity analysis of the examined compounds based on their partition coefficients. According to the computed lipophilicity parameters (the logP values), the studied compounds can be divided into two groups based on variation in their structural particularities, i.e., depending on the kind of substituent that may influence the lipophilic character of these compounds. The first group consists of the compounds C13, C8, and C3 in the first subgroup, with C12, C7, and C2 in the second subgroup. Similarly, the second group is divided into two smaller visible subgroups, with C15, C10, and C5 in the first subgroup and C14, C9, C4, C11, C6, and C1 in the second subgroup.
The cluster analysis shown in Figure 4 confirms the greatest similarity in lipophilic character of derivatives C10 and C5.
In the next steps of this study, the calculated logP values using different computational methods, summarized in Table 1, were correlated with the experimental lipophilicity parameter RM0, presented in Table 2. The correlation coefficients and the equations obtained for the correlation of the RM0 parameter with the logP values obtained by means of an appropriate program are summarized in Table 10, which shows the strongest relationships between logP and RM0.
In order to obtain the linear equations listed in Table 10, RP-TLC was chosen as the method for experimental determination of the lipophilicity parameter of the tested diquinothiazines 115. Using this technique, the relative lipophilicity parameter RM0 was obtained. It was then converted to the relative parameter logPTLC (as described in Section 2 and Section 4). For all investigated derivatives, in a wide range of organic modifier concentrations in the mobile phase, high values of correlation coefficients (r = 0.99) made it possible to determine the lipophilicity parameter RM0 by extrapolation (Table 2). The analysis of the chromatographic parameter RM0 of the examined angularly condensed N-dialkyloaminoalkyldiquinothiazines 115 presented in Table 2 shows that they are relatively smaller compared to the theoretical logP values obtained by means of different computer software and their mean value (logPaverage). Generally, a difference in the obtained RM0 within one unit to logP was noted in the case of all chromatographic descriptors. This fact confirms that calculation approaches are useful only for the rapid prediction of logP values during the first stage of study of new drug candidates, after which they should be complemented with experimental data such as chromatographic descriptors. Of all the algorithms, iLogP showed the biggest similarity to the chromatographic value, especially for compounds 2, 7, 11, 12, 14, and 15. This shows the potential utility of iLogP for the rapid estimation of the lipophilic character of these compounds.
In the next stage, satisfactory relationships between the RM0 values and b (slope) values of the linear equations allowed the calculation of the next lipophilicity parameter, namely C0. The results of these data confirm that all compounds studied belong to a congeneric group. However, as shown in Figure 5, the dissimilarity of the C0 values with the partition coefficients as well as the chromatographic parameter (RM0) indicates that both parameters cannot be fully replaced by C0.
The linear equations obtained between the theoretical partition coefficients and chromatographic parameters of lipophilicity of the tested compounds 115 (Table 10) confirmed the usefulness of iLOGP, XLOGP3, WLOGP, MLOGP, LOGP, miLOGP, logPaverage, and SILICOS-IT for predicting the experimental value of this parameter in the form of RM0.
Analyzing the experimentally obtained values of the logPTLC parameter, it can be observed that in each of the three groups of isomeric angularly condensed diquinothiazines, the lowest lipophilicity was found in derivatives with a dimethylaminopropyl substituent (2, 7, and 12), while the highest was found in derivatives with an N-methylpiperidine substituent (5, 10, and 15). Comparing the three tested groups of isomeric angularly condensed diquinothiazines, it can be concluded that the compounds with the 7-substituted diquino[3,2-b;3′,4′-e]thiazine structure, 1–5, are characterized by the lowest lipophilicity, while those with the 7-substituted diquino[3,2-b;6′,5′-e]thiazine structure, 6–10, show the highest values of the logPTLC parameter.
Selected molecular descriptors were also calculated for all tested diquinothiazines. The values are presented in Table 5. As shown by the results in Table 5, the same values of these parameters were obtained for isomeric angularly condensed diquinothiazines with the same substituents. The results obtained correlated with the experimentally obtained RM0 parameter (Table 11). Good correlations were obtained for all descriptors (r ≥ 0.8). Correlations in subgroups of individual isomers resulted in improved correlations.

In Silico ADME Prediction

Potential drug candidates usually have similar physicochemical properties. Therefore, based on analysis of the simple molecular properties of existing drugs and/or drug candidates, simple filters defining acceptable limits for these properties are being developed. In 2007, Lipinski proposed the ‘Rule of Five’ [46], the most famous drug similarity filter which defines four rules that determine whether a molecule is well absorbed or not after oral administration. These parameters include molecular weight (MW ≤ 500), octanol/water partition coefficient (ClogP ≤ 5), number of hydrogen bond donors (HBD ≤ 5), and number of hydrogen bond acceptors (HBA ≤ 10). If the relationship violates two or more bases, it may not be active when administered orally. Data in the literature indicate that 85.4% of FDA-approved drugs meet the rule of five [47]. The results obtained on the SwissADME platform show that all tested diquinothiazines meet Lipiński’s rule.
Since then, various principles have been developed similar to the Rule of Five. For example, after parsing a file of 6304 molecules in a database, Ghose et al. found that over 80% of compounds from the Comprehensive Medicinal Chemistry (CMC) database meet the following qualification ranges: AlogP between −0.4 and 5.6, MW between 160 and 480, molar refractive index from 40 to 130, and integer atoms from 20 to 70 [48]. From the current group of 15 tested diquinothiazines, according to the SwissADME platform, derivatives with piperidine (4, 9, and 14) and N-methylpiperidine (5, 10, and 15) substituents do not meet this rule because their molar refractive index is greater than 130. However, the value of this parameter obtained with the program Chem3D is less than 130 (Table 5).
All tested substances, 115, also meet Veber and Egan’s rules. According to Veber’s rule, compounds that satisfy only the following two criteria will most likely have good oral bioavailability in rats: have 10 or fewer rotatable bonds and a polar surface area equal to or less than 140 Å2 (or 12 or fewer hydrogen bond donors and acceptors). This may be because a reduced polar area correlates better with an increased permeation rate than lipophilicity, and an increased number of rotatable bonds has a negative effect on the permeation rate [49,50]. In turn, Egan’s rule is based on two descriptors: PSA and AlogP98.
It has been shown that molecular weight, although often used in passive absorption models, is unnecessary because it is already a component of both PSA and AlogP98. Following extensive model validation on hundreds of known orally administered drugs, “drug-like” molecules that have been tested for Caco-2 cell permeability have shown a good rate of successful predictions (74–92%) [51]. However, none of the tested substances meet Muegge’s rule due to the calculated XLOGP value being greater than 5.
The bioavailability score calculated using the SwissADME platform for all tested substances was 0.55. Substances with a bioavailability score ≥ 0.55 are considered very well absorbed by the body. A bioavailability score of 0.55 and a TPSA score between 60 and 140 Å2 indicate optimal absorption [52].
When developing new drugs, it is also important to determine and evaluate ADME properties (absorption, distribution, metabolism, and excretion). Poor ADME properties are a significant cause of failure and increased drug design costs [53]. It was considered acceptable to determine these parameters using in silico methods at the early stages of substance evaluation for drug suitability. This allows, without incurring high research costs, to eliminate compounds with undesirable ADME properties at an early stage of research [54]. For the discussed substances, selected parameters responsible for absorption (water solubility, Caco-2 permeability, intestinal absorption, and skin permeability), distribution (VDss, unbound fraction, BBB permeability, and CNS permeability), and excretion and toxicity (total clearance, max. tolerated dose, oral rat acute toxicity, oral rat chronic toxicity, Tetrahymena pyriformis toxicity, and minnow toxicity) were also determined in silico. The pkCSM platform was used to determine these parameters. All 15 tested compounds are characterized by poor solubility in water due to their chemical structure (five six-membered fused rings and no hydrophilic groups). The obtained parameter values range from −5.871 for diquinothiazine 1 to −3.965 for diquinothiazine 13. Caco-2 permeability obtained using the pkCSM platform is given as the logarithm of the apparent permeability coefficient (log Papp). Compounds with a predicted log Papp at 10−6 cm/s greater than 0.9 are considered to have high Caco-2 permeability. According to the analysis results, all the synthesized compounds show high Caco-2 cell permeability (1.003–1.226). The calculated absorption values in the human intestine show that all compounds have a similar, very high probability of intestinal absorption (92.24–96.51%). The lowest values were obtained for derivatives with the 7-substituted diquino[3,2-b;3′,4′-e]thiazine structure, i.e., 15. For the other two groups of isomers, the values were similar, which shows the influence of the shape of the molecule on this parameter. The calculated skin permeability values are in the range of −2.744 to −2.694 and indicate poor permeability (Table 7). Almost all tested angularly condensed diquinothiazines can be substrates for P-glycoprotein (except for substance 1), and all of them are predicted to inhibit P-glycoproteins I and II (Table S1).
VDss is a pharmacokinetic parameter that determines half-life and represents the degree of drug distribution in tissues. The predicted log VDss values ranged from 0.86 to 1.37, which means that the compounds have a high constant distribution in plasma and tissues. The unbound fraction in plasma (Fu) is also an important pharmacokinetic parameter which determines the amount of drug that is “free” in the plasma, and therefore the fraction capable of diffusing from the plasma into tissues. From the results obtained for the tested diquinothiazines 115, the unbound fraction had low values ranging from 0.189 to 0.268. The BBB permeability parameter was also determined using the pkCSM platform. The blood–brain barrier controls the transfer of essential substances needed for a proper functioning brain and participates in the removal of cellular metabolites and toxins. The expected permeability through the blood–brain barrier for the tested substances ranged from 0.357 to 0.572, which means that they will pass through the BBB barrier. All tested diquinothiazines may also be able to penetrate the central nervous system because the logPS parameter values obtained for them were in the range of −1.304 to −1.503, and substances with logPS > −2 are considered to penetrate the CNS, while those with logPS < −3 do not penetrate the CNS (Table 8).
The possible interaction of the tested diquinothiazines with cytochrome P450 was also calculated. The obtained results show that all tested substances may be CYP1A2 inhibitors, and almost all of them may be CYP2D6 (except substances 3 and 4) and CYP3A4 (except substance 2) inhibitors; however, none of the tested compounds may be CYP2C9 inhibitors (Table S1). Although the abundance of CYP3A4 is poor, it contributes to nearly 50% of drug metabolism, so drug–drug interactions need to be considered in the subsequent development of compounds.
The excretion and toxicity parameters determined in silico are presented in Table 9. Total clearance of the tested diquinothiazines ranged from 0.526 to 0.832. Total clearance measures the efficiency of drug elimination from the entire body and is useful in determining the rate of drug dosing. The maximum tolerated dose value for the tested diquinothiazines was calculated in the range from 0.188 to 0.672. A dose equal to or less than 0.477 is considered low, and such values were obtained for diquinothiazines 610, while dose values higher than 0.477, which are considered high, were obtained for the remaining compounds tested. The values of the parameters of oral rat acute toxicity and chronic toxicity were calculated in the range from 2.272 to 3.184 and from 0.449 to 1.096, respectively. The determined values of the Tetrahymena pyriformis toxicity parameter were above −0.5 (0.291–0.331), which means that the compounds may be toxic, while the minnow toxicity parameter values ranged from −0.989 to 1.944. Values of this parameter of less than −0.3 may indicate the toxicity of the substance and were calculated for diquinothiazines 4, 8, 9, 13, and 14. It was also checked whether the tested substances could be substrates of Organic Cation Transporter 2 (Renal OCT2). Renal OCT2 is a renal uptake transporter that plays an important role in renal clearance and the clearance of drugs and endogenous compounds. Diquinothiazines 8 and 9 may have this effect (Table S2). However, none of the substances discussed should show AMES toxicity or irritate the skin, but they may be hepatotoxic. Calculations using the pkCSM program also showed that the tested substances may be hERG II inhibitors but should not show any inhibition against hERG I (Table S1).
Absorption from the gastrointestinal tract and access to the brain are two pharmacokinetic parameters important to estimate at different stages of the drug discovery process. For this purpose, an estimated method of penetration into the brain or intestines has been developed—boiled egg [55]. It is an accurate predictive model that works by calculating the lipophilicity and polarity of small molecules. At the same time, penetration into both the brain and intestines is predicted based on the same two physicochemical descriptors. This computational method was also used to analyze the tested diquinothiazines 115 (Figure 6). Molecules located in the white region are expected to be passively absorbed in the gastrointestinal tract, while molecules in the “yolk” region are expected to passively cross the BBB. It can be observed that none of the analyzed compounds can be passively absorbed in the intestines, but all of them can pass through the BBB. This shows that the dosing compliance may be poor due to the inability of all the compounds to be absorbed through the intestines. Additionally, all tested compounds can become substrates for p-glycoprotein (blue points).

4. Materials and Methods

4.1. Solvents and Reference Standards

The studied compounds, i.e., 7-substituted diquino[3,2-b;3′,4′-e]thiazines 15, 7-substituted diquino[3,2-b;6′,5′-e]thiazines 610, and 14-substituted diquino[3,2-b;8′,7′-e]thiazines 1115, were previously designed and synthesized [24]. The chemical structures of the examined compounds are shown in Figure 1. The summary formulas, molecular weights, melting points, and physical appearances of the tested diquinothiazines 1–15 are listed in Table S3.
Analytical grade acetone (POCh, Gliwice, Poland) and TRIS (tris(hydroxymethyl)aminomethane, Fluka, Buchs, Switzerland) were used in RP-TLC studies as the mobile phase components. To prepare the calibration curve, the reference standards of five chemical compounds with the described lipophilicity parameter (logPlit) were used: acetanilide (I, 1.21 [30], POCh, Gliwice, Poland), benzoic acid (II, 1.87 [31], POCh, Gliwice, Poland), benzophenone (III, 3.18 [31], Fluka, Buchs, Switzerland), anthracene (IV, 4.45 [31], POCh, Gliwice, Poland), and p,p′-DDT (1,1,1-trichloro-2,2-bis(4-chlorophenyl, Merck, Darmstadt, Germany), V, 6.38 [32]).

4.2. RP-TLC Analysis

All chromatographic experiments were performed on a silica gel 60 RP-18 F254S (10 cm × 10 cm) RP-TLC plate (Merck, Darmstadt, Germany). The mobile phase was prepared by mixing the respective amount of aqueous TRIS (tris(hydroxymethyl)aminomethane) buffer, pH = 7.4 (ionic strength 0.2 M), to meet physiological conditions and acetone in a range from 50 to 80% (v/v) in increments of 5%. Ethanol (96%) was used to prepare the solutions of angularly condensed diquinothiazines 115 and the standards IV at a concentration of 2.0 mg/mL. Next, solutions (2 mL) of the analyzed compounds were applied to the plates 5 mm apart and 10 mm from the lower edge and sides of the plates. Before plate development, the chromatographic chambers were saturated with the mobile phase for 0.5 h. After the development of the plates and drying in a stream of air, the chromatograms were observed in ultraviolet light at 254 nm. Each chromatographic experiment was run in triplicate, and then mean values of the retardation factor, i.e., RF values, were used to determine the RMW values of the compounds.
The RM values calculated from the experimental RF values by using Equation (1) were linearly dependent on the concentration of acetone:
R M = log 1 R F 1
The RM0 values were obtained by extrapolating to zero acetone concentration by using Equation (2):
R M = R M 0 + b C
where C is the volume fraction of the organic modifier in the mobile phase and b is the change in the RM value due to the 1% increase in the organic modifier in the mobile phase (associated with the specific hydrophobic surface area) [56].
Next, the chromatographic data in the form of RM0 and coefficient b were used to calculate the lipophilicity parameter C0 [56] from the following Equation (3):
C 0 = R M 0 b
The C0 parameter is interpreted as the molecular hydrophobicity per unit of the specific hydrophobic surface, i.e., corresponding to the specific hydrophobic surface area of the substance in contact with the stationary phase.

4.3. In Silico Calculated Descriptors

Theoretical values of the partition coefficients (logP) for each compound were determined using various Internet servers such as SwissADME [26], pkCSM [27], Molinspiration [28], and the ChemDraw program [29], including iLOGP, XLOGP3, WLOGP, MLOGP, SILCOS-IT, LogP, logP, and milogP. The logPaverage is an arithmetic mean of the values predicted by the algorithms. The software used are based on different methods of logP calculation [52,57,58]. The following lipophilicity descriptors were obtained using the web tool SwissADME: iLOGP, a physics-based method that relies on free energies of solvation in n-octanol and water calculated by the generalized-born and solvent accessible surface area (GB/SA) model developed by Daina and coworkers; WLOGP, an atomistic method stationed on fragments and topological descriptors; XLOGP3, an atomistic method including corrective factors and a knowledge-based library; MLOGP, based on topological indices and the linear relationship of structure–logP; and SILICOS-IT, a hybrid fragment/topological approach based on 27 fragments and 7 topological descriptors [26,52,57,58]. The next method used for the calculation of miLogP is a fragment-based approach developed by Molinspiration, and it calculates log P from the sum of group or fragment contributions and correction factors [28]. To compute all other descriptors, such as LogPa and logPc, the pkCSM and ChemDraw programs were used, respectively [27,29].
The molecular descriptors and the ADME parameters were also obtained using the SwissADME and pkCSM platform [26,27].
The cluster analysis procedure was conducted using Statistica software v. 13.3. Calculations were performed on Euclidean distances and a single linkage distance [59].

5. Conclusions

In the search for the most promising drug-like substances, the key parameter that describes the pharmacokinetic behavior of a drug is lipophilicity. Therefore, in the present work, fifteen newly synthesized anticancer, angularly condensed diquinothiazines with pharmacophore dialkylaminoalkyl substituents were studied with the goal of predicting selected pharmacokinetic parameters influencing their ADMET properties, including the lipophilicity descriptors logP and RM0. Lipophilicity indices were determined experimentally by the RP-TLC method, with RP18 as the stationary phase and an acetone–TRIS buffer (pH 7.4) mixture as the mobile phase, as well as theoretically using computer approaches. In addition to this, the selected ADMET parameters were determined in silico using the SwissADME and pkCSM platforms and then correlated with the experimental lipophilicity parameters of the examined compounds. The results of the lipophilicity study confirm that the applied software can be useful for the rapid prediction of logP values during the first stage of study of the examined drug candidates, and next they should be complemented with experimental data such as chromatographic descriptors. Of all the algorithms used, iLogP showed the biggest similarity to the chromatographic value (RM0), especially for compounds 2, 7, 12, 14, and 15. In addition to this, it was stated that both the SwissADME and pkCSM web tools are good sources of a wide range of ADMET parameters that describe the pharmacokinetic profiles of the studied compounds and can be fast and low-cost tools in the evaluation of drug candidates during the early stages of the development process.
The satisfactory linear correlation equations of the chromatographic parameter of lipophilicity (RM0) with the theoretical logP values and other physicochemical properties of the tested compounds, such as the molar mass, molar volume, molar refractive index, and surface area, and with the results of other pharmacokinetics parameters that are key in the ADMET profile demonstrate the utility of the computational approaches in estimating the physicochemical properties of new drug candidates in the development process.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ph17060725/s1, Table S1: The absorption and metabolism descriptors for compounds 115; Table S2: The excretion and toxicity descriptors for compounds 115; Figure S1: The bioavailability radars for compounds 115; Table S3: The summary formulas, molecular weights, melting points, and appearance of the tested diquinothiazines 115.

Author Contributions

Conceptualization, M.J.; methodology, M.J. and D.K.; software, M.J., D.K. and B.M.-M.; validation, M.J. and M.D.; formal analysis, M.J. and D.K.; investigation, M.J. and D.K.; writing—original draft preparation, M.J. and M.D.; writing—review and editing, M.J., M.D., D.K. and B.M.-M.; visualization, B.M.-M.; supervision, M.J.; project administration, M.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Medical University of Silesia in Katowice, grant number BNW-1-013/N/3/F.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data from the research described in the manuscript are available from the authors.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Vitaku, E.; Smith, D.T.; Njardarson, J.T. Analysis of the structural diversity, substitution patterns, and frequency of nitrogen heterocycles among US FDA approved pharmaceuticals. J. Med. Chem. 2014, 57, 10257–10274. [Google Scholar] [CrossRef]
  2. Ohlow, M.J.; Moosmann, B. Foundation Review: Phenothiazine: The seven lives of pharmacology’s first lead structure. Drug Discov. Today 2011, 16, 119–131. [Google Scholar] [CrossRef]
  3. Edinoff, A.N.; Armistead, G.; Rosa, C.A.; Anderson, A.; Patil, R.; Cornett, E.M.; Murnane, K.S.; Kaye, A.M.; Kaye, A.D. Phenothiazines and their evolving roles in clinical practice: A narrative review. Health Psychol. Res. 2022, 10, 38930. [Google Scholar] [CrossRef]
  4. Posso, M.C.; Domingues, F.C.; Ferreira, S.; Silvestre, S. Development of phenothiazine hybrids with potential medicinal interest: A Review. Molecules 2022, 27, 276. [Google Scholar] [CrossRef]
  5. Pluta, K.; Morak-Młodawska, B.; Jeleń, M. Recent progress in biological activities of synthesized phenothiazines. Eur. J. Med. Chem. 2011, 46, 3179–3189. [Google Scholar] [CrossRef]
  6. Spengler, G.; Csonka, Á.; Molnár, J.; Amaral, L. The anticancer activity of the old neuroleptic phenothiazine-type drug thioridazine. Anticancer Res. 2016, 36, 5701–5706. [Google Scholar] [CrossRef]
  7. Zubas, A.; Ghinet, A.; Farce, A.; Dubois, J.; Bîcu, E. Phenothiazineand carbazole-cyanochalcones as dual inhibitors of tubulin polymerization and human farnesyltransferase. Pharmaceuticals 2023, 16, 888. [Google Scholar] [CrossRef]
  8. Mosnaim, A.D.; Ranade, V.V.; Wolf, M.E.; Puente, J.; Valenzuela, A.M. Phenothiazine molecule provides the basic chemical structure for various classes of pharmacotherapeutic agents. Am. J. Ther. 2006, 13, 261–273. [Google Scholar] [CrossRef]
  9. Chuang, S.T.; Papp, H.; Kuczmog, A.; Eells, R.; Condor Capcha, J.M.; Shehadeh, L.A.; Jakab, F.; Buchwald, P. Methylene blue is a nonspecific protein–protein interaction inhibitor with potential for repurposing as an antiviral for COVID-19. Pharmaceuticals 2022, 15, 621. [Google Scholar] [CrossRef]
  10. Boyd-Kimball, D.; Gonczy, K.; Lewis, B.; Mason, T.; Siliko, N.; Wolfe, J. Classics in chemical neuroscience: Chlorpromazine. ACS Chem. Neurosci. 2019, 10, 79–88. [Google Scholar] [CrossRef]
  11. Ginex, T.; Vazquez, J.; Gilbert, E.; Herrero, E.; Luque, F.J. Lipophilicity in drug design: An overview of lipophilicity descriptors in 3D-QSAR studies. Future Med. Chem. 2019, 11, 1177–1193. [Google Scholar] [CrossRef]
  12. Johnson, T.W.; Dress, K.R.; Edwards, M. Using the Golden Triangle to optimize clearance and oral absorption. Bioorg. Med. Chem. Lett. 2009, 19, 5560–5564. [Google Scholar] [CrossRef]
  13. Johnson, T.W.; Gallego, R.A.; Edwards, M.P. Lipophilic efficiency as an important metric in drug design. J. Med. Chem. 2018, 61, 6401–6420. [Google Scholar] [CrossRef]
  14. Ditzinger, F.; Price, D.J.; Ilie, A.R.; Köhl, N.J.; Jankovic, S.; Tsakiridou, G.; Aleandri, S.; Kalantzi, L.; Holm, R.; Nair, A.; et al. Lipophilicity and hydrophobicity considerations in bioenabling oral formulations approaches—A PEARRL review. J. Pharm. Pharmacol. 2019, 71, 464–482. [Google Scholar] [CrossRef]
  15. Leeson, P.D.; Davis, A.M. Time-related differences in the physical property profiles of oral drugs. J. Med. Chem. 2004, 47, 6338–6348. [Google Scholar] [CrossRef]
  16. Markovic, M.; Ben-Shabat, S.; Keinan, S.; Aponick, A.; Zimmermann, E.M.; Dahan, A. Lipidic prodrug approach for improved oral drug delivery and therapy. Med. Res. Rev. 2019, 39, 579–607. [Google Scholar] [CrossRef]
  17. Di, L.; Kerns, E.H. Profiling drug-like properties in discovery research. Curr. Opin. Chem. Biol. 2003, 7, 402–408. [Google Scholar] [CrossRef]
  18. Dulsat, J.; López-Nieto, B.; Estrada-Tejedor, R.; Borrell, J.I. Evaluation of free online ADMET tools for academic or small biotech environments. Molecules 2023, 28, 776. [Google Scholar] [CrossRef]
  19. Rutkowska, E.; Pająk, K.; Jóźwiak, K. Lipophilicity—Methods of determination and its role in medicinal chemistry. Acta Pol. Pharm. Drug Res. 2013, 70, 3–18. [Google Scholar]
  20. OECD. Test No. 117: Partition Coefficient (n-Octanol/Water), HPLC Method; OECD Publishing: Paris, France, 2022. [Google Scholar]
  21. Roman, I.P.; Mastromichali, A.; Tyrovola, K.; Canals, A.; Psillakis, E. Rapid determination of octanol-water partition coefficient using vortex-assisted liquid-liquid microextraction. J. Chromatogr. A 2014, 1330, 1–5. [Google Scholar] [CrossRef]
  22. Eadsforth, C.V.; Moser, P. Assessment of reverse-phase chromatographic methods for determining partition coefficients. Chemosphere 1983, 12, 1459–1475. [Google Scholar] [CrossRef]
  23. Soares, J.X.; Santos, Á.; Fernandes, C.; Pinto, M.M.M. Liquid chromatography on the different methods for the determination of lipophilicity: An essential analytical tool in medicinal chemistry. Chemosensors 2022, 10, 340. [Google Scholar] [CrossRef]
  24. Jeleń, M.; Pluta, K.; Latocha, M.; Morak-Młodawska, B.; Suwińska, K.; Kuśmierz, D. Evaluation of angularly condensed diquinothiazines as potential anticancer agents. Bioorg. Chem. 2019, 87, 810–820. [Google Scholar] [CrossRef]
  25. Skonieczna, M.; Kasprzycka, A.; Jeleń, M.; Morak-Młodawska, B. Tri- and pentacyclic azaphenothiazine as pro-apoptotic agents in lung carcinoma with a protective potential to healthy cell lines. Molecules 2022, 27, 5255. [Google Scholar] [CrossRef]
  26. Available online: http://swissadme.ch (accessed on 10 January 2024).
  27. Available online: https://biosig.unimelb.edu.au/pkcsm/prediction (accessed on 5 January 2024).
  28. Available online: https://www.molinspiration.com/services/logp.html (accessed on 5 January 2024).
  29. ChemDraw: ChemDraw Ultra, Version 22.2.0; CambridgeSoft, PerkinElmer Informatics: Austin, TX, USA, 2022.
  30. Bodor, N.; Gabanyi, Z.; Wong, C.K. A new method for the estimation of partition coefficient. J. Am. Chem. Soc. 1989, 111, 3783–3786. [Google Scholar] [CrossRef]
  31. Mannhold, R.; Cruciani, G.; Dross, K.; Rekker, R. Multivariate analysis of experimental and computational descriptors of molecular lipophilicity. J. Comput. Mol. Des. 1998, 12, 573–581. [Google Scholar] [CrossRef]
  32. Brooke, D.; Dobbs, J.; Williams, N. Octanol:water partition coefficients (P): Measurement, estimation, and interpretation, particularly for chemicals with P > 105. Ecotoxicol. Environ. Saf. 1986, 11, 251–260. [Google Scholar] [CrossRef]
  33. Azman, M.; Sabri, A.H.; Anjani, Q.K.; Mustaffa, M.F.; Hamid, K.A. Intestinal absorption study: Challenges and absorption enhancement strategies in improving oral drug delivery. Pharmaceuticals 2022, 15, 975. [Google Scholar] [CrossRef]
  34. Tsakovska, I.; Pajeva, I.; Al Sharif, M.; Alov, P.; Fioravanzo, E.; Kovarich, S.; Worth, A.P.; Richarz, A.N.; Yang, C.; Mostrag-Szlichtyng, A.; et al. Quantitative structure-skin permeability relationships. Toxicology 2017, 387, 27–42. [Google Scholar] [CrossRef]
  35. Gombar, V.K.; Hall, S.D. Quantitative structure−activity relationship models of clinical pharmacokinetics: Clearance and volume of distribution. J. Chem. Inf. Model. 2013, 53, 948–957. [Google Scholar] [CrossRef]
  36. Poulin, P.; Theil, F.P. Prediction of pharmacokinetics prior to in vivo studies. 1. Mechanism-based prediction of volume of distribution. J. Pharm. Sci. 2002, 91, 129–156. [Google Scholar] [CrossRef]
  37. Mulpuru, V.; Mishra, N. In silico prediction of fraction unbound in human plasma from chemical fingerprint using automated machine learning. ACS Omega 2021, 10, 6791–6797. [Google Scholar] [CrossRef]
  38. Sweeney, M.D.; Zhao, Z.; Montagne, A.; Nelson, A.R.; Zlokovic, B.V. Blood-Brain Barrier: From physiology to disease and back. Physiol. Rev. 2019, 99, 21–78. [Google Scholar] [CrossRef]
  39. Horde, G.W.; Gupta, V. Drug Clearance. 20 June 2023. In StatPearls [Internet]; StatPearls Publishing: Treasure Island, FL, USA, 2024. [Google Scholar] [PubMed]
  40. Das, A.; Kumar, S.; Persoons, L.; Daelemans, D.; Schols, D.; Alici, H.; Tahtaci, H.; Karki, S.S. Synthesis, in silico ADME, molecular docking and in vitro cytotoxicity evaluation of stilbene linked 1,2,3-triazoles. Heliyon 2021, 7, e05893. [Google Scholar] [CrossRef]
  41. Goetz, G.H.; Shalaeva, M. Leveraging chromatography based physicochemical properties for efficient drug design. ADMET DMPK 2018, 6, 85. [Google Scholar] [CrossRef]
  42. Kaliszan, R. QSRR: Quantitative Structure-(Chromatographic) Retention Relationships. Chem. Rev. 2007, 107, 3212–3246. [Google Scholar] [CrossRef]
  43. Kempińska, D.; Chmiel, T.; Kot-Wasik, A.; Mróz, A.; Mazerska, Z.; Namie’snik, J. State of the art and prospects of methods for determination of lipophilicity of chemical compounds. Trends Anal. Chem. 2019, 113, 54–73. [Google Scholar] [CrossRef]
  44. Di, L.; Kerns, E.H. Drug-Like Properties; Academic Press: New York, NY, USA, 2008. [Google Scholar]
  45. Petrauskas, A.A.; Kolovanov, E.A. ACD/Log P method description. Perspect. Drug Discov. Des. 2000, 19, 99–116. [Google Scholar] [CrossRef]
  46. Lipinski, C.A.; Lombardo, F.; Dominy, B.W.; Feeney, P.J. Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv. Drug Deliv. Rev. 1997, 23, 3–25. [Google Scholar] [CrossRef]
  47. Law, V.; Knox, C.; Djoumbou, Y.; Jewison, T.; Guo, A.C.; Liu, Y.; Maciejewski, A.; Arndt, D.; Wilson, M.; Neveu, V.; et al. DrugBank 4.0: Shedding new light on drug metabolism. Nucleic Acids Res. 2014, 42, D1091–D1097. [Google Scholar] [CrossRef]
  48. Ghose, A.K.; Viswanadhan, V.N.; Wendoloski, J.J. A knowledge-based approach in designing combinatorial or medicinal chemistry libraries for drug discovery. A qualitative and quantitative characterization of known drug databases. J. Comb. Chem. 1999, 1, 55–68. [Google Scholar] [CrossRef]
  49. Veber, D.F.; Johnson, S.R.; Cheng, H.-Y.; Smith, B.R.; Ward, K.W.; Kopple, K.D. Molecular properties that influence the oral bioavailability of drug candidates. J. Med. Chem. 2002, 45, 2615–2623. [Google Scholar] [CrossRef]
  50. Egan, W.J.; Merz, K.M.; Baldwin, J.J. Prediction of drug absorption using multivariate statistics. J. Med. Chem. 2000, 43, 3867–3877. [Google Scholar] [CrossRef]
  51. Martin, Y.C. A Bioavailability Score. J. Med. Chem. 2005, 48, 3164–3170. [Google Scholar] [CrossRef]
  52. Daina, A.; Michielin, O.; Zoete, V. SwissADME: A free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules. Sci. Rep. 2017, 7, 42717. [Google Scholar] [CrossRef]
  53. Lobo, S. Is there enough focus on lipophilicity in drug discovery? Expert Opin. Drug Discov. 2020, 15, 261–263. [Google Scholar] [CrossRef]
  54. Pantaleao, S.Q.; Fernandes, P.O.; Goncalves, J.E.; Maltarollo, V.G.; Honorio, K.M. Recent advances in the prediction of pharmacokinetics properties in drug design studies: A review. Chem. Med. Chem. 2022, 17, e202100542. [Google Scholar] [CrossRef]
  55. Daina, A.; Zoete, V. A BOILED-Egg to predict gastrointestinal absorption and brain penetration of small molecules. Chem. Med. Chem. 2016, 11, 1117–1121. [Google Scholar] [CrossRef]
  56. Šegan, S.; Penjišević, J.; Šukalović, V.; Andrić, D.; Milojković-Opsenica, D.; Kostić-Rajačić, S. Investigation of lipophilicity and pharmacokinetic properties of 2-(methoxy)phenylpiperazine dopamine D2 ligands. J. Chromatogr. B Anal. Technol. Biomed. Life Sci. 2019, 1124, 146–153. [Google Scholar] [CrossRef]
  57. Moriguchi, I.; Shuichi, H.; Liu, Q.; Nakagome, I.; Matsushita, Y. Simple method of calculating octanol/water partition coefficient. Chem. Pharm. Bull. 1992, 40, 127–130. [Google Scholar] [CrossRef]
  58. Moriguchi, I.; Shuichi, H.; Nakagome, I.; Hirano, H. Comparison of reliability of log P values for drugs calculated by several methods. Chem. Pharm. Bull. 1994, 42, 976–978. [Google Scholar] [CrossRef]
  59. Stanisz, A. Przystępny Kurs Statystyki z Zastosowaniem STATISTICA PL na Przykładach Medycyny; Analizy Wielowymiarowe; Stat Soft Polska: Kraków, Poland, 2007; Volume 3. [Google Scholar]
Figure 1. Structures of angularly condensed N-dialkyloaminoalkylodiquinothiazines 115.
Figure 1. Structures of angularly condensed N-dialkyloaminoalkylodiquinothiazines 115.
Pharmaceuticals 17 00725 g001
Figure 2. Comparison of the theoretical lipophilicity parameters (the logP values) of the examined angularly condensed N-dialkyloaminoalkyldiquinothiazines 115.
Figure 2. Comparison of the theoretical lipophilicity parameters (the logP values) of the examined angularly condensed N-dialkyloaminoalkyldiquinothiazines 115.
Pharmaceuticals 17 00725 g002
Figure 3. Dendrogram of partition coefficients of the examined angularly condensed N-dialkyloaminoalkyldiquinothiazines 115.
Figure 3. Dendrogram of partition coefficients of the examined angularly condensed N-dialkyloaminoalkyldiquinothiazines 115.
Pharmaceuticals 17 00725 g003
Figure 4. Dendrogram of the examined angularly condensed N-dialkyloaminoalkyldiquinothiazines 115 based on their partition coefficients.
Figure 4. Dendrogram of the examined angularly condensed N-dialkyloaminoalkyldiquinothiazines 115 based on their partition coefficients.
Pharmaceuticals 17 00725 g004
Figure 5. Dendrogram of theoretical (logPaverge) and chromatographic parameters (RM0, C0) of examined angularly condensed N-dialkyloaminoalkyldiquinothiazines 115.
Figure 5. Dendrogram of theoretical (logPaverge) and chromatographic parameters (RM0, C0) of examined angularly condensed N-dialkyloaminoalkyldiquinothiazines 115.
Pharmaceuticals 17 00725 g005
Figure 6. The boiled egg representation of the intestinal absorption and permeation through the blood–brain barrier for diquinothiazines 115.
Figure 6. The boiled egg representation of the intestinal absorption and permeation through the blood–brain barrier for diquinothiazines 115.
Pharmaceuticals 17 00725 g006
Table 1. The computed lipophilicity parameters (logPcalcd) for diquinothiazines 115 using the following Internet databases: SwissADME [26], pkCSM a [27], Molinspiration b [28], and ChemDraw c [29].
Table 1. The computed lipophilicity parameters (logPcalcd) for diquinothiazines 115 using the following Internet databases: SwissADME [26], pkCSM a [27], Molinspiration b [28], and ChemDraw c [29].
No. of
Compound
iLOGPXLOGP3WLOGPMLOGPSILICOS-ITLogP amilogP blogP clogPaverage
14.215.395.354.244.665.735.615.175.05 (±0.60)
23.985.014.964.034.275.345.134.604.67 (±0.52)
34.095.144.724.244.515.485.264.814.91 (±0.61)
44.175.505.114.454.745.875.775.225.10 (±0.61)
54.285.935.504.654.846.266.035.315.35 (±0.71)
64.085.395.354.244.665.735.445.175.01 (±0.61)
73.835.014.964.034.275.344.964.604.63 (±0.53)
84.065.144.724.244.515.485.094.814.76 (±0.48)
94.015.505.114.454.745.875.605.225.06 (±0.63)
104.285.935.504.654.846.265.865.315.33 (±0.69)
113.995.395.354.244.665.735.145.174.96 (±0.60)
123.675.014.964.034.275.344.664.604.57 (±0.55)
134.125.144.724.244.515.484.794.814.73 (±0.45)
143.775.505.114.454.745.875.305.225.00 (±0.66)
153.975.935.504.654.846.265.555.315.25 (±0.74)
Table 2. Data for linear correlation (RM = RM0 + bC) for compounds 115.
Table 2. Data for linear correlation (RM = RM0 + bC) for compounds 115.
No. of CompoundRM0brC0
13.45−4.770.99320.7233
23.01−4.090.99610.7359
33.05−4.230.99160.7210
43.48−4.890.99530.7117
53.64−5.140.99800.7082
63.68−4.850.99510.7588
73.31−4.340.99380.7627
83.48−4.600.99410.7565
93.69−4.880.99620.7562
103.83−5.040.99620.7599
113.60−4.750.99270.7579
123.17−4.250.99350.7458
133.23−4.210.99070.7672
143.61−4.690.99520.7697
153.82−5.010.99420.7625
Table 3. The RM0 and logPlit values and b (slope) and r (correlation coefficient) values of the equation RM = RM0+ bC for standards IV.
Table 3. The RM0 and logPlit values and b (slope) and r (correlation coefficient) values of the equation RM = RM0+ bC for standards IV.
Lipophilicity
Parameters
Standards
IIIIIIIVV
logPlit1.21 [30]1.87 [31]3.18 [31]4.45 [31]6.38 [32]
RM00.781.162.513.334.69
-b0.01620.02470.03280.04120.0564
r0.99230.99370.99710.99820.9977
logPTLC1.211.703.434.496.24
Table 4. The experimental lipophilicity parameters (logPTLC values) for compounds 115.
Table 4. The experimental lipophilicity parameters (logPTLC values) for compounds 115.
No. of CompoundlogPTLCNo. of CompoundlogPTLCNo. of CompoundlogPTLC
14.6464.94114.84
24.0874.46124.28
34.1384.68134.36
44.6894.95144.85
54.89105.13155.12
Table 5. The molecular descriptors for compounds 115.
Table 5. The molecular descriptors for compounds 115.
No. of
Compound
Molar Mass (M) [g/mol]Molar Volume (VM) [cm3]Molar Refractivity (RefM) [cm3/mol]Surface Area [Å]
1400.55322.1121.077175.12
2386.52305.6116.446168.75
3398.54304.6119.121174.11
4412.56322.4123.722180.48
5426.59342.3128.258186.84
6400.55322.1121.077175.12
7386.52305.6116.446168.75
8398.54304.6119.121174.11
9412.56322.4123.722180.48
10426.59342.3128.258186.84
11400.55322.1121.077175.12
12386.52305.6116.446168.75
13398.54304.6119.121174.11
14412.56322.4123.722180.48
15426.59342.3128.258186.84
Table 6. The drug-likeness and ADME properties predicted by in silico studies using SwissADME.
Table 6. The drug-likeness and ADME properties predicted by in silico studies using SwissADME.
Predicted ParameterCompound No.
1, 6, 112, 7, 123, 8, 134, 9, 145, 10, 15
Physicochemical Properties
Num. heavy atoms2928293031
Num. arom. heavy atoms2020202020
Hydrogen bond acceptors33333
Hydrogen bond donors00000
Number of rotatable bonds54333
Topological polar surface area [Å2]57.5657.5657.5657.5657.56
Drug-Likeness Prediction
Rule of Lipinski+++++
Rule of Ghose+++−(MR>130)−(MR>130)
Rule of Veber+++++
Rule of Egan+++++
Rule of Muegge−(XLOGP3>5)−(XLOGP3>5)−(XLOGP3>5)−(XLOGP3>5)−(XLOGP3>5)
Bioavailability
Bioactivity score0.550.550.550.550.55
Table 7. The absorption descriptors for compounds 115.
Table 7. The absorption descriptors for compounds 115.
No. of CompoundWater Solubility
[log mol/L]
Caco-2 Permeability [log Papp in 10−6 cm/s]Intestinal Absorption [% Absorbed]Skin Permeability [log Kp]
1−5.8711.00392.241−2.697
2−5.7831.02892.931−2.701
3−4.5511.14592.725−2.733
4−4.6601.14392.337−2.733
5−5.8201.02192.309−2.711
6−5.7931.03295.220−2.694
7−5.7371.05796.401−2.700
8−4.5171.22694.349−2.723
9−4.6521.22493.960−2.723
10−5.8291.04995.779−2.712
11−5.3731.04395.329−2.691
12−5.2631.06896.510−2.697
13−3.9651.21193.469−2.744
14−4.0611.20993.080−2.743
15−5.3291.06095.888−2.709
Table 8. The distribution descriptors for compounds 115.
Table 8. The distribution descriptors for compounds 115.
No. of
Compound
VDss
[log L/kg]
Unbound Fraction [Fu]BBB Permeability [log BB]CNS Permeability [log PS]
10.9860.2590.478−1.464
20.8640.2530.483−1.402
30.8680.1930.424−1.394
40.9230.1890.437−1.381
51.0410.2550.535−1.422
61.1780.2680.537−1.493
71.0620.2620.484−1.478
81.2690.2060.559−1.318
91.3280.2000.572−1.304
101.2440.2570.536−1.498
111.3320.2670.387−1.475
121.2000.2610.357−1.483
131.1700.2000.425−1.370
141.2220.1960.438−1.356
151.3680.2620.410−1.503
Table 9. The excretion and toxicity for compounds 115.
Table 9. The excretion and toxicity for compounds 115.
No. of
Compound
Total Clearance [log ml/min/kg]Max. Tolerated Dose
[log mg/kg/day]
Oral Rat Acute Toxicity
[mol/kg]
Oral Rat Chronic Toxicity
[log mg/kg bw/day]
Tetrahymena pyriformis Toxicity [log μg/L] Minnow Toxicity [log mM]
10.6720.6722.2770.4490.2991.983
20.5260.6472.2720.4900.3001.944
30.8240.5372.7340.8270.291−0.279
40.7790.5502.7560.8480.291−0.396
50.5990.5682.3630.5420.2911.736
60.7140.2932.3610.7730.3270.423
70.5840.2582.3120.8070.3310.597
80.8320.2113.1571.0740.300−0.872
90.7870.2283.1841.0960.299−0.989
100.6570.1882.4130.6060.3040.389
110.7070.5942.5130.5900.2990.542
120.5620.5572.4850.6100.3000.460
130.8260.5432.9491.0200.292−0.848
140.7810.5572.9721.0410.291−0.965
150.6340.4832.5880.6630.2910.252
Table 10. Linear correlations between the partition coefficients and the chromatographic parameters of lipophilicity of the tested compounds 115 (p < 0.05).
Table 10. Linear correlations between the partition coefficients and the chromatographic parameters of lipophilicity of the tested compounds 115 (p < 0.05).
CompoundsLipophilicity ParameterEquationr
15iLOGPiLOGP = 0.3877 RM0 + 2.85660.9403
115XLOGP3XLOGP3 = 1.0578 RM0 + 1.72330.8426
15 XLOGP3 = 1.1972RM0 + 1.41220.9384
610 XLOGP3 = 1.6538 RM0 − 0.55650.9427
1115 XLOGP3 = 1.2383 RM0 + 1.07730.9577
115WLOGPWLOGP = 0.8350 RM0 + 2.23040.7683
15 WLOGP = 0.9794 RM0 + 1.87050.8868
1115 WLOGP = 0.9985 RM0 + 1.64730.8919
115MLOGPMLOGP = 0.6500 RM0 + 2.06640.7840
610 MLOGP = 1.0545 RM0 + 0.52790.9101
1115 MLOGP = 0.7511 RM0 + 1.70350.8795
115SILCOS-ITSILCOS-IT = 0.6630 RM0 + 2.30340.8486
15 SILCOS-IT = 0.7442 RM0 + 2.12870.9374
610 SILCOS-IT = 1.0799 RM0 + 0.71840.9892
1115 SILCOS-IT = 0.7588 RM0 + 1.95880.9432
115LogPLogP = 1.0697 RM0 + 2.0243 0.8484
15 LogP = 1.2120 RM0 + 1.70500.9460
610 LogP = 1.6739 RM0 − 0.28660.9500
1115 LogP = 1.2499 RM0 + 1.37880.9624
610 logP = 1.4693 RM0 − 0.26440.9877
1115 logP = 1.0652 RM0 + 1.30860.9714
115milogPmilogP = 0.9755 RM0 + 1.96100.6499
15 milogP = 1.2903 RM0 + 1.26840.9813
610 milogP = 1.7585 RM0 − 0.93710.9725
1115 milogP = 1.3037 RM0 + 0.54310.9867
115logPaveragelogPaverage = 0.7369 RM0 + 2.40170.7854
15 logPaverage = 0.8307 RM0 + 2.25300.9292
610 logPaverage = 1.3129 RM0 + 0.23430.9791
1115 logPaverage = 0.9300 RM0 + 1.66010.9819
Table 11. Linear correlations between experimentally determined RM0 and the molecular descriptors’ values (p < 0.05).
Table 11. Linear correlations between experimentally determined RM0 and the molecular descriptors’ values (p < 0.05).
CompoundsMolecular DescriptorEquationr
115Molar massM = 42.779 RM0 + 256.5080.8000
15 M = 47.859 RM0 + 245.7730.8807
610 M = 68.869 RM0 + 157.1610.9216
1115 M = 49.570 RM0 + 232.1520.8999
115Molar volumeVM = 45.782 RM0 + 160.7110.8444
15 VM = 52.381 RM0 + 145.1800.9518
610 VM = 69.514 RM0 + 69.2880.9182
1115 VM = 53.970 RM0 + 131.2620.9674
115Molar refractivityRefM = 13.202 RM0 + 75.9130.8311
15 RefM = 14.871 RM0 + 72.2630.9213
610 RefM = 20.942 RM0 + 46.3750.9434
1115 RefM = 15.363 RM0 + 68.1700.9389
115Surface areaSurface area = 19.960 RM0 + 109.8820.8017
15 Surface area = 21.667 RM0 + 105.0000.8831
610 Surface area = 31.142 RM0 + 65.0100.9229
1115 Surface area = 22.437 RM0 + 98.8430.9021
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Klimoszek, D.; Jeleń, M.; Dołowy, M.; Morak-Młodawska, B. Study of the Lipophilicity and ADMET Parameters of New Anticancer Diquinothiazines with Pharmacophore Substituents. Pharmaceuticals 2024, 17, 725. https://doi.org/10.3390/ph17060725

AMA Style

Klimoszek D, Jeleń M, Dołowy M, Morak-Młodawska B. Study of the Lipophilicity and ADMET Parameters of New Anticancer Diquinothiazines with Pharmacophore Substituents. Pharmaceuticals. 2024; 17(6):725. https://doi.org/10.3390/ph17060725

Chicago/Turabian Style

Klimoszek, Daria, Małgorzata Jeleń, Małgorzata Dołowy, and Beata Morak-Młodawska. 2024. "Study of the Lipophilicity and ADMET Parameters of New Anticancer Diquinothiazines with Pharmacophore Substituents" Pharmaceuticals 17, no. 6: 725. https://doi.org/10.3390/ph17060725

APA Style

Klimoszek, D., Jeleń, M., Dołowy, M., & Morak-Młodawska, B. (2024). Study of the Lipophilicity and ADMET Parameters of New Anticancer Diquinothiazines with Pharmacophore Substituents. Pharmaceuticals, 17(6), 725. https://doi.org/10.3390/ph17060725

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