NAD+-dependent DNA Ligase (Rv3014c) from Mycobacterium tuberculosis. Crystal structure of the adenylation domain and identification of novel inhibitors

DNA ligases utilize either ATP or NAD+ as cofactors to catalyze the formation of phosphodiester bonds in nicked DNA. Those utilizing NAD+ are attractive drug targets because of the unique cofactor requirement for ligase activity. We report here the crystal structure of the adenylation domain of the Mycobacterium tuberculosis NAD+-dependent ligase with bound AMP. The adenosine nucleoside moiety of AMP adopts a syn-conformation. The structure also captures a new spatial disposition between the two subdomains of the adenylation domain. Based on the crystal structure and an in-house compound library, we have identified a novel class of inhibitors for the enzyme using in silico docking calculations. The glycosyl ureide-based inhibitors were able to distinguish between NAD+- and ATP-dependent ligases as evidenced by in vitro assays using T4 ligase and human DNA ligase I. Moreover, assays involving an Escherichia coli strain harboring a temperature-sensitive ligase mutant and a ligase-deficient Salmonella typhimurium strain suggested that the bactericidal activity of the inhibitors is due to inhibition of the essential ligase enzyme. The results can be used as the basis for rational design of novel antibacterial agents.     

Mutational Analysis of Bacterial NAD^＋-dependent DNA Ligase： Role of Motif IV in Ligation Catalysis

The bacterial DNA ligase as a multiple domain protein is involved in DNA replication, repair and recombination. Its catalysis of ligation can be divided into three steps. To delineate the roles of amino acid residues in motif IV in ligation catalysis, site-directed mutants were constructed in a bacterial NAD^＋-dependent DNA ligase from Thermus sp. TAK16D. It was shown that four conserved residues （D286, G287, V289 and K291） in motif IV had significant roles on the overall ligation. Under single turnover conditions, the observed apparent rates of D286E, G287A, V289I and K291R mutants were clearly reduced compared with that of WT ligase on both match and mismatch nicked substrates. The effects of D286E mutation on overall ligation may not only be ascribed to the third step. The G287A mutation has a major effect on the second step. The effects of V289I and K291R mutation on overall ligation are not on the third step, perhaps other aspects, such as conformation change of ligase protein in ligation catalysis, are involved. Moreover, the amino acid substitutions of above four residues were more sensitive on mismatch nicked substrate, indicating an enhanced ligation fidelity.     

NAD+-dependent DNA ligases of Mycobacterium tuberculosis and Streptomyces coelicolor

Sequencing of the genomes of Mycobacterium tuberculosis H37Rv and Streptomyces coelicolor A3(2) identified putative genes for an NAD(+)-dependent DNA ligase. We have cloned both open reading frames and overexpressed the protein products in Escherichia coli. In vitro biochemical assays confirm that each of these proteins encodes a functional DNA ligase that uses NAD(+) as its cofactor. Expression of either protein is able to complement E. coli GR501, which carries a temperature-sensitive mutation in ligA. Thus, in vitro and in vivo analyses confirm predictions that ligA genes from M. tuberculosis and S. coelicolor are NAD(+)-dependent DNA ligases.     

PPE protein (Rv3425) from DNA segment RD11 of Mycobacterium tuberculosis: a novel immunodominant antigen of Mycobacterium tuberculosis induces humoral and cellular immune responses in mice

Subtractive DNA hybridization of pathogenic M. bovis and BCG, and comparative genome-wide DNA microarray analysis of M. tuberculosis H37Rv and BCG identified several RD, designated as RD1 to RD16, between M. tuberculosis and M. bovis on the one hand and BCG on the other. These regions cover 108 ORF of M. tuberculosis H37Rv, and are deleted from all 13 BCG sub-strains currently used as anti-tuberculosis vaccines in different parts of the world. In this study, we evaluated cellular and humoral immune response in C57BL/6 mice immunized with the PPE protein Rv3425, encoded by an ORF found in RD11 of M. tuberculosis. Rv3425 protein induced an increased Th1/Th2 type immune response in mice, characterized by an elevated concentration of IFN-gamma in antigen stimulated splenocyte culture and a strong IgG(1) antibody response. These results provide evidence on the immunogenicity of the PPE protein Rv3425 which, together with its reported immunodominant characteristics, imply that it may be a candidate for development of a vaccine for the control of TB.     

In vitro and in vivo identification of a novel cytotoxic T lymphocyte epitope from Rv3425 of Mycobacterium tuberculosis

The identification of novel cytotoxic T lymphocyte (CTL) epitopes is important to analysis of the involvement of CD8(+) T cells in Mycobacterium tuberculosis infection as well as to the development of peptide vaccines. In this study, a novel CTL epitope from region of difference 11 encoded antigen Rv3425 was identified. Epitopes were predicted by the reversal immunology approach. Rv3425-p118 (LIASNVAGV) was identified as having relatively strong binding affinity and stability towards the HLA-A*0201 molecule. Peripheral blood mononuclear cells pulsed by this peptide were able to release interferon-γ in healthy donors (HLA-A*02(+) purified protein derivative(+)). In cytotoxicity assays in vitro and in vivo, Rv3425-p118 induced CTLs to specifically lyse the target cells. Therefore, this epitope could provide a subunit component for designing vaccines against Mycobacterium tuberculosis.     

Structural and Biochemical Insight into the Mechanism of Rv2837c from Mycobacterium tuberculosis as a c-di-NMP Phosphodiesterase

The intracellular infections of Mycobacterium tuberculosis, which is the causative agent of tuberculosis, are regulated by many cyclic dinucleotide signaling. Rv2837c from M. tuberculosis is a soluble, stand-alone DHH-DHHA1 domain phosphodiesterase that down-regulates c-di-AMP through catalytic degradation and plays an important role in M. tuberculosis infections. Here, we report the crystal structure of Rv2837c (2.0 Å), and its complex with hydrolysis intermediate 5′-pApA (2.35 Å). Our structures indicate that both DHH and DHHA1 domains are essential for c-di-AMP degradation. Further structural analysis shows that Rv2837c does not distinguish adenine from guanine, which explains why Rv2837c hydrolyzes all linear dinucleotides with almost the same efficiency. We observed that Rv2837c degraded other c-di-NMPs at a lower rate than it did on c-di-AMP. Nevertheless, our data also showed that Rv2837c significantly decreases concentrations of both c-di-AMP and c-di-GMP in vivo. Our results suggest that beside its major role in c-di-AMP degradation Rv2837c could also regulate c-di-GMP signaling pathways in bacterial cell.     

Crystal structure of calcium dodecin (Rv0379), from Mycobacterium tuberculosis with a unique calcium-binding site

In eukaryotes, calcium-binding proteins play a pivotal role in diverse cellular processes, and recent findings suggest similar roles for bacterial proteins at different stages in their life cycle. Here, we report the crystal structure of calcium dodecin, Rv0379, from Mycobacterium tuberculosis with a dodecameric oligomeric assembly and a unique calcium-binding motif. Structure and sequence analysis were used to identify orthologs of Rv0379 with different ligand-binding specificity.     

Inhibitor binding studies of Mycobacterium tuberculosis MraY (Rv2156c): Insights from molecular modeling,docking, and simulation studies

Abstract

Tuberculosis (TB) is a contagious disease caused by Mycobacterium tuberculosis (M.tb) or tubercule bacillus, and H37Rv is the most studied clinical strain. The recent development of resistance to existing drugs is a global health-care challenge to control and cure TB. Hence, there is a critical need to discover new drug targets in M.tb. The members of peptidoglycan biosynthesis pathway are attractive target proteins for antibacterial drug development. We have performed in silico analysis of M.tb MraY (Rv2156c) integral membrane protein and constructed the three-dimensional (3D) structure model of M.tb MraY based on homology modeling method. The validated model was complexed with antibiotic muraymycin D2 (MD2) and was used to generate structure-based pharmacophore model (e-pharmacophore). High-throughput virtual screening (HTVS) of Asinex database and molecular docking of hits was performed to identify the potential inhibitors based on their mode of interactions with the key residues involved in M.tb MraY–MD2 binding. The validation of these molecules was performed using molecular dynamics (MD) simulations for two best identified hit molecules complexed with M.tb MraY in the lipid bilayer, dipalmitoylphosphatidyl-choline (DPPC) membrane. The results indicated the stability of the complexes formed and retained non-bonding interactions similar to MD2. These findings may help in the design of new inhibitors to M.tb MraY involved in peptidoglycan biosynthesis.     

Biochemical characterization and novel inhibitor identification of Mycobacterium tuberculosis Endonuclease VIII 2 (Rv3297)

Nei2 (Rv3297) is a DNA Base Excision Repair (BER) glycosylase that is essential for survival of Mycobacterium tuberculosis in primates. We show that MtbNei2 is a bifunctional glycosylase that specifically acts on oxidized pyrimidine-containing single-stranded, double-stranded, 5’/3’ fork and bubble DNA substrates. MtbNei2 possesses Uracil DNA glycosylase activity unlike E. coli Nei. Mutational studies demonstrate that Pro2 and Glu3 located in the active site are essential for glycosylase activity of MtbNei2. Mutational analysis demonstrated that an unstructured C-terminal zinc finger domain that was important for activity in E. coli Nei and Fpg, was not required for the glycosylase activity of MtbNei2. Lastly, we screened the NCI natural product compound database and identified three natural product inhibitors with IC50 values ranging between 41.8 μM-92.7 μM against MtbNei2 in in vitro inhibition assays. Surface Plasmon Resonance (SPR) experiments showed that the binding affinity of the best inhibitor, NSC31867, was 74 nM. The present results set the stage for exploiting this important target in developing new therapeutic strategies that target Mycobacterial BER.     

Crystal structure of a putative pyridoxine 5′‐phosphate oxidase (Rv2607) from Mycobacterium tuberculosis

The three‐dimensional structure of Rv2607, a putative pyridoxine 5′‐phosphate oxidase (PNPOx) from Mycobacterium tuberculosis, has been determined by X‐ray crystallography to 2.5 Å resolution. Rv2607 has a core domain similar to known PNPOx structures with a flavin mononucleotide (FMN) cofactor. Electron density for two FMN at the dimer interface is weak despite the bright yellow color of the protein solution and crystal. The shape and size of the putative binding pocket is markedly different from that of members of the PNPOx family, which may indicate some significant changes in the FMN binding mode of this protein relative to members of the family. Proteins 2006. © 2005 Wiley‐Liss, Inc.     

Adenylation-dependent conformation and unfolding pathways of the NAD+-dependent DNA ligase from the thermophile Thermus scotoductus

In the last few years, an increased attention has been focused on NAD(+)-dependent DNA ligases. This is mostly due to their potential use as antibiotic targets, because effective inhibition of these essential enzymes would result in the death of the bacterium. However, development of an efficient drug requires that the conformational modifications involved in the catalysis of NAD(+)-dependent DNA ligases are understood. From this perspective, we have investigated the conformational changes occurring in the thermophilic Thermus scotoductus NAD(+)-DNA ligase upon adenylation, as well as the effect of cofactor binding on protein resistance to thermal and chemical (guanidine hydrochloride) denaturation. Our results indicate that cofactor binding induces conformational rearrangement within the active site and promotes a compaction of the enzyme. These data support an induced "open-closure" process upon adenylation, leading to the formation of the catalytically active enzyme that is able to bind DNA. These conformational changes are likely to be associated with the protein function, preventing the formation of nonproductive complexes between deadenylated ligases and DNA. In addition, enzyme adenylation significantly increases resistance of the protein to thermal denaturation and GdmCl-induced unfolding, establishing a thermodynamic link between ligand binding and increased conformational stability. Finally, chemical unfolding of deadenylated and adenylated enzyme is accompanied by accumulation of at least two equilibrium intermediates, the molten globule and premolten globule states. Maximal populations of these intermediates are shifted toward higher GdmCl concentrations in the case of the adenylated ligase. These data provide further insights into the properties of partially folded intermediates.     

Direct evidence for a glutamate switch necessary for substrate recognition: crystal structures of lysine epsilon-aminotransferase (Rv3290c) from Mycobacterium tuberculosis H37Rv

Lysine epsilon-aminotransferase (LAT) is a PLP-dependent enzyme that is highly up-regulated in nutrient-starved tuberculosis models. It catalyzes an overall reaction involving the transfer of the epsilon-amino group of L-lysine to alpha-ketoglutarate to yield L-glutamate and alpha-aminoadipate-delta-semialdehyde. We have cloned and characterized the enzyme from Mycobacterium tuberculosisH37Rv. We report here the crystal structures of the enzyme, the first from any source, in the unliganded form, external aldimine with L-lysine, with bound PMP and with its C5 substrate alpha-ketoglutarate. In addition to interaction details in the active site, the structures reveal a Glu243 "switch" through which the enzyme changes substrate specificities. The unique substrate L-lysine is recognized specifically when Glu243 maintains a salt-bridge with Arg422. On the other hand, the binding of the common C5 substrates L-glutamate and alpha-ketoglutarate is enabled when Glu243 switches away and unshields Arg422. The structures reported here, sequence conservation and earlier mutational studies suggest that the "glutamate switch" is an elegant solution devised by a subgroup of fold type I aminotransferases for recognition of structurally diverse substrates in the same binding site and provides for reaction specificity.     

Docking- and pharmacophore-based virtual screening for the identification of novel Mycobacterium tuberculosis protein tyrosine phosphatase B (MptpB) inhibitor with a thiobarbiturate scaffold

Mycobacterium tuberculosis (Mtb) protein tyrosine phosphatase B (MptpB) is an important virulence factor for Mtb that contributes to survival of the bacteria in macrophages. The absence of a human ortholog makes MptpB an attractive target for new therapeutics to treat tuberculosis. MptpB inhibitors could be an effective treatment to overcome emerging TB drug resistance. Adopting a structure-based virtual screening strategy, we successfully identified thiobarbiturate-based drug-like MptpB inhibitor 15 with an IC₅₀ of 22.4 μM, and as a non-competitive inhibitor with a K_i of 24.7 μM. Importantly, not only did it exhibit moderate cell membrane permeability, compound 15 also displayed potent inhibition of intracellular TB growth in the macrophage, making it an excellent lead compound for anti-TB drug discovery. To the best of our knowledge, this novel thiobarbiturate is the first class of MptpB inhibitor reported so far that leveraged docking- and pharmacophore-based virtual screening approaches. The results of preliminary structure-activity relationship demonstrated that compound 15 identified herein was not a singleton and may inspire the design of novel selective and drug-like MptpB inhibitors.     

Arcyptera fusca and Arcyptera tornosi repetitive DNA families: whole-comparative genomic hybridization (W-CGH) as a novel approach to the study of satellite DNA libraries

Whole-comparative genomic hybridization (W-CGH) has been used to exemplify a simple methodology which allows identifying and mapping whole genome differences for highly repetitive DNA sequences between two related species of unknown genomic background. The use of this technique to the species binomy Arcyptera fusca/Arcyptera tornosi has allowed the identification of different DNA families mainly concentrated within the para-/peri-centromeric and distal heterochromatic regions of different chromosomes, which are differentially expanded in both genomes. Additionally, W-CGH allowed chromosome mapping of particular euchromatic regions immersed in the chromosome arms which have been affected by processes of DNA amplification and losses. A molecular approach was also conducted to analyse satellite DNA families in these species. We have found three different families showing an unequal representation in both species. Two of these families showed a centromeric location (EcoRV-390CEN and Sau3A-419CEN), whereas the last one was located at distal heterochromatic regions (Sau3A-197TEL). As A. fusca is a widely distributed species represented in most European high mountains, whereas A. tornosi is an endemic species represented in the Iberian Peninsula, the differences and resemblances reported here offer a good basis to support a close evolutionary relationship between both of the actually isolated species. Finally, W-CGH allowed identification of an asynchronic pattern of heterochromatin condensation through early prophase (characteristic in both species) which is uncommon or probably has been poorly analysed within classical early condensing chromosome domains through meiosis. The congruence of the obtained cytological and molecular results is analysed in light of the ancestral genome relationship between both species.     

A novel interaction partner for the C-terminus of Arabidopsis thaliana plasma membrane H+-ATPase (AHA1 isoform): site and mechanism of action on H+-ATPase activity differ from those of 14-3-3 proteins

Using the two-hybrid technique we identified a novel protein whose N-terminal 88 amino acids (aa) interact with the C-terminal regulatory domain of the plasma membrane (PM) H+-ATPase from Arabidopsis thaliana (aa 847-949 of isoform AHA1). The corresponding gene has been named Ppi1 for Proton pump interactor 1. The encoded protein is 612 aa long and rich in charged and polar residues, except for the extreme C-terminus, where it presents a hydrophobic stretch of 24 aa. Several genes in the A. thaliana genome and many ESTs from different plant species share significant similarity (50-70% at the aa level over stretches of 200-600 aa) to Ppi1. The PPI1 N-terminus, expressed in bacteria as a fusion protein with either GST or a His-tag, binds the PM H+-ATPase in overlay experiments. The same fusion proteins and the entire coding region fused to GST stimulate H+-ATPase activity. The effect of the His-tagged peptide is synergistic with that of fusicoccin (FC) and of tryptic removal of a C-terminal 10 kDa fragment. The His-tagged peptide binds also the trypsinised H+-ATPase. Altogether these results indicate that PPI1 N-terminus is able to modulate the PM H+-ATPase activity by binding to a site different from the 14-3-3 binding site and is located upstream of the trypsin cleavage site.     

Synthesis and enantiopreferential DNA-binding profile of late 3d transition metal R- and S-enantiomeric complexes derived from N,N-bis-(1-benzyl-2-ethoxyethane): Validation of R-enantiomer of copper(II) complex as a human topoisomerase II inhibitor

To evaluate the biological preference of chiral drug candidates for molecular target DNA, new potential metal‐based chemotherapeutic agents 1 , 1a , 1b , 2 , 2a , 2b , 3 , 3a , 3b of late 3d transition metals Ni(II), Cu(II), and Zn(II), respectively, derived from (R)‐ and (S)‐2‐amino‐2‐phenylethanol with  CH₂ CH₂ linker were synthesized and thoroughly characterized. Interaction studies of 1 , 1a , 1b , 2 , 2a , 2b , 3 , 3a , 3b with calf thymus DNA in Tris buffer were studied by electronic absorption titrations, luminescence titrations, cyclic voltammetry, and circular dichroism. The results reveal that the extent of DNA binding of R‐enantiomer of copper 1a was highest in comparison to rest of the complexes via electrostatic interaction mode. The nuclease activity of 1a , 1b with supercoiled pBR322 DNA was further examined by gel electrophoresis, which reveals that complex 1a exhibits a remarkable DNA cleavage activity (concentration dependent) with pBR322DNA, and the cleavage activity of both enantiomers of complex 1 was significantly enhanced in the presence of activators. The activating efficiency follows the order Asc > H₂O₂ > MPA for 1a , and reverse order was observed for 1b , because of the differences in enantioselectivity and conformation. Further, it was observed that cleavage reaction involves singlet oxygen species and superoxide radicals via oxidative cleavage mechanism. In addition, complex 1a exhibits significant inhibitory effects on the topoisomerase II (topo II) activity at a very low concentration ∼24 μM, which suggest that complex 1a is indeed catalytic inhibitor or (poison) of human topo II. Chirality 2011 © 2011 Wiley‐Liss, Inc.