Protein targets for structure-based anti-Mycobacterium tuberculosis drug discovery

Inhibitors of Aspartate Transcarbamoylase Inhibit Mycobacterium tuberculosis Growth

Aspartate transcarbamoylase (ATCase) plays a key role in the second step of de novo pyrimidine biosynthesis in eukaryotes and has been proposed to be a target to suppress cell proliferation in E. coli, human cells and the malarial parasite. We hypothesized that a library of ATCase inhibitors developed for malarial ATCase (PfATCase) may also contain inhibitors of the tubercular ATCase and provide a similar inhibition of cellular proliferation. Of the 70 compounds screened, 10 showed single-digit micromolar inhibition in an in vitro activity assay and were tested for their effect on M. tuberculosis cell growth in culture. The most promising compound demonstrated a MIC90 of 4 μM. A model of MtbATCase was generated using the experimental coordinates of PfATCase. In silico docking experiments showed this compound can occupy a similar allosteric pocket on MtbATCase to that seen on PfATCase, explaining the observed species selectivity seen for this compound series.

Synthesis of Novel Derivatives of 5,6,7,8-Tetrahydroquinazolines Using α-Aminoamidines and In Silico Screening of Their Biological Activity

α-Aminoamidines are promising reagents for the synthesis of a diverse family of pyrimidine ring derivatives. Here, we demonstrate the use of α-aminoamidines for the synthesis of a new series of 5,6,7,8-tetrahydroquinazolines by their reaction with bis-benzylidene cyclohexanones. The reaction occurs in mild conditions and is characterized by excellent yields. It has easy workup, as compared to the existing methods of tetrahydroquinazoline preparation. Newly synthesized derivatives of 5,6,7,8-tetrahydroquinazoline bear protecting groups at the C2-tert-butyl moiety of a quinazoline ring, which can be easily cleaved, opening up further opportunities for their functionalization. Moreover, molecular docking studies indicate that the synthesized compounds reveal high binding affinity toward some essential enzymes of Mycobacterial tuberculosis, such as dihydrofolate reductase (DHFR), pantothenate kinase (MtPanK), and FAD-containing oxidoreductase DprE1 (MtDprE1), so that they may be promising candidates for the molecular design and the development of new antitubercular agents against multidrug-resistant strains of the Tubercle bacillus. Finally, the high inhibition activity of the synthesized compounds was also predicted against β-glucosidase, suggesting a novel tetrahydroquinazoline scaffold for the treatment of diabetes.

An Outline of the Latest Crystallographic Studies on Inhibitor-Enzyme Complexes for the Design and Development of New Therapeutics against Tuberculosis

The elucidation of the structure of enzymes and their complexes with ligands continues to provide invaluable insights for the development of drugs against many diseases, including bacterial infections. After nearly three decades since the World Health Organization’s (WHO) declaration of tuberculosis (TB) as a global health emergency, Mycobacterium tuberculosis (Mtb) continues to claim millions of lives, remaining among the leading causes of death worldwide. In the last years, several efforts have been devoted to shortening and improving treatment outcomes, and to overcoming the increasing resistance phenomenon. The structural elucidation of enzyme-ligand complexes is fundamental to identify hot-spots, define possible interaction sites, and elaborate strategies to develop optimized molecules with high affinity. This review offers a critical and comprehensive overview of the most recent structural information on traditional and emerging mycobacterial enzymatic targets. A selection of more than twenty enzymes is here discussed, with a special emphasis on the analysis of their binding sites, the definition of the structure–activity relationships (SARs) of their inhibitors, and the study of their main intermolecular interactions. This work corroborates the potential of structural studies, substantiating their relevance in future anti-mycobacterial drug discovery and development efforts.

Tuberculosis: An Overview of the Immunogenic Response, Disease Progression, and Medicinal Chemistry Efforts in the Last Decade toward the Development of Potential Drugs for Extensively Drug-Resistant Tuberculosis Strains

Tuberculosis (TB) is a slow growing, potentially debilitating disease that has plagued humanity for centuries and has claimed numerous lives across the globe. Concerted efforts by researchers have culminated in the development of various strategies to combat this malady. This review aims to raise awareness of the rapidly increasing incidences of multidrug-resistant (MDR) and extensively drug-resistant (XDR) tuberculosis, highlighting the significant modifications that were introduced in the TB treatment regimen over the past decade. A description of the role of pathogen-host immune mechanisms together with strategies for prevention of the disease is discussed. The struggle to develop novel drug therapies has continued in an effort to reduce the treatment duration, improve patient compliance and outcomes, and circumvent TB resistance mechanisms. Herein, we give an overview of the extensive medicinal chemistry efforts made during the past decade toward the discovery of new chemotypes, which are potentially active against TB-resistant strains.

Probing the dual inhibitory mechanisms of novel thiophenecarboxamide derivatives against Mycobacterium tuberculosis PyrG and PanK: an insight from biomolecular modeling study

The growing occurrence of multidrug-resistant (MDR) and extensively drug-resistant (XDR) Mycobacterium tuberculosis (Mtb) strains underscores an urgent need for new antibiotics. The development of more bioactive antibiotics against drug-resistant organisms with a different mode of action could be a game-changer for the cure and eradication of tuberculosis (TB). Pantothenate Kinase (PanK) and CTP synthetase (PyrG) are both essential for RNA, DNA, and Lipids biosynthesis pathways. Given the extensive knowledge on these biosynthesis pathways inhibition of Mtb growth and survival, these enzymes present a fascinating opportunity for anti-mycobacterial drug discovery. Recently, it was experimentally established that the active metabolite 11426026 of compound 7947882 (a prodrug activated by EthA monooxygenase, 5-methyl-N-(4-nitrophenyl) thiophene-2-carboxamide) inhibits the activities of PyrG and PanK to indicate novel multitarget therapy aimed at discontinuing Mtb growth. However, the molecular mechanisms of their selective inhibition remain subtle. In this work, molecular dynamics simulations were employed to investigate the inhibitory mechanism as well as the selectivity impact of the active metabolite inhibitor of these enzymes. Computational modeling of the studied protein-ligand systems reveals that the active metabolite can potentially inhibit both PanK and PyrG, thereby creating a pathway as a double target approach in tuberculosis treatment.Communicated by Ramaswamy H. Sarma.

Structural basis for the broad substrate specificity of two acyl-CoA dehydrogenases FadE5 from mycobacteria

FadE, an acyl-CoA dehydrogenase, introduces unsaturation to carbon chains in lipid metabolism pathways. Here, we report that FadE5 from Mycobacterium tuberculosis (MtbFadE5) and Mycobacterium smegmatis (MsFadE5) play roles in drug resistance and exhibit broad specificity for linear acyl-CoA substrates but have a preference for those with long carbon chains. Here, the structures of MsFadE5 and MtbFadE5, in the presence and absence of substrates, have been determined. These reveal the molecular basis for the broad substrate specificity of these enzymes. FadE5 interacts with the CoA region of the substrate through a large number of hydrogen bonds and an unusual π-π stacking interaction, allowing these enzymes to accept both short- and long-chain substrates. Residues in the substrate binding cavity reorient their side chains to accommodate substrates of various lengths. Longer carbon-chain substrates make more numerous hydrophobic interactions with the enzyme compared with the shorter-chain substrates, resulting in a preference for this type of substrate.

Using a Heat Diffusion Model to Detect Potential Drug Resistance Genes of Mycobacterium tuberculosis

BACKGROUND

Tuberculosis (TB), caused by Mycobacterium tuberculosis (Mtb), is one of the oldest known and most dangerous diseases. Although the spread of TB was controlled in the early 20th century using antibiotics and vaccines, TB has again become a threat because of increased drug resistance. There is still a lack of effective treatment regimens for a person who is already infected with multidrug-resistant Mtb (MDR-Mtb) or extensively drug-resistant Mtb (XDRMtb). In the past decades, many research groups have explored the drug resistance profiles of Mtb based on sequence data by GWAS, which identified some mutations that were significantly linked with drug resistance, and attempted to explain the resistance mechanisms. However, they mainly focused on several significant mutations in drug targets (e.g. rpoB, katG). Some genes which are potentially associated with drug resistance may be overlooked by the GWAS analysis.

OBJECTIVE

In this article, our motivation is to detect potential drug resistance genes of Mtb using a heat diffusion model.

METHODS

All sequencing data, which contained 127 samples of Mtb, i.e. 34 ethambutol-, 65 isoniazid-, 53 rifampicin- and 45 streptomycin-resistant strains. The raw sequence data were preprocessed using Trimmomatic software and aligned to the Mtb H37Rv reference genome using Bowtie2. From the resulting alignments, SAMtools and VarScan were used to filter sequences and call SNPs. The GWAS was performed by the PLINK package to obtain the significant SNPs, which were mapped to genes. The P-values of genes calculated by GWAS were transferred into a heat vector. The heat vector and the Mtb protein-protein interactions (PPI) derived from the STRING database were inputted into the heat diffusion model to obtain significant subnetworks by HotNet2. Finally, the most significant (P < 0.05) subnetworks associated with different phenotypes were obtained. To verify the change of binding energy between the drug and target before and after mutation, the method of molecular dynamics simulation was performed using the AMBER software.

RESULTS

We identified significant subnetworks in rifampicin-resistant samples. Excitingly, we found rpoB and rpoC, which are drug targets of rifampicin. From the protein structure of rpoB, the mutation location was extremely close to the drug binding site, with a distance of only 3.97 Å. Molecular dynamics simulation revealed that the binding energy of rpoB and rifampicin decreased after D435V mutation. To a large extent, this mutation can influence the affinity of drug-target binding. In addition, topA and pyrG were reported to be linked with drug resistance, and might be new TB drug targets. Other genes that have not yet been reported are worth further study.

CONCLUSION

Using a heat diffusion model in combination with GWAS results and protein-protein interactions, the significantly mutated subnetworks in rifampicin-resistant samples were found. The subnetwork not only contained the known targets of rifampicin (rpoB, rpoC), but also included topA and pyrG, which are potentially associated with drug resistance. Together, these results offer deeper insights into drug resistance of Mtb, and provides potential drug targets for finding new antituberculosis drugs.

Predictive Binding Affinity of Plant-Derived Natural Products Towards the Protein Kinase G Enzyme of Mycobacterium tuberculosis (MtPknG)

Tuberculosis (TB), caused by Mycobacterium tuberculosis, is a growing public health concern worldwide, especially with the emerging challenge of drug resistance to the current drugs. Efforts to discover and develop novel, more effective, and safer anti-TB drugs are urgently needed. Products from natural sources, such as medicinal plants, have played an important role in traditional medicine and continue to provide some inspiring templates for the design of new drugs. Protein kinase G, produced by M. tuberculosis (MtPKnG), is a serine/threonine kinase, that has been reported to prevent phagosome-lysosome fusion and help prolong M. tuberculosis survival within the host’s macrophages. Here, we used an in silico, target-based approach (docking) to predict the interactions between MtPknG and 84 chemical constituents from two medicinal plants (Pelargonium reniforme and Pelargonium sidoides) that have a well-documented historical use as natural remedies for TB. Docking scores for ligands towards the target protein were calculated using AutoDock Vina as the predicted binding free energies. Ten flavonoids present in the aerial parts of P. reniforme and/or P.sidoides showed docking scores ranging from −11.1 to −13.2 kcal/mol. Upon calculation of all ligand efficiency indices, we observed that the (−ΔG/MW) ligand efficiency index for flavonoids (4), (5) and (7) was similar to the one obtained for the AX20017 control. When taking all compounds into account, we observed that the best (−ΔG/MW) efficiency index was obtained for coumaric acid, coumaraldehyde, p-hydroxyphenyl acetic acid and p-hydroxybenzyl alcohol. We found that methyl gallate and myricetin had ligand efficiency indices superior and equal to the AX20017 control efficiency, respectively. It remains to be seen if any of the compounds screened in this study exert an effect in M. tuberculosis-infected macrophages.

Investigation of the anti-TB potential of selected propolis constituents using a molecular docking approach

Human tuberculosis (TB), caused by Mycobacterium tuberculosis, is the leading bacterial killer disease worldwide and new anti-TB drugs are urgently needed. Natural remedies have long played an important role in medicine and continue to provide some inspiring templates for drug design. Propolis, a substance naturally-produced by bees upon collection of plant resins, is used in folk medicine for its beneficial anti-TB activity. In this study, we used a molecular docking approach to investigate the interactions between selected propolis constituents and four ‘druggable’ proteins involved in vital physiological functions in M. tuberculosis, namely MtPanK, MtDprE1, MtPknB and MtKasA. The docking score for ligands towards each protein was calculated to estimate the binding free energy, with the best docking score (lowest energy value) indicating the highest predicted ligand/protein affinity. Specific interactions were also explored to understand the nature of intermolecular bonds between the most active ligands and the protein binding site residues. The lignan (+)-sesamin displayed the best docking score towards MtDprE1 (−10.7 kcal/mol) while the prenylated flavonoid isonymphaeol D docked strongly with MtKasA (−9.7 kcal/mol). Both compounds showed docking scores superior to the control inhibitors and represent potentially interesting scaffolds for further in vitro biological evaluation and anti-TB drug design.

Recent advancements in the development of anti-tuberculosis drugs

Modern chemotherapy has significantly improved patient outcomes against drug-sensitive tuberculosis. However, the rapid emergence of drug-resistant tuberculosis, together with the bacterium's ability to persist and remain latent present a major public health challenge. To overcome this problem, research into novel anti-tuberculosis targets and drug candidates is thus of paramount importance. This review article provides an overview of tuberculosis highlighting the recent advances and tools that are employed in the field of anti-tuberculosis drug discovery. The predominant focus is on anti-tuberculosis agents that are currently in the pipeline, i.e. clinical trials.

Bioinformatics Identification of Drug Resistance-Associated Gene Pairs in Mycobacterium tuberculosis

Tuberculosis is a chronic infectious disease caused by Mycobacterium tuberculosis (Mtb). Due to the extensive use of anti-tuberculosis drugs and the development of mutations, the emergence and spread of multidrug-resistant tuberculosis is recognized as one of the most dangerous threats to global tuberculosis control. Some single mutations have been identified to be significantly linked with drug resistance. However, the prior research did not take gene-gene interactions into account, and the emergence of transmissible drug resistance is connected with multiple genetic mutations. In this study we use the bioinformatics software GBOOST (The Hong Kong University, Clear Water Bay, Kowloon, Hong Kong, China) to calculate the interactions of Single Nucleotide Polymorphism (SNP) pairs and identify gene pairs associated with drug resistance. A large part of the non-synonymous mutations in the drug target genes that were included in the screened gene pairs were confirmed by previous reports, which lent sound solid credits to the effectiveness of our method. Notably, most of the identified gene pairs containing drug targets also comprise Pro-Pro-Glu (PPE) family proteins, suggesting that PPE family proteins play important roles in the drug resistance of Mtb. Therefore, this study provides deeper insights into the mechanisms underlying anti-tuberculosis drug resistance, and the present method is useful for exploring the drug resistance mechanisms for other microorganisms.

Fast assessment of structural models of ion channels based on their predicted current-voltage characteristics

Computational prediction of protein structures is a difficult task, which involves fast and accurate evaluation of candidate model structures. We propose to enhance single-model quality assessment with a functionality evaluation phase for proteins whose quantitative functional characteristics are known. In particular, this idea can be applied to evaluation of structural models of ion channels, whose main function - conducting ions - can be quantitatively measured with the patch-clamp technique providing the current-voltage characteristics. The study was performed on a set of KcsA channel models obtained from complete and incomplete contact maps. A fast continuous electrodiffusion model was used for calculating the current-voltage characteristics of structural models. We found that the computed charge selectivity and total current were sensitive to structural and electrostatic quality of models. In practical terms, we show that evaluating predicted conductance values is an appropriate method to eliminate models with an occluded pore or with multiple erroneously created pores. Moreover, filtering models on the basis of their predicted charge selectivity results in a substantial enrichment of the candidate set in highly accurate models. Tests on three other ion channels indicate that, in addition to being a proof of the concept, our function-oriented single-model quality assessment method can be directly applied to evaluation of structural models of some classes of protein channels. Finally, our work raises an important question whether a computational validation of functionality should be included in the evaluation process of structural models, whenever possible.

Structural Annotation of the Mycobacterium tuberculosis Proteome

Efforts from the TB Structural Genomics Consortium together with those of tuberculosis structural biologists worldwide have led to the determination of about 350 structures, making up nearly a tenth of the pathogen's proteome. Given that knowledge of protein structures is essential to obtaining a high-resolution understanding of the underlying biology, it is desirable to have a structural view of the entire proteome. Indeed, structure prediction methods have advanced sufficiently to allow structural models of many more proteins to be built based on homology modeling and fold recognition strategies. By means of these approaches, structural models for about 2,877 proteins, making up nearly 70% of the Mycobacterium tuberculosis proteome, are available. Knowledge from bioinformatics has made significant inroads into an improved annotation of the M. tuberculosis genome and in the prediction of key protein players that interact in vital pathways, some of which are unique to the organism. Functional inferences have been made for a large number of proteins based on fold-function associations. More importantly, ligand-binding pockets of the proteins are identified and scanned against a large database, leading to binding site-based ligand associations and hence structure-based function annotation. Near proteome-wide structural models provide a global perspective of the fold distribution in the genome. New insights about the folds that predominate in the genome, as well as the fold combinations that make up multidomain proteins, are also obtained. This chapter describes the structural proteome, functional inferences drawn from it, and its applications in drug discovery.

Rv2969c, essential for optimal growth in Mycobacterium tuberculosis, is a DsbA-like enzyme that interacts with VKOR-derived peptides and has atypical features of DsbA-like disulfide oxidases

The bacterial disulfide machinery is an attractive molecular target for developing new antibacterials because it is required for the production of multiple virulence factors. The archetypal disulfide oxidase proteins in Escherichia coli (Ec) are DsbA and DsbB, which together form a functional unit: DsbA introduces disulfides into folding proteins and DsbB reoxidizes DsbA to maintain it in the active form. In Mycobacterium tuberculosis (Mtb), no DsbB homologue is encoded but a functionally similar but structurally divergent protein, MtbVKOR, has been identified. Here, the Mtb protein Rv2969c is investigated and it is shown that it is the DsbA-like partner protein of MtbVKOR. It is found that it has the characteristic redox features of a DsbA-like protein: a highly acidic catalytic cysteine, a highly oxidizing potential and a destabilizing active-site disulfide bond. Rv2969c also has peptide-oxidizing activity and recognizes peptide segments derived from the periplasmic loops of MtbVKOR. Unlike the archetypal EcDsbA enzyme, Rv2969c has little or no activity in disulfide-reducing and disulfide-isomerase assays. The crystal structure of Rv2969c reveals a canonical DsbA fold comprising a thioredoxin domain with an embedded helical domain. However, Rv2969c diverges considerably from other DsbAs, including having an additional C-terminal helix (H8) that may restrain the mobility of the catalytic helix H1. The enzyme is also characterized by a very shallow hydrophobic binding surface and a negative electrostatic surface potential surrounding the catalytic cysteine. The structure of Rv2969c was also used to model the structure of a paralogous DsbA-like domain of the Ser/Thr protein kinase PknE. Together, these results show that Rv2969c is a DsbA-like protein with unique properties and a limited substrate-binding specificity.

A novel clade of unique eukaryotic ribonucleotide reductase R2 subunits is exclusive to apicomplexan parasites

Apicomplexa are protist parasites of tremendous medical and economic importance, causing millions of deaths and billions of dollars in losses each year. Apicomplexan-related diseases may be controlled via inhibition of essential enzymes. Ribonucleotide reductase (RNR) provides the only de novo means of synthesizing deoxyribonucleotides, essential precursors for DNA replication and repair. RNR has long been the target of antibacterial and antiviral therapeutics. However, targeting this ubiquitous protein in eukaryotic pathogens may be problematic unless these proteins differ significantly from that of their respective host. The typical eukaryotic RNR enzymes belong to class Ia, and the holoenzyme consists minimally of two R1 and two R2 subunits (α₂β₂). We generated a comparative, annotated, structure-based, multiple-sequence alignment of R2 subunits, identified a clade of R2 subunits unique to Apicomplexa, and determined its phylogenetic position. Our analyses revealed that the apicomplexan-specific sequences share characteristics with both class I R2 and R2lox proteins. The putative radical-harboring residue, essential for the reduction reaction by class Ia R2-containing holoenzymes, was not conserved within this group. Phylogenetic analyses suggest that class Ia subunits are not monophyletic and consistently placed the apicomplexan-specific clade sister to the remaining class Ia eukaryote R2 subunits. Our research suggests that the novel apicomplexan R2 subunit may be a promising candidate for chemotherapeutic-induced inhibition as it differs greatly from known eukaryotic host RNRs and may be specifically targeted.