Probing mechanisms of resistance to the tuberculosis drug isoniazid: Conformational changes caused by inhibition of InhA, the enoyl reductase from Mycobacterium tuberculosis

In silico analyses of isoniazid and streptomycin resistance-associated mutations in Mycobacterium tuberculosis

Multi-drug resistant tuberculosis is categorised by the World Health Organisation (WHO) as a public health crisis. In silico techniques were used to probe the structural basis of Mycobacterium tuberculosis resistance to isoniazid and streptomycin. Isoniazid resistance-associated mutations in InhA were predicted to reduce the binding affinity of NADH to InhA, without affecting INH-NAD (competitive-inhibitor) binding. Perturbation of the mutated residues was predicted (with the AlloSigMA server) to modulate the free energy of allosteric modulation of key binding site residues F41, F149, Y158 and W222. These results suggest that allosteric modulation of the protein structure may be key to the mechanism by which isoniazid resistance-associated mutations act. Mutations in the methyltransferase glucose-inhibited division gene B (GidB) are associated with streptomycin resistance. Molecular docking was carried out to predict the structure of the GidB bound to its substrate (s-adenosyl methionine). The effects of streptomycin resistance-associated mutations in GidB on protein stability and substrate binding were predicted (using SDM and mCSM-lig). All GidB mutants were predicted to disfavour SAM binding.

Investigations into Isoniazid Treated Mycobacterium tuberculosis by Electrospray Mass Spectrometry Reveals New Insights into Its Lipid Composition

Many of the earlier studies involving the effect of isoniazid (INH) treatment have solely focused on the fatty acyl (FA) category of Mycobacterium tuberculosis (MTB) lipids. This motivated us with the major interest to examine the impact of INH on various other categories of MTB lipids. Towards this, we chose to interpret our mass spectral data (LC-ESI-MS) by a standalone software, MS-LAMP, in which "Mtb LipidDB" was integrated. Analysis by MS-LAMP revealed that INH treatment can alter the composition of "glycerolipids (GLs)" and "glycerophospholipids (GPLs)" categories of MTB lipids, in addition to the variations to FA category. Interpretation by "MycoMass" database yielded similar results as that of Mtb LipidDB, except that significant alterations to polyketides (PKs) category also were observed. Probing biosynthetic pathways of certain key lipids belonging to any of GLs, GPLs, and PKs categories can be attractive target(s) for drug discovery or can be useful to identify means to overcome drug resistance or to obtain insights into the causal factors of virulence. To the best of our knowledge, this is the first report hinting at the influence of INH on GLs, GPLs, and PKs of MTB.

Examining the role of protein structural dynamics in drug resistance in Mycobacterium tuberculosis

2D-IR spectroscopy reveals a role for protein structural dynamics in antimicrobial-resistance.

Antimicrobial resistance represents a growing global health problem. The emergence of novel resistance mechanisms necessitates the development of alternative approaches to investigate the molecular fundamentals of resistance, leading ultimately to new strategies for counteracting them. To gain deeper insight into antibiotic–target interactions, the binding of the frontline anti-tuberculosis drug isoniazid (INH) to a target enzyme, InhA, from Mycobacterium tuberculosis was studied using ultrafast two-dimensional infrared (2D-IR) spectroscopy and molecular simulations. Comparing wild-type InhA with a series of single point mutations, it was found that binding of the INH–NAD inhibitor to susceptible forms of the enzyme increased the vibrational coupling between residues located in the Rossmann fold co-factor binding site of InhA and suppressed dynamic fluctuations of the enzyme structure. The effect correlated with biochemical assay data, being reduced in the INH-resistant S94A mutant and absent in the biochemically-inactive P193A control. Molecular dynamics simulations and calculations of inter–residue couplings indicate that the changes in coupling and dynamics are not localised to the co-factor binding site, but permeate much of the protein. We thus propose that the resistant S94A mutation circumvents subtle changes in global structural dynamics caused by INH upon binding to the wild-type enzyme that may impact upon the formation of important protein–protein complexes in the fatty acid synthase pathway of M. tuberculosis.

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.

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.

A novel three-component system-based regulatory model for D-xylose sensing and transport in Clostridium beijerinckii

D-Xylose is the most abundant fermentable pentose in nature and can serve as a carbon source for many bacterial species. Since D-xylose constitutes the major component of hemicellulose, its metabolism is important for lignocellulosic biomass utilization. Here, we report a six-protein module for D-xylose signaling, uptake and regulation in solvent-producing Clostridium beijerinckii. This module consists of a novel 'three-component system' (a putative periplasmic ABC transporter substrate-binding protein XylFII and a two-component system LytS/YesN) and an ABC-type D-xylose transporter XylFGH. Interestingly, we demonstrate that, although XylFII harbors a transmembrane domain, it is not involved in D-xylose transport. Instead, XylFII acts as a signal sensor to assist the response of LytS/YesN to extracellular D-xylose, thus enabling LytS/YesN to directly activate the transcription of the adjacent xylFGH genes and thereby promote the uptake of D-xylose. To our knowledge, XylFII is a novel single transmembrane sensor that assists two-component system to respond to extracellular sugar molecules. Also of significance, this 'three-component system' is widely distributed in Firmicutes, indicating that it may play a broad role in this bacterial phylum. The results reported here provide new insights into the regulatory mechanism of D-xylose sensing and transport in bacteria.

On the chemical mechanism of succinic semialdehyde dehydrogenase (GabD1) from Mycobacterium tuberculosis

Succinic semialdehyde dehydrogenases (SSADHs) are ubiquitous enzymes that catalyze the NAD(P)+-coupled oxidation of succinic semialdehyde (SSA) to succinate, the last step of the γ-aminobutyrate shunt. Mycobacterium tuberculosis encodes two paralogous SSADHs (gabD1 and gabD2). Here, we describe the first mechanistic characterization of GabD1, using steady-state kinetics, pH-rate profiles, ¹H NMR, and kinetic isotope effects. Our results confirmed SSA and NADP+ as substrates and demonstrated that a divalent metal, such as Mg²+, linearizes the time course. pH-rate studies failed to identify any ionizable groups with pK(a) between 5.5 and 10 involved in substrate binding or rate-limiting chemistry. Primary deuterium, solvent and multiple kinetic isotope effects revealed that nucleophilic addition to SSA is very fast, followed by a modestly rate-limiting hydride transfer and fast thioester hydrolysis. Proton inventory studies revealed that a single proton is associated with the solvent-sensitive rate-limiting step. Together, these results suggest that product dissociation and/or conformational changes linked to it are rate-limiting. Using structural information for the human homolog enzyme and ¹H NMR, we further established that nucleophilic attack takes place at the Si face of SSA, generating a thiohemiacetal with S stereochemistry. Deuteride transfer to the Pro-R position in NADP+ generates the thioester intermediate and [4A-²H, 4B-¹H] NADPH. A chemical mechanism based on these data and the structural information available is proposed.

Mining flexible-receptor docking experiments to select promising protein receptor snapshots

Background

Molecular docking simulation is the Rational Drug Design (RDD) step that investigates the affinity between protein receptors and ligands. Typically, molecular docking algorithms consider receptors as rigid bodies. Receptors are, however, intrinsically flexible in the cellular environment. The use of a time series of receptor conformations is an approach to explore its flexibility in molecular docking computer simulations, but it is extensively time-consuming. Hence, selection of the most promising conformations can accelerate docking experiments and, consequently, the RDD efforts.

Results

We previously docked four ligands (NADH, TCL, PIF and ETH) to 3,100 conformations of the InhA receptor from M. tuberculosis. Based on the receptor residues-ligand distances we preprocessed all docking results to generate appropriate input to mine data. Data preprocessing was done by calculating the shortest interatomic distances between the ligand and the receptor’s residues for each docking result. They were the predictive attributes. The target attribute was the estimated free-energy of binding (FEB) value calculated by the AutodDock3.0.5 software. The mining inputs were submitted to the M5P model tree algorithm. It resulted in short and understandable trees. On the basis of the correlation values, for NADH, TCL and PIF we obtained more than 95% correlation while for ETH, only about 60%. Post processing the generated model trees for each of its linear models (LMs), we calculated the average FEB for their associated instances. From these values we considered a LM as representative if its average FEB was smaller than or equal the average FEB of the test set. The instances in the selected LMs were considered the most promising snapshots. It totalized 1,521, 1,780, 2,085 and 902 snapshots, for NADH, TCL, PIF and ETH respectively.

Conclusions

By post processing the generated model trees we were able to propose a criterion of selection of linear models which, in turn, is capable of selecting a set of promising receptor conformations. As future work we intend to go further and use these results to elaborate a strategy to preprocess the receptors 3-D spatial conformation in order to predict FEB values. Besides, we intend to select other compounds, among the million catalogued, that may be promising as new drug candidates for our particular protein receptor target.

A novel interaction linking the FAS-II and phthiocerol dimycocerosate (PDIM) biosynthetic pathways

The fatty acid biosynthesis (FAS-II) pathway in Mycobacterium tuberculosis generates long chain fatty acids that serve as the precursors to mycolic acids, essential components of the mycobacterial cell wall. Enzymes in the FAS-II pathway are thought to form one or more noncovalent multi-enzyme complexes within the cell, and a bacterial two-hybrid screen was used to search for missing components of the pathway and to furnish additional data on interactions involving these enzymes in vivo. Using the FAS-II beta-ketoacyl synthase, KasA, as bait, an extensive bacterial two-hybrid screen of a M. tuberculosis genome fragment library unexpectedly revealed a novel interaction between KasA and PpsB as well as PpsD, two polyketide modules involved in the biosynthesis of the virulence lipid phthiocerol dimycocerosate (PDIM). Sequence analysis revealed that KasA interacts with PpsB and PpsD in the region of the acyl carrier domain of each protein, raising the possibility that lipids could be transferred between the FAS-II and PDIM biosynthetic pathways. Subsequent studies utilizing purified proteins and radiolabeled lipids revealed that fatty acids loaded onto PpsB were transferred to KasA and also incorporated into long chain fatty acids synthesized using a Mycobacterium smegmatis lysate. These data suggest that in addition to producing PDIMs, the growing phthiocerol product can also be shuttled into the FAS-II pathway via KasA as an entry point for further elongation. Interactions between these biosynthetic pathways may exist as a simple means to increase mycobacterial lipid diversity, enhancing functionality and the overall complexity of the cell wall.

Arylamine N-acetyltransferases in mycobacteria

Polymorphic Human arylamine N-acetyltransferase (NAT2) inactivates the anti-tubercular drug isoniazid by acetyltransfer from acetylCoA. There are active NAT proteins encoded by homologous genes in mycobacteria including M. tuberculosis, M. bovis BCG, M. smegmatis and M. marinum. Crystallographic structures of NATs from M. smegmatis and M. marinum, as native enzymes and with isoniazid bound share a similar fold with the first NAT structure, Salmonella typhimurium NAT. There are three approximately equal domains and an active site essential catalytic triad of cysteine, histidine and aspartate in the first two domains. An acetyl group from acetylCoA is transferred to cysteine and then to the acetyl acceptor e.g. isoniazid. M. marinum NAT binds CoA in a more open mode compared with CoA binding to human NAT2. The structure of mycobacterial NAT may promote its role in synthesis of cell wall lipids, identified through gene deletion studies. NAT protein is essential for survival of M. bovis BCG in macrophage as are the proteins encoded by other genes in the same gene cluster (hsaA-D). HsaA-D degrade cholesterol, essential for mycobacterial survival inside macrophage. Nat expression remains to be fully understood but is co-ordinated with hsaA-D and other stress response genes in mycobacteria. Amide synthase genes in the streptomyces are also nat homologues. The amide synthases are predicted to catalyse intramolecular amide bond formation and creation of cyclic molecules, e.g. geldanamycin. Lack of conservation of the CoA binding cleft residues of M. marinum NAT suggests the amide synthase reaction mechanism does not involve a soluble CoA intermediate during amide formation and ring closure.