SANCDB: an update on South African natural compounds and their readily available analogs

Artificial intelligence and cheminformatics tools: a contribution to the drug development and chemical science

In the ever-evolving field of drug discovery, the integration of Artificial Intelligence (AI) and Machine Learning (ML) with cheminformatics has proven to be a powerful combination. Cheminformatics, which combines the principles of computer science and chemistry, is used to extract chemical information and search compound databases, while the application of AI and ML allows for the identification of potential hit compounds, optimization of synthesis routes, and prediction of drug efficacy and toxicity. This collaborative approach has led to the discovery, preclinical evaluations and approval of over 70 drugs in recent years. To aid researchers in the pursuit of new drugs, this article presents a comprehensive list of databases, datasets, predictive and generative models, scoring functions and web platforms that have been launched between 2021 and 2022. These resources provide a wealth of information and tools for computer-assisted drug development, and are a valuable asset for those working in the field of cheminformatics. Overall, the integration of AI, ML and cheminformatics has greatly advanced the drug discovery process and continues to hold great potential for the future. As new resources and technologies become available, we can expect to see even more groundbreaking discoveries and advancements in these fields.Communicated by Ramaswamy H. Sarma.

A critical assessment of bioactive compounds databases

Compound databases (DBs) are essential tools for drug discovery. The number of DBs in public domain is increasing, so it is important to analyze these DBs. In this article, the main characteristics of 64 DBs will be presented. The methodological strategy used was a literature search. To analyze the characteristics obtained in the review, the DBs were categorized into two subsections: Open Access and Commercial DBs. Open access includes generalist DBs (containing compounds of diverse origins), DBs with specific applicability, DBs exclusive to natural products and those containing compounds with specific pharmacological action. The literature review showed that there are challenges to making these repositories available, such as standardizing information curation practices and funding to maintain and sustain them.

Natural Perylenequinone Compounds as Potent Inhibitors of Schistosoma mansoni Glutathione S-Transferase

The existing treatment strategy for Schistosomiasis centers on praziquantel, a single drug, but its effectiveness is limited due to resistance and lack of preventive benefits. Thus, there is an urgent need for novel antischistosomal agents. Schistosoma glutathione S-transferase (GST) is an essential parasite enzyme, with a high potential for targeted drug discovery. In this study, we conducted a screening of compounds possessing antihelminth properties, focusing on their interaction with the Schistosoma mansoni glutathione S-transferase (SmGST) protein. We demonstrated the unique nature of SmGST in comparison to human GST. Evolutionary analysis indicated its close relationship with other parasitic worms, setting it apart from free-living worms such as C. elegans. Through an assessment of binding pockets and subsequent protein-ligand docking, we identified Scutiaquinone A and Scutiaquinone B, both naturally derived Perylenequinones, as robust binders to SmGST. These compounds have exhibited effectiveness against similar parasites and offer promising potential as antischistosomal agents.

Investigating the role of Sterol C24-Methyl transferase mutation on drug resistance in leishmaniasis and identifying potential inhibitors

Leishmaniasis is a fatal disease caused by the leishmania parasite. For the survival of the leishmania parasite, Sterol C24-Methyl Transferase (SMT) is essential which is an enzyme of the ergosterol pathway. SMT protein mutation is responsible for Amphotericin-B drug resistance in Leishmania, which is the main treatment for visceral leishmaniasis. Amphotericin-B resistance is caused by three mutated residues V131I, V321I and F72C. The underlying mechanisms and structural changes in SMT enzymes responsible for resistance due to mutation are still not well understood. In the current study, the potential mechanism of resistance due to these mutations and the structure variation of wild and mutant SMT proteins were investigated through molecular dynamics simulations and molecular docking analysis. The results showed that AmB established strong bonding interaction with wild SMT as compare to mutants SMT. The binding energy calculation showed that binding energy of AmB with mutants SMT increases as compare to the wild SMT. Further structural based virtual screening was carried out to design potential inhibitors for the mutant SMT. On the basis of structural-based virtual screening four inhibitors (SANC01057, SANC00882, SANC00414, SANC01047) were computationally identified as potential mutant SMT (F72C) inhibitors. This work provides valuable information for improved management of drug resistant Leishmaniasis.Communicated by Ramaswamy H. Sarma.

Identification of Drug Targets and Their Inhibitors in Yersinia pestis Strain 91001 through Subtractive Genomics, Machine Learning, and MD Simulation Approaches

Yersinia pestis, the causative agent of plague, is a Gram-negative bacterium. If the plague is not properly treated it can cause rapid death of the host. Bubonic, pneumonic, and septicemic are the three types of plague described. Bubonic plague can progress to septicemic plague, if not diagnosed and treated on time. The mortality rate of pneumonic and septicemic plague is quite high. The symptom-defining disease is the bubo, which is a painful lymph node swelling. Almost 50% of bubonic plague leads to sepsis and death if not treated immediately with antibiotics. The host immune response is slow as compared to other bacterial infections. Clinical isolates of Yersinia pestis revealed resistance to many antibiotics such as tetracycline, spectinomycin, kanamycin, streptomycin, minocycline, chloramphenicol, and sulfonamides. Drug discovery is a time-consuming process. It always takes ten to fifteen years to bring a single drug to the market. In this regard, in silico subtractive proteomics is an accurate, rapid, and cost-effective approach for the discovery of drug targets. An ideal drug target must be essential to the pathogen's survival and must be absent in the host. Machine learning approaches are more accurate as compared to traditional virtual screening. In this study, k-nearest neighbor (kNN) and support vector machine (SVM) were used to predict the active hits against the beta-ketoacyl-ACP synthase III drug target predicted by the subtractive genomics approach. Among the 1012 compounds of the South African Natural Products database, 11 hits were predicted as active. Further, the active hits were docked against the active site of beta-ketoacyl-ACP synthase III. Out of the total 11 active hits, the 3 lowest docking score hits that showed strong interaction with the drug target were shortlisted along with the standard drug and were simulated for 100 ns. The MD simulation revealed that all the shortlisted compounds display stable behavior and the compounds formed stable complexes with the drug target. These compounds may have the potential to inhibit the beta-ketoacyl-ACP synthase III drug target and can help to combat Yersinia pestis-related infections. The dataset and the source codes are freely available on GitHub.

Multi-dimensional structural footprint identification for the design of potential scaffolds targeting METTL3 in cancer treatment from natural compounds

CONTEXT

[Formula: see text]-adenosine-methyltransferase (METTL3) is the catalytic domain of the 'writer' proteins which is involved in the post modifications of [Formula: see text]-methyladinosine ([Formula: see text]). Though its activities are essential in many biological processes, it has been implicated in several types of cancer. Thus, drug developers and researchers are relentlessly in search of small molecule inhibitors that can ameliorate the oncogenic activities of METTL3. Currently, STM2457 is a potent, highly selective inhibitor of METTL3 but is yet to be approved.

METHODS

In this study, we employed structure-based virtual screening through consensus docking by using AutoDock Vina in PyRx interface and Glide virtual screening workflow of Schrodinger Glide. Thermodynamics via MM-PBSA calculations was further used to rank the compounds based on their total free binding energies. All atom molecular dynamics simulations were performed using AMBER 18 package. FF14SB force fields and Antechamber were used to parameterize the protein and compounds respectively. Post analysis of generated trajectories was analyzed with CPPTRAJ and PTRAJ modules incorporated in the AMBER package while Discovery studio and UCSF Chimera were used for visualization, and origin data tool used to plot all graphs.

RESULTS

Three compounds with total free binding energies higher than STM2457 were selected for extended molecular dynamics simulations. The compounds, SANCDB0370, SANCDB0867, and SANCDB1033, exhibited stability and deeper penetration into the hydrophobic core of the protein. They engaged in relatively stronger intermolecular interactions involving hydrogen bonds with resultant increase in stability, reduced flexibility, and decrease in the surface area of the protein available for solvent interactions suggesting an induced folding of the catalytic domain. Furthermore, in silico pharmacokinetics and physicochemical analysis of the compounds revealed good properties suggesting these compounds could serve as promising MEETL3 entry inhibitors upon modifications and optimizations as presented by natural compounds. Further biochemical testing and experimentations would aid in the discovery of effective inhibitors against the berserk activities of METTL3.

Identification of promising high-affinity inhibitors of SARS-CoV-2 main protease from African Natural Products Databases by Virtual Screening

With the rapid spread of the new severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the pathogen agent of COVID-19 pandemic created a serious threat to global public health, requiring the most urgent research for potential therapeutic agents. The availability of genomic data of SARS-CoV-2 and efforts to determine the protein structure of the virus facilitated the identification of potent inhibitors by using structure-based approach and bioinformatics tools. Many pharmaceuticals have been proposed for the treatment of COVID-19, although their effectiveness has not been assessed yet. However, it is important to find out new-targeted drugs to overcome the resistance concern. Several viral proteins such as proteases, polymerases or structural proteins have been considered as potential therapeutic targets. But the virus target must be essential for host invasion match some drugability criterion. In this Work, we selected the highly validated pharmacological target main protease Mpro and we performed high throughput virtual screening of African Natural Products Databases such as NANPDB, EANPDB, AfroDb, and SANCDB to identify the most potent inhibitors with the best pharmacological properties. In total, 8753 natural compounds were virtually screened by AutoDock vina against the main protease of SARS-CoV-2. Two hundred and five (205) compounds showed high-affinity scores (less than - 10.0 Kcal/mol), while fifty-eight (58) filtered through Lipinski's rules showed better affinity than known Mpro inhibitors (i.e., ABBV-744, Onalespib, Daunorubicin, Alpha-ketoamide, Perampanel, Carprefen, Celecoxib, Alprazolam, Trovafloxacin, Sarafloxacin, Ethyl biscoumacetate…). Those promising compounds could be considered for further investigations toward the developpement of SARS-CoV-2 drug development.

Discovery of Isojacareubin as a covalent inhibitor of SARS-CoV-2 main protease using structural and experimental approaches

The ongoing pandemic with the emergence of immune evasion potential and, particularly, the current omicron subvariants intensified the situation further. Although vaccines are available, the immune evasion capabilities of the recent variants demand further efficient therapeutic choices to control the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic. Hence, considering the necessity of the small molecule inhibitor, we target the main protease (3CLpro), which is an appealing target for the development of antiviral drugs against SARS-CoV-2. High-throughput molecular in silico screening of South African natural compounds database reported Isojacareubin and Glabranin as the potential inhibitors for the main protease. The calculated docking scores were reported to be -8.47 and -8.03 kcal/mol, respectively. Moreover, the structural dynamic assessment reported that Isojacareubin in complex with 3CLpro exhibit a more stable dynamic behavior than Glabranin. Inhibition assay indicated that Isojacareubin could inhibit SARS-CoV-2 3CLpro in a time- and dose-dependent manner, with half maximal inhibitory concentration values of 16.00 ± 1.35 μM (60 min incubation). Next, the covalent binding sites of Isojacareubin on SARS-CoV-2 3CLpro was identified by biomass spectrometry, which reported that Isojacareubin can covalently bind to thiols or Cysteine through Michael addition. To evaluate the inactivation potency of Isojacareubin, the inactivation kinetics was further investigated. The inactivation kinetic curves were plotted according to various concentrations with gradient-ascending incubation times. The KI value of Isojacareubin was determined as 30.71 μM, whereas the Kinact value was calculated as 0.054 min-1 . These results suggest that Isojacareubin is a covalent inhibitor of SARS-CoV-2 3CLpro .

Ultra-Large-Scale Screening of Natural Compounds and Free Energy Calculations Revealed Potential Inhibitors for the Receptor-Binding Domain (RBD) of SARS-CoV-2

The emergence of immune-evading variants of SARS-CoV-2 further aggravated the ongoing pandemic. Despite the deployments of various vaccines, the acquired mutations are capable of escaping both natural and vaccine-induced immune responses. Therefore, further investigation is needed to design a decisive pharmacological treatment that could efficiently block the entry of this virus into cells. Hence, the current study used structure-based methods to target the RBD of the recombinant variant (Deltacron) of SARS-CoV-2, which was used as a model variant. From the virtual drug screenings of various databases, a total of four hits were identified as potential lead molecules. Key residues were blocked by these molecules with favorable structural dynamic features. The binding free energies further validated the potentials of these molecules. The TBE for MNP was calculated to be -32.86 ± 0.10 kcal/mol, for SANC00222 the TBE was -23.41 ± 0.15 kcal/mol, for Liriodenine the TBE was -34.29 ± 0.07 kcal/mol, while for Carviolin the TBE was calculated to be -27.67 ± 0.12 kcal/mol. Moreover, each complex demonstrated distinct internal motion and a free energy profile, indicating a different strategy for the interaction with and inhibition of the RBD. In conclusion, the current study demands further in vivo and in vitro validation for the possible usage of these compounds as potential drugs against SARS-CoV-2 and its variants.

In-silico investigation of a novel inhibitors against the antibiotic-resistant Neisseria gonorrhoeae bacteria

Antibiotics are drugs that are used to treat or prevent bacterial infections. They work by either killing or stopping bacteria from spreading. Nevertheless, it appeared in the last decade, Antibiotic-resistant bacteria are bacteria resistant to antibiotics and cannot be controlled or killed by them. In the presence of an antibiotic, they can live and even reproduce. The Neisseria gonorrhoeae bacteria is appearing to be a multidrug-resistant pathogen. Many factors contribute to antibiotic resistance, including unfettered access to antimicrobials, incorrect drug selection, misuse, and low-quality antibiotics. Here, we investigated in-silico docking screening and analysis for ten natural marine fungus extracted compounds. The resulted data were examined for the best binding affinity, toxicity, and chemical interactions. The most superior compound was elipyrone A with six hydrogen bonds, −8.5 of binding affinity, and preferable results in the SWISS-ADME examination. It is well known that “Declining corporate investment and a lack of innovation in the development of new antibiotics are weakening efforts to battle drug-resistant illnesses,” according to the World Health Organization (WHO). So, we extended our effort to predict a new natural compound to overcome the resistance of this bacteria.

Anti-trypanosomatid structure-based drug design - lessons learned from targeting the folate pathway

INTRODUCTION

Trypanosomatidic parasitic infections of humans and animals caused by Trypanosoma brucei, Trypanosoma cruzi, and Leishmania species pose a significant health and economic burden in developing countries. There are few effective and accessible treatments for these diseases, and the existing therapies suffer from problems such as parasite resistance and side effects. Structure-based drug design (SBDD) is one of the strategies that has been applied to discover new compounds targeting trypanosomatid-borne diseases.

AREAS COVERED

We review the current literature (mostly over the last 5 years, searched in PubMed database on Nov 11^th 2021) on the application of structure-based drug design approaches to identify new anti-trypanosomatidic compounds that interfere with a validated target biochemical pathway, the trypanosomatid folate pathway.

EXPERT OPINION

The application of structure-based drug design approaches to perturb the trypanosomatid folate pathway has successfully provided many new inhibitors with good selectivity profiles, most of which are natural products or their derivatives or have scaffolds of known drugs. However, the inhibitory effect against the target protein(s) often does not translate to anti-parasitic activity. Further progress is hampered by our incomplete understanding of parasite biology and biochemistry, which is necessary to complement SBDD in a multiparameter optimization approach to discovering selective anti-parasitic drugs.

Learning about allosteric drugs and ways to design them

While the accelerating quest for precision medicine requires new individually targeting and selective drugs, and the ability to work with so-called undruggable targets, the realm of allosteric drugs meeting this need remains largely uncharted. Generalizing the observations on two major drug targets with widely observed inherent allostery, GPCRs and kinases, we describe and discuss basic allosteric modes of action that are universally applicable in all types of structures and functions. Using examples of Class A GPCRs and CMGC protein kinases, we show how Allosteric Signalling and Probing Fingerprints can be used to identify potential allosteric sites and reveal effector-leads that may serve as a starting point for the development of allosteric drugs targeting these regulatory sites. A set of distinct characteristics of allosteric ligands was established, which highlights the versatility of their design and make them advantageous before their orthosteric counterparts in personalized medicine. We argue that rational design of allosteric drugs should begin with the search for latent sites or design of non-natural binding sites followed by fragment-based design of allosteric ligands and by the mutual adjustment of the site-ligand pair in order to achieve required effects. On the basis of the perturbative nature and reversibility of allosteric communication, we propose a generic protocol for computational design of allosteric effectors, enabling also the allosteric tuning of biologics, in obtaining allosteric control over protein functions.

Targeting the RBD of Omicron Variant (B.1.1.529) with Medicinal Phytocompounds to Abrogate the Binding of Spike Glycoprotein with the hACE2 Using Computational Molecular Search and Simulation Approach

Simple Summary

The current study based on virtual drugs screening and simulations identified novel drugs to target the RBD of the spike protein from Omicron variant of SARS-CoV-2. Using molecular modeling tools to search for a good binding drugs we identified SANC00944, SANC01032, SANC00992, and SANC00317 from South African natural compounds database as potential inhibitor of the Spike-ACE2 complex. In sum, this study will help in the design and discovery of novel drug therapeutics, which may be used against the emerging Omicron variant of SARS-CoV-2.

Abstract

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus continues to inflict chaos globally. The emergence of a novel Omicron variant (B.1.1.529) in South Africa harbors 30 mutations in the spike protein. The variant is distinguished from other variants of concern (VOCs) with an increased (15) number of mutations in the receptor-binding domain (RBD) and suggests higher chances of causing reinfections. Initial reports also claimed that this variant escapes all the neutralizing antibodies, thus demanding a novel strategy against it. Thus, in this study, we performed a computational molecular screening against the RBD of the Omicron (B.1.1.529) variant and assessed the binding affinity of potent drugs against the RBD. The multi-steps screening of the South African Natural Compounds Database (SANCDB) revealed four medicinal compounds as excellent (potential) anti-viral agents against the Omicron variant, namely SANC00944, SANC01032, SANC00992, and SANC00317. The simulation analysis of these compounds in complex with the RBD demonstrated stable dynamics and structural compactness. Moreover, the residual flexibility analysis revealed that the flexibility of three loops required for interaction with hACE2 has been reduced by the binding of these drugs. The post-simulation validation of these compounds such as binding free energy, in silico bioactivity, and dissociation constant prediction validated the anti-viral potency of these compounds. The total binding free energy (TBFE) for the SANC01032–RBD complex was reported to be −46.54 kcal/mol; for the SANC01032–RBD complex, the TBFE was −41.88 kcal/mol; for the SANC00992–RBD complex the TBFE was −29.05 kcal/mol, while for the SANC00317–RBD complex the TBFE was −31.03 kcal/mol. The results showed the inhibition potential of these compounds by targeting the RBD. In conclusion, this study will help in the design and discovery of novel drug therapeutics, which may be used against the emerging Omicron variant of SARS-CoV-2.

Rescuing the Host Immune System by Targeting the Immune Evasion Complex ORF8-IRF3 in SARS-CoV-2 Infection with Natural Products Using Molecular Modeling Approaches

The perennial emergence of SARS-CoV-2 and its new variants causing upper respiratory complexities since December 2019 has aggravated the pandemic situation around the world. SARS-CoV-2 encodes several proteins among which ORF8 is a novel factor that is unique to SARS-CoV-2 only and is reported to help the virus in disease severity and immune evasion. ORF8-IRF3 complex induces endoplasmic reticulum stress, thus helps in the evasion of immune response. Consequently, targeting the ORF8-IRF3 complex is considered as a prime target for the discovery of novel drugs against SARS-CoV-2. In this regard, computational methods are of great interest to fast track the identification and development of novel drugs. Virtual screening of South African Natural Compounds Database (SANCDB), followed by docking and molecular dynamics (MD) simulation analysis, were performed to determine novel natural compounds. Computational molecular search and rescoring of the SANCDB database followed by induced-fit docking (IFD) protocol identified Quercetin 3-O-(6″-galloyl)-beta-D-galactopyranoside (SANC00850), Tribuloside (SANC01050), and Rutin (SANC00867) are the best scoring compounds. Structural-dynamic properties assessment revealed that these three compounds have stable dynamics, compactness, and a higher number of hydrogen bonds. For validation, we used MM/GBSA, in silico bioactivity estimation and dissociation constant (KD) approaches, which revealed that these compounds are the more potent inhibitors of the ORF8-IRF3 complex and would rescue the host immune system potentially. These compounds need further in vitro and in vivo validations to be used as therapeutics against SARS-CoV-2 to rescue the host immune system during COVID-19 infection.

Novel dynamic residue network analysis approaches to study allosteric modulation: SARS-CoV-2 Mpro and its evolutionary mutations as a case study

The rational search for allosteric modulators and the allosteric mechanisms of these modulators in the presence of mutations is a relatively unexplored field. Here, we established novel in silico approaches and applied them to SARS-CoV-2 main protease (Mpro) as a case study. First, we identified six potential allosteric modulators. Then, we focused on understanding the allosteric effects of these modulators on each of its protomers. We introduced a new combinatorial approach and dynamic residue network (DRN) analysis algorithms to examine patterns of change and conservation of critical nodes, according to five independent criteria of network centrality. We observed highly conserved network hubs for each averaged DRN metric on the basis of their existence in both protomers in the absence and presence of all ligands (persistent hubs). We also detected ligand specific signal changes. Using eigencentrality (EC) persistent hubs and ligand introduced hubs we identified a residue communication path connecting the allosteric binding site to the catalytic site. Finally, we examined the effects of the mutations on the behavior of the protein in the presence of selected potential allosteric modulators and investigated the ligand stability. One crucial outcome was to show that EC centrality hubs form an allosteric communication path between the allosteric ligand binding site to the active site going through the interface residues of domains I and II; and this path was either weakened or lost in the presence of some of the mutations. Overall, the results revealed crucial aspects that need to be considered in rational computational drug discovery.

Computational Applications in Secondary Metabolite Discovery (CAiSMD): an online workshop

We report the major conclusions of the online open-access workshop "Computational Applications in Secondary Metabolite Discovery (CAiSMD)" that took place from 08 to 10 March 2021. Invited speakers from academia and industry and about 200 registered participants from five continents (Africa, Asia, Europe, South America, and North America) took part in the workshop. The workshop highlighted the potential applications of computational methodologies in the search for secondary metabolites (SMs) or natural products (NPs) as potential drugs and drug leads. During 3 days, the participants of this online workshop received an overview of modern computer-based approaches for exploring NP discovery in the "omics" age. The invited experts gave keynote lectures, trained participants in hands-on sessions, and held round table discussions. This was followed by oral presentations with much interaction between the speakers and the audience. Selected applicants (early-career scientists) were offered the opportunity to give oral presentations (15 min) and present posters in the form of flash presentations (5 min) upon submission of an abstract. The final program available on the workshop website ( https://caismd.indiayouth.info/ ) comprised of 4 keynote lectures (KLs), 12 oral presentations (OPs), 2 round table discussions (RTDs), and 5 hands-on sessions (HSs). This meeting report also references internet resources for computational biology in the area of secondary metabolites that are of use outside of the workshop areas and will constitute a long-term valuable source for the community. The workshop concluded with an online survey form to be completed by speakers and participants for the goal of improving any subsequent editions.