Ligand-Based Virtual Screening for Discovery of Indole Derivatives as Potent DNA Gyrase ATPase Inhibitors Active against Mycobacterium tuberculosis and Hit Validation by Biological Assays.

Mycobacterium tuberculosis is the single most important global infectious disease killer and a World Health Organization critical priority pathogen for development of new antimicrobials. M. tuberculosis DNA gyrase is a validated target for anti-TB agents, but those in current use target DNA breakage-reunion, rather than the ATPase activity of the GyrB subunit. Here, virtual screening, subsequently validated by whole-cell and enzyme inhibition assays, was applied to identify candidate compounds that inhibit M. tuberculosis GyrB ATPase activity from the Specs compound library. This approach yielded six compounds: four carbazole derivatives (1, 2, 3, and 8), the benzoindole derivative 11, and the indole derivative 14. Carbazole derivatives can be considered a new scaffold for M. tuberculosis DNA gyrase ATPase inhibitors. IC50 values of compounds 8, 11, and 14 (0.26, 0.56, and 0.08 µM, respectively) for inhibition of M. tuberculosis DNA gyrase ATPase activity are 5-fold, 2-fold, and 16-fold better than the known DNA gyrase ATPase inhibitor novobiocin. MIC values of these compounds against growth of M. tuberculosis H37Ra are 25.0, 3.1, and 6.2 µg/mL, respectively, superior to novobiocin (MIC > 100.0 µg/mL). Molecular dynamics simulations of models of docked GyrB:inhibitor complexes suggest that hydrogen bond interactions with GyrB Asp79 are crucial for high-affinity binding of compounds 8, 11, and 14 to M. tuberculosis GyrB for inhibition of ATPase activity. These data demonstrate that virtual screening can identify known and new scaffolds that inhibit both M. tuberculosis DNA gyrase ATPase activity in vitro and growth of M. tuberculosis bacteria.

Elucidating specific interactions for designing novel pyrrolamide derivatives as potential GyrB inhibitors based on ab initio fragment molecular orbital calculations.

Tuberculosis (TB), the second leading infectious killer, causes serious public health problems worldwide. To develop novel anti-TB agents, many biochemical studies have targeted the subunit B of DNA gyrase (GyrB), which captures a second DNA segment and responses for ATP hydrolysis. Here, we investigated specific interactions between GyrB residues and existing pyrrolamide derivatives at an electronic level using ab initio fragment molecular orbital (FMO) calculations and designed potent inhibitors against GyrB. The evaluated binding affinities between GyrB and pyrrolamides were confirmed to be consistent with the IC50 values obtained from previous experiments. Thus, we employed the most potent pyrrolamide (compound 1) as a lead compound and proposed novel pyrrolamide derivatives. The specific interactions between GyrB and these derivatives were investigated using molecular mechanic optimizations and FMO calculations. The results revealed that our proposed derivatives had strong hydrogen bonds with Asp79 and Arg141 and exhibited electrostatic interactions with Glu56 and Ile84 of GyrB. In addition, the binding affinity between GyrB and compound 1 was enhanced significantly by the replacement at the R3 site of compound 1. The present results may provide structural concepts for the rational design of potent GyrB inhibitors as anti-TB agents.Communicated by Ramaswamy H. Sarma.

Bioisosteric Design Identifies Inhibitors of Mycobacterium tuberculosis DNA Gyrase ATPase Activity.

Mutations in DNA gyrase confer resistance to fluoroquinolones, second-line antibiotics for Mycobacterium tuberculosis infections. Identification of new agents that inhibit M. tuberculosis DNA gyrase ATPase activity is one strategy to overcome this. Here, bioisosteric designs using known inhibitors as templates were employed to define novel inhibitors of M. tuberculosis DNA gyrase ATPase activity. This yielded the modified compound R3-13 with improved drug-likeness compared to the template inhibitor that acted as a promising ATPase inhibitor against M. tuberculosis DNA gyrase. Utilization of compound R3-13 as a virtual screening template, supported by subsequent biological assays, identified seven further M. tuberculosis DNA gyrase ATPase inhibitors with IC50 values in the range of 0.42-3.59 µM. The most active compound 1 showed an IC50 value of 0.42 µM, 3-fold better than the comparator ATPase inhibitor novobiocin (1.27 µM). Compound 1 showed noncytotoxicity to Caco-2 cells at concentrations up to 76-fold higher than its IC50 value. Molecular dynamics simulations followed by decomposition energy calculations identified that compound 1 occupies the binding pocket utilized by the adenosine group of the ATP analogue AMPPNP in the M. tuberculosis DNA gyrase GyrB subunit. The most prominent contribution to the binding of compound 1 to M. tuberculosis GyrB subunit is made by residue Asp79, which forms two hydrogen bonds with the OH group of this compound and also participates in the binding of AMPPNP. Compound 1 represents a potential new scaffold for further exploration and optimization as a M. tuberculosis DNA gyrase ATPase inhibitor and candidate anti-tuberculosis agent.

Structure-based drug design of novel M. tuberculosis InhA inhibitors based on fragment molecular orbital calculations.

2-trans enoyl-acyl carrier protein reductase (InhA) is a promising target for developing novel chemotherapy agents for tuberculosis, and their inhibitory effects on InhA activity were widely investigated by the physicochemical experiments. However, the reason for the wide range of their inhibitory effects induced by similar agents was not explained by only the difference in their chemical structures. In our previous molecular simulations, a series of heteroaryl benzamide derivatives were selected as candidate inhibitors against InhA, and their binding properties with InhA were investigated to propose novel derivatives with higher binding affinity to InhA. In the present study, we extended the simulations for a series of 4-hydroxy-2-pyridone derivatives to search widely for more potent inhibitors against InhA. Using ab initio fragment molecular orbital (FMO) calculations, we elucidated the specific interactions between InhA residues and the derivatives at an electronic level and highlighted key interactions between InhA and the derivatives. The FMO results clearly indicated that the most potent inhibitor has strong hydrogen bonds with the backbones of Tyr158, Thr196, and NADH of InhA. This finding may provide informative structural concepts for designing novel 4-hydroxy-2-pyridone derivatives with higher binding affinity to InhA. Our previous and present molecular simulations could provide important guidelines for the rational design of more potent InhA inhibitors.

Identification of Potent DNA Gyrase Inhibitors Active against Mycobacterium tuberculosis.

Mycobacterium tuberculosis DNA gyrase manipulates the DNA topology using controlled breakage and religation of DNA driven by ATP hydrolysis. DNA gyrase has been validated as the enzyme target of fluoroquinolones (FQs), second-line antibiotics used for the treatment of multidrug-resistant tuberculosis. Mutations around the DNA gyrase DNA-binding site result in the emergence of FQ resistance in M. tuberculosis; inhibition of DNA gyrase ATPase activity is one strategy to overcome this. Here, virtual screening, subsequently validated by biological assays, was applied to select candidate inhibitors of the M. tuberculosis DNA gyrase ATPase activity from the Specs compound library (www.specs.net). Thirty compounds were identified and selected as hits for in vitro biological assays, of which two compounds, G24 and G26, inhibited the growth of M. tuberculosis H37Rv with a minimal inhibitory concentration of 12.5 µg/mL. The two compounds inhibited DNA gyrase ATPase activity with IC50 values of 2.69 and 2.46 µM, respectively, suggesting this to be the likely basis of their antitubercular activity. Models of complexes of compounds G24 and G26 bound to the M. tuberculosis DNA gyrase ATP-binding site, generated by molecular dynamics simulations followed by pharmacophore mapping analysis, showed hydrophobic interactions of inhibitor hydrophobic headgroups and electrostatic and hydrogen bond interactions of the polar tails, which are likely to be important for their inhibition. Decreasing compound lipophilicity by increasing the polarity of these tails then presents a likely route to improving the solubility and activity. Thus, compounds G24 and G26 provide attractive starting templates for the optimization of antitubercular agents that act by targeting DNA gyrase.

Discovery of New and Potent InhA Inhibitors as Antituberculosis Agents: Structure-Based Virtual Screening Validated by Biological Assays and X-ray Crystallography.

The enoyl-acyl carrier protein reductase InhA of Mycobacterium tuberculosis is an attractive, validated target for antituberculosis drug development. Moreover, direct inhibitors of InhA remain effective against InhA variants with mutations associated with isoniazid resistance, offering the potential for activity against MDR isolates. Here, structure-based virtual screening supported by biological assays was applied to identify novel InhA inhibitors as potential antituberculosis agents. High-speed Glide SP docking was initially performed against two conformations of InhA differing in the orientation of the active site Tyr158. The resulting hits were filtered for drug-likeness based on Lipinski's rule and avoidance of PAINS-like properties and finally subjected to Glide XP docking to improve accuracy. Sixteen compounds were identified and selected for in vitro biological assays, of which two (compounds 1 and 7) showed MIC of 12.5 and 25 µg/mL against M. tuberculosis H37Rv, respectively. Inhibition assays against purified recombinant InhA determined IC50 values for these compounds of 0.38 and 0.22 µM, respectively. A crystal structure of the most potent compound, compound 7, bound to InhA revealed the inhibitor to occupy a hydrophobic pocket implicated in binding the aliphatic portions of InhA substrates but distant from the NADH cofactor, i.e., in a site distinct from those occupied by the great majority of known InhA inhibitors. This compound provides an attractive starting template for ligand optimization aimed at discovery of new and effective compounds against M. tuberculosis that act by targeting InhA.

Specific interactions between 2-trans enoyl-acyl carrier protein reductase and its ligand: Protein-ligand docking and ab initio fragment molecular orbital calculations.

2-trans enoyl-acyl carrier protein reductase (InhA) has been identified as a promising target for the development of novel chemotherapy for tuberculosis. In the present study, a series of heteroaryl benzamide derivatives were selected as potent inhibitors against InhA, and their binding properties with InhA were investigated at atomic and electronic levels by ab initio molecular simulations based on protein-ligand docking, classical molecular mechanics optimizations and ab initio fragment molecular orbital (FMO) calculations. The results evaluated by FMO highlight some key interactions between InhA and the derivatives, indicating that the most potent derivative has strong hydrogen bonds with the Met98 side chain of InhA and strong electrostatic interactions with the nicotinamide adenine dinucleotide cofactor. These findings provide informative structural concepts for designing novel heteroaryl benzamide derivatives with higher binding affinity to InhA.