- Research article
- Open Access
Assessment of lipophilicity of newly synthesized celecoxib analogues using reversed-phase HPLC
BMC Chemistryvolume 13, Article number: 84 (2019)
Lipophilicity is a physicochemical property of an essential importance in medicinal chemistry; therefore, fast and reliable measurement of lipophilicity will affect greatly the drug discovery process.
A series of N-benzenesulfonamide-1H-pyrazoles, oximes and hydrazones as celecoxib analogues was investigated with regard to their retention behavior using reversed-phase high performance liquid chromatography (RP-HPLC). The mobile phases employed for this study consist of a mixture of water and methanol in different proportions. In addition, the stationary phase utilized for this separation was C18 silanized silica gel and using 200 nm as a detection wavelength. The retention behavior of the investigated compounds was determined based on practical determination of log k at different concentrations of methanol (as an organic modifier) in the mobile phase ranging from 60 to 80%. It was observed that the retention of these compounds (expressed as log k) decreased in a linear manner with increasing the concentration of methanol. High correlation coefficients (more than 0.90 in most cases) were obtained for the relationship between the volume fraction of the organic solvent and the retention values represented as log kw. A comparative evaluation was carried out between chromatographically-obtained lipophilicity parameters (represented as lipophilicity chromatographic index log kw or isocratic chromatographic hydrophobicity index, \(\varphi\) 0) and the computationally calculated log P values (miLogP, ALOGP, ACD/Labs and ALOGPs).
It was found that a good correlation exists between the experimental and computed log P values. In the future, these results can give a deep insight about the anti-inflammatory and analgesic activity of the newly synthesized compounds.
Lipophilicity is a paramount descriptor which represents an essential part in design of new medications with required biological action . It is also used in quantitative structure activity relationship (QSAR) investigations. The IUPAC (International Union of Pure and Applied Chemistry) defined it as “a physicochemical property which describes a partitioning equilibrium of solute molecules between water and an immiscible organic solvent” . This property is one of the leading considerations as it has an important effect on pharmacodynamic and pharmacokinetic properties of drugs . It suggests how the ADME (Absorption, Distribution, Metabolism, and Elimination) features of medications will proceed, in addition to its influence on their toxicological profile . For its vital role on drug discovery and design process, the estimation of lipophilicity and how to regulate it has turned out to be routine procedure in drug development field . Lipophilicity is also important in demonstrating the destiny of a compound in the environment, where it influences the bioavailability and bio-concentration in the food chain via sorption from water, and soil or dregs, which makes it an essential issue in danger evaluation and controlling of risky resources .
Lipophilicity is a fundamental molecular property which is defined as the logarithm of the octanol–water partition coefficient (log POW), which is, in turn, expressed as the non-ionized drug spread ratio between the two phases of octanol and water . log P measurements by different experimental protocols are illustrated in detail by the Organization for Economic Cooperation and Development (OECD) guidelines as test No. 107, which is Shake Flask method , and Test No. 123, which is the Slow Stirring method . These traditional procedures are time-consuming, limited to extremely pure compound and need specific reagent to be performed. Hence, they have recently been replaced by more adaptable simpler methodologies which are chromatographic techniques (the OECD Test No. 117 ) that provide coherent results in the same log P range.
In the last few years, a great interest has been given for RP-HPLC as a tool for lipophilicity estimation and for characterization of pharmacological activity of molecules [11,12,13]. This is attributed to the close relationship between the retention performance of molecules in reversed phase chromatographic system and its lipophilicity . The major advantages in using HPLC as a superior strategy are: smaller sample amount, high sensitivity to impurities, reproducibility and accuracy in addition to short determination time. The majority of the chromatographic systems used for measuring the lipophilicity rely on hydrophobic reversed-phase silica-based stationary phase (C8 or C18) .
In parallel to HPLC methods, computational methodologies have also been employed. Their wide use is based on considering the industrial requests as being simple, low cost, fast, and consistent approach to deliver information for fast screening of compounds under focus .
Celecoxib,4-[5-(4-methylphenyl)-3-(trifluoromethyl)pyrazol-1yl]benzenesulfonamide  is a selective inhibitor of cyclooxygenase-2 (COX-2) enzyme that is prescribed as a non-steroidal anti-inflammatory, analgesic, and antipyretic drug . Being a selective cyclooxygenase-2 inhibitor, celecoxib seems to show fewer side effects than non-selective non-steroidal anti-inflammatory agents . Unlike celecoxib, these non-selective agents act as inhibitors of both COX-1 and COX-2 enzymes which increases the possibility of inducing gastric ulcers, gastric bleeding and suppressed renal functions . On the other hand, celecoxib has been alerted by FDA for its adverse side effects on cardiovascular system, and that represents the main motivation for discovering novel compounds with COX-2 inhibitory activity .
The principal aim of the present investigation is to develop an appropriate and proficient strategy for estimating the lipophilicity of a series of recently discovered N-benzenesulfonamide-1H-pyrazoles, oximes and hydrazones as celecoxib analogues  by RP-HPLC method and comparing these results with those obtained from computational programs. All the compounds in this series are weakly basic in nature (pKa of celecoxib is 11.1) . The investigated compounds are named 11a, 11b, 11c, 11d, 16a, 16b, 17a, and 17b in addition to celecoxib. The chemical structures of compounds under investigation are depicted in Fig. 1. Additional file 1 included H-NMR for the newly synthesized compounds.
Reagents and chemicals
Celecoxib of purity 99.44% was kindly provided by Amoun Pharmaceutical Co. (Cairo, Egypt). The synthesis and characterization of the studied compounds (celecoxib analogues) were previously reported . HPLC grade methanol was purchased from Sigma-Aldrich (Seelze, Germany). Double-distilled water was utilized throughout the study and was filtered through a Millipore membrane filter of bore size 0.45 μm (Additional file 1).
Instrumentation and chromatographic conditions
For chromatographic analysis, stock solutions (200 μg/mL) of each of the investigated samples were prepared in methanol. HPLC measurements were performed utilizing Shimadzu Prominence HPLC system (Shimadzu Corporation, Japan) consisting of LC-20 AD pump, DGU-20 A5 degasser, CBM-20A interface, and SPD-20A UV–Vis detector with a 20 µL injection loop. The column utilized for studying the chromatographic behavior of the analytes was the most frequently used for lipophilicity investigations: Waters symmetry® C18 column [250 mm × 4.6 mm (i.d.), 5-μm particle size]. The flow rate was 1.2 mL/min at room temperature. The retention behaviour of the analytes was investigated as a function of mobile phase composition ranging from 60–80% methanol and 40–20% water. The HPLC analyses were carried out at room temperature under isocratic conditions. Membrane filters (0.45 μm) from Sartorius (Goettingen, Germany) were used for filtration of samples. The absorbance of the analytes during a chromatographic run was performed at 200 nm. Each experiment was run in duplicate. Chromatographic data were collected using LabSolutions CS software. H1NMR spectra were recorded using Bruker NMR spectrometer (Bruker GmbH, Germany),
Computational methods were performed using CS ChemDraw Ultra software (Cambridge Soft Corporation, Cambridge, MA, USA) running under Windows 7 operating system.
Computational lipophilicity (clog P) was calculated by the Advanced Chemistry Development ACD/Labs online service (https://ilab.acdlabs.com/iLab2/), Molinspiration online service (miLOGP, performed online at http://www.molinspiration.com/) and ALOGPS 2.1 online service (ALOGPS, performed online at http://www.vcclab.org/lab/alogps/).
Results and discussion
It was shown from previous reports that the chromatographic systems using methanol: water as mobile phases were better for log k modeling than those containing acetonitrile as an organic modifier . Moreover, it was found that the most preferred packing material for reversed-phase columns in lipophilicity assessment is C18 silanized silica gel. Chromatographic techniques including HPLC are based on the determination of the retention parameters . The indices of lipophilicity measured by HPLC are obtained from the logarithm of the retention factor (log k), which is calculated by the following equation:
where; tR is the retention time of the solute and t0 is the retention time of unretained species .
log k is proved experimentally to be in a linear relationship to the volume fraction of the organic solvent in the mobile phase (\(\varphi\)) according to the following linear regression equation :
where; log kw is the intercept of the linear regression curve which reflects the lipophilicity chromatographic index and denotes the factor of retention for pure water as eluent, and the slope (S) is extensively related to the strength of solvent or with specific hydrophobic surface area of solutes, and \(\varphi\) is the portion of organic modifier volume.
Estimating log kw directly is often very difficult, and is nearly non-feasible, because it can produce elongated retention time and also extreme peak broadening. Hence, estimating k with various proportions of water-organic solvent combination as mobile phases is favored, then extrapolating the correlation between log k and percentage of organic modifier specifies the value of log k upon utilizing pure water as mobile phase . log k was estimated for each of the new celecoxib analogues at different proportions of methanol. Consequently, to obtain log kw, an extrapolation to 100% of pure water was carried out. Then, linear regression analysis was performed . Linear extrapolation is favored using definite ratios of organic modifier consistent with the lipophilicity level of the solutes .
The retention data were expressed as logarithm of the retention factor (log k) as shown in Table 1. It is clear that the lower the volume fraction of the organic modifier, the lower the elution power of the mobile phase and the higher the value of log k. The highest log k values were obtained when 60% v/v methanol was utilized.
The chromatographic data obtained for the nine investigated compounds together with the results of linear regression analysis are listed in Table 2 and represented in Fig. 2. For all examined drugs, high values of correlation coefficients were achieved (r > 0.9098), with small values of standard error of estimate, which proves the high significance of Eq. (2) for determination of lipophilicity. S is related to the specific hydrophobic surface area of solutes and it relies on the solute and the chromatographic system. A good correlation exists between log kw and slope with r ≈ − 0.9716 which suggests a similar chromatographic retention mechanism for this homologous series of studied compounds. This relationship can be represented by the following equation:
(r = 0.9716, R2 = 0.9440, F = 118.12, p value < 0.001, Syx = 0.2822, n = 9), R2 is the determination coefficient, F is the value of test F- Snedecora, Syx is standard error of estimate.
Another substitute to the index of lipophilicity is the chromatographic index of hydrophobicity \(\varphi\)0 which was first given by Valkò et al. . It is calculated simply by the equation: \(\varphi\)0 = − log kw/S. \(\upvarphi_0\) denotes the volume fraction of the organic modifier in the mobile phase at which equal partitioning of the solute between the mobile and stationary phases is obtained, that is, k = 1, log k = 0 . The value of \(\varphi\)0 is specific for each compound; it relies solely on the pH, organic modifier and temperature. It does not rely on the column type and length, flow rate and the composition of the mobile phase with which the retention determination are made. Moreover, \(\varphi\)0 can be accurately estimated as it possesses a solid physical meaning, that is, the concentration of the organic modifier in the mobile phase where the retention time is accurately double the dead time. This is considered a distinguishable advantage over extrapolating of log k values to pure water mobile phase . \(\varphi\)0 values calculated for the studied compounds are listed in Table 2.
Eight compounds of the N-benzenesulfonamide-1H-pyrazoles, oximes and hydrazones series were investigated in this study as celecoxib analogues. They possess similar structural features with different substituents that significantly impact their physicochemical, pharmacokinetic, and pharmacodynamic properties. RP-HPLC was utilized to estimate their lipophilicity (conveyed as log kw or \(\varphi\)0). In this study, C18 column was utilized as the stationary phase which is consisting of silica gel amended by a lipophilic C18 hydrocarbon chain while the mobile phase was more polar than the stationary phase. It is well-known that polar molecules have lower retention times and hence lower retention factors than non-polar molecules which are strongly attracted to the stationary phase hydrocarbon groups as explained by the Van der Waals forces .
For the first group of the studied compounds (aryl sulfones), in relation to the retention data, it can be observed that 11a is more polar than 11b. The methyl group substituent in 11b reduces its polarity. It is also noted that 11c is more polar than 11d due to the methyl substituent in the aryl moiety of the sulfone group in 11d. These observations coincide very well with the data obtained (Tables 1 and 2). Similarly, in the second and third groups of the studied compounds (oximes and hydrazones, respectively), due to the methyl substituent, compound 16a is more polar than 16b and compound 17a is more polar than compound 17b.
The linear regression parameters between log k values and methanol volume fraction are abridged in Table 2.
Correlation between log Kw and φ0 and computed log P values
The following step in this work is to compare chromatographically-obtained lipophilicity parameters; log kw or \(\varphi\)0 of the studied compounds with clog P calculated using different computational methods (that utilize different theoretical approaches, principles of which were described by Mornar et al. ) such as ACD/Labs, miLOGP, ALOGPS, and APLOGP. A very good correlation exists between log kw and computed log P values (r in all cases ≥ 0.8283) as represented in Table 3, which confirms the reliability of HPLC technique as a tool for lipophilicity estimation. Although a fair correlation exists between \(\varphi\)0 and miLogP, a good correlation still presents between \(\varphi\)0 and the other computed log P values (Table 3), which is a further proof for the validity of this approach for determination of lipophilicity. All the investigated compounds purity was confirmed by the H1NMR spectrometer (Table 4).
For further investigation of the differences among the computational and experimental lipophilicity, the lipophilicity values were arranged in a matrix of dimensions nine (compounds) × six (2 computational and four experimental lipophilicity (log kw or \(\varphi\)0)), then principal component analysis PCA was performed on the whole matrix.
The first principal component PC1 accounts for 86.7% of the data variation, whereas PC2 explains 8.1% of the data variation, respectively (94.8% total). As given in Fig. 3, the variation in the lipophilicity of the investigated analogues can be mainly explained by PC1, while PC2 explains differences among chromatographic and computational techniques, which strongly justify the utilization of chromatographic techniques for lipophilicity estimation as a promising substitute to experimental methods based on Shake-Flask method.
In this work, determination of lipophilicity of celecoxib and its newly synthesized analogues was carried out by RP-HPLC using a C18 silanized silica gel stationary phase. This experimental study proved to be simple, reliable and does not need sophisticated or complicated methodologies which show the importance of this method for pharmaceutical applications. Lipophilicity is an essential physicochemical property that affects pharmacokinetics and pharmacodynamics as well as toxicity of drug molecules. Moreover, it is well known that this essential property is important for the in vivo distribution of organic compounds by influencing their solubility, oral absorption, cell uptake, blood–brain penetration and metabolism and to rationalizing a number of biological events as membrane penetration and permeability. Under suitable isocratic chromatographic conditions, extrapolated retention factors are in good correlation with the computationally calculated log P obtained using different software programs. Hence, RP-HPLC is an easily applicable technique for determination of lipophilicity and the chromatographic retention.
Availability of data and materials
The data material are available whenever its needed from authors.
reversed-phase high performance liquid chromatography
- log P OW :
octanol–water partition coefficient
- \(\varphi\) 0 :
the chromatographic index of hydrophobicity
Marciniec K, Bafeltowska J, Maślankiewicz MJ, Buszman E, Boryczka S (2016) Determination of the lipophilicity of quinolinesulfonamides by reversed-phase HPLC and theoretical calculations. J Liq Chromatogr Relat Technol 39(15):702–709
Duffus JH, Nordberg M, Templeton DM (2007) Glossary of terms used in toxicology, (IUPAC Recommendations 2007). Pure Appl Chem 79(7):1153–1344
Ayouni L, Cazorla G, Chaillou D, Herbreteau B, Rudaz S, Lantéri P et al (2005) Fast determination of lipophilicity by HPLC. Chromatographia 62(5):251–255
Giaginis C, Tsantili-Kakoulidou A (2007) Current state of the art in HPLC methodology for lipophilicity assessment of basic drugs. A review. J Liq Chromatogr Relat Technol 31(1):79–96
Santos Á, Soares JX, Cravo S, Tiritan ME, Reis S, Afonso C et al (2018) Lipophilicity assessement in drug discovery: experimental and theoretical methods applied to xanthone derivatives. J Chromatogr B. 1072(Supplement C):182–192
Andrić F, Héberger K (2015) Towards better understanding of lipophilicity: assessment of in silico and chromatographic logP measures for pharmaceutically important compounds by nonparametric rankings. J Pharm Biomed Anal 115(Supplement C):183–191
Lipinski CA, Lombardo F, Dominy BW, Feeney PJ (1997) Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv Drug Deliv Rev 23:3–25
OECD (1995) Test No. 107, Partition coefficient(n-octanol/water), Shake flask method. OECD, Paris
OECD (2005) Test No. 123, Partition coefficient(n-octanol/water), Slow-stirring method. OECD, Paris
OECD (2004) Test No. 117, Partition coefficient(n-octanol/water), HPLC method. OECD, Paris
Hawry AM, Popiołek LP, Hawrył MA, Ryszard HS, Niejedli MA (2015) Chromatographic and calculation methods for analysis of the lipophilicity of newly synthesized thiosemicarbazides and their cyclic analogues 1,2,4-triazol-3-thiones. J Brazil Chem Soc 26(8):1617–1624
Rudraraju AV, Amoyaw PNA, Hubin TJ, Khan MOF (2014) Determination of log P values of new cyclen based antimalarial drug leads using RP-HPLC. Pharmazie 69:655–662
Marciniec K, Bafeltowska J, Maślankiewicz M, Buszman E, Boryczka S (2016) Determination of the lipophilicity of quinoline sulfonamides by reversed-phase HPLC and theoretical calculations. J Liqu Chromatogr Relat Technol 39(1):702–709
Tuzimski T, Sztanke K (2005) 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. JPC J Plan Chromatogr Mod TLC. 18(104):274–281
Moffat ACO, David M, Widdop B, Watts J (2011) Clarke’s analysis of drugs and poisons, 4th edn. Pharmaceutical press, London
Sweetman SC (2014) Martindale: the complete drug reference, 38th edn. Pharmaceutical Press, London
Miles S (2007) Celecoxib.* xPharm: the comprehensive pharmacology reference. Elsevier, New York, pp 1–7
Bennett JS, Daugherty A, Herrington D, Greenland P, Roberts H, Taubert KA (2005) The use of nonsteroidal anti-inflammatory drugs (NSAIDs). A science advisory from the American Heart Association. Circulation 111(13):1713–1716
Abdel-Aziz HA, Al-Rashood KA, ElTahir KEH, Suddek GM (2014) Synthesis of N-benzenesulfonamide-1H-pyrazoles bearing arylsulfonyl moiety: novel celecoxib analogs as potent anti-inflammatory agents. Eur J Med Chem. 80:416–422
Snyder LR, Kirkland JJ, Dolan JW (2011) Introduction to modern liquid chromatography. Wiley, Hoboken
Haber P, Baczek T, Kaliszan R, Snyder LR, Dolan JW, Wehr CT (2000) Computer simulation for the simultaneous optimization of any two variables and any chromatographic procedure. J Chromatogr Sci 38(9):386–392
Sima IA, Kot-Wasik A, Wasik A, Namieśnik J, Sârbu C (2017) Assessment of lipophilicity indices derived from retention behavior of antioxidant compounds in RP-HPLC. Molecules 22(4):550
Volná T, Motyka K, Hlaváč J (2016) Potential of high-performance liquid chromatography for distribution coefficient determination of 3-hydroxyquinolin-4(1H)-one derivatives. Chromatographia 79(17):1153–1163
Valkó K, Slégel P (1993) New chromatographic hydrophobicity index (\(\varphi\) 0) based on the slope and the intercept of the log k′ versus organic phase concentration plot. J Chromatogr A 631(1):49–61
Jevrić LR, Karadžić MŽ, Mandić AI, Podunavac Kuzmanović SO, Kovačević SZ, Nikolić AR et al (2017) Lipophilicity estimation and characterization of selected steroid derivatives of biomedical importance applying RP HPLC. J Pharm Biomed Anal 134:27–35
Mornar A, Damic A, Nigovic B (2011) Pharmacokinetic parameters of statin drugs characterized by reversed phase high-performance liquid chromatography. Anal Lett 44:1009–1020
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.