Abstract

The lipophilicity of ten ruthenium(II)-arene complexes was assessed by reversed-phase thin-layer chromatography (RP-TLC) on octadecyl silica stationary phase. The binary solvent systems composed of water and acetonitrile were used as mobile phase in order to determine chromatographic descriptors for lipophilicity estimation. Octanol-water partition coefficient, , of tested complexes was experimentally determined using twenty-eight standard solutes which were analyzed under the same chromatographic conditions as target substances. In addition, ab initio density functional theory (DFT) computational approach was employed to calculate values from the differences in Gibbs’ free solvation energies of the solute transfer from n-octanol to water. A good overall agreement between DFT calculated and experimentally determined values was established (R2 = 0.8024–0.9658).

1. Introduction

Apart from being important in material science and catalysis, metal ions and their complexes play a significant role in the vital functions of living organisms. Numerous applications of metal-based compounds as both therapeutic and diagnostic agents as well as mineral supplements were studied in scope of relatively young but rapidly developing research discipline named medicinal inorganic chemistry [1]. Metal complexes have a long history of use as medicines with cytostatic, antirheumatic, or anti-inflammatory properties. Further, they have been used in treatment of cardiac and many other diseases, while the anticancer activity of metal complexes, especially of cisplatin related drugs, is of the greatest importance [2].

In addition to the widely used platinum-based chemotherapeutic drugs such as cisplatin or carboplatin, numerous non-platinum-based compounds were investigated as anticancer agents [1, 2]. Among them various ruthenium complexes attracted recent attention. Their anticancer activity, as well as clinical toxicity, is clearly distinct from platinum complexes [35]. Diversity of modes of action that involves both extra- and intracellular processes was achieved through interactions of ruthenium complexes with plasma proteins, extracellular matrix collagens, actins on the cell surface, regulatory enzymes in the plasma membranes or cytoplasm, and DNA in the cell nucleus [6]. In addition, some ruthenium complexes exhibit greater efficacy against cancer metastasis than against primary tumors by modulating adhesion, migration, invasion, proteolytic degradation of extracellular matrix, and new blood vessel formation. Over the last thirty years the interest for various ruthenium complexes, their synthesis and potential anticancer activity is constantly increasing, among them ruthenium(II)-arene compounds are being in the focus of research.

Our previous studies have been focused on syntheses, characterization, and cytotoxic activity of series of ruthenium(II)-arene compounds [79]. It has been shown that the presence of arene ligand is crucial. It is involved in several steps: it stabilizes the +2 ruthenium oxidation state, affects the cell uptake, influences the interactions with potential intracellular targets, and provides satisfactory lipophilicity needed to cross the cell membrane.

Since the lipophilicity is one of the major parameters affecting important biological processes that follow drug intake such as adsorption, passage through membranes, drug-receptor interactions, metabolism, and toxicity of molecules [10], we have decided to put an emphasis on the determination and estimation of lipophilicity of several ruthenium(II)-arene complexes in the scope of the present work.

Octanol-water partition coefficient is a widely accepted measure of lipophilicity and can be determined for various compounds including metal complexes in several ways. According to original “shake-flask” method, is defined as a concentration ratio of compound distributed between n-octanol and aqueous phase. However, this time-consuming method with many experimental limitations has been replaced with chromatographic methods such as reversed-phase thin-layer chromatography (RP-TLC) [11], high-performance liquid chromatography (HPLC) in both, isocratic and gradient elution modes [12, 13], microemulsion electrokinetic chromatography (MEEKC) [14], and immobilized artificial membrane chromatography (IAM) [15]. Also, many software applications such as ALOGP, KOWIN, and CLOGP that incorporate different calculation approaches and algorithms have been frequently applied. However, these commonly used programs for estimation are not suitable in the case of complex compounds simply because the appropriate input of the central metal atom is usually missing, therefore resulting in poor predictions [16].

While considerable effort has been devoted to studying the of organic compounds, to the best of our knowledge, only few papers are dealing with the measure of of various complexes mostly by the means of the shake-flask method, HPLC in various reversed-phase modalities, and MEKC [16], probably because of the fact that complex compounds exhibit a very intricate behavior in liquid-liquid partitioning systems including HPLC that requires a delicate approach [17].

However, despite its numerous advantages and long history of successful application for lipophilicity determination of small organic molecules, modern high-performance thin-layer chromatography is rarely used for determination of complex compounds, except for a few cases [18].

On the other hand, modern computational approach to the estimation of lipophilicity of metal complexes mostly relies on quantitative structure property relationship studies (QSPR) that establish quantitative models based on experimentally determined data and ab initio calculated molecular descriptors, such as the case with numerous platinum-based complexes [19, 20]. However, complete, reliable, and comparable data which correlate the with relevant properties and molecular descriptors are needed. Molecular descriptors often reflect complex and multiple physical interactions, and model should include as many descriptors as possible, often leading to difficult and ambiguous interpretation of the results. Another approach based on ab-initio calculations is to estimate the free solvation energy change for the solute transfer from n-octanol to pure aqueous phase. Modern theoretical methods in quantum chemistry, such as density functional theory (DFT), possess great predictive power and in conjunction with continuum solvation models are proven to be very reliable in the determination of the free energy of solvation [21]. It is noteworthy to mention that advantages of DFT are particularly important for transition metal compounds [22, 23], although they are not always considered as innocent systems. Continuum solvation models are efficient tools for studying solvent effects on molecular structure, spectra, and energetics [2427] and have been used with success for determination of partition coefficients [28], pKa values [29], redox potentials [30], and so forth. Improvement of solvent models is active area in research in computational chemistry, and models are continually being improved and new versions and models frequently appear [31]. From a great variety of solvent models, in this work the universal model proposed by Marenich et al, based on Density (SMD) [32] is used due to the proven accuracy for first-principle calculation of solvation energies [33].

In the present work the attention was focused on the application of RP-TLC as a simple, fast, and reliable tool to determine values of ten ruthenium(II)-arene complexes with potential anticancer and antiproliferative activity. Moreover, there is a serious lack of information in the present literature about lipophilicity assessment of complex compounds by the means of TLC. In order to assure accuracy of the proposed methodology we have employed ab-initio DFT computational approach to calculate values from the differences in Gibbs’ free solvation energies of the solute transfer from n-octanol to water and to compare estimated values with experimentally determined ones.

Being simple, fast and reliable RP-TLC provides retention data in the form of and corresponding values ((1) and (2)) that can be further used to derive several chromatographic descriptors for lipophilicity estimation: , , , and PC1. Consider the following: Parameters and represent migration distance of the solute and solvent front, respectively. Consider the following:

The first parameter represents the intercept in the linear relation between retention of a solute, , and volume fraction of organic-mobile-phase modifier, (3). As a value extrapolated to 0% v/v of organic solvent it accounts for partitioning of a solute between pure water and the nonpolar stationary phase [34].

Slope, , indicates the rate at which the solubility of the solute increases with changes in the mobile-phase composition and it is related to the specific hydrophobic surface area of a solute molecule [35]. Parameter , introduced by Bieganowska et al. [36], represents the ratio of and and is considered as a concentration of organic modifier in the mobile phase for which the distribution of the solute between the two phases is equal.

The last descriptor, the first principal component, PC1, is derived from principal component analysis (PCA), that is, principal component regression (PCR), multivariate chemometric methods often applied on chromatographic data. It has been demonstrated that the scores of the principal components (usually the PC1 scores are sufficient) are better correlated with since PCA combines all chromatographic data in one single feature [37].

2. Experimental

2.1. Materials

All reagents and solvents, including standard substances, were of analytical or HPLC purity grade. They are purchased from commercial suppliers: Aldrich (Milwaukee, WI, USA), Fluka (Buchs, Switzerland), and Merck (Darmstadt, Germany), and used as received. Ruthenium(II)-arene complex compounds (Figure 1) were synthesized and characterized as already described [7, 8].

In the present work a set of 28 standard solutes, mainly mono- and polysubstituted phenols, aromatic carboxylic acids, ketones, amines, esters, and few polyaromatic hydrocarbons, of known values (Table 1) were chosen for calibration. Experimentally determined values of were taken from the KOWWIN software (EPI Suite, v. 4.11 US EPA). The set of standard compounds was compiled considering that hydrogen bonds, for both proton accepting and proton donating, as well as dipolar interactions, among the overall interactions that solutes exhibit in a chromatographic environment, should be present in a way as to ensure the ability of the final model to describe the behavior of a diverse set of compounds. The optimal range of values was considered broad enough to provide a reliable regression performance (from 1 to 5 units). Dissociation constants for ionogenic compounds were collected from several sources and presented along with the calculated degree of ionization at in Table 1.

2.2. Chromatographic Procedure

For all chromatographic experiments solutions of standard substances, as well as of studied ruthenium complex compounds, were prepared by dissolving appropriate amount of substance in acetone in concentration of 0.1 mg/mL. Commercially available octadecyl modified silica aluminum sheets (Art. number 5559, Merck, Darmstadt, Germany) were cut into  cm plates. The plates were manually spotted with approximate volume of 1.0 μL of freshly prepared solutions, at 5 mm distance from each edge of the plate. Mixtures of acetonitrile and water were used as mobile phases, with organic component content increasing in the range of 30–60% v/v, with increment of 5%. Other organic modifiers such as methanol and acetone have been tested, but some of the ruthenium complexes were either moved through the system nonselectively or completely retained at the starting points. They are, therefore, discarded from the further study. The chromatograms were developed using horizontal developing chamber (CAMAG, Mutenz, Switzerland). The solvent migration distance was about 4.5 cm. The plates were visually inspected under the UV light (254 nm) and each zone was clearly marked and its distance was manually measured. All measurements were done in triplicate and average values were used in further calculations. All experiments were performed at ambient temperature (°C).

2.3. Calculation of Chromatographic Descriptors

The following chromatographic descriptors: , , and , were calculated based on the retention data (Table 2). and were obtained as the intercept and slope, respectively, according to (3). The hydrophobicity parameter was calculated as the ratio . However, this parameter did not show promising retention— models, and was excluded from the further study. All necessary calculations were performed using the Data analysis tool-pack (Microsoft Excel 2010). The obtained models are summarized in Table 3 followed by accompanying statistics.

2.4. Modeling and Validation of TLC Retention— Relationship

Calibration models based on individual values as well as extrapolated and descriptors were established using standard set of compounds (Table 1) and ordinary least squares procedure (Data analysis tool-pack, Microsoft Excel 2010). Principal component regression was employed instead of separate calculation of PC1 score values, using singular value decomposition (SVD) algorithm, followed by random subset cross-validation (five splits and one iteration) procedure. All calculations were made with PLS Tool Box v. 7.2., (Eigenvectors Research Inc.) for MATLAB R2011a (Mathworks Inc.). Quality of obtained models was assessed by the means of Pearson’s and accompanying statistics.

2.5. DFT Computational Details

The common logarithm of is calculated as where is the standard state solvation free energy of a given complex in octanol (oct) or in water (water) at  K. The standard-state solvation free energy is defined as the free energy of transfer from the gas phase to the condensed phase under standard state conditions. Because the gas-phase free energies are calculated with respect to a standard state of 1 atm, a correction factor of (1.894 kcal mol−1 at 298 K) needs to be added to convert it into the standard state of 1 mol dm−3.

All the DFT calculations have been carried out with the Gaussian 09, revision C.01 electronic structure program suite [38]. The local density approximation characterized by the Vosko-Wilk-Nusair (SVWN5) parametrization [39] has been used for the gas phase geometry optimizations, proven to be accurate for geometries of Werner type complexes [4042]. The ruthenium cation was described by LANL2DZ relativistic effective core potentials (ECP) that replace 28 core electrons with a nonlocal effective potential and associated basis set for remaining electrons [43]. 6-31G+(d,p) basis set [44] was used for all other atoms. For the determination of , M06L density functional [45, 46], with the ECP LANL2TZ basis set [47] for ruthenium, and 6-311+G(2d,p) for other atoms were used. The free energies of solvation were calculated using continuum solvation model based on density (SMD) [32]. With SMD, the 298 K solvation Gibbs energy is defined as the difference between the solvent and gas electronic energies [48], necessitating corresponding gas-phase calculation.

3. Results and Discussion

For all studied compounds a good linearity between retention constants and the volume fraction of acetonitrile was obtained resulting in well-established chromatographic descriptors ( and ) (Table 3), even in the case of strongly polar or ionizable compounds as well as highly hydrophobic solutes. Although was initially calculated, this parameter was later excluded from further consideration, because of statistically significant, but poor, correlation with values (, ).

All calibration models were obtained using entire set of standard compounds. Accompanying statistics and equations are summarized in Table 4. Depending on the particular model, some of standard compounds were identified as outliers, with borderline statistical significance. However, all data were kept because of lack of reasonable justification for outlier removal. The best model was achieved using parameter , while direct calibration of retention in 55% acetonitrile-water system (model no. 4) showed the worst correlation . There is obvious deterioration in statistical performance of calibration models based on values with an increase of acetonitrile content in mobile phase; however, with exception of the model number 4, no model can be statistically justified as being preferred over the others. Use of direct calibration is more convenient compared to calculation of , , and PCR since it does not require additional experimental and computational work. However, and can be assessed for both highly lipophilic and polar compounds that exhibit measurable retention in different range of organic modifier volume ratios, in contrast to direct calibration based on isocratic chromatographic conditions. Because of its physical meaning, parameter might be the most suitable.

The calculated gas-phase geometries of all complexes under investigation are in excellent agreement with available X-ray crystal structures [7, 8] (Table 5).

Chromatographically determined values were obtained using retention data (Table 2) and established calibration models (Table 3). Data are summarized in Table 6. All estimation methods (models 1–7) give mutually coherent data, with no statistical difference in between. Also there is a good overall agreement between DFT calculated and experimentally determined values .

Studied ruthenium(II)-arene complexes exhibit unusually high lipophilicity, in the range of 1–4 log units, compared with reported values of different series of platinum(II) and platinum(IV) complexes [19, 20] as well as with some Ru-(-arene) compounds of similar structure that show significantly low octanol-water distribution coefficients [49]. However, all mentioned compounds have been positively charged, which might explain their low lipophilicity in comparison with neutral complexes investigated in current study.

4. Conclusion

The lipophilicity of ten ruthenium(II)-arene complexes with potential anticancer and antiproliferative activity has been determined by means of RPTLC on RP-18 silica as stationary and binary acetonitrile-water solvent systems as mobile phase. Based on retention data corresponding chromatographic descriptors for lipophilicity assessment have been obtained. Most of the experimental findings described above have been confirmed by DFT free energy calculations of complexes in octanol and water as solvents, using solvation model based on density (SMD). As results are promising, it can be considered as reliable tool for prediction of lipophilicity and rational design of coordination compounds with desired properties.

Conflict of Interests

The authors declare that they have no conflict of interests.

Acknowledgment

This work has been supported by the Ministry of Education, Science and Technological Development of Serbia, Grant no. 172017.