Molecular modeling and dynamics studies of cytidylate kinase from Mycobacterium tuberculosis H37Rv

Resistance related metabolic pathways for drug target identification in Mycobacterium tuberculosis

BACKGROUND

Increasing resistance to anti-tuberculosis drugs has driven the need for developing new drugs. Resources such as the tropical disease research (TDR) target database and AssessDrugTarget can help to prioritize putative drug targets. Hower, these resources do not necessarily map to metabolic pathways and the targets are not involved in dormancy. In this study, we specifically identify drug resistance pathways to allow known drug resistant mutations in one target to be offset by inhibiting another enzyme of the same metabolic pathway. One of the putative targets, Rv1712, was analysed by modelling its three dimensional structure and docking potential inhibitors.

RESULTS

We mapped 18 TB drug resistance gene products to 15 metabolic pathways critical for mycobacterial growth and latent TB by screening publicly available microarray data. Nine putative targets, Rv1712, Rv2984, Rv2194, Rv1311, Rv1305, Rv2195, Rv1622c, Rv1456c and Rv2421c, were found to be essential, to lack a close human homolog, and to share >67 % sequence identity and >87 % query coverage with mycobacterial orthologs. A structural model was generated for Rv1712, subjected to molecular dynamic simulation, and identified 10 compounds with affinities better than that for the ligand cytidine-5'-monophosphate (C5P). Each compound formed more interactions with the protein than C5P.

CONCLUSIONS

We focused on metabolic pathways associated with bacterial drug resistance and proteins unique to pathogenic bacteria to identify novel putative drug targets. The ten compounds identified in this study should be considered for experimental studies to validate their potential as inhibitors of Rv1712.

Nucleotide Metabolism and DNA Replication

The development and application of a highly versatile suite of tools for mycobacterial genetics, coupled with widespread use of "omics" approaches to elucidate the structure, function, and regulation of mycobacterial proteins, has led to spectacular advances in our understanding of the metabolism and physiology of mycobacteria. In this article, we provide an update on nucleotide metabolism and DNA replication in mycobacteria, highlighting key findings from the past 10 to 15 years. In the first section, we focus on nucleotide metabolism, ranging from the biosynthesis, salvage, and interconversion of purine and pyrimidine ribonucleotides to the formation of deoxyribonucleotides. The second part of the article is devoted to DNA replication, with a focus on replication initiation and elongation, as well as DNA unwinding. We provide an overview of replication fidelity and mutation rates in mycobacteria and summarize evidence suggesting that DNA replication occurs during states of low metabolic activity, and conclude by suggesting directions for future research to address key outstanding questions. Although this article focuses primarily on observations from Mycobacterium tuberculosis, it is interspersed, where appropriate, with insights from, and comparisons with, other mycobacterial species as well as better characterized bacterial models such as Escherichia coli. Finally, a common theme underlying almost all studies of mycobacterial metabolism is the potential to identify and validate functions or pathways that can be exploited for tuberculosis drug discovery. In this context, we have specifically highlighted those processes in mycobacterial DNA replication that might satisfy this critical requirement.

Kinetic mechanism and energetics of binding of phosphoryl group acceptors to Mycobacterium tuberculosis cytidine monophosphate kinase

Cytidine monophosphate kinase from Mycobacterium tuberculosis (MtCMK) likely plays a role in supplying precursors for nucleic acid synthesis. MtCMK catalyzes the ATP-dependent phosphoryl group transfer preferentially to CMP and dCMP. Initial velocity studies and Isothermal titration calorimetry (ITC) measurements showed that MtCMK follows a random-order mechanism of substrate (CMP and ATP) binding, and an ordered mechanism for product release, in which ADP is released first followed by CDP. The thermodynamic signatures of CMP and CDP binding to MtCMK showed favorable enthalpy and unfavorable entropy, and ATP binding was characterized by favorable changes in enthalpy and entropy. The contribution of linked protonation events to the energetics of MtCMK:phosphoryl group acceptor binary complex formation suggested a net gain of protons. Values for the pKa of a likely chemical group involved in proton exchange and for the intrinsic binding enthalpy were calculated. The Asp187 side chain of MtCMK is suggested as the likely candidate for the protonation event. Data on thermodynamics of binary complex formation were collected to evaluate the contribution of 2'-OH group to intermolecular interactions. The data are discussed in light of functional and structural comparisons between CMP/dCMP kinases and UMP/CMP ones.

Theoretical study of structure, bonding, and electronic behavior of novel sandwich compounds M₃(C6R6)₂ (M = Ni, Pd, Pt; R = H, F)

The correlations between the structural and electronic properties of the monolayer clusters M₃ (where M = Ni, Pd, Pt) and the sandwich complexes M₃(C₆R₆)₂ (where M = Ni, Pd, Pt; R = H, F) were studied by performing quantum-chemical calculations. All of the sandwich complexes are strongly donating and backdonating metal-ligand bonding structures. The influence of the ligand as well as significant variations in the M-C, M-M, and C-C bond lengths and binding energies were examined to obtain a qualitative and quantitative picture of the intramolecular interactions in C₆R₆-M₃. Our theoretical investigations show that the binding energies of these sandwich complexes gradually decrease from Ni to Pt as well as from H to F, which can be explained via the frontier orbitals of the clusters M₃ and C₆R₆.

Binding dynamics and energetic insight into the molecular forces driving nucleotide binding by guanylate kinase

Plasmodium deoxyguanylate pathways are an attractive area of investigation for future metabolic and drug discovery studies due to their unique substrate specificities. We investigated the energetic contribution to guanylate kinase substrate binding and the forces underlying ligand recognition. In the range from 20 to 35°C, the thermodynamic profiles displayed marked decrease in binding enthalpy, while the free energy of binding showed little changes. GMP produced a large binding heat capacity change of -356 cal mol(-1) K(-1), indicating considerable conformational changes upon ligand binding. Interestingly, the calculated ΔCp was -32 cal mol(-1) K(-1), indicating that the accessible surface area is not the central change in substrate binding, and that other entropic forces, including conformational changes, are more predominant. The thermodynamic signature for GMP is inconsistent with rigid-body association, while dGMP showed more or less rigid-body association. These binding profiles explain the poor catalytic efficiency and low affinity for dGMP compared with GMP. At low temperature, the ligands bind to the receptor site under the effect of hydrophobic forces. Interestingly, by increasing the temperature, the entropic forces gradually vanish and proceed to a nonfavorable contribution, and the interaction occurs mainly through bonding, electrostatic forces, and van der Waals interactions.

Exploring the molecular basis for selective binding of Mycobacterium tuberculosis Asp kinase toward its natural substrates and feedback inhibitors: a docking and molecular dynamics study

Tuberculosis (TB) is still a major public health problem, compounded by the human immunodeficiency virus (HIV)-TB co-infection and recent emergence of multidrug-resistant (MDR) and extensively drug resistant (XDR)-TB. In this context, aspartokinase of mycobacterium tuberculosis has drawn attention for designing novel anti-TB drugs. Asp kinase is an enzyme responsible for the synthesis of 4-phospho-L-aspartate from L-aspartate and involved in the branched biosynthetic pathway leading to the synthesis of amino acids lysine, threonine, methionine and isoleucine. An intermediate of lysine biosynthetic branch, mesodiaminopimelate is also a component of the peptidoglycan which is a component of bacterial cell wall. To interfere with the production of all these amino acids and cell wall, it is possible to inhibit Asp kinase activity. This can be achieved using Asp kinase inhibitors. In order to design novel Asp kinase inhibitors as effective anti-TB drugs, it is necessary to have an understanding of the binding sites of Asp kinase. As no crystal structure of the enzyme has yet been published, we built a homology model of Asp kinase using the crystallized Asp kinase from M. Jannaschii, as template structures (2HMF and 3C1M). After the molecular dynamics refinement, the optimized homology model was assessed as a reliable structure by PROCHECK, ERRAT, WHAT-IF, PROSA2003 and VERIFY-3D. The results of molecular docking studies with natural substrates, products and feedback inhibitors are in agreement with the published data and showed that ACT domain plays an important role in binding to ligands. Based on the docking conformations, pharmacophore model can be developed by probing the common features of ligands. By analyzing the results, ACT domain architecture, certain key residues that are responsible for binding to feedback inhibitors and natural substrates were identified. This would be very helpful in understanding the blockade mechanism of Asp kinase and providing insights into rational design of novel Asp kinase inhibitors for M.tuberculosis.

Molecular modeling and dynamics simulation of human cyclin-dependent kinase 3 complexed with inhibitors

The complex CDK3-cyclin is involved in the control of the progression of G0. While the mechanisms governing early and late G1 progression are well understood, very little is known about the G0-G1 transition. Human CDK3 is closely related to cyclin-dependent kinase 2 (CDK2). Since there is no crystallographic structure of human CDK3, this work describes for the first time a molecular model of human CDK3 complexed with several inhibitors. Comparison of the binary complexes with different inhibitors strongly indicates that those inhibitors should inhibit CDK3 as well as CDK2.

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.