Determination of absolute configuration and binding efficacy of benzimidazole-based FabI inhibitors through the support of electronic circular dichroism and MM-GBSA techniques.

We have previously reported benzimidazole-based compounds to be potent inhibitors of FabI for Francisella tularensis (FtFabI), making them promising antimicrobial hits. Optically active enantiomers exhibit markedly differing affinities toward FtFabI. The IC50 of benzimidazole (-)-1 is ∼100× lower than the (+)-enantiomer, with similar results for the 2 enantiomers. Determining the absolute configuration for these optical compounds and elucidating their binding modes is important for further design. Electronic circular dichroism (ECD) quantum calculations have become important in determining absolute configurations of optical compounds. We determined the absolute configuration of (-)/(+)-1 and (-)/(+)-2 by comparing experimental spectra and theoretical density functional theory (DFT) simulations of ECD spectra at the B3LYP/6-311+G(2d, p) level using Gaussian09. Comparison of experimental and calculated ECD spectra indicates that the S configuration corresponds to the (-)-rotation for both compounds 1 and 2, while the R configuration corresponds to the (+)-rotation. Further, molecular dynamics simulations and MM-GBSA binding energy calculations for these two pairs of enantiomers with FtFabI show much tighter binding MM-GBSA free energies for S-1 and S-2 than for their enantiomers, R-1 and R-2, consistent with the S configuration being the more active one, and with the ECD determination of the S configuration corresponding to (-) and the R configuration corresponding to (+). Thus, our computational studies allow us to assign (-) to (S)- and (+) to (R)- for compounds 1 and 2, and to further evaluate structural changes to improve efficacy.

Evaluating thermodynamic integration performance of the new amber molecular dynamics package and assess potential halogen bonds of enoyl-ACP reductase (FabI) benzimidazole inhibitors.

Thermodynamic integration (TI) can provide accurate binding free energy insights in a lead optimization program, but its high computational expense has limited its usage. In the effort of developing an efficient and accurate TI protocol for FabI inhibitors lead optimization program, we carefully compared TI with different Amber molecular dynamics (MD) engines (sander and pmemd), MD simulation lengths, the number of intermediate states and transformation steps, and the Lennard-Jones and Coulomb Softcore potentials parameters in the one-step TI, using eleven benzimidazole inhibitors in complex with Francisella tularensis enoyl acyl reductase (FtFabI). To our knowledge, this is the first study to extensively test the new AMBER MD engine, pmemd, on TI and compare the parameters of the Softcore potentials in the one-step TI in a protein-ligand binding system. The best performing model, the one-step pmemd TI, using 6 intermediate states and 1 ns MD simulations, provides better agreement with experimental results (RMSD = 0.52 kcal/mol) than the best performing implicit solvent method, QM/MM-GBSA from our previous study (RMSD = 3.00 kcal/mol), while maintaining similar efficiency. Briefly, we show the optimized TI protocol to be highly accurate and affordable for the FtFabI system. This approach can be implemented in a larger scale benzimidazole scaffold lead optimization against FtFabI. Lastly, the TI results here also provide structure-activity relationship insights, and suggest the parahalogen in benzimidazole compounds might form a weak halogen bond with FabI, which is a well-known halogen bond favoring enzyme.

Comparison of radii sets, entropy, QM methods, and sampling on MM-PBSA, MM-GBSA, and QM/MM-GBSA ligand binding energies of F. tularensis enoyl-ACP reductase (FabI).

To validate a method for predicting the binding affinities of FabI inhibitors, three implicit solvent methods, MM-PBSA, MM-GBSA, and QM/MM-GBSA were carefully compared using 16 benzimidazole inhibitors in complex with Francisella tularensis FabI. The data suggests that the prediction results are sensitive to radii sets, GB methods, QM Hamiltonians, sampling protocols, and simulation length, if only one simulation trajectory is used for each ligand. In this case, QM/MM-GBSA using 6 ns MD simulation trajectories together with GB(neck2) , PM3, and the mbondi2 radii set, generate the closest agreement with experimental values (r(2) = 0.88). However, if the three implicit solvent methods are averaged from six 1 ns MD simulations for each ligand (called "multiple independent sampling"), the prediction results are relatively insensitive to all the tested parameters. Moreover, MM/GBSA together with GB(HCT) and mbondi, using 600 frames extracted evenly from six 0.25 ns MD simulations, can also provide accurate prediction to experimental values (r(2) = 0.84). Therefore, the multiple independent sampling method can be more efficient than a single, long simulation method. Since future scaffold expansions may significantly change the benzimidazole's physiochemical properties (charges, etc.) and possibly binding modes, which may affect the sensitivities of various parameters, the relatively insensitive "multiple independent sampling method" may avoid the need of an entirely new validation study. Moreover, due to large fluctuating entropy values, (QM/)MM-P(G)BSA were limited to inhibitors' relative affinity prediction, but not the absolute affinity. The developed protocol will support an ongoing benzimidazole lead optimization program.

Structural and biological evaluation of a novel series of benzimidazole inhibitors of Francisella tularensis enoyl-ACP reductase (FabI).

Francisella tularensis, the causative agent of tularemia, presents a significant biological threat and is a Category A priority pathogen due to its potential for weaponization. The bacterial FASII pathway is a viable target for the development of novel antibacterial agents treating Gram-negative infections. Here we report the advancement of a promising series of benzimidazole FabI (enoyl-ACP reductase) inhibitors to a second-generation using a systematic, structure-guided lead optimization strategy, and the determination of several co-crystal structures that confirm the binding mode of designed inhibitors. These compounds display an improved low nanomolar enzymatic activity as well as promising low microgram/mL antibacterial activity against both F. tularensis and Staphylococcus aureus and its methicillin-resistant strain (MRSA). The improvements in activity accompanying structural modifications lead to a better understanding of the relationship between the chemical structure and biological activity that encompasses both enzymatic and whole-cell activity.

Hit identification and optimization in virtual screening: practical recommendations based on a critical literature analysis.

A critical analysis of virtual screening results published between 2007 and 2011 was performed. The activity of reported hit compounds from over 400 studies was compared to their hit identification criteria. Hit rates and ligand efficiencies were calculated to assist in these analyses, and the results were compared with factors such as the size of the virtual library and the number of compounds tested. A series of promiscuity, druglike, and ADMET filters were applied to the reported hits to assess the quality of compounds reported, and a careful analysis of a subset of the studies that presented hit optimization was performed. These data allowed us to make several practical recommendations with respect to selection of compounds for experimental testing, definition of hit identification criteria, and general virtual screening hit criteria to allow for realistic hit optimization. A key recommendation is the use of size-targeted ligand efficiency values as hit identification criteria.

Discovery of a novel and potent class of F. tularensis enoyl-reductase (FabI) inhibitors by molecular shape and electrostatic matching.

Enoyl-acyl carrier protein (ACP) reductase, FabI, is a key enzyme in the bacterial fatty acid biosynthesis pathway (FAS II). FabI is an NADH-dependent oxidoreductase that acts to reduce enoyl-ACP substrates in a final step of the pathway. The absence of this enzyme in humans makes it an attractive target for the development of new antibacterial agents. FabI is known to be unresponsive to structure-based design efforts due to a high degree of induced fit and a mobile flexible loop encompassing the active site. Here we discuss the development, validation, and careful application of a ligand-based virtual screen used for the identification of novel inhibitors of the Francisella tularensis FabI target. In this study, four known classes of FabI inhibitors were used as templates for virtual screens that involved molecular shape and electrostatic matching. The program ROCS was used to search a high-throughput screening library for compounds that matched any of the four molecular shape queries. Matching compounds were further refined using the program EON, which compares and scores compounds by matching electrostatic properties. Using these techniques, 50 compounds were selected, ordered, and tested. The tested compounds possessed novel chemical scaffolds when compared to the input query compounds. Several hits with low micromolar activity were identified and follow-up scaffold-based searches resulted in the identification of a lead series with submicromolar enzyme inhibition, high ligand efficiency, and a novel scaffold. Additionally, one of the most active compounds showed promising whole-cell antibacterial activity against several Gram-positive and Gram-negative species, including the target pathogen. The results of a preliminary structure-activity relationship analysis are presented.

Characterization of the Escherichia coli ClpY (HslU) substrate recognition site in the ClpYQ (HslUV) protease using the yeast two-hybrid system.

In Escherichia coli, ClpYQ (HslUV) is a two-component ATP-dependent protease in which ClpQ is the peptidase subunit and ClpY is the ATPase and the substrate-binding subunit. The ATP-dependent proteolysis is mediated by substrate recognition in the ClpYQ complex. ClpY has three domains, N, I, and C, and these domains are discrete and exhibit different binding preferences. In vivo, ClpYQ targets SulA, RcsA, RpoH, and TraJ molecules. In this study, ClpY was analyzed to identify the molecular determinants required for the binding of its natural protein substrates. Using yeast two-hybrid analysis, we showed that domain I of ClpY contains the residues responsible for recognition of its natural substrates, while domain C is necessary to engage ClpQ. Moreover, the specific residues that lie in the amino acid (aa) 137 to 150 (loop 1) and aa 175 to 209 (loop 2) double loops in domain I of ClpY were shown to be necessary for natural substrate interaction. Additionally, the two-hybrid system, together with random PCR mutagenesis, allowed the isolation of ClpY mutants that displayed a range of binding activities with SulA, including a mutant with no SulA binding trait. Subsequently, via methyl methanesulfonate tests and cpsB::lacZ assays with, e.g., SulA and RcsA as targets, we concluded that aa 175 to 209 of loop 2 are involved in the tethering of natural substrates, and it is likely that both loops, aa 137 to 150 and aa 175 to 209, of ClpY domain I may assist in the delivery of substrates into the inner core for ultimate degradation by ClpQ.