Skip to main content
  • Research article
  • Open access
  • Published:

Titanocene / cyclodextrin supramolecular systems: a theoretical approach



Recently, various metallocenes were synthesized and analyzed by biological activity point of view (such as antiproliferative properties): ruthenocenes, cobaltoceniums, titanocenes, zirconocenes, vanadocenes, niobocenes, molibdocenes etc. Two main disadvantages of metallocenes are the poor hydrosolubility and the hydrolytic instability. These problems could be resolved in two ways: synthetically modifying the structure or finding new formulations with enhanced properties. The aqueous solubility of metallocenes with cytostatic activities could be enhanced by molecular encapsulation in cyclodextrins, as well as the hydrolytic instability of these compounds could be reduced.


This study presents a theoretical approach on the nanoencapsulation of a series of titanocenes with cytotoxic activity in α-, β-, and γ-cyclodextrin. The HyperChem 5.11 package was used for building and molecular modelling of titanocene and cyclodextrin structures, as well as for titanocene/cyclodextrin complex optimization. For titanocene/cyclodextrin complex optimization experiments, the titanocene and cyclodextrin structures in minimal energy conformations were set up at various distances and positions between molecules (molecular mechanics functionality, MM+). The best interaction between titanocene structures and cyclodextrins was obtained in the case of β- and γ-cyclodextrin, having the hydrophobic moieties oriented to the secondary face of cyclodextrin. The hydrophobicity of titanocenes (logP) correlate with the titanocene-cyclodextrin interaction parameters, especially with the titanocene-cyclodextrin interaction energy; the compatible geometry and the interaction energy denote that the titanocene/β- and γ-cyclodextrin complex can be achieved. Valuable quantitative structure-activity relationships (QSARs) were also obtained in the titanocene class by using the same logP as the main parameter for the in vitro cytotoxic activity against HeLa, K562, and Fem-x cell lines.


According to our theoretical study, the titanocene/cyclodextrin inclusion compounds can be obtained (high interaction energy; the encapsulation is energetically favourable). Further, the most hydrophobic compounds are better encapsulated in β- and γ-cyclodextrin molecules and are more stable (from energetically point of view) in comparison with α-cyclodextrin case. This study suggests that the titanocene / β- and γ-cyclodextrin complexes (or synthetically modified cyclodextrins with higher water solubility) could be experimentally synthesized and could have enhanced cytotoxic activity and even lower toxicity.


Cancer is a generic name comprises a great number of medical affections, having various locations and symptoms [15]. Even this disease is studied more than fifty years, the cause and action mechanisms are not completely elucidated [3, 6, 7]. Chemotherapy is widely used in order to cure this disease, by using various cytostatic or cytotoxic compounds: alkylating agents, antimetabolites, hormones, immunostimulating agents, antibiotics, alkaloids, all of them with higher toxicity [1, 8].

Organometallic compounds is an important class used in chemotherapy and the main groups studied are metallocenes (compounds which contains two cyclopentadienyl anions bound to a metal centre in the oxidation state II), ruthenium-, osmium-, iridium half-sandwich complexes, rhenium organometallics, metal N-heterocyclic carbene complexes, metal carbonyl complexes, or miscellaneous organometallic compounds [9]. The actual trend in cancer treatment is to replace some of the more toxic drugs such as cisplatin with less toxic compounds. Organometallic compounds are widely studied from cytotoxic point of view. Ferrocene was one of the first organometallic compounds from the first group evaluated for its antiproliferative properties [9, 10]. Ferrocene derivatives were obtained as antimalarial or cytostatic drugs and drug candidates [9]. Recently, similar metallocenes were synthesized and analyzed by biological activity point of view: ruthenocenes, cobaltoceniums, titanocenes, zirconocenes, vanadocenes, niobocenes, molibdocenes etc. [923]. The titanocene compounds are promising such drugs, but the hydrolytic instability and slightly water solubility conduct to a lower cytotoxic activity (approximately ten fold lower than cisplatin). Further, in the titanocene series cytotoxic activity against HeLa, K562, and Fem-x cell lines increases with the overall hydrophobicity of compounds [16, 17, 22]. On the other hand, increasing the hydrophobicity of titanocenes conducts to a more lower water solubility and reducing the transport capacity in aqueous layers (even the transport capacity through lipid layers are increased) [16, 17]. Despite of the resemblance of titanocene dichloride derivatives with cisplatin, seems that the mode of action as anti-cancerigene is different: binding to DNA and apoptosis of the cancer cell for the cisplatin and binding to DNA phosphate group, with additional interaction stabilizing the binding to DNA, for titanocenes. Two main problems of the titanocene dihalides are the poor hydrosolubility and the hydrolytic instability [22].

These problems could be resolved in two ways: synthetically modifying the titanocene structure (laborious, other physico-chemical and biological analyses needed) or finding new formulations with enhanced properties [11, 13, 15, 19, 2427]. Natural or chemically modified cyclodextrins (cyclic oligosaccharides with hydrophobic inner cavities and hydrosolubilizing outer groups) are widely used for protection, enhancing water solubility, and controlled release properties of bioactive compounds [2834]. The aqueous solubility of metallocenes (i.e. titanocenes) with cytostatic or cytotoxic activities could be enhanced by molecular encapsulation in cyclodextrins, as well as the hydrolytic instability of these compounds could be reduced (by reducing the access of water molecules to the metallocene halide reaction centre) [13, 19, 25, 27].

This study presents a theoretical approach on the molecular encapsulation of a series of titanocenes with cytostatic activity in α-, β-, and γ-cyclodextrin, in order to obtain supramolecular systems with enhanced stability and bioavailability. Further, a quantitative structure-biological activity relationships (QSAR) studies were performed in order to evaluate the main parameters which influencing the in vitro cytostatic activity.

Results and discussion

Quantitative structure-activity relationships (QSARs)

Molecular modelling and conformational analysis of titanocenes were performed by using the default parameters of MM+ functionality from the HyperChem 5.11 package [35]. This method is appropriate especially for organic molecules, but in our organometallic series the “inner” metal did not have a major influence on the overall geometry of compounds. Titanocenes with cytotoxic activity indicate a reduced number of stable conformations, especially for structures with a lower flexible bonds (codes 01TC, 02TC, 08TC to 10TC); structures having dimethyl-vinyl-silyl-ethyl or trivinyl-silyl-ethyl moieties at the Si, or pyridinium-methyl moieties at the cyclopentadienyl rings have a higher number of stable conformers; the most stable conformations are partially superimposed, especially at the titanocene skeleton. The most stable conformations have all flexible substituents (Y, R, and R’ from the general structure, Table 1 and Scheme 1) oriented close to a pseudoplan formed by Ti and the gravity centres of the two cyclopentadiene rings (Figure 1) [see Additional file 1.

Table 1 Titanocene structures (see Scheme 1) and in vitro cytotoxic activities (a)
Scheme 1
scheme 1

General structure of titanocene compounds.

Figure 1
figure 1

Superimposed minimum energy conformations of titanocene structures (OX, OY, and OZ view).

In order to evaluate the importance of titanocene structure/conformation to the overall cytotoxic activity, the most important steric, electronic, and hydrophobic descriptors were evaluated by using QSAR Properties functionality. Molecular van der Waals surface and volume, refractivity, and polarizability for titanocenes in minimum energy conformations were evaluated from sterical descriptor class (Table 2), but some intercorrelation exists, especially in the case of SvdWVvdW and RfPol pairs (intercorrelation coefficient of 0.87 and 1.0, respectively, Table 3). This intercorrelation is poor for SvdW/VvdW and Rf/Pol. Hydration energy has no correlation with descriptors from the other classes, but the transport parameter, logP, correlates with Rf and Pol parameters with higher intercorrelation coefficients (inverse correlation of −0.85 and −0.87, respectively, Table 3).

Table 2 Values of the structural descriptors (a) for minimum energy conformations of titanocenes
Table 3 Intercorrelational matrix for titanocene structural descriptors

Taking into account these intercorrelations, only logP, Rf, and Pol reveals statistically significant quantitative structure – cytotoxic activity relationships (QSARs) [36, 37], which could drive molecular changes in this titanocene class in order to obtain new compounds with higher biological activity. In the case of HeLa cell line, in vitro cytotoxic activity of titanocenes correlates with the logP parameter (Eq. 1), having a correlation coefficient of 0.80, which is higher than in the other cases; the correlation coefficient for the case of K562 is little bit lower (Eq. 2), but statistically significant, as well as in the case of Fem-x cell line (Eq. 3). A significant correlation was also obtained if the polarizability, Pol, was used as independent parameter (Eq. 4); for example, an inverse correlation exists between cytotoxic activity against K562 cells (expressed as pIC 50 ) and polarizability (r = 0.80), but this parameter is intercorrelated with the hydrophobic one. As a result, if the polarizability is lower, the hydrophobicity is increased and the transport, the interaction of titanocene with the hydrophobic cell membrane is favoured, and further the cytotoxic activity is increased. The predictive accuracy was evaluated by cross-validation method (leave-half-out method), all models having cross-validation correlation coefficients q2 cv > 0.75 [see Additional file 1. The predicted pA activities, calculated with equations Eq. 1 to Eq. 4, as well as the differences between experimental and predicted activities, are presented in Table 4 [see Additional file 1.

Table 4 QSAR results for cytotoxic activity of titanocenes against HeLa, K562, and Fem-x cell lines (experimental activitiesA and p A, predicted activitiesp A ( pred .) , and the differences between experimental and predicted activities, Δp A)

QSAR results (Table 4 and Eqs. 1–4) indicate very good predictions of cytotoxic activity against HeLa and Fem-x cell lines, especially in the case of compounds 01TC-11TC (|ΔpA| < 0.11, except compound 09TC for HeLa cytotoxic activity; this compound contains Ge as the second metal atom, in comparison with the other titanocene compounds, which contain only Si instead). Higher differences were observed in the case of pyridinium-methyl derivatives which have no Si as the second metal atom; the presence of Si atom in the structure stabilizes the overall titanocene structure (compounds 23TC, 24TC, and 26TC). Good results were obtained for prediction of the cytotoxic activity against K562 cell lines for models obtained with logP and Rf as structural parameters (Eqs. 2 and 4). The difference between experimental and predicted activity was lower than 0.2 in almost all titanocene cases (except compound 18TC, having trivinyl-silyl moiety and the highest Rf value) (Table 4). However, a correlation between predicted activities against K562 cell lines obtained with Eqs. 2 and 4 exists (r = 0.67); this result is sustained by the logPRf intercorrelation coefficient which are involved in Eqs. 2 and 4 (Table 3).

p A 1 i = 3.40 ± 0.17 + 0.69 ± 0.17 log P i

n = 11; r = 0.80; F = 16; q2 cv = 0.75

p A 2 i = 3.54 ± 0.23 + 0.69 ± 0.30 l o g P i

n = 8; r = 0.70; F = 5.5; q2 cv = 0.88

p A 3 i = 3.63 ± 0.12 + 0.39 ± 0.16 l o g P i

n = 7; r = 0.74; F = 5.9; q2 cv = 0.80

p A 2 i = 5.27 ± 0.37 0.014 ± 0.004 R f i

n = 11; r = 0.80; F = 10.5; q2 cv = 0.75

Geometry optimization of titanocene / cyclodextrin supramolecular systems

It is known that the cytotoxic activity of titanocene dichloride is different from cisplatin action (the last being more toxic, but with higher cytotoxic activity): titanocenes conduct to adduct with DNA and prevent the replication and/or transcription, resulting in cell death [22]. Thus, the transport of titanocene to the DNA is very important, but the lower water solubility and instability (hydrolysis of chloride ligands) reduces the access to the target - biomacromolecules; on the other hand, if the hydrophobicity is increased (the QSAR study indicates a higher in vitro biological activity for more hydrophobic titanocenes), the titanocene transfer in aqueous layer is decreased (especially in the case of the in vivo experiments). Thus, the increase of the hydrophobicity in order to enhance the in vitro cytotoxic activity conducts to a less water solubility and a harder transport of the titanocenes to the target. The water solubility of hydrophobic molecules could be realized by molecular encapsulation in matrices such as cyclodextrins. Among the increasing of water solubility, the advantage of this procedure is to protect of easier hydrolysable titanocenes against degradation and controlled release to the target (i.e. DNA); further, the higher hydrophobicity of modified titanocenes conduct to a better interaction with the hydrophobic inner cavity of cyclodextrins.

Fully geometry optimization of titanocene / cyclodextrin supramolecular systems or docking of organometallic compounds in cyclodextrins [38, 39] could provide information on the stability of complexes and suggest chemical modifications for new titanocenes with higher cytotoxic activity. In our theoretical experiments of fully geometry optimization of titanocene / cyclodextrin supramolecular systems in vacuum (the default HyperChem molecular mechanics MM+ force field was used), only the interaction of the hydrophobic moiety of titanocene with the inner cavity of cyclodextrins from the secondary face was efficient (higher stability of the complex). The main interactions which stabilize the titanocene / cyclodextrin complex in vacuum were bond stretching energy, the angle bending energy, the torsional energy, and the energy arising from van der Waals interactions of non-bonded pairs of atoms. Due to the fact that all titanocene compounds from the studied series are chemically similar molecules, we have an internal consistency from the force field point of view [15, 4049].

In the case of α-cyclodextrin, the interaction of titanocenes with the secondary face is poor, even with structures having thin hydrophobic chains (e.g. 03TC, 18TC, 23TC, 24TC, and 26TC; Figure 2) [see Additional file 1]; the worst interaction was observed for the case of titanocene oriented to the primary face of α-cyclodextrin or with the titanocene moiety in the front of both cyclodextrin rings. Theoretical interaction energy between titanocenes and α-cyclodextrin was in the range of 11.9 kcal/mole to 19.7 kcal/mole (Table 5). The best interactions of titanocenes with α-cyclodextrin were obtained in the case of structures having compounds monosubstituted at the Si (02TC and 10TC, E int. 19.7 kcal/mole and 17.0 kcal/mole, respectively) or having flexible chain connected to this atom (03TC and 18TC, E int. 18.0 kcal/mole and 17.0 kcal/mole, respectively). Lower interaction energies between titanocenes and α-cyclodextrin were obtained for compounds di-substituted at the Si atom with short chain moieties (methyl or vinyl) and for those substituted at the cyclopentadienyl rings with pyridinium-methyl moieties (Table 5). These low interaction energies are especially due to the geometric compatibilities, which are close to the limits and the interactions with the hydrophobic inner cavity of α-cyclodextrin are reduced. This is the reason that the hydrophobicity (expressed as logP) is not important (r < 0.5) in relation with the titanocene / α-cyclodextrin interaction energy, while the hydration energy (E hydr. ) weak correlates with this interaction energy (Eq. 5).

E i n t . , a C D , i = 15.54 ± 0.57 + 0.356 ± 0.137 E h y d r . , a C D , i
Figure 2
figure 2

Titanocene (code 03TC) / α-cyclodextrin supramolecular system (theoretically modeled by MM+).

Table 5 Energies (resulted from the MM+ molecular modeling and titanocene/cyclodextrin optimization experiments) for cyclodextrins (E CD , α-, β-, and γ-cyclodextrin, codes aCD, bCD, and gCD), titanocenes (E TC . , codes xTC, where x = 0103, 0811, 18, 23, 24, and 26), the sum of titanocene and cyclodextrin energies, with no interaction (E TC .+ CD ), the energies of the TC-CD complex (E TC .- CD complex ), and the TC-CD interaction energies (E int . ), determined as the difference between the TC+CD energy, with no interaction, and the energy of the TC-CD complex

n = 11; r = 0.656; F = 6.8

The superior homologues, β- and γ-cyclodextrins, are more proper for molecular encapsulation of the studied titanocenes (Figures 3 and 4). In all β- and γ-cyclodextrin complexes the mono- or disubstituted silyl moiety, the bulky moiety (dimethyl-vinyl-silyl or trivinyl-silyl substituents) from the edge of the corresponding substituent, or the pyridinium-methyl substituent from the cyclopentadienyl ring are almost completely encapsulated in the cyclodextrin cavities [see Additional file 1]. In these cases the hydrophobicity of titanocenes (quantified by the logarithm of the 1-octanol/water partition coefficient) was important for the molecular encapsulation process and the interaction energy correlates with logP in a statistically significant way (Eq. 6 and Eq. 7). Thus, better interactions were obtained with the titanocenes having pyridinium-methyl substituents on the cyclopentadienyl rings (codes 24TC and 26TC), the interaction energy being 23.7-24.3 kcal/mole and 27.3-28.1 kcal/mole for β- and γ-cyclodextrin, respectively (Table 5). Higher interaction energies were obtained also for titanocenes having monosubstituted silyl moieties, 02TC and 10TC, especially in the case of γ-cyclodextrin complexes (28.4 kcal/mole and 28.6 kcal/mole, respectively; Table 5). Furthermore, a statistically significant correlation between titanocene / β- and γ-cyclodextrin interaction energies was observed (r = 0.82) [see Additional file 1].

E int . , b C D , i = 16.24 ± 1.64 + 4.459 ± 1.687 l o g P b C D , i
Figure 3
figure 3

Titanocene (code 23TC) / β-cyclodextrin supramolecular system (theoretically modeled by MM+).

Figure 4
figure 4

Titanocene (code 01TC) / γ-cyclodextrin supramolecular system (theoretically modeled by MM+).

n = 11; r = 0.661; F = 7.0

E i n t . , g C D , i = 21.46 ± 1.99 + 4.347 ± 2.047 l o g P g C D , i

n = 11; r = 0.578; F = 4.1

The stability of the titanocene / cyclodextrin supramolecular system is observed from Figure 5 for β- and γ-cyclodextrin cases: the better stability is obtained in the case of γ-cyclodextrin, followed by β-cyclodextrin case, and finally for α-cyclodextrin complex after more than double number of iterations. This variation of the energy interaction reveals the possibility of forming the titanocene/cyclodextrin complexes and the stability of these supramolecular systems [see Additional files 1 and 2].

Figure 5
figure 5

Interaction energy vs . Number of cycles from MM+ titanocene/cyclodextrin geometry optimization experiments for 23TC.


In our study the importance of the hydrophobic parameter (logP – the logarithm of the octanol/water partition coefficient) in both of in vitro cytotoxic activity against HeLa, K562, and Fem-x cell lines, as well as in the cyclodextrin nanoencapsulation of the titanocene compounds were demonstrated.

Our theoretical studies demonstrates that the molecular encapsulation of titanocenes in natural cyclodextrins could resolves some of cytotoxic titanocene disadvantages: a more hydrophobic titanocene, which have higher cytotoxic activity (the in vitro cytotoxic activity against HeLa, K562, and Fem-x cell lines is increased with the logP of the titanocene compound; r > 0.7), is better encapsulated in cyclodextrins and could be transported through the aqueous layers (cyclodextrin inclusion compounds are water soluble) and protected against hydrolysis. Our theoretical study on the titanocene / cyclodextrin complexes indicate that the hydrophobic biologically active compounds (with higher cytotoxic activity) are better encapsulated in β- and γ-cyclodextrin; the highest titanocene / β- and γ-cyclodextrin complex interaction energies of ~24 kcal/mole and ~28 kcal/mole was obtained in vacuum for bis(pyridinium-methyl)-silyl derivatives, respectively. This study indicate that the titanocene compound can be controlled released to the target from the cyclodextrin complex; this complex allow to transfer the more hydrophobic titanocene through lipid layers, to increase the concentration of bioactive compound to the DNA phosphoester groups and further to form the titanocene-DNA adducts. Due to this process the inhibition of DNA transcription and/or replication appears (cytotoxicity).

According to our theoretical study, these titanocene/cyclodextrin inclusion compounds can be obtained (the encapsulation process is energetically favourable for β- and γ-cyclodextrin complexes). Further, the most hydrophobic compounds are better encapsulated in β- and γ-cyclodextrin molecules and are more stable (from energetically point of view) in comparison with α-cyclodextrin case. This study suggests that the titanocene/β- and γ-cyclodextrin complexes (or synthetically modified cyclodextrins with higher water solubility) could be experimentally synthesized and could have enhanced cytotoxic activity and even lower toxicity.


Titanocene structure selection and cytotoxic activity

Titanocenes with potential anticarcinogenic properties were recently synthesized by Gómez-Ruiz et al.[16, 17] and Potter et al.[22] and have structural variability at the cyclopentadienyl moieties. All these new titanocene structures (eleven organometallic compounds) with cytotoxic activity against cervical carcinoma cell line HeLa, human myelogenous leukemia cell line K562, and human malignant carcinoma cell line Fem-x were considered in this theoretical study. Cytotoxic activity was expressed as the logarithm of the inverse inhibitory concentration 50%, pIC 50 = log(1/IC 50 ); pA 1 , pA 2 , pA 3 were used for cytotoxic activity against HeLa, K562, and Fem-x, respectively (Table 1).

Molecular modelling

Molecular modelling of titanocene molecules as well as α-, β-, and γ-cyclodextrins was performed by using the molecular mechanics MM+ functionality from the HyperChem 5.11. The MM+ molecular mechanics force field with a RMS of 0.005 kcal/mole, a number of maximum cycles to limit the search directions of fifteen times the number of atoms, and a Polak-Ribiere algorithm (a gradient method using one-dimensional searches) were used in the molecular modelling process. Bond dipole was used to calculate all nonbonded electrostatic interactions. In the MM+ calculations potential energy depends on bond lengths, bond angles, torsion angles, and nonbonded interactions (van der Waals forces, electrostatic interactions, and hydrogen bonds) [50, 51].

Conformational analysis

In order to find the most stable conformation even for titanocenes or cyclodextrins, a conformational analysis by using Conformational Search functionality (HyperChem 5.11) was performed. In titanocene structures only some side chains have flexible bonds. On the other hand, the flexible bonds in cyclodextrins were only those corresponding to the hydroxymethyl from C5 position of glucopyranose unit; the flexible rings were all glucopyranose rings and the corresponding macrocyclic ring. The following conditions were set up for conformational search: variation of the flexible torsion angles ±60º ÷ ±180º, energy criterion for acceptance of the conformation 4 kcal/mole above minimum, all conformations with atomic distances lower than 0.5 Å, and differences between torsion angles lower than 15º were not considered as well as conformations with energy differences lower than 0.05 kcal/mole (duplicates); the maximum number of optimization and iterative calculations was 1000 and maximum 20 conformations were retained. The hydrogen atoms were neglected.

Geometry optimization of titanocene / cyclodextrin supramolecular systems

The geometry optimization of titanocene (the most stable conformations) / cyclodextrin (α-, β-, and γ-cyclodextrin) complexes was realized by using the molecular mechanics interactions of the host-guest molecules in vacuum. The titanocene and cyclodextrin structures in minimal energy conformations were set up at distances of ~8Å between the gravity centres of the host-guest molecules, and the titanocene structure was oriented with the hydrophobic side chain in front of the primary (A) or secondary (B) face of cyclodextrin (the principal axis corresponding to the biocompound or side chain moiety was perpendicular to the A or B plan of cyclodextrin). The complex was modeled in absence of water molecules by using the same MM+ functionality and the interaction was stopped when the RMS gradient was lower than 0.005 kcal/mole. The titanocene-cyclodextrin interaction energy was evaluated as the difference between the overall energies of these two molecules and the energy of the complex.

Structural parameters, correlations, and QSARs

The main molecular descriptors of titanocenes were evaluated by using QSAR Properties functionality from the HyperChem 5.11 package. The following descriptors were calculated and were used as structural parameters for obtaining correlations with titanocene-cyclodextrin interaction energy or quantitative structure-activity relationships (QSARs): van der Waals molecular surface (SvdW, Å2; van der Waals surface area was carried out by an approximate developed by Still and co-workers [52, 53]), van der Waals molecular volume (VvdW, Å3; the grid method described by Bodor et al. [54] was used for van der Waals volume calculation. The QSAR Properties functionality uses the atomic radii of Gavezotti [55] for this method), hydration energy (E hydr , kcal/mole; the method of calculation is developed by Ooi et al. [56]. The calculation is based on exposed surface area, and employs the surface area as computed by the approximate method, weighted by atom type), logarithm of the octanol/water partition coefficient (logP; it was calculated by means of atomic contributions [57, 58]. Metal atoms were excepted from calculation; this approximation did not significantly affect the final results due to the “hidden” positions of these metal atoms), refractivity (Rf, Å3; the refractivity is estimated by the same method as logP, presented by Ghose and Crippen [59]), and polarizability (Pol, Å3; this parameter was estimated according to an additive method of Miller [60], where the different increments are associated with different atom types).

Monolinear correlations (Eqs. 8a and 8b) in theoretical titanocene-cyclodextrin interaction experiments and QSARs were evaluated by using interaction energy between titanocene and cyclodextrin molecules (E int. ) or the above mentioned pA cytotoxic activity and structural parameters (P).

E i n t . , i = a 0 + b P i
p A i = a 0 + b P i


  1. Avendaño C, Menéndez JC: Medicinal chemistry of anticancer drugs. 2008, Amsterdam: Elsevier, 1-3.

    Book  Google Scholar 

  2. Quintans-Júnior L, Fagundes da Rocha R, Freitas Caregnato F, Fonseca Moreira JC, Amaral da Silva F, Antunes de Souza Araújo A, Almeida dos Santos JP, Santos Melo M, Pergentino de Sousa D, Rigoldi Bonjardim L, Pens Gelain D: Antinociceptive Action and Redox Properties of Citronellal, an Essential Oil Present in Lemongrass. J Med Food. 2011, 14: 630-639. 10.1089/jmf.2010.0125.

    Article  Google Scholar 

  3. Clapp RW, Jacobs MM, Loechler EL: Environmental and occupational causes of cancer. 2007, Lowell: The Lowell Center for Sustainable Production, University of Massachusetts, 24-28.

    Google Scholar 

  4. Vogelstein B, Kinzler KW: The Genetic Basis of Human Cancer. 2002, New York: McGraw-Hill, Medical Publication Division

    Google Scholar 

  5. Tong R, Cheng J: Anticancer polymeric nanomedicines. Polym Rev. 2007, 47: 345-381. 10.1080/15583720701455079.

    Article  CAS  Google Scholar 

  6. Andersen M, Kiel P, Larsen H, Maxild J: Mutagenic action of aromatic epoxy resins. Nature. 1978, 276: 391-392. 10.1038/276391a0.

    Article  CAS  Google Scholar 

  7. Chaudhary A, Pandeya SN, Kumar P, Sharma PP, Gupta S, Soni N, Verma KK, Bhardwaj G: Combretastatin A-4 analogs as anticancer agents. Mini Rev Med Chem. 2007, 7: 1186-1205. 10.2174/138955707782795647.

    Article  CAS  Google Scholar 

  8. Cragg GM, Kingston DGI, Newman DJ: Anticancer agents from natural products. 2005, Boca Raton: CRC Press

    Google Scholar 

  9. Gasser G, Ott I, Metzler-Nolte N: Organometallic anticancer compounds. J Med Chem. 2011, 54: 3-25. 10.1021/jm100020w.

    Article  CAS  Google Scholar 

  10. Vessières A, Plamont M-A, Cabestaing C, Claffey J, Dieckmann S, Hogan M, Müller-Bunz H, Strohfeldt K, Tacke M: Proliferative and anti-proliferative effects of titanium- and iron-based metallocene anti-cancer drugs. J Organomet Chem. 2009, 694: 874-879. 10.1016/j.jorganchem.2008.11.071.

    Article  Google Scholar 

  11. Ashton PR, Balzani V, Clemente-León M, Colonna B, Credi A, Jayaraman N, Raymo FM, Stoddart JF, Venturi M: Ferrocene-containing carbohydrate dendrimers. Chem Eur J. 2002, 8: 673-684. 10.1002/1521-3765(20020201)8:3<673::AID-CHEM673>3.0.CO;2-D.

    Article  CAS  Google Scholar 

  12. Campbell KS, Dillon CT, Smith SV, Harding MM: Radiotracer studies of the antitumor metallocene molybdocene dichloride with biomolecules. Polyhedron. 2007, 26: 456-459. 10.1016/j.poly.2006.07.004.

    Article  CAS  Google Scholar 

  13. Casas-Solvas JM, Ortiz-Salmeron E, Fernandez I, Garcia-Fuentes L, Santoyo-Gonzalez F, Vargas-Berenguel A: Ferrocene-β-cyclodextrin conjugates: synthesis, supramolecular behavior, and use as electrochemical sensors. Chem Eur J. 2009, 15: 8146-8162. 10.1002/chem.200900593.

    Article  CAS  Google Scholar 

  14. Chohan ZH, Sumrra SH, Youssoufi MH, Hadda TB: Metal based biologically active compounds: Design, synthesis, and antibacterial/antifungal/cytotoxic properties of triazole-derived Schiff bases and their oxovanadium(IV) complexes. Eur J Med Chem. 2010, 45: 2739-2747. 10.1016/j.ejmech.2010.02.053.

    Article  CAS  Google Scholar 

  15. Fey N: Organometallic molecular modelling - the computational chemistry of metallocenes: a review. J Chem Technol Biotechnol. 1999, 74: 852-862. 10.1002/(SICI)1097-4660(199909)74:9<852::AID-JCTB131>3.0.CO;2-T.

    Article  CAS  Google Scholar 

  16. Gómez-Ruiz S, Kaluđerović GN, Polo-Cerón D, Prashar S, Fajardo M, Zižak Z, Juranić ZD, Sabo TJ: Study of the cytotoxic activity of alkenyl-substituted ansa-titanocene complexes. Inorg Chem Commun. 2007, 10: 748-752. 10.1016/j.inoche.2007.03.016.

    Article  Google Scholar 

  17. Gómez-Ruiz S, Kaluđerović GN, Prashar S, Polo-Cerón D, Fajardo M, Zižak Z, Sabo TJ, Juranić ZD: Cytotoxic studies of substituted titanocene and ansa-titanocene anticancer drugs. J Inorg Biochem. 2008, 102: 1558-1570. 10.1016/j.jinorgbio.2008.02.001.

    Article  Google Scholar 

  18. Immel TA, Martin JT, Dürr CJ, Groth U, Huhn T: Dimethyl titanocene Y: a valuable precursor for libraries of cytotoxic titanocene derivatives. J Inorg Biochem. 2010, 104: 863-867. 10.1016/j.jinorgbio.2010.04.003.

    Article  CAS  Google Scholar 

  19. Liu Y, Zhong R-Q, Zhang H-Y, Song H-B: A unique tetramer of 4: 5 β-cyclodextrin-ferrocene in the solid state. Chem Commun. 2005, 17: 2211-2213. 10.1039/b418220k.

    Article  Google Scholar 

  20. Lu Z, Lu C, Ren X, Meng Q: New metallocene-bridged cyclodextrin dimer: A stable derivative of the antitumor drug titanocene dichloride and its potent cytotoxity against human breast cancer (MCF-7) cells. J Organomet Chem. 2006, 691: 5895-5899. 10.1016/j.jorganchem.2006.09.052.

    Article  CAS  Google Scholar 

  21. Napoli M, Saturnino C, Sirignano E, Popolo A, Pinto A, Longo P: Synthesis, characterization and cytotoxicity studies of methoxy alkyl substituted metallocenes. Eur J Med Chem. 2011, 46: 122-128. 10.1016/j.ejmech.2010.10.021.

    Article  CAS  Google Scholar 

  22. Potter GD, Baird MC, Cole SPC: A new series of titanocene dichloride derivatives bearing cyclic alkylammonium groups: Assessment of their cytotoxic properties. J Organomet Chem. 2007, 692: 3508-3518. 10.1016/j.jorganchem.2007.04.024.

    Article  CAS  Google Scholar 

  23. Wallis D, Claffey J, Gleeson B, Hogan M, Müller-Bunz H, Tacke M: Novel zirconocene anticancer drugs?. J Organomet Chem. 2009, 694: 828-833. 10.1016/j.jorganchem.2008.08.020.

    Article  CAS  Google Scholar 

  24. Braga SS, Almeida Paz FA, Pillinger M, Seixas JD, Romão CC, Gonçalves IS: Structural studies of β-cyclodextrin and permethylated β-cyclodextrin inclusion compounds of cyclopentadienyl metal carbonyl complexes. Eur J Inorg Chem. 2006, 8: 1662-1669. 10.1002/ejic.200501006.

    Article  Google Scholar 

  25. Morales A, Weber RT, Melendez E: Spectroscopic and thermal characterization of the host-guest interactions between α-, β- and γ -cyclodextrins and vanadocene dichloride. Appl Organomet Chem. 2008, 22: 440-450. 10.1002/aoc.1420.

    Article  CAS  Google Scholar 

  26. Singh R, Bharti N, Madan J, Hiremath SN: Characterization of cyclodextrin inclusion complexes - a review. J Pharm Sci Technol. 2010, 2: 171-183.

    CAS  Google Scholar 

  27. Sokolov VI: Cyclodextrin-metallocene inclusion complexes. Supramolecular stereochemistry. Edited by: Siegel JS. 1995, Dordrecht: Kluwer Academic Publishers, 239-245. 10.1007/978-94-011-0353-4_30.

    Chapter  Google Scholar 

  28. Brewster ME, Loftsson T: Cyclodextrins as pharmaceutical solubilizers. Adv Drug Deliv Rev. 2007, 59: 645-666. 10.1016/j.addr.2007.05.012.

    Article  CAS  Google Scholar 

  29. Hădărugă DI, Hădărugă NG, Bandur GN, Isengard H-D: Water content of flavonoid/cyclodextrin nanoparticles: relationship with the structural descriptors of biologically active compounds. Food Chem. 2012, 132: 1651-1659. 10.1016/j.foodchem.2011.06.004.

    Article  Google Scholar 

  30. Hădărugă DI, Hădărugă NG, Bandur GN, Riviş A, Costescu C, Ordodi V, Ardelean A: Berberis vulgaris extract/β-cyclodextrin nanoparticles: synthesis and characterization. Rev Chim. 2010, 61: 669-675. 10.1002/recl.19420610908.

    Article  Google Scholar 

  31. Hădărugă DI, Hădărugă NG, Butnaru G, Tatu C, Gruia A: Bioactive microparticles (10): thermal and oxidative stability of nicotine and its complex with β-cyclodextrin. J Incl Phenom Macrocycl Chem. 2010, 68: 155-164. 10.1007/s10847-010-9761-0.

    Article  Google Scholar 

  32. Hădărugă NG, Hădărugă DI, Păunescu V, Tatu C, Ordodi VL, Bandur GN, Lupea AX: Bioactive nanoparticles (6). thermal stability of linoleic acid / α- and β-cyclodextrin complexes. Food Chem. 2006, 99: 500-508. 10.1016/j.foodchem.2005.08.012.

    Article  Google Scholar 

  33. Szejtli J, Szente L: Elimination of bitter, disgusting tastes of drugs and foods by cyclodextrins. Eur J Pharm Biopharm. 2005, 61: 115-125. 10.1016/j.ejpb.2005.05.006.

    Article  CAS  Google Scholar 

  34. Serafini MR, Menezes PP, Costa LP, Lima CM, Quintans LJ, Cardoso JC, Matos JR, Soares-Sobrinho JL, Grangeiro S, Nunes PS, et al: Interaction of p-cymene with beta-cyclodextrin. J Therm Anal Calorim. 2012, 109: 951-955. 10.1007/s10973-011-1736-x.

    Article  CAS  Google Scholar 

  35. HyperChemTM release 5.11 Professional for Windows. 1999, FL, USA: Hypercube, Inc, Gainsville,,

  36. Hansch C, Maloney PP, Fujita T, Muir RM: Correlation of biological activity of phenoxyacetic acids with Hammett substituent constants and partition coefficients. Nature. 1962, 194: 178-180. 10.1038/194178b0.

    Article  CAS  Google Scholar 

  37. Hansch C, Fujita T: p-σ-π analysis. A method for the correlation of biological activity and chemical structure. J Am Chem Soc. 1964, 86: 1616-1626. 10.1021/ja01062a035.

    Article  CAS  Google Scholar 

  38. Menger FM, Sherrod MJ: Docking calculations on ferrocene complexation with cyclodextrins. J Am Chem Soc. 1988, 110: 8606-8611. 10.1021/ja00234a005.

    Article  CAS  Google Scholar 

  39. Sherrod MJ: Exploration of cyclomalto-oligosaccharide (cyclodextrin) chemistry with molecular mechanics: Docking calculations on the complexation of ferrocenes with cyclodextrins. Carbohydr Res. 1989, 192: 17-32. 10.1016/0008-6215(89)85161-4.

    Article  CAS  Google Scholar 

  40. Duchamp DJ: Molecular mechanics and crystal structure analysis in drug design. Computer-Assisted Drug Design. Volume 112. Edited by: Olson EC, Christoffersen RE. 1979, Washington, DC: American Chemical Society, 79-102. 10.1021/bk-1979-0112.ch003.

    Chapter  Google Scholar 

  41. Ulrich B, Allinger NL: Molecular mechanics. 1982, Washington: American Chemical Society

    Google Scholar 

  42. Leach AR: Molecular modelling. Principles and Applications. 2001, Harlow: Pearson Education Limited

    Google Scholar 

  43. Hinchliffe A: Modelling molecular structures. 2000, Chichester: John Wiley & Sons, Ltd.

    Google Scholar 

  44. Bowen JP, Allinger NL: Molecular Mechanics: The Art and Science of Parameterization. Reviews in Computational Chemistry. Volume 2. Edited by: Lipkowitz KB, Boyd DB. 2007, Hoboken, NJ: John Wiley & Sons, Inc, 10.1002/9780470125793.ch3.

    Google Scholar 

  45. Allinger NL, Zhou X, Bergsma J: Molecular mechanics parameters. J Mol Struct (THEOCHEM). 1994, 312: 69-83. 10.1016/S0166-1280(09)80008-0.

    Article  Google Scholar 

  46. Jianu C: Synthesis of nonionic-anionic colloidal systems based on alkaline and ammonium β-nonylphenol polyethyleneoxy (n = 3-20) propionates / dodecylbenzenesulfonates with prospects for food hygiene. Chem Cent J. 2012, 95: 10.1186/1752-153X-6-95.

    Google Scholar 

  47. Doman TN, Landis CR, Bosnich B: Molecular mechanics force fields for linear metallocenes. J Am Chem Soc. 1992, 114: 7264-7272. 10.1021/ja00044a042.

    Article  CAS  Google Scholar 

  48. Doman TN, Lollis TK, Bosnich B: Molecular mechanics force fields for bent metallocenes of the type [M(Cp)2Cl2]. J Am Chem Soc. 1995, 117: 1352-1368. 10.1021/ja00109a020.

    Article  CAS  Google Scholar 

  49. Timofeeva TV, Lii J-H, Allinger NL: Molecular mechanics explanation of the metallocene bent sandwich structure. J Am Chem Soc. 1995, 117: 7452-7459. 10.1021/ja00133a018.

    Article  CAS  Google Scholar 

  50. Comba P, Hambley TW: Molecular Modeling of Inorganic Compounds. 2001, Weinheim: Wiley-VCH

    Google Scholar 

  51. Yao S, Shoji T, Iwamoto Y, Kamei E: Consideration of an activity of the metallocene catalyst by using molecular mechanics, molecular dynamics and QSAR. Comput Theor Polym Sci. 1999, 9: 41-46. 10.1016/S1089-3156(98)00051-8.

    Article  CAS  Google Scholar 

  52. Hasel W, Hendrickson TF, Still WC: A rapid approximation to the solvent accessible surface areas of atoms. Tetrahedron Comput Meth. 1988, 1: 103-116. 10.1016/0898-5529(88)90015-2.

    Article  CAS  Google Scholar 

  53. Still WC, Tempczyk A, Hawley RC, Hendrickson T: Semianalytical treatment of solvation for molecular mechanics and dynamics. J Am Chem Soc. 1990, 112: 6127-6129. 10.1021/ja00172a038.

    Article  CAS  Google Scholar 

  54. Bodor N, Gabanyi Z, Wong CK: A new method for the estimation of partition coefficient. J Am Chem Soc. 1989, 111: 3783-10.1021/ja00193a003.

    Article  CAS  Google Scholar 

  55. Gavezzotti A: The calculation of molecular volumes and the use of volume analysis in the investigation of structured media and of solid-state organic reactivity. J Am Chem Soc. 1983, 105: 5220-5225. 10.1021/ja00354a007.

    Article  CAS  Google Scholar 

  56. Ooi T, Oobatake M, Némethy G, Scheraga HA: Accessible surface areas as a measure of the thermodynamic parameters of hydration of peptides. Proc Natl Acad Sci. 1987, 84: 3086-3090. 10.1073/pnas.84.10.3086.

    Article  CAS  Google Scholar 

  57. Ghose AK, Pritchett A, Crippen GM: Atomic physicochemical parameters for three dimensional structure directed quantitative structure-activity relationships III: modeling hydrophobic interaction. J Comput Chem. 1988, 9: 80-90. 10.1002/jcc.540090111.

    Article  CAS  Google Scholar 

  58. Viswanadhan VN, Ghose AK, Revankar GR, Robins RK: Atomic physicochemical parameters for three dimensional structure directed quantitative structure-activity relationships. 4. Additional parameters for hydrophobic and dispersive interactions and their application for an automated superposition of certain naturally occurring nucleoside antibiotics. J Chem Inf Comp Sci. 1989, 29: 163-172. 10.1021/ci00063a006.

    Article  CAS  Google Scholar 

  59. Ghose AK, Crippen GM: Atomic physicochemical parameters for three-dimensional-structure-directed quantitative structure-activity relationships. 2. Modeling dispersive and hydrophobic interactions. J Chem Inf Comp Sci. 1987, 27: 21-35. 10.1021/ci00053a005.

    Article  CAS  Google Scholar 

  60. Miller KJ: Additivity methods in molecular polarizability. J Am Chem Soc. 1990, 112: 8533-8542. 10.1021/ja00179a044.

    Article  CAS  Google Scholar 

Download references


This work was supported by Ministry of Education, Research, Youth, and Sports from Romania, PN2_ID_PCCE_140/2008. Authors are grateful to Professor Mircea Mracec (“Coriolan Drăgulescu” Institute of Chemistry, Timişoara, Romania) for permission to use the HyperChem 5.11 molecular modeling package.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Daniel I Hădărugă.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

NGH and DIH carried out all theoretical experiments and prepared the final manuscript. AR and ZG help to discuss the theoretical results. All authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: Molecular modelling and complex optimization of titanocene / cyclodextrin systems. Description: In this additional file the molecular modelling of titanocenes (A1), the quantitative structure-activity relationships (QSARs), which contain correlations between experimental and predicted activities, as well as the cross-validation data (leave-half-out method) for QSARs (A2), and complex optimization of titanocene/cyclodextrin complexes (starting positions and the most stable titanocene / cyclodextrin supramolecular systems; variation of the titanocene / cyclodextrin interaction energy in the complexation process; titanocene / cyclodextrin interaction energies correlations) (A3) are presented. These data supports all considerations presented in the manuscript. (DOC 4 MB)


Additional file 2: Titanocene / cyclodextrin interaction energies. Description: In this additional file the interaction energy versus the number of cycles data from MM+ optimization experiments for titanocene / cyclodextrin complexation process is presented; this supports all considerations presented in the manuscript. (XLS 932 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Riviş, A., Hădărugă, N.G., Gârban, Z. et al. Titanocene / cyclodextrin supramolecular systems: a theoretical approach. Chemistry Central Journal 6, 129 (2012).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: