Identification of potential inhibitors of Mtb InhA: a pharmacoinformatics approach

The emergence of superbugs of multi-drug resistant (MDR/RR) and extensively drug-resistant (XDR) Mycobacterium tuberculosis (Mtb) strains at a faster rate is posing a serious threat to Tuberculosis (TB) control worldwide. Mtb enoyl-acyl carrier protein reductase (InhA) is a well-established target of the front-line anti-TB prodrug Isoniazid (INH), which requires activation by Catalase-peroxidase enzyme (KatG) in order to inhibit InhA enzyme, that is crucial for the biosynthesis of the mycobacterial cell wall. Currently, due to widespread resistance to this drug, it is necessary to identify new clinical candidates that directly inhibit InhA enzyme and do not require activation by KatG, thereby circumventing most of the resistance mechanisms. In the present study, high-throughput virtual screening of ASINEX database was carried out to identify potential direct inhibitors of Mtb InhA. Best twenty compounds with good binding energies ranging between -12.36 and -9.27 kcal/mol were selected as promising virtual screening hits. These molecules were subjected to ADME study followed by toxicity prediction. Finally, four top-ranked molecules which are structurally diverse and possess best binding affinity than the co-crystalized ligand have been chosen for MD simulation studies followed by MM-GBSA analysis to validate and ensure the stability of hits in the active site of the enzyme. Based on the 100 ns MD simulation studies and binding free energy estimates, three hit molecules B244, B369, and B310 could be considered as potential inhibitors for Mtb InhA, which are likely to be potent against INH-resistant Mtb strains after successful experimental validation.Communicated by Ramaswamy H. Sarma.

Rethinking the MtInhA tertiary and quaternary structure flexibility: a molecular dynamics view

Flexibility and function are related properties in the study of protein dynamics. Flexibility reflects in the conformational potential of proteins and thus in their functionalities. The presence of interactions between protein-ligands and protein-protein complexes, substrates, and environmental changes can alter protein plasticity, acting from the rearrangement of the side chains of amino acids to the folding/unfolding of large structural motifs. To evaluate the effects of the flexibility in protein systems, we defined the enzyme 2-trans-enoyl-ACP (CoA) reductase from Mycobacterium tuberculosis, or MtInhA, as our target system. MtInhA is biologically active as a tetramer in solution; however, computational studies commonly use the monomer justifying the independence of its active sites due to their distances. However, differences in flexibility between tertiary and quaternary structures could present impact on the size of the active site, influencing the drug discovery process. In this study, we investigated the influence of flexibility restrictions in A- and B-loops of the MtInhA in order to suggest a monomeric structure that describes the conformational behavior of the tetrameric system. Overall, we observed that simulations where restrictions were applied to the A- and B-loops present a more similar behavior to the native structure when compared to unrestricted simulations. Therefore, our work presents a monomeric model of MtInhA, which has conformational characteristics of the biologically active structure. Thus, the data obtained in this work can be applied to the MtInhA system for the generation of more reliable flexible models for molecular docking experiments, and also for the performance of longer simulations by molecular dynamics and with a lower computational cost.

Inhibition of Mycobacterium tuberculosis InhA by 3-nitropropanoic acid

3-Nitropropanoic acid (3NP), a bioactive fungal natural product, was previously demonstrated to inhibit growth of Mycobacterium tuberculosis. Here we demonstrate that 3NP inhibits the 2-trans-enoyl-acyl carrier protein reductase (InhA) from Mycobacterium tuberculosis with an IC₅₀ value of 71 μM, and present the crystal structure of the ternary InhA-NAD⁺ -3NP complex. The complex contains the InhA substrate-binding loop in an ordered, open conformation with Tyr158, a catalytically important residue whose orientation defines different InhA substrate/inhibitor complex conformations, in the "out" position. 3NP occupies a hydrophobic binding site adjacent to the NAD⁺ cofactor and close to that utilized by the diphenyl ether triclosan, but binds predominantly via electrostatic and water-mediated hydrogen-bonding interactions with the protein backbone and NAD⁺ cofactor. The identified mode of 3NP binding provides opportunities to improve inhibitory activity toward InhA.

Mycobacterium enoyl acyl carrier protein reductase (InhA): A key target for antitubercular drug discovery

Enoyl acyl carrier protein reductase (InhA) is a key enzyme involved in fatty acid synthesis mainly mycolic acid biosynthesis that is a part of NADH dependent acyl carrier protein reductase family. The aim of the present literature is to underline the different scaffolds or enzyme inhibitors that inhibit mycolic acid biosynthesis mainly cell wall synthesis by inhibiting enzyme InhA. Various scaffolds were identified based on the screening technologies like high throughput screening, encoded library technology, fragment-based screening. The compounds studied include indirect inhibitors (Isoniazid, Ethionamide, Prothionamide) and direct inhibitors (Triclosan/Diphenyl ethers, Pyrrolidine Carboxamides, Pyrroles, Acetamides, Thiadiazoles, Triazoles) with better efficacy against drug resistance. Out of the several scaffolds studied, pyrrolidine carboxamides were found to be the best molecules targeting InhA having good bioavailability properties and better MIC. This review provides with a detailed information, analysis, structure activity relationship and useful insight on various scaffolds as InhA inhibitors.

Design, synthesis, in silico and in vitro evaluation of novel diphenyl ether derivatives as potential antitubercular agents

Diphenyl ether derivatives inhibit mycobacterial cell wall synthesis by inhibiting an enzyme, enoyl-acyl carrier protein reductase (InhA), which catalyses the last step in the fatty acid synthesis cycle of genus Mycobacterium. To select and validate a protein crystal structure of enoyl-acyl carrier protein reductase of Mycobacterium tuberculosis for designing inhibitors using molecular modelling, a cross-docking and correlation study was performed. A series of novel 1-(3-(3-hydroxy-4-phenoxyphenyl)-5-phenyl-4,5-dihydro-1H-pyrazol-1-yl) ethan-1-ones were synthesized from this model and screened for their antitubercular activity against M. tuberculosis H37Rv. Compound PYN-8 showed good antitubercular activity on M. tuberculosis H37Rv (MIC = 4-7 µM) and Mycobacterium bovis (% inhibition at 10 µM = 95.91%). Cytotoxicity of all the synthesized derivatives was assessed using various cell lines, and they were found to be safe. Structure of PYN-8 was also confirmed by single-crystal X-ray diffraction. The molecular modelling studies also corroborated the biological activity of the compounds. Further, in silico findings revealed that all these tested compounds exhibited good ADME properties and drug likeness and thus may be considered as potential candidates for further drug development.

Identifying the Structural Features of Diphenyl Ether Analogues for InhA Inhibition: A 2D and 3D QSAR Based Study

Background:

Enoyl acyl carrier protein reductase (InhA) is a validated target for Mycobacterium. It is an enzyme which is associated with the biosynthesis of mycolic acids in type II fatty acid synthase system. Mycobacterial cell wall majorly comprises mycolic acids, which are responsible for virulence of the microorganism. Several diphenyl ether derivatives have been known to be direct inhibitors of InhA.

Objective:

In the present work, a Quantitative Structure Activity Relationship (QSAR) study was performed to identify the structural features of diphenyl ether analogues which contribute to InhA inhibitory activity in a favourable way.

Method:

Both 2D and 3D QSAR models were built and compared. Several fingerprint based 2D QSAR models were generated and their relationship with the structural features was studied. Models which corroborated the inhibitory activity of the molecules with their structural features were selected and studied in detail.

Results:

A 2D-QSAR model, with dendritic fingerprints having regression coefficient, for test set molecules Q2 =0.8132 and for the training set molecules, R2 =0.9607 was obtained. Additionally, an atom-based 3D-QSAR model with Q2 =0.7697 and R2 =0.9159 was also constructed.

Conclusion:

The data reported by various models provides guidance for the designing of structurally new diphenyl ether inhibitors with potential activity against InhA of M. tuberculosis.

An integrated computational approach of molecular dynamics simulations, receptor binding studies and pharmacophore mapping analysis in search of potent inhibitors against tuberculosis

Tuberculosis is an infectious chronic disease caused by obligate pathogen Mycobacterium tuberculosis that affects millions of people worldwide. Although many first and second line drugs are available for its treatment, but their irrational use has adversely lead to the emerging cases of multiple drug resistant and extensively drug-resistant tuberculosis. Therefore, there is an intense need to develop novel potent analogues for its treatment. This has prompted us to develop potent analogues against TB. The Mycobacterium tuberculosis genome provides us with number of validated targets to combat against TB. Study of Mtb genome disclosed six epoxide hydrolases (A to F) which convert harmful epoxide into diols and act as a potential drug target for rational drug design. Our current strategy is to develop such analogues which inhibits epoxide hydrolase enzyme present in Mtb genome. To achieve this, we adopted an integrated computational approach involving QSAR, pharmacophore mapping, molecular docking and molecular dynamics simulation studies. The approach envisaged vital information about the role of molecular descriptors, essential pharmacophoric features and binding energy for compounds to bind into the active site of epoxide hydrolase. Molecular docking analysis revealed that analogues exhibited significant binding to Mtb epoxide hydrolase. Further, three docked complexes 2s, 37s and 15s with high, moderate and low docking scores respectively were selected for molecular dynamics simulation studies. RMSD analysis revealed that all complexes are stable with average RMSD below 2 Å throughout the 10 ns simulations. The B-factor analysis showed that the active site residues of epoxide hydrolase are flexible enough to interact with inhibitor. Moreover, to confirm the binding of these urea derivatives, MM-GBSA binding energy analysis were performed. The calculations showed that 37s has more binding affinity (ΔGtotal = -52.24 kcal/mol) towards epoxide hydrolase compared to 2s (ΔGtotal = -51.70 kcal/mol) and 15s (ΔGtotal = -49.97 kcal/mol). The structural features inferred in our study may provide the future directions to the scientists towards the discovery of new chemical entity exhibiting anti-TB property.

Cell wall permeability assisted virtual screening to identify potential direct InhA inhibitors of Mycobacterium tuberculosis and their biological evaluation

The arising cases of isoniazid-resistance have motivated research interests toward new class of molecules known as direct InhA inhibitors. Here, a combine approach of shape-based pharmacophore and descriptor-based 2D QSAR was used to identify the potential direct InhA inhibitors. The approach is duly assisted with in vitro testing and molecular dynamics simulations. A combination of empirical parameters was derived to use as a filter for cell wall permeability while 2D QSAR was used as another filter to predict the biological activity. Both filters were applied to prioritize the molecules for biological evaluation against anti-TB activity. It led to 6 potential molecules which showed > 90% inhibition of H37Rv strain of Mycobacterium tuberculosis in BACTEC assay. Further, MMGBSA binding free energy of identified molecules was compared with available highly potent molecule, 5-hexyl-2-(2-methylphenoxy) phenol (IC₅₀ = 5nM) using molecular dynamics simulations. It showed two molecules with comparatively higher affinity toward InhA as compared to potent molecule. It indicated the candidature of identified molecules to be further considered in anti-TB drug development pipeline.

Recent Advances and Structural Features of Enoyl-ACP Reductase Inhibitors of Mycobacterium tuberculosis

Mycobacterium tuberculosis enoyl-ACP reductase (InhA) has been validated as a promising target for antitubercular agents. Isoniazid (INH), the most prescribed drug to treat tuberculosis (TB), inhibits a NADH-dependent InhA that provides precursors of mycolic acids, which are components of the mycobacterial cell wall. It is a pro-drug that needs activation to form the inhibitory INH-NAD adduct by KatG coding for catalase-peroxidase. The INH resistance of M. tuberculosis is caused by mutations in KatG, which may lead to multidrug-resistant TB (MDR-TB). Hence, there is a need for new drugs that can combat MDR-TB. The rationale for the development of new drugs to combat MDR-TB strains is the design of InhA inhibitors that can bypass bioactivation by KatG. In the present review, special attention was paid to discuss the chemical nature and recent developments of direct InhA inhibitors. The InhA inhibitors reported here have significant inhibitory effects against Mtb InhA. The diphenyl ether derivatives have shown slow onset, a tight-binding mechanism, and high affinity at the InhA active site. However, some of the diphenyl ethers have significant in vitro efficacy, which fails to transform into in vivo efficacy. Among the InhA inhibitors, 4-hydroxy-2-pyridones have emerged as a new chemical class with significant InhA inhibitory activity and better pharmacokinetic parameters when compared to diphenyl ethers.

Aqueous Molecular Dynamics Simulations of the M. tuberculosis Enoyl-ACP Reductase-NADH System and Its Complex with a Substrate Mimic or Diphenyl Ethers Inhibitors

Molecular dynamics (MD) simulations of 12 aqueous systems of the NADH-dependent enoyl-ACP reductase from Mycobacterium tuberculosis (InhA) were carried out for up to 20–40 ns using the GROMACS 4.5 package. Simulations of the holoenzyme, holoenzyme-substrate, and 10 holoenzyme-inhibitor complexes were conducted in order to gain more insight about the secondary structure motifs of the InhA substrate-binding pocket. We monitored the lifetime of the main intermolecular interactions: hydrogen bonds and hydrophobic contacts. Our MD simulations demonstrate the importance of evaluating the conformational changes that occur close to the active site of the enzyme-cofactor complex before and after binding of the ligand and the influence of the water molecules. Moreover, the protein-inhibitor total steric (E_LJ) and electrostatic (E_C) interaction energies, related to Gly96 and Tyr158, are able to explain 80% of the biological response variance according to the best linear equation, pK_i = 7.772 − 0.1885 × Gly96 + 0.0517 × Tyr158 (R² = 0.80; n = 10), where interactions with Gly96, mainly electrostatic, increase the biological response, while those with Tyr158 decrease. These results will help to understand the structure-activity relationships and to design new and more potent anti-TB drugs.