High-Throughput Virtual Screening of Compounds with Electrophilic Fragments for New Potential Covalent Inhibitors of Bacterial Proteins †

: The search for new antibacterial drugs has continued to be an urgent matter. One of the approaches is the development of covalent inhibitors using biochemoinformatics at the initial stages. In this work, structures of a few plant-derived substances with electrophilic unsaturated carbonyl and structures of small synthetic compounds suitable for fragment-based drug discovery (FBDD) with -CH 2 -Br group were selected as ligands for sets of structures of bacterial proteins. The theoretical assessment was carried out using the Autodock Vina program for calculation and FYTdock for the organization of the process and the analysis of results. Natural Ixerine D as well as synthetic 4-(4-(2-bromoethyl)piperazin-1-yl)-7-nitrobenzofurazan demonstrated the most promising results as potential Cys-targeted inhibitors.


Introduction
Natural antibiotics, their derivatives, and synthetic antimicrobials are the primary tools to treat bacterial infections. A number of bacteria have developed resistance to many or even all such currently available drugs, hindering the treatment of these diseases. Therefore, the search for new antibacterial drugs does not cease to be an urgent scientific task. One of the approaches to creating therapeutic agents is the development of covalent inhibitors. The initial stages of any modern drug design company include the use of modern methods of biochemoinformatics, in particular, molecular docking, i.e., computing of ligand-protein complexes with an assessment of their geometry and affinity [1][2][3].
In this work,~20 structures of plant-derived electrophilic substances from the Pubchem database were selected as ligands as well as some structures suitable for the fragmentbased ligand design approach (FBDD-Fragment-based drug discovery) containing the electrophilic fragment -CH 2 -Br. The electrophilic nature of the phytochemicals and the fragments provide a possibility of covalent medication of nucleophilic atoms of Cys and His residues in proteins. In this work, the possibility was additionally evaluated in silico using the Autodock Vina program for docking simulations and FYTdock [4] to organize, run, and analyze the docking results.
Approximately 2900 protein structures of Mycobacterium tuberculosis and 500 random protein structures of some other bacterial species were selected to create a library of bacterial protein structures from the Protein Data Bank. Docking results were initially collected and processed using FYTdock software as an Excel spreadsheet showing binding energies, amino acid environment, protein ligand-amino acid interactions, protein-ligand complexes.
The result was taken into account if the value of E bind was no more than −6.0 kcal/mol and the distance from the electrophilic fragment of the ligand structure to the sulfur atom of the thiol group of cysteine residues in the protein-ligand complexes obtained in silico did not exceed each 0.45 nm (distance criterion). For a graphical representation of the result, the Biovia program was used.

Results
Fortunately, it was found that Ixerin D (Pubchem database number CID101553163) from dandelion, a common and widely grown plant, demonstrated a number of interactions with high affinity and the location of its electrophilic fragment within 0.4 nm from the sulfur atom of the cysteine of Mtb proteins, lipoyl synthase, inosine monophosphate dehydrogenase, and beta-ketoacyl-acyl carrier protein synthase III from Mycobacterium tuberculosis ( Table 1 and Figure 1).  This compound is a metabolite of the common dandelion (Taraxacum officinale) and is probably of low toxicity to humans due to the use of parts of this plant as food or medicine by humans and some animals. The beta-ketoacyl-acyl carrier protein synthase III is very important for fatty acid biosynthesis and for the normal life cycle of Mtb [5]. Such calculated and theoretical data indicate the possibility of a favorable outcome of the This compound is a metabolite of the common dandelion (Taraxacum officinale) and is probably of low toxicity to humans due to the use of parts of this plant as food or medicine by humans and some animals. The beta-ketoacyl-acyl carrier protein synthase III is very important for fatty acid biosynthesis and for the normal life cycle of Mtb [5]. Such calculated and theoretical data indicate the possibility of a favorable outcome of the biological testing of Ixerin D, and it can be obtained from a natural source, which does not make it necessary to develop a scheme for its chemical synthesis.
For the synthetic ligand 4-(4-(2-bromoethyl)piperazin-1-yl)-7-nitrobenzofurazan, compiled using the FBDD approach, the -CH 2 -Br fragment was found to be located close to the cysteine sulfur atom in Sortase B from Staphylococcus aureus, a human pathogen [6], E. coli Gsp amidase, which regulates the redox state of E. coli cells [7], β-lactamase S70C BlaC from Mycobacterium tuberculosis, which contributes to the development of the bacteria natural resistance to β-lactam antibiotics [8] (Table 2). It is important to note that, despite the localization of the electrophilic fragment of the ligand near cysteine, definitely these ligand-receptor interactions can be hindered due to geometrical features and energetically favorable location of the ligand in the protein. The 4-(4-(2-bromoethyl)piperazin-1-yl)-7-nitrobenzofurazan showed interactions with high affinity and favorable location of the -CH 2 -Br fragment for covalent binding with the following bacterial proteins: E. coli bifunctional glutathionylspermidine synthetase/amidase and E. coli K-12 methionine aminopeptidase with binding energies of −8.3 kcal/mol and −7.0 kcal/mol, respectively (Table 3 and Figure 2). These enzymes are potential new drug targets, and inhibitors of these enzymes may be useful as prototypes of new antibacterial agents [9,10]. Table 3. The proteins from protein-ligand complexes where an electrophilic fragment of 4-(4-(2bromoethyl)piperazin-1-yl)-7-nitrobenzofurazan located within 0.4 nm from the sulfur atom of a cysteine residue and their binding energies. The compound 2-bromo-1-(4-(nitrobenzofurazan-4-yl)piperazin-1-yl)ethanone had the orientation of the -CH 2 -Br fragment near the C37L/C151T/C442A histidine of the CYP51 triple mutant from Mycobacterium tuberculosis, Microcin-processing metalloprotease TldD/E from E. coli, CYP134A1 with a closed-loop substrate binding from Bacillus subtilis and other bacterial proteins with binding energies from −6.7 to −9.0 kcal/mol (Table 4 and Figure 3). The compound 2-bromo-1-(4-(nitrobenzofurazan-4-yl)piperazin-1-yl)ethanone had the orientation of the -CH2-Br fragment near the C37L/C151T/C442A histidine of the CYP51 triple mutant from Mycobacterium tuberculosis, Microcin-processing metalloprotease TldD/E from E. coli, CYP134A1 with a closed-loop substrate binding from Bacillus subtilis and other bacterial proteins with binding energies from −6.7 to −9.0 kcal/mol (Table 4 and Figure 3).

Conclusions
Based on in silico molecular docking, natural compound Ixerin D from common dandelion, as well as synthetic ligands, 4-(4-(2-bromoethyl)piperazin-1-yl)-7nitrobenzofurazan, 2-bromo-1-(4-(nitrobenzofurazan-4-yl)piperazin-1-yl)ethenone, and 2-bromo-N-(4-bromophenyl)acetamide, fragments for structures of new covalent molecular tools or drugs, were identified to be able to covalently modified inhibitors of various bacterial proteins with localization of their electrophilic fragments within 0.4 nm from functional amino acid fragments. Residues and binding energy in the range from −6.0 to −10.7 kcal/mol. Thus, the results substantiate perspectives of experimental studies of these ligands as potential antibacterial agents or molecular tools with covalent modifier properties.

Conclusions
Based on in silico molecular docking, natural compound Ixerin D from common dandelion, as well as synthetic ligands, 4-(4-(2-bromoethyl)piperazin-1-yl)-7-nitrobenzofurazan, 2-bromo-1-(4-(nitrobenzofurazan-4-yl)piperazin-1-yl)ethenone, and 2-bromo-N-(4-bromophenyl) acetamide, fragments for structures of new covalent molecular tools or drugs, were identified to be able to covalently modified inhibitors of various bacterial proteins with localization of their electrophilic fragments within 0.4 nm from functional amino acid fragments. Residues and binding energy in the range from −6.0 to −10.7 kcal/mol. Thus, the results substantiate perspectives of experimental studies of these ligands as potential antibacterial agents or molecular tools with covalent modifier properties.