- Open Access
Prioritization of anti-malarial hits from nature: chemo-informatic profiling of natural products with in vitro antiplasmodial activities and currently registered anti-malarial drugs
Malaria Journal volume 15, Article number: 50 (2016)
A large number of natural products have shown in vitro antiplasmodial activities. Early identification and prioritization of these natural products with potential for novel mechanism of action, desirable pharmacokinetics and likelihood for development into drugs is advantageous. Chemo-informatic profiling of these natural products were conducted and compared to currently registered anti-malarial drugs (CRAD).
Natural products with in vitro antiplasmodial activities (NAA) were compiled from various sources. These natural products were sub-divided into four groups based on inhibitory concentration (IC50). Key molecular descriptors and physicochemical properties were computed for these compounds and analysis of variance used to assess statistical significance amongst the sets of compounds. Molecular similarity analysis, estimation of drug-likeness, in silico pharmacokinetic profiling, and exploration of structure–activity landscape were also carried out on these sets of compounds.
A total of 1040 natural products were selected and a total of 13 molecular descriptors were analysed. Significant differences were observed among the sub-groups of NAA and CRAD for at least 11 of the molecular descriptors, including number of hydrogen bond donors and acceptors, molecular weight, polar and hydrophobic surface areas, chiral centres, oxygen and nitrogen atoms, and shape index. The remaining molecular descriptors, including clogP, number of rotatable bonds and number of aromatic rings, did not show any significant difference when comparing the two compound sets. Molecular similarity and chemical space analysis identified natural products that were structurally diverse from CRAD. Prediction of the pharmacokinetic properties and drug-likeness of these natural products identified over 50 % with desirable drug-like properties. Nearly 70 % of all natural products were identified as potentially promiscuous compounds. Structure–activity landscape analysis highlighted compound pairs that form ‘activity cliffs’. In all, prioritization strategies for the NAA were proposed.
Chemo-informatic profiling of NAA and CRAD have produced a wealth of information that may guide decisions and facilitate anti-malarial drug development from natural products. Articulation of the information provided within an interactive data-mining environment led to a prioritized list of NAA.
Malaria is a major health burden in several developing countries and imposes a huge strain on health systems, particularly in Africa where over 50 % of the malaria deaths occur in children under the age of five and pregnant women [1, 2]. The emergence of drug-resistant Plasmodium falciparum strains, a major causative organism of malaria, has led to increasing numbers of fatal cases [1, 2]. Consequently there is an urgent need to discover or design new anti-malarial drugs with mechanism of actions that will circumvent the current resistance profile of P. falciparum. Arguably, natural products from plants (phytochemicals) have been the most consistent and successful source and template of anti-malarial drugs . Starting with quinine (from the bark of Cinchona)  to artemisinin (from Artemisia annua) , natural plant products have provided an invaluable armament against malaria infection. More promising is the fact that the literature revealed an increasing number of natural products, from ethnomedicine in malaria-endemic regions, with good in vitro and/or in vivo antiplasmodial activities [6–9]. Yet, many of these natural products have not made it into or made much progress down the anti-malarial drug development pipeline [10, 11]. Therefore there is a dire need to begin the process of identifying natural products with potential for anti-malarial drug discovery.
Perhaps a pertinent question to ask in light of the high cost, long duration and high failure rate of drug discovery [12–15] is: Should natural products with in vitro antiplasmodial activities (NAA) that are most likely to be successfully developed into anti-malarial drug candidates be prioritized? One approach that may help prioritize such natural products is chemo-informatics profiling. Chemo-informatics integrates chemical information with biological information  and translate such information into knowledge that could be used to assist decision making in the area of compound prioritization, selection, optimization, and ultimately clinical development . Chemo-informatics profiling of natural products with antiplasmodial activities, hereafter referred to as NAA, may allow researchers to prioritize and select NAA for the next stage of anti-malarial drug development. Conceivably, what may be more informative is to carry out such profiling with reference to currently registered anti-malarial drugs (CRAD). Currently registered anti-malarial drugs, which have successfully passed through all the drug development hurdles, have molecular descriptors and physicochemical properties to which NAA need to conform or, more interestingly, deviate.
Chemo-informatic profiling of NAA and CRAD was conducted. Specific approaches used included comparison of key molecular descriptors and physicochemical properties of NAA and CRAD, molecular similarity/diversity analysis, exploration of structure–activity landscape, estimation of drug-likeness, bioavailability, and toxicity profile. Literature search showed no report of chemo-informatic profiling of NAA in comparison with CRAD. The results from this analysis may provide insight into the important molecular features that define the reported in vitro antiplasmodial activities, potential for good bioavailability (which is essential for in vivo assay), toxicity liabilities and structural-activity relationships that may prioritize promising NAA. Such knowledge may expedite the progress of NAA for anti-malarial drug development.
Data collection and preparation of dataset
The trivial name, source, chemical class and antiplasmodial activity (IC50) of NAA were retrieved from published articles [6–9], MSc and PhD theses, textbook chapters, collaborative drug discovery (CDD) database , ChEMBL and PubChem (see Additional file 1). A total of 1040 NAA were selected based on availability of chemical structure and bioactivity data (C50) (see Additional file 2). The selected NAA were sub-divided into four categories based on a normalized IC50 (pIC50): highly active (HA) with IC50 less than 1 µM (pIC50 >0), active (A) with IC50 equal or greater than 1 µM but less than 5 µM (pIC50 ≤ 0 > −0.7), moderately active (MA) with IC50 equal or greater than 5 µM but less than 10 µM (pIC50 ≤ −0.7 but > −1) and low active(N) with IC50 equal or greater than 10 µM (pIC50 ≤ −1). The chemical structures of these NAA were downloaded from PubChem and ChEMBL databases in two dimensional (2D) SDF format. The NAA that were not found in public chemical databases were drawn with GchemPaint chemical structure editor for Linux and exported in molfile format. All chemical structures were combined and duplicates removed according to InChIKey generated by Open Babel . CRAD were retrieved from ChEMBL (name and smiles format). Three dimensional (3D) structures were generated for all the compounds from either 2D SDF, smiles or molfile formats (using builder module), corrected and minimized (using MMFF94 force field) with Molecular Operating Environment (MOE) 2013 software .
Calculation of molecular descriptors and physicochemical properties
The QuSAR module of the MOE package  was employed to calculate structure-related 2D molecular descriptors. Other physicochemical properties (e.g., ligand efficiency (LE), number chemical functional groups) were computed with ICM Chemist Pro (v3.7) from Molsoft Inc. and DataWarrior  running on a Linux platform on a Dell Vostro 2520 computer. Boxplots of the molecular descriptors and physicochemical properties were plotted for NAA (HA, A, MA, N) and CRAD using DataWarrior. The mean of the molecular descriptors and physicochemical properties for NAA (HA, A, MA, N) and CRAD were compared and statistical differences assessed with analysis of variance (ANOVA) with significance set at p < 0.05. Furthermore, the association between the in vitro antiplasmodial activities (pIC50) of NAA and the molecular descriptors and physicochemical properties were assessed using Spearman correlation coefficient (r).
Molecular similarity/diversity analysis
MoSS most common substructure (MoSSMCSS)
The KNIME (Konstanz Information Miner) workflow , shown in Fig. 1, was used to compute and visualize the molecular similarity, based on most common substructure (MoSSMCSS) , amongst the NAA and between NAA and CRAD. Molecular similarity within the dataset was visualized with a heat map.
Chemical space analysis (ChemGPS-NP)
The coordinates of chemical space occupied by NAA relative to those occupied by CRAD were explored with ChemGPS-NP . A text file containing the smiles and identifier for the compounds in the dataset was submitted to the ChemGPS-NP web service. The result, with four principal components added for each compound, was downloaded as a text file. The text file was then opened in DataWarrior and the four principal components plotted on a 3D scatter plot (markers were sized relative to one of the principal components).
Analysis of structure–activity landscape (identifying activity cliffs)
In recent time, structure–activity landscape analysis was used to visualize the relationship between molecular structure and bio-activities, especially in large activity screening datasets . The structure–activity landscape analysis was conducted on NAA using DataWarrior . Similarity/activity cliff analysis was chosen from the Chemistry menu. SkelSpheres descriptors were used as similarity criteria for arranging molecules on a 2D self-organizing maps (SOM). The activity column containing the numeric value of IC50 of the NAA was selected to calculate the structure–activity landscape index (SALI). For any couple of molecules, the SALI value reflects how much activity is gained with a small modification of the chemical structure. The hit status, i.e., NAA sub-group (A, HA, MA, N) was chosen as the Identifier column. Similarity limit was set at 0.8 (Tanimoto coefficient). Visualization based on similarity relationships were created (markers were coloured by activity).
Prediction of pharmacokinetic properties and drug-likeness of NAA and CRAD
Four models were used to characterize the bio-availability profile of NAA and CRAD namely: Lipinski’s rule of five , Egan Egg , golden triangle  and a model proposed by Veber et al. . For the Lipinski’s rule of five, a 3D plot was constructed in DataWarrior for NAA and CRAD using the number of hydrogen bond donors (HBD) and acceptors (HBA), calculated logarithm of partition coefficient (clogP) and molecular weight (MW). Markers on the plot were sized by molecular weight. The filters for the descriptors (in Lipinski’s rule of five) were then adjusted, in the DataWarrrior software , to the cut-off proposed by the Lipinski’s model (HBD <5, HBA <10, clogP <5, MW <500). The proportion of NAA and CRAD that fell within these limits were retained.
ICM Chemist Pro (v3.7) from Molsoft Inc. was used to spot compounds with potential for good, borderline and poor absorption based on the Egan Egg model . A plot of the clogP versus polar surface area (PSA) was also generated, using DataWarrior, to interactively visualize the proportion of the compounds in each sub-group of NAA and CRAD that fall within the Egan Egg (PSA <131.6 Å2, clogP < 5.88).
Johnson et al.  proposed a golden triangle, with a base between −2 and 5 for logD (pH 7.4) and peak at 500 Daltons, which enclose compounds with potential for good absorption and low clearance. This golden triangle was superimposed on a plot of logD (pH 7.4) versus MW for the sub-groups of NAA and CRAD. Compounds enclosed within the golden triangles were enumerated and recorded for each sub-group of NAA and CRAD.
For the model proposed by Veber et al. , a plot of PSA against number of rotatable bonds (NRB) for the sub-groups of NAA and CRAD was plotted in DataWarrior. The DataWarrior filters were then set to the cut-off proposed by the model (PSA [140 Å2] and NRB (10)). The proportion of compounds, in sub-groups of NAA and CRAD, within the area bounded by this cut-off were then tallied and noted.
The toxicity risk of NAA and CRAD was assessed with DataWarrior (from open molecules). The toxicity assessment was based on a search for reported toxic substructures or ‘toxicophores’ in the NAA. The collection of toxicophores in the software used for the prediction of toxicity (DataWarrior) were obtained by shredding compounds in the Registry of Toxic Effects of Chemical Substances (RTECS®). The toxicity information appearing in the Registry is derived from reports of the toxic effects of chemical substances. Registry of Toxic Effects of Chemical Substances (RTECS®) consists of tabulations of the lowest dose reported to have caused the listed toxic effect in the designated species (including mammalian cells) by the designated route of administration. The result was presented as a Table with each compound tag with none, low and high risk for mutagenic, tumorigenic and irritant toxicity class. Nasty or reactive chemical functions identified by the DataWarrior from each compound were also inserted into the table of results.
Fragment-based drug-likeness and quantitative estimate of drug-likeness(QED) were calculated with DataWarrior and ICM Chemist Pro (v3.7), respectively.
Frequent hitters (promiscuous compounds)
MedChem rules , a set of 275 rules, was also used to identify compounds that may interfere with bioassays to produce ostensible activity (false positives). ICM Chemist Pro (v3.7) software from Molsoft Inc. was used to flag such compounds.
In this study, the approach to prioritizing the selected 1040 published NAA included calculation of molecular descriptors, calculation of LE metrics, assessment of structural similarity, overview of structural-activity landscape, pharmacokinetic profiling, toxicity profiling, and identification of frequent hitters. All these were done in relation to CRAD.
Description of dataset
A description of the dataset used in this study is presented here. A total of 1040 NAA and 27 CRAD were included in the dataset. With regards to in vitro antiplasmodial bioactivity (IC50) described earlier in Methods, NAA consist of 21 % HA compounds, 28 % A compounds, 20 % MA, and 31 % N compounds. Concerning chemical class, a trend was observed between the prevalence of chemical class and bioactivity of NAA. The three most prevalent chemical class of compounds for each bioactivity category of NAA were: alkaloids (40 %), diterpenes (9 %) and flavonoids (9 %) for HA; alkaloids (34 %), terpenes (14 %) and flavonoids (13.5 %) for A; flavonoids (27 %), triterpenes (19 %) and alkaloids (15 %) for MA; triterpenes (30 %), alkaloids (25 %) and flavonoids (19 %) for N. The proportion of alkaloid compounds in NAA seems to decrease as the antiplasmodial potency decreases (i.e., most prevalent in HA bioactivity category). In contrast, percentage of flavonoids and the terpenoids (terpenes, diterpenes and triterpenes) appear to increase as the antiplasmodial activity decreases. This suggests that alkaloids may be an important chemical class for antiplasmodial activity.
Molecular descriptors and physicochemical properties of NAA and CRAD
In this section, the distribution and summary statistics of selected molecular descriptors and physicochemical properties of NAA and CRAD were determined and assessed as a contribution towards the prioritization of NAA.
Number of hydrogen bond acceptors and donors
HBA and HBD are essential to mediate interactions between compounds and biochemical macromolecules. They are also determinants of oral absorption of compounds . The results (Fig. 2, panel a) showed that the average number of HBA was similar for the A, HA and MA sub-groups of NAAs (median = 6). However, HBA was significantly lower in CRAD (median, 4) and low active sub-group of NAA (N) (median = 5 p < 0.05) in comparison to other sub-groups of NAA. In all cases the numbers of HBD (Fig. 2, panel b) were lower than the HBA. The HBD was significantly higher in the HA and MA than CRAD (p < 0.05). This observation aligned with other studies that report predominance of HBA and HBD in natural products over synthetic compounds [31–33]. The preponderance of strongly electronegative atoms (especially oxygen) in NAA may be responsible for the higher numbers of HBD and HBA [31–34]. In addition, there was very small negative correlation between the number of HBA (r = −0.20) and HBD (r = −0.10) and antiplasmodial activities (IC50) of NAA. This suggests the need for HBD and HBA for bioactivities. In conclusion, the results showed that the number of HBA and HBD was higher in the NAA compared to CRAD and showed slight association with in vitro antiplasmodial activities.
Total molecular weight (TMW)
A low MW (less than 500 Daltons) is usually favoured because of its impact on oral absorption of compounds . However, some natural products with MW >500 Daltons have been absorbed via the biological membrane . From the results, the average computed total molecular weight (TMW) was similar for HA (462.47), A (430.93) and MA (447.26) but significantly different (p value < 0.05) from N and CRAD (Fig. 2, panel c). This result attests to the greater structural complexity of NAA over CRAD. A slight negative correlation (r = −0.303) was observed between TMW and antiplasmodial activities of the NAA i.e., most active (HA) showed the highest TMW. Although it has been demonstrated that the median MW of oral drugs has increased substantially over the past years with about 2 % having MW >500 Daltons , a TMW below 500 may be desirable for hit compounds that are yet to be optimized into lead compounds.
Calculated logarithm of partition coefficient (clogP)
The clogP of a compound is a vital consideration for oral absorption  ] and it also influences binding to biological targets . From the results (Fig. 2, panel d), there was no significant difference in the clogP of CRAD and NAA (HA, A, MA, N). Given that the average clogP value observed for all the CRAD and up to 70 % of NAA were lower than five, as prescribed by Lipinski’s rule of five , it is expected that these NAA may show good oral absorption profile comparable to the CRAD. However Lipinski’s rule of five has been reported to have exceptions, predominantly in natural products [26, 38]. It is therefore expected that NAA with clogP > 5 may also have suitable oral absorption.
In addition, clogP have been reported to be predictive of bioactivity . Yet, correlation between clogP and the activity profiles of the NAA was negligible (r = −0.05). The inability of clogP to discriminate among the activity sub-groups of NAA highlights the contributions of other molecular properties to the observed antiplasmodial activity among the NAA. In view of the prospect for development into drug candidates, studies have shown that the clogP of molecules that make it all the way to commercialization have remained in the same range (2.6) for a number of years and that optimal range for clogP between one and three is vital for desirable physicochemical properties . From the results, though the average clogP of the NAA was similar to CRAD, it was however higher than the 2.6 reported for commercialized drugs [39, 40]. Given that the logP of hits and leads increase as they move down the drug development pipeline  and that high logP [>5] comes with some liabilities, e.g., receptor promiscuity [42, 43], poor metabolic clearance  and toxicity [44, 45], it is desirable to select NAA with minimal clogP values to ensure that eventual anti-malarial leads and drug candidates have clogP less than four  or less than five .
Number of rotatable bonds
A minimal NRB is desired for a NAA that may be prioritized and selected for development. The results (see Additional file 3) showed no significant difference (p > 0.05) in the average NRB between CRAD and the NAA. Although the NRB of HA (median, 3) was lower than that detected for A (median, 4) and MA (median, 4), no significant correlation was observed between NRB and bioactivity of the NAA (r = 0.03). NRB and PSA have been reported as determinants of oral absorption in a rat model  and specifically, compounds which meet the criterion of ten or fewer NRB are predicted to have good oral bioavailability . Glaxo SmithKline (GSK) has also devised the rule of seven rotatable bonds: “Less than seven rotatable bonds are essential for good bioavailablity” . Based on Veber et al. criteria , the results revealed that 92 % of HA had less than ten NRB, which is similar to 96.2 % observed in CRAD. Using the GSK criterion, a higher percentage of HA (81.27 %) had less than seven NRB in comparison to CRAD (70.3 %). This result attests to the potential of NAA to have good oral absorption.
Total polar surface area
Total polar surface area (TPSA) has also been implicated as a predictor for gastro-intestinal tract (GIT) penetration by many investigators [47, 48], and may be a key factor to consider during prioritization of NAA for lead development. The TPSA of CRAD was significantly lower (p < 0.05) than the TPSA of sub-groups of NAA (see Additional file 3). Among the latter, the highest TPSA was observed in MA followed by the HA, A and N groups. There was no significant correlation between TPSA and the activity of the NAA (r = 0.13). According to Veber et al. , the upper limit for TPSA for a molecule to penetrate the GIT is around 140 Å2.
All of the CRAD, up to 75 % of the HA and over 75 % of A compounds fell below the upper limit for the TPSA. This implies that most of the active NAA are expected to have good oral absorption, which is particularly relevant during the next stage of in vivo anti-malarial assessment. The high number of polar groups (hydroxyl and carbonyl groups) may be responsible for the high TPSA observed in some of the NAA. In addition, TPSA may also determine the extent of plasma protein binding . Albumin is involved in binding of mainly polar compounds , and high TPSA resulting from ionization of polar groups (e.g., acidic groups) may increase plasma protein binding. The higher a compound is bound to plasma protein the lower the proportion of the compound free for therapeutic effect. It is likely that some NAA that have good activity in in vitro assays may show poor activity in in vivo assays due to high plasma protein binding resulting from high TPSA. Optimization strategies that are geared towards reduction of TPSA (e.g., methylation of hydroxyl groups)may improve the oral absorption, reduce plasma protein binding and consequently improve bioactivity in vivo.
Van der Waals hydrophobic surface areas of hydrophobic atoms (vsa_hyd)
The Van der Waals hydrophobic surface area of hydrophobic atoms (vsa_hyd) measures the level of hydrophobicity of compounds. The vsa_hyd of HA and MA of NAA were significantly higher (p < 0.05) than that of CRAD (see Additional file 3). The vsa_hyd of N was lower than the vsa_hyd of CRAD. The vsa_hyd of the active NAA (A, HA, MA) were significantly (p < 0.05) higher than the low active NAA (N) suggesting that hydrophobicity may be vital for bioactivity. This observation was corroborated by the slight negative correlation (r = −0.231) seen between vsa_hyd and IC50 of NAA.
Hydrophobicity determines many biological processes, such as transport, distribution, metabolism, and molecular interactions of biological molecules. It is reported that the binding affinity and drug efficacy can be optimized and increased by incorporating hydrophobic groups [51, 52]. It is estimated that addition of a methyl group will lead to a 3.5-fold increase in binding constant . Therefore, moderately active NAA (MA) may be optimized as described above to improve bioactivity. Notably, improved hydrophobic interactions may also increase incidence of side effects and toxicity . With regard to absorption of bioactive compounds via biological membrane, hydrophobicity is also a key factor in various absorption models, e.g., Lipinski’s rule of five . Poor absorption or permeation is more likely for compounds with low hydrophobicity and the NAA that possess a higher vsa_hyd (hydrophobicity) than CRAD may be expected to show good absorption and permeation.
Overall, this brings to fore the pivotal role of hydrophobicity in achieving delicate balance of desirable activity, low toxicity and good absorption. Optimization of hydrophobicity of NAA towards a reference point, as observed for CRAD and as reported for marketed drugs, may be desirable to ensure successful development of these compounds.
Number of chiral centres
Chiral centres (asymmetric) are tetrahedral atoms (usually carbons) that have four different substituents . Compounds that have chiral centres are optically active and rotate the plane of polarized light to the left (levorotatory) or to the right (dextrorotatory) . Such optically active pairs are referred to as enantiomers .
The average number of chiral centres for HA and MA were significantly higher (p < 0.05) than CRAD (see Additional file 3). Among the NAA, the average number of chiral centres showed slight negative correlation (r = −0.2) with antiplasmodial activity (IC50); HA had the highest average number of chiral centres and N had the least. This result suggests that a high number of chiral centres may be essential for antiplasmodial activity. This may be because a high number of chiral centres increase flexibility of compounds and their tendency for more interaction with binding site of macromolecules. Moreover, the higher number of chiral centres observed in the NAA compared to CRAD, which have been previously observed between natural products and synthetic compounds , suggests that more compounds in HA, A and MA may have enantiomers. Given that enantiomeric molecules may interact in a different mode with biological receptors, binding affinities can differ between enantiomers . In clinical settings, enantiomers of chiral drugs can have decreased, had no, or even adverse effects [57–59]. Therefore, it is imperative to elucidate which of the enantiomers may be responsible for the observed antiplasmodial activities of the NAA and separate such enantiomers. Although the technology for the separation and analysis of chiral compounds has greatly advanced [59, 60], this may not yet be available in many laboratories in least developed countries where malaria is endemic and where these NAA may be sourced. In light of this limitation it may be necessary to prioritize NAAs without chiral centres (15 % of HA have no chiral centre) for preclinical development. Another option is to create racemic mixture (containing two enantiomers) provided the safety and efficacy can be justified as required by the Food Drug Administration of the United States of America . In conclusion, most of the compounds in NAA had higher number of chiral centres than those in CRAD and the number of chiral centres seems to correlate with antiplasmodial activities of NAA.
Number of oxygen atoms
The NAA had significantly (p value < 0.05) greater number of oxygen atoms than the CRAD (see Additional file 3). This may be because the NAA (consisting of alkaloids, terpenoids and flavonoids) are rich in oxygen atom as earlier reported for natural products [31–33, 62]. There was no significant difference (p > 0.05) in the number of oxygen (nO) among the sub-groups of NAA and little correlation of nO with bioactivity (r = −0.09). Although oxygen atoms, particularly sp2 -oxygen atom is important for ligands to form hydrogen bonds with receptors/enzymes , the low correlation with bioactivity [IC50] observed in the results suggests that the number of oxygen atoms may not be the sole contributing factor to the reported antiplasmodial activities of the NAA.
Number of nitrogen atoms and amine functional groups
The average number of nitrogen atoms was significantly higher (p value < 0.05) in CRAD than NAA (see Additional file 3), in accord with previous observation for natural products [31–33], although the most prevalent chemical class in the NAA was alkaloids (which usually contain nitrogen atoms). However, the other prevalent chemical classes, terpenoids and flavonoids, have no or low number of nitrogen atoms, which may explain the lower average number of nitrogen atoms observed in NAA compared to CRAD. Among the NAA, a slight negative correlation (r = −0.17) was observed between the number of nitrogen atoms and bioactivity. This suggests that the nitrogen atom may be relevant for anti-malarial activities .
Most of the nitrogen atoms were present as amine functional groups. The amine groups were predominant in CRAD and some of the compounds in HA sub-group of NAA. The presence of amine groups (or basic nitrogen) in nearly all CRAD highlights its importance for anti-malarial activity [64, 65]. However, the absence of the amine group in some of the HA and A compounds suggest that other functional groups, present within NAA compounds, are also vital for good antiplasmodial activities. These non-amine functional groups may have the potential to be the basis of a new series of anti-malarial compounds. A typical example is artemisinin-based drugs, which are endoperoxides.
Number of aromatic rings
Limited numbers of aromatic rings have been shown to improve developability of hits to drug candidates . The median number of aromatic rings for CRAD and all the sub-groups of NAA were between one and two. The average number of aromatic rings for CRAD was similar to the majority of the NAA (see Additional file 3). There was no correlation between the number of aromatic rings and bioactivity of NAA (r = −0.03).
Increase in aromatic ring count correlates with decreased aqueous solubility (with attendant poor absorption) [41, 62, 66, 67], increased plasma protein binding (leading to low clearance and therapeutic efficacy) , increased potential for inhibition of enzyme [with attendant toxicity and drug interactions] , increased mean observed hERG activity [hERG toxicity]  and decreased chance of developability of compounds into marketed drugs [66, 69, 70]. On these premise, NAA with minimal aromatic rings are expected to possess good aqueous solubility that is essential for oral absorption and exhibit reduced plasma protein binding that may lead to a greater free fraction of such compounds in systemic circulation, particularly during in vivo assay. In addition, low potential for enzyme inhibition and minimal toxicity [including hERG toxicity] as well as potential to be developed into successful anti-malarial drug candidates may also be expected from these NAA. It is however worth noting that a high number of aromatic rings may be desirable. This is because increasing number of aromatic rings decreases entropy [random movement] of molecules and favours binding of compounds to biological targets [66, 71, 72]. Increasing the number of aromatic rings during lead optimization may be favoured as a means to increase potency. Another reason why aromatic rings are favoured in drug design programmes is because of the well-established synthetic methods to make aryl–aryl links , which has made it attractive to design and synthesize compounds with increased number of aromatic rings in most combinatorial libraries .
Overall, limiting the number of aromatic rings during optimization of these NAA will make them broadly more developable and more ‘drug-like’ despite the likelihood of aromatic rings to increase potency and be readily amenable to synthesis and transformation .
Shape index, computed with Datawarrior, is a parameter that estimates the 3D shapes of compounds. Shape index less than 0.5 suggests presence of 3D (non-flat or spherical) scaffolds while shape index greater than 0.5 is for flat scaffolds. The average shape index of CRAD (0.53) was significantly (p < 0.05) higher than that of HA (0.43), A (0.45) and MA (0.46) but not significantly (p = 0.34) different from N (0.52) (see Additional file 3). This result suggests that CRAD and N contain compounds with flat scaffolds while HA, A and MA contain compounds with spherical or non-flat scaffolds. Other authors [74, 75] have reported the predominance of non-flat compounds in natural products as observed for NAA. Among the sub-groups of NAA, only the low active (N) NAA showed mean shape index >0.5, indicating the prevalence of compounds with flat scaffolds. There was some positive correlation between the shape index and in vitro antiplasmodial activity (r = 0.319). This suggests that non-flat scaffolds may be essential for antiplasmodial activity.
Generally, the results showed significant difference (p < 0.05) between the shape indices of NAA and CRAD. The presence of non-flat scaffolds seems to be essential for antiplasmodial activity amongst the NAA.
Synthetic feasibility [rsynth]
This parameter estimates how feasible it is to synthesize the compounds, with 1 being the most synthetically feasible and 0 the least synthetically feasible. The results (see Additional file 3) showed that all the sub-groups of NAA showed mean synthetic feasibility (rsynth) values that were significantly lower (p < 0.05) than that of CRAD (0.659). Among the NAA, there was very little positive correlation (r = 0.12) between rsynth and bioactivity (IC50). These results suggest that the CRAD and NAA N may be relatively easier to synthesize. The prevalence of flat and low MW (low complexity) compounds in N may be responsible for this observation. Synthetic feasibility and cost of synthesis may have significant impact on the development and eventual cost of drugs. This is particular relevant for neglected diseases in low-to-middle income countries where low-cost drugs are desired. Natural products have a high number of chiral centres that require advanced chemical synthetic techniques and chiral separation technology. This is evident in artemisinin, a recent drug of choice for malaria, which is still being sourced from the plant, A. annua, because its chemical synthesis, although achieved in 1983, is too expensive for commercialization .
Ligand efficiency metrics
The calculation of the binding efficiency metrics for HA, A, MA, and N sub-groups of NAA are described as well as the contribution to, and impact of these metrics, on prioritization and selection of NAA to take forward into anti-malarial drug development.
Ligand efficiency assesses the contribution of heavy atoms in or MW of a compound to potency or binding affinity of such a compound i.e., potency or binding affinity per heavy atom/molecular weight, given by Eq. (1) .
where ΔG = −RT ln (IC50/2), R = gas constant and T = absolute temperature. The unit of LE is kcal/mol/non hydrogen atom or heavy atom.
Figure 3, panel a showed that the mean LE for HA (0.48) was significantly higher than that of A (0.44) and MA (0.42) but not N (0.45). The exemplar values for LE should be greater than 0.3 kcal per mole per heavy atom  and the results (Fig. 3, panel a) showed that the mean LE values for most of the NAA (A, HA, MA, N) were within these exemplar values. This indicates that a good proportion of NAA (up to 80 %) have desired LE (potency at the right weight). It is particularly important to identify compounds with low weight and low potency (may be present within the MA sub-group of NAA) because it has been reported that such compounds have ‘room’ for optimization to increase potency and pharmacokinetic properties without the risk of losing LE [77, 79, 80]. A downside of LE is that it does not take lipophilicity, which is an important determinant of binding and/or potency, into account in its estimation of efficiency of binding or potency . Ligand lipophilicity efficiency (LLE) however provides a link between potency/binding affinity and lipophilicity.
Ligand lipophilicity efficiency
LLE measures how efficiently a ligand/compound exploits its lipophilicity to bind to target protein or create its potency [Eq. (2)]. In other words, it evaluates how well compounds improve potency while maintaining low lipophilicity .
The results (Fig. 3, panel b) showed that HA displayed a significantly higher average LLE value (HA, 6.05) than A (4.63), MA (4.14), and N (4.12). The ideal value for LLE has been reported to be greater than five  and only HA had a mean LLE value above five in spite of the similar clogP value with A and MA sub-groups of NAA. This suggests that the hydrophobic region of compounds in HA may be in such orientation that ensures optimal interaction with biological targets that brought about the observed bioactivity. Since the ultimate goal is to have compounds with good potency at the minimal lipophilicity, the results therefore suggests that HA consists of compounds that may be good starting points for anti-malarial drug development. Other studies have used LLE as criteria to find compounds suitable as starting points for optimization and drug development . On the other hand, compounds with low potency and low lipophilicity (that may be found in the MA sub-group of NAA) have also been reported as good starting points for drug development [83, 84]. This is because such compounds have big ‘lipophilicity room’ that are generally ‘filled’ during optimization towards improved potency. Monitoring the LLE of a compound collection during optimization will also allow medicinal chemists to track the efficiency of each lipophilic addition made towards improved potency.
One limitation of LLE is that it does not account for molecular size [heavy atom or MW]. A binding efficiency metric that combine lipophilicity, molecular size and potency is the ligand efficiency dependent lipophilicity index [LELP] .
Ligand efficiency dependent lipophilicity index
LELP is calculated using Eq. (3). LELP has been shown to reliably identify fragments, lead-like and drug-like compounds [77, 85]. Moreover, LELP was a better predictor of pharmacokinetic liabilities than LLE . The ideal LELP values have been stated to be between -10 and 10 for acceptable leads .
Looking at the results (Fig. 3 panel c), average LELP obtained for HA (8.93) was significantly lower than that of A (11.46) and MA (11.81) but somewhat similar to that of N (9.09). In addition, only HA and N showed mean LELP value within the ideal range. In so much as it has been reported that compounds with LELP values outside the exemplar range may not proceed far in the drug development pipeline , it is anticipated that NAA that fall outside the ideal LELP range may have lower chance of success in the anti-malarial drug development process. Moreover, lead optimization strategies should aim to increase LE or reduce logP in order to bring elevated LELP values within the desired range. In addition, monitoring LELP will help to control essential physicochemical properties that will maintain desirable potency and pharmacokinetic profile during optimization [80, 84].
In conclusion, ranking and selection of NAA from the initial list of potential antiplasmodial hits is a critical step in successful anti-malarial drug discovery [87, 88]. Given the influence of logP and molecular size [heavy atoms or MW] on potency and pharmacokinetic properties, the use of binding efficiency indices [LE, LLE, LELP] as a guiding criteria is important not only for hit selection, but also for lead generation and optimization [80, 84]. A plot of LELP against LLE (Fig. 4), previously described by Tarcsay et al. , may give medicinal chemists an idea of where the NAA compounds are in terms of these parameters and guide the optimization process to get the compounds to the desired region (as shown in Fig. 4). A key consideration is to be aware of the optimizing strategies that can increase potency and keep LE more or less constant or within exemplar limits.
Molecular similarity between CRAD and NAA
Molecular similarity analysis was conducted on NAA and CRAD. The objectives of this analysis were: firstly, to examine the extent of molecular diversity within the NAA and secondly to identify NAA that are structurally similar to and diverse from CRAD. Highly active NAA that are structurally diverse from CRAD may be potentially new anti-malarial agents with novel mechanism of actions.
The result of the molecular similarity assessment based on most common substructure, using the MoSSMCSS algorithm  is presented as a heat map (Fig. 5). The heat map showed a larger area of low similarity (lower values of Tanimoto coefficient) amongst the compounds. This is indicative of the substructural diversity amongst the NAA. A closer look at the heat map, using similarity viewer in KNIME, revealed Tanimoto coefficient in the range 0.1–0.7 between NAA and CRAD. This suggests that most of the NAA are structurally diverse from CRAD.
A similar observation was made when ChemGPS-NP  was used to compare the spatial coordinates and volume of chemical space occupied by NAA relative to CRAD. The result (Fig. 6) shows a plot of the first four dimensions: aromaticity, lipophilicity and flexibility, representing PC2, PC3 and PC4, respectively, plotted on the x, y and z axes. PC1, which represented size, was indicated by the size of the markers on the plot. The compounds in CRAD were identified by coloured markers. A cursory look at Fig. 6 shows that the CRAD did not form a tight cluster but were dispersed within the chemical space. A wider dispersion was however observed for NAA particularly in the positive direction of PC1 (size) A closer look revealed that NAA were bigger in size (PC1) and more aromatic (PC2) than the CRAD. Additionally, the dispersal of NAA and CRAD appears to be similar along the PC4 (flexibility) but slightly different along the PC3 (lipophilicity) with the NAA tending towards less lipophilicity. The proximity of some of the NAA to some of the CRAD may be a sign of their desirable drug-like properties and their amenability as a starting point in anti-malarial drug development. Although other studies have revealed that natural products occupy, in comparison to drug-like compounds, unique regions of property space [32, 89, 90], this result showed that some of the NAA occupy similar chemical spaces as CRAD, which also consist of natural product based anti-malarial drugs.
Based on the ‘similarity property principle’ , NAA that are diverse from CRAD may have different bioactive profiles or mechanism of action. Structurally diverse NAA may have different spatial orientation within biological targets in Plasmodium with resultant unique molecular interactions . This is particularly relevant in light of recent resistance of Plasmodium to CRAD. Compounds that may have different interactions pattern with targets implicated in Plasmodium resistance to CRAD or that interact with new targets are highly desirable for anti-malarial drug development. Moreover, such structurally diverse active antiplasmodial compounds provide insight into new chemical groups required for anti-malarial activities.
Conversely, NAA that are structurally similar to CRAD may have similar pharmacokinetic properties and may be drug-like. Since pharmacokinetic properties have been implicated as a major determinant of compound success or attrition during drug development [93, 94], NAA that are structurally similar to CRAD may suffer less attrition going through anti-malarial drug development. In addition, the chemical space analysis (Fig. 6), using the CRAD as reference or signpost, allows the visualization of relative position of NAA in chemical space compared to CRAD. Such view may enable the medicinal chemist to identify NAA that are within or outside the desired region (i.e., space occupied by CRAD). It may also assist the medicinal chemist to recognize the necessary properties to optimize and the extent of optimization required to move the NAA towards the desirable drug-like region or ‘sweet spot’ .
Overall, molecular diversity from CRAD may be indicative of new mechanism of action and potential for circumvention of current drug resistance while molecular similarity to CRAD may be indicative of favourable drug-like profile.
Structure–activity landscape: identifying activity cliffs
Exploration of structure–activity landscape represents a core aspect of medicinal chemistry . Activity cliff has been defined as pair of structurally similar compounds with large difference in bioactivity/potency  and has been of interest to the medicinal and computational chemist for a long time [97, 98]. To identify a pair of NAA that display activity cliffs, structure–activity similarity analysis was conducted with DataWarrior .
The result is presented as a SOM (Fig. 7) that display the relative position, in a 2D space, of all the NAA. Similar compounds are connected with a line and the markers are coloured by antiplasmodial activity (IC50) of the compounds from green (active NAA (≤5 uM)) to dark blue (inactive NAA (≥45 uM)). Clusters of green markers connected with lines (one of the cluster is marked ‘A’) representing similar NAA with similar antiplasmodial activity were observed within the landscape. This group of NAA form the ‘smooth region’ of the structure–activity landscape where minor changes in molecular structure usually lead to small change in bioactivity. This collection of NAA may be particularly appealing because it will allow the medicinal chemist to rationalize chemical substitutions that will improve pharmacokinetic parameters without sacrificing potency or bioactivity. This group of compounds are also amenable to quantitative structural activity relationship (QSAR) modelling because their structure–activity property aligns with the assumption of statistical modelling .
A closer look at Fig. 7 also reveals few clusters, one is marked ‘B’, which contains green markers (active NAA) connected to red and blue markers (inactive NAA). These clusters of NAA are structurally similar but have diverse bioactivity representing typical activity cliffs (see Additional file 4). Two examples of pairs of NAA that displayed activity cliff are shown in Fig. 8.
A plot of molecular similarity against activity similarity between pairs of compounds with markers coloured by fold change in normalized activity value [Δp(IC50)], provides another visualization for identification of the pairs of NAA that exhibit activity cliff (see Additional file 5). Markers that fall within activity similarity <0.8 (low activity similarity) and molecular similarity >0.8 (high molecular similarity) as well as high Δp(IC50) value signify pairs of NAA that show an activity cliff. Conversely, markers with activity similarity >0.8 (high activity similarity) and molecular similarity >0.8 (high molecular similarity) as well as low Δp(IC50) value signify pairs of NAA that are in the smooth region of the structure–activity landscape. The smooth region consists of compounds that adhere to the similarity principle and are pliable and desirable for QSAR models, as mentioned earlier.
Although Maggiora  proposed that activity cliffs may be responsible for the inefficient performance of many quantitative QSAR models , activity cliffs help pinpoint regions of the activity landscape that contain maximum information for structure activity relationship [SAR] studies. This is because this allows the medicinal chemist to identify the subtle molecular difference between a pair of compounds responsible for a dramatic shift in bioactivity. The rich SAR information from activity cliffs have been used in many drug discovery studies [99–101].
Pharmacokinetic profiling and drug-likeness of NAA
The pharmacokinetic profile and drug-likeness of NAA were assessed in silico using the models discussed below:
Lipinski’s rule of five
Lipinski’s rule of five got its name from the cut-off values for each of the four parameters that define the potential of a drug candidate for good absorption: the molecule has less than five HBD and less than ten HBA, its MW is below 500, and its LogP is less than five .
These four parameters were calculated for NAA as well as CRAD and plotted on a 3D graph with markers sized by MW (see Additional file 6). It was observed that most of the CRAD occupied space that was within the Lipinski’s rule of five, while most of the NAA (HA, A, MA) were dispersed away from Lipinski’s rule of five space. This suggests that the NAA (HA (47 %), A (52 %), MA (48 %), and N (71 %)) within Lipinski’s rule of five space may likely have good passive absorption. It was noted that a greater proportion of the N NAA showed propensity for good absorption. This suggests that though violation of Lipinski’s rule of five confers absorption liabilities it may be associated with good bioactivity among the NAA. Although limitations of Lipinski’s rule of five to predict the absorption of natural products and molecules that are actively transported have been reported [102–105], it is desirable that NAA that fall within Lipinski’s rule of five space be prioritized to lower attrition rates during anti-malarial drug development and increase the chance of new anti-malarial drugs reaching the market [26, 81].
Although formulation and drug delivery strategies have been developed to improve the absorption of compounds that violate Lipinski’s rule of five (i.e., poorly absorbed drugs) [106–110], it is important to be aware of the cost of such technology and its impact on the eventual market price of an anti-malarial drug. This is in view of the economically disadvantaged population in malaria-endemic regions of the world that require these anti-malarial drugs.
Veber et al. model
Veber et al.  suggested NRB and PSA of a compound as determinants of oral absorption. A plot of NRB and PSA was generated for NAA and CRAD (see Additional file 7, panel A). Veber et al. found that the majority of compounds with good oral bioavailability in rats had fewer than 10 rotatable bonds and PSA less than 140 Å2.
All the CRAD, except one with 11 NRB (halofantrine), were within the desirable area (red rectangle). A large number of the NAA [A (81 %), HA (76 %)and MA (80 %)] were also present within this region suggesting that these compounds may be well absorbed orally. The NAA that were dispersed outside the desired area and along the NRB axis have long aliphatic chains in their structure, while those dispersed along the PSA axis contain high number of hydroxyl and carbonyl groups. Lead optimization strategies may change single aliphatic bonds to double bonds to reduce rotation and polar groups may be methylated to reduce PSA to improve oral absorption of such compounds. However, the effect of such modification on the bioactivity of the compounds needs to be monitored.
Egan Egg plot [passive gut absorption]
Similar to Lipinski’s rule of five, Egan et al.  used statistical analysis to correlate passive intestinal absorption with PSA and clogP. A plot of PSA against clogP for NAA and CRAD is shown in Additional file 7, panel B. The ellipsoidal area of the plot (aka Egan Egg) enclose compounds that are expected to have good passive gut absorption. Compounds that fall outside the outer Egan Egg are predicted to have poor passive gut absorption, but may be absorbed by active transport processes. The results showed that most of the NAA, like the CRAD, fall within and at the border of the Egan Egg, suggesting that they may be well passively absorbed. A closer look revealed that 50 % of compounds in HA may have good oral absorption, 15 % may have borderline oral absorption while 45 % may show poor oral absorption. Slightly similar distribution in proportion was observed for compounds in A and MA. This model not only identifies NAA, especially HA and A, that may have poor absorption but the implicated physicochemical properties (PSA or clogP) may be identified and noted as one of the parameters to be addressed during lead optimization.
Golden triangle model
An analysis of Caco-2 permeability data for more than 16,000 compounds and human liver microsome clearance data for about 47,000 on a plot of distribution constant (logD) versus MW showed that compounds with good permeability and low clearance are concentrated within a triangular-shaped area (golden triangle) . A similar plot, logD (at pH 7.4) versus MW, was generated for the compounds in the dataset (see Additional file 7, panel C). The results showed a small proportion of the NAA (as well as CRAD) within the golden triangle (defined by: base of triangle is logD between −2 and 5 and apex is at MW 500). Using DataWarrior to select compounds within the triangle, it was observed that the following proportions of NAA were predicted to have tendency for good permeability and low clearance: A (35 %), HA (25 %) and MA (33 %). Approximately 33 % of the compounds in CRAD were predicted by this model to possess propensity for good permeability and low clearance. The bioavailability and clearance [half-life] data provided in the drug bank  for the CRAD predicted to possess good permeability and low clearance were explored to ascertain the consistency of the prediction. The data provided by drug bank for the identified compounds in CRAD align with some of the predictions by this model.
Fragment-based drug-likeness of NAA and CRAD, calculated with DataWarrior , is presented as boxplots (Fig. 9, panel a). A higher proportion of the compounds in CRAD and HA, and a lower proportion of A, MA and N were in the positive region of the drug-likeness score (drug-likeness score >0). Comparison of the mean drug-likeness score of NAA to CRAD showed that the HA (mean = 1.389) was slightly lower while A, MA and N were significantly lower than CRAD.
This approach to assessment of drug-likeness is based on presence of substructure fragments which are frequently present in commercial drugs. A positive value shows that the query molecule contains predominantly fragments that are frequently present in commercial drugs. HA, among the NAA, contains a higher proportion of compounds that may be drug-like and may have greater chance of success during anti-malarial drug development.
Quantitative estimate of drug-likeness
Bickerton et al.  recommended a new metric, the QED, to estimate drug-likeness of hits, leads or drug candidates. The QED value (‘desirability functions’) is mapped onto a scale between 0 and 1, where a desirability of 1 signifies an ideal value of the drug-like property and a desirability of 0 relates to a completely intolerable outcome.
The QED calculated for CRAD and NAA (sub-groups), using ICM Chemist Pro (Molsoft), is summarized as a boxplot (Fig. 9, panel b). As expected, CRAD showed a distribution that tends towards 1 (ideal value of the drug-like property) with an average QED value of 0.602. The average QED value for HA, A and MA were significantly lower than that of CRAD. On the other hand, N showed a similar drug-like score to CRAD. This is a contrast to what was observed using the fragment-based drug-likeness score which shows that HA had similar drug-likeness score to CRAD.
This observation may be due to the different approach used to estimate the QED value. While the fragment-based drug-likeness is based on presence of substructure fragments in commercial drugs, QED is generated from eight properties commonly used to define drug-likeness: MW, logP, HBA, HBD, PSA, number of aromatic groups [AROM], NRB and ALERTS [the number of matches to undesirable functionalities] . The distributions of these eight properties for a set of oral drugs were conducted and a desirability of 1 was assigned to the property values of oral drugs that occur most commonly, and 0 to property values that are not observed.
The low QED value estimated for HA and A thus signify that these compounds have low similarity to bulk of oral drugs and may have reduced chance of success during drug development. On the other hand, the low similarity of HA and A to CRAD and bulk of oral drugs may also be an indication of their structural peculiarity or novelty (as observed during molecular similarity analysis). These compounds may be the starting point of new anti-malarial drugs with unique mechanism of action.
Toxicity potential assessment
The potential for toxicity from the NAA was assessed by checking for the presence of reactive chemical groups and potential to cause tumour, irritation and mutagenesis. Assessment of tumourigenic, irritant and mutagenic risk seeks to identify compounds with possibility to cause tumour, irritation and mutation in vivo. Figure 10, a visualization of the results of the toxicity assessment, showed that a greater proportion of NAA (>80 %) showed no risk for tumourigenic, irritant and mutagenic potential in comparison to CRAD (<60 %). This may attest to the low toxicity of these NAA, as previously observed for other natural products , and their potential as a source of new and safe anti-malarial drug candidates. However, it is worth noting that drug approvals are based on rigorous benefit-risk assessment [114, 115] and NAA, with high risk of the assessed toxicity parameters, may be considered for anti-malarial drug development provided the benefits from such compounds outweigh the potential risk.
In the context of drug development, reactive groups are usually electrophiles or free radicals that may react readily [via covalent binding] with nucleophilic components, such as DNA and proteins within the biological system [116, 117]. Reactive groups were observed in 33.3 % of CRAD and that was significantly lower than observed for the NAA (47.49 % in HA, 47.37 % in A, 55.4 % in MA, and 44.58 % in N). Reactive groups identified include peroxo, oxiran/aziridine, allyl/benzyl chloride, 2-halo-enone, 3-halo-enone, and quaternary ammonium. Reactive groups present within some drugs or formed after metabolism of such drugs within the biological system, have been implicated for unexpected toxicities of drugs that become apparent only after the launch of such drug entities [117–119].
Overall, these results showed that NAA have lower potential for toxicity in comparison to CRAD. Although the presence of reactive groups, tumourigenic, irritant and mutagenic risk is indicative of toxicity risk, these toxicity risk alerts are by no means to be taken as fully reliable predictors of toxicity. Nor should the absence of these toxicity risk alerts be a confirmatory indication that a compound will be completely free of toxicity. Nonetheless, in silico toxicity assessment of NAA early during hit profiling allows the de-prioritization of compounds that may have unexpected toxicity issues. In addition, implicated chemical groups may be replaced with other groups while retaining the biological activity of the compound through scaffold hopping. This may help to design and bring safer anti-malarial drugs to the market.
Assessing promiscuity of NAA: Eli Lilly MedChem rules
Promiscuous compounds or frequent hitters in NAA may be false positives from antiplasmodial assays . Identifying and flagging such compounds will guide selection of NAA hits, preferably excluding such frequent hitters, for the next stage of anti-malarial drug development.
Applying the Eli Lilly MedChem rules on the NAA and CRAD identified promiscuous or reactive compounds in both datasets. Approximately 63 % of the compounds in CRAD failed the rules with the predominate reasons for failure being the presence of ‘peroxide’ and ‘para quinone’ groups . The former captures direct oxidants, which may be the artemisinin derivatives present in the CRAD, while the latter identify para-positioned quinones that have high redox potential. Although these compounds in CRAD are registered drugs, they represent chemotypes that would have been triaged out of screening sets due to their high reactivity.
Promiscuity indicates a compound with non-specific mode of action and may lead to increased incidence of toxicity . In the case of NAA, 82 % of the compounds in HA failed the rules, 79 % compounds failed in A, 78 % in MA, while 64 % failed in N. The predominant reason for failure in all the sub-groups of the NAA was ‘presence of catechol’ group which implies the presence of dihydroxybenzene groups . Such groups, which have high redox capacity and are highly promiscuous, are predominant in natural products from plants, e.g., flavonoids, anthraquinones, terpenoids compounds, etc. However, there was no significant (p > 0.05) difference in the proportion of compounds with reactive groups between CRAD and sub-groups of NAA (except MA) (see Additional file 3). No correlation (r = −0.031) was observed between the number of reactive compounds and the bioactivity (IC50) of the compounds (see Additional file 3). These observations either underscore the aggressiveness of the Eli Lilly MedChem rules or attest to the presence of natural products and natural product-derived compounds in CRAD. One may argue that the rule should be used with caution but it is worthwhile as it will flag compounds with bioactivity to be looked at more closely. This will guard against the use of false positives or non-specific NAA as potential starting points for anti-malarial drug development.
Prioritization of NAA: integration of chemo-informatic profiling data
Using the information provided in this study, especially within an interactive data-mining environment such as DataWarrior software , it is possible to provide an answer to the key question posed: Should NAA that are most likely to be successfully developed into anti-malarial drug candidates be prioritized? One way to provide an affirmative response to this question was to adjust the filters of the molecular descriptors or/and physicochemical properties within DataWarrior to the reported exemplar values for all the calculated molecular descriptors and physicochemical properties. The result of this operation will be a prioritized list of NAA that possess the desired molecular descriptors or/and physicochemical properties. About 28 % of the NAA fell within these exemplar values or limits. However, the proportion of NAA prioritized may vary based on the selected molecular descriptors and physicochemical properties and the choice of exemplar values. In addition, the use of filters will permanently remove seemingly interesting compounds that do not fall within the range set on the filters.
Another approach used to get a prioritized list of NAA was to generate a consensus scoring function for each compound in NAA. The cut-off or exemplar values for each calculated molecular descriptor and physicochemical property were used as the benchmark to score each compound (molecular descriptor and physicochemical property that were not significantly different between the two datasets, i.e., CRAD and NAA, were excluded). Compounds with the values of molecular descriptors and physicochemical properties outside the desirable range were penalized with a score of −1 while compounds within the desirable range were rewarded with a score of 1. The average of these scores for all the molecular descriptors and physicochemical properties discussed in this study were taken as the consensus score. A consensus score of 1 suggests that such compounds have all molecular descriptors and physicochemical properties within the acceptable range and may be prioritized for the next stage of pre-clinical anti-malarial drug development. The prioritized list of the NAA with their consensus score [list was sorted by the consensus score, highest to lowest] is shown in Additional file 8. As expected, the consensus scoring showed high values (0.5–1.0) for over 90 % of the CRAD, which have successfully passed through- anti-malarial drug development. Prospective NAA that showed high score (e.g., above 0.5) share similar drug-like properties with CRAD and may stand a greater chance of successfully passing to development to become anti-malarial drugs.
Overall, compounds within NAA with consensus score close to 1 (i.e., those that fall on the positive side of the cut-off values of the various chemo-informatic properties assessed) may have greater chance of success during anti-malarial pre-clinical drug development. Alternatively, visualization of all metrics in the context of CRAD can aid prioritization and selection of NAA for downstream anti-malarial drug discovery (Fig. 11). The graph allows visualization of the profile of NAA in comparison to CRAD based on selected molecular descriptors. It is expected that NAA that have similar profile to CRAD (i.e., have similar values of selected molecular descriptors) should be more drug-like and therefore be prioritized for the next stage of anti-malarial drug development.
Chemo-informatic profiling of NAA and CRAD has led to development of prioritization strategies and prioritized lists of at least 1000 compounds that may guide decisions and facilitate anti-malarial drug development from natural products with antiplasmodial activities. This prioritized list includes structurally diverse NAA that may encompass new biologically relevant chemical space and could be developed into anti-malarial drug candidates with possible different mechanism of action. Structure–activity landscape analysis revealed NAA pairs that form activity cliffs, which are particular relevant for SAR studies. Finally, this study was able to identify NAA with desired drug-like properties and toxicity liabilities as well as promiscuous compounds or ‘frequent hitters’ among the NAA.
natural products with in vitro antiplasmodial activities
highly active natural products with in vitro antiplasmodial activities
active natural products with in vitro antiplasmodial activities
moderately active natural products with in vitro antiplasmodial activities
low active natural products with in vitro antiplasmodial activities
currently registered anti-malarial drugs
hydrogen bond donor
hydrogen bond acceptor
total molecular weight
number of rotatable bonds
total polar surface area
Van der Waals hydrophobic surface area of hydrophobic atoms
ligand lipophilicity efficiency
ligand efficiency dependent lipophilicity
Konstanz information miner
quantitative structural activity relationship
structural activity relationship
quantitative estimate of drug-likeness
gastro instestinal tract
MoSS most common substructure
structure activity landscape index
WHO. World malaria report 2011. Geneva: World Health Organization; 2011.
WHO. Global Malaria Action Plan. World Health Organization, Geneva. Available online at: http://www.rbm.who.int/keyfacts.html.Downloaded January 2011.
Wells TN. Natural products as starting points for future anti-malarial therapies: going back to our roots. Malar J. 2011;10(Suppl 1):S3.
Phillipson JD. Phytochemistry and medicinal plants. Phytochemistry. 2001;56:237–43.
ElSohly HN, Croom E, ElSohly M. Analysis of the antimalarial sesquiterpene artemisinin in Artemisia annua by high-performance liquid chromatography (HPLC) with postcolumn derivatization and ultraviolet detection. Pharm Res. 1987;4:258–60.
Batista R, Silva Ade J Jr, de Oliveira AB. Plant-derived antimalarial agents: new leads and efficient phytomedicines. Part II. Non-alkaloidal natural products. Molecules. 2009;14:3037–72.
Frederich M, Tits M, Angenot L. Potential antimalarial activity of indole alkaloids. Trans R Soc Trop Med Hyg. 2008;102:11–9.
Kaur K, Jain M, Kaur T, Jain R. Antimalarials from nature. Bioorg Med Chem. 2009;17:3229–56.
Nogueira CR, Lopes LM. Antiplasmodial natural products. Molecules. 2011;16:2146–90.
Harvey A. Strategies for discovering drugs from previously unexplored natural products. Drug Discov Today. 2000;5:294–300.
Butler MS, Fontaine F, Cooper MA. Natural product libraries: assembly, maintenance, and screening. Planta Med. 2014;80:1161–70.
Bains W. Failure rates in drug discovery and development-will we ever get any better? Drug Discovery World. 2004;5:9–19.
Adams CP, Brantner VV. Estimating the cost of new drug development: is it really 802 million dollars? Health Aff (Millwood). 2006;25:420–8.
Barrenho E, Miraldo M, Smith P. The determinants of attrition in drug development: a duration analysis. Discussion paper 2013/4, Health Management Group, Imperial College Business School, (https://spiral.imperial.ac.uk/bitstream/10044/1/12204/4/Barrenho%202013-04.pdf).
Herper M. The cost of creating a new drug now $5 billion, pushing big pharma to change. http://www.forbes.com/sites/matthewherper/2013/08/11/how-the-staggering-cost-of-inventing-new-drugs-is-shaping-the-future-of-medicine/. Accessed January 2015.
Mannhold R, Kubinyi H, Folkers G, Oprea TI. Chemoinformatics in drug discovery. John Wiley; 2006.
Jacoby E, Schuffenhauer A, Popov M, Azzaoui K, Vangrevelinghe E, Priestle J, et al. Molecular informatics as an enabling in silico technology platform for drug discovery. CHIMIA Int J Chem. 2004;58:577–84.
Hohman M, Gregory K, Chibale K, Smith PJ, Ekins S, Bunin B. Novel web-based tools combining chemistry informatics, biology and social networks for drug discovery. Drug Discov Today. 2009;14:261–70.
O’Boyle NM, Banck M, James CA, Morley C, Vandermeersch T, Hutchison GR. Open babel: an open chemical toolbox. J Cheminf. 2011;3:33.
Moe M. Chemical Computing Group. Quebec: Montreal; 2006.
Sander T, Freyss J, von Korff M, Rufener C. DataWarrior: an open-source program for chemistry aware data visualization and analysis. J Chem Inf Model. 2015;55:460–73.
Berthold MR, Cebron N, Dill F, Gabriel TR, Kötter T, Meinl T. KNIME-the Konstanz information miner: version 2.0 and beyond. AcM SIGKDD Explor Newsl. 2009;11:26–31.
Moss: a program for molecular substructure mining. Proceedings of the 1st international workshop on open source data mining: frequent pattern mining implementations: ACM; 2005.
Larsson J, Gottfries J, Muresan S, Backlund A. ChemGPS-NP: tuned for navigation in biologically relevant chemical space. J Nat Prod. 2007;70:789–94.
Bajorath J, Peltason L, Wawer M, Guha R, Lajiness MS, Van Drie JH. Navigating structure–activity landscapes. Drug Discov Today. 2009;14:698–705.
Lipinski CA. Lead-and drug-like compounds: the rule-of-five revolution. Drug Discov Today Technol. 2004;1:337–41.
Egan WJ, Merz KM, Baldwin JJ. Prediction of drug absorption using multivariate statistics. J Med Chem. 2000;43:3867–77.
Johnson TW, Dress KR, Edwards M. Using the golden triangle to optimize clearance and oral absorption. Bioorg Med Chem Lett. 2009;19:5560–4.
Veber DF, Johnson SR, Cheng H, Smith BR, Ward KW, Kopple KD. Molecular properties that influence the oral bioavailability of drug candidates. J Med Chem. 2002;45:2615–23.
Bruns RF, Watson IA. Rules for identifying potentially reactive or promiscuous compounds. J Med Chem. 2012;55:9763–72.
Singh SB, Culberson JC. Chemical space and the difference between natural products and synthetics. Natural product chemistry for drug discovery. Cambridge: RSC publishing; 2010.
Lachance H, Wetzel S, Kumar K, Waldmann H. Charting, navigating, and populating natural product chemical space for drug discovery. J Med Chem. 2012;55:5989–6001.
Vasilevich NI, Kombarov RV, Genis DV, Kirpichenok MA. Lessons from natural products chemistry can offer novel approaches for synthetic chemistry in drug discovery. J Med Chem. 2012;55:7003–9.
Grabowski K, Schneider G. Properties and architecture of drugs and natural products revisited. Curr Chem Biol. 2007;1:115–27.
Tan DS. Current progress in natural product-like libraries for discovery screening. Comb Chem High Throughput Screen. 2004;7:631–43.
Wermuth CG. The practice of medicinal chemistry. Amsterdam: Academic Press; 2011.
Gleeson MP, Hersey A, Montanari D, Overington J. Probing the links between in vitro potency, ADMET and physicochemical parameters. Nat Rev Drug Discov. 2011;10:197–208.
Lipinski CA. Poor aqueous solubility-an industry wide problem in ADME screening. Am Pharm Rev. 2002;5:82–5.
Leeson PD, St-Gallay SA. The influence of the ‘organizational factor’ on compound quality in drug discovery. Nat Rev Drug Discovery. 2011;10:749–65.
Proudfoot JR. The evolution of synthetic oral drug properties. Bioorg Med Chem Lett. 2005;15:1087–90.
Arnott JA, Planey SL. The influence of lipophilicity in drug discovery and design. Expert Opin Drug Discov. 2012;7:863–75.
Azzaoui K, Hamon J, Faller B, Whitebread S, Jacoby E, Bender A, et al. Modeling promiscuity based on in vitro safety pharmacology profiling data. ChemMedChem. 2007;2:874–80.
Bender A, Scheiber J, Glick M, Davies JW, Azzaoui K, Hamon J, et al. Analysis of pharmacology data and the prediction of adverse drug reactions and off-target effects from chemical structure. ChemMedChem. 2007;2:861–73.
Hughes JD, Blagg J, Price DA, Bailey S, DeCrescenzo GA, Devraj RV, et al. Physiochemical drug properties associated with in vivo toxicological outcomes. Bioorg Med Chem Lett. 2008;18:4872–5.
Greene N, Aleo MD, Louise-May S, Price DA, Will Y. Using an in vitro cytotoxicity assay to aid in compound selection for in vivo safety studies. Bioorg Med Chem Lett. 2010;20:5308–12.
Gleeson MP. Generation of a set of simple, interpretable ADMET rules of thumb. J Med Chem. 2008;51:817–34.
Leach AR, Hann MM, Burrows JN, Griffen EJ. Fragment screening: an introduction. Mol BioSyst. 2006;2:429–46.
Feng M. Assessment of blood-brain barrier penetration: in silico, in vitro and in vivo. Curr Drug Metab. 2002;3:647–57.
Ghafourian T, Amin Z. QSAR models for the prediction of plasma protein binding. Bioimpacts. 2013;3:21.
Colmenarejo G. In silico prediction of drug-binding strengths to human serum albumin. Med Res Rev. 2003;23:275–301.
DiMagno SG, Sun H. The strength of weak interactions: aromatic fluorine in drug design. Curr Top Med Chem. 2006;6:1473–82.
Qian SB, Waldron L, Choudhary N, Klevit RE, Chazin WJ, Patterson C. Engineering a ubiquitin ligase reveals conformational flexibility required for ubiquitin transfer. J Biol Chem. 2009;284:26797–802.
Bissantz C, Kuhn B, Stahl M. A medicinal chemist’s guide to molecular interactions. J Med Chem. 2010;53:5061–84.
Lipton SA. Turning down, but not off. Nature. 2004;428:473.
Sheldon RA. Chirotechnology: industrial synthesis of optically active compounds. New York: CRC Press; 1993.
Brooks WH, Guida WC, Daniel KG. The significance of chirality in drug design and development. Curr Top Med Chem. 2011;11:760–70.
McConathy J, Owens MJ. Stereochemistry in Drug Action. Prim Care Companion J Clin Psychiatry. 2003;5:70–3.
Smith SW. Chiral toxicology: it’s the same thing…only different. Toxicol Sci. 2009;110:4–30.
Sekhon BS. Enantioseparation of chiral drugs—an overview. Int J Pharm Technol Res. 2010;2:1584–94.
Rentsch KM. The importance of stereoselective determination of drugs in the clinical laboratory. J Biochem Biophys Methods. 2002;54:1–9.
De Camp WH. The FDA perspective on the development of stereoisomers. Chirality. 1989;1:2–6.
Feher M, Schmidt JM. Property distributions: differences between drugs, natural products, and molecules from combinatorial chemistry. J Chem Inf Comput Sci. 2003;43:218–27.
Kubinyi H. Hydrogen bonding: the last mystery in drug design. Pharmacokinetic Optimization in Drug Research. 2001:513–24.
de Villiers KA, Marques HM, Egan TJ. The crystal structure of halofantrine–ferriprotoporphyrin IX and the mechanism of action of arylmethanol antimalarials. J Inorg Biochem. 2008;102:1660–7.
Ramu K, Baker JK. Synthesis, characterization, and antimalarial activity of the glucuronides of the hydroxylated metabolites of arteether. J Med Chem. 1995;38:1911–21.
Ritchie TJ, Macdonald SJ. The impact of aromatic ring count on compound developability–are too many aromatic rings a liability in drug design? Drug Discov Today. 2009;14:1011–20.
Yusof I, Segall MD. Considering the impact of ‘drug-like’properties on the chance of success. Drug Discov Today. 2013;18:659–66.
Blagg J. Structure-activity relationships for in vitro and in vivo toxicity. Annu Rep Med Chem. 2006;41:353.
Hann MM, Leach AR, Harper G. Molecular complexity and its impact on the probability of finding leads for drug discovery. J Chem Inf Comput Sci. 2001;41:856–64.
Ward SE, Beswick P. What does the aromatic ring number mean for drug design? Expert Opin Drug Discov. 2014;9:995–1003.
Levitt M, Perutz MF. Aromatic rings act as hydrogen bond acceptors. J Mol Biol. 1988;201:751–4.
Perutz M. The role of aromatic rings as hydrogen-bond acceptors in molecular recognition. Philos Transact Royal Soc London Ser A Phys Eng Sci. 1993;345:105–12.
Buss AD, Butler MS. Natural product chemistry for drug discovery. Royal Society of Chemistry. 2010.
López-Vallejo F, Giulianotti MA, Houghten RA, Medina-Franco JL. Expanding the medicinally relevant chemical space with compound libraries. Drug Discov Today. 2012;17:718–26.
Lovering F, Bikker J, Humblet C. Escape from flatland: increasing saturation as an approach to improving clinical success. J Med Chem. 2009;52:6752–6.
Hommel M. The future of artemisinins: natural, synthetic or recombinant? J Biol. 2008;7:38.
Planey SL, Kumar R. Lipophilicity indices for drug development. J Appl Biopharm Pharmacokinet. 2013;1:31–6.
Hopkins AL, Groom CR, Alex A. Ligand efficiency: a useful metric for lead selection. Drug Discov Today. 2004;9:430–1.
Carr RA, Congreve M, Murray CW, Rees DC. Fragment-based lead discovery: leads by design. Drug Discov Today. 2005;10:987–92.
Mortenson PN, Murray CW. Assessing the lipophilicity of fragments and early hits. J Comput Aided Mol Des. 2011;25:663–7.
Leeson PD, Springthorpe B. The influence of drug-like concepts on decision-making in medicinal chemistry. Nat Rev Drug Discov. 2007;6:881–90.
Mowbray CE, Burt C, Corbau R, Gayton S, Hawes M, Perros M, et al. Pyrazole NNRTIs 4: selection of UK-453,061 (lersivirine) as a development candidate. Bioorg Med Chem Lett. 2009;19:5857–60.
Hann MM, Keserü GM. Finding the sweet spot: the role of nature and nurture in medicinal chemistry. Nat Rev Drug Discovery. 2012;11:355–65.
Keserü GM, Makara GM. The influence of lead discovery strategies on the properties of drug candidates. Nat Rev Drug Discovery. 2009;8:203–12.
Tarcsay A, Nyíri K, Keserű GM. Impact of lipophilic efficiency on compound quality. J Med Chem. 2012;55:1252–60.
Wager TT, Chandrasekaran RY, Hou X, Troutman MD, Verhoest PR, Villalobos A, et al. Defining desirable central nervous system drug space through the alignment of molecular properties, in vitro ADME, and safety attributes. ACS Chemical Neuroscience. 2010;1:420–34.
Di L, Kerns EH. Profiling drug-like properties in discovery research. Curr Opin Chem Biol. 2003;7:402–8.
Hughes J, Rees S, Kalindjian S, Philpott K. Principles of early drug discovery. Br J Pharmacol. 2011;162:1239–49.
Rosén J, Gottfries J, Muresan S, Backlund A, Oprea TI. Novel chemical space exploration via natural products. J Med Chem. 2009;52:1953–62.
Mander L, Liu H. Comprehensive natural products, vol 2 : Chemistry and Biology, Elsevier; 2010, ISBN: 978-0-08-045382-8.
Maldonado AG, Doucet J, Petitjean M, Fan B. Molecular similarity and diversity in chemoinformatics: from theory to applications. Mol Divers. 2006;10:39–79.
Bon RS, Waldmann H. Bioactivity-guided navigation of chemical space. Acc Chem Res. 2010;43:1103–14.
Kennedy T. Managing the drug discovery/development interface. Drug Discov Today. 1997;2:436–44.
Hou T, Xu X. ADME evaluation in drug discovery. Mol Model Annu. 2002;8:337–49.
Guha R. Exploring structure–activity data using the landscape paradigm. Wiley Interdiscip Rev Comput Mol Sci. 2012;2:829–41.
Hu Y, Stumpfe D, Bajorath J. Advancing the activity cliff concept. F1000Res. 2013;2.
Lajiness M. Evaluation of the performance of dissimilarity selection methodology. QSAR: Rational approaches to the design of bioactive compounds; 1991. p. 201–4.
Maggiora GM. On outliers and activity cliffs why QSAR often disappoints. J Chem Inf Model. 2006;46:1535.
Wawer M, Peltason L, Weskamp N, Teckentrup A, Bajorath J. Structure-activity relationship anatomy by network-like similarity graphs and local structure-activity relationship indices. J Med Chem. 2008;51:6075–84.
Vogt M, Huang Y, Bajorath J. From activity cliffs to activity ridges: informative data structures for SAR analysis. J Chem Inf Model. 2011;51:1848–56.
Stumpfe D, Bajorath J. Exploring activity cliffs in medicinal chemistry: miniperspective. J Med Chem. 2012;55:2932–42.
Keller TH, Pichota A, Yin Z. A practical view of ‘druggability’. Curr Opin Chem Biol. 2006;10:357–61.
Ganesan A. The impact of natural products upon modern drug discovery. Curr Opin Chem Biol. 2008;12:306–17.
Bauer RA, Wurst JM, Tan DS. Expanding the range of ‘druggable’targets with natural product-based libraries: an academic perspective. Curr Opin Chem Biol. 2010;14:308–14.
Doak BC, Over B, Giordanetto F, Kihlberg J. Oral druggable space beyond the rule of 5: insights from drugs and clinical candidates. Chem Biol. 2014;21:1115–42.
Rogge MC, Taft DR. Preclinical drug development. CRC Press; 2005.
Shaikh M, Derle ND, Bhamber R. Permeability enhancement techniques for poorly permeable drugs: a review. J Appl Pharm Sci. 2012;2:34–9.
Gupta S, Kesarla R, Omri A. Formulation strategies to improve the bioavailability of poorly absorbed drugs with special emphasis on self-emulsifying systems. ISRN pharm. 2013;2013.
Basavaraj S, Betageri GV. Can formulation and drug delivery reduce attrition during drug discovery and development—review of feasibility, benefits and challenges. Acta Pharmaceutica Sinica B. 2014;4:3–17.
Mahapatra AK, Murthy PN. Self-emulsifying drug delivery systems (SEDDS): an update from formulation development to therapeutic strategies. Self. 2014;6:546–68.
Wishart DS, Knox C, Guo AC, Shrivastava S, Hassanali M, Stothard P, et al. DrugBank: a comprehensive resource for in silico drug discovery and exploration. Nucleic Acids Res. 2006;34:D668–72.
Bickerton GR, Paolini GV, Besnard J, Muresan S, Hopkins AL. Quantifying the chemical beauty of drugs. Nat Chem. 2012;4:90–8.
Hoult J, Paya M. Pharmacological and biochemical actions of simple coumarins: natural products with therapeutic potential. Gen Pharmacol. 1996;27:713–22.
Garrison LP Jr, Towse A, Bresnahan BW. Assessing a structured, quantitative health outcomes approach to drug risk-benefit analysis. Health Aff (Millwood). 2007;26:684–95.
Eichler H, Pignatti F, Flamion B, Leufkens H, Breckenridge A. Balancing early market access to new drugs with the need for benefit/risk data: a mounting dilemma. Nat Rev Drug Discov. 2008;7:818–26.
Chen WG, Zhang C, Avery MJ, Fouda HG. Reactive metabolite screen for reducing candidate attrition in drug discovery. Biological Reactive Intermediates VI. US: Springer; 2001. p. 521–4.
Attia SM. Deleterious effects of reactive metabolites. Oxid Med Cell Longev. 2010;3:238–53.
Nelson SD. Mechanisms of the formation and disposition of reactive metabolites that can cause acute liver injury. Drug Metab Rev. 1995;27:147–77.
Prakash C, Sharma R, Gleave M, Nedderman A. In vitro screening techniques for reactive metabolites for minimizing bioactivation potential in drug discovery. Curr Drug Metab. 2008;9:952–64.
Baell JB, Holloway GA. New substructure filters for removal of pan assay interference compounds (PAINS) from screening libraries and for their exclusion in bioassays. J Med Chem. 2010;53:2719–40.
Pearce BC, Sofia MJ, Good AC, Drexler DM, Stock DA. An empirical process for the design of high-throughput screening deck filters. J Chem Inf Model. 2006;46:1060–8.
SAE, SM and AC conceived the project concept. ESA designed all experiments and carried out the analysis. JS, SM and AC contributed to data analysis. AC provided funding for the project. All authors read and approved the final manuscript.
This work was supported by the South African Research Chairs Initiative of the Department of Science and Technology and National Research Foundation of South Africa.
The authors declare that they have no competing interests.
About this article
Cite this article
Egieyeh, S.A., Syce, J., Malan, S.F. et al. Prioritization of anti-malarial hits from nature: chemo-informatic profiling of natural products with in vitro antiplasmodial activities and currently registered anti-malarial drugs. Malar J 15, 50 (2016). https://doi.org/10.1186/s12936-016-1087-y
- Natural products
- Anti-malarial drugs