Cheminformatic models based on machine learning for pyruvate kinase inhibitors of Leishmania mexicana
© Jamal and Scaria; licensee BioMed Central Ltd. 2013
Received: 18 March 2013
Accepted: 1 October 2013
Published: 19 November 2013
Leishmaniasis is a neglected tropical disease which affects approx. 12 million individuals worldwide and caused by parasite Leishmania. The current drugs used in the treatment of Leishmaniasis are highly toxic and has seen widespread emergence of drug resistant strains which necessitates the need for the development of new therapeutic options. The high throughput screen data available has made it possible to generate computational predictive models which have the ability to assess the active scaffolds in a chemical library followed by its ADME/toxicity properties in the biological trials.
In the present study, we have used publicly available, high-throughput screen datasets of chemical moieties which have been adjudged to target the pyruvate kinase enzyme of L. mexicana (LmPK). The machine learning approach was used to create computational models capable of predicting the biological activity of novel antileishmanial compounds. Further, we evaluated the molecules using the substructure based approach to identify the common substructures contributing to their activity.
We generated computational models based on machine learning methods and evaluated the performance of these models based on various statistical figures of merit. Random forest based approach was determined to be the most sensitive, better accuracy as well as ROC. We further added a substructure based approach to analyze the molecules to identify potentially enriched substructures in the active dataset. We believe that the models developed in the present study would lead to reduction in cost and length of clinical studies and hence newer drugs would appear faster in the market providing better healthcare options to the patients.
Leishmaniasis is a tropical disease affecting 12 million people worldwide, with approximately ~2 million (1.5 million incidences of cutaneous leishmaniasis and 500,000 visceral leishmaniasis) new people getting infected each year . It is considered as one of the world’s most neglected disease given its strong association with poverty and limited resources invested in new tools for diagnosis, treatment, and control . Among tropical diseases, leishmaniasis ranks second as a causative factor in mortality and fourth in morbidity and has been reported to occur in as much as 88 countries. It affects massive populations in most tropical and subtropical regions resulting in a huge number of deaths. The disease has become a major threat to one-third of the world population with more than 90% of the cases arising out of India, Bangladesh, Sudan, South Sudan, Brazil and Ethiopia [3-5]. Leishmaniasis is caused by a Trypanosomatid protozoan parasite belonging to the genus Leishmania, which infect both human and domestic animals, resulting in significant social and economic losses, especially in developing nations . The infection spreads through the bite of the phlebotomine sandflies which injects the promastigotes into the host . Approximately 21 of 30 species cause infections in humans and include L. donovani complex with three species (L. donovani, L. infantum, and L. chagasi); the L. mexicana complex with four main species (L. mexicana, L. amazonensis, and L. venezuelensis); L. tropica; L. major; L. aethiopica; and the subgenus Viannia with four main species (L. (V.) braziliensis, L. (V.) guyanensis, L. (V.) panamensis, and L. (V.) peruviana) . L. Mexicana causes both cutaneous and diffused cutaneous types of infection . The disease is considered as a major constraint to economic development with symptoms ranging from self-healing ulcers to highly disfiguring lesions and serious, often lethal visceral diseases which affect the haemopoetic organs .
The therapy of Leishmaniasis has been quite a challenge given the fact that the commonly used drugs available for treatment are characterized by high toxicity, high costs, limited activity and considerable possibility of drug resistance [11, 12]. The first line drugs used in the therapy are antimonial compounds such as sodium stibogluconate and meglumine antimoniate which form the traditional therapy for leishmaniasis. They are administered through the parenteral route and have severe side effects. In case of failure of the first line drugs, second line drugs are used, which include pentamidine (Lomidine) and amphotericin B (Fungizone). However, both these drugs are also associated with high levels of toxicity and side effects. Similarly, the new drug, Miltefosine (Impavido) prescribed for visceral and cutaneous leishmaniasis has also been identified to cause adverse effects [13, 14]. It has a long residence time which may contribute to the selection of resistant parasites, limiting its applicability. Miltefosine stays for a longer duration time in circulation approximately 150 hours which may lead to the development of resistance owing to which the parasite spreads rapidly . The current approach based on chemotherapy relies on a handful of drugs which are limited by factors such as high costs, toxicity, difficult routes of administration, and less efficacy . Keeping in mind all these factors, it is necessary to develop reasonably priced, secure, and effective antileishmanial vaccines for the acceptable therapy of leishmaniasis.
In Leishmania sp., sugar uptake and gluconeogenesis are essential to synthesize hexose-phosphates necessary for the production of glycoconjugates and the intracellular polysaccharide mannan, which form essential components for both replication and virulence of the parasite . Therefore, glycolytic enzymes are extremely pertinent for the growth and infectivity of the parasite.
Glycolysis pathway and enzymes in the pathway has been extensively reviewed as a potential drug target candidates . One of the well characterized enzymes in the pathway is Pyruvate kinase, which also has been extensively studied as a candidate drug target. Trypanosomatids entirely depend on the carbon sources available inside the host to meet their energy requirements and the only source of ATP generation is glycolysis as they lack Kreb’s cycle. Pyruvate kinase plays an important role in carbohydrate and amino acid metabolism and catalyses the last step in glycolysis to produce ATP and pyruvate kinase. Several reports have exploited the features of glycolytic enzymes on the basis of the enzyme’s structure further leading to the utilization of these features for the design of specific inhibitors [19-21]. Leishmania mexicana is known to encode for two copies of the enzyme, organized tandem to one another . The crystal structure of the enzyme has also been elucidated recently . The recent availability of high throughput screens for drug discovery of neglected diseases has motivated us to create predictive models based on molecular properties and machine learning approaches [23-26]. Recently a large dataset of high-throughput screens have been made available in public domain for Leishmania Mexicana Pyruvate Kinase and forms the baseline for the present study.
In the present study, we have used a computational strategy to create predictive classification models from the high-throughput assay which target pyruvate kinase enzyme from L. mexicana (LmPK). We have further analyzed chemical substructures to find enriched bioactive molecules using Maximum Common Substructure (MCS) approach and we also show that machine learning based cheminformatic modeling can create predictive models with high accuracy which can be effectively used to screen large molecular databases in silico, thus drastically reducing the cost of finding hits for drug discovery.
Bioassay and data sources
The assay used in the current study targets pyruvate kinase from Leishmania mexicana (LmPK). The datasets for the assay have been deposited at PubChem, a database collecting information on small molecules and datasets on high throughput biological assays and maintained by the National Centre for Biotechnology Information (NCBI) . The assay [AID: 1721] consisted of a total of 292,740 compounds capable of inhibiting the enzyme, pyruvate kinase, derived from Leishmania mexicana. Compounds were characterized based on a PubChem activity score. The compounds that had an activity score between 40 and 100 were defined as active (N = 1,087) and the compounds having an activity score of 0 were defined as inactive (N = 289,657). All the compounds having activity score between 1 and 39 were considered to have inconclusive activity and were not included in our analysis in order to avoid uncertainty in the predictive ability of the models.
Dataset pre-processing and calculation of molecular descriptors
The chemical structure of each of the molecules was downloaded in the Structural Data Format (SDF) from PubChem. These structures were imported into the molecular descriptor generator and visualization software PowerMV . PowerMV generates 2D molecular descriptors and is freely available. The dataset was split into smaller SDF files using SplitSDFiles Perl script available from Mayachem tools . A total of 179 descriptors were computed for the molecules. These descriptors encompassed different categories and included 147 descriptors which were pharmacophore fingerprints, 24 descriptors which were weighted burden numbers and 8 which belonged to property descriptors. The bit-string fingerprint attributes of only one value (all 0’s or all 1’s) all across the molecules were removed to reduce the dimensionality of the dataset. The full set of compounds were randomly divided into 20% independent test set and 80% training cum validation set using a bespoke Perl script. We used 5-fold cross validation in our study.
Machine Learning methods and implementation
Machine learning is a scientific discipline that broadly refers to a collection of algorithms and computational methods for predictive learning from tagged datasets . In cheminformatics, such methods have been extensively used to predict molecular properties, or biological activities. Generally, molecular datasets are tagged on the basis of their activity; say active/inactive and binary classification based on a set of molecular descriptions could be attempted. We have earlier shown that such an approach could accurately predict the activities in diverse sets of datasets with activities as diverse as anti-tubercular [23, 24] molecules, anti-malarial molecules  and RNA-binders . Similar tagging and learning could be attempted for multiple classes, rather than binary sets and have been extensively reviewed in . Multiple algorithms and implementations have been used in the area previously; nevertheless we attempted four popular classifier algorithms, that is, Naïve Bayes, Random Forest, J48 and SMO. All four methods have been previously determined to be quite efficient in terms of both computation time and classification accuracies. The Naive Bayes classifier is based on the Bayesian theorem, which assumes that for a given target value, the description of each predictor is independent of the other predictors. The final prediction is obtained by considering all descriptor-based properties . Random Forest algorithm is based on decision trees, where each tree is independently constructed and each node is split using the best among a subset of predictors randomly chosen at the node. It is the most accurate classifier and produces most precise results for all the datasets . J48 is a version of an earlier algorithm, the very popular C4.5, developed by J. Ross Quinlan and employs a tree pruning approach which produces fewer but more easily interpreted results. The J48 algorithm chooses one attribute of the data and splits the set of samples into subsets, one for every value of the attribute. The attribute having the maximum information gain is chosen to make the decision . Sequential Minimization Optimization (SMO) algorithm developed by John Platt in 1998 is widely used for training support vector machines. SMO, an iterative algorithm, breaks up the quadratic programming (QP) optimization problem into smaller problems which are then solved analytically. The SMO algorithm is simple, easy to use and faster in comparison to the standard SVM training algorithm .
Cost sensitive classification
One of the key issues that needs to be taken into consideration while using machine learning technique on a highly imbalanced dataset is the cost of misclassification. This is an important issue because standard classifiers presume equal weighing for all the classes and thus are unable to handle imbalanced data . The use of cost-sensitive classifiers can abrogate this issue and minimize misclassification errors. In cost sensitive learning, misclassification costs are used in which molecules are predicted to have the class with lower expected cost .
In the present study, we have used Weka (Waikato Environment for Knowledge Analysis), a collection of machine learning algorithms, for data mining tasks . Weka uses a confusion matrix consisting of four sections: True Positives (TP) for correctly classified actives; False Positives (FP) for inactive classified as actives; True Negatives (TN) inactive classified as inactive and false negatives (FN) for active compounds incorrectly classified as inactive. One of the most important points to be well thought-out during the development of classifiers is the fact that the false negatives are considered to be more important than the false positives. Consequently, we can minimize the % of false negatives at the expense of increasing false positives. To keep a check on the rate of false positives, an upper limit of 20% is set on the false positives. In Weka, there are no rules to set misclassification cost. It exclusively depends on the base classifier used .
Statistical measures for evaluation of cheminformatics models
A variety of measures were used to evaluate the performance of models such as sensitivity, specificity, accuracy, and BCR. Sensitivity (TP/(TP + FN)) is the proportion of positively labeled molecules predicted correctly. Specificity (TN/(TN + FP)) is the percentage of negatively labeled instances predicted as negative. Accuracy ((TP + TN)/(TP + TN + FP + FN) * 100) is the percentage coverage of correct predictions. Balanced Classification Rate (BCR) (½. (Sensitivity + Specificity)) is the mean of sensitivity and specificity which introduces a balance amid the classification rate of the two classes. Matthews Correlation Coefficient (MCC) is regarded as the balanced measure that measures the quality of a binary classification. We also evaluated the models based on the Receiver Operating Characteristic (ROC) curve which is the plot between the true positive rate and false positive rate.
Evaluation of enriched substructures
We used a hierarchical clustering algorithm ‘LibMCS’ , available from Chemaxon to find out potentially enriched molecular substructures in bioactive molecules . The maximum common substructure (MCS) based approach retrieves and compares the substructure common to a group of molecules. The MCS size which corresponds to the number of constituent atoms was set to an empirical threshold of 14 atoms in this study. The scaffolds obtained were then used to search for similar molecules in active and inactive datasets using the ‘jcsearch’ algorithm available from Chemaxon . The Chi-square test and the associated p-value were used for the further evaluation of substructures and to test the significance of enrichment, respectively. We calculated the enrichment factor and used a threshold of 2 to prioritize the molecules for further analysis. Using vROCS (release 3.1.2) , we further performed a molecular alignment of the selected scaffolds with the molecules in the active dataset and visualized the alignment in VIDA software  available from Open Eye Scientific Software .
Results and discussion
Significantly enriched substructures in the active dataset
0.00E + 00
0.00E + 0
The recent years have seen a wealth of information being available in public domain on molecular structures and biological assays of very small subsets of known small molecule repertoire using high-throughput screening platforms. The major challenge within the field pertains to the assigning potential biological activities to molecules so that they could be studied in detail. On an average less than a percent of the molecule library screened show some biological activity. Given the large costs associated with setting up screens for biological activities, it would be economically not plausible to exhaustively screen large parts of the known small molecule repertoire known to mankind. This problem becomes acute in cases of neglected tropical diseases. The challenge therefore would be to effectively mine large libraries using computational tools so they could be effectively prioritized for experimental screening for their biological activities. This necessitates the creation of highly accurate computational methods capable of predicting biological activities in silico. In the present study, we use machine learning as an approach to build highly accurate predictive models for bio-activity against pyruvate kinase on Leishmania species. We show how high-throughput experimental datasets on a diverse set of molecules could potentially be used to build highly accurate predictive models. These models could potentially be used to mine and annotate large molecular datasets and prioritize molecules for biological activity screening experiments and could contribute significantly to the ongoing efforts for drug discovery for neglected tropical diseases.
Leishmaniasis is one of the major neglected tropical diseases in recent years, killing close to 100,000 individuals worldwide annually, mostly in the tropical and sub-tropical countries. The disease is majorly distributed in the tropical and sub-tropical regions. Though efficient treatment regimens are available for its therapy, the drugs used are largely toxic. In addition, wide-spread drug resistance has been reported in several regions, adding to the urgency for the discovery of novel, efficient and less-toxic molecules with anti-leishmanial activity. In this study, we have employed computational strategy via machine learning approach to create predictive models for classification of molecules to discover new therapeutic compounds for leishmaniasis. The approach will help in quick search of large libraries of chemical structures in order to pick potential hits which are most probable to bind to a drug target. Additionally, we have used a substructure based approach to explore potentially enriched substructures in the active dataset of molecules. We show that accurate models for mining large datasets could be built based on high throughput assays available in public domain using machine learning approaches. We have previously reported similar approaches for mining molecules with anti-tubercular activities, and suggest that Random Forest based learning can systematically learn from bioassay datasets with high accuracies. We argue that automated approaches based on Random Forest based algorithms could be implemented on a large-scale to learn from bioassay datasets and automatically annotate molecules from PubChem for biological activities. We hope that such approaches could accelerate the process and efficiency of screening for discovery of novel molecules with specific biological activities not just for tropical diseases but others as well.
The predictive model generated by Weka, a stepwise manual and the scripts to be used for preprocessing of the dataset are available online at http://vinodscaria.rnabiology.org/2C4C/models.
The authors thank Dr S Ramachandran and Dr Sridhar Sivasubbu for reviewing the manuscript and for scientific suggestions. We also acknowledge the scientific suggestions and support and insights provided by Ms Vinita Periwal, Institute of Genomics and Integrative Biology (CSIR-IGIB). The computation was supported by C-DAC India through the Garuda Grid and supported by the National Knowledge Network (NKN), India. Authors acknowledge the support provided by Mrs Mangala N, Janaki Chintalapati and Dr. Subrata Chattopadhyay from CDAC for setting up the cheminformatics test bed on Garuda. This work was funded by the Council of Scientific and Industrial Research (CSIR), India for funding through the Open Source Drug Discovery Project (HCP001).
Open Source Drug Discovery Consortium.
CSIR Open Source Drug Discovery Unit, Anusandhan Bhavan, New Delhi 110001, India.
- World Health Organization: http://www.who.int/leishmaniasis/burden/magnitude/burden_magnitude/en/index.html,
- Alvar J, Yactayo S, Bern C: Leishmaniasis and poverty. Trends Parasitol. 2006, 22 (12): 552-10.1016/j.pt.2006.09.004.View ArticlePubMedGoogle Scholar
- Yamey G, Torreele E: The world’s most neglected diseases. Br Med J. 2002, 325: 176-177. 10.1136/bmj.325.7357.176.View ArticleGoogle Scholar
- Mathers CD, Ezzati M, Lopez AD: Measuring the burden of neglected tropical diseases: the global burden of disease framework. PLoS Negl Trop Dis. 2007, 1: e114-10.1371/journal.pntd.0000114.PubMed CentralView ArticlePubMedGoogle Scholar
- Alvar J, Ve’lez ID, Bern C, Herrero M, Desjeux P: Leishmaniasis worldwide and global estimates of its incidence. PLoS ONE. 2012, 7 (5): e35671-10.1371/journal.pone.0035671.PubMed CentralView ArticlePubMedGoogle Scholar
- World Health Organization: Thirteenth Programme Report, UNDP/World Bank/World Health Organization Programme for Research and Training in Tropical Diseases. 1997, Geneva: World Health OrganizationGoogle Scholar
- Hepburn NC: Cutaneous Leishmaniasis: an overview. J Postgrad Med. 2003, 49: 50-10.4103/0022-3859.928.View ArticlePubMedGoogle Scholar
- Banuls AL, Hide M, Prugnolle F: Leishmania and the leishmaniases: a parasite genetic update and advances in taxonomy, epidemiology and pathogenicity in humans. Adv Parasitol. 2007, 64: 1-109.View ArticlePubMedGoogle Scholar
- Piscopo TV, Mallia AC: Leishmaniasis. Postgrad Med J. 2007, 83: 649-657.PubMed CentralView ArticlePubMedGoogle Scholar
- Kato H, Gomez EA, Cáceres AG, Uezato H, Mimori T: Molecular epidemiology for vector research on leishmaniasis. Int J Environ Res Public Health. 2010, 7: 814-826. 10.3390/ijerph7030814.PubMed CentralView ArticlePubMedGoogle Scholar
- Reports of the World Health Organization. 2011, http://www.who.int/leishmaniasis/burden/en/,
- Croft SL, Coombs GH: Leishmaniasis- current chemotherapy and recent advances in the search for novel drugs. TRENDS in Parasitol. 2003, 19 (11): 503-508.View ArticleGoogle Scholar
- Kedzierski L, Sakthianandeswaren A, Curtis JM, Andrews PC, Junk PC, Kedzierska K: Leishmaniasis: current treatment and prospects for new drugs and vaccines. Curr Med Chem. 2009, 16 (5): 599-614. 10.2174/092986709787458489.View ArticlePubMedGoogle Scholar
- Lukasz K: Leishmaniasis vaccine: where are We today?. J Glob Infect Dis. 2010, 2 (2): 177-185. 10.4103/0974-777X.62881.View ArticleGoogle Scholar
- Maltezou HC: Drug resistance in visceral leishmaniasis. J Biomed Biotechnol. 2010, 2010: 617521-PubMed CentralView ArticlePubMedGoogle Scholar
- Kedzierski L, Zhu Y, Handman E: Leishmania vaccines: progress and problems. Parasitology. 2006, 133: S87-S112. 10.1017/S0031182006001831.View ArticlePubMedGoogle Scholar
- Chawla B, Madhubala R: Drug targets in Leishmania. J Parasit Dis. 2010, 34 (1): 1-13. 10.1007/s12639-010-0006-3.PubMed CentralView ArticlePubMedGoogle Scholar
- Rigden DJ, Phillips SE, Michels PA, Fothergill-Gilmore LA: The structure of pyruvate kinase from Leishmania mexicana reveals details of the allosteric transition and unusual effector specificity. J Mol Biol. 1999, 293 (3): 745-749. 10.1006/jmbi.1999.3170.View ArticleGoogle Scholar
- Fothergill-Gillmore LA, Rigden DJ, Michels PA, Phillips SE: The structure of pyruvate kinase from Leishmania mexicana reveals details of the allosteric transition and unusual effector specificity. J Mol Biol. 1999, 291 (3): 615-635. 10.1006/jmbi.1999.2918.View ArticleGoogle Scholar
- Hannaert V, Yernaux C, Rigden DJ, Fothergill-Gilmore LA, Opperdoes FR, Michels PA: The putative effector-binding site of Leishmania mexicana pyruvate kinase studied by site-directed mutagenesis. FEBS Lett. 2002, 514 (2-3): 255-259.View ArticlePubMedGoogle Scholar
- Verlinde CL, Hannaert V, Blonski C, Willson M, Périé JJ, Fothergill-Gilmore LA, Opperdoes FR, Gelb MH, Hol WG, Michels PA: Glycolysis as a target for the design of new anti-trypanosome drugs. Drug Resist Updat. 2001, 4 (1): 50-65. 10.1054/drup.2000.0177.View ArticlePubMedGoogle Scholar
- Morgan HP, McNae IW, Hsin KY, Michels PA, Fothergill-Gilmore LA, Walkinshaw MD: An improved strategy for the crystallization of Leishmania mexicana pyruvate kinase. Acta Crytallogr. 2010, 66: 215-218.Google Scholar
- Periwal V, Rajappan JK, Jaleel AU, Scaria V: Predictive models for anti-tubercular molecules using machine learning on high-throughput biological screening datasets. BMC Res Notes. 2011, 4: 504-10.1186/1756-0500-4-504.PubMed CentralView ArticlePubMedGoogle Scholar
- Periwal V, Kishtapuram S, Scaria V: Computational models for in-vitro antitubercular activity of molecules based on high-throughput chemical biology screening datasets. BMC Pharmacol. 2012, 12: 1-10.1186/1471-2210-12-1.PubMed CentralView ArticlePubMedGoogle Scholar
- Jamal S, Periwal V, Scaria V: Predictive modeling of anti-malarial molecules inhibiting apicoplast formation. BMC Bioinforma. 2013, 14: 55-10.1186/1471-2105-14-55.View ArticleGoogle Scholar
- Jamal S, Periwal V, Scaria V: Computational analysis and predictive modeling of small molecule modulators of microRNA. J Cheminformatics. 2012, 4: 16-10.1186/1758-2946-4-16.View ArticleGoogle Scholar
- Wang Y, Xiao J, Suzek TO, Zhang J, Wang J, Bryant SH: PubChem: a public information system for analyzing bioactivities of small molecules. Nucleic Acids Res. 2009, 37: W623-W633. 10.1093/nar/gkp456.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu K, Feng J, Young SS: PowerMV: a software environment for molecular viewing, descriptor generation, data analysis and hit evaluation. J Chem Inf Model. 2005, 45: 515-522. 10.1021/ci049847v.View ArticlePubMedGoogle Scholar
- Sud M: MayaChemTools. 2010, http://www.mayachemtools.org/.Google Scholar
- Melville JL, Burke EK, Hirst JD: Machine learning in virtual screening. Comb Chem High Throughput Screen. 2009, 12: 332-343. 10.2174/138620709788167980.View ArticlePubMedGoogle Scholar
- Lars J, Jasmin S, Peer B: Literature mining for the biologist: from information retrieval to biological discovery. Nat Rev Genet. 2006, 7 (2): 119-129. 10.1038/nrg1768.View ArticleGoogle Scholar
- Friedman N, Geiger D, GoldSzmidt M: Bayesian Net Classif Machine Learn. 1997, 29: 131-163.View ArticleGoogle Scholar
- Breiman L: Random forests. Machine Learn. 2001, 45: 5-32. 10.1023/A:1010933404324.View ArticleGoogle Scholar
- Quinlan JR: C4.5: programs for machine learning. 1993, San Francisco: Morgan Kaufmann PublishersGoogle Scholar
- Cortes C, Vapnik V: Support vector networks. Mach Learn. 1995, 20: 273-297.Google Scholar
- Japkowicz N: Proceedings of the International Conference on Artificial Intelligence. The Class Imbalance Problem: Significance and Strategies. 2000Google Scholar
- Elkan C: The foundations of cost-sensitive learning. Proceed Seventeenth Internl Joint Confer Artific Intellig. 2001, 2: 973-978.Google Scholar
- Bouckaert RR, Frank E, Hall MA, Holmes G, Pfahringer B, Reutemann P: Weka -experiences with a java open-source project. J Mach Learn Res. 2010, 10: 2533-2541.Google Scholar
- Schierz AC: Virtual screening of bioassay data. J Cheminform. 2009, 22: 1-21.Google Scholar
- Chemaxon: 2008, Budapest H. Library MCS, version 0.7. 2008
- Chemaxon: Budapest H. Jcsearch version 5.8.2
- vROCS, release 3.1.2, OpenEye Scientific Software, Inc: 2010, NM, USA: Santa Fe, http://www.eyesopen.com,
- VIDA, version 4.1.1, OpenEye Scientific Software, Inc: 2010, NM, USA: Santa Fe, http://www.eyesopen.com,
- OpenEye Scientific Software, Inc: 2010, NM, USA: Santa Fe, http://www.eyesopen.com,
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.