Virtualizing the p-ANAPL Library: A Step towards Drug Discovery from African Medicinal Plants

Abstract

Background

Natural products play a key role in drug discovery programs, both serving as drugs and as templates for the synthesis of drugs, even though the quantities and availabilities of samples for screening are often limitted.

Experimental approach

A current collection of physical samples of > 500 compound derived from African medicinal plants aimed at screening for drug discovery has been made by donations from several researchers from across the continent to be directly available for drug discovery programs. A virtual library of 3D structures of compounds has been generated and Lipinski’s “Rule of Five” has been used to evaluate likely oral availability of the samples.

Results

A majority of the compound samples are made of flavonoids and about two thirds (2/3) are compliant to the “Rule of Five”. The pharmacological profiles of thirty six (36) selected compounds in the collection have been discussed.

Conclusions and implications

The p-ANAPL library is the largest physical collection of natural products derived from African medicinal plants directly available for screening purposes. The virtual library is also available and could be employed in virtual screening campaigns.

About the Authors:

Fidele Ntie-Kang

Contributed equally to this work with: Fidele Ntie-Kang, Pascal Amoa Onguéné

Affiliation: Department of Chemistry, Chemical and Bioactivity Information Centre, Faculty of Science, University of Buea, Buea, Cameroon

Pascal Amoa Onguéné

Contributed equally to this work with: Fidele Ntie-Kang, Pascal Amoa Onguéné

Affiliation: Department of Chemistry, Faculty of Science, University of Douala, Douala, Cameroon

Ghislain W. Fotso

Affiliation: Department of Organic Chemistry, University of Yaoundé I, Yaoundé, Cameroon

Kerstin Andrae-Marobela

* E-mail: [email protected] (LM); [email protected] (KM)

Affiliation: Department of Biological Sciences, University of Botswana, Gaborone, Botswana

Merhatibeb Bezabih

Affiliation: Department of Chemistry, University of Botswana, Gaborone, Botswana

Jean Claude Ndom

Affiliation: Department of Chemistry, Faculty of Science, University of Douala, Douala, Cameroon

Bonaventure T. Ngadjui

Affiliation: Department of Pharmaceutical Sciences and Traditional Pharmacopeia, Faculty of Medicine and Biomedical Sciences, University of Yaoundé I, Yaoundé, Cameroon

Abiodun O. Ogundaini

Affiliation: Department of Pharmaceutical Chemistry, Obafemi Awolowo University, Ile-Ife, Nigeria

Berhanu M. Abegaz

Affiliation: Academy of Sciences, Nairobi, Kenya

Luc Mbaze Meva’a

* E-mail: [email protected] (LM); [email protected] (KM)

Affiliation: Department of Chemistry, Faculty of Science, University of Douala, Douala, Cameroon

Introduction

The African continent is highly endowed with diverse vegetation types constituting tropical rainforests, coastal and alpine forests, savannahs, woodlands and scrublands which is a reservoir for diverse natural product (NP) classes [1], [2]. Natural products have always served as sources of inspiration for the design of new drugs and/or as drugs themselves [3]–[5]. In addition, it has been verified that natural products from Africa have an enormous potential for drugs [6]. The measured activities of compounds derived from plants of African origin have been reported in at least 2000 publications and it is believed that the active principles for the development of drugs against some of the neglected tropical diseases may be found in the African flora [5], [6], since some of these plants have been used extensively in the treatment of some of these diseases in African Traditional Medicine (ATM) [7]. In addition, the potential of natural products derived from African medicinal plants has been demonstrated by a number of recent review articles [6], [8]–[16]. In order to properly exploit these data for drug development, African researchers have initiated efforts towards the collection of compound samples which should be readily available for screening. It is important to note that although natural products play important/evolving roles in drug discovery [17], huge pharmaceutical companies have reduced emphasis on natural products in terms of new drug development for the last two decades [17]–[19]. There are many reasons for this decreasing interest is the strategy of drug discovery employed by pharmaceutical companies, of which the most important is the time factor involved in the search for NP lead compounds. The trend has been drastically changing from a laborious search for potent NP leads (bioassay-guided isolation of natural products from crude extracts, without guaranteeing reasonable quantities of isolates for testing) to the rapid high-throughput screening of molecular target-based pure compound chemical libraries. The chemical libraries are often generated, to a large extent, using combinatorial chemistry. However, streamlined screening procedures and enhanced organism sourcing mechanisms are among the new technologies which have been put in place in order to enhance natural product drug discovery in an industrial setting [20]. In order to considerably reduce the number of compound samples required for testing in bio-assays in drug discovery projects, computer-aided methods like docking and virtual screening (VS) have been employed. Computer-aided drug design (CADD) often involves VS of large compound datasets and the availability of such is vital for drug discovery protocols. Thus, the development of virtual libraries for this purpose is of utmost importance in a modern drug discovery program [21]–[23]. In this paper, we demonstrate that a “drug-like” and diverse library, constituting > 500 natural products collected from diverse regions within the continent of Africa (for which samples are directly available for screening), might be both a significant starting point and a useful tool in the drug discovery process, if virtual and biological screenings are employed.

What is p-ANAPL? The pan-African natural products library (p-ANAPL) is a consortium of natural product collections isolated from African biota and owned by scientists and/or groups of scientists working in African institutions. The p-ANAPL project was established in April 2009 by a group of natural product scientists from Botswana, Cameroon, Ethiopia, Kenya, Sudan and Tanzania, to build up the collection of natural products and to create an enhanced and efficient platform for their biological screening, and thereby promote the discovery of useful products from them. The p-ANAPL library is associated with the Center for Scientific Research, Indigenous Knowledge and Innovation (Cesriki), in the University of Botswana.

Materials and Methods

Collection method of the physical library

The compounds in the p-ANAPL collection were mainly obtained through the Network for Analytical and Bioassay Services in Africa (NABSA) which is based at the Department of Chemistry at the University of Botswana. NABSA was initiated in 1992 to share facilities with other African scientist working in facility-constrained environments. NABSA promotes and continues to do so, short-term visits by African researchers and postgraduate students to the University of Botswana and providing access to spectroscopic services (High Field NMR, High Resolution MS, Chiroptical Spectroscopy) for less privileged researchers and institutions. NABSA therefore promoted intra-African cooperation amongst several countries including Cameroon, Nigeria, Kenya, Ethiopia, Tanzania, the Democratic Republic of Congo and Zimbabwe. It is the policy of NABSA that visiting researchers who isolate and elucidate the structures of natural products voluntarily deposit specimens to establish a collection of natural products with full entitlement to any intellectual property (IP) that may arise from subsequent assay of these substances. It is mostly these compounds that formed the initial basis for p-ANAPL.

The first p-ANAPL consortium meeting in Gaborone, Botswana in 2009, defined its main goal as systematically organising compounds in a well-characterized natural product collection with the potential to expand and to use for bioactivity screening purposes. One of the co-authors of this paper, G. W. F., undertook the task of characterization and cataloguing all compounds forming the basis of this paper during a post-doctoral term at the University of Botswana from 2011 to 2013. The challenges and lessons learnt from this exercise will be published elsewhere.

Generation of 3D models and in silico determination of molecular descriptors

The 3D structures of the compounds had been sketched and energy minimisation subsequently using the graphical user interface (GUI) of the MOE software [24] running on a Linux workstation with a 3.5GHz Intel Core2 Duo processor. The 3D structures were generated using the builder module of MOE and energy minimization was subsequently carried out using the MMFF94 force field [25] until a gradient of 0.01 kcal/mol was reached. The 3D structures of the compounds were then saved as.mol2 files subsequently included into a MOE database (.mdb) file and converted to other file formats (.sdf,.mol,.mol2 and.ldb), which are suitable for use in several virtual screening workflow protocols. The molar weight (MW), number of rotatable bonds (NRB), lipophilicity parameter (log P), number of hydrogen bond acceptors (HBA), number of hydrogen bond donors (HBD) and number of Lipinski violations were calculated using the molecular descriptor calculator included in the QuSAR module of the MOE package [24]. A further treatment was carried out to ensure that the protanation states of the compounds were correct with respect to biological pH (strong acids were deprotonated, strong bases protonated and metal atoms removed). A maximum of 10 tautomers were generated for each molecule in the dataset.

Results and Discussion

Unique natural compounds within the library

The compounds in the p-ANAPL collection comprise about thirty (30) different classes of compounds that reflect a substantial chemical diversity even in a relatively small collection. It shows the potential of African natural products in terms of molecular diversity which has been noted to be often limited in large compound libraries used in high throughput screening (HTS) format [21], [26]. Molecular diversity goes hand in hand with bioactivity diversity, which is an important consideration for drug discovery. A recent study based on a database of 197 201 natural products (NP) from plants, animals and microorganisms of mainly Chinese origin revealed that there was a large overlap between natural products and FDA-approved drugs, in terms of chemical space, indicating the potential of NPs to become lead compounds [27]. An analysis of NP-target networks suggested a high degree of polypharmacology associated with NPs. This seems to confirm comparisons by Feher and Schmidt, showing that overall, NPs are more similar to drugs than compounds obtained from combinatorial chemistry [28]. Close to half of the p-ANAPL compounds are flavonoids (mono-, di-, tri- prenylated and geranylated flavones, flavanones, chalcones; homoisoflavonoids, rotenoids, biflavonoids, bichalcones, etc.), a very versatile class of NPs, which have been associated with a relatively large margin of safety for therapeutic use in humans and which have recently been reported to display diverse anti-HIV bioactivities [29]. More than two hundred (200) flavonoids, such as proropensin (1) or Dinklagin B (2), isolated primarily from Dorstenia species [30], [31], Figure 4, are available in good quantities in the p-ANAPL collection. The p-ANAPL collection, furthermore, contains about forty (40) homoisoflavonoids isolated mainly from Hyacinthaceae species such as Scilla nervosa [32], which constitutes a quite unique contingent and the largest collection of compounds representing this class. An initial assessment of commercial availability of p-ANAPL compounds from major suppliers, PubChem entries and specifically dedicated webpages for supplies of chemicals revealed that approximately half of the compounds are currently not commercially available. Amongst them are compounds with very unique properties. Examples are the phenyl anthraquinones, Gaboroquinone A and B (3 and 4), with interesting stereochemistry and good antiplasmodial (IC50 = 4.2 μg mL−1 against Plasmodium falciparum (NF54)) and anti-trypanosomal (IC50 = 5.1 μg mL−1 against Trypanosoma b. rhodesiense) activities [33].

Furthermore, the p-ANAPL collection contains a series of sulphate derivatives of 6′-O-sulphated phenylanthraquinones, which are unique as natural products with sulphate functional groups are considered to be very rare in plants [34]. The p-ANAPL also contains the largest collections of isofuranonapthoquinones. Broussoflavonol B (5), contained in the p-ANAPL collection, is a recently identified novel anti-cancer agent which has been shown to inhibit growth of estrogen-negative breast cancer MDA-MB-231 cells, a particularly more malignant and aggressive form of cancer cells which account for about one third of breast cancers [35]. Additionally, the collection contains julocrotine (6), an alkaloid, which has been shown to be effective against Leishmania (L.) amazonensis amastigotes, the causative agent for cutaneous leishmaniasis (IC50 = 19.8 μM) [36], and a series of compounds dorsmanin A to I, isolated from the genus Dorstenia, which showed versatile antimicrobial activities [37].

A substantial number of NP databases have been established in the last years [21], [38]–[43], which provided valuable information for assessments of chemical space occupation for virtual screening purposes. To a much lesser extent non-commercial actual NP collections with proper characterization and traceability of samples do exist for the African continent. An example of a similar effort is the Natural Product Discovery Institute (NPDI, formerly Merck U.S. Natural Products Library and Schering-Plough Legacy Culture Collection, www.npdi-us.org) and the Developmental Therapeutics Program (DTP) collection of the National Cancer Institute/National Institutes of Health (NCI/NIH) in the US (http://dtp.nci.nih.gov/branches/npb/repository.html), with only a partial coverage of African natural resources. Both collections are not based on the African continent. Close to half of the p-ANAPL compounds are available in amounts of more than 10 mg which allows biological screening to some extent. Therefore, the p-ANAPL collection is unique in the sense that it has brought together NPs from African researchers that remain under the control of African institutions. We hope that this first virtual characterization of our collection triggers interest in that collection and, above all, serves as a springboard to extend p-ANAPL to other African research institutions to transform it into a true pan-African collection.

The chemical composition of the p-ANAPL library

A pie chart, showing the distribution of the five hundred and thirty four (534) samples currently lodged in the p-ANAPL library is shown in Fig. 1. As previously mentioned, a majority of the compound samples are flavonoids (flavanones, biflavonoid, isoflavones, homoisoflavonoids, and chalcones), constituting 47.19% of compound samples. This is followed by terpenoids (12.54%), quinones (11.61%), monoaromatic compounds (7.87%), alkaloids (5.43%), coumarins (3.74%), steroids (2.81%) and xanthones (2.43%). The remaining compound types each represented < 2% of the total sample collection.

[Figure omitted. See PDF.]

Figure 1. Pie chart showing the distribution by compound type.

https://doi.org/10.1371/journal.pone.0090655.g001

Assessment of its “drug-likeness” and “lead-likeness” potential

Lipinski’s “Rule of Five” [44] (ro5) is a very useful tool for assessing chemical compound libraries to be used in drug discovery programs. This “rule” was derived from chemical libraries from the World Dug Index (WDI), as a criterion to evaluate likely oral bioavailability [44]–[45]. However, the highly valuable class of NPs was initially omitted, since Lipinski had postulated that the ro5 was not respected by NPs. NP libraries have however been previously analysed comparatively using the ro5 in order to have a rough idea of the extent of “drug-likeness” of a compound library to be used in virtual screening [21], [41]–[43]. It is on these grounds that Lipinski’s criteria [44] are often used for the evaluation of “drug-likeness” of compounds within the designed libraries. Thus, Lipinski’s ro5 is often regarded as a useful filter for the elimination of compounds not likely to be orally available in the early stages of drug discovery protocols [45]. In summary, Lipinski’s ro5 defines a “drug-like” molecule as one with high likelihood to be orally available, for which the molar weight (MW) ≤ 500 Daltons (Da), the logarithm of the octan-1-ol/water partition coefficient (log P(o/w) or log P) ≤ 5, the number of hydrogen bond acceptors (HBA) ≤ 10 and the number of hydrogen bond donors (HBD) ≤ 5. An additional rule for the number of rotatable bonds (NRB ≤ 5) is often included to the ro5. An evaluation of “lead-likeness” is often carried out using more stringent criteria defined by Oprea et al [46]–[49]. The “Rule of 3.5” for “lead-like” compounds is defined as: 150 ≤ MW ≤ 350; log P(o/w) ≤ 4; HBD ≤ 3; HBA ≤ 6). The criteria for “fragment-like” libraries is defined by Verdonk et al [50]. The criteria (also referred to as the “Rule of 2.5”) are such that MW ≤ 250; –2 ≤ log P(o/w) ≤ 3; HBD < 3; HBA < 6; NRB < 3. A pairwaise scatter plot for the physico-chemical parameters defining the ro5 are shown in Fig. 2. These plots show that the regions of highest population density of points fall within the Lipinski compliance regions (LCR).

[Figure omitted. See PDF.]

Figure 2. Pair-wise distribution of Lipinski parameters for the p-ANAPL library; (A) log P against MW, (B) HBA against MW, (C) HBD against MW and (D) NRB against MW.

https://doi.org/10.1371/journal.pone.0090655.g002

Fig. 3 shows the distribution of violations and compliance of Lipinski parameters in the p-ANAPL library, while the maxima, minima and mean values for each parameter are shown in Table 1, along with those of the “drug-like”, “lead-like” and “fragment-like” subsets, defined following the aforementioned criteria. An analysis of the total library showed that about 67% of the compounds had no Lipinski violations and almost 85% of the compounds showed < 2 violations (Fig 3A). The ‘drug-likeness’ is further highlighted by an analysis of individual parameters (MW, log P(o/w), HBA and HBD). The MW values showed a Gaussian distribution with a peak value located between 301 and 400 Da, Fig. 3B. This interval corresponded to about 34% of the compounds currently included in the p-ANAPL database. In addition, only about 14% of all the compounds showed MW falling outside the recommended range for the ro5 (> 500 Da). The mean values for MW ( =  370 Da), further highlighted the “drug-like” nature of our library, the mean MW of a typical drug being equal to 310 Da [51]. The log P(o/w) distribution curve (Fig. 3C) showed a very steep Gaussian shape with a peak centred at 2.5 log P(o/w) units, with only 9.5% of the compounds having log P(o/w) > 5. Only six compounds (compounds 8 to 13) showed log P(o/w) > 8 units. These constitute stilbenoid and long chain alkyl esters of some pentacyclic triterpenoids (Fig. 5). The distribution curves for HBA and HBD respectively rose rapidly to maxima of 6 acceptors (corresponding to 17.66% of the compounds) and 2 donors (corresponding to 18.96% of the compounds). The two graphs then fell to 24 acceptors and 13 donors respectively (Fig. 3D,E). Moreover, only 10.41% of the compounds had HBA > 10 and only 12.45% of the compounds showed HBD > 5. The graph of the NRB showed two cusps; at 1 and 3 RBs (Fig. 3F) and rapidly fell to 21 rotatable bond (RBs), with only 14.68% of the compounds having NRB > 5. The mean values of all the Lipinski parameters indicate a high probability of finding ‘drug-like’ and ‘lead-like’ compounds within the p-ANAPL library (Table 1).

[Figure omitted. See PDF.]

Figure 3. Distribution of Lipinski parameters for the p-ANAPL library; (A) Bar chart showing the number of Lipinski violations, (B) Bar chart showing the MW parameter, (C) Plot of the lipophilicity parameter, (D) Plot of the HBA parameter, (E) Plot of the HBD parameter, and (F) Plot of the NRB parameter.

https://doi.org/10.1371/journal.pone.0090655.g003

[Figure omitted. See PDF.]

Figure 4. Chemical structures of selected bioactive natural products from the p-ANAPL library (1 to 7).

https://doi.org/10.1371/journal.pone.0090655.g004

[Figure omitted. See PDF.]

Figure 5. Chemical structures of natural products from the p-ANAPL library with log P > 8 units.

https://doi.org/10.1371/journal.pone.0090655.g005

[Figure omitted. See PDF.]

Table 1. Summary of physico-chemical properties (often used to predict ‘drug-likeness’) of the compounds within the p-ANAPL library.

https://doi.org/10.1371/journal.pone.0090655.t001

A discussion of pharmacological profiles of selected compounds

The chemical structures of selected compounds with interesting pharmacological profiles from the p-ANAPL library have been shown in Fig. 6. Esters of the monoaromatic mandelic acid (14) have been involved in the modulation of the enantioselectivity of lipases via controlled immobilization on glutaraldehyde supports [52]. In addition to previously known biological activities, Zofou et al. recently demonstrated the anti-malarial properties of the flavonoids quercitrin (15) and quercetin (16), derived from Dacryodes edulis (Burseraceae), a plant currently used in the treatment of malaria and fevers in West Cameroon [53]. Compound 15 exhibited IC50 values of 5.96 and 2.26 μg mL−1, against the 3D7 and Dd2 strains of Plasmodium falciparum respectively, while compound 16 exhibited IC50 values of 6.07 and 5.91 μg mL−1, against 3D7 and Dd2 respectively. The anthraquinones emodin (18) and aloe emodin (20) have demonstrated interesting anticancer properties [54]–[57], while chrysophanol (19) has shown the ability to induce necrosis through the production of ROS and alteration of ATP levels in J5 human liver cancer cells [58], in addition to its antidiabetic [59], anti-inflammatory [60] and antimicrobial [61]–[62] properties. Moreover, chrysophanol-8-O-glucoside, has shown antiplatelet and anticoagulant activities [63]. It is also proven that physcion (17) and emodin (18) exhibit antibacterial properties against three Bacillus species, emodin exhibiting MICs in the range 0.5 – 2.0 μg mL−1 [64]. Both compounds inhibited Pseudomonas aeruginosa, emodin being more effective, showing an MIC of 70 μg mL−1 [64].

[Figure omitted. See PDF.]

Figure 6. Chemical structures of selected bioactive natural products from the p-ANAPL library (14 to 36).

https://doi.org/10.1371/journal.pone.0090655.g006

Shikimic acid (21) is known to play a key role in aromatic biosynthesis [65]. This has been exploited in the design of modified shikimic acid derivatives as potential anti-tubercular agents, inhibiting Mycobacterium tuberculosis shikimate kinase enzyme [66]. The quinone knipholone (22) has been derived from Kniphophia foliosa (Asphodelaceae) and is known to exhibit anti-malarial activity [67]. The antiplasmodial activity of knipholone is IC50 = 0.38 μM compared to Chloroquine 0.09 μM [68]. The chalcones bartericin A (23), stipulin (24), kanzonol B (26) and 4-hydroxylonchocarpin (27) have also shown anti-malarial activities [69]. while caffein (25) is a common psychostimulant. Caffeine was also reported to inhibit HIV-1 transduction of non-dividing cells [70]. Knipholone anthrone (28) has shown anti-cancer activity, with IC50 of 0.5 to 3.3 μM against leukaemic and melanocyte cancer cell lines [71]. Ellagic acid (29), isolated from the leaves of Alchornea cordifolia (Euphorbiaceae) by Banzouzi et al. [72], also showed good activity against Plasmodium berghei in mice with an ED50 in the range of 0.2 – 0.151 μg mL−1. The coumarin umbelliferone (30) has exhibited analgesic and anti-inflammatory [73], antihyperglycemic [74] and antioxidant [75]–[76] activities. The derivatives of umbelliferone have proven to be potent 5α-reductase type 1 inhibitors and provided the basis for further development for the treatment of human endocrine disorders associated with overproduction of DHT by 5α-reductase type 1 [77]. Aesculetin (31) induces apoptosis through a ROS-mediated mitochondrial dysfunction pathway in human cervical cancer cells [78], in addition to its antimicrobial [79], and anti-inflammatory [80] activities.

Yohimbine (32) is an alkaloid extracted from Pausinystalia johimbe, a plant commonly used to treat erectile dysfunction in ATM in West and Central Africa [2]. Both yohimbine [81] and its hydrochloride [82] have proven to be potent in the treatment of erectile dysfunction by preferential blockade of presynaptic α-adrenoceptors in rabbits [83]. Starke et al. demonstrated that yohimbine is more potent in blocking the presynaptic than the postsynaptic α-adrenoceptors of the artery. Cannabidiol (33) is a non-psychotropic component of Cannabis with possible therapeutic use as an anti-inflammatory drug [84]. It is a neuroprotective antioxidant [85], an oral anti-arthritic therapeutic in murine collagen-induced arthritis [86], which induces anxiety and psychotic-like symptoms in healthy subjects [87]. The plant alkaloid piperine (34) exhibits several biological activities, including inhibition of ethidium bromide efflux in Mycobacterium smegmatis [88], inhibition of Rv1258c, a putative multidrug efflux pump of Mycobacterium tuberculosis [89], selective inhibition of CYP3A4 [90], immunomodulatory and antitumor activities [91], inhibition of human P-glycoprotein [92], anti-inflammatory and antiarthritic effects [93]–[94], and inhibition of mammosphere formation [95]. Compound 35 (p-hydroxycinnamic acid) is known to exhibit anti-malarial activity [96], stimulate bone formation/inhibit bone resorption [97] and act as a natural mediator for laccase oxidation of recalcitrant compounds [98]. Reynosin (36) has exhibited neuronal toxicity protection in Parkinson’s disease models [99].

Accessibility of the compounds

The virtual compound library is available as supplementary data accompanying this publication (Dataset S1). Updated versions of this dataset will be subsequently be available on request for non-commercial purposes through the corresponding author L. M. M. of this article. The 3D structures have been generated based on the chemical structure reported in the literature, and treated as previously described [41]–[43]. All requests for compounds samples for biological assays should be formally addressed to the corresponding author K. A. M. of this article.

Conclusions

We have recently reported the development of natural product virtual libraries for African medicinal plants [41]–[43]. However, to the best of our knowledge, the p-ANAPL library constitutes the largest collection of physical samples of NPs derived from African medicinal plants, which is available for biological screening. In addition, the virtual library is provided in several file formats (mol2, sdf, mdb, and ldb), which are readable by several drug discovery software. These could be useful in virtual screening campaigns. Concerning the physical samples, the purity of each compound was tested by measurement of melting points, and confirmed to be > 95% pure before including in the database. The stability has been assured by keeping compounds at below 0°C in freezers. The uniqueness of the p-ANAPL library also lies in the fact that it is the largest collection of NPs with physical samples, specifically derived from African medicinal plants.

Supporting Information

[Figure omitted. See PDF.]

Dataset S1.

Electronic Supplementary Information (ESI) available: [Low energy 3D structures of compounds currently included in the p-ANAPL library for virtual screening, along with the generated tautomers, drug-like, lead-like and fragment-like subsets]. Correspondence and requests for materials should be addressed to LMM and KAM for the virtual library and physical samples respectively.

https://doi.org/10.1371/journal.pone.0090655.s001

(RAR)

Acknowledgments

We wish to sincerely acknowledge the colleagues and post-graduate students: J. Bipa, S. Famuyiwa, D. Kapche, S. Khalid, J. J. Magadula, L. K. Mdee, J. Mutanyatta, N. Ngamga, M. Ngameni, and O. Shetonde who have contributed compounds to the p-ANAPL collection.

Author Contributions

Conceived and designed the experiments: LMM KAM BMA BTN AOO. Performed the experiments: FNK PAO GWF MB JCN. Analyzed the data: GWF FNK PAO MB. Contributed reagents/materials/analysis tools: LMM KAM BMA BTN AOO FNK PAO GWF BTN. Wrote the paper: FNK GWF KAM LMM BMA AOO.

Citation: Ntie-Kang F, Amoa Onguéné P, Fotso GW, Andrae-Marobela K, Bezabih M, Ndom JC, et al. (2014) Virtualizing the p-ANAPL Library: A Step towards Drug Discovery from African Medicinal Plants. PLoS ONE 9(3): e90655. https://doi.org/10.1371/journal.pone.0090655

References

1. Vasisht K, Kumar V (2002) Trade and Production of Herbal Medicines and Natural Health Products. ICS-UNIDO, Trieste, pp 3.

2. Vasisht K, Kumar V (2004) Compendium of Medidincal and aromatic plants. Volume 1: Africa. ICS-UNIDO, Trieste, pp 23–56.

3. Harvey AL (2008) Natural products in drug discovery. Drug Discov Today 13: 894–901.

4. Newman DJ (2008) Natural products as leads to potential drugs: an old process or the new hope for drug discovery? J Med Chem 51: 2589–2599.

5. Efange SMN (2002) Natural products: a continuing source of inspiration for the medicinal chemist. In M.M. Iwu & J.C. Wootton (Eds.), Advances in Phytomedicine, Vol. 1, Ethnomedicine and Drug Discovery (61–69). Amsterdam, The Netherlands: Elsevier Science.

6. Hostettmann K, Marston A, Ndjoko K, Wolfender JL (2000) The potential of African plants as a source of drugs. Curr Org Chem 4: 973–1010.

7. Chibale K, Davies-Coleman M, Masimirembwa C (2012). Drug discovery in Africa: impacts of genomics, natural products, traditional medicines, insights into medicinal chemistry, and technology platforms in pursuit of new drugs. Springer.

8. Titanji VPK, Zofou D, Ngemenya M (2008) The antimalarial potential of medicinal plants used for the treatment of malaria in Cameroonian folk medicine. Afr J Trad CAM 5(3): 302–321.

9. Schwikkard S, van Heerden FR (2002) Antimalarial activity of plant metabolites. Nat Prod Rep 19: 675–692.

10. Kuete V, Efferth T (2010) Cameroonian medicinal plants: pharmacology and derived natural products. Front Pharmacol 1: 123.

11. Ntie-Kang F, Lifongo LL, Mbaze LM, Ekwelle N, Owono LCO, et al. (2013) Cameroonian medicinal plants: a bioactivity versus ethnobotanical survey and chemotaxonomic classification. BMC Complement Altern Med 13(1): 147.

12. Magadula JJ, Erasto P (2009) Bioactive natural products derived from the East African flora. Nat Prod Rep 26: 1535–1554.

13. Onguéné PA, Ntie-Kang F, Lifongo LL, Ndom JC, Sippl W, et al. (2013) The potential of anti-malarial compounds derived from African medicinal plants. Part I: A pharmacological evaluation of alkaloids and terpenoids. Malar J 2013 12: 449.

14. Ntie-Kang F, Onguéné PA, Lifongo LL, Ndom JC, Sippl W, et al.. (2013) The potential of anti-malarial compounds derived from African medicinal plants. Part II: A pharmacological evaluation of non alkaloids and non terpenoids. Malar J submitted.

15. Kuete V (2010) Potential of Cameroonian plants and derived products against microbial infections: a review. Planta Med 76: 1479–1491.

16. Zofou D, Ntie-Kang F, Sippl W, Efange SMN (2013) Bioactive natural products derived from the Central African flora against neglected tropical diseases and HIV. Nat Prod Rep 30: 1098–1120.

17. Koehn FE, Carter GT (2005) The evolving role of natural products in drug discovery. Nat Rev Drug Discov 4: 206–220.

18. Li JWH, Vederas JC (2009) Drug discovery and natural products: end of an era or an endless frontier? Science 325: 161–165.

19. Lam KS (2007) New aspects of natural products in drug discovery. Trends Microbiol 15: 279–289.

20. Pan L, Chai HB, Kinghorn AD (2013) Discovery of new anticancer agents from higher plants. Front. Biosci. (Schol Ed.) 4: 142–156.

21. Quinn RJ, Carroll AR, Pham MB, Baron P, Palframan ME, et al. (2008) J Nat Prod. 71: 464–468.

22. Irwin JJ, Shoichet BK (2005) ZINC – A free database of commercially available compounds for virtual screening. J Chem Inf Model 45(1): 177–182.

23. Irwin JJ, Sterling T, Mysinger MM, Bolstad ES, Coleman RG (2012) ZINC: A free tool to discover chemistry for biology. J Chem Inf Model 52: 1757–1768.

24. Chemical Computing Group Inc: Molecular Operating Environment Software. Montreal; 2010.

25. Halgren TA (1996) Merck molecular forcefield. J Comput Chem 17: 490–641.

26. Dobson CM (2004) Chemical space and biology. Nature 432(7019): 824–828.

27. Gu J, Gui Y, Yuan G, Lu HZ, Xu X (2013) Use of natural products as chemical library for drug discovery and network pharmacology. PLOS ONE 8(4): e62839.

28. Feher M, Schmidt JM (2003) Property distributions: differences between drugs, natural products, and molecules from combinatorial chemistry. J Chem Inf Comput Sci 43: 218–227.

29. Andrae-Marobela K, Fotso WG, Okatch H, Majinda RRT (2013) Polyphenols: A diverse class of multi-target anti HIV-agents. Current Drug Metabolism 14: 392–413.

30. Abegaz BM, Ngadjui BT, Dongo E, Bezabih MT (2000) Chemistry of the genus Dorstenia. Current Contents in Organic Chemistry 4: 1079–1090.

31. Abegaz BM, Ngadjui BT, Dongo E, Ngameni B, Nindi N, Bezabih MT (2002) Chalcones and other constituents of Dorstenia proropens and D. zenkeri. Phytochemistry 59(8): 877.

32. Silayo A, Ngadjui BT, Abegaz BM (1999) Homoisoflavonoids and stilbenes from the bulbs of Scilla nervosa.. Phytochemistry 52: 947–955.

33. Abegaz BM, Bezabih M, Msuta T, Brun R, Menche D, et al. (2002) Gaboroquinones A and B and 4’-O-demethylknipholone-4’-O-beta-D-glucopyranoside, phenylanthraquinones from the roots of Bulbine frutescens. J Nat Prod 65(8): 1117–1121.

34. Mutanyatta J, Bezabih M, Abegaz BM, Dreyer M, Brun R, et al. (2005) The first 6’-O-sulfated phenylanthraquinones: Isolation from Bulbine frutescens, structural elucidation, enantiomeric purity, and partial synthesis. Tetrahedron 61: 8475–8484.

35. Guo MX, Wang ML, Deng H, Zhang XT, Wang ZY (2013) A novel anticancer agent Broussoflavonol B downregulates estrogen receptor (ER)-alpha36 expression and inhibits growth of ER-negative breast cancer MDA-MB-231 cells. Eur J Pharmacol 714: 56–64.

36. Guimarães LR, Rodrigues AP, Marinho PS, Muller AH, Guilhon GM, et al. (2010) Activity of the julocrotine, a glutarimide alkaloide from Croton pullei var. glabrior, on Leishmania (L.) amazonensis. Parasitol Res 107(5): 1075–1081.

37. Mbaveng AT, Kuete V, Ngameni B, Beng VP, Ngadjui BT, et al. (2012) Antimicrobial activities of the methanol extract and compounds from the twigs of Dorstenia mannii (Moraceae). BMC Complement Altern Med 12: 83.

38. Yongye AB, Waddell J, Medicna-Franco JL (2012) Molecular scaffold analyisis of natural product databases in the public domain. Chem Biol Drug Des 80(5): 717–724.

39. Fullbeck M, Michalsky E, Dunkel M, Preissner R (2006) Natural products: sources and databases. Nat Prod Rep 23: 347–356.

40. Dunkel M, Fullbeck M, Neumann S, Preissner R (2006) SuperNatural: a searchable database of available natural compounds. Nucleic Acids Res 34: D678–D683.

41. Ntie-Kang F, Mbah JA, Mbaze LM, Lifongo LL, Scharfe M, et al. (2013) CamMedNP: Building the Cameroonian 3D structural natural products database for virtual screening. BMC Complement Altern Med 13: 88.

42. Ntie-Kang F, Onguéné PA, Scharfe M, Owono LCO, Megnassan E, et al. (2014) ConMedNP: A natural product library from Central African medicinal plants for drug discovery. RSC Adv 2014 4: 409–419.

43. Ntie-Kang F, Zofou D, Babiaka SB, Meudom R, Scharfe M, et al. (2013) AfroDb: A select highly potent and diverse natural product library from African medicinal plants. PLOS ONE 8: e78085.

44. Lipinski CA, Lombardo F, Dominy BW, Feeney PJ (2001) Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv Drug Delivery Rev 46: 3–26.

45. Lipinski CA (2000) Drug-like properties and the causes of poor solubility and poor permeability. J Pharmacol Toxicol Methods 44: 253–249.

46. Teague SJ, Davis AM, Leeson PD, Opea TI (1999) The design of leadlike combinatorial libraries. Angew Chem, Int Ed 38: 3743–3748.

47. Hann MM, Oprea TI (2004) Pursuing the leadlikeness concept in pharmaceutical research. Curr Opin Chem Biol 8: 255–263.

48. Oprea TI (2002) Current trends in lead discovery: are we looking for the appropriate properties J Comput-Aided Mol Des. 16: 325–334.

49. Schneider G (2002) Trends in virtual computational library design. Curr Med Chem 9: 2095–2102.

50. Verdonk ML, Cole JC, Hartshorn ML, Murray CW, Taylor RD (2003) Improved protein-ligand docking using GOLD. Proteins 52: 609–623.

51. Khanna V, Ranganathan S (2011) Structural diversity of biologically interesting datasets: a scaffold analysis approach. J Cheminform 3: 30.

52. Palomo JM, Fernandez-Lorente G, Mateo C, Ortiz C, Fernandez-Lafuente R, et al. (2002) Modulation of the enantioselectivity of lipases via controlled immobilization and medium engineering: hydrolytic resolution of mandelic acid esters. Enzyme and Microbial Technology 31: 775–783.

53. Zofou D, Tematio EL, Ntie-Kang F, Tene M, Ngemenya MN, et al. (2013) New antimalarial hits from Dacryodes edulis (Burseraceae) - Part I: isolation, in vitro activity, in silico “drug-likeness” and pharmacokinetic profiles. PLOS ONE 8(11): e79544.

54. Pecere T, Gazzola MV, Mucignat C, Parolin C, Dalla F, et al. (2000) Aloe-emodin is a new type of anticancer agent with selective activity against neuroectodermal tumors. Cancer Res 60: 2800–2804.

55. Kuo PL, Lin TC, Lin CC (2002) The antiproliferative activity of aloe-emodin is through p53-dependent and p21-dependent apoptotic pathway in human hepatoma cell lines. Life Sci 71: 1879–1892.

56. Cárdenas C, Quesada AR, Medina MA (2006) Evaluation of the anti-angiogenic effect of aloe-emodin. Cell Mol Life Sci 63: 3083–3089.

57. Chen YY, Chiang SY, Lin JG, Ma YS, Liao CL, et al. (2010) Emodin, aloe-emodin and rhein inhibit migration and invasion in human tongue cancer SCC-4 cells through the inhibition of gene expression of matrix metalloproteinase-9. Int J Oncology 2010 36: 1113–1120.

58. Lu CC, Yang JS, Huang AC, Hsia TC, Chou ST, et al. (2010) Chrysophanol induces necrosis through the production of ROS and alteration of ATP levels in J5 human liver cancer cells. Mol Nutr Food Res 54(7): 967–976.

59. Lee MS, Sohn CB (2008) Anti-diabetic properties of chrysophanol and its glucoside from rhubarb rhizome. Biol Pharm Bull 31(11): 2154–2157.

60. Kim SJ, Kim MC, Lee BJ, Park DH, Hong SH, et al. (2010) Anti-Inflammatory activity of chrysophanol through the suppression of NF-κB/caspase-1 activation in vitro and in vivo. Molecules 15: 6436–6451.

61. Coopoosamy RM, Magwa ML (2006) Antibacterial activity of aloe emodin and aloin A isolated from Aloe excelsa. Afr J Biotech 5(11): 1092–1094.

62. García-Sosa K, Villarreal-Alvarez N, Lübben P, Peña-Rodríguez LM (2006) Chrysophanol, an antimicrobial anthraquinone from the root extract of Colubrina greggii. J Mex Chem Soc 50(2): 76–78.

63. Seo EJ, Ngoc TM, Lee SM, Kim YS, Jung YS (2012) Chrysophanol-8-O-glucoside, an anthraquinone derivative in rhubarb, has antiplatelet and anticoagulant activities. J Pharmacol Sci 118: 245–254.

64. Basu S, Ghosh A, Hazra B (2005) Evaluation of the antibacterial activity of Ventilago madraspatana Gaertn., Rubia cordifolia Linn. and Lantana camara Linn.: isolation of emodin and physcion as active antibacterial agents. Phytother Res 19: 888–894.

65. Davis BD (1951) Aromatic biosynthesis. I. The role of shikimic acid. J Biol Chem 191: 315–325.

66. Blanco B, Prado V, Lence E, Otero JM, Garcia-Doval C, et al. (2013) Mycobacterium tuberculosis shikimate kinase inhibitors: design and simulation studies of the catalytic turnover. J Am Chem Soc 135(33): 12366–12376.

67. Dagne E, Steglich W (1984) Knipholone: a unique antraquinone derivative from Kniphofia foliosa. Phytochemistry 23: 1729–1731.

68. Bringmann G, Menche D, Bezabih M, Abegaz BM, Kaminsky R (1999) Antiplasmodial activity of knipholone and related natural phenylanthraquinones. Planta Med 65: 757–758.

69. Ngameni B, Watchueng J, Boyom FF, Keumedjio F, Ngadjui BT, et al. (2007) Antimalarial prenylated chalcones from the twigs of Dorstenia barteri var. subtriangularis. Arkivoc xiii: 116–123.

70. Daniel R, Marusich E, Argyris E, Zhao RY, Skalka AM, et al. (2005) Caffeine inhibits human immunodeficiency virus type 1 transduction of nondividing cells. J Virol 79(4): 2058–2065.

71. Habtemariam S (2010) Knipholone anthrone from Kniphofia foliosa induces a rapid onset of necrotic cell death in cancer cells. Fitoterapia 81: 1013–1019.

72. Banzouzi JT, Prado R, Menan H, Valentin A, Roumestan C, Mallie M, et al. (2002) Studies on medicinal plants of Ivory coast: Investigation of an active constituent. J Ethnopharmacol 81(3): 399–341.

73. Lino CS, Taveira ML, Viana GSB, Matos FJA (1997) Analgesic and antiinflammatory activities of Justicia pectoralis Jacq and its main constituents: coumarin and umbelliferone. Phytother Res 11: 211–215.

74. Ramesh B, Pugalendi KV (2006) Antihyperglycemic effect of umbelliferone in streptozotocin-diabetic rats. J Med Food 9(4): 562–566.

75. Ramesh B, Pugalendi KV (2006) Antioxidant role of umbelliferone in STZ-diabetic rats. Life Sciences 79(3): 306–310.

76. Singh R, Singh B, Singh S, Kumar N, Kumar S, et al. (2010) Umbelliferone – An antioxidant isolated from Acacia nilotica (L.) Willd. Ex. Del. Food Chemistry 120(3): 825–830.

77. Fan GJ, Mar W, Park MK, Choi EW, Kim K, et al. (2001) A novel class of inhibitors for steroid 5α-Reductase: synthesis and evaluation of umbelliferone derivatives. Bioorg Med Chem Lett 11(7): 2361–2363.

78. Yang J, Xiao YL, He XR, Qiu GF, Hu XM (2010) Aesculetin-induced apoptosis through a ROS-mediated mitochondrial dysfunction pathway in human cervical cancer cells. J Asian Nat Prod Res 12(3): 185–193.

79. Jurd L, Corse J, King Jr AD, Bayne H, Mihara K (1971) Antimicrobial properties of 6,7-dihydroxy-,7,8-dihydroxy-,6-hydroxy- and 8-hydroxycoumarins. Phytochemistry 10(2): 2971–2974.

80. Silván AM, Abad MJ, Bermejo P, Sollhuber M, Villar A (1996) Antiinflammatory activity of coumarins from Santolina oblongifolia. J Nat Prod 59(12): 1183–1185.

81. Morales A, Condra M, Owen JA, Surridge DH, Fenemore J, et al. (1987) Is yohimbine effective in the treatment of organic impotence? Results of a controlled trial. J Urol 137(6): 1168–1172.

82. Susset JG, Tessier CD, Wincze J, Bansal S, Malhotra C, Schwacha MG (1989) Effect of yohimbine hydrochloride on erectile impotence: a double-blind study. J Urol 141(6): 1360–1363.

83. Starke K, Borowski E, Endo T (1975) Preferential blockade of presynaptic α-adrenoceptors by yohimbine. Eur J Pharmacol 34(2): 385–388.

84. Bisogno T, Hanuš L, De Petrocellis L, Tchilibon S, Ponde DE, Brandi I, et al. (2001) Molecular targets for cannabidiol and its synthetic analogues: effect on vanilloid VR1 receptors and on the cellular uptake and enzymatic hydrolysis of anandamide. Brit J Pharmacol 134: 845–852.

85. Hampson AJ, Grimaldi M, Axelrod J, Wink D (1998) Cannabidiol and (-)Δ9-tetrahydrocannabinol are neuroprotective antioxidants. Proc Natl Acad Sci USA 95: 8268–8273.

86. Malfait AM, Gallily R, Sumariwalla PF, Malik AS, Andreakos E, et al. (2000) The nonpsychoactive cannabis constituent cannabidiol is an oral anti-arthritic therapeutic in murine collagen-induced arthritis. Proc Natl Acad Sci USA 97: 9561–9566.

87. Zuardi AW, Crippa JAS, Hallak JEC, Moreira FA, Guimarães FS (2006) Cannabidiol, a Cannabis sativa constituent, as an antipsychotic drug. Braz J Med Biol Res 39: 421–429.

88. Jin J, Zhang J, Guo N, Feng H, Li L, et al. (2011) The plant alkaloid piperine as a potential inhibitor of ethidium bromide efflux in Mycobacterium smegmatis. J Med Microbiol 60: 223–229.

89. Sharma S, Kumar M, Sharma S, Nargotra A, Koul S, Khan IA (2010) Piperine as an inhibitor of Rv1258c, a putative multidrug efflux pump of Mycobacterium tuberculosis. J Antimicrob Chemother 65, 1694–1701.

90. Volak LP, Ghirmai S, Cashman JR, Court MH (2008) Curcuminoids inhibit multiple human cytochromes P450, UDP-glucuronosyltransferase, and sulfotransferase enzymes, whereas piperine is a relatively selective CYP3A4 inhibitor. Drug Metab Dispos 36(8): 1594–1605.

91. Sunila ES, Kuttan G (2004) Immunomodulatory and antitumor activity of Piper longum Linn. and piperine. J Ethnopharmacol 90: 339–346.

92. Bhardwaj RK, Glaeser H, Becquemont L, Klotz U, Gupta SK, et al. (2002) Piperine, a major constituent of black pepper, inhibits human P-glycoprotein and CYP3A4. J Pharmacol Exp Ther 302: 645–650.

93. Bang JS, Oh DH, Choi HM, Sur BJ, Lim SJ, et al. (2009) Anti-inflammatory and antiarthritic effects of piperine in human interleukin 1β-stimulated fibroblast-like synoviocytes and in rat arthritis models. Arthritis Res Ther 11: R49.

94. Mujumdar M, Dhuley JN, Deshmukh VK, Raman PH, Naik SR (1990) Anti-inflammatory activity of piperine. Jpn J Med Sci Biol 43(3): 95–100.

95. Kakarala M, Brenner DE, Khorkaya H, Cheng C, Tazi K, et al. (2010) Targeting breast stem cells with the cancer preventive compounds curcumin and piperine. Breast Cancer Res. Treat 122(3): 777–785.

96. Zofou D, Tene M, Tane P, Titanji VPK (2012) Antimalarial drug interactions of compounds isolated from Kigelia africana (Bignoniaceae) and their synergism with artemether, against the multidrug-resistant W2mef Plasmodium falciparum strain. Parasitol Res 110: 539–544.

97. Lai YL, Yamaguchi M (2006) Phytocomponent p-hydroxycinnamic acid stimulates bone formation and inhibits bone resorption in rat femoral tissues in vitro. Mol Cell Biochem 292: 45–52.

98. Camarero S, Cañas AI, Nousiainen P, Record E, Lomascolo A, et al. (2008) P-hydroxycinnamic acids as natural mediators for laccase oxidation of recalcitrant compounds. Environ Sci Technol 42 (17): 6703–6709.

99. Ham A, Kim DW, Kim KH, Lee SJ, Oh KB, et al. (2013) Reynosin protects against neuronal toxicity in dopamine-induced SH-SY5Y cells and 6-hydroxydopamine-lesioned rats as models of Parkinson’s disease: Reciprocal up-regulation of E6-AP and down-regulation of α-synuclein. Brain Res 1524: 54–61.