Buruli ulcer (BU) is caused by Mycobacterium ulcerans and is predominant in both tropical and subtropical regions. The neglected debilitating disease is characterized by chronic necrotizing skin lesions attributed to a mycolactone, which is a macrolide toxin secreted by M. ulcerans. The preferred treatment is surgical excision of the lesions followed by a prolonged combination antibiotic therapy using existing drugs such as rifampicin and streptomycin or clarithromycin. These antibiotics appear not to be adequately potent and efficacious against persistent and late stage ulcers. In addition, emerging drug resistance to treatment poses great challenges. There is a need to identify novel natural product-derived lead compounds, which are potent and efficacious for the treatment of Buruli ulcer. Natural products present a rich diversity of chemical compounds with proven activity against various infectious diseases, and therefore, are considered in this study. This study sought to computationally predict natural product-derived lead compounds with the potential to be developed further into potent drugs with better therapeutic efficacy than the existing anti-buruli ulcer compounds. The three-dimensional (3D) structure of Isocitrate lyase (ICL) of Mycobacterium ulcerans was generated using homology modeling and was further scrutinized with molecular dynamics simulations. A library consisting of 885 compounds retrieved from the AfroDb database was virtually screened against the validated ICL model using AutoDock Vina. AfroDb is a compendium of “drug-like” and structurally diverse 3D structures of natural products originating from different geographical regions in Africa. The molecular docking with the ICL model was validated by computing a Receiver Operating Characteristic (ROC) curve with a reasonably good Area Under the Curve (AUC) value of 0.89375. Twenty hit compounds, which docked firmly within the active site pocket of the ICL receptor, were assessed via in silico bioactivity and pharmacological profiling. The three compounds, which emerged as potential novel leads, comprise ZINC38143792 (Euscaphic acid), ZINC95485880, and ZINC95486305 with reasonable binding energies (high affinity) of −8.6, −8.6, and −8.8 kcal/mol, respectively. Euscaphic acid has been reported to show minimal inhibition against a drug-sensitive strain of M. tuberculosis. The other two leads were both predicted to possess dermatological activity while one was antibacterial. The leads have shown promising results pertaining to efficacy, toxicity, pharmacokinetic, and safety. These leads can be experimentally characterized to assess their anti-mycobacterial activity and their scaffolds may serve as rich skeletons for developing anti-buruli ulcer drugs.
15 September 2021 in Informatics in Medicine Unlocked
The emergence of drug resistance against the known hookworm drugs namely albendazole and mebendazole and their reduced efficacies necessitate the need for new drugs. Chemically diverse natural products present plausible templates to augment hookworm drug discovery. The present work utilized pharmacoinformatics techniques to predict African natural compounds ZINC95486082, ZINC95486052 and euphohelionon as potential inhibitory molecules of the hookworm Necator americanus β tubulin gene. A library of 3390 compounds was screened against a homology-modelled structure of β tubulin. The docking results obtained from AutoDock Vina was validated with an acceptable area under the curve (AUC) of 0.714 computed from the receiver operating characteristic (ROC) curve. The three selected compounds had favourable binding affinities and were predicted to form no interactions with the resistance-associated mutations Phe167, Glu198 and Phe200. The compounds were predicted as anthelmintics using a Bayesian-based technique and were pharmacologically profiled to be druglike. Further molecular dynamics simulations and MM-PBSA calculations showed the compounds as promising anthelmintic drug leads. Novel critical residues comprising Leu246, Asn247 and Asn256 were also predicted for binding. Euphohelionon was selected as a template for the de novo fragment-based design of five compounds labelled A1, A2, A3, A4 and A5; with four of them having SAscore values below 6, denoting easy synthesis. All the five de novo molecules docked firmly in the binding pocket of the β tubulin with no binding interactions with the three known resistance mutation residues. Binding energies of −8.2, −7.6, −7.3, −7.2 and −6.8 kcal/mol were obtained for A1, A2, A3, A4 and A5, respectively. The identified compounds can serve as treasure troves from which future potent anthelmintics can be designed. The current study strives to assuage the hookworm disease burden, especially making available molecules with the potential to circumvent the chemoresistance.
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.
Objective: Falciparum malaria is a major global health concern, affecting more than half of the world’s population and causing over half a million deaths annually. Red cell invasion is a crucial step in the parasite’s life cycle, where the parasite invade human erythrocytes to sustain infection and ensure survival. Two parasite proteins, Apical Membrane Antigen 1 (AMA-1) and Rhoptry Neck Protein 2 (RON2), are involved in tight junction formation, which is an essential step in parasite invasion of the red blood cell. Targeting the AMA-1 and RON2 interaction with inhibitors halts the formation of the tight junction, thereby preventing parasite invasion, which is detrimental to parasite survival. This study leverages machine learning (ML) to predict potential small molecule inhibitors of the AMA-1–RON2 interaction, providing putative antimalaria compounds for further chemotherapeutic exploration. Method: Data was retrieved from the PubChem database (AID 720542), comprising 364,447 inhibitors and non-inhibitors of the AMA-1–RON2 interaction. The data was processed by computing Morgan fingerprints and divided into training and testing with an 80:20 ratio, and the classes in the training data were balanced using the Synthetic Minority Oversampling Technique. Five ML models developed comprised Random Forest (RF), Gradient Boost Machines (GBMs), CatBoost (CB), AdaBoost (AB) and Support Vector Machine (SVM). The performances of the models were evaluated using accuracy, F1 score, and receiver operating characteristic—area under the curve (ROC-AUC) and validated using held-out data and a y-randomization test. An applicability domain analysis was carried out using the Tanimoto distance with a threshold set at 0.04 to ascertain the sample space where the models predict with confidence. Results: The GBMs model emerged as the best, achieving 89% accuracy and a ROC-AUC of 92%. CB and RF had accuracies of 88% and 87%, and ROC-AUC scores of 93% and 91%, respectively. Conclusions: Experimentally validated inhibitors of the AMA-1–RON2 interaction could serve as starting blocks for the next-generation antimalarial drugs. The models were deployed as a web-based application, known as PLASMOpred.