Search

Approximately 60 % of men globally over the age of 50 experience a diminished quality of life as a result of benign prostatic hyperplasia (BPH). Consequently, while the national direct cost for managing BPH is hovering around $4 billion in the USA, the individual direct cost for medications is estimated at $1536 and $425 annually for the USA and Ghana, respectively. Due to the chemotherapeutic drawbacks in BPH treatment, a concerted effort is urgently needed to find new chemotypes with novel mechanism of action to combat the disease. This research aims to employ computational techniques to identify new compounds from ethno-pharmacological plants from Ghana for treating BPH. Altogether, 250 natural products from Croton membranaceus, Heliotropium angiospermum, Annona muricata, Serenoa repens, Cissus quadrangularis, Tribulus terrestrials, Vernonia amagdalina and Momordica foetida from Ghanaian origin were virtually screened against 5alpha reductase 2 (5aR2), a primary drug target implicated in the parthenogenesis of BPH. Out of the 254 compounds docked, three compounds (A, B and C) showed a binding energy lower than testosterone (-10.4 kcal/mol), the main substrate of the target and comparative to finasteride (-11.7 kcal/mol), one of the drugs currently used for treating the disease. Molecular docking and dynamics simulations studies revealed Leu224 to be critical for ligand binding and complex stability. The pharmacological and physicochemical profiles show the compounds to possess good pharmacodynamics with negligible toxicities. While the compounds A, C, G, H, I and J were, respectively, predicted to treat prostate disorders with Pa (0.909, 0.221, 0.743, 0.729, 0.978, and 0.443) and Pi (0.003, 0.132, 0.005,0.005, 0.002, 0.016), that of A, B, D, G, H, I, and J were found to be associated with the inhibition of 5alpha reductase with Pa (0.840, 0.102, 0.097, 0.377, 0.123, 0.909, and 0.201) and Pi (0.001, 0.075, 0.089, 0.003, 0.036, 0.001, and 0.004), respectively. These lead-like compounds possessing plausible chemotherapeutic moieties with putative 5aR2 inhibitory potential need further experimental validation via in vitro studies.

This review aims to highlight the structures of ADAR proteins that have been crucial in the discernment of their functions and are relevant to future therapeutic development. ADAR proteins can correct or diversify genetic information, underscoring their pivotal contribution to protein diversity and the sophistication of neuronal networks. ADAR proteins have numerous functions in RNA editing independent roles and through the mechanisms of A-I RNA editing that continue to be revealed. Provided is a detailed examination of the ADAR family members—ADAR1, ADAR2, and ADAR3—each characterized by distinct isoforms that offer both structural diversity and functional variability, significantly affecting RNA editing mechanisms and exhibiting tissue-specific regulatory patterns, highlighting their shared features, such as double-stranded RNA binding domains (dsRBD) and a catalytic deaminase domain (CDD). Moreover, it explores ADARs’ extensive roles in immunity, RNA interference, and disease modulation, demonstrating their ambivalent nature in both the advancement and inhibition of diseases. Through this comprehensive analysis, the review seeks to underline the potential of targeting ADAR proteins in therapeutic strategies, urging continued investigation into their biological mechanisms and health implications.

Monkeypox (now Mpox), a zoonotic disease caused by the monkeypox virus (MPXV) is an emerging threat to global health. In the time span of only six months, from May to October 2022, the number of MPXV cases breached 80,000 and many of the outbreaks occurred in locations that had never previously reported MPXV. Currently there are no FDA-approved MPXV-specific vaccines or treatments, therefore, finding drugs to combat MPXV is of utmost importance. The A42R profilin-like protein of the MPXV is involved in cell development and motility making it a critical drug target. A42R protein is highly conserved across orthopoxviruses, thus A42R inhibitors may work for other family members. This study sought to identify potential A42R inhibitors for MPXV treatment using computational approaches. The energy minimized 3D structure of the A42R profilin-like protein (PDB ID: 4QWO) underwent virtual screening using a library of 36,366 compounds from Traditional Chinese Medicine (TCM), AfroDb, and PubChem databases as well as known inhibitor tecovirimat via AutoDock Vina. A total of seven compounds comprising PubChem CID: 11371962, ZINC000000899909, ZINC000001632866, ZINC000015151344, ZINC000013378519, ZINC000000086470, and ZINC000095486204, predicted to have favorable binding were shortlisted. Molecular docking suggested that all seven proposed compounds have higher binding affinities to A42R (–7.2 to –8.3 kcal/mol) than tecovirimat (–6.7 kcal/mol). This was corroborated by MM/PBSA calculations, with tecovirimat demonstrating the highest binding free energy of –68.694 kJ/mol (lowest binding affinity) compared to the seven shortlisted compounds that ranged from –73.252 to –97.140 kJ/mol. Furthermore, the 7 compounds in complex with A42R demonstrated higher stability than the A42R-tecovirimat complex when subjected to 100 ns molecular dynamics simulations. The protein-ligand interaction maps generated using LigPlot+ suggested that residues Met1, Glu3, Trp4, Ile7, Arg127, Val128, Thr131, and Asn133 are important for binding. These seven compounds were adequately profiled to be potential antivirals via PASS predictions and structural similarity searches. All seven potential lead compounds were scored Pa > Pi for antiviral activity while ZINC000001632866 and ZINC000015151344 were predicted as poxvirus inhibitors with Pa values of 0.315 and 0.215, and Pi values of 0.052 and 0.136, respectively. Further experimental validations of the identified lead compounds are required to corroborate their predicted activity. These seven identified compounds represent solid footing for development of antivirals against MPXV and other orthopoxviruses.

Cervical cancer is the fourth most diagnosed cancer and the fourth leading cause of cancer death in women globally. Its onset and progression have been attributed to high-risk human papillomavirus (HPV) types, especially 16 and 18, while the Epstein–Barr virus (EBV) is believed to also significantly contribute to cervical cancer growth. The E6 protein associated with high-risk HPV strains, such as HPV16 and HPV18, is known for its role in promoting cervical cancer and other anogenital cancers. E6 proteins contribute to the malignant transformation of infected cells by targeting and degrading tumor suppressor proteins, especially p53. On the other hand, EBV nuclear antigen 1 (EBNA1) plays a crucial role in the maintenance and replication of the EBV genome in infected cells. EBNA1 is believed to increase HPV E6 and E7 levels, as well as c-MYC, and BIRC5 cellular genes in the HeLa cell line, implying that HPV/EBV co-infection accelerates cervical cancer onset and growth. Thus, the E6 and EBNA1 antigens of HPV and EBV, respectively, are attractive targets for cervical cancer immunotherapy. This study, therefore, virtually screened for potential drug candidates with good binding affinity to all three oncoviral proteins, HPV16 E6, HPV18 E6, and EBNA1. The compounds were further subjected to ADMET profiling, biological activity predictions, molecular dynamics (MD) simulations, and molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) calculations. A total of six compounds comprising ZINC000013380012, ZINC000070454124, ZINC000014588133, ZINC000085568136, ZINC000095909247, and ZINC000085597263 demonstrated very strong affinity (≤−60 kJ/mol) to the three oncoviral proteins (EBNA1, HPV16 E6, and HPV18 E6) after being subjected to docking, MD, and MM/PBSA. These compounds demonstrated relatively stronger binding than the controls used, inhibitors of EBNA1 (VK-1727) and HPV E6 (baicalein and gossypetin). Biological activity predictions also corroborated their antineoplastic, p53-enhancing, Pin1 inhibitory, and JAK2 inhibitory activities. Further experimental testing is required to validate the ability of the shortlisted compounds to silence the insidious effects of HPV E6 and EBNA1 proteins in cervical cancers.

The prognosis of mixed-lineage leukemia (MLL) has remained a significant health concern, especially for infants. The minimal treatments available for this aggressive type of leukemia has been an ongoing problem. Chromosomal translocations of the KMT2A gene are known as MLL, which expresses MLL fusion proteins. A protein called menin is an important oncogenic cofactor for these MLL fusion proteins, thus providing a new avenue for treatments against this subset of acute leukemias. In this study, we report results using the structure-based drug design (SBDD) approach to discover potential novel MLL-mediated leukemia inhibitors from natural products against menin. The three-dimensional (3D) protein model was derived from Protein Databank (Protein ID: 4GQ4), and EasyModeller 4.0 and I-TASSER were used to fix missing residues during rebuilding. Out of the ten protein models generated (five from EasyModeller and I-TASSER each), one model was selected. The selected model demonstrated the most reasonable quality and had 75.5% of residues in the most favored regions, 18.3% of residues in additionally allowed regions, 3.3% of residues in generously allowed regions, and 2.9% of residues in disallowed regions. A ligand library containing 25,131 ligands from a Chinese database was virtually screened using AutoDock Vina, in addition to three known menin inhibitors. The top 10 compounds including ZINC000103526876, ZINC000095913861, ZINC000095912705, ZINC000085530497, ZINC000095912718, ZINC000070451048, ZINC000085530488, ZINC000095912706, ZINC000103580868, and ZINC000103584057 had binding energies of −11.0, −10.7, −10.6, −10.2, −10.2, −9.9, −9.9, −9.9, −9.9, and −9.9 kcal/mol, respectively. To confirm the stability of the menin–ligand complexes and the binding mechanisms, molecular dynamics simulations including molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) computations were performed. The amino acid residues that were found to be potentially crucial in ligand binding included Phe243, Met283, Cys246, Tyr281, Ala247, Ser160, Asn287, Asp185, Ser183, Tyr328, Asn249, His186, Leu182, Ile248, and Pro250. MI-2-2 and PubChem CIDs 71777742 and 36294 were shown to possess anti-menin properties; thus, this justifies a need to experimentally determine the activity of the identified compounds. The compounds identified herein were found to have good pharmacological profiles and had negligible toxicity. Additionally, these compounds were predicted as antileukemic, antineoplastic, chemopreventive, and apoptotic agents. The 10 natural compounds can be further explored as potential novel agents for the effective treatment of MLL-mediated leukemia.

The mortalities and morbidities of leishmaniasis are high and the disease is under reported globally. The absence of vaccines coupled with chemotherapeutic challenges including chemoresistance, scarcity and toxicity have made the fight against leishmaniasis an arduous one. Furthermore, the treatment options currently available for leishmaniasis are long and sometimes require hospitalization. There is therefore the need to explore novel pathways to identify new compounds with alternative mechanisms of action. A pharmacophore-based screening was employed in identifying new potential inhibitors with unique scaffolds targeting Leishmania donovani sterol methyltransferase (LdSMT), a key enzyme for ergosterol biosynthesis. To accomplish this, 22,26-azasterol, a known inhibitor of this target and five other derivatives with IC50 less than 10 μM were used to generate a robust 3D pharmacophore model via LigandScout with a score of 0.9144. The validated model was used as a query to screen a library of 69034 natural products obtained from the InterBioScreen Limited. Compounds with pharmacophore fit scores above 50 were docked against the modelled structure of LdSMT. Altogether, ten molecules with binding energies between −7 and −11 kcal/mol were identified as potential bioactive molecules. The molecular dynamics simulation and molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) calculations reinforced the results from the docking studies suggesting the selected hits bind effectively at the active sites of the target protein. The compounds were observed to bind in the S-adenosine-L-homocysteine binding pocket of the modelled LdSMT with Trp208 and Val330 predicted as key residues critical for ligand binding. Prediction of biological activity with probability of activity (Pa) greater than probability of inactivity (Pi) revealed that seven compounds (STOCKIN-54848, STOCKIN-89115, STOCKIN-68720, STOCKIN-44724, STOCKIN-76694, STOCKIN-47277 and STOCKIN-95708) possessed antileishmanial properties. STOCKIN-89115, STOCKIN-68720, STOCKIN-44724, and STOCKIN-47277 were predicted to be membrane permeability inhibitors, while all ten hit compounds possessed antineoplastic activity. The compounds have the propensity of disrupting ergosterol biosynthesis leading to the suppression of growth in Leishmania donovani. The compounds were predicted to have good absorption, distribution, metabolism, excretion and toxicity profiles, hence their potential antileishmanial activity can be exploited upon experimental corroboration.

Adenosine deaminase acting on RNA 2 (ADAR2) is an important enzyme involved in RNA editing processes, particularly in the conversion of adenosine to inosine in RNA molecules. Dysregulation of ADAR2 activity has been implicated in various diseases, including neurological disorders (including schizophrenia), inflammatory disorders, viral infections, and cancers. Therefore, targeting ADAR2 with small molecules presents a promising therapeutic strategy for modulating RNA editing and potentially treating associated pathologies. However, there are limited compounds that effectively inhibit ADAR2 reactions. This study therefore employed computational approaches to virtually screen natural compounds from the traditional Chinese medicine (TCM) library. The shortlisted compounds demonstrated a stronger binding affinity to the ADAR2 (2CR with binding energies ranging from −7.8 to −12.9 kcal/mol. Further subjecting the top ADAR2–ligand complexes to molecular dynamics simulations and molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) calculations revealed that five potential hit compounds comprising ZINC000014637370, ZINC000085593577, ZINC000042890265, ZINC000039183320, and ZINC000101100339 had favorable binding free energies of −174.911, −137.369, −117.236, −67.023, and −64.913 kJ/mol, respectively, with the human ADAR2 protein. Residues Lys350, Cys377, Glu396, Cys451, Arg455, Ser486, Gln488, and Arg510 were also predicted to be crucial in ligand recognition and binding. This finding will provide valuable insights into the molecular interactions between ADAR2 and small molecules, aiding in the design of future ADAR2 inhibitors with potential therapeutic applications. The potential lead compounds were also profiled to have insignificant toxicities. A structural similarity search via DrugBank revealed that ZINC000039183320 and ZINC000014637370 were similar to naringin and naringenin, which are known adenosine deaminase (ADA) inhibitors. These potential novel ADAR2 inhibitors identified herein may be beneficial in treating several neurological disorders, cancers, viral infections, and inflammatory disorders caused by ADAR2 after experimental validation.

The therapeutic challenges pertaining to leishmaniasis due to reported chemoresistance and toxicity necessitate the need to explore novel pathways to identify plausible inhibitory molecules. Leishmania donovani 24-sterol methyltransferase (LdSMT) is vital for the synthesis of ergosterols, the main constituents of Leishmania cellular membranes. So far, mammals have not been shown to possess SMT or ergosterols, making the pathway a prime candidate for drug discovery. The structural model of LdSMT was elucidated using homology modeling to identify potential novel 24-SMT inhibitors via virtual screening, scaffold hopping, and de-novo fragment-based design. Altogether, six potential novel inhibitors were identified with binding energies ranging from −7.0 to −8.4 kcal/mol with e-LEA3D using 22,26-azasterol and S1–S4 obtained from scaffold hopping via the ChEMBL, DrugBank, PubChem, ChemSpider, and ZINC15 databases. These ligands showed comparable binding energy to 22,26-azasterol (−7.6 kcal/mol), the main inhibitor of LdSMT. Moreover, all the compounds had plausible ligand efficiency-dependent lipophilicity (LELP) scores above 3. The binding mechanism identified Tyr92 to be critical for binding, and this was corroborated via molecular dynamics simulations and molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) calculations. The ligand A1 was predicted to possess antileishmanial properties with a probability of activity (Pa) of 0.362 and a probability of inactivity (Pi) of 0.066, while A5 and A6 possessed dermatological properties with Pa values of 0.205 and 0.249 and Pi values of 0.162 and 0.120, respectively. Structural similarity search via DrugBank identified vabicaserin, daledalin, zanapezil, imipramine, and cefradine with antileishmanial properties suggesting that the de-novo compounds could be explored as potential antileishmanial agents.

Altered RNA editing has been linked to several neurodevelopmental disorders, including autism spectrum disorder (ASD) and intellectual disability, in addition to depression, schizophrenia, some cancers, viral infections and autoimmune disorders. The human ADAR2 is a potential therapeutic target for managing these various disorders due to its crucial role in adenosine to inosine editing. This study applied consensus scoring to rank potential ADAR2 inhibitors after performing molecular docking with AutoDock Vina and Glide (Maestro), using a library of 35,161 compounds obtained from traditional Chinese medicine. A total of 47 compounds were predicted to be good binders of the human ADAR2 and had insignificant toxicity concerns. Molecular dynamics (MD) simulations, including the molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) procedure, also emphasized the binding of the shortlisted compounds. The potential compounds had plausible binding free energies ranging from −81.304 to −1068.26 kJ/mol from the MM/PBSA calculations. ZINC000085511995, a naphthoquinone had more negative binding free energy (−1068.26 kJ/mol) than inositol hexakisphosphate (IHP) [−873.873 kJ/mol], an agonist and a strong binder of ADAR2. The potential displacement of IHP by ZINC000085511995 in the IHP binding site of ADAR2 could be explored for possible deactivation of ADAR2. Bayesian-based biological activity prediction corroborates the neuropharmacological, antineoplastic and antiviral activity of the potential lead compounds. All the potential lead compounds, except ZINC000014612330 and ZINC000013462928, were predicted to be inhibitors of various deaminases. The potential lead compounds also had probability of activity (Pa) > 0.442 and probability of inactivity (Pi) < 0.116 values for treating acute neurologic disorders, except for ZINC000085996580 and ZINC000013462928. Pursuing these compounds for their anti-ADAR2 activities holds a promising future, especially against neurological disorders, some cancers and viral infections caused by RNA viruses. Molecular interaction, hydrogen bond and per-residue decomposition analyses predicted Arg400, Arg401, Lys519, Trp687, Glu689, and Lys690 as hot-spot residues in the ADAR2 IHP binding site. Most of the top compounds were observed to have naphthoquinone, indole, furanocoumarin or benzofuran moieties. Serotonin and tryptophan, which are beneficial in digestive regulation, improving sleep cycle and mood, are indole derivatives. These chemical series may have the potential to treat neurological disorders, prion diseases, some cancers, specific viral infections, metabolic disorders and eating disorders through the disruption of ADAR2 pathways. A total of nine potential lead compounds were shortlisted as plausible modulators of ADAR2.

The Ebola virus (EBOV) is still highly infectious and causes severe hemorrhagic fevers in primates. However, there are no regulatorily approved drugs against the Ebola virus disease (EVD). The highly virulent and lethal nature of EVD highlights the need to develop therapeutic agents. Viral protein 40 kDa (VP40), the most abundantly expressed protein during infection, coordinates the assembly, budding, and release of viral particles into the host cell. It also regulates viral transcription and RNA replication. This study sought to identify small molecules that could potentially inhibit the VP40 protein by targeting the N-terminal domain using an in silico approach. The statistical quality of AutoDock Vina’s capacity to discriminate between inhibitors and decoys was determined, and an area under the curve of the receiver operating characteristic (AUC-ROC) curve of 0.791 was obtained. A total of 29,519 natural-product-derived compounds from Chinese and African sources as well as 2738 approved drugs were successfully screened against VP40. Using a threshold of −8 kcal/mol, a total of 7, 11, 163, and 30 compounds from the AfroDb, Northern African Natural Products Database (NANPDB), traditional Chinese medicine (TCM), and approved drugs libraries, respectively, were obtained after molecular docking. A biological activity prediction of the lead compounds suggested their potential antiviral properties. In addition, random-forest- and support-vector-machine-based algorithms predicted the compounds to be anti-Ebola with IC50 values in the micromolar range (less than 25 μM). A total of 42 natural-product-derived compounds were identified as potential EBOV inhibitors with desirable ADMET profiles, comprising 1, 2, and 39 compounds from NANPDB (2-hydroxyseneganolide), AfroDb (ZINC000034518176 and ZINC000095485942), and TCM, respectively. A total of 23 approved drugs, including doramectin, glecaprevir, velpatasvir, ledipasvir, avermectin B1, nafarelin acetate, danoprevir, eltrombopag, lanatoside C, and glycyrrhizin, among others, were also predicted to have potential anti-EBOV activity and can be further explored so that they may be repurposed for EVD treatment. Molecular dynamics simulations coupled with molecular mechanics Poisson–Boltzmann surface area calculations corroborated the stability and good binding affinities of the complexes (−46.97 to −118.9 kJ/mol). The potential lead compounds may have the potential to be developed as anti-EBOV drugs after experimental testing.

The recent outlook of leishmaniasis as a global public health concern coupled with the reportage of resistance and lack of efficacy of most antileishmanial drugs calls for a concerted effort to find new leads. The study combined In silico and in vitro approaches to identify novel potential synthetic small-molecule inhibitors targeting the Leishmania donovani sterol methyltransferase (LdSMT). The LdSMT enzyme in the ergosterol biosynthetic pathway is required for the parasite’s membrane fluidity, distribution of membrane proteins, and control of the cell cycle. The lack of LdSMT homologue in the human host and its conserved nature among all Leishmania parasites makes it a viable target for future antileishmanial drugs. Initially, six known inhibitors of LdSMT with IC50 < 10 μM were used to generate a pharmacophore model with a score of 0.9144 using LigandScout. The validated model was used to screen a synthetic library of 95,630 compounds obtained from InterBioScreen limited. Twenty compounds with pharmacophore fit scores above 50 were docked against the modelled three-dimensional structure of LdSMT using AutoDock Vina. Consequently, nine compounds with binding energies ranging from −7.5 to −8.7 kcal/mol were identified as potential hit molecules. Three compounds comprising STOCK6S-06707, STOCK6S-84928, and STOCK6S-65920 with respective binding energies of −8.7, −8.2, and −8.0 kcal/mol, lower than 22,26-azasterol (−7.6 kcal/mol), a known LdSMT inhibitor, were selected as plausible lead molecules. Molecular dynamics simulation studies and molecular mechanics Poisson–Boltzmann surface area calculations showed that the residues Asp25 and Trp208 were critical for ligand binding. The compounds were also predicted to have antileishmanial activity with reasonable pharmacological and toxicity profiles. When the antileishmanial activity of the three hits was evaluated in vitro against the promastigotes of L. donovani, mean half-maximal inhibitory concentrations (IC50) of 21.9 ± 1.5 μM (STOCK6S-06707), 23.5 ± 1.1 μM (STOCK6S-84928), and 118.3 ± 5.8 μM (STOCK6S-65920) were obtained. Furthermore, STOCK6S-84928 and STOCK6S-65920 inhibited the growth of Trypanosoma brucei, with IC50 of 14.3 ± 2.0 μM and 18.1 ± 1.4 μM, respectively. The identified compounds could be optimised to develop potent antileishmanial therapeutic agents.

The effect of Ebola virus disease (EVD) is fatal and devastating, necessitating several efforts to identify potent biotherapeutic molecules. This review seeks to provide perspectives on complementing existing work on Ebola virus (EBOV) by discussing the role of machine learning (ML) techniques in the prediction of small molecule inhibitors of EBOV. Different ML algorithms have been used to predict anti-EBOV compounds, including Bayesian, support vector machine, and random forest algorithms, which present strong models with credible outcomes. The use of deep learning models for predicting anti-EBOV molecules is underutilized; therefore, we discuss how such models could be leveraged to develop fast, efficient, robust, and novel algorithms to aid in the discovery of anti-EBOV drugs. We further discuss the deep neural network as a plausible ML algorithm for predicting anti-EBOV compounds. We also summarize the plethora of data sources necessary for ML predictions in the form of systematic and comprehensive high-dimensional data. With ongoing efforts to eradicate EVD, the application of artificial intelligence-based ML to EBOV drug discovery research can promote data-driven decision making and may help to reduce the high attrition rates of compounds in the drug development pipeline.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a serious global challenge requiring urgent and permanent therapeutic solutions. These solutions can only be engineered if the patterns and rate of mutations of the virus can be elucidated. Predicting mutations and the structure of proteins based on these mutations have become necessary for early drug and vaccine design purposes in anticipation of future viral mutations. The amino acid composition (AAC) of proteomes and individual viral proteins provide avenues for exploitation since AACs have been previously used to predict structure, shape and evolutionary rates. Herein, the frequency of amino acid residues found in 1637 complete proteomes belonging to 11 SARS-CoV-2 variants/lineages were analyzed. Leucine is the most abundant amino acid residue in the SARS-CoV-2 with an average AAC of 9.658% while tryptophan had the least abundance of 1.11%. The AAC and ranking of lysine and glycine varied in the proteome. For some variants, glycine had higher frequency and AAC than lysine and vice versa in other variants. Tryptophan was also observed to be the most intolerant to mutation in the various proteomes for the variants used. A correlogram revealed a very strong correlation of 0.999992 between B.1.525 (Eta) and B.1.526 (Iota) variants. Furthermore, isoleucine and threonine were observed to have a very strong negative correlation of −0.912, while cysteine and isoleucine had a very strong positive correlation of 0.835 at p < 0.001. Shapiro-Wilk normality test revealed that AAC values for all the amino acid residues except methionine showed no evidence of non-normality at p < 0.05. Thus, AACs of SARS-CoV-2 variants can be predicted using probability and z-scores. AACs may be beneficial in classifying viral strains, predicting viral disease types, members of protein families, protein interactions and for diagnostic purposes. They may also be used as a feature along with other crucial factors in machine-learning based algorithms to predict viral mutations. These mutation-predicting algorithms may help in developing effective therapeutics and vaccines for SARS-CoV-2.

Broadly neutralizing antibodies (bNAbs) are potent in neutralizing a wide range of HIV strains. VRC01 is a CD4-binding-site (CD4-bs) class of bNAbs that binds to the conserved CD4-binding region of HIV-1 envelope (env) protein. Natural products that mimic VRC01 bNAbs by interacting with the conserved CD4-binding regions may serve as a new generation of HIV-1 entry inhibitors by being broadly reactive and potently neutralizing. This study aimed to identify compounds that mimic VRC01 by interacting with the CD4-bs of HIV-1 gp120 and thereby inhibiting viral entry into target cells. Libraries of purchasable natural products were virtually screened against clade A/E recombinant 93TH057 (PDB: 3NGB) and clade B (PDB ID: 3J70) HIV-1 env protein. Protein–ligand interaction profiling from molecular docking and dynamics simulations showed that the compounds had intermolecular hydrogen and hydrophobic interactions with conserved amino acid residues on the CD4-binding site of recombinant clade A/E and clade B HIV-1 gp120. Four potential lead compounds, NP-005114, NP-008297, NP-007422, and NP-007382, were used for cell-based antiviral infectivity inhibition assay using clade B (HXB2) env pseudotype virus (PV). The four compounds inhibited the entry of HIV HXB2 pseudotype viruses into target cells at 50% inhibitory concentrations (IC50) of 15.2 µM (9.7 µg/mL), 10.1 µM (7.5 µg/mL), 16.2 µM (12.7 µg/mL), and 21.6 µM (12.9 µg/mL), respectively. The interaction of these compounds with critical residues of the CD4-binding site of more than one clade of HIV gp120 and inhibition of HIV-1 entry into the target cell demonstrate the possibility of a new class of HIV entry inhibitors.

Leishmania is a parasitic protozoon responsible for the neglected tropical disease Leishmaniasis. Approximately, 350 million people are susceptible and close to 70,000 death cases globally are reported annually. The lack of effective leishmanicides, the emergence of drug resistance and toxicity concerns necessitate the pursuit for effective antileishmanial drugs. Natural compounds serve as reservoirs for discovering new drugs due to their chemical diversity. Hardwickiic acid (HA) isolated from the stembark of Croton sylvaticus was evaluated for its leishmanicidal potential against Leishmania donovani and L. major promastigotes. The susceptibility of the promastigotes to HA was determined using the 3-[4,5-dimethylthiazol-2-yl]-2,5-diphenyltetrazolium bromide/phenazine methosulfate colorimetric assay with Amphotericin B serving as positive control. HA showed a significant antileishmanial activity on L. donovani promastigotes with an IC50 value of 31.57± 0.06 µM with respect to the control drug, amphotericin B with IC50 of 3.35 ± 0.14 µM). The cytotoxic activity was observed to be CC50 = 247.83 ± 6.32 µM against 29.99 ± 2.82 µM for curcumin, the control, resulting in a selectivity index of SI = 7.85. Molecular modeling, docking and dynamics simulations of selected drug targets corroborated the observed antileishmanial activity of HA. Novel insights into the mechanisms of binding were obtained for trypanothione reductase (TR), pteridine reductase 1 (PTR1), and glutamate cysteine ligase (GCL). The binding affinity of HA to the drug targets LmGCL, LmPTR1, LdTR, LmTR, LdGCL, and LdPTR1 were obtained as -8.0, -7.8, -7.6, -7.5, -7.4 and -7.1 kcal/mol, respectively. The role of Lys16, Ser111, and Arg17 as critical residues required for binding to LdPTR1 was reinforced. HA was predicted as a Caspase-3 stimulant and Caspase-8 stimulant, implying a possible role in apoptosis, which was shown experimentally that HA induced parasite death by loss of membrane integrity. HA was also predicted as antileishmanial molecule corroborating the experimental activity. Therefore, HA is a promising antileishmanial molecule worthy of further development as a biotherapeutic agent.

Organometallic compounds are molecules that contain at least one metal-carbon bond. Due to resistance of the Plasmodium parasite to traditional organic antimalarials, the use of organometallic compounds has become widely adopted in antimalarial drug discovery. Ferroquine, which was developed due to the emergence of chloroquine resistance, is currently the most advanced organometallic antimalarial drug and has paved the way for the development of new organometallic antimalarials. In this review, a general overview of organometallic antimalarial compounds and their antimalarial activity in comparison to purely organic antimalarials are presented. Furthermore, recent developments in the field are discussed, and future applications of this emerging class of therapeutics in antimalarial drug discovery are suggested.

Onchocerciasis is the leading cause of blindness and severe skin lesions which remain a major public health problem, especially in tropical areas. The widespread use of antibiotics and the long duration required for effective treatment continues to add to the increasing global menace of multi-resistant pathogens. Onchocerca volvulus harbors the endosymbiont bacteria Wolbachia, essential for the normal development of embryos, larvae and long-term survival of the adult worm, O. volvulus. We report here results of using structure-based drug design (SBDD) approach aimed at identifying potential novel Wolbachia inhibitors from natural products against the Wolbachia surface protein (WSP). The protein sequence of the WSP with UniProtKB identifier Q0RAI4 was used to model the three-dimensional (3D) structure via homology modelling techniques using three different structure-building algorithms implemented in Modeller, I-TASSER and Robetta. Out of the 15 generated models of WSP, one was selected as the most reasonable quality model which had 82, 15.5, 1.9 and 0.5% of the amino acid residues in the most favored regions, additionally allowed regions, generously allowed regions and disallowed regions, respectively, based on the Ramachandran plot. High throughput virtual screening was performed via Autodock Vina with a library comprising 42,883 natural products from African and Chinese databases, including 23 identified anti-Onchocerca inhibitors. The top six compounds comprising ZINC000095913861, ZINC000095486235, ZINC000035941652, NANPDB4566, acetylaleuritolic acid and rhemannic acid had binding energies of −12.7, −11.1, −11.0, −11, −10.3 and −9.5 kcal/mol, respectively. Molecular dynamics simulations including molecular mechanics Poisson-Boltzmann (MMPBSA) calculations reinforced the stability of the ligand-WSP complexes and plausible binding mechanisms. The residues Arg45, Tyr135, Tyr148 and Phe195 were predicted as potential novel critical residues required for ligand binding in pocket 1. Acetylaleuritolic acid and rhemannic acid (lantedene A) have previously been shown to possess anti-onchocercal activity. This warrants the need to evaluate the anti-WSP activity of the identified molecules. The study suggests the exploitation of compounds which target both pockets 1 and 2, by investigating their potential for effective depletion of Wolbachia. These compounds were predicted to possess reasonably good pharmacological profiles with insignificant toxicity and as drug-like. The compounds were computed to possess biological activity including antibacterial, antiparasitic, anthelmintic and anti-rickettsials. The six natural products are potential novel antiwolbachial agents with insignificant toxicities which can be explored further as filaricides for onchocerciasis.

Ebola virus (EBOV) is one of the most lethal pathogens that can infect humans. The Ebola viral protein VP35 (EBOV VP35) inhibits host IFN-α/β production by interfering with host immune responses to viral invasion and is thus considered as a plausible drug target. The aim of this study was to identify potential novel lead compounds against EBOV VP35 using computational techniques in drug discovery. The 3D structure of the EBOV VP35 with PDB ID: 3FKE was used for molecular docking studies. An integrated library of 7675 African natural product was pre-filtered using ADMET risk, with a threshold of 7 and, as a result, 1470 ligands were obtained for the downstream molecular docking using AutoDock Vina, after an energy minimization of the protein via GROMACS. Five known inhibitors, namely, amodiaquine, chloroquine, gossypetin, taxifolin and EGCG were used as standard control compounds for this study. The area under the curve (AUC) value, evaluating the docking protocol obtained from the receiver operating characteristic (ROC) curve, generated was 0.72, which was considered to be acceptable. The four identified potential lead compounds of NANPDB4048, NANPDB2412, ZINC000095486250 and NANPDB2476 had binding affinities of −8.2, −8.2, −8.1 and −8.0 kcal/mol, respectively, and were predicted to possess desirable antiviral activity including the inhibition of RNA synthesis and membrane permeability, with the probable activity (Pa) being greater than the probable inactivity (Pi) values. The predicted anti-EBOV inhibition efficiency values (IC50), found using a random forest classifier, ranged from 3.35 to 11.99 μM, while the Ki values ranged from 0.97 to 1.37 μM. The compounds NANPDB4048 and NANPDB2412 had the lowest binding energy of −8.2 kcal/mol, implying a higher binding affinity to EBOV VP35 which was greater than those of the known inhibitors. The compounds were predicted to possess a low toxicity risk and to possess reasonably good pharmacological profiles. Molecular dynamics (MD) simulations of the protein–ligand complexes, lasting 50 ns, and molecular mechanisms Poisson-Boltzmann surface area (MM-PBSA) calculations corroborated the binding affinities of the identified compounds and identified novel critical interacting residues. The antiviral potential of the molecules could be confirmed experimentally, while the scaffolds could be optimized for the design of future novel anti-EBOV chemotherapeutics.

Severely ill coronavirus disease 2019 (COVID-19) patients show elevated concentrations of pro-inflammatory cytokines, a situation commonly known as a cytokine storm. The p38 MAPK receptor is considered a plausible therapeutic target because of its involvement in the platelet activation processes leading to inflammation. This study aimed to identify potential natural product-derived inhibitory molecules against the p38α MAPK receptor to mitigate the eliciting of pro-inflammatory cytokines using computational techniques. The 3D X-ray structure of the receptor with PDB ID 3ZS5 was energy minimized using GROMACS and used for molecular docking via AutoDock Vina. The molecular docking was validated with an acceptable area under the curve (AUC) of 0.704, which was computed from the receiver operating characteristic (ROC) curve. A compendium of 38,271 natural products originating from Africa and China together with eleven known p38 MAPK inhibitors were screened against the receptor. Four potential lead compounds ZINC1691180, ZINC5519433, ZINC4520996 and ZINC5733756 were identified. The compounds formed strong intermolecular bonds with critical residues Val38, Ala51, Lys53, Thr106, Leu108, Met109 and Phe169. Additionally, they exhibited appreciably low binding energies which were corroborated via molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) calculations. The compounds were also predicted to have plausible pharmacological profiles with insignificant toxicity. The molecules were also predicted to be anti-inflammatory, kinase inhibitors, antiviral, platelet aggregation inhibitors, and immunosuppressive, with probable activity (Pa) greater than probable inactivity (Pi). ZINC5733756 is structurally similar to estradiol with a Tanimoto coefficient value of 0.73, which exhibits anti-inflammatory activity by targeting the activation of Nrf2. Similarly, ZINC1691180 has been reported to elicit anti-inflammatory activity in vitro. The compounds may serve as scaffolds for the design of potential biotherapeutic molecules against the cytokine storm associated with COVID-19.

The huge burden of leishmaniasis caused by the trypanosomatid protozoan parasite Leishmania is well known. This illness was included in the list of neglected tropical diseases targeted for elimination by the World Health Organization. However, the increasing evidence of resistance to existing antimonial drugs has made the eradication of the disease difficult to achieve, thus warranting the search for new drug targets. We report here studies that used computational methods to identify inhibitors of receptors from natural products. The cell division cycle-2-related kinase 12 (CRK12) receptor is a plausible drug target against Leishmania donovani. This study modelled the 3D molecular structure of the L. donovani CRK12 (LdCRK12) and screened for small molecules with potential inhibitory activity from African flora. An integrated library of 7722 African natural product-derived compounds and known inhibitors were screened against the LdCRK12 using AutoDock Vina after performing energy minimization with GROMACS 2018. Four natural products, namely sesamin (NANPDB1649), methyl ellagic acid (NANPDB1406), stylopine (NANPDB2581), and sennecicannabine (NANPDB6446) were found to be potential LdCRK12 inhibitory molecules. The molecular docking studies revealed two compounds NANPDB1406 and NANPDB2581 with binding affinities of −9.5 and −9.2 kcal/mol, respectively, against LdCRK12 which were higher than those of the known inhibitors and drugs, including GSK3186899, amphotericin B, miltefosine, and paromomycin. All the four compounds were predicted to have inhibitory constant (Ki) values ranging from 0.108 to 0.587 μM. NANPDB2581, NANPDB1649 and NANPDB1406 were also predicted as antileishmanial with Pa and Pi values of 0.415 and 0.043, 0.391 and 0.052, and 0.351 and 0.071, respectively. Molecular dynamics simulations coupled with molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) computations reinforced their good binding mechanisms. Most compounds were observed to bind in the ATP binding pocket of the kinase domain. Lys488 was predicted as a key residue critical for ligand binding in the ATP binding pocket of the LdCRK12. The molecules were pharmacologically profiled as druglike with inconsequential toxicity. The identified molecules have scaffolds that could form the backbone for fragment-based drug design of novel leishmanicides but warrant further studies to evaluate their therapeutic potential.

The coronavirus disease 2019 (COVID-19) pandemic caused by the severe acute respiratory syndrome virus 2 (SARS-CoV-2) has impacted negatively on public health and socioeconomic status, globally. Although, there are currently no specific drugs approved, several existing drugs are being repurposed, but their successful outcomes are not guaranteed. Therefore, the search for novel therapeutics remains a priority. We screened for inhibitors of the SARS-CoV-2 main protease and the receptor-binding domain of the spike protein from an integrated library of African natural products, compounds generated from machine learning studies and antiviral drugs using AutoDock Vina. The binding mechanisms between the compounds and the proteins were characterized using LigPlot+ and molecular dynamics simulations techniques. The biological activities of the hit compounds were also predicted using a Bayesian-based approach. Six potential bioactive molecules NANPDB2245, NANPDB2403, fusidic acid, ZINC000095486008, ZINC0000556656943 and ZINC001645993538 were identified, all of which had plausible binding mechanisms with both viral receptors. Molecular dynamics simulations, including molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) computations revealed stable protein-ligand complexes with all the compounds having acceptable free binding energies de novo inhibitors. The compounds are proposed as worthy of further in vitro assaying and as scaffolds for the development of novel SARS-CoV-2 therapeutic molecules.

Ulcers due to infections with Mycobacterium ulcerans are characterized by complete lack of wound healing processes, painless, an underlying bed of host dead cells and undermined edges due to necrosis. Mycolactone, a macrolide produced by the mycobacterium, is believed to be the toxin responsible. Of interest and relevance is the knowledge that Buruli ulcer (BU) patients remember experiencing trauma previously at the site of the ulcers, suggesting an impairment of wound healing processes, the plausible effect due to the toxin. Wound healing processes involve activation of the blood platelets to release the contents of the dense granules mainly serotonin, calcium ions, and ADP/ATP by exocytosis into the bloodstream. The serotonin release results in attracting more platelets and mast cells to the wound site, with the mast cells also undergoing degranulation, releasing compounds into the bloodstream by exocytosis. Recent work has identified interference in the co-translational translocation of many secreted proteins via the endoplasmic reticulum and cell death involving Wiskott-Aldrich syndrome protein (WASP), Sec61, and angiotensin II receptors (AT2R). We hypothesized that mycolactone by being lipophilic, passively crosses cell membranes and binds to key proteins that are involved in exocytosis by platelets and mast cells, thus inhibiting the initiation of wound healing processes. Based on this, molecular docking studies were performed with mycolactone against key soluble n-ethylmaleimide-sensitive factor attachment protein receptor (SNARE) proteins and regulators, namely Vesicle-associated membrane protein (VAMP8), Synaptosomal-associated protein (SNAP23, syntaxin 11, Munc13-4 (its isoform Munc13-1 was used), and Munc18b; and also against known mycolactone targets (Sec61, AT2R, and WASP). Munc18b was shown to be a plausible mycolactone target after the molecular docking studies with binding affinity of −8.5 kcal/mol. Structural studies and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) binding energy calculations of the mycolactone and Munc18b complex was done with 100 ns molecular dynamics simulations using GROMACS. Mycolactone binds strongly to Munc18b with an average binding energy of −247.571 ± 37.471 kJ/mol, and its presence elicits changes in the structural conformation of the protein. Analysis of the binding interactions also shows that mycolactone interacts with Arg405, which is an important residue of Munc18b, whose mutation could result in impaired granule exocytosis. These findings consolidate the possibility that Munc18b could be a target of mycolactone. The implication of the interaction can be experimentally evaluated to further understand its role in granule exocytosis impairment in Buruli ulcer.