Exploring Structure-Activity Relationships with Three-Dimensional Matched Molecular Pairs-A Review

S. C, Gangadhar S, A. S, Murali M, Shivamallu C, Achar RR, Silina E, Stupin V, Manturova N, Shati AA, Alfaifi MY, Elbehairi SEI, Kollur SP, Amruthesh KN. Exploration of Type III effector Xanthomonas outer protein Q (XopQ) inhibitor from Picrasma quassioides as an antibacterial agent using chemoinformatics analysis

The present study was focused on exploring the efficient inhibitors of closed state (form) of type III effector Xanthomonas outer protein Q (XopQ) (PDB: 4P5F) from the 44 phytochemicals of Picrasma quassioides using cutting-edge computational analysis. Among them, Kumudine B showed excellent binding energy (-11.0 kcal/mol), followed by Picrasamide A, Quassidine I and Quassidine J with the targeted closed state of XopQ protein compared to the reference standard drug (Streptomycin). The molecular dynamics (MD) simulations performed at 300 ns validated the stability of top lead ligands (Kumudine B, Picrasamide A, and Quassidine I)-bound XopQ protein complex with slightly lower fluctuation than Streptomycin. The MM-PBSA calculation confirmed the strong interactions of top lead ligands (Kumudine B and QuassidineI) with XopQ protein, as they offered the least binding energy. The results of absorption, distribution, metabolism, excretion, and toxicity (ADMET) analysis confirmed that Quassidine I, Kumudine B and Picrasamide A were found to qualify most of the drug-likeness rules with excellent bioavailability scores compared to Streptomycin. Results of the computational studies suggested that Kumudine B, Picrasamide A, and Quassidine I could be considered potential compounds to design novel antibacterial drugs against X. oryzae infection. Further in vitro and in vivo antibacterial activities of Kumudine B, Picrasamide A, and Quassidine I are required to confirm their therapeutic potentiality in controlling the X. oryzae infection.

Camptothecin structure simplification elaborated new imidazo[2,1-b]quinazoline derivative as a human topoisomerase I inhibitor with efficacy against bone cancer cells and colon adenocarcinoma

Camptothecin is a pentacyclic natural alkaloid that inhibits the hTop1 enzyme involved in DNA transcription and cancer cell growth. Camptothecin structure pitfalls prompted us to design new congeners using a structure simplification strategy to reduce the ring extension number from pentacyclic to tetracyclic while maintaining potential stacking of the new compounds with the DNA base pairs at the Top1-mediated cleavage complex and aqueous solubility, as well as minimizing compound-liver toxicity. The principal axis of this study was the verification of hTop1 inhibiting activity as a possible mechanism of action and the elaboration of new simplified inhibitors with improved pharmacodynamic and pharmacokinetic profiling using three structure panels (A-C) of (isoquinolinoimidazoquinazoline), (imidazoquinazoline), and (imidazoisoquinoline), respectively. DNA relaxation assay identified five compounds as hTop1 inhibitors belonging to the imidazoisoquinolines 3a,b, the imidazoquinazolines 12, and the isoquinolinoimidazoquinazolines 7a,b. In an MTT cytotoxicity assay against different cancer cell lines, compound 12 was the most potent against HOS bone cancer cells (IC50 = 1.47 μM). At the same time, the other inhibitors had no detectable activity against any cancer cell type. Compound (12) demonstrated great penetrating power in the HOS cancer cells' 3D-multicellular tumor spheroid model. Bioinformatics research of the hTop1 gene revealed that the TP53 cell proliferative gene is in the network of hTop1. The finding is confirmed empirically using the gene expression assay that proved the increase in p53 expression. The impact of structure simplification on compound 12 profile, characterized by the absence of acute oral liver toxicity when compared to Doxorubicin as a standard inhibitor, the lethal dose measured on Swiss Albino female mice and reported at LD50 = 250 mg/kg, and therapeutic significance in reducing colon adenocarcinoma tumor volume by 75.36 % after five weeks of treatment with compound 12. The molecular docking solutions of the active CPT-based derivative 12 and the inactive congener 14 into the active site of hTop1 and the activity cliffing of such MMP directed us to recommend the addition of HBD and HBA variables to compound 12 imidazoquinazoline core scaffold to enhance the potency via hydrogen bond formation with the major groove amino acids (Asp533, Lys532) as well as maintaining the hydrogen bond with the minor groove amino acid Arg364.

Secondary metabolites of Trichoderma spp. as EGFR tyrosine kinase inhibitors: Evaluation of anticancer efficacy through computational approach

The present study explores the epidermal growth factor receptor (EGFR) tyrosine kinase inhibition efficacy of secondary metabolites in Trichoderma spp. through molecular docking, molecular dynamics (MD) simulation and MM-PBSA approach. The result of molecular docking confirmed that out of 200 metabolites screened, three metabolites such as Harzianelactone A, Pretrichodermamide G and Aspochalasin M, potentially bound with the active binding site of EGFR tyrosine kinase domain(PDB ID: 1M17) with a threshold docking score of ≤- 9.0 kcal/mol when compared with the standard EGFR inhibitor (Erlotinib). The MD simulation was run to investigate the potential for stable complex formation in EGFR tyrosine kinase domain-unbound/lead metabolite (Aspochalasin M)-bound/standard inhibitor (Erlotinib)-bound complex. The MD simulation analysis at 100 ns revealed that Aspochalasin M formed the stable complex with EGFR. Besides, the in silico predication of pharmacokinetic properties further confirmed that Aspochalasin M qualified the drug-likeness rules with no harmful side effects (viz., hERG toxicity, hepatotoxicity and skin sensitization), non-mutagenicity and favourable logBB value. Moreover, the BOILED-Egg model predicted that Aspochalasin M showed a higher gastrointestinal absorption with improved bioavailability when administered orally and removed from the central nervous system (CNS). The results of the computational studies concluded that Aspochalasin M possessed significant efficacy in binding EGFR's active sites compared to the known standard inhibitor (Erlotinib). Therefore, Aspochalasin M can be used as a possible anticancer drug candidate and further in vitro and in vivo experimental validation of Aspochalasin M of Trichoderma spp. are required to determine its anticancer potential.

Visualization of Topological Pharmacophore Space with Graph Edit Distance

A topological pharmacophore (TP) is a chemical graph-based pharmacophore representation, where nodes are pharmacophoric features (PF) and edges are topological distances between PFs. Previously proposed sparse pharmacophore graphs (SPhGs) for TPs were shown to be effective in identifying structurally different active compounds while maintaining the interpretability of the graphs. However, one limitation of using SPhGs as queries is that many structurally similar SPhGs can be identified from a set of active compounds, requiring the classification and visualization of SPhGs, followed by an understanding of the pharmacophore hypotheses. In this study, we propose a scheme for SPhG analysis based on dimensionality reduction techniques with the graph edit distance (GED) metric. This metric enables measuring similarities among SPhGs in a quantitative manner. The visualization of SPhGs, which themselves are the graphs shared by active compounds, can help us understand the pharmacophore hypotheses as well as the data set. As a proof-of-concept study, we generated two-dimensional SPhG-maps using three dimensionality reduction techniques for six biological targets. A comparison with other pharmacophore representations was also conducted. We demonstrated knowledge extraction (interpretation of the data set) from the generated maps. Our findings include a suitable mapping algorithm as well as a pharmacophore hypothesis analysis procedure using an SPhG-map.

Performance of Relative Binding Free Energy Calculations on an Automatically Generated Dataset of Halogen-Deshalogen Matched Molecular Pairs

In this study, we generated a matched molecular pair dataset of halogen/deshalogen compounds with reliable binding affinity data and structural binding mode information from public databases. The workflow includes automated system preparation and setup of free energy perturbation relative binding free energy calculations. We demonstrate the suitability of these datasets to investigate the performance of molecular mechanics force fields and molecular simulation algorithms for the purpose of in silico affinity predictions in lead optimization. Our datasets of a total of 115 matched molecular pairs show highly accurate binding free energy predictions with an average error of <1 kcal/mol despite the semi-automated calculation scheme. We quantify the accuracy of the optimized potential for liquid simulations (OPLS) force field to predict the effect of halogen addition to compounds, a commonly employed chemical modification in the design of drug-like molecules.

Structure-based molecular modeling in SAR analysis and lead optimization

In silico methods like molecular docking and pharmacophore modeling are established strategies in lead identification. Their successful application for finding new active molecules for a target is reported by a plethora of studies. However, once a potential lead is identified, lead optimization, with the focus on improving potency, selectivity, or pharmacokinetic parameters of a parent compound, is a much more complex task. Even though in silico molecular modeling methods could contribute a lot of time and cost-saving by rationally filtering synthetic optimization options, they are employed less widely in this stage of research. In this review, we highlight studies that have successfully used computer-aided SAR analysis in lead optimization and want to showcase sound methodology and easily accessible in silico tools for this purpose.

Discovery of potent and selective inhibitors of the Escherichia coli M1-aminopeptidase via multicomponent solid-phase synthesis of tetrazole-peptidomimetics

The Escherichia coli neutral M1-aminopeptidase (ePepN) is a novel target identified for the development of antimicrobials. Here we describe a solid-phase multicomponent approach which enabled the discovery of potent ePepN inhibitors. The on-resin protocol, developed in the frame of the Distributed Drug Discovery (D3) program, comprises the implementation of parallel Ugi-azide four-component reactions with resin-bound amino acids, thus leading to the rapid preparation of a focused library of tetrazole-peptidomimetics (TPMs) suitable for biological screening. By dose-response studies, three compounds were identified as potent and selective ePepN inhibitors, as little inhibitory effect was exhibited for the porcine ortholog aminopeptidase. The study allowed for the identification of the key structural features required for a high ePepN inhibitory activity. The most potent and selective inhibitor (TPM 11) showed a non-competitive inhibition profile of ePepN. We predicted that both diastereomers of compound TPM 11 bind to a site distinct from that occupied by the substrate. Theoretical models suggested that TPM 11 has an alternative inhibition mechanism that doesn't involve Zn coordination. On the other hand, the activity landscape analysis provided a rationale for our findings. Of note, compound TMP 2 showed in vitro antibacterial activity against Escherichia coli. Furthermore, none of the three identified inhibitors is a potent haemolytic agent, and only two compounds showed moderate cytotoxic activity toward the murine myeloma P3X63Ag cells. These results point to promising compounds for the future development of rationally designed TPMs as antibacterial agents.