Effects of the magnesium and chloride ions and shikimate on the structure of shikimate kinase from Mycobacterium tuberculosis

An overview of mechanism and chemical inhibitors of shikimate kinase

Tuberculosis is a highly infectious disease other than HIV/AIDS and it is one of the top ten causes of death worldwide. Resistance development in the bacteria occurs because of genetic alterations, and the molecular insights suggest that the accumulation of mutation in the individual drug target genes is the primary mechanism of multi-drug resistant tuberculosis. Chorismate is an essential structural fragment for the synthesis of aromatic amino acids and synthesized biochemically by a number of bacteria, including Mycobacterium tuberculosis, utilizing the shikimate pathway. This shikimate kinase is the newer possible target for the generation of novel antitubercular drug because this pathway is expressed only in mycobacterium and not in Mammals. The discovery and development of shikimate kinase inhibitors provide an opportunity for the development of novel selective medications. Multiple shikimate kinase inhibitors have been identified via insilico virtual screening and related protein-ligand interactions along with their in-vitro studies. These inhibitors bind to the active site in a similar fashion to shikimate. In the current review, we present an overview of the biology and chemistry of the shikimate kinase protein and its inhibitors, with special emphasis on the various active scaffold against the enzyme. A variety of chemically diversified synthetic scaffolds including Benzothiazoles, Oxadiazoles, Thiobarbiturates, Naphthoquinones, Thiazoleacetonitriles, Hybridized Pyrazolone derivatives, Orthologous biological macromolecule derivatives, Manzamine Alkaloids derivatives, Dipeptide inhibitor, and Chalcones are discussed in detail. These derivatives bind to the specific target appropriately proving their potential ability through different binding interactions and effectively explored as an effective and selective Sk inhibitor.Communicated by Ramaswamy H. Sarma.

Shikimate Kinase Inhibitors: An Update on Promising Strategy against Mycobacterium tuberculosis

Humanity has been battling with tuberculosis (TB) for a long period, and despite the availability of drugs well-known to act against the deadly microbe, the menace is still very far from reaching its end. Moreover, problems related to TB chemotherapy, such as lengthy treatment periods leading to poor patient compliance, increasing drug resistance, and association with another deadlier disease HIV-AIDS, make the situation alarming, thereby pressing the need for the discovery of new potent drugs urgently. Therefore, a drug target that is essential for survival and exclusive to M. tuberculosis presents a promising platform to explore novel molecules against the microorganism for better pathogen clearance with minimal toxicity. The shikimate pathway that leads to the synthesis of essential aromatic amino acids is one such attractive target. Shikimate kinase, the fifth enzyme of this pathway, converts shikimate to shikimate-3-phosphate by using ATP as a cosubstrate. Targeting shikimate kinase could be an effective strategy in light of its essentiality and absence of any homologue in mammals. This review discusses different strategies adopted for discovering novel compounds or scaffolds targeting M. tuberculosis shikimate kinase (MtSK) in vitro. The application of substrate analogues, their structure, and ligand-based approach for screening a library of anti-mycobacterial compounds, marine-derived molecules, and commercially available libraries have yielded promising MtSK inhibitors exhibiting micro-molar activities. To develop these leads into future drugs with minimum off-target effects on the host microenvironment, the molecules need to be structurally optimized for improved activities against enzymes and whole-cell organisms.

Identification of novel scaffold using ligand and structure based approach targeting shikimate kinase

Tuberculosis (TB) remains a major global health problem. It causes ill-health among millions of people each year and rank as the second leading cause of death from an infectious disease worldwide, after the human immunodeficiency virus (HIV). Shikimate kinase is one of the major enzymes targeted for TB. Most approaches to overcome TB were based on synthesis and screening of a known compounds to obtain a few representatives with desired potency. In this study, we have applied a virtual screening approach which combines ligand- and structure-based approaches to screen a large library of compounds as a starting point for the identification of new scaffolds for the development of shikimate kinase inhibitors. The combined approach has identified 2 new scaffolds as potential inhibitors of shikimate kinase. To prove the approach, few of the molecules and their derivatives, a total of 17 compounds, were synthesized. The compounds were tested for biological activity and shows moderate activity against shikimate kinase. The shikimate kinase enzyme inhibition study reveals that the compounds showed inhibition (IC₅₀) at concentrations of 50 µg/mL (Compounds 21, 22, 24, 25, 26, 27, 30, 32, 34) and 25 µg/mL (14, 19, 23, 31, 33).

How Docking Programs Work

Protein-ligand docking simulations are of central interest for computer-aided drug design. Docking is also of pivotal importance to understand the structural basis for protein-ligand binding affinity. In the last decades, we have seen an explosion in the number of three-dimensional structures of protein-ligand complexes available at the Protein Data Bank. These structures gave further support for the development and validation of in silico approaches to address the binding of small molecules to proteins. As a result, we have now dozens of open source programs and web servers to carry out molecular docking simulations. The development of the docking programs and the success of such simulations called the attention of a broad spectrum of researchers not necessarily familiar with computer simulations. In this scenario, it is essential for those involved in experimental studies of protein-ligand interactions and biophysical techniques to have a glimpse of the basics of the protein-ligand docking simulations. Applications of protein-ligand docking simulations to drug development and discovery were able to identify hits, inhibitors, and even drugs. In the present chapter, we cover the fundamental ideas behind protein-ligand docking programs for non-specialists, which may benefit from such knowledge when studying molecular recognition mechanism.

Docking with GemDock

GEMDOCK is a protein-ligand docking software that makes use of an elegant biologically inspired computational methodology based on the differential evolution algorithm. As any docking program, GEMDOCK has two major features to predict the binding of a small-molecule ligand to the binding site of a protein target: the search algorithm and the scoring function to evaluate the generated poses. The GEMDOCK scoring function uses a piecewise potential energy function integrated into the differential evolutionary algorithm. GEMDOCK has been applied to a wide range of protein systems with docking accuracy similar to other docking programs such as Molegro Virtual Docker, AutoDock4, and AutoDock Vina. In this chapter, we explain how to carry out protein-ligand docking simulations with GEMDOCK. We focus this tutorial on the protein target cyclin-dependent kinase 2.

Electrostatic Energy in Protein-Ligand Complexes

Computational analysis of protein-ligand interactions is of pivotal importance for drug design. Assessment of ligand binding energy allows us to have a glimpse of the potential of a small organic molecule as a ligand to the binding site of a protein target. Considering scoring functions available in docking programs such as AutoDock4, AutoDock Vina, and Molegro Virtual Docker, we could say that they all rely on equations that sum each type of protein-ligand interactions to model the binding affinity. Most of the scoring functions consider electrostatic interactions involving the protein and the ligand. In this chapter, we present the main physics concepts necessary to understand electrostatics interactions relevant to molecular recognition of a ligand by the binding pocket of a protein target. Moreover, we analyze the electrostatic potential energy for an ensemble of structures to highlight the main features related to the importance of this interaction for binding affinity.

SAnDReS: A Computational Tool for Docking

Since the early 1980s, we have witnessed considerable progress in the development and application of docking programs to assess protein-ligand interactions. Most of these applications had as a goal the identification of potential new binders to protein targets. Another remarkable progress is taking place in the determination of the structures of protein-ligand complexes, mostly using X-ray diffraction crystallography. Considering these developments, we have a favorable scenario for the creation of a computational tool that integrates into one workflow all steps involved in molecular docking simulations. We had these goals in mind when we developed the program SAnDReS. This program allows the integration of all computational features related to modern docking studies into one workflow. SAnDReS not only carries out docking simulations but also evaluates several docking protocols allowing the selection of the best approach for a given protein system. SAnDReS is a free and open-source (GNU General Public License) computational environment for running docking simulations. Here, we describe the combination of SAnDReS and AutoDock4 for protein-ligand docking simulations. AutoDock4 is a free program that has been applied to over a thousand receptor-ligand docking simulations. The dataset described in this chapter is available for downloading at https://github.com/azevedolab/sandres.

Van der Waals Potential in Protein Complexes

Van der Waals forces are determinants of the formation of protein-ligand complexes. Physical models based on the Lennard-Jones potential can estimate van der Waals interactions with considerable accuracy and with a computational complexity that allows its application to molecular docking simulations and virtual screening of large databases of small organic molecules. Several empirical scoring functions used to evaluate protein-ligand interactions approximate van der Waals interactions with the Lennard-Jones potential. In this chapter, we present the main concepts necessary to understand van der Waals interactions relevant to molecular recognition of a ligand by the binding pocket of a protein target. We describe the Lennard-Jones potential and its application to calculate potential energy for an ensemble of structures to highlight the main features related to the importance of this interaction for binding affinity.

Homology Modeling of Protein Targets with MODELLER

Homology modeling is a computational approach to generate three-dimensional structures of protein targets when experimental data about similar proteins are available. Although experimental methods such as X-ray crystallography and nuclear magnetic resonance spectroscopy successfully solved the structures of nearly 150,000 macromolecules, there is still a gap in our structural knowledge. We can fulfill this gap with computational methodologies. Our goal in this chapter is to explain how to perform homology modeling of protein targets for drug development. We choose as a homology modeling tool the program MODELLER. To illustrate its use, we describe how to model the structure of human cyclin-dependent kinase 3 using MODELLER. We explain the modeling procedure of CDK3 apoenzyme and the structure of this enzyme in complex with roscovitine.

Hydrogen Bonds in Protein-Ligand Complexes

Fast and reliable evaluation of the hydrogen bond potential energy has a significant impact in the drug design and development since it allows the assessment of large databases of organic molecules in virtual screening projects focused on a protein of interest. Semi-empirical force fields implemented in molecular docking programs make it possible the evaluation of protein-ligand binding affinity where the hydrogen bond potential is a common term used in the calculation. In this chapter, we describe the concepts behind the programs used to predict hydrogen bond potential energy employing semi-empirical force fields as the ones available in the programs AMBER, AutoDock4, TreeDock, and ReplicOpter. We described here the 12-10 potential and applied it to evaluate the binding affinity for an ensemble of crystallographic structures for which experimental data about binding affinity are available.

Slow-Binding Inhibition of Mycobacterium tuberculosis Shikimate Kinase by Manzamine Alkaloids

Tuberculosis represents a significant public health crisis. There is an urgent need for novel molecular scaffolds against this pathogen. We screened a small library of marine-derived compounds against shikimate kinase from Mycobacterium tuberculosis ( MtSK), a promising target for antitubercular drug development. Six manzamines previously shown to be active against M. tuberculosis were characterized as MtSK inhibitors: manzamine A (1), 8-hydroxymanzamine A (2), manzamine E (3), manzamine F (4), 6-deoxymanzamine X (5), and 6-cyclohexamidomanzamine A (6). All six showed mixed noncompetitive inhibition of MtSK. The lowest K_I values were obtained for 6 across all MtSK-substrate complexes. Time-dependent analyses revealed two-step, slow-binding inhibition. The behavior of 1 was typical; initial formation of an enzyme-inhibitor complex (EI) obeyed an apparent K_I of ∼30 μM with forward ( k₅) and reverse ( k₆) rate constants for isomerization to an EI* complex of 0.18 and 0.08 min^-1, respectively. In contrast, 6 showed a lower K_I for the initial encounter complex (∼1.5 μM), substantially faster isomerization to EI* ( k₅ = 0.91 min^-1), and slower back conversion of EI* to EI ( k₆ = 0.04 min^-1). Thus, the overall inhibition constants, K_I*, for 1 and 6 were 10 and 0.06 μM, respectively. These findings were consistent with docking predictions of a favorable binding mode and a second, less tightly bound pose for 6 at MtSK. Our results suggest that manzamines, in particular 6, constitute a new scaffold from which drug candidates with novel mechanisms of action could be designed for the treatment of tuberculosis by targeting MtSK.

Computational Approaches in Multitarget Drug Discovery

Current therapeutic strategies entail identifying and characterizing a single protein receptor whose inhibition is likely to result in the successful treatment of a disease of interest, and testing experimentally large libraries of small molecule compounds "in vitro" and "in vivo" to identify promising inhibitors in model systems and determine if the findings are extensible to humans. This highly complex process is largely based on tests, errors, risk, time, and intensive costs. The virtual computational study of compounds simulates situations predicting possible drug linkages with multiple protein target atomic structures, taking into account the dynamic protein inhibitor, and can help identify inhibitors efficiently, particularly for complex drug-resistant diseases. Some discussions will relate to the potential benefits of this approach, using HIV-1 and Plasmodium falciparum infections as examples. Some authors have proposed a virtual drug discovery that not only identifies efficient inhibitors but also helps to minimize side effects and toxicity, thus increasing the likelihood of successful therapies. This chapter discusses concepts and research of bioactive multitargets related to toxicology.

Benzothiazole Derivative as a Novel Mycobacterium tuberculosis Shikimate Kinase Inhibitor: Identification and Elucidation of Its Allosteric Mode of Inhibition

Mycobacterium tuberculosis shikimate kinase (Mtb-SK) is a key enzyme involved in the biosynthesis of aromatic amino acids through the shikimate pathway. Since it is proven to be essential for the survival of the microbe and is absent from mammals, it is a promising target for anti-TB drug discovery. In this study, a combined approach of in silico similarity search and pharmacophore building using already reported inhibitors was used to screen a procured library of 20,000 compounds of the commercially available ChemBridge database. From the in silico screening, 15 hits were identified, and these hits were evaluated in vitro for Mtb-SK enzyme inhibition. Two compounds presented significant enzyme inhibition with IC50 values of 10.69 ± 0.9 and 46.22 ± 1.2 μM. The best hit was then evaluated for the in vitro mode of inhibition where it came out to be an uncompetitive and noncompetitive inhibitor with respect to shikimate (SKM) and ATP, respectively, suggesting its binding at an allosteric site. Potential binding sites of Mtb-SK were identified which confirmed the presence of an allosteric binding pocket apart from the ATP and SKM binding sites. The docking simulations were performed at this pocket in order to find the mode of binding of the best hit in the presence of substrates and the products of the enzymatic reaction. Molecular dynamics (MD) simulations elucidated the probability of inhibitor binding at the allosteric site in the presence of ADP and shikimate-3-phosphate (S-3-P), that is, after the formation of products of the reaction. The inhibitor binding may prevent the release of the product from Mtb-SK, thereby inhibiting its activity. The binding stability and the key residue interactions of the inhibitor to this product complex were also revealed by the MD simulations. Residues ARG43, ILE45, and PHE57 were identified as crucial that were involved in interactions with the best hit. This is the first report of an allosteric binding site of Mtb-SK, which could largely address the selectivity issue associated with kinase inhibitors.

Structure of shikimate kinase, an in vivo essential metabolic enzyme in the nosocomial pathogen Acinetobacter baumannii, in complex with shikimate

Acinetobacter baumannii is an opportunistic Gram-negative pathogen that is an important cause of healthcare-associated infections exhibiting high mortality rates. Clinical isolates of multidrug-resistant (MDR) and extremely drug-resistant (XDR) A. baumannii strains are increasingly being observed. Compounding this concern is the dearth of new antibacterial agents in late-stage development that are effective against MDR and XDR A. baumannii. As part of an effort to address these concerns, two genes (aroA and aroC) of the shikimate pathway have previously been determined to be essential for the growth and survival of A. baumannii during host infection (i.e. to be essential in vivo). This study expands upon these results by demonstrating that the A. baumannii aroK gene, encoding shikimate kinase (SK), is also essential in vivo in a rat soft-tissue infection model. The crystal structure of A. baumannii SK in complex with the substrate shikimate and a sulfate ion that mimics the binding interactions expected for the β-phosphate of ATP was then determined to 1.91 Å resolution and the enzyme kinetics were characterized. The flexible shikimate-binding domain and LID region are compared with the analogous regions in other SK crystal structures. The impact of structural differences and sequence divergence between SKs from pathogenic bacteria that may influence antibiotic-development efforts is discussed.

"Superchiral" Spectroscopy: Detection of Protein Higher Order Hierarchical Structure with Chiral Plasmonic Nanostructures

Optical spectroscopic methods do not routinely provide information on higher order hierarchical structure (tertiary/quaternary) of biological macromolecules and assemblies. This necessitates the use of time-consuming and material intensive techniques, such as protein crystallography, NMR, and electron microscopy. Here we demonstrate a spectroscopic phenomenon, superchiral polarimetry, which can rapidly characterize ligand-induced changes in protein higher order (tertiary/quaternary) structure at the picogram level, which is undetectable using conventional CD spectroscopy. This is achieved by utilizing the enhanced sensitivity of superchiral evanescent fields to mesoscale chiral structure.

Selective Mycobacterium tuberculosis Shikimate Kinase Inhibitors as Potential Antibacterials

Owing to the persistence of tuberculosis (TB) as well as the emergence of multidrug-resistant and extensively drug-resistant (XDR) forms of the disease, the development of new antitubercular drugs is crucial. Developing inhibitors of shikimate kinase (SK) in the shikimate pathway will provide a selective target for antitubercular agents. Many studies have used in silico technology to identify compounds that are anticipated to interact with and inhibit SK. To a much more limited extent, SK inhibition has been evaluated by in vitro methods with purified enzyme. Currently, there are no data on in vivo activity of Mycobacterium tuberculosis shikimate kinase (MtSK) inhibitors available in the literature. In this review, we present a summary of the progress of SK inhibitor discovery and evaluation with particular attention toward development of new antitubercular agents.

The mode of action of recombinant Mycobacterium tuberculosis shikimate kinase: kinetics and thermodynamics analyses

Tuberculosis remains as one of the main cause of mortality worldwide due to a single infectious agent, Mycobacterium tuberculosis. The aroK-encoded M. tuberculosis Shikimate Kinase (MtSK), shown to be essential for survival of bacilli, catalyzes the phosphoryl transfer from ATP to the carbon-3 hydroxyl group of shikimate (SKH), yielding shikimate-3-phosphate and ADP. Here we present purification to homogeneity, and oligomeric state determination of recombinant MtSK. Biochemical and biophysical data suggest that the chemical reaction catalyzed by monomeric MtSK follows a rapid-equilibrium random order of substrate binding, and ordered product release. Isothermal titration calorimetry (ITC) for binding of ligands to MtSK provided thermodynamic signatures of non-covalent interactions to each process. A comparison of steady-state kinetics parameters and equilibrium dissociation constant value determined by ITC showed that ATP binding does not increase the affinity of MtSK for SKH. We suggest that MtSK would more appropriately be described as an aroL-encoded type II shikimate kinase. Our manuscript also gives thermodynamic description of SKH binding to MtSK and data for the number of protons exchanged during this bimolecular interaction. The negative value for the change in constant pressure heat capacity (ΔCp) and molecular homology model building suggest a pronounced contribution of desolvation of non-polar groups upon binary complex formation. Thermodynamic parameters were deconvoluted into hydrophobic and vibrational contributions upon MtSK:SKH binary complex formation. Data for the number of protons exchanged during this bimolecular interaction are interpreted in light of a structural model to try to propose the likely amino acid side chains that are the proton donors to bulk solvent following MtSK:SKH complex formation.

New tuberculosis drug development: targeting the shikimate pathway

BACKGROUND

Tuberculosis (TB) remains a leading cause of morbidity and mortality worldwide, yet no new drugs have been developed in the last 40 years.

OBJECTIVE

The exceedingly lengthy TB chemotherapy and the increasing emergence of drug resistance complicated by HIV co-infection call for the development of new TB drugs. These problems are further compounded by a poor understanding of the biology of persister bacteria.

METHODS

New molecular tools have offered insights into potential new drug targets, particularly the enzymes of the shikimate pathway, which is the focus of this review.

RESULTS/CONCLUSION

Shikimate pathway enzymes, especially shikimate kinase, may offer attractive targets for new TB drug and vaccine development.

Comprehensive structural classification of ligand-binding motifs in proteins

Comprehensive knowledge of protein-ligand interactions should provide a useful basis for annotating protein functions, studying protein evolution, engineering enzymatic activity, and designing drugs. To investigate the diversity and universality of ligand-binding sites in protein structures, we conducted the all-against-all atomic-level structural comparison of over 180,000 ligand-binding sites found in all the known structures in the Protein Data Bank by using a recently developed database search and alignment algorithm. By applying a hybrid top-down-bottom-up clustering analysis to the comparison results, we determined approximately 3000 well-defined structural motifs of ligand-binding sites. Apart from a handful of exceptions, most structural motifs were found to be confined within single families or superfamilies, and to be associated with particular ligands. Furthermore, we analyzed the components of the similarity network and enumerated more than 4000 pairs of structural motifs that were shared across different protein folds.

Structural studies of shikimate dehydrogenase from Bacillus anthracis complexed with cofactor NADP

Bacillus anthracis has been employed as an agent of bioterrorism, with high mortality, despite anti-microbial treatment, which strongly indicates the need of new drugs to treat anthrax. Shikimate pathway is a seven step biosynthetic route which generates chorismic acid from phosphoenol pyruvate and erythrose-4-phosphate. Chorismic acid is the major branch point in the synthesis of aromatic amino acids, ubiquinone, and secondary metabolites. The shikimate pathway is essential for many pathological organisms, whereas it is absent in mammals. Therefore, these enzymes are potential targets for the development of nontoxic antimicrobial agents and herbicides and have been submitted to intensive structural studies. The forth enzyme of this pathway is responsible for the conversion of dehydroshikimate to shikimate in the presence of NADP. In order to pave the way for structural and functional efforts toward development of new antimicrobials we describe the molecular modeling of shikimate dehydrogenase from Bacillus anthracis complexed with the cofactor NADP. This study was able to identify the main residues of the NADP binding site responsible for ligand affinities. This structural study can be used in the design of more specific drugs against infectious diseases.

Molecular modeling and dynamics studies of Shikimate Kinase from Bacillus anthracis

Bacillus anthracis has been used as weapon in bioterrorist activities, with high mortality, despite anti-microbial treatment, which strongly indicates a need of new drugs to treat anthrax. Shikimate Pathway is a seven-step biosynthetic route which generates chorismic acid. The shikimate pathway is essential for many pathological organisms, whereas it is absent in mammals. Therefore, these enzymes are potential targets for the development of non-toxic anti-microbial agents and herbicides and have been submitted to intensive structural studies. Shikimate Kinase is the fifth enzyme of shikimate pathway and catalyzes the specific phosphorylation of the 3-hydroxyl group of shikimate using ATP as a co-substrate, resulting in shikimate-3-phosphate and ADP. The present work describes for the first time a structural model for the Shikimate Kinase from B. anthracis using molecular modeling approach and molecular dynamics simulations. This study was able to identify the main residues of the ATP-binding and the shikimate pockets responsible for ligand affinities. Analysis of the molecular dynamics simulations indicates the structural features responsible for the stability of the structure. This study may help in the identification of new inhibitors for this enzyme.