Crystal structure and function of the isoniazid target of Mycobacterium tuberculosis

Identification of some novel pyrazolo[1,5-a]pyrimidine derivatives as InhA inhibitors through pharmacophore-based virtual screening and molecular docking

The InhA inhibitors play key role in mycolic acid synthesis by preventing the fatty acid biosynthesis pathway. In this present article, Pharmacophore modelling and molecular docking study followed by in silico virtual screening could be considered as effective strategy to identify newer enoyl-ACP reductase inhibitors. Pyrrolidine carboxamide derivatives were opted to generate pharmacophore models using HypoGen algorithm in Discovery studio 2.1. Further it was employed to screen Zinc and Minimaybridge databases to identify and design newer potent hit molecules. The retrieved newer hits were further evaluated for their drug likeliness and docked against enoyl acyl carrier protein reductase. Here, novel pyrazolo[1,5-a]pyrimidine analogues were designed and synthesized with good yields. Structural elucidation of synthesized final molecules was perform through IR, MASS, ¹H-NMR, ¹³C-NMR spectroscopy and further tested for its in vitro anti-tubercular activity against H37Rv strain using Microplate Alamar blue assay (MABA) method. Most of the synthesized compounds displayed strong anti-tubercular activities. Further, these potent compounds were gauged for MDR-TB, XDR-TB and cytotoxic study.

Screening of a Novel Fragment Library with Functional Complexity against Mycobacterium tuberculosis InhA

Our findings reported herein provide support for the benefits of including functional group complexity (FGC) within fragments when screening against protein targets such as Mycobacterium tuberculosis InhA. We show that InhA fragment actives with FGC maintained their binding pose during elaboration. Furthermore, weak fragment hits with functional group handles also allowed for facile fragment elaboration to afford novel and potent InhA inhibitors with good ligand efficiency metrics for optimization.

New active leads for tuberculosis booster drugs by structure-based drug discovery

The transcriptional repressor EthR from Mycobacterium tuberculosis, a member of the TetR family of prokaryotic homo-dimeric transcription factors, controls the expression of the mycobacterial mono-oxygenase EthA. EthA is responsible for the bio-activation of the second-line tuberculosis pro-drug ethionamide, and consequently EthR inhibitors boost drug efficacy. Here, we present a comprehensive in silico structure-based screening protocol that led to the identification of a number of novel scaffolds of EthR inhibitors in subsequent biophysical screening by thermal shift assay. Growth inhibition assays demonstrated that five of the twenty biophysical hits were capable of boosting ethionamide activity in vitro, with the best novel scaffold displaying an EC₅₀ of 34 μM. In addition, the co-crystal structures of EthR with four new ligands at resolution ranging from 2.1 to 1.4 Å confirm the binding and inactivation mode, and will enable future lead development.

An overview on crystal structures of InhA protein: Apo-form, in complex with its natural ligands and inhibitors

The enoyl-ACP reductase InhA from the mycobacterial fatty acid biosynthesis pathway has become a target of interest for the development of new anti-tubercular drugs. This protein has been identified as essential for the survival of Mycobacterium tuberculosis, the causative agent of tuberculosis, and as the main target of two pro-drugs: isoniazid, the frontline anti-tubercular drug, and ethionamide, a second-line medicine. Since most cases of resistance to isoniazid and ethionamide result from mutations in the mycobacterial activating enzyme (KatG for isoniazid and EthA for ethionamide), research of direct InhA inhibitors, avoiding the activation step, has emerged as a promising strategy for combating tuberculosis. Thereby, InhA is drawing much attention and its three-dimensional structure has been particularly studied. A better understanding of key sites of interactions responsible for InhA inhibition arises thus as an essential tool for the rational design of new potent inhibitors. In this paper, we propose an overview of the 80 available crystal structures of wild-type and mutant InhA, in its apo form, in complex with its cofactor, with an analogue of its natural ligands (C16 fatty acid and/or NADH) or with inhibitors. We will first discuss structural and mechanistic aspects in order to highlight key features of the protein before delivering thorough inventory of structures of InhA in the presence of synthetic ligands to underline the key interactions implicated in high affinity inhibition.

Molecular modeling studies and anti-TB activity of trisubstituted indolizine analogues; molecular docking and dynamic inputs

A series of trisubstituted indolizine analogues has been designed as a result of a fragment-based approach to target the inhibition of mycobacterial enoyl-acyl carrier protein reductase. Anti-tuberculosis (TB) screening of the characterized compounds by a resazurin microplate assay method revealed that ethyl group at second position of indolizine nucleus exhibited activity against susceptible and multidrug-resistant strains of Mycobacterium tuberculosis at concentration of 5.5 and 11.3 μg/mL, respectively. A molecular docking study was also conducted to evaluate the stability of the active compounds, and compound with ethyl substitution at second position of indolizine nucleus showed the highest free binding energy of ΔG -24.11 (kcal/mol), a low clash score of 3.04, and high lipo score of -13.33. Indolizine analog with ethyl substitution at second position demonstrated Molecular Mechanics/Generalized Born Surface Area (-23.85 kcal/mol). Two molecular dynamics studies were computed (100 ps and 50 ns) to calculate the relationship between the potential and kinetic energies of the active anti-TB compound with time and temperature. The discovery of this lead may have a positive impact on anti-TB drug discovery.

Disruption of key NADH-binding pocket residues of the Mycobacterium tuberculosis InhA affects DD-CoA binding ability

Tuberculosis (TB) is a global health problem that affects over 10 million people. There is an urgent need to develop novel antimicrobial therapies to combat TB. To achieve this, a thorough understanding of key validated drug targets is required. The enoyl reductase InhA, responsible for synthesis of essential mycolic acids in the mycobacterial cell wall, is the target for the frontline anti-TB drug isoniazid. To better understand the activity of this protein a series of mutants, targeted to the NADH co-factor binding pocket were created. Residues P193 and W222 comprise a series of hydrophobic residues surrounding the cofactor binding site and mutation of both residues negatively affect InhA function. Construction of an M155A mutant of InhA results in increased affinity for NADH and DD-CoA turnover but with a reduction in V_max for DD-CoA, impairing overall activity. This suggests that NADH-binding geometry of InhA likely permits long-range interactions between residues in the NADH-binding pocket to facilitate substrate turnover in the DD-CoA binding region of the protein. Understanding the precise details of substrate binding and turnover in InhA and how this may affect protein-protein interactions may facilitate the development of improved inhibitors enabling the development of novel anti-TB drugs.

EthA/R-Independent Killing of Mycobacterium tuberculosis by Ethionamide

Ethionamide (ETH) is part of the drug arsenal available to treat multi-drug resistant tuberculosis. The current paradigm of this pro-drug activation involves the mycobacterial enzyme EthA and the transcriptional repressor, EthR. However, several lines of evidence suggest the involvement of additional players. The ethA/R locus was deleted in Mycobacterium bovis BCG and three Mycobacterium tuberculosis (MTB) strains. While complete resistance to ETH was observed with BCG ethA/R KO, drug susceptibility and dose-dependent killing were retained in the ethA/R KO MTB mutants, suggesting the existence of an alternative pathway of ETH bio-activation in MTB. We further demonstrated that this alternative pathway is EthR-independent, whereby re-introduction of ethR in ethA/R KO MTB did not lead to increased resistance to ETH. Consistently, ethA KO MTB (with intact ethR expression) displayed similar ETH susceptibility profile as their ethA/R KO counterparts. To identify the alternative ETH bio-activator, spontaneous ETH-resistant mutants were obtained from ethA/R KO MTB and whole genome sequencing identified single nucleotide polymorphisms in mshA, involved in mycothiol biosynthesis and previously linked to ETH resistance. Deletion of mshA in ethA/R KO MTB led to complete ETH resistance, supporting that the role of MshA in ETH killing is EthA/R-independent. Furthermore mshA single KO MTB displayed levels of ETH resistance similar or greater than those obtained with ethA/R KO strains, supporting that mshA is as critical as ethA/R for ETH killing efficacy.

Multiplex detection of extensively drug resistant tuberculosis using binary deoxyribozyme sensors

Current diagnostic tools for Mycobacterium tuberculosis (Mtb) have many disadvantages including low sensitivity, slow turnaround times, or high cost. Accurate, easy to use, and inexpensive point of care molecular diagnostic tests are urgently needed for the analysis of multidrug resistant (MDR) and extensively drug resistant (XDR) Mtb strains that emerge globally as a public health threat. In this study, we established proof-of-concept for a novel diagnostic platform (TB-DzT) for Mtb detection and the identification of drug resistant mutants using binary deoxyribozyme sensors (BiDz). TB-DzT combines a multiplex PCR with single nucleotide polymorphism (SNP) detection using highly selective BiDz sensors targeting loci associated with species typing and resistance to rifampin, isoniazid and fluoroquinolone antibiotics. Using the TB-DzT assay, we demonstrated accurate detection of Mtb and 5 mutations associated with resistance to three anti-TB drugs in clinical isolates. The assay also enables detection of a minority population of drug resistant Mtb, a clinically relevant scenario referred to as heteroresistance. Additionally, we show that TB-DzT can detect the presence of unknown mutations at target loci using combinatorial BiDz sensors. This diagnostic platform provides the foundation for the development of cost-effective, accurate and sensitive alternatives for molecular diagnostics of MDR- and XDR-TB.

Evaluating the Contribution of Transition-State Destabilization to Changes in the Residence Time of Triazole-Based InhA Inhibitors

A critical goal of lead compound selection and optimization is to maximize target engagement while minimizing off-target binding. Since target engagement is a function of both the thermodynamics and kinetics of drug-target interactions, it follows that the structures of both the ground states and transition states on the binding reaction coordinate are needed to rationally modulate the lifetime of the drug-target complex. Previously, we predicted the structure of the rate-limiting transition state that controlled the time-dependent inhibition of the enoyl-ACP reductase InhA. This led to the discovery of a triazole-containing diphenyl ether with an increased residence time on InhA due to transition-state destabilization rather than ground-state stabilization. In the present work, we evaluate the inhibition of InhA by 14 triazole-based diphenyl ethers and use a combination of enzyme kinetics and X-ray crystallography to generate a structure-kinetic relationship for time-dependent binding. We show that the triazole motif slows the rate of formation for the final drug-target complex by up to 3 orders of magnitude. In addition, we identify a novel inhibitor with a residence time on InhA of 220 min, which is 3.5-fold longer than that of the INH-NAD adduct formed by the tuberculosis drug, isoniazid. This study provides a clear example in which the lifetime of the drug-target complex is controlled by interactions in the transition state for inhibitor binding rather than the ground state of the enzyme-inhibitor complex, and demonstrates the important role that on-rates can play in drug-target residence time.

Discovery of cofactor-specific, bactericidal Mycobacterium tuberculosis InhA inhibitors using DNA-encoded library technology

Millions of individuals are infected with and die from tuberculosis (TB) each year, and multidrug-resistant (MDR) strains of TB are increasingly prevalent. As such, there is an urgent need to identify novel drugs to treat TB infections. Current frontline therapies include the drug isoniazid, which inhibits the essential NADH-dependent enoyl-acyl-carrier protein (ACP) reductase, InhA. To inhibit InhA, isoniazid must be activated by the catalase-peroxidase KatG. Isoniazid resistance is linked primarily to mutations in the katG gene. Discovery of InhA inhibitors that do not require KatG activation is crucial to combat MDR TB. Multiple discovery efforts have been made against InhA in recent years. Until recently, despite achieving high potency against the enzyme, these efforts have been thwarted by lack of cellular activity. We describe here the use of DNA-encoded X-Chem (DEX) screening, combined with selection of appropriate physical properties, to identify multiple classes of InhA inhibitors with cell-based activity. The utilization of DEX screening allowed the interrogation of very large compound libraries (10¹¹ unique small molecules) against multiple forms of the InhA enzyme in a multiplexed format. Comparison of the enriched library members across various screening conditions allowed the identification of cofactor-specific inhibitors of InhA that do not require activation by KatG, many of which had bactericidal activity in cell-based assays.

Drug resistance mechanisms and novel drug targets for tuberculosis therapy

Drug-resistant tuberculosis (TB) poses a significant challenge to the successful treatment and control of TB worldwide. Resistance to anti-TB drugs has existed since the beginning of the chemotherapy era. New insights into the resistant mechanisms of anti-TB drugs have been provided. Better understanding of drug resistance mechanisms helps in the development of new tools for the rapid diagnosis of drug-resistant TB. There is also a pressing need in the development of new drugs with novel targets to improve the current treatment of TB and to prevent the emergence of drug resistance in Mycobacterium tuberculosis. This review summarizes the anti-TB drug resistance mechanisms, furnishes some possible novel drug targets in the development of new agents for TB therapy and discusses the usefulness using known targets to develop new anti-TB drugs. Whole genome sequencing is currently an advanced technology to uncover drug resistance mechanisms in M. tuberculosis. However, further research is required to unravel the significance of some newly discovered gene mutations in their contribution to drug resistance.

A new 'golden age' for the antitubercular target InhA

The increasing prevalence of multidrug-resistant strains of Mycobacterium tuberculosis is the main contributing factor in unfavorable outcomes in the treatment of tuberculosis. Studies suggest that direct inhibitors of InhA, an enoyl-ACP-reductase, might yield promising clinical candidates that can be developed into new antitubercular drugs. In this review, we describe the application of different hit-identification strategies to InhA, which clearly illustrate the druggability of its active site through distinct binding mechanisms. We further characterize four classes of InhA inhibitors that show novel binding modes, and provide evidence of their successful target engagement as well as their in vivo activity.

Thiolactomycin-Based Inhibitors of Bacterial β-Ketoacyl-ACP Synthases with in Vivo Activity

β-Ketoacyl-ACP synthases (KAS) are key enzymes involved in the type II bacterial fatty acid biosynthesis (FASII) pathway and are putative targets for antibacterial discovery. Several natural product KAS inhibitors have previously been reported, including thiolactomycin (TLM), which is produced by Nocardia spp. Here we describe the synthesis and characterization of optically pure 5R-thiolactomycin (TLM) analogues that show improved whole cell activity against bacterial strains including methicillin-resistant Staphylococcus aureus (MRSA) and priority pathogens such as Francisella tularensis and Burkholderia pseudomallei. In addition, we identify TLM analogues with in vivo efficacy against MRSA and Klebsiella pneumoniae in animal models of infection.

N-Benzyl-4-((heteroaryl)methyl)benzamides: A New Class of Direct NADH-Dependent 2-trans Enoyl-Acyl Carrier Protein Reductase (InhA) Inhibitors with Antitubercular Activity

Isoniazid (INH) remains one of the cornerstones of antitubercular chemotherapy for drug-sensitive strains of M. tuberculosis bacteria. However, the increasing prevalence of multidrug-resistant (MDR) and extensively drug-resistant (XDR) strains containing mutations in the KatG enzyme, which is responsible for the activation of INH into its antitubercular form, have rendered this drug of little or no use in many cases of drug-resistant tuberculosis. Presented herein is a novel family of antitubercular direct NADH-dependent 2-trans enoyl-acyl carrier protein reductase (InhA) inhibitors based on an N-benzyl-4-((heteroaryl)methyl)benzamide template; unlike INH, these do not require prior activation by KatG. Given their direct InhA target engagement, these compounds should be able to circumvent KatG-related resistance in the clinic. The lead molecules were shown to be potent inhibitors of InhA and showed activity against M. tuberculosis bacteria. This new family of inhibitors was found to be chemically tractable, as exemplified by the facile synthesis of analogues and the establishment of structure-activity relationships. Furthermore, a co-crystal structure of the initial hit with the enzyme is disclosed, providing valuable information toward the design of new InhA inhibitors for the treatment of MDR/XDR tuberculosis.

Resistance to Isoniazid and Ethionamide in Mycobacterium tuberculosis: Genes, Mutations, and Causalities

Isoniazid (INH) is the cornerstone of tuberculosis (TB) chemotherapy, used for both treatment and prophylaxis of TB. The antimycobacterial activity of INH was discovered in 1952, and almost as soon as its activity was published, the first INH-resistant Mycobacterium tuberculosis strains were reported. INH and its structural analog and second-line anti-TB drug ethionamide (ETH) are pro-drugs. INH is activated by the catalase-peroxidase KatG, while ETH is activated by the monooxygenase EthA. The resulting active species reacts with NAD+ to form an INH-NAD or ETH-NAD adduct, which inhibits the enoyl ACP reductase InhA, leading to mycolic acid biosynthesis inhibition and mycobacterial cell death. The major mechanism of INH resistance is mutation in katG, encoding the activator of INH. One specific KatG variant, S315T, is found in 94% of INH-resistant clinical isolates. The second mechanism of INH resistance is a mutation in the promoter region of inhA (c-15t), which results in inhA overexpression and leads to titration of the drug. Mutations in the inhA open reading frame and promoter region are also the major mechanism of resistance to ETH, found more often in ETH-resistant clinical isolates than mutations in the activator of ETH. Other mechanisms of resistance to INH and ETH include expression changes of the drugs' activators, redox alteration, drug inactivation, and efflux pump activation. In this article, we describe each known mechanism of resistance to INH and ETH and its importance in M. tuberculosis clinical isolates.

Molecular Basis of Drug Resistance in Mycobacterium tuberculosis

In this chapter we review the molecular mechanisms of drug resistance to the major first- and second-line antibiotics used to treat tuberculosis.

Molecular Dynamics Assisted Mechanistic Study of Isoniazid-Resistance against Mycobacterium tuberculosis InhA

Examination of InhA mutants I16T, I21V, I47T, S94A, and I95P showed that direct and water mediated H-bond interactions between NADH and binding site residues reduced drastically. It allowed conformational flexibility to NADH, particularly at the pyrophosphate region, leading to weakening of its binding at dinucleotide binding site. The highly scattered distribution of pyrophosphate dihedral angles and chi1 side chain dihedral angles of corresponding active site residues therein confirmed weak bonding between InhA and NADH. The average direct and water mediated bridged H-bond interactions between NADH and mutants were observed weaker as compared to wild type. Further, estimated NADH binding free energy in mutants supported the observed weakening of InhA-NADH interactions. Similarly, per residue contribution to NADH binding was also found little less as compared to corresponding residues in wild type. This investigation clearly depicted and supported the effect of mutations on NADH binding and can be accounted for isoniazid resistance as suggested by previous biochemical and mutagenic studies. Further, structural analysis of InhA provided the crucial points to enhance the NADH binding affinity towards InhA mutants in the presence of direct InhA inhibitors to combat isoniazid drug resistance. This combination could be a potential alternative for treatment of drug resistant tuberculosis.

Computer-Aided Design of Orally Bioavailable Pyrrolidine Carboxamide Inhibitors of Enoyl-Acyl Carrier Protein Reductase of Mycobacterium tuberculosis with Favorable Pharmacokinetic Profiles

We have carried out a computational structure-based design of new potent pyrrolidine carboxamide (PCAMs) inhibitors of enoyl-acyl carrier protein reductase (InhA) of Mycobacterium tuberculosis (MTb). Three-dimensional (3D) models of InhA-PCAMx complexes were prepared by in situ modification of the crystal structure of InhA-PCAM1 (Protein Data Bank (PDB) entry code: 4U0J), the reference compound of a training set of 20 PCAMs with known experimental inhibitory potencies (IC50(exp)). First, we built a gas phase quantitative structure-activity relationships (QSAR) model, linearly correlating the computed enthalpy of the InhA-PCAM complex formation and the IC50(exp). Further, taking into account the solvent effect and loss of inhibitor entropy upon enzyme binding led to a QSAR model with a superior linear correlation between computed Gibbs free energies (ΔΔGcom) of InhA-PCAM complex formation and IC50(exp) (pIC50(exp) = -0.1552·ΔΔGcom + 5.0448, R² = 0.94), which was further validated with a 3D-QSAR pharmacophore model generation (PH4). Structural information from the models guided us in designing of a virtual combinatorial library (VL) of more than 17 million PCAMs. The VL was adsorption, distribution, metabolism and excretion (ADME) focused and reduced down to 1.6 million drug like orally bioavailable analogues and PH4 in silico screened to identify new potent PCAMs with predicted IC50(pre) reaching up to 5 nM. Combining molecular modeling and PH4 in silico screening of the VL resulted in the proposed novel potent antituberculotic agent candidates with favorable pharmacokinetic profiles.

Dissecting the Structural Elements for the Activation of β-Ketoacyl-(Acyl Carrier Protein) Reductase from Vibrio cholerae

β-Ketoacyl-(acyl carrier protein) reductase (FabG) catalyzes the key reductive reaction in the elongation cycle of fatty acid synthesis (FAS), which is a vital metabolic pathway in bacteria and a promising target for new antibiotic development. The activation of the enzyme is usually linked to the formation of a catalytic triad and cofactor binding, and crystal structures of FabG from different organisms have been captured in either the active or inactive conformation. However, the structural elements which enable activation of FabG require further exploration. Here we report the findings of structural, enzymatic, and binding studies of the FabG protein found in the causative agent of cholera, Vibrio cholerae (vcFabG). vcFabG exists predominantly as a dimer in solution and is able to self-associate to form tetramers, which is the state seen in the crystal structure. The formation of the tetramer may be promoted by the presence of the cofactor NADP(H). The transition between the dimeric and tetrameric states of vcFabG is related to changes in the conformations of the α4/α5 helices on the dimer-dimer interface. Two glycine residues adjacent to the dimer interface (G92 and G141) are identified to be the hinge for the conformational changes, while the catalytic tyrosine (Y155) and a glutamine residue that forms hydrogen bonds to both loop β4-α4 and loop β5-α5 (Q152) stabilize the active conformation. The functions of the aforementioned residues were confirmed by binding and enzymatic assays for the corresponding mutants.

IMPORTANCE

This paper describes the results of structural, enzymatic, and binding studies of FabG from Vibrio cholerae (vcFabG). In this work, we dissected the structural elements responsible for the activation of vcFabG. The structural information provided here is essential for the development of antibiotics specifically targeting bacterial FabG, especially for the multidrug-resistant strains of V. cholerae.

A user-friendly web portal for analyzing conformational changes in structures of Mycobacterium tuberculosis

Initiation of the Tuberculosis Structural Consortium has resulted in the expansion of the Mycobacterium tuberculosis (MTB) protein structural database. Currently, 969 experimentally solved structures are available for 354 MTB proteins. This includes multiple crystal structures for a given protein under different functional conditions, such as the presence of different ligands or mutations. In depth analysis of the multiple structures reveal that subtle differences exist in conformations of a given protein under varied conditions. Therefore, it is immensely important to understand the conformational differences between the multiple structures of a given protein in order to select the most suitable structure for molecular docking and structure-based drug designing. Here, we introduce a web portal ( http://bmi.icmr.org.in/mtbsd/torsion.php ) that we developed to provide comparative data on the ensemble of available structures of MTB proteins, such as Cα root means square deviation (RMSD), sequence identity, presence of mutations and torsion angles. Additionally, torsion angles were used to perform principal component analysis (PCA) to identify the conformational differences between the structures. Additionally, we present a few case studies to demonstrate this database. Graphical Abstract Conformational changes seen in the structures of the enoyl-ACP reductase protein encoded by the Mycobacterial gene inhA.

Anti-tubercular drug development: computational strategies to identify potential compounds

InhA is an attractive target to combat tuberculosis (TB), which is targeted by many pro-drugs (isoniazid, etc.) and drugs such as triclosan. However, triclosan is less useful as an antitubercular drug due to its low bioavailability and therefore, in order to overcome this difficulty, many derivatives of triclosan were prepared. Here, we have combined various computational techniques to virtually screen out four potential triclosan derivatives. Molecular docking methods have been employed to screen out 32 out of 62 triclosan derivatives considering the mode of binding and the top re-rank scores. A comparative study on the chemical properties of triclosan and some of its derivatives has been performed using density functional theory (DFT) calculations. DFT based global reactivity descriptors (GRD), such as hardness, chemical potential, chemical softness, electrophilicity index, Fukui function, and local philicity calculated at the optimized geometries were used to investigate the usefulness of these descriptors for understanding the reactive nature and sites of the molecules. QSAR equations were built using these descriptors considering these 32 compounds. Four common compounds showing the best correlation and the best docking scores were considered for the ADMET property calculations and their dynamical movements have been studied using molecular dynamics simulations. Our results showed that these four compounds are chemically more active than triclosan and have the potential to inhibit the Mycobacterium tuberculosis enoyl acyl carrier protein reductase. This work shows that combination of different computational techniques may help to screen out potential drug candidates from a list of possible ones.

An Effective Approach for Clustering InhA Molecular Dynamics Trajectory Using Substrate-Binding Cavity Features

Protein receptor conformations, obtained from molecular dynamics (MD) simulations, have become a promising treatment of its explicit flexibility in molecular docking experiments applied to drug discovery and development. However, incorporating the entire ensemble of MD conformations in docking experiments to screen large candidate compound libraries is currently an unfeasible task. Clustering algorithms have been widely used as a means to reduce such ensembles to a manageable size. Most studies investigate different algorithms using pairwise Root-Mean Square Deviation (RMSD) values for all, or part of the MD conformations. Nevertheless, the RMSD only may not be the most appropriate gauge to cluster conformations when the target receptor has a plastic active site, since they are influenced by changes that occur on other parts of the structure. Hence, we have applied two partitioning methods (k-means and k-medoids) and four agglomerative hierarchical methods (Complete linkage, Ward's, Unweighted Pair Group Method and Weighted Pair Group Method) to analyze and compare the quality of partitions between a data set composed of properties from an enzyme receptor substrate-binding cavity and two data sets created using different RMSD approaches. Ensembles of representative MD conformations were generated by selecting a medoid of each group from all partitions analyzed. We investigated the performance of our new method for evaluating binding conformation of drug candidates to the InhA enzyme, which were performed by cross-docking experiments between a 20 ns MD trajectory and 20 different ligands. Statistical analyses showed that the novel ensemble, which is represented by only 0.48% of the MD conformations, was able to reproduce 75% of all dynamic behaviors within the binding cavity for the docking experiments performed. Moreover, this new approach not only outperforms the other two RMSD-clustering solutions, but it also shows to be a promising strategy to distill biologically relevant information from MD trajectories, especially for docking purposes.

Crystal structure of the enoyl-ACP reductase of Mycobacterium tuberculosis (InhA) in the apo-form and in complex with the active metabolite of isoniazid pre-formed by a biomimetic approach

InhA is an enoyl-ACP reductase of Mycobacterium tuberculosis implicated in the biosynthesis of mycolic acids, essential constituents of the mycobacterial cell wall. To date, this enzyme is considered as a promising target for the discovery of novel antitubercular drugs. In this work, we describe the first crystal structure of the apo form of the wild-type InhA at 1.80Å resolution as well as the crystal structure of InhA in complex with the synthetic metabolite of the antitubercular drug isoniazid refined to 1.40Å. This metabolite, synthesized in the absence of InhA, is able to displace and replace the cofactor NADH in the enzyme active site. This work provides a unique opportunity to enlighten the structural adaptation of apo-InhA to the binding of the NADH cofactor or of the isoniazid adduct. In addition, a differential scanning fluorimetry study of InhA, in the apo-form as well as in the presence of NAD(+), NADH and INH-NADH was performed showing that binding of the INH-NADH adduct had a strong stabilizing effect.

Tuberculosis Drug Development: History and Evolution of the Mechanism-Based Paradigm

Modern tuberculosis (TB) chemotherapy is widely viewed as a crowning triumph of anti-infectives research. However, only one new TB drug has entered clinical practice in the past 40 years while drug resistance threatens to further destabilize the pandemic. Here, we review a brief history of TB drug development, focusing on the evolution of mechanism(s)-of-action studies and key conceptual barriers to rational, mechanism-based drugs.

Elucidating structural basis of benzofuran pyrrolidine pyrazole derivatives for enhancing potency against both the InhA enzyme and intact M. tuberculosis cells: a combined MD simulations and 3D-QSAR study

The structural concept for enhancing both IC₅₀ and MIC₉₀ activities summarized from MD simulations and CoMSIA results.

Biological evaluation of potent triclosan-derived inhibitors of the enoyl-acyl carrier protein reductase InhA in drug-sensitive and drug-resistant strains of Mycobacterium tuberculosis

New triclosan (TRC) analogues were evaluated for their activity against the enoyl-acyl carrier protein reductase InhA in Mycobacterium tuberculosis (Mtb). TRC is a well-known inhibitor of InhA, and specific modifications to its positions 5 and 4' afforded 27 derivatives; of these compounds, seven derivatives showed improved potency over that of TRC. These analogues were active against both drug-susceptible and drug-resistant Mtb strains. The most active compound in this series, 4-(n-butyl)-1,2,3-triazolyl TRC derivative 3, had an MIC value of 0.6 μg mL(-1) (1.5 μM) against wild-type Mtb. At a concentration equal to its MIC, this compound inhibited purified InhA by 98 %, and showed an IC50 value of 90 nM. Compound 3 and the 5-methylisoxazole-modified TRC 14 were able to inhibit the biosynthesis of mycolic acids. Furthermore, mc(2) 4914, an Mtb strain overexpressing inhA, was found to be less susceptible to compounds 3 and 14, supporting the notion that InhA is the likely molecular target of the TRC derivatives presented herein.

Key Structures and Interactions for Binding of Mycobacterium tuberculosis Protein Kinase B Inhibitors from Molecular Dynamics Simulation

Substituted aminopyrimidine inhibitors have recently been introduced as antituberculosis agents. These inhibitors show impressive activity against protein kinase B, a Ser/Thr protein kinase that is essential for cell growth of M. tuberculosis. However, up to now, X-ray structures of the protein kinase B enzyme complexes with the substituted aminopyrimidine inhibitors are currently unavailable. Consequently, structural details of their binding modes are questionable, prohibiting the structural-based design of more potent protein kinase B inhibitors in the future. Here, molecular dynamics simulations, in conjunction with molecular mechanics/Poisson-Boltzmann surface area binding free-energy analysis, were employed to gain insight into the complex structures of the protein kinase B inhibitors and their binding energetics. The complex structures obtained by the molecular dynamics simulations show binding free energies in good agreement with experiment. The detailed analysis of molecular dynamics results shows that Glu93, Val95, and Leu17 are key residues responsible to the binding of the protein kinase B inhibitors. The aminopyrazole group and the pyrimidine core are the crucial moieties of substituted aminopyrimidine inhibitors for interaction with the key residues. Our results provide a structural concept that can be used as a guide for the future design of protein kinase B inhibitors with highly increased antagonistic activity.

Elucidating the structural basis of diphenyl ether derivatives as highly potent enoyl-ACP reductase inhibitors through molecular dynamics simulations and 3D-QSAR study

Diphenyl ether derivatives are good candidates for anti-tuberculosis agents that display a promising potency for inhibition of InhA, an essential enoyl-acyl carrier protein (ACP) reductase involved in fatty acid biosynthesis pathways in Mycobacterium tuberculosis. In this work, key structural features for the inhibition were identified by 3D-QSAR CoMSIA models, constructed based on available experimental binding properties of diphenyl ether inhibitors, and a set of four representative compounds was subjected to MD simulations of inhibitor-InhA complexes for the calculation of binding free energies. The results show that bulky groups are required for the R1 substituent on the phenyl A ring of the inhibitors to favor a hydrophobic pocket formed by residues Phe149, Met155, Pro156, Ala157, Tyr158, Pro193, Met199, Val203, Leu207, Ile215, and Leu218. Small substituents with a hydrophilic property are required at the R3 and R4 positions of the inhibitor phenyl B rings to form hydrogen bonds with the backbones of Gly96 and Met98, respectively. For the R2 substituent, small substituents with simultaneous hydrophilic or hydrophobic properties are required to favor the interaction with the pyrophosphate moiety of NAD(+) and the methyl side chain of Ala198, respectively. The reported data provide structural guidance for the design of new and potent diphenyl ether-based inhibitors with high inhibitory activities against M. tuberculosis InhA.

A structural characterization of the isoniazid Mycobacterium tuberculosis drug target, Rv2971, in its unliganded form

Aldo-keto reductases (AKR) are a large superfamily of NADPH-dependent oxidoreductases and play a role in detoxification of toxic metabolites. Rv2971, an AKR in Mycobacterium tuberculosis, has recently been identified as a target of isoniazid, a key first-line drug against tuberculosis. Here, the cloning, expression, purification, crystallization and structural characterization of Rv2971 are described. To gain insight into its function, the crystal structure of Rv2971 was successfully determined to 1.60 Å resolution in its unliganded form. The structure exhibits a TIM-barrel fold typical of AKRs, revealing structural characteristics essential for function and substrate specificities, allowing a structural comparison between Rv2971 and other mycobacterial AKRs.

Rational design of InhA inhibitors in the class of diphenyl ether derivatives as potential anti-tubercular agents using molecular dynamics simulations

A series of diphenyl ether derivatives were developed and showed promising potency for inhibiting InhA, an essential enoyl acyl carrier protein reductase involved in mycolic acid biosynthesis, leading to the lysis of Mycobacterium tuberculosis. To understand the structural basis of diphenyl ether derivatives for designing more potent inhibitors, molecular dynamics (MD) simulations were performed. Based on the obtained results, the dynamic behaviour in terms of flexibility, binding free energy, binding energy decomposition, conformation, and the inhibitor-enzyme interaction of diphenyl ether inhibitors were elucidated. Phe149, Tyr158, Met161, Met199, Val203 and NAD+ are the key residues for binding of diphenyl ether inhibitors in the InhA binding pocket. Our results could provide the structural concept to design new diphenyl ether inhibitors with better enzyme inhibitory activity against M. tuberculosis InhA. The present work facilitates the design of new and potentially more effective anti-tuberculosis agents.

Rational design of broad spectrum antibacterial activity based on a clinically relevant enoyl-acyl carrier protein (ACP) reductase inhibitor

Determining the molecular basis for target selectivity is of particular importance in drug discovery. The ideal antibiotic should be active against a broad spectrum of pathogenic organisms with a minimal effect on human targets. CG400549, a Staphylococcus-specific 2-pyridone compound that inhibits the enoyl-acyl carrier protein reductase (FabI), has recently been shown to possess human efficacy for the treatment of methicillin-resistant Staphylococcus aureus infections, which constitute a serious threat to human health. In this study, we solved the structures of three different FabI homologues in complex with several pyridone inhibitors, including CG400549. Based on these structures, we rationalize the 65-fold reduced affinity of CG400549 toward Escherichia coli versus S. aureus FabI and implement concepts to improve the spectrum of antibacterial activity. The identification of different conformational states along the reaction coordinate of the enzymatic hydride transfer provides an elegant visual depiction of the relationship between catalysis and inhibition, which facilitates rational inhibitor design. Ultimately, we developed the novel 4-pyridone-based FabI inhibitor PT166 that retained favorable pharmacokinetics and efficacy in a mouse model of S. aureus infection with extended activity against Gram-negative and mycobacterial organisms.

Time-dependent diaryl ether inhibitors of InhA: structure-activity relationship studies of enzyme inhibition, antibacterial activity, and in vivo efficacy

The diaryl ethers are a novel class of antituberculosis drug candidates that inhibit InhA, the enoyl-ACP reductase involved in the fatty acid biosynthesis (FASII) pathway, and have antibacterial activity against both drug-sensitive and drug-resistant strains of Mycobacterium tuberculosis. In the present work, we demonstrate that two time-dependent B-ring modified diaryl ether InhA inhibitors have antibacterial activity in a mouse model of TB infection when delivered by intraperitoneal injection. We propose that the efficacy of these compounds is related to their residence time on the enzyme, and to identify structural features that modulate drug-target residence time in this system, we have explored the inhibition of InhA by a series of B-ring modified analogues. Seven ortho-substituted compounds were found to be time-dependent inhibitors of InhA, where the slow step leading to the final enzyme-inhibitor complex (EI*) is thought to correlate with closure and ordering of the InhA substrate binding loop. A detailed mechanistic understanding of the molecular basis for residence time in this system will facilitate the development of InhA inhibitors with improved in vivo activity.

Crystal structures and kinetic properties of enoyl-acyl carrier protein reductase I from Candidatus Liberibacter asiaticus

Huanglongbing (HLB) is a destructive citrus disease. The leading cause of HLB is Candidatus Liberibacter asiaticus. Fatty acid biosynthesis is essential for bacterial viability and has been validated as a target for the discovery of novel antibacterial agents. Enoyl-acyl carrier protein reductase (also called ENR or FabI and a product of the fabI gene) is an enzyme required in a critical step of bacterial fatty acid biosynthesis and has attracted attention as a target of novel antimicrobial agents. We determined the crystal structures of FabI from Ca. L. asiaticus in its apoform as well as in complex with b-nicotinamide adenine dinucleotide (NAD) at 1.7 and 2.7 Å resolution, respectively, to facilitate the design and screening of small molecule inhibitors of FabI. The monomeric ClFabI is highly similar to other known FabI structures as expected; however, unlike the typical tetramer, ClFabI exists as a hexamer in crystal, whereas as dimer in solution, on the other hand, the substrate binding loop which always disordered in apoform FabI structures is ordered in apo-ClFabI. Interestingly, the structure of ClFabI undergoes remarkable conformational change in the substrate-binding loop in the presence of NAD. We conclude that the signature sequence motif of FabI can be considered as Gly-(Xaa)5-Ser-(Xaa)n-Val-Tyr-(Xaa)6-Lys-(Xaa)n-Thr instead of Tyr-(Xaa)6-Lys. We have further identified isoniazid as a competitive inhibitor with NADH.

N-substituted 2-isonicotinoylhydrazinecarboxamides--new antimycobacterial active molecules

This report presents a new modification of the isoniazid (INH) structure linked with different anilines via a carbonyl group obtained by two synthetic procedures and with N-substituted 5-(pyridine-4-yl)-1,3,4-oxadiazole-2-amines prepared by their cyclisation. All synthesised derivatives were characterised by IR, NMR, MS and elemental analyses and were evaluated in vitro for their antimycobacterial activity against Mycobacterium tuberculosis H₃₇Rv, Mycobacterium avium 330/88, Mycobacterium kansasii 235/80 and one clinical isolated strain of M. kansasii 6509/96. 2-Isonicotinoyl-N-(4-octylphenyl)hydrazinecarboxamide displayed an in vitro efficacy comparable to that of INH for M. tuberculosis with minimum inhibitory concentrations (MICs) of 1–2 μM. Among the halogenated derivatives, the best anti-tuberculosis activity was found for 2-isonicotinoyl-N-(2,4,6-trichlorophenyl)hydrazinecarboxamide (MIC = 4 μM). In silico modelling on the enoyl-acyl carrier protein reductase InhA confirmed that longer alkyl substituents are advantageous for the interactions and affinity to InhA. Most of the hydrazinecarboxamides, especially those derived from 4-alkylanilines, exhibited significant activity against INH-resistant nontuberculous mycobacteria.

Context-based preprocessing of molecular docking data

Background

Data preprocessing is a major step in data mining. In data preprocessing, several known techniques can be applied, or new ones developed, to improve data quality such that the mining results become more accurate and intelligible. Bioinformatics is one area with a high demand for generation of comprehensive models from large datasets. In this article, we propose a context-based data preprocessing approach to mine data from molecular docking simulation results. The test cases used a fully-flexible receptor (FFR) model of Mycobacterium tuberculosis InhA enzyme (FFR_InhA) and four different ligands.

Results

We generated an initial set of attributes as well as their respective instances. To improve this initial set, we applied two selection strategies. The first was based on our context-based approach while the second used the CFS (Correlation-based Feature Selection) machine learning algorithm. Additionally, we produced an extra dataset containing features selected by combining our context strategy and the CFS algorithm. To demonstrate the effectiveness of the proposed method, we evaluated its performance based on various predictive (RMSE, MAE, Correlation, and Nodes) and context (Precision, Recall and FScore) measures.

Conclusions

Statistical analysis of the results shows that the proposed context-based data preprocessing approach significantly improves predictive and context measures and outperforms the CFS algorithm. Context-based data preprocessing improves mining results by producing superior interpretable models, which makes it well-suited for practical applications in molecular docking simulations using FFR models.

Encoded library technology as a source of hits for the discovery and lead optimization of a potent and selective class of bactericidal direct inhibitors of Mycobacterium tuberculosis InhA

Tuberculosis (TB) is one of the world's oldest and deadliest diseases, killing a person every 20 s. InhA, the enoyl-ACP reductase from Mycobacterium tuberculosis, is the target of the frontline antitubercular drug isoniazid (INH). Compounds that directly target InhA and do not require activation by mycobacterial catalase peroxidase KatG are promising candidates for treating infections caused by INH resistant strains. The application of the encoded library technology (ELT) to the discovery of direct InhA inhibitors yielded compound 7 endowed with good enzymatic potency but with low antitubercular potency. This work reports the hit identification, the selected strategy for potency optimization, the structure-activity relationships of a hundred analogues synthesized, and the results of the in vivo efficacy studies performed with the lead compound 65.

An altered Mycobacterium tuberculosis metabolome induced by katG mutations resulting in isoniazid resistance

The most common form of drug resistance found in tuberculosis (TB)-positive clinical samples is monoresistance to isoniazid. Various genomics and proteomics studies to date have investigated this phenomenon; however, the exact mechanisms relating to how this occurs, as well as the implications of this on the TB-causing organisms function and structure, are only partly understood. Considering this, we followed a metabolomics research approach to identify potential new metabolic pathways and metabolite markers, which when interpreted in context would give a holistic explanation for many of the phenotypic characteristics associated with a katG mutation and the resulting isoniazid resistance in Mycobacterium tuberculosis. In order to achieve these objectives, gas chromatography-time of flight mass spectrometry (GCxGC-TOFMS)-generated metabolite profiles from two isoniazid-resistant strains were compared to a wild-type parent strain. Principal component analyses showed clear differentiation between the groups, and the metabolites best describing the separation between these groups were identified. It is clear from the data that due to a mutation in the katG gene encoding catalase, the isoniazid-resistant strains experience increased susceptibility to oxidative stress and have consequently adapted to this by upregulating the synthesis of a number of compounds involved in (i) increased uptake and use of alkanes and fatty acids as a source of carbon and energy and (ii) the synthesis of a number of compounds directly involved in reducing oxidative stress, including an ascorbic acid degradation pathway, which to date hasn't been proposed to exist in these organisms.