Recent advances in quantitative structure-activity relationship models of antimalarial drugs

In silico modeling revealed phytomolecules derived from Cymbopogon citratus (DC.) leaf extract as promising candidates for malaria therapy

The emergence of varying levels of resistance to currently available antimalarial drugs significantly threatens global health. This factor heightens the urgency to explore bioactive compounds from natural products with a view to discovering and developing newer antimalarial drugs with novel mode of actions. Therefore, we evaluated the inhibitory effects of sixteen phytocompounds from Cymbopogon citratus leaf extract against Plasmodium falciparum drug targets such as P. falciparum circumsporozoite protein (PfCSP), P. falciparum merozoite surface protein 1 (PfMSP1) and P. falciparum erythrocyte membrane protein 1 (PfEMP1). In silico approaches including molecular docking, pharmacophore modeling and 3D-QSAR were adopted to analyze the inhibitory activity of the compounds under consideration. The molecular docking results indicated that a compound swertiajaponin from C. citratus exhibited a higher binding affinity (-7.8 kcal/mol) to PfMSP1 as against the standard artesunate-amodiaquine (-6.6 kcal/mol). Swertiajaponin also formed strong hydrogen bond interactions with LYS29, CYS30, TYR34, ASN52, GLY55 and CYS28 amino acid residues. In addition, quercetin another compound from C. citratus exhibited significant binding energies -6.8 and -8.3 kcal/mol with PfCSP and PfEMP1, respectively but slightly lower than the standard artemether-lumefantrine with binding energies of -7.4 kcal/mol against PfCSP and -8.7 kcal/mol against PfEMP1. Overall, the present study provides evidence that swertiajaponin and other phytomolecules from C. citratus have modulatory properties toward P. falciparum drug targets and thus may warrant further exploration in early drug discovery efforts against malaria. Furthermore, these findings lend credence to the folkloric use of C. citratus for malaria treatment.Communicated by Ramaswamy H. Sarma.

Structure-activity relationships andz interindividual variability of drug responses: pharmacogenomics with antimicrobial drugs as a paradigm

Adverse drug reactions represent a major health burden because they cause notable patient morbidity and mortality. From this viewpoint, several strategies have been developed to prevent or reduce adverse drug reactions. One such strategy is the use of pharmacogenomics. Interindividual variability in drug response and adverse effects is mainly attributable to genetic variation in enzymes such as sulfotransferases and cytochrome P450s. The current narrative review discusses the relationship between the structure and activity of drugs. Specifically, the activity of drugs can be increased and/or their adverse effects can be reduced by altering specific positions in their structures.

A comprehension on structure guided alignment dependent 3D-QSAR modelling, and molecular dynamics simulation on 2,4-thiazolidinediones as aldose reductase inhibitors for the management of diabetic complications

Aldose reductase is an oxo-reductase enzyme belonging to the aldo-keto reductase class. Compounds having thiazolidine-2,4-dione scaffold are reported as potential aldose reductase inhibitors for diabetic complications. The present work uses structure-guided alignment-dependent Gaussian field- and atom-based 3D-QSAR on a dataset of 84 molecules. 3D-QSAR studies on two sets of dataset alignment have been carried out to understand the favourable and unfavourable structural features influencing the affinity of these inhibitors towards the enzyme. Using common pharmacophore hypotheses, the five-point pharmacophores for aldose reductase favourable features were generated. The molecular dynamics simulations (up to 100 ns) were performed for the potent molecule from each alignment set (compounds 24 and 65) compared to reference standard tolrestat and epalrestat to study target-ligand complexes' binding energy and stability. Compound 65 was most stable with better interactions in the aldose reductase binding pocket than tolrestat. The MM-PBSA study suggests compound 65 possessed better binding energy than reference standard tolrestat, i.e. -87.437 ± 19.728 and -73.424 ± 12.502 kJ/mol, respectively. The generated 3D-QSAR models provide information about structure-activity relationships and ligand-target binding energy. Target-specific stability data from MD simulation would be helpful for rational compound design with better aldose reductase activity.Communicated by Ramaswamy H. Sarma.

Unsupervised machine learning, QSAR modelling and web tool development for streamlining the lead identification process of antimalarial flavonoids

Identification of lead compounds with the traditional laboratory approach is expensive and time-consuming. Nowadays, in silico techniques have emerged as a promising approach for lead identification. In this study, we aim to develop robust and predictive 2D-QSAR models to identify lead flavonoids by predicting the IC₅₀ against Plasmodium falciparum. We applied machine learning algorithms (Principal component analysis followed by K-means clustering) and Pearson correlation analysis to select 9 molecular descriptors (MDs) for model building. We selected and validated the three best QSAR models after execution of multiple linear regression (MLR) 100 times with different combinations of MDs. The developed models have fulfilled the five principles for QSAR models as specified by the Organization for Economic Co-operation and Development. The outcome of the study is a reliable and sustainable in silico method of IC₅₀ (Mean ± SD) prediction that will positively impact the antimalarial drug development process by reducing the money and time required to identify potential antimalarial lead compounds from the class of flavonoids. We also developed a web tool (JazQSAR, https://etflin.com/news/4) to offer an easily accessible platform for the developed QSAR models.

The importance of selecting crystal form for triazole fungicide tebuconazole to enhance its botryticidal activity

The growing evidences of resistant fungi stimulate fully understanding tebuconazole regarding its crystal structure on fungicidal activity. In this study, the crystal structures of six technical tebuconazoles (BX, HH, JP, QZ, SJ, and YT) were characterized by using high-resolution X-ray powder diffraction and three-dimensional crystal structure modeling. A structure-activity relationship of the tebuconazoles on the susceptible (HLS and YJS) or resistant (XHR) Botrytis cinerea isolates was analyzed, the differential tricarboxylic acid (TCA) cycle metabolism was determined, and molecular docking with sterol 14α-demethylase (CYP51) was performed. The results showed that tebuconazole existed in three types of crystal forms: an overlapping-pair conformation, a side-by-side-pair conformation, and a parallel-pair conformation. QZ with the parallel-pair conformation and the minimum crystal cell volume exhibited a higher activity and a lower resistant level. XHR possessed a higher content of TCA cycle metabolites and phosphate than YJS, but the exposure to QZ significantly reduced the contents of citrate, isocitrate, α-ketoglutarate and oxaloacetate in XHR, as did the exposure to other technical tebuconazoles. Moreover, the point mutations F487L, G464S, and G443S altered the binding properties of chiral stereoscopic R-QZ with CYP51 protein. Especially the G443S mutation promoted a weak linking of R-QZ with LEU380 and TYR126, and greatly slashed the binding action at lower docking score. In conclusion, our results evidenced an efficient crystal conformation of tebuconazole to improve botryticidal activity and a potential adaptability of B. cinerea to tebuconazole inhibition in TCA cycle metabolism and CYP51 protein mutation.

Systematic review on the application of machine learning to quantitative structure-activity relationship modeling against Plasmodium falciparum

Malaria accounts for over two million deaths globally. To flatten this curve, there is a need to develop new and high potent drugs against Plasmodium falciparum. Some major challenges include the dearth of suitable animal models for anti-P. falciparum assays, resistance to first-line drugs, lack of vaccines and the complex life cycle of Plasmodium. Gladly, newer approaches to antimalarial drug discovery have emerged due to the release of large datasets by pharmaceutical companies. This review provides insights into these new approaches to drug discovery covering different machine learning tools, which enhance the development of new compounds. It provides a systematic review on the use and prospects of machine learning in predicting, classifying and clustering IC₅₀ values of bioactive compounds against P. falciparum. The authors identified many machine learning tools yet to be applied for this purpose. However, Random Forest and Support Vector Machines have been extensively applied though on a limited dataset of compounds.

First total synthesis of hoshinoamide A

Hoshinoamides A, B and C, linear lipopeptides, were isolated from the marine cyanobacterium Caldora penicillata, with potent antiplasmodial activity against chloroquine-sensitive Plasmodium falciparum. Herein, we describe the first total synthesis of hoshinoamide A by the combination of liquid-phase and solid-phase peptide synthesis. Liquid-phase synthesis is to improve the coupling yield of ʟ-Val³ and N-Me-ᴅ-Phe². Connecting other amino acids efficiency and convergence is achieved by solid-state synthesis. Our synthetic strategy could synthesize the target peptide in high yield with good purity.

2D-quantitative structure-activity relationships model using PLS method for anti-malarial activities of anti-haemozoin compounds

Background

Emergence of cross-resistance to current anti-malarial drugs has led to an urgent need for identification of potential compounds with novel modes of action and anti-malarial activity against the resistant strains. One of the most promising therapeutic targets of anti-malarial agents related to food vacuole of malaria parasite is haemozoin, a product formed by the parasite through haemoglobin degradation.

Methods

With this in mind, this study developed two-dimensional-quantitative structure–activity relationships (QSAR) models of a series of 21 haemozoin inhibitors to explore the useful physicochemical parameters of the active compounds for estimation of anti-malarial activities. The 2D-QSAR model with good statistical quality using partial least square method was generated after removing the outliers.

Results

Five two-dimensional descriptors of the training set were selected: atom count (a_ICM); adjacency and distance matrix descriptor (GCUT_SLOGP_2: the third GCUT descriptor using atomic contribution to logP); average total charge sum (h_pavgQ) in pKa prediction (pH = 7); a very low negative partial charge, including aromatic carbons which have a heteroatom-substitution in “ortho” position (PEOE_VSA-0) and molecular descriptor (rsynth: estimating the synthesizability of molecules as the fraction of heavy atoms that can be traced back to starting material fragments resulting from retrosynthetic rules), respectively. The model suggests that the anti-malarial activity of haemozoin inhibitors increases with molecules that have higher average total charge sum in pKa prediction (pH = 7). QSAR model also highlights that the descriptor using atomic contribution to logP or the distance matrix descriptor (GCUT_SLOGP_2), and structural component of the molecules, including topological descriptors does make for better anti-malarial activity.

Conclusions

The model is capable of predicting the anti-malarial activities of anti-haemozoin compounds. In addition, the selected molecular descriptors in this QSAR model are helpful in designing more efficient compounds against the P. falciparum 3D7A strain.