- Research article
- Open Access
On the ease of predicting the thermodynamic properties of beta-cyclodextrin inclusion complexes
- Andreas Steffen^{1} and
- Joannis Apostolakis^{2}Email author
- Received: 18 September 2007
- Accepted: 15 November 2007
- Published: 15 November 2007
Abstract
Background
In this study we investigated the predictability of three thermodynamic quantities related to complex formation. As a model system we chose the host-guest complexes of β-cyclodextrin (β-CD) with different guest molecules. A training dataset comprised of 176 β-CD guest molecules with experimentally determined thermodynamic quantities was taken from the literature. We compared the performance of three different statistical regression methods – principal component regression (PCR), partial least squares regression (PLSR), and support vector machine regression combined with forward feature selection (SVMR/FSS) – with respect to their ability to generate predictive quantitative structure property relationship (QSPR) models for ΔG°, ΔH° and ΔS° on the basis of computed molecular descriptors.
Results
We found that SVMR/FFS marginally outperforms PLSR and PCR in the prediction of ΔG°, with PLSR performing slightly better than PCR. PLSR and PCR proved to be more stable in a nested cross-validation protocol. Whereas ΔG° can be predicted in good agreement with experimental values, none of the methods led to comparably good predictive models for ΔH°. In using the methods outlined in this study, we found that ΔS° appears almost unpredictable. In order to understand the differences in the ease of predicting the quantities, we performed a detailed analysis. As a result we can show that free energies are less sensitive (than enthalpy or entropy) to the small structural variations of guest molecules. This property, as well as the lower sensitivity of ΔG° to experimental conditions, are possible explanations for its greater predictability.
Conclusion
This study shows that the ease of predicting ΔG° cannot be explained by the predictability of either ΔH° or ΔS°. Our analysis suggests that the poor predictability of TΔS° and, to a lesser extent, ΔH° has to do with a stronger dependence of these quantities on the structural details of the complex and only to a lesser extent on experimental error.
Keywords
- Support Vector Regression
- Molecular Descriptor
- Partial Little Square Regression
- Outer Loop
- Principal Component Regression
Background
Cyclodextrins (CDs) are cyclic oligomers of α-D-Glucose, which can be categorised into four types: α-, β-, γ- and δ-CDs, corresponding to 6, 7, 8 or 9 α-D-glucose units. The shape of CDs has been described as torus- or doughnut-like, reflecting the existence of a cavity within the molecule. The exterior region of the CD, which is populated with hydroxyl groups, is hydrophilic, whereas the cavity is dominated by hydrophobic interactions, enabling CDs to form relatively strong complexes with hydrophobic guests. The hydrophobicity of the CDs' exterior ensures their water solubility, a property that results in their being potential solubilisers [1].
The size of the molecules that can be bound by a particular CD is related to the size of the cavity. α-CDs bind alkyl-chains of various lengths, whilst benzene, for instance, is seen to be too large. β-CDs cavities can complex more bulky molecules, such as adamantane, naphthalene or various benzene derivatives. γ-CDs can bind annelated ring systems and even buckyballs up to C_{60}. The ability of CDs to bind molecules of particular sizes has been termed 'size recognition' [1, 2].
β-CDs in particular have proven to be in high demand in the pharmaceutical industry because their cavities seem to be almost 'predestined' to bind drug-size molecules. Several formulations are on the marketplace in which β-CDs are applied as solubilisers of insoluble drugs [3–7]. Other industrial applications have been reported in the food industry, where CDs have been used to protect flavours or vitamins from oxidation [8]. It is therefore felt that an ability to predict the thermodynamic properties of this particular host-guest system could potentially have a direct impact on the development of novel pharmaceutical formulations.
Some time has been spent on studying and predicting the binding free energies (ΔG°) of CD inclusion complexes using computational methods [9, 10]. Amongst these statistical methods, those based on multiple regression [11, 12] or neural nets [13] have particularly proven to lead to robust prediction models.
In this study we investigate the predictability of the thermodynamic quantities relevant to complex formation, i.e. the free energy change ΔG°, the change of enthalpy ΔH° and the change of entropy ΔS°. The study is combined with a detailed performance comparison of three different types of statistical regression methods, namely principal components regression (PCR) [14], partial least squares regression (PLSR) [14] and support vector regression with forward feature selection (SVMR/FFS) [15, 16]. Whereas the first two methods are well established in the field of cheminformatics, the latter is a relatively new machine learning technique, which has been successfully applied in recent research projects [17, 18]. We have also reported its application as a valuable tool for increasing hit rates in similarity based virtual screenings [19].
Our study shows that the ease of predicting ΔG° cannot be explained by the predictability of either ΔH° or ΔS°. We shall discuss this finding in the context of a concise analysis of the experimental accuracy of thermodynamic data.
Results and discussion
Comparison of the regression methods for ten-fold cross validation. The maximal q^{2} values are reported for each thermodynamic parameter.
ΔG ° | ΔH ° | TΔS ° | |
---|---|---|---|
q^{2}(max) | q^{2}(max) | Q^{2}(max) | |
PCR | 0.71 | 0.54 | 0.35 |
PLSR | 0.74 | 0.53 | 0.31 |
SVMR/FFS | 0.89 | 0.75 | 0.63 |
Comparison of the regression methods
where σ^{2}(...) = variance of the respective quantity in brackets, Δy is the deviation between prediction and experimental value, and y is the quantity being predicted.
On applying PCR to predict ΔG°, ΔH° and TΔS°, the highest cross-validation q^{2} values obtained are 0.71, 0.54 and 0.35 respectively (Table 1). PLSR leads to models with maximal q^{2} values for the three properties of 0.74, 0.53 and 0.31, respectively. The highest q^{2} values obtained are for SVMR/FFS, namely, 0.89, 0.75 and 0.63 respectively.
In order to validate accurately the statistical models we performed the nested cross-validation protocol as described by Ruschhaupt et al. [20]. This type of validation gives an accurate estimate of the reliability of predicting external data. The method consists of inner and outer cross-validation loops. In the inner loop 10-fold cross-validation is used to identify the optimal settings for the complete algorithm that are then used for the outer loop. In the outer loop the performance of the best model obtained from the inner loop is tested on unseen data. Thus the test-sets in the outer loop do not in any way influence the training, with the performance of the learner on these test sets being an objective measure of its expected performance on similar data. Compared to the typical validation by external test-sets, this approach has an advantage of being less dependent on the partitioning into test and training sets, as each data point is part of the test-set exactly once. For each regression method this procedure was performed three times resulting in nine different models and prediction assessments. A more detailed description of nested cross validation is outlined in the methods section.
Comparison of the regression methods for nested cross-validation (ΔG°). Shown are the maximal q^{2} in the inner loop (q^{2}(max)-inner loop), the maximal q^{2} in the outer loop (q^{2}(max)-outer loop) and the q^{2} of the outer loop predicted by the model with the maximal q^{2} in the inner loop (q^{2} (inner loop-max)-outer loop).
Regression method | q ^{2} (max)-inner loop | q ^{2} (max)-outerloop | q ^{2} (inner loop-max)-outer loop |
---|---|---|---|
PCR | 0.71 ± 0.03 | 0.7 ± 0.03 | 0.69 ± 0.03 |
PLSR | 0.71 ± 0.03 | 0.7 ± 0.01 | 0.69 ± 0.03 |
SVMR/FFS | 0.87 ± 0.03 | 0.74 ± 0.01 | 0.71 ± 0.03 |
Comparison of the regression methods for nested cross-validation (ΔH°) with respect to q^{2}. Shown are the maximal q^{2} in the inner loop (q^{2}(max)-inner loop), the maximal q^{2} in the outer loop (q^{2}(max)-outer loop) and the q^{2} of the outer loop predicted by the model with the maximal q^{2} in the inner loop (q^{2}(inner loop-max)-outer loop).
Regression method | q ^{2} (max)-inner loop | q ^{2} (max)-outerloop | q ^{2} (inner loop-max)-outer loop |
---|---|---|---|
PCR | 0.49 ± 0.07 | 0.48 ± 0.02 | 0.48 ± 0.02 |
PLSR | 0.5 ± 0.08 | 0.47 ± 0.02 | 0.47 ± 0.02 |
SVMR/FFS | 0.73 ± 0.08 | 0.43 ± 0.02 | 0.4 ± 0.08 |
Comparison of the regression methods for nested cross-validation (TΔS°) with respect to q^{2}. Shown are the maximal q^{2} in the inner loop (q^{2}(max)-inner loop), the maximal q^{2} in the outer loop (q^{2}(max)-outer loop) and the q^{2} of the outer loop predicted by the model with the maximal q^{2} in the inner loop (q^{2}(inner loop-max)-outer loop).
Regression method | q ^{2} (max)-inner loop | q ^{2} (max)-outerloop | q ^{2} (inner loop-max)-outer loop |
---|---|---|---|
PCR | 0.33 ± 0.09 | 0.3 ± 0.05 | 0.3 ± 0.03 |
PLSR | 0.32 ± 0.1 | 0.29 ± 0.04 | 0.29 ± 0.04 |
SVMR/FFS | 0.64 ± 0.08 | 0.26 ± 0.04 | 0.21 ± -0.09 |
Predictability of different thermodynamic quantities
We also analyzed differences in the thermodynamic properties of structurally related guest molecules in order to obtain a more detailed picture of the reasons behind the poor predictability of ΔH° and TΔS°. To extend the experimental data-basis of our analysis, we integrated additional data if multiple measurements for a guest molecule were listed in Rekharsky's review [23]. For those compounds for which we had independent data from other published studies, we calculated the standard deviations for ΔG°, ΔH° and TΔS°, and averaged these over all compounds; the respective values obtained being 1.8 kJmol^{-1}, 2.1 kJmol^{-1}, and 2.7 kJmol^{-1}. These values describe the average absolute error in the experimental determination of the thermodynamic parameters. Interestingly, the magnitudes are entirely consistent with those obtained using the usual practice of determining entropy changes, that is, from the difference between the measured change in enthalpy and the measured change in the binding free energy. If we assume independent errors in the two latter quantities, we can calculate the expected error in the entropy change by means of the law of error propagation – the root of the sum of the squares of the errors in enthalpy and free energy is 2.8 kJmol^{-1}.
It is noteworthy that the magnitudes of the experimental errors found here are higher than those generally reported. This is mainly because our data include systematic errors arising from the compilation of data from different laboratories, whose experimental protocols most likely differ. The error values certainly agree with the predictability of the three quantities. However, the error is rather low when compared to the overall spread of the corresponding quantities: the overall standard deviations of the thermodynamic parameters in our dataset are 5.3 kJmol^{-1}, 9.6 kJmol^{-1}, and 8.5 kJmol^{-1} for ΔG°, ΔH° and TΔS° respectively. The average root mean square errors of the predicted to the experimental values obtained with SVMR/FFS are 2.8 kJmol^{-1} (ΔG°), 7.5 kJmol^{-1}(ΔH°) and 7.4 kJmol^{-1} (TΔS°). While the prediction of ΔG° appears to be limited mainly by experimental error (prediction error 2.8 kJmol^{-1} compared to the experimental error 1.8 kJmol^{-1}), ΔH° and TΔS° are clearly poorly predicted, which cannot be explained by the slightly higher values of the experimental error alone.
In addition, we attempted a nearest-neighbour prediction of ΔG°, ΔH° and TΔS° using the graph-based similarity of the molecules [31]. This method is independent of the E-Dragon descriptors [29] and regression methods. For each molecule within the dataset the three thermodynamic quantities were predicted to be equal to those of the most similar compound within the set. We obtained squared correlation coefficients r^{2} values of 0.50 for ΔG°, 0.47 for ΔH° and 0.29 for TΔS°. Except for a loss of accuracy in the prediction of ΔG°, the results are very similar to those obtained from the regression-based prediction. The main trend in the predictability of the thermodynamic quantities observed in the regression analysis can also be observed in this analysis, where again TΔS° proves to be the least predictable thermodynamic parameter.
This analysis indicates that the lower ability to predict TΔS° – and to a lesser extent, ΔH° – for different ligands has to do with the more complex dependence of TΔS° on even small structural changes in the ligand. This explanation is also consistent with the empirical observation of enthalpy-entropy compensation. The relative insensitivity of ΔG° to small structural changes compared to the other two quantities would lead to the compensation effects in enthalpy and entropy according to ΔG° = ΔH°-TΔS°, and conversely, given entropy-enthalpy compensation, changes in entropy would lead to smaller changes in free energy.
Conclusion
In this study we investigated the predictability of three important thermodynamic quantities, namely, the free energy of binding, heat of formation and the entropy change upon binding. To this end, we chose β-cyclodextrin with its ligands, a very well-studied system for which there is a large amount of high quality binding data available. We were able to show that free energies of binding can be reliably predicted by means of simple, readily available molecular descriptors with all three of the linear regression methods studied. The SVMR/FFS method has the advantage of leading to a (partly) interpretable model with comparably few descriptors. However, in the application of SVMR/FFS it is important to perform a nested cross-validation in order to obtain a realistic impression of its generalisation ability. The predictability of ΔG° obviously cannot be compared directly to that of ΔH°, as the latter is reproduced with significantly lower accuracy by the models analyzed. We found that TΔS° appears almost unpredictable, with an analysis of our results in the context of further data from the literature suggesting that its poor predictability – and, to a lesser extent, that of ΔH° – is explained by a stronger dependence of those quantities on the complex's structural details, and to a lesser extent on the wider experimental error. This would also explain the well-documented empirical finding of entropy-enthalpy compensation. In this sense our conclusion is in disagreement with that of Sharp, which suggested that entropy-enthalpy compensation is most likely due to lower accuracy in the experimental determination of binding enthalpy and entropy.
Methods
Dataset choice and preparation of the molecules
We assembled a dataset consisting of 176 β-CD ligands (see additional file 1; Table 1). These molecules are a subset of those collected by Rekharsky et al [23]. We applied the following selection criteria:
• The availability of experimental data derived from either calorimetric (cal) or UV-spectroscopic measurements;
• The availability of ΔG°, ΔH° and TΔS° data;
• The exclusion of all guest molecules whose data deviated from measurements of other groups.
We drew two-dimensional Lewis structures of the molecules with ISIS-Draw and exported them as MDL MOL files [24]. The protonation state of each molecule was manually set according to the pH value at which the measurement was performed. When no pH data were available, a reasonable state was set. To generate three dimensional low-energy structures from the mol file we used CORINA [25]. We then converted the structures to SD-files [26]. Finally, all structures were manually inspected and, when needed, corrected.
Calculation and processing of molecular descriptors
We calculated molecular descriptors for all molecules using the web service E-Dragon, which is part of the Virtual Computational Chemistry Laboratory [27, 28]. E-Dragon can calculate up to 1,666 different molecular descriptors [29], which are grouped into different categories ranging from simple atom-type descriptors or fragment counts to more sophisticated topological, geometrical or quantum chemical descriptors. In order to prevent numerical problems and to ensure the avoidance of any bias in the descriptor space, we normalised all descriptor values to a range between -1 and +1.
Regression Methods
The statistical methods used in this work have been employed on numerous occasions elsewhere. For PCR and Partial PLSR the R-package PLS was used [14]. The support vector machine regression was performed with LIBSVM, which was developed by Chang et al [30].
Principal component regression
In PCR a multiple linear regression is performed on principal components. Principal components are linear combinations of the descriptors in the data matrix and explain their variance. They are derived from the covariance matrix of the calculated descriptors. The number of principal components corresponds to the data matrix rank. Its maximal value is the minimum of the number of data points (i.e. molecules) and descriptors.
The first principal component of a data matrix points in the direction that maximizes the variance of the descriptors and corresponds to the largest eigenvalue of the covariance matrix. The second corresponds to the second largest eigenvalue and points in the direction that maximizes the variance and is orthogonal to the first principal component, and so on for the remaining principal components. The PCR model is generated on a subset of the components. The subset is built by selecting the components in the order of their ability to explain the variance in the dependent variable, i.e. in this study, the thermodynamic properties.
Partial Least Squares Regression
PLSR is very similar to PCR, however, while the covariance matrix of the data is used to generate the principal components, in PLSR the principal components are derived from the cross-covariance between the data matrix and the dependent variables. Hence, while in PCR the eigenvectors of the data covariance matrix are used to span the solution space, in PLSR the directions of maximal covariance between data and the dependent variables are used.
Support vector machine regression combined with forward feature selection
The theoretical background of support vector machine regression has been described in detail by Drucker et al [15]. Support vector regression is a straightforward variant of support vector machine (SVM) classification [16]. In classification problems SVMs find the hyper plane that separates positive examples from negative examples with a maximum margin (where the margin is defined as the distance of the closest data point from the separating hyper plane). In this way a statistical model is produced that only depends on a subset of the training data, namely those data points that are close enough to influence the size of the margin and the orientation of the hyper plane. These are the most difficult examples in the training set. They are termed 'support vectors' because they define the orientation of the separating plane. In support vector regression (SVMR) the same effect (namely that the final model depends only on a subset of the data) is achieved by the use of a so-called 'ε-insensitive cost function', which during model optimization ignores errors up to a defined threshold. This means that any training data being predicted by the current model with an accuracy of up to ε can be neglected.
In this work we added a 'forward feature selection procedure', which is in some respects similar to the component extension in PCR and PLSR. Forward feature selection increases the learning performance and the interpretability of the regression model as only descriptors are selected that significantly improve the SVMR model. The selection of descriptors produces combinatorial explosion if all possible subsets of all the available descriptors have to be considered. This, of course, is not feasible if the number of descriptors is too large. To overcome this problem, forward feature selection uses the following greedy heuristic, that is, for each single descriptor a support vector regression model is trained with tenfold cross-validation. The descriptor leading to the model with the highest q^{2} value is selected as the start descriptor. The procedure is repeated to find the next descriptor to form the best pair, triplet etc., iteratively expanding the model by one descriptor until q^{2} reaches a maximum, at which point the final model is obtained.
Validation by means of nested cross-validation
In order to validate whether our model generation procedures can lead to a predictive model that provides reliable output, we performed a nested three-way cross-validation protocol as proposed by Ruschhaupt et al [20] for each of the regression methods. To this end, we first split the training set into three equally sized subsets by randomly assigning training dataset molecules to one of the subsets (S_{1} and S_{2} consist of 59 molecules, S_{3} consists of 58 molecules). We then generated three validation sets, each as a combination of two subsets (V_{1}→ S_{1} and S_{2}, V_{2}→ S_{1} and S_{3}, V_{3}→ S_{2} and S_{3}), such that each of the validation sets could be used as a training-set for predicting the binding energies of the remaining subset (test set) that is not included in the respective training set. A ten-fold cross validation is used within the validation set (e.g. V_{1}) to identify the optimal model (e.g. the number of components in PCR/PLS or descriptors in SVMR), meaning that V_{1} is again separated into ten subsets that are used for cross validation in a ten-fold loop. The model resulting from this inner loop is then used to predict the test set (e.g. S_{3}). This is performed once for every validation set/test set combination, three times in total. Thus the reported results on the test sets from the outer loop do not in any way influence the choice of model parameters and are comparable to independent test set validations. The advantage of this approach over independent test set validation, is that every data point is predicted once in one of the three test sets, thus reducing the effects of test set choice.
Calculation of molecular similarity and clustering of the molecules
To cluster the dataset molecules and for the nearest-neighbour prediction, we calculated all pair-wise molecular similarities using our in-house similarity tool GMA [31]. The molecular similarity was calculated on the basis of a graph-based alignment. The better the molecular graphs (i.e. the topology and the atom types) of two molecules can be matched, the greater the similarity between these two molecules (1 = identical, 0 = dissimilar). On the basis of these similarities we performed complete-linkage hierarchical clustering. The cluster tree was cut off at a similarity threshold of 0.7. Hence, within one cluster only those molecules that exhibit a similarity of 0.7 or higher are grouped (see additional file 1, Table 2).
Declarations
Acknowledgements
We are grateful to Deutsche Forschungsgemeinschaft for funding part of this work (grants AP-101/1 and KA-1804/1). The authors wish to thank Joachim Büch for technical support, and Thomas Lengauer for critically reading an early version of the manuscript.
Authors’ Affiliations
References
- Wenz G: Cyclodextrins As Building-Blocks For Supramolecular Structures And Functional Units. Angew Chem-Int Edit Engl. 1994, 33 (8): 803-822. 10.1002/anie.199408031.View ArticleGoogle Scholar
- Muller A, Wenz G: Thickness recognition of bolaamphiphiles by alpha-cyclodextrin. Chemistry. 2007, 13 (8): 2218-2223. 10.1002/chem.200600764.View ArticleGoogle Scholar
- Davis ME, Brewster ME: Cyclodextrin-based pharmaceutics: past, present and future. Nat Rev Drug Discov. 2004, 3: 1023-1035. 10.1038/nrd1576.View ArticleGoogle Scholar
- Fugen G, Cuijing L: The preparation of inclusion compound of diclofenac sodium-[beta]- cyclodextrin. Chinese Pharmaceutical Journal. 1998, 33 (3): 153-Google Scholar
- Kang JC, Kumar V, Yang D, Chowdhury PR, Hohl RJ: Cyclodextrin complexation: influence on the solubility, stability, and cytotoxicity of camptothecin, an antineoplastic agent. Eur J Pharm Sci. 2002, 15 (2): 163-170. 10.1016/S0928-0987(01)00214-7.View ArticleGoogle Scholar
- Stanton J, Vincent P: Tumor necrosis factor receptor 2 . 2001, USA , Nuvelo, Inc.Google Scholar
- Stuerzebecher CS, Witt W, Raduechel B, Skuballa W, Vorbrueggen H: Prostacyclins, their analogs or prostaglandins and thromboxane antagonists for treatment of thrombotic and thromboembolic syndromes. 1996, USA , Schering AktiengeselleschaftGoogle Scholar
- Szejtli J, Bolla-Pusztai E, Szabo P, Ferenczy T: Enhancement of stability and biological effect on cholecalciferol by beta-cyclodextrin complexation. Pharmazie. 1980, 35 (12): 779-787.Google Scholar
- Connors KA: Prediction of binding constants of alpha-cyclodextrin complexes. J Pharm Sci. 1996, 85 (8): 796-802. 10.1021/js960167j.View ArticleGoogle Scholar
- Lipkowitz KB: Applications of computational chemistry to the study of cyclodextrins. Chem Rev. 1998, 98 (5): 1829-1873. 10.1021/cr9700179.View ArticleGoogle Scholar
- Suzuki T: A nonlinear group contribution method for predicting the free energies of inclusion complexation of organic molecules with alpha- and beta-cyclodextrins. J Chem Inf Comput Sci. 2001, 41 (5): 1266-1273. 10.1021/ci010295f.View ArticleGoogle Scholar
- Katritzky AR, Fara DC, Yang HF, Karelson M, Suzuki T, Solov'ev VP, Varnek A: Quantitative structure-property relationship modeling of beta-cyclodextrin complexation free energies. J Chem Inf Comput Sci. 2004, 44 (2): 529-541. 10.1021/ci034190j.View ArticleGoogle Scholar
- Liu L, Guo QX: Wavelet neural network and its application to the inclusion of beta-cyclodextrin with benzene derivatives. Journal Of Chemical Information And Computer Sciences. 1999, 39 (1): 133-138. 10.1021/ci980097x.View ArticleGoogle Scholar
- Mevik RWBH: pls: Partial Least Squares Regression (PLSR) and Principal ComponentRegression (PCR). 2007, [http://mevik.net/work/software/pls.html]R package version 2.0-1,Google Scholar
- Drucker H, Burges CJC, Kaufman L, Smola A, Vapnik V: Support Vector Regression Machines. 1996, MIT Press, 9: 155-161.Google Scholar
- Cortes C, Vapnik V: Support-Vector Networks. Machine Learning. 1995, 20 (3): 273-297.Google Scholar
- Briem H, Gunther J: Classifying "kinase inhibitor-likeness" by using machine-learning methods. Chembiochem. 2005, 6 (3): 558-566. 10.1002/cbic.200400109.View ArticleGoogle Scholar
- Liu X, Lu WC, Jin SL, Li YW, Chen NY: Support vector regression applied to materials optimization of sialon ceramics. Chemometr Intell Lab Sys. 2006, 82 (1-2): 8-14. 10.1016/j.chemolab.2005.08.011.View ArticleGoogle Scholar
- Jorissen RN, Gilson MK: Virtual screening of molecular databases using a Support Vector Machine. J Chem Inf Model. 2005, 45 (3): 549-561. 10.1021/ci049641u.View ArticleGoogle Scholar
- Ruschhaupt M, Huber W, Poustka A, Mansmann U: A Compendium to Ensure Computational Reproducibility in High-Dimensional Classification Tasks. Statistical Applications in Genetics and Molecular Biology. 2004Google Scholar
- Sharp K: Entropy-enthalpy compensation: fact or artifact?. Protein Sci. 2001, 10 (3): 661-667. 10.1110/ps.37801.View ArticleGoogle Scholar
- Ross PD, Rekharsky MV: Thermodynamics of hydrogen bond and hydrophobic interactions in cyclodextrin complexes. Biophys J. 1996, 71 (4): 2144-2154.View ArticleGoogle Scholar
- Rekharsky MV, Inoue Y: Complexation Thermodynamics of Cyclodextrins. Chem Rev. 1998, 98 (5): 1875-1918. 10.1021/cr970015o.View ArticleGoogle Scholar
- MDL Information Systems I: MDL ISIS/Draw. MDL Information Systems, Inc., 1990-2002, 2.5Google Scholar
- Sadowski J, Gasteiger J: From Atoms and Bonds to Three-Dimensional Atomic Coordinates: Automatic Model Builders. Chem Rev. 1993, 93: 2567-2581. 10.1021/cr00023a012.View ArticleGoogle Scholar
- Dalby A, Nourse JG, Hounshell WD, Gushurst AKI, Grier DL, Leland BA, Laufer J: Description of several chemical structure file formats used by computer programs developed at Molecular Design Limited. J Chem Inf Comput Sci. 1992, 32 (3): 244-255.View ArticleGoogle Scholar
- Tetko IV, Gasteiger J, Todeschini R, Mauri A, Livingstone D, Ertl P, Palyulin V, Radchenko E, Zefirov NS, Makarenko AS, Tanchuk VY, Prokopenko VV: Virtual computational chemistry laboratory - design and description. J Comput-Aided Mol Des. 2005, 19 (6): 453-463. 10.1007/s10822-005-8694-y.View ArticleGoogle Scholar
- Tetko IV: Computing chemistry on the web. Drug Discov Today. 2005, 10 (22): 1497-1500. 10.1016/S1359-6446(05)03584-1.View ArticleGoogle Scholar
- Todeschini R, Consonni V: Handbook of Molecular Descriptors. Methods and Principles in Medicinal Chemistry. Edited by: Mannhold R, Kubyini H, Timmermann H. 2000, Weinheim , Wiley-VCH, 11:Google Scholar
- Chang CC, Lin CJ: LIBSVM : a library for support vector machines . 2001Google Scholar
- Marialke J, Korner R, Tietze S, Apostolakis J: Graph-Based Molecular Alignment (GMA). J Chem Inf Model. 2007, 47 (2): 591-601. 10.1021/ci600387r.View ArticleGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.