Challenges in natural product-based drug discovery assisted with in silico-based methods.

The application of traditional medicine by humans for the treatment of ailments as well as improving the quality of life far outdates recorded history. To date, a significant percentage of humans, especially those living in developing/underprivileged communities still rely on traditional medicine for primary healthcare needs. In silico-based methods have been shown to play a pivotal role in modern pharmaceutical drug discovery processes. The application of these methods in identifying natural product (NP)-based hits has been successful. This is very much observed in many research set-ups that use rationally in silico-based methods in combination with experimental validation techniques. The combination has rendered the use of in silico-based approaches even more popular and successful in the investigation of NPs. However, identifying and proposing novel NP-based hits for experimental validation comes with several challenges such as the availability of compounds by suppliers, the huge task of separating pure compounds from complex mixtures, the quantity of samples available from the natural source to be tested, not to mention the potential ecological impact if the natural source is exhausted. Because most peer-reviewed publications are biased towards "positive results", these challenges are generally not discussed in publications. In this review, we highlight and discuss these challenges. The idea is to give interested scientists in this field of research an idea of what they can come across or should be expecting as well as prompting them on how to avoid or fix these issues.

An Analysis Based on Molecular Docking and Molecular Dynamics Simulation Study of Bromelain as Anti-SARS-CoV-2 Variants.

The rapid spread of a novel coronavirus known as SARS-CoV-2 has compelled the entire world to seek ways to weaken this virus, prevent its spread and also eliminate it. However, no drug has been approved to treat COVID-19. Furthermore, the receptor-binding domain (RBD) on this viral spike protein, as well as several other important parts of this virus, have recently undergone mutations, resulting in new virus variants. While no treatment is currently available, a naturally derived molecule with known antiviral properties could be used as a potential treatment. Bromelain is an enzyme found in the fruit and stem of pineapples. This substance has been shown to have a broad antiviral activity. In this article, we analyse the ability of bromelain to counteract various variants of the SARS-CoV-2 by targeting bromelain binding on the side of this viral interaction with human angiotensin-converting enzyme 2 (hACE2) using molecular docking and molecular dynamics simulation approaches. We have succeeded in making three-dimensional configurations of various RBD variants using protein modelling. Bromelain exhibited good binding affinity toward various variants of RBDs and binds right at the binding site between RBDs and hACE2. This result is also presented in the modelling between Bromelain, RBD, and hACE2. The molecular dynamics (MD) simulations study revealed significant stability of the bromelain and RBD proteins separately up to 100 ns with an RMSD value of 2 Å. Furthermore, despite increases in RMSD and changes in Rog values of complexes, which are likely due to some destabilized interactions between bromelain and RBD proteins, two proteins in each complex remained bonded, and the site where the two proteins bind remained unchanged. This finding indicated that bromelain could have an inhibitory effect on different SARS-CoV-2 variants, paving the way for a new SARS-CoV-2 inhibitor drug. However, more in vitro and in vivo research on this potential mechanism of action is required.

Pharmacoinformatic Investigation of Medicinal Plants from East Africa.

Medicinal plants have widely been used in the traditional treatment of ailments and have been proven effective. Their contribution still holds an important place in modern drug discovery due to their chemical, and biological diversities. However, the poor documentation of traditional medicine, in developing African countries for instance, can lead to the loss of knowledge related to such practices. In this study, we present the Eastern Africa Natural Products Database (EANPDB) containing the structural and bioactivity information of 1870 unique molecules isolated from about 300 source species from the Eastern African region. This represents the largest collection of natural products (NPs) from this geographical region, covering literature data of the period from 1962 to 2019. The computed physicochemical properties and toxicity profiles of each compound have been included. A comparative analysis of some physico-chemical properties like molecular weight, H-bond donor/acceptor, logPo/w , etc. as well scaffold diversity analysis has been carried out with other published NP databases. EANPDB was combined with the previously published Northern African Natural Products Database (NANPDB), to form a merger African Natural Products Database (ANPDB), containing ∼6500 unique molecules isolated from about 1000 source species (freely available at http://african-compounds.org). As a case study, latrunculins A and B isolated from the sponge Negombata magnifica (Podospongiidae) with previously reported antitumour activities, were identified via substructure searching as molecules to be explored as putative binders of histone deacetylases (HDACs).

The potential of anti-malarial compounds derived from African medicinal plants: a review of pharmacological evaluations from 2013 to 2019.

BACKGROUND: African Traditional Medicine (ATM) is used for the healthcare of about 80% of the rural populations of the continent of Africa. The practices of ATM make use of plant-products, which are known to contain plant-based secondary metabolites or natural products (NPs), likely to play key roles in drug discovery, particularly as lead compounds. For various reasons, including resistance of strains of Plasmodium to known anti-malarial drugs, local African populations often resort to plant-based treatments and/or a combination of this and standard anti-malarial regimens. Emphasis has been laid in this review to present the anti-malarial virtue of the most recently published phytochemicals or natural products, which have been tested by in vitro and in vivo assays. METHODS: The data was based on the current version of the African Compound Libraries, which are constantly being updated based on inputs from journal articles and student theses (M.Sc/Ph.D) from African University libraries. Emphasis was laid on data published after 2012. In order to carry out the original data collection, currently being included in the African Compounds Database, individual journal websites were queried using the country names in Africa as search terms. Over 40,000 articles "hits" were originally retrieved, then reduced to about 9000 articles. The retained articles/theses was further queried with the search terms "malaria", "malarial", "plasmodium", "plasmodial" and a combination of them, resulting in over 500 articles. Those including compounds with anti-malarial activities for which the measured activities fell within the established cut off values numbered 55, which were all cited in the review as relevant references. RESULTS AND DISCUSSION: Pure compounds derived from African medicinal plants with demonstrated anti-malarial/antiplasmodial properties with activities ranging from "very active" to "weakly active" have been discussed. The majority of the 187 natural products were terpenoids (30%), followed by flavonoids (22%), alkaloids (19%) and quinones (15%), with each of the other compound classes being less than 5% of the entire compound collection. It was also observed that most of the plant species from which the compounds were identified were of the families Rubiaceae, Meliaceae and Asphodelaceae. The review is intended to continue laying the groundwork for an African-based anti-malarial drug discovery project.

Design of selective histone deacetylase inhibitors: rethinking classical pharmacophore.

For two decades, a classical pharmacophore model comprising a zinc binding group, a linker and a cap group, has been used for the development of histone deacetylase (HDAC) inhibitors. However, some of the recently reported selective HDAC inhibitors targeting additional, usually subtype specific, cavities in the binding pocket show supplementary features which do not fit this classical pharmacophore. We, therefore, propose an extended pharmacophore model, which can describe almost all currently known HDAC inhibitors. This pharmacophore consists of six pharmacophoric features and should be helpful for the classification and design of selective HDAC inhibitors.

Identification of PKMYT1 inhibitors by screening the GSK published protein kinase inhibitor set I and II.

As a member of the Wee-kinase family protein kinase PKMYT1 is involved in G2/M checkpoint regulation of the cell cycle. Recently, a peptide microarray approach led to the identification of a small peptide; EFS247-259 as substrate of PKMYT1, which allowed for subsequent development of an activity assay. The developed activity assay was used to characterize the PKMYT1 catalyzed phosphorylation of EFS247-259. For the first time kinetic parameters for PKMYT1, namely Km, Km, ATP and vmax were determined. The optimized assay was used to screen the published protein kinase inhibitor sets (PKIS I and II), two sets of small molecule ATP-competitive kinase inhibitors reported by GlaxoSmithKline. We identified ten inhibitors, providing different scaffolds. The inhibitors were further characterized by using binding assay, activity and functional assay. In addition, docking studies were carried out in order to rationalize the observed biological activities. The derived results provide the basis for further chemical optimization of PKMYT1 inhibitors and for further analysis of PKMYT1 as target for anti-cancer therapy.

Compounds from African Medicinal Plants with Activities Against Selected Parasitic Diseases: Schistosomiasis, Trypanosomiasis and Leishmaniasis.

Parasitic diseases continue to represent a threat on a global scale, particularly among the poorest countries in the world. This is particularly because of the absence of vaccines, and in some cases, resistance against available drugs, currently being used for their treatment. In this review emphasis is laid on natural products and scaffolds from African medicinal plants (AMPs) for lead drug discovery and possible further development of drugs for the treatment of parasitic diseases. In the discussion, emphasis has been laid on alkaloids, terpenoids, quinones, flavonoids and narrower compound classes of compounds with micromolar range activities against Schistosoma, Trypanosoma and Leishmania species. In each subparagraph, emphasis is laid on the compound subclasses with most promising in vitro and/or in vivo activities of plant extracts and isolated compounds. Suggestions for future drug development from African medicinal plants have also been provided. This review covering 167 references, including 82 compounds, provides information published within two decades (1997-2017).

Identification of Bichalcones as Sirtuin Inhibitors by Virtual Screening and In Vitro Testing.

Sirtuins are nicotinamide adenine dinucleotide (NAD⁺)-dependent class III histone deacetylases, which have been linked to the pathogenesis of numerous diseases, including HIV, metabolic disorders, neurodegeneration and cancer. Docking of the virtual pan-African natural products library (p-ANAPL), followed by in vitro testing, resulted in the identification of two inhibitors of sirtuin 1, 2 and 3 (sirt1-3). Two bichalcones, known as rhuschalcone IV (8) and an analogue of rhuschalcone I (9), previously isolated from the medicinal plant Rhus pyroides, were shown to be active in the in vitro assay. The rhuschalcone I analogue (9) showed the best activity against sirt1, with an IC50 value of 40.8 µM. Based on the docking experiments, suggestions for improving the biological activities of the newly identified hit compounds have been provided.

Identification of Highly Potent Protein Kinase†C-Related Kinase†1 Inhibitors by Virtual Screening, Binding Free Energy Rescoring, and in†vitro Testing.

Despite the considerable interest in protein kinase C-related kinase 1 (PRK1) as a target in cancer research, there is still a lack of PRK1 inhibitors with suitable selectivity profiles and physicochemical properties. To identify new PRK1 inhibitors we applied a virtual screening approach, which combines ensemble docking, minimization of the protein-ligand complex, binding free energy calculations, and application of quantitative structure-activity relationship (QSAR) models for predicting in vitro activity. The developed approach was then applied in a prospective manner to screen available libraries of kinase inhibitors from Selleck and GlaxoSmithKline (GSK). Compounds that showed favorable prediction were then tested in vitro for PRK1 inhibition. Some of the hits were found to inhibit PRK1 in the low-nanomolar range. Three in vitro hits were additionally tested in a mass-spectrometry-based cellular kinase profiling assay to examine selectivity. Our findings show that nanomolar and drug-like inhibitors can be identified by the virtual screening approach presented herein. The identified inhibitors are valuable tools for gaining a better understanding of PRK1 inhibition, and the identified hits can serve as starting points for further chemical optimization.

Pharmacophore modeling and in silico toxicity assessment of potential anticancer agents from African medicinal plants.

Molecular modeling has been employed in the search for lead compounds of chemotherapy to fight cancer. In this study, pharmacophore models have been generated and validated for use in virtual screening protocols for eight known anticancer drug targets, including tyrosine kinase, protein kinase B ß, cyclin-dependent kinase, protein farnesyltransferase, human protein kinase, glycogen synthase kinase, and indoleamine 2,3-dioxygenase 1. Pharmacophore models were validated through receiver operating characteristic and Güner-Henry scoring methods, indicating that several of the models generated could be useful for the identification of potential anticancer agents from natural product databases. The validated pharmacophore models were used as three-dimensional search queries for virtual screening of the newly developed AfroCancer database (~400 compounds from African medicinal plants), along with the Naturally Occurring Plant-based Anticancer Compound-Activity-Target dataset (comprising ~1,500 published naturally occurring plant-based compounds from around the world). Additionally, an in silico assessment of toxicity of the two datasets was carried out by the use of 88 toxicity end points predicted by the Lhasa's expert knowledge-based system (Derek), showing that only an insignificant proportion of the promising anticancer agents would be likely showing high toxicity profiles. A diversity study of the two datasets, carried out using the analysis of principal components from the most important physicochemical properties often used to access drug-likeness of compound datasets, showed that the two datasets do not occupy the same chemical space.

The Chemistry and Biological Activities of Natural Products from Northern African Plant Families: From Taccaceae to Zygophyllaceae.

Traditional medicinal practices have a profound influence on the daily lives of people living in developing countries, particularly in Africa, since the populations cannot generally afford the cost of Western medicines. We have undertaken to investigate the correlation between the uses of plants in Traditional African medicine and the biological activities of the derived natural products, with the aim to validate the use of traditional medicine in Northern African communities. The literature is covered for the period 1959-2015 and part III of this review series focuses on plant families with names beginning with letters T to Z. The authors have focused on curating data from journals in natural products and phytomedicine. Within each journal home page, a query search based on country name was conducted. All articles "hits" were then verified, one at a time, that the species was harvested within the Northern African geographical regions. The current data partly constitutes the bases for the development of the Northern African natural compounds database. The review discusses 284 plant-based natural compounds from 34 species and 11 families. It was observed that the ethnobotanical uses of less than 40 % of the plant species surveyed correlated with the bioactivities of compounds identified.

Substituted 2-(2-aminopyrimidin-4-yl)pyridine-4-carboxylates as potent inhibitors of JumonjiC domain-containing histone demethylases.

BACKGROUND: Aberrant expression of iron(II)- and 2-oxoglutarate-dependent JumonjiC histone demethylases has been linked to cancer. Potent demethylase inhibitors are drug candidates and biochemical tools to elucidate the functional impact of demethylase inhibition. METHODS & RESULTS: Virtual screening identified a novel lead scaffold against JMJD2A with low-micromolar potency in vitro. Analogs were acquired from commercial sources respectively synthesized in feedback with biological testing. Optimized compounds were transformed into cell-permeable prodrugs. A cocrystal x-ray structure revealed the mode of binding of these compounds as competitive to 2-oxoglutarate and confirmed kinetic experiments. Selectivity studies revealed a preference for JMJD2A and JARID1A over JMJD3. CONCLUSION: Virtual screening and rational structural optimization led to a novel scaffold for highly potent and selective JMJD2A inhibitors.

Fluorescence-based screening assays for the NAD⁺-dependent histone deacetylase smSirt2 from Schistosoma mansoni.

Sirtuins are NAD(+)-dependent histone deacetylases (HDACs) that cleave off acetyl but also other acyl groups from the ϵ-amino group of lysines in histones and other substrate proteins. Five sirtuin isoforms are encoded in the genome of the parasitic pathogen Schistosoma mansoni. During its life cycle, S. mansoni undergoes drastic changes in phenotype that are associated with epigenetic modifications. Previous work showed strong effects of hSirt2 inhibitors on both worm life span and reproduction. Thus, we postulate smSirt2 as a new antiparasite target. We report both the optimization of a homogeneous fluorescence-based assay and the development of a new heterogeneous fluorescence-based assay to determine smSirt2 activity. The homogeneous assay uses a coumarin-labeled acetyl lysine derivative, and the heterogeneous version is using a biotinylated and fluorescence-labeled oligopeptide. Magnetic streptavidin-coated beads allow higher substrate loading per well than streptavidin-coated microtiter plates and make it possible to screen for inhibitors of either smSirt2 or its human isoform (hSirt2) for selectivity studies. We also present hits from a pilot screen with inhibitors showing an IC50 lower than 50 µM. Binding of the hits to their targets is rationalized by docking studies using a homology model of smSirt2.

Molecular modeling of potential anticancer agents from African medicinal plants.

Naturally occurring anticancer compounds represent about half of the chemotherapeutic drugs which have been put in the market against cancer until date. Computer-based or in silico virtual screening methods are often used in lead/hit discovery protocols. In this study, the "drug-likeness" of ~400 compounds from African medicinal plants that have shown in vitro and/or in vivo anticancer, cytotoxic, and antiproliferative activities has been explored. To verify potential binding to anticancer drug targets, the interactions between the compounds and 14 selected targets have been analyzed by in silico modeling. Docking and binding affinity calculations were carried out, in comparison with known anticancer agents comprising ~1,500 published naturally occurring plant-based compounds from around the world. The results reveal that African medicinal plants could represent a good starting point for the discovery of anticancer drugs. The small data set generated (named AfroCancer) has been made available for research groups working on virtual screening.

The potential of anti-malarial compounds derived from African medicinal plants, part II: a pharmacological evaluation of non-alkaloids and non-terpenoids.

Malaria is currently a public health concern in many countries in the world due to various factors which are not yet under check. Drug discovery projects targeting malaria often resort to natural sources in the search for lead compounds. A survey of the literature has led to a summary of the major findings regarding plant-derived compounds from African flora, which have shown anti-malarial/antiplasmodial activities, tested by in vitro and in vivo assays. Considerations have been given to compounds with activities ranging from "very active" to "weakly active", leading to >500 chemical structures, mainly alkaloids, terpenoids, flavonoids, coumarins, phenolics, polyacetylenes, xanthones, quinones, steroids and lignans. However, only the compounds that showed anti-malarial activity, from "very active" to "moderately active", are discussed in this review.

Virtual screening of PRK1 inhibitors: ensemble docking, rescoring using binding free energy calculation and QSAR model development.

Protein kinase C Related Kinase 1 (PRK1) has been shown to be involved in the regulation of androgen receptor signaling and has been identified as a novel potential drug target for prostate cancer therapy. Since there is no PRK1 crystal structure available to date, multiple PRK1 homology models were generated in order to address the protein flexibility. An in-house library of compounds tested on PRK1 was docked into the ATP binding site of the generated models. In most cases a correct pose of the inhibitors could be identified by ensemble docking, while there is still a challenge of finding a reasonable scoring function that is able to rank compounds according to their biological activity. We estimated the binding free energy for our data set of structurally diverse PRK1 inhibitors using the MM-PB(GB)SA and QM/MM-GBSA methods. The obtained results demonstrate that a correlation between calculated binding free energies and experimental IC50 values was found to be usually higher than using docking scores. Furthermore, the developed approach was tested on a set of diverse PRK1 inhibitors taken from literature, which resulted in a significant correlation. The developed method is computationally inexpensive and can be applied as a postdocking filter in virtual screening as well as for optimization of PRK1 inhibitors in order to prioritize compounds for further biological characterization.

Synthesis, biological evaluation and molecular docking studies of benzyloxyacetohydroxamic acids as LpxC inhibitors.

The inhibition of the UDP-3-O-[(R)-3-hydroxymyristoyl]-N-acetylglucosamine deacetylase (LpxC) represents a promising strategy to combat infections caused by multidrug-resistant Gram-negative bacteria. In order to elucidate the functional groups being important for the inhibition of LpxC, the structure of our previously reported hydroxamic acid 4 should be systematically varied. Therefore, a series of benzyloxyacetohydroxamic acids was prepared, of which the diphenylacetylene derivatives 28 (Ki=95nM) and 21 (Ki=66nM) were the most potent inhibitors of Escherichia coli LpxC. These compounds could be synthesized in a stereoselective manner employing a Sharpless asymmetric dihydroxylation and a Sonogashira coupling in the key steps. The obtained structure-activity relationships could be rationalized by molecular docking studies.

The potential of anti-malarial compounds derived from African medicinal plants, part III: an in silico evaluation of drug metabolism and pharmacokinetics profiling.

BACKGROUND: Malaria is an endemic disease affecting many countries in Tropical regions. In the search for compound hits for the design and/or development of new drugs against the disease, many research teams have resorted to African medicinal plants in order to identify lead compounds. Three-dimensional molecular models were generated for anti-malarial compounds of African origin (from 'weakly' active to 'highly' active), which were identified from literature sources. Selected computed molecular descriptors related to absorption, distribution, metabolism, excretion and toxicity (ADMET) of the phytochemicals have been analysed and compared with those of known drugs in order to access the 'drug-likeness' of these compounds. RESULTS: In the present study, more than 500 anti-malarial compounds identified from 131 distinct medicinal plant species belonging to 44 plant families from the African flora have been considered. On the basis of Lipinski's 'Rule of Five', about 70% of the compounds were predicted to be orally bioavailable, while on the basis of Jorgensen's 'Rule of Three', a corresponding >80% were compliant. An overall drug-likeness parameter indicated that approximately 55% of the compounds could be potential leads for the development of drugs. CONCLUSIONS: From the above analyses, it could be estimated that >50% of the compounds exhibiting anti-plasmodial/anti-malarial activities, derived from the African flora, could be starting points for drug discovery against malaria. The 3D models of the compounds have been included as an accompanying file and could be employed in virtual screening.

The potential of anti-malarial compounds derived from African medicinal plants, part I: a pharmacological evaluation of alkaloids and terpenoids.

Traditional medicine caters for about 80% of the health care needs of many rural populations around the world, especially in developing countries. In addition, plant-derived compounds have played key roles in drug discovery. Malaria is currently a public health concern in many countries in the world due to factors such as chemotherapy faced by resistance, poor hygienic conditions, poorly managed vector control programmes and no approved vaccines. In this review, an attempt has been made to assess the value of African medicinal plants for drug discovery by discussing the anti-malarial virtue of the derived phytochemicals that have been tested by in vitro and in vivo assays. This survey was focused on pure compounds derived from African flora which have exhibited anti-malarial properties with activities ranging from "very active" to "weakly active". However, only the compounds which showed anti-malarial activities from "very active" to "moderately active" are discussed in this review. The activity of 278 compounds, mainly alkaloids, terpenoids, flavonoids, coumarines, phenolics, polyacetylenes, xanthones, quinones, steroids, and lignans have been discussed. The first part of this review series covers the activity of 171 compounds belonging to the alkaloid and terpenoid classes. Data available in the literature indicated that African flora hold an enormous potential for the development of phytomedicines for malaria.

AfroDb: a select highly potent and diverse natural product library from African medicinal plants.

Computer-aided drug design (CADD) often involves virtual screening (VS) of large compound datasets and the availability of such is vital for drug discovery protocols. We assess the bioactivity and "drug-likeness" of a relatively small but structurally diverse dataset (containing >1,000 compounds) from African medicinal plants, which have been tested and proven a wide range of biological activities. The geographical regions of collection of the medicinal plants cover the entire continent of Africa, based on data from literature sources and information from traditional healers. For each isolated compound, the three dimensional (3D) structure has been used to calculate physico-chemical properties used in the prediction of oral bioavailability on the basis of Lipinski's "Rule of Five". A comparative analysis has been carried out with the "drug-like", "lead-like", and "fragment-like" subsets, as well as with the Dictionary of Natural Products. A diversity analysis has been carried out in comparison with the ChemBridge diverse database. Furthermore, descriptors related to absorption, distribution, metabolism, excretion and toxicity (ADMET) have been used to predict the pharmacokinetic profile of the compounds within the dataset. Our results prove that drug discovery, beginning with natural products from the African flora, could be highly promising. The 3D structures are available and could be useful for virtual screening and natural product lead generation programs.