Prospects of Indole derivatives as methyl transfer inhibitors: antimicrobial resistance managers
BMC Pharmacology and Toxicology volume 21, Article number: 33 (2020)
It is prudent that novel classes of antibiotics be urgently developed to manage the WHO prioritized multi-drug resistant (MDR) pathogens posing an unprecedented medical crisis. Simultaneously, multiple essential proteins have to be targeted to prevent easy resistance development.
An integration of structure-based virtual screening and ligand-based virtual screening was employed to explore the antimicrobial properties of indole derivatives from a compound database.
Whole-genome sequences of the target pathogens were aligned exploiting DNA alignment potential of MAUVE to identify putative common lead target proteins. S-adenosyl methionine (SAM) biosynthesizing MetK was taken as the lead target and various literature searches revealed that SAM is a critical metabolite. Furthermore, SAM utilizing CobA involved in the B12 biosynthesis pathway, Dam in the regulation of replication and protein expression, and TrmD in methylation of tRNA were also taken as drug targets. The ligand library of 715 indole derivatives chosen based on kinase inhibition potential of indoles was created from which 102 were pursued based on ADME/T scores. Among these, 5 potential inhibitors of MetK in N. gonorrhoeae were further expanded to molecular docking studies in MetK proteins of all nine pathogens among which 3 derivatives exhibited inhibition potential. These 3 upon docking in other SAM utilizing enzymes, CobA, Dam, and TrmD gave 2 potential compounds with multiple targets. Further, docking with human MetK homolog also showed probable inhibitory effects however SAM requirements can be replenished from external sources since SAM transporters are present in humans.
We believe these molecules 3-[(4-hydroxyphenyl)methyl]-6-(1H-indol-3-ylmethyl)piperazine-2,5-dione (ZINC04899565) and 1-[(3S)-3-[5-(1H-indol-3-ylmethyl)-1,3,4-oxadiazol-2-yl]pyrrolidin-1-yl]ethanone (ZINC49171024) could be a starting point to help develop broad-spectrum antibiotics against infections caused by N. gonorrhoeae, A. baumannii, C. coli, K. pneumoniae, E. faecium, H. pylori, P. aeruginosa, S. aureus and S. typhi.
The finding and development of antibiotics was a milestone in medical sciences that prevented fatality from simple infections. Unfortunately, the emergence of antibiotic-resistant strains among these pathogens appears to be inevitable as selective pressure for survival . The most alarming is the prevalence of resistance even in the last resort antibiotic colistin  that has added a serious challenge to the current antibiotic crisis.
The cost of developing a prescription drug estimated by The Tufts Center for the Study of Drug Development as published in the Journal of Health Economics in March 2016, is a massive 2.558 billion dollars . Huge research costs and numerous failures in various stages of drug development have lowered the interests of commercial pharmaceutical companies in drug discovery research. The rapid increase in drug resistance among pathogens and the excessive time and cost parameters required to develop a drug demand a robust and faster method of drug discovery. This is where computational strategies come into play, efficiently assisting drug discovery and development with the available in vitro techniques .
Computer-aided drug design (CADD) approaches were applied in this study to find the probable drug targets and discover potential lead candidates against these. SAM is a critical metabolite involved in several biochemical reactions in bacteria. MetK, a SAM producer and various SAM utilizers including DNA adenosine methylase (Dam), Uroporphyrinogen-III methyltransferase (CobA), and tRNA (guanine-N (1)-)-methyltransferase (TrmD) were taken as drug targets in this study. Dam is responsible for DNA replication and mRNA transcription which methylates adenine in DNA of bacteria in contrary to human cytosine. TrmD is responsible for proper reading of codons that prevents + 1 frameshift reading error thus involved in proper peptide elongation. CobA is responsible for corrin ring contraction in vitamin B12 synthesis, an important cofactor for membrane biosynthesis. Thus, all the genes/proteins involved from DNA replication to peptide elongation, and even membrane biosynthesis were targeted to discover new lead candidates, simultaneously preventing easy resistance buildup in these targets.
Selection of MDR strains and obtaining their genomic sequences
Nine prioritized pathogens by WHO  as ‘critical’ and ‘high’ against whom new antimicrobials are sought were taken as the reference organisms. The whole-genome sequences of these organisms published in NCBI were taken for whole-genome alignment. The genomic sequences of the 9 selected pathogens were downloaded from NCBI FTP site in the annotated gbk format.
Multiple sequence alignment (MSA)
MSA was performed using the progressive Mauve algorithm in MAUVE, a multiple sequence alignment software. The genomic regions common to all the aligned sequences were searched for, in MAUVE via visual observation of locally collinear blocks (LCBs) denoted by certain color codes. LCBs represent homologous regions of sequence shared by two or more of the genomes under study without rearrangement .
Alignment of amino acid sequences of the lead proteins
Clustal Omega was used to align the amino acid sequences of S-adenosyl methionine synthase (MetK) for all 9 selected pathogens including Acinetobacter baumannii strain 1656–2 , Campylobacter coli 15–537,360 , Enterococcus faecium Aus0004 , Helicobacter pylori 2017 , Klebsiella pneumoniae subsp. pneumoniae 1084 , Neisseria gonorrhoeae FA 1090 (https://www.ncbi.nlm.nih.gov/nuccore/AE004969) P. aeruginosa B136–33 (https://www.ncbi.nlm.nih.gov/nuccore/NC_020912.1) Staphylococcus aureus 04–02981 (https://www.ncbi.nlm.nih.gov/nuccore/NC_017340.1) and Salmonella typhimurium .
Gene essentiality analysis
The common genes obtained from MAUVE alignment were looked for their essentiality in DEG and OGEE, databases of essential genes.
Obtaining the dockable crystal structures of the target proteins
The X-ray diffraction structures of S-adenosyl methionine synthase, MetK from N. gonorrhoeae (PDB id: 5T8S) ; cobA from P. aeruginosa (PDB id: 2YBQ)  and that of TrmD from P. aeruginosa (PDB id: 5WYQ)  were obtained from Protein Data Bank. For those whose crystal structures were not available in RCSB PDB, homology modeling tools including Phyre2, RaptorX, ps2v2, Swiss-model, and CPHmodel were used to predict their tertiary structures and the best structures were selected based upon the completeness and Z-scores of the predicted structures using Prosa-server.
Preparation of ligand database
In the present work, both ligand-based (LBVS) and structure-based virtual screening (SBVS) was performed. LBVS was done because similar compounds exhibit similar Physico-chemical and biological properties so a broad chemical database with structural diversity would offer an ideal solution for effective lead discovery. In this study, a ligand database containing 715 indole derivatives including marine indoles  was prepared from ZINC database .
Protein and ligand preparation
SBVS was performed based on the common gene in all nine pathogens, MetK, and the metabolite that it produces, SAM which is further utilized in methylation reactions. Prior to molecular docking, the proteins and ligands were prepared for efficient and more accurate docking results. Protein preparation was done by deleting water, adding hydrogen atoms, merging non-polar bonds, and computing Gasteiger charges in AutoDockTools (http://mgltools.scripps.edu/). Similarly, ligand preparation was done in Openbabel GUI  available in PyRx interface by adding hydrogens, energy minimization and converted to pdbqt file format, a useable file format for docking afterward.
Setting reference values for docking
The native ligands were removed from each of the target proteins in Discovery Studio Visualizer 2017 and docked back in their binding sites, a process called re-docking. The binding energy thus calculated was taken as a reference value for identifying potential leads based on their binding energy in the respective binding pockets of the target proteins.
Binding sites prediction
The ligand-binding sites in the target proteins were identified from the RCSB protein data bank. For the homology modeled 3D-structures, 3DLigandSite, a web server that superimposes the ligands bound to the structures similar to the query and thus predicts the binding site , was used to predict the ligand-binding sites. All the amino acid residues around the binding site (Supplementary Table 2) are marked to create a site for molecular docking.
Molecular docking, rescoring and clustering analysis of docked poses
Docking was carried out using AutoDockVina in a virtual screening software, PyRx against the target proteins with the selected ligand database. The conformation with the lowest docked energy was chosen after the docking interactions since, the higher the negative binding energy value, the stronger is the binding of the ligand in the target .
The rescoring of docked poses was done by using the python implementation of NNScore 1.0  in combination with a consensus of the top 24 scoring networks.
AuPosSOM (Automatic analysis of Poses using SOM)  was used for the clustering analysis and to differentiate active compounds from inactive ones. AuPosSOM default parameters were used. The tree was visualized using PhyloWidget .
Protein-ligand interaction visualization
The 2D and 3D protein-ligand interaction for the lead compound was observed and analyzed using Discovery Studio Visualizer 2017 and ligplot+.
The toxic profiles and drug-likeness based on the binding energies were predicted using the OSIRIS program . OSIRIS calculates various drug relevant properties like molecular weight, cLogP, cLogS, Druglikeness, and toxicities like mutagenicity, tumorigenicity, reproductive effects and irritant effects in the lead molecules based on functional groups present in their structures .
Drug target identification
The MAUVE result showed MetK as one of the probable therapeutic targets and was taken as a reference on our search for other therapeutic targets (Fig. 1).
Manual curing of the alignment revealed twenty-four common genes found in most of the strains with diverse roles. Most of the genes (Supplementary Table 1) were ribosomal proteins (14 proteins) and some were involved in ATP synthesis (3 proteins), DNA directed RNA polymerase (2 proteins), chaperones (2 proteins), elongation factor (1 protein), protein translocator (1 protein), involved in thiol assimilation (1 protein) which were not pursued further due to lack of required computational resources to work on these.
Amino acid alignment
The active binding sites in MetK were found same (conserved) in all the pathogens under study upon amino acid sequence alignment using clustal Omega (Fig. 2) making this a common target in all these pathogens. Further analysis of gene essentiality for pathogen survival could reveal whether this is the lead target protein.
Search in DEG and OGEE for gene essentiality of the genes under study in the target organisms revealed metK as essential in H. pylori, P. aeruginosa, S. aureus and S. typhi. Similarly, dam was found to be essential in S. enterica subsp. enterica serovar Typhimurium; cobA was reported as non-essential in P. aeruginosa; and trmD was essential in P. aeruginosa, S. aureus, S. typhi whereas non-essential in H. pylori and Acinetobacter sp. However, these databases do not make use of the interrelation between the genes to record gene essentiality. So the genes mentioned as non-essential here could still be essential when the activity of one is inhibited. Since our works were concerned with multiple targets simultaneously and all these genes under study are SAM utilizers, they were taken for further study despite being non-essential in some instances.
Ligand database preparation
The close structural proximity of the indole ring to the adenosyl moiety of SAM (Fig. 3) pushes indole derivatives to be probable candidates against SAM binding pocket of MetK. Thus, indole derivatives were presumed as the potential ligand sources for virtual screening. A total of 715 indole derivatives were taken from ZINC database among which only 102 showed the drug-like properties based on ADMET parameters (Fig. 4) and used for molecular docking studies.
Molecular docking results
One hundred two ligands that passed ADME/T tests were subjected to molecular docking against the MetK protein of N. gonorrhoeae (PDB id: 5T8S). Fifty three exhibited higher binding energy in the SAM binding pocket of MetK than its native ligand SAM (Supplementary Table 3). Top 3 ligands ZINC04899565 (B.E = − 11.0 Kcal/mol), ZINC06096559 (B.E = − 10.9 Kcal/mol) and ZINC19909549 (B.E = − 10.9 Kcal/mol) with higher B. E than that of SAM (B.E = − 8.7 Kcal/mol) were further investigated to look for their binding modes into the catalytic site of MetK. The common aminoacid residues involved for these three top ligands were Lys274, Gly120, Ile103, Ile307, Phe235, Lys168, Ser191, Asp121, and Asp243 which could have contributed in stronger binding affinities in the binding pocket of MetK (Fig. 5).
NNScore, a neural network based scoring function was then used to re-rank the small-molecule ligands which resulted in 12 positive hits (good binders) as potential inhibitors of MetK in N. gonorrhoeae (Supplementary Table 3). ZINC04899565 was still among the top binders.
We further used a contact activity relationship (CAR) analysis to overcome the limitations of the scoring functions used for docking. Aupossom analyses all the docking poses in multiple conformations given by the docking algorithm and discriminates active and inactive compounds using only mean protein contacts’ footprints calculation. The 12 ligands along with SAM were clustered into 10 different groups with varied scores. The score is determined by the combination between contact specificity and contact intensity of the ligands with various atoms of the protein molecule. The plot (Fig. 6) shows the ligands in leaves 0, 3, 4 and 5 as the most active ones. ZINC01494627 from cluster 0, ZINC14824027 and ZINC49171024 from cluster 3, ZINC04899565 from cluster 4 and ZINC15219763 from cluster 5 can thus be concluded as the potential MetK inhibitors of N. gonorrhoeae.
CAR results showed that these 12 potential leads were distributed in 10 different clusters (0 to 9) with different protein-contact footprints represented in the clustering tree (Fig. 7). The different clusters indicate the differences in the interacting residues with the protein. Cluster 0 contains 1 compound along with SAM (native ligand) which interacts predominantly with A41, E56, Z99, K274, D166 and I237 residues. Cluster 5 contains 1 compound that interacts predominantly with P16, D166, I237 and G242 residues and so on. The similar binding residues of the ligands with that of SAM could have contributed to their inhibition potential of MetK. CAR analysis can thus be used to cluster the compounds as per the binding residues of the protein.
Those 5 screened ligands on further docking against the SAM binding pocket of MetK of all other pathogens under study resulted in 3 potential leads (Table 1). These 3 lead molecules were further docked against all other protein targets CobA, Dam and TrmD (Tables 2, 3, 4) to assess their inhibition potential in multiple targets resulting 2 potential lead candidates (Table 5).
Cross-reactivity with human homologs
Questions on cross-reactivity with human homologs of S-adenosylmethionine synthase (Uniprot id: Q00266) could be raised as it has more than 50% structural similarity with that of bacteria (Table 6). Unfortunately, both leads were potential inhibitors of its human homolog as well (Supplementary Table 4).
Nevertheless, MetK inhibitors of humans could still be used as antibacterial therapeutics because of the presence of SAM transporters in humans  and the SAM requirements can be replenished from external sources. The lack of crystal structures of SAM transporters in humans constrained the molecular docking studies of possible inhibition of the transport system.
Human homologs for other target proteins were not considered because of their dissimilarity with humans.
Drug target identification
Lead molecules with multiple target proteins in a single pathogen and with a common target in multiple pathogens are highly sought  to develop a broad-spectrum drug. This strategy has been designed for preventing easy resistance development against these new drugs. Developing resistance in multiple targets at once could be evolutionarily challenging for any pathogens and probably impossible for the bacteria to survive against such developed drugs. Hence, genome-level sequence alignment of major pathogens could give common new lead target proteins for the screening of lead inhibitor molecules of these proteins, MetK being the common target in this study (Fig. 1). From this, the respective metabolic pathway or other target proteins could be identified based on protein-protein interactions.
MetK codes the formation of SAM from ATP and methionine as substrates. SAM is utilized by three major metabolic pathways: transmethylation, transsulfuration, and polyamine synthesis making SAM an important molecule in normal cell functioning and survival . In addition, SAM is a primary methyl donor in multiple reactions including corrin ring methylation , RNA methylation , and DNA methylation  thus, these steps were also taken as lead targets.
Quorum sensing is one of the major causes of resistance in pathogens that utilize autoinducers which inturn utilize SAM as a substrate . Also, it controls biofilm formation and virulence in bacteria . The literature search further verified the essentiality of MetK in many different pathogens [34, 35]. Also, the lack of reports about SAM transporters in any of the mentioned target pathogens makes this a better target. So, including MetK, other SAM utilizing proteins, namely CobA, Dam, and TrmD were taken as potential drug targets against which virtual screening of compounds was done.
Indoles as potential drugs
The adenosyl moiety of SAM and ATP binding domain present in kinases  probably suggests kinase inhibitors as potential inhibitors of these proteins that biosynthesize or utilize SAM. Protein kinase inhibitors represent an important class of targeted therapeutic agents, particularly as anticancer drugs . Several indoles have been reported to possess kinase inhibition potentials [38,39,40,41]. Also, Indole is reported to be metabolized by human cytochrome P450 2A6  and the source of indole could be from tryptophan metabolism by gut microflora. This indicates that indole could be easily transported in humans through gut suggesting indoles are metabolized in humans thus indicating these could not pose toxicity . In addition, indole has been suggested to be pharmaceutical scaffolds for drug development . Also, its metabolized derivative Indirubin has been reported to be kinase inhibitor .
The primary reason for lead molecules not being able to pass the clinical trials is their inability to reach the target and perform its predicted function, and also the toxicity issues . Thus, ADMET and pharmacokinetic properties evaluation in the early stages of drug discovery seem to be a wiser choice. The toxic profiles and drug-likeness were predicted by OSIRIS using various parameters. Various physicochemical and drug relevant properties such as Molecular weight, cLogP, cLogS, number of hydrogen bond donors and acceptors, topological polar surface area, number of rotatable bonds and drug-likeness were analyzed for each of the lead molecules. The parameters (Fig. 4) were set based on Lipinski’s rule of five which predicts whether a chemical compound has chemical and physical properties that would make it likely to be an orally active drug in humans. The rule states that most orally active drugs will have molecular weight ≤ 500, logP ≤5, hydrogen bond donors ≤5, and hydrogen bond acceptors ≤10. Also, the aqueous solubility of a compound measured as logS significantly affects its absorption and distribution characteristics. Low solubility usually goes along with a bad absorption so poorly soluble compounds should be avoided. The number of rotatable bonds determines the flexibility of compounds and can predict oral bioavailability. It has been reported that 10 or fewer rotatable bonds in a molecule indicated good oral bioavailability . Polar surface area (PSA) of a molecule can predict membrane permeability including crossing of the blood-brain barrier. Most known drugs have PSA values less than 120 Å .
Apart from indole, ZINC04899565 has a benzene and a 2,5-diketopiperazine ring. Out of all the naturally occurring peptide antibiotics, the 2,5-diketopiperazine rings are among the most numerous. Cycloserine diketopiperazine active against Mycobacterium tuberculosis, bicyclomycin active against gram-negative bacteria, avrainvillamide which contains a 3-alkylidene-3H-indole-1-oxide function active even against multi drug-resistant bacteria are some instances . Moreover, 2,5-diketopiperazines have also been reported to inhibit quorum sensing in certain gram-negative bacteria thus blocking cell-to-cell communication and restraining the virulence as well . Structurally, the 2,5-diketopiperazines are the smallest possible cyclic peptides, which are peptidoimmetic in nature resembling a constrained protein beta turn. They have two cis-amide bonds thus possessing 2H-bond acceptors and donors each. Although they contain conformationally constrained heterocyclic scaffolds, they are flexible since the six-membered ring can exist either in a flat or a slight puckered boat conformation. Moreover, these are stable to proteolysis . All these features support them bind to a wide range of enzymes and receptors, and their good bioavailability and resistance to enzymatic degradation make them excellent drug candidates. Thus, ZINC04899565 has the potential to be a broad-spectrum antimicrobial based on this study and could be pursued further.
ZINC49171024 is an indole with a pyrrolidine and a benzimidazole ring. Pyrrolidine moiety containing compounds have been reported as antimicrobials and fungals . The strong lipophilic properties of benzimidazoles contribute in producing antimicrobial effects . All these features indicate the high probability of these two compounds to be developed as broad-spectrum antimicrobials.
CADD approaches were used in this study to discover potential methyltransferase inhibitory activities of indole derivatives. Multiple protein targets were subjected to molecular docking that resulted ZINC04899565 and ZINC49171024 as probable therapeutic drug candidates with multi-target potential and probable antimicrobial resistance managers since these multiple targets would be troublesome for the pathogens to easily evolve resistance as these are really critical in its survival and mutating these targets could be more lethal than survival.
Availability of data and materials
The 3D structures of the proteins and ligands used in the current study are available from RCSB PDB (https://www.rcsb.org/) and ZINC (http://zinc15.docking.org/) databases respectively. Molecular docking can be performed using PyRx (https://pyrx.sourceforge.io/). The library of screened ligands, 3D structures of homology modeled proteins and any other information regarding the current work will be made available by the corresponding author on reasonable request.
Multi Drug Resistant
World Health Organization
Absorption, Digestion, Metabolism, Excretion and Toxicity
Computer Aided Drug Design
Database of Essential Genes
Online GEne Essentiality database
Contact Activity Relationship
Tello A, Austin B, Telfer TC. Selective pressure of antibiotic pollution on bacteria of importance to public health. Environ Health Perspect. 2012;120(8):1100–6. https://doi.org/10.1289/ehp.1104650.
Caniaux I, van Belkum A, Zambardi G, Poirel L, Gros MF. MCR: modern colistin resistance. Eur J Clin Microbiol Infect Dis, 2017;36:415–420. https://doi.org/10.1007/s10096-016-2846-y.
University, T. Tufts CSDD assessment of cost to develop and win marketing approval for a new drug now published. In: Retrieved from development. Tufts: Center for the Study of Drug website; 2016. http://csdd.tufts.edu/news/complete_story/tufts_csdd_rd_cost_study_now_published.
Sliwoski G, Kothiwale S, Meiler J, Lowe EW. Computational methods in drug discovery. Pharmacol Rev. 2014;66(1):334–95. https://doi.org/10.1124/pr.112.007336.
Lawe-Davies O, Bennett S. WHO publishes list of bacteria for which new antibiotics are urgently needed: WHO; 2017. p. 1. Retrieved from http://www.who.int/mediacentre/news/releases/2017/bacteria-antibiotics-needed/en/.
Darling ACE, Mau B, Blattner FR, Perna NT. Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 2004;14(7):1394–403. https://doi.org/10.1101/gr.2289704.
Park JY, Kim S, Kim SM, Cha SH, Lim SK, Kim J. Complete genome sequence of multidrug-resistant Acinetobacter baumannii strain 1656-2 Which Forms Sturdy Biofilm. J Bacteriology. 2011;193:6393–4. https://doi.org/10.1128/JB.06109-11.
Pearson BM, Rokney A, Crossman LC, Miller WG, Wain J, van Vliet AHM. Complete genome sequence of the Campylobacter coli clinical isolate 15-537360. Genome Announcements. 2013;1(6):e01056–13. https://doi.org/10.1371/journal.
Lam MMC, Seemann T, Bulach DM, Gladman SL, Chen H, Haring V, Stineara TP. Comparative analysis of the first complete Enterococcus faecium genome. J Bacteriol. 2012;194(9):2334–41. https://doi.org/10.1128/JB.00259-12.
Avasthi TS, Devi SH, Taylor TD, Kumar N, Baddam R, Kondo S, Ahmed N. Genomes of two chronological isolates (Helicobacter pylori 2017 and 2018) of the west African Helicobacter pylori strain 908 obtained from a single patient. J Bacteriol. 2011;193:3385–6. https://doi.org/10.1128/JB.05006-11.
Lin AC, Liao TL, Lin YC, Lai YC, Lu MC, Chen YT. Complete genome sequence of Klebsiella pneumoniae 1084, a hypermucoviscosity-negative K1 clinical strain. J Bacteriol. 2012;194:6316. https://doi.org/10.1128/JB.01548-12.
Silva C, Calva E, Puente JL, Zaidi MB, Vinuesa P. Complete genome sequence of Salmonella enterica Serovar Typhimurium strain YU15 (sequence type 19) harboring the Salmonella Genomic Island 1 and virulence plasmid pSTV. Genome Announcements. 2016;4(2):e00252–16. https://doi.org/10.1128/genomeA.00252-16.
Dranow DM, Delker SL, Lorimer DD, Edwards TE. Crystal structure of a S-adenosylmethionine synthase from Neisseria gonorrhoeae with bound S-adenosylmethionine. AMP, Pyrophosphate, Phosphate, and Magnesium. n.d.. https://doi.org/10.2210/PDB5T8S/PDB.
Storbeck, S, Saha, S, Krausze, J, Klink, BU, Heinz, DW, Layer, G. (2011). Crystal structure of the heme d 1 biosynthesis enzyme NirE in complex with its substrate reveals new insights into the catalytic mechanism of S-adenosyl-L-methionine-dependent uroporphyrinogen III methyltransferases. J Biol Chem https://doi.org/10.1074/jbc. M111.239855.
Jaroensuk, J, Wong, YH, Zhong, W, Liew, CW, Maenpuen, S, Sahili, AE, Fuangthong, M. (2019). Crystal structure and catalytic mechanism of the essential m 1 G37 tRNA methyltransferase TrmD from Pseudomonas aeruginosa. https://doi.org/10.1261/rna.066746.118.
Netz N, Opatz T. Marine indole alkaloids. Marine Drugs. 2015;13(8):4814–914. https://doi.org/10.3390/md13084814.
Irwin JJ, Shoichet BK. ZINC - a free database of commercially available compounds for virtual screening. J Chem Inf Model. 2005;45(1):177–82. https://doi.org/10.1021/ci049714+.
O’Boyle NM, Banck M, James CA, Morley C, Vandermeersch T, Hutchison GR. Open babel: an open chemical toolbox. Journal of Cheminformatics. 2011;3(10). https://doi.org/10.1186/1758-2946-3-33.
Wass MN, Kelley LA, Sternberg MJE. 3DLigandSite: predicting ligand-binding sites using similar structures. Nucleic Acids Res. 2010;38(SUPPL. 2). https://doi.org/10.1093/nar/gkq406.
Dallakyan S, Olson AJ. Small-molecule library screening by docking with PyRx. Methods Mol Biol. 2015;1263:243–50. https://doi.org/10.1007/978-1-4939-2269-7_19.
Durrant JD, McCammon JA. NNScore: a neural-network-based scoring function for the characterization of protein-ligand complexes. J Chem Inf Model. 2010. https://doi.org/10.1021/ci100244v.
Bouvier G, Evrard-todeschi N, Girault JP, Bertho G. Automatic clustering of docking poses in virtual screening process using self-organizing map. Bioinformatics. 2009. https://doi.org/10.1093/bioinformatics/btp623.
Jordan GE, Piel WH. PhyloWidget: web-based visualizations for the tree of life. Bioinformatics. 2008. https://doi.org/10.1093/bioinformatics/btn235.
Nisha CM, Kumar A, Nair P, Gupta N, Silakari C, Tripathi T, Kumar A. Molecular docking and in silico admet study reveals acylguanidine 7a as a potential inhibitor of β -secretase. Adv Bioinforma. 2016;2016. https://doi.org/10.1155/2016/9258578.
Sander T, Freyss J, Von Korff M, Rufener C. DataWarrior: an open-source program for chemistry aware data visualization and analysis. J Chem Inf Model. 2015;55(2):460–73. https://doi.org/10.1021/ci500588j.
Agrimi G, Di Noia MA, Marobbio CMT, Fiermonte G, Lasorsa FM, Palmieri F. Identification of the human mitochondrial S-adenosylmethionine transporter: bacterial expression, reconstitution, functional characterization and tissue distribution. Biochem J. 2004;379:183–90. https://doi.org/10.1042/bj20031664.
Bolognesi L, M. Polypharmacology in a single drug: multitarget drugs. Curr Med Chem. 2013;20(13):1639–45. https://doi.org/10.2174/0929867311320130004.
Lu, S. C. (2000). S-Adenosylmethionine. Int J Biochem Cell Biol, 32(4), 391–395. https://doi.org/10.1016/S1357-2725(99)00139-9.
Roper, J. M., Raux, E., Brindley, A. A., Schubert, H. L., Gharbia, S. E., Shah, H. N., & Warren, M. J. (2000). The enigma of cobalamin (vitamin B12) biosynthesis in Porphyromonas gingivalis: identification and characterization of a functional corrin pathway. J Biol Chem, 275(51), 40316–40323. https://doi.org/10.1074/jbc. M007146200.
Motorin Y, Helm M. RNA nucleotide methylation. Wiley Interdisciplinary Reviews: RNA. 2011;2:611–31. https://doi.org/10.1002/wrna.79.
Sánchez-Romero MA, Cota I, Casadesús J. DNA methylation in bacteria: from the methyl group to the methylome. Curr Opin Microbiol. 2015;25:9–16. https://doi.org/10.1016/j.mib.2015.03.004.
Wei Y, Perez LJ, Ng WL, Semmelhack MF, Bassler BL. Mechanism of Vibrio cholerae autoinducer-1 biosynthesis. ACS Chem Biol. 2011;6(4):356–65. https://doi.org/10.1021/cb1003652.
Rutherford ST, Bassler BL. Bacterial quorum sensing: its role in virulence and possibilities for its control. Cold Spring Harbor Perspectives in Medicine, Vol. 2012;2. https://doi.org/10.1101/cshperspect.a012427.
Walker DJF, Heap JT, Winzer K, Minton NP. A genetic assay for gene essentiality in Clostridium. Anaerobe. 2016;42:40–3. https://doi.org/10.1016/j.anaerobe.2016.07.007.
Wei Y, Newman EB. Studies on the role of the metK gene product of Escherichia coli K-12. Mol Microbiol. 2002;43(6):1651–6. https://doi.org/10.1046/j.1365-2958.2002.02856.x.
Villamor, J G, Kaschani, F, Colby, T, Oeljeklaus, J, Zhao, D, Kaiser, M, van der Hoorn, RAL (2013). Profiling protein kinases and other ATP binding proteins in Arabidopsis using acyl-ATP probes. Mol Cell Proteomics, 12(9), 2481–2496. https://doi.org/10.1074/mcp. M112.026278.
Grant SK. Therapeutic protein kinase inhibitors. Cell Mol Life Sci. 2009;66:1163–77. https://doi.org/10.1007/s00018-008-8539-7.
Chowdhury S, Sessions EH, Pocas JR, Grant W, Schröter T, Lin L, Feng Y. Discovery and optimization of indoles and 7-azaindoles as rho kinase (ROCK) inhibitors (part-I). Bioorganic and Medicinal Chemistry Letters. 2011. https://doi.org/10.1016/j.bmcl.2011.09.083.
Kiliç Z, Işgör YG, Olgen S. Evaluation of new indole and bromoindole derivatives as pp60(c-Src) tyrosine kinase inhibitors. Chem Biol Drug Des. 2009;74(4):397–404. https://doi.org/10.1111/j.1747-0285.2009.00876.x.
Ölgen S, Kiliç-Kurt Z, Şener F, Işgör YG, Çoban T. Evaluation of novel aminomethyl indole derivatives as Src kinase inhibitors and antioxidant agents. Chemotherapy. 2011;57(1):1–6. https://doi.org/10.1159/000317764.
Polychronopoulos P, Magiatis P, Skaltsounis AL, Myrianthopoulos V, Mikros E, Tarricone A, Meijer L. Structural basis for the synthesis of Indirubins as potent and selective inhibitors of glycogen synthase Kinase-3 and Cyclin-dependent kinases. J Med Chem. 2004. https://doi.org/10.1021/jm031016d.
Gillam EMJ, Notley LM, Cai H, De Voss JJ, Guengerich FP. Oxidation of indole by cytochrome P450 enzymes. Biochemistry. 2000;39(45):13817–24. https://doi.org/10.1021/bi001229u.
Banoglu E, Jha GG, King RS. Hepatic microsomal metabolism of indole to indoxyl, a precursor of indoxyl sulfate. Eur J Drug Metab Pharmacokinet. 2001;26(4):235–40. https://doi.org/10.1007/BF03226377.
Suenkel B, Fischer F, Steegborn C. Inhibition of the human deacylase Sirtuin 5 by the indole GW5074. Bioorganic and Medicinal Chemistry Letters. 2013;23(1):143–6. https://doi.org/10.1016/j.bmcl.2012.10.136.
Wu ZL, Aryal P, Lozach O, Meijer L, Guengerich FP. Biosynthesis of new indigoid inhibitors of protein kinases using recombinant cytochrome P450 2A6. Chem Biodivers. 2005;2(1):51–65. https://doi.org/10.1002/cbdv.200490166.
Hughes JP, Rees S, Kalindjian SB, Philpott KL. Principles of early drug discovery. Br J Pharmacol. 2011;162(6):1239–49. https://doi.org/10.1111/j.1476-5381.2010.01127.x.
Veber DF, Johnson SR, Cheng HY, Smith BR, Ward KW, Kopple KD. Molecular properties that influence the oral bioavailability of drug candidates. J Med Chem. 2002. https://doi.org/10.1021/jm020017n.
Ertl P, Rohde B, Selzer P. Fast calculation of molecular polar surface area as a sum of fragment-based contributions and its application to the prediction of drug transport properties. J Med Chem. 2000. https://doi.org/10.1021/jm000942e.
Borthwick AD. 2,5-diketopiperazines: synthesis, reactions, medicinal chemistry, and bioactive natural products. Chem Rev. 2012. https://doi.org/10.1021/cr200398y.
Scoffone VC, Chiarelli LR, Makarov V, Brackman G, Israyilova A, Azzalin A, Buroni S. Discovery of new diketopiperazines inhibiting Burkholderia cenocepacia quorum sensing in vitro and in vivo. Sci Rep. 2016. https://doi.org/10.1038/srep32487.
Muneer S, Memon S, Pahnwar QK, Bhatti AA, Khokhar TS. Synthesis and investigation of antimicrobial properties of pyrrolidine appended calix  arene. Journal of Analytical Science and Technology. 2017. https://doi.org/10.1186/s40543-017-0111-3.
Alasmary FAS, Snelling AM, Zain ME, Alafeefy AM, Awaad AS, Karodia N. Synthesis and evaluation of selected benzimidazole derivatives as potential antimicrobial agents. Molecules. 2015. https://doi.org/10.3390/molecules200815206.
The authors would like to acknowledge Central Department of Biotechnology, Tribhuvan University for providing space to perform lab works.
Not applicable (Self-funded by the authors).
Ethics approval and consent to participate
Consent for publication
“The authors declare that they have no competing interests”.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Table 1. Putative lead target proteins identified as common, by whole genome sequence alignment of WHO priority list pathogens  and M. tuberculosis.
Supplementary Table 2. Predicted binding sites in MetK of Helicobacter pylori in 3DligandSite predicted via Phyre2 webserver.
Supplementary Table 3. 53 indole derivatives with higher B.E in MetK of Neisseria gonorrhoeae than SAM with their NNScore values.
Supplementary Table 4. Docking results of two lead molecules in human homolog of MetK, Mat2B (RCSB id: 4KTT).
About this article
Cite this article
Tha, S., Shakya, S., Malla, R. et al. Prospects of Indole derivatives as methyl transfer inhibitors: antimicrobial resistance managers. BMC Pharmacol Toxicol 21, 33 (2020). https://doi.org/10.1186/s40360-020-00402-9