- Open Access
Mechanistic analysis elucidating the relationship between Lys96 mutation in Mycobacterium tuberculosispyrazinamidase enzyme and pyrazinamide susceptibility
© Vats et al.; licensee BioMed Central Ltd. 2015
- Published: 21 January 2015
- drug resistance
- molecular dynamics simulations
Pyrazinamide (PZA) is one of the most effective first line treatments against tuberculosis disease. The drug generally has bacteriostatic action. It also acts on bacterial spores which eliminates the chances of resurfacing of the infection. However, in recent years there has been a major increase in the occurrence of drug resistant bacterial strains. Resistance against PZA is caused by mutations in pyrazinamidase (PncA) protein which is the activator of the prodrug PZA. In the present study, we have tried to gain insights into the mechanism by which resistance develops due to K96R mutation occurring in the PncA catalytic region
The binding cavity analysis showed an increase of 762.3 Å3 in the volume of the mutant protein. Docking studies revealed that PZA has a greater binding affinity for the native protein in comparison to the mutant protein. Molecular dynamics simulations further showed the role of flap region, which is present in PncA protein, in development of resistance to the drug.
The residues of flap region acquire more flexibility in mutant form of protein and thus move away from the active site. This leads to weak binding of the drug to the target residues which might interfere with the activation of the drug to functional form thereby giving rise to drug resistant bacterial strains.
One of the major concerns being faced by the world today is tuberculosis disease (TB). About one third of the world population is suffering from the infection caused by Mycobacterium tuberculosis . It was estimated, that in 2012, around 8.6 million people developed TB and approximately 1.3 million people died because of this disease . The emergence of drug-resistant bacterial strains is one of the main causes for the current spread of TB . MDR TB or multi drug resistant tuberculosis is caused when resistance is developed for two or more first line drugs. Treatment of drug-resistant tuberculosis is hampered by poor efficacy and high toxicity of the second-line drugs . Continuous efforts are being made worldwide to find out some novel protein targets against which new effective drugs can be designed [5, 6].
Various new approaches have come up for the designing of drugs for such complex diseases. One of the strategies includes the development of small molecule which could either inhibit multiple sites on the same target or inhibit altogether different protein targets. Some previous work done in our lab to check pathologies associated with the Alzheimer's disease illustrates this well. In one of the experiments two small molecular compounds were proposed, which had inhibitory activity against two of the important targets associated with Alzheimer's- β-amyloid cleavage enzyme and acetylcholine esterase (AChE) . In another work, compounds having inhibitory activity against two different sites of AChE enzyme were identified . Molecular dynamics (MD) simulations have further made it easier to investigate the dynamic behavior of proteins using in silico techniques. It drastically narrows down the investment of time, cost and effort by limiting the number of cases for which experimental verification needs to be carried out . MD simulations have made it possible to study conformational changes that might occur in the three dimensional structural of proteins upon any change in the sequence. Therefore the effect of mutations occurring in a protein can be elucidated in detail using this computational approach. This can further be used to investigate the mechanism that leads to resistance of drugs in certain diseases. This insight into the mode of action of drug with respect to changes that occur in the three dimensional structure of protein can then be used to design drugs effective against both wild type and mutant isoforms of the target associated with the disease.
Pyrazinamide (PZA), a derivative of nicotinamide is one of the most imperative first-line drug treatments against tuberculosis . PZA is significantly used in MDR tuberculosis in combination with isoniazid, rifampicin and ethambutol in regimens. The most potent action of the drug is against the semi-dormant bacilli in an acidic environment, which cannot be treated with most other drugs and thus helps in shortening the chemotherapy period .
However, PZA is a pro-drug and acts only when metabolized to pyrazonic acid (POA). Pyrazonic acid is derived by the action of an enzyme called pyrazinamidase, encoded by pncA gene. Mycobacterial pyrazinamidase (MtPncA) is constitutively expressed in the cytoplasm and plays a pivotal role in the activation of pyrazinamide. Only after conversion of PZA into pyrazinoic acid by this enzyme the drug becomes potent. The bactericidal effect of this drug is attributed to destabilization of membrane potential and changes in the transport function [13–15].
Emergence of drug resistant strains of tuberculosis is one of the major concerns being faced by both developed and developing countries. Several mutations have been identified in Mtb proteins making them resistant to current drug therapies. Mutations in pncA play a major role in PZA resistance in M. tuberculosis [16–20]. The loss of PZA activity has been attributed to non-synonymous mutations in the gene pncA. Studies have shown that mutations, based on their location, were able to either decrease the activity several folds or confer resistance to PZA [19, 21–23].
Several studies have been performed to analyze the resistance mechanism of the drug. One such study analyzed three mutations namely D8G, S104R and C138Y. Study involved analysis of binding pocket, molecular docking studies and molecular dynamics simulations. This helped in understanding the interaction pattern of the drug with the enzyme and three-dimensional (3D) conformational behaviour of native and mutant PncA. The analysis predicted that the inhibition of the activity of the drug is due to change in conformation of the enzyme. According to the results, the binding cavity of the enzyme becomes rigid due to these mutations, which hampers the binding of the conversion of PZA to its active form. The major reason for rigidity is excessive hydrogen bonding between PZA binding cavity residues and their neighboring residues .
The 3D-structure of PncA protein is available in PDB (PDB ID: 3PLI). The structure was determined at a resolution of 2.2Å using PhPncA by molecular replacement . The protein consists of 185 amino acid residues as revealed by X- ray crystallography. Six stranded parallel beta sheets and four alpha helices pack together to form alpha/beta domain. Site directed mutagenesis based experiments revealed Asp 8, Lys 96 and Cys138 to be the pivotal residues involved in catalysis and Asp 49, His 57 and His 71 as the metal ion binding residues. Mutations at these residues greatly affect the function of PncA enzyme [19, 21, 25, 26]. These active residues are mostly present in the alpha- 3 helix which lies in the N-terminal region of the protein .
One of the most prominent characteristic of the mutations identified from the data obtained from clinical isolates is their diversity. These mutations range from missense mutations to insertions and deletions and are spread throughout the pncA gene [16, 27, 28]. Since the drug is metabolized in acidic conditions, it is cumbersome to obtain reliable interpretation of data on resistance of M. tuberculosis to the drug. Therefore, PCR amplification and DNA sequencing of pncA is an approach of interest as it allows rapid identification of the pncA mutations that may be involved in PZA resistance. This information can be used to study the correlation between structure and function i.e. to study the effects of mutation on the function of the enzyme. Mutational analysis thus can be used to assess PZA susceptibility.
In this study, a novel mutation Lys-96-Arg (K96R) has been considered. This mutation was revealed by deep sequencing via NGS (Next Generation sequencing) in the pncA gene of 26 random Multi Drug Resistant Tuberculosis patients of West African origin residing in USA. 9 out of 26 isolates were carrying this mutation making them resistant to pyrazinamide and in one case to all 1st and 2nd line drugs. This mutation lies in the active center of the protein and thereby is of utmost importance. Considering the ratio 9/26, the mutation seems to be quite prevalent. Since, tuberculosis is one of the most widespread disease, counteracting the causes of drug resistance via computational studies is a good choice .
Here, we have tried to answer the questions recently raised by some researches. First, what are the structural changes produced by this mutation? Second, what is the mechanism by which resistance to this drug develops due to this mutation? In the present work, we studied the relationship between structural and physico-chemical characteristics of PncA, with K96R mutation in the active site, and the enzymatic function to understand PZA resistance in M. tuberculosis. Critical understanding of these mechanisms allows the development of robust and efficient molecular diagnostic tests and provides a platform for the development of new drugs. Further, knowledge of these mechanisms can help in developing precautionary measures to curb the development of resistance .
In this study, we have performed several analyses to understand the mechanism of resistance including docking studies, binding pocket analysis and molecular dynamics of both native and mutant PncA. The MD simulations are known to provide the dynamic conformational changes in the interaction of protein and ligand . The results obtained inferred that the drug resistance could be due to the alteration in 3-D conformation of the protein by mobility of the binding residues in the mutant. The comparison of dynamic properties of native and mutant protein was done on the basis of molecular flexibility of the binding residues, the network of bonds and the fluctuations observed in them. The main objective of our work was to understand the molecular mechanism of PncA resistance at molecular level which causes PZA resistance in TB patients.
Structure of protein and ligand
We retrieved the native PncA protein structure (PDB: 3PL1) from Protein Data Bank (PDB) [25, 32]. The pdb structure contained a Fe2+ ion surrounded by one aspartate and three histidines in the substrate binding cavity along with three water molecules. The three most important functional residues include Asp 8, Lys 96 and Asp 138. In order to perform further analysis the crystallized water molecules were removed using Accelyrs Viewerlite 5.0 . The 3-D structure of the drug molecule PZA was taken from drug bank in .sdf format . Since the mutant structure was not available in the database, a point mutation was introduced in silico at 96th position in the native structure of the protein to create the mutant form.
Binding cavity analysis
The size of the binding cavity is an important measure in mutational studies. To compare the binding cavities of native and mutant proteins, we used CASTp server . CASTp algorithm is based upon recent developments in Computational Geometry. It has certain advantages which include, 1) analytical identification of pockets and, 2) precisely defined boundary between solvent and pocket, and 3) calculated parameters does not make use of dot surface or grid points and are rotationally invariant. The server uses weighted Delaunay triangulation and alpha complex for shape measurements. It generates information about surface accessible pockets and interior inaccessible cavities for a given protein. It calculates the volume and area of each cavity in solvent accessible surface (SA, Richards' surface) and molecular surface (MS, Connolly's surface). Number, area and circumference of mouth pockets are also calculated. It requires the protein .pdb structure and the results are displayed in Jmol. We used the default values for all the calculations using this server. A value of 1.4 Å was used as probe radius. The results were validated using another program, MOLE 2.0 .
Protein ligand interaction
Receptor-ligand interaction analysis was performed for native and mutant PncA with PZA. In order to carry out this docking study, we used the Glide module of Schrodinger [37, 38]. The protein structure was prepared using Schrodinger's protein preparation wizard [39, 40]. In this process hydrogen were added, bond lengths were optimized, disulfide bonds were created, terminal residues were capped and selenomethionines were converted to methionine. Flexible docking protocols were used in which both protein and ligand were kept flexible. PZA was prepared using LigPrep . It generated all possible chiral, stereochemical and ionization variants of the ligand.
A grid was generated around the active site residues of the protein molecule using Glide module of Schrodinger. The grid was generated around Asp8, Phe13, Asp49, His51, His57, Phe58, Gly60, Trp68, His71, Lys96, Ile133, Ala134, Thr135, and Cys138. Extra precision docking was then performed for native and mutant proteins with PZA.
Total ligand-receptor interaction energy along with its components like van der Waals energy and electrostatic energy was also calculated. The electrostatic energy gives a measure of how a molecule interacts with another nearby molecule. The two molecules interact with each other through their mutual electrostatic interactions. The negative value of electrostatic energy indicates better interaction and vice versa. Since hydrogen bonds and van der Waals contacts depend on complementary surfaces, these surfaces must be able to pack closely together, creating many contact points allowing charged atoms to make electrostatic bonds. Thus, polar interactions and van der Waals contribute to the dynamic stability of the ligand-receptor complex .
Molecular dynamics simulations of the docked complexes
In order to investigate the mechanism of resistance in the docked complexes, molecular dynamic study was performed in the presence of an explicit fully hydrated solvent model using triclinic boundary. The simulations were performed using Desmond Molecular Dynamics module of Schrodinger, with Optimized Potentials for Liquid Simulations (OPLS) all-atom force field 2005 [46, 47]. The complexes were prepared before simulation by addition of hydrogen atoms followed by optimization, capping of end terminals, conversion of selenomethionine to methionine and generation of disulphide bonds using the protein preparation wizard. Prepared protein-ligand complexes were then solvated with SPC water model in a triclinic periodic boundary box. To prevent interaction of the protein complex with its own periodic image, the distance between the complex and the box wall was kept 10 Å. Energy of the prepared systems was minimized for 5000 steps using steepest descent method or until a gradient threshold of 25kcal/mol/Å was achieved. It was followed by L-BFGS (Low-memory Broyden- Fletcher- Goldfarb Shanno quasi-Newtonian minimiser) until a convergence threshold of 1 kcal/mol/Å was met. For system equilibration, the default parameters in Desmond were applied. The equilibrated systems were then used for simulations at a temperature of 300 K and a constant pressure of 1atm, with a time step of 2fs. For handling long range electrostatic interactions Smooth Particle Mesh Ewald method was used whereas Cutoff method was selected to define the short range electrostatic interactions. A cut-off of 9 Å radius was used.
Continuous emergence of new mutations associated with drug resistance in Mycobacterium tuberculosis is posing a great danger to the conventional therapies. The methods need to be revised according to the current situation, which requires understanding of the mechanisms of resistance. In order to identify one such mechanism we have tried to analyze the mutation K96R occurring in the active site of Mycobacterium tuberculosis. Various analyses have been done in order to understand the mechanism at molecular level.
Preparation of mutant protein
Since the structure of mutant protein was not available in the database, in silico mutation was created at position 96 by replacing lysine with arginine. In order to allow the protein to rearrange the backbone and the side chains of the mutated residue and the nearby amino acids so as to reach the conformational state with no steric clashes and minimum internal energy, it was subjected to MD simulation for 15 ns. RMSD of each frame in reference to the first frame was plotted against the simulation time to comment on the equilibrium state of the mutated protein. Initially an RMSD of 4Å was observed for the first 6 ns after which an almost stable trajectory was obtained. The coordinates of all the frames within the most stable time frame, i.e. from 6 to 15 ns, were averaged to compute a structure representing the equilibrated state of the mutated protein. This representative average structure was used for further study.
Binding cavity analysis
The change in the cavity size was further validated using MOLE 2.0. A significant increase in the active site volume was observed for the mutant protein with this program as well.
Interaction analysis between receptor and ligand molecule
Scores from docking analysis of native and mutant protein
Vander waal Energy
All the three factors, namely, volume of binding pocket, docking score and interaction energy, typically used here to explain the reason for resistance to drug supported the idea of weak binding of the ligand with the protein in mutated form, which was hampering its activation. However, increased number of hydrogen bonds in mutant complex was suspicious. Therefore, in order to study the conformational changes occurring in the complexes in dynamic environment similar to in vivo conditions, we performed molecular dynamics simulation studies. Calculation and analysis of parameters like RMSD, radius of gyration and RMSF was done to understand the behavior of binding of the ligand with the enzyme.
Investigation of flexibility of native and mutant PncA enzyme
Some changes were observed in the interaction pattern of the native and mutant complexes post molecular dynamics studies. In native complex, the number of hydrogen bonds increased to two. The oxygen atom of PZA was making two hydrogen bonds with nitrogen of Ala 134 and Cys 138 with bond lengths of 3.33Å and 2.89Å respectively (Figure 2B). The hydrophobic interactions however reduced. PZA was now making hydrophobic contacts with Asp 8, Phe 13, Val 21, Ile 133 and His 137 of native PncA (Figure 3B). On the other hand, mutant complex was now making only one hydrogen bond with Ala 134 (Figure 2D). The bond length for the corresponding hydrogen bond was observed to be 3.06Å. Seven amino acid residues of mutant proteins, namely Asp 8, Gln 10, Phe 13, Asp 49, Ile 133, His 137 and Cys 138 were now interacting hydrophobically with PZA (Figure 3D).
RMSF values of binding cavity residues
All these analyses brought the discussion to the point that in mutated form the flap region of the protein becomes more flexible and moves away from the active site. This results in a much bigger active cavity. Once the ligand, PZA comes and binds to the active residues, flap is not able to hold it in bound state thereby preventing the activation of the prodrug.
Previously studies have been performed on PncA mutations namely, D8G, S104R and C138Y. The analysis showed that the binding cavity volume decreased and the cavity became rigid. This did not permit proper alignment of PZA inside the cavity and hampered hydrogen bond formation. Our study included a novel mutation. The mutation being present in the flap region, affected the volume of the cavity in contrast to the previous mutations. The mutation led to a substantial increase in the binding cavity. This prohibited the enzyme from holding the drug properly and therefore PZA could not take its active form .
This study is intended to present a comprehensive analysis of the drug resistance mechanism caused by K96R mutation in PncA enzyme of M. tuberculosis. Here we have performed binding pocket analysis, molecular docking and molecular dynamics simulation studies. The binding cavity analysis showed a significant increase in the volume of the active cavity from 551.9 to 1314.2 Å3 upon mutation. This fact was further supported by RMSF, radius of gyration and simulations studies. We showed that the flap region becomes highly flexible in mutated protein and moves away from the active cavity and is no longer able to hold the bound drug in its place. This prevents the activation of the prodrug to this functional form. And hence the patient with this K96R mutation develops resistance to PZA drug. This study gives insight into the mechanism by which drug resistance is acquired and therefore can be helpful while designing new drugs for tuberculosis treatment. Using the analysis, new classes of drugs can be produced which have a better conformation such that the binding to important residues is not affected by the changing volumes. Ligands should possess functional groups which permit the establishment of strongly interacting hydrogen bonds with the binding residues.
AG is thankful to Jawaharlal Nehru University for usage of all computational facilities. AG would like thank University Grants Commission, India for the Faculty Recharge position.
AG would like to acknowledge financial support from Jawaharlal Nehru University towards publication of this article.
This article has been published as part of BMC Genomics Volume 16 Supplement 2, 2015: Selected articles from the Thirteenth Asia Pacific Bioinformatics Conference (APBC 2015): Genomics. The full contents of the supplement are available online at http://0-www.biomedcentral.com.brum.beds.ac.uk/bmcgenomics/supplements/16/S2
- Onyango RO: State of the Globe: Tracking Tuberculosis is the Test of Time. J Glob Infect Dis. 2011, 3 (1): 1-3. 10.4103/0974-777X.77287.PubMed CentralView ArticlePubMedGoogle Scholar
- McQuade Billingsley K, Smith N, Shirley R, Achieng L, Keiser P: A quality assessment tool for tuberculosis control activities in resource limited settings. Tuberculosis. 2011, 91 (Suppl 1): S49-53.View ArticlePubMedGoogle Scholar
- Heifets LB, Cangelosi GA: Drug susceptibility testing of Mycobacterium tuberculosis: a neglected problem at the turn of the century. The international journal of tuberculosis and lung disease: the official journal of the International Union against Tuberculosis and Lung Disease. 1999, 3 (7): 564-581.Google Scholar
- Dooley KE, Mitnick CD, Ann DeGroote M, Obuku E, Belitsky V, Hamilton CD, Makhene M, Shah S, Brust JC, Durakovic N, et al: Old drugs, new purpose: retooling existing drugs for optimized treatment of resistant tuberculosis. Clinical infectious diseases: an official publication of the Infectious Diseases Society of America. 2012, 55 (4): 572-581. 10.1093/cid/cis487.View ArticleGoogle Scholar
- Dhiman H, Dhanjal JK, Sharma S, Chacko S, Grover S, Grover A: Resisting resistant Mycobacterium tuberculosis naturally: mechanistic insights into the inhibition of the parasite's sole signal peptidase Leader peptidase B. Biochem Biophys Res Commun. 2013, 433 (4): 552-557. 10.1016/j.bbrc.2013.03.013.View ArticlePubMedGoogle Scholar
- Dhanjal JK, Grover S, Sharma S, Singh A, Grover A: Structural insights into mode of actions of novel natural Mycobacterium protein tyrosine phosphatase B inhibitors. BMC Genomics. 2014, 15 (Suppl 1): S3-10.1186/1471-2164-15-S1-S3.View ArticleGoogle Scholar
- Goyal M, Dhanjal JK, Goyal S, Tyagi C, Hamid R, Grover A: Development of dual inhibitors against Alzheimer's disease using fragment-based QSAR and molecular docking. Biomed Res Int. 2014, 2014: 979606-PubMed CentralView ArticlePubMedGoogle Scholar
- Goyal M, Grover S, Dhanjal JK, Goyal S, Tyagi C, Grover A: Molecular modelling studies on flavonoid derivatives as dual site inhibitors of human acetyl cholinesterase using 3D-QSAR, pharmacophore and high throughput screening approaches. Medicinal Chemistry Research. 2014, 23 (4): 2122-2132. 10.1007/s00044-013-0810-2.View ArticleGoogle Scholar
- Soni S, Tyagi C, Grover A, Goswami SK: Molecular modeling and molecular dynamics simulations based structural analysis of the SG2NA protein variants. BMC Res Notes. 2014, 7 (4): 446-PubMed CentralView ArticlePubMedGoogle Scholar
- Wright A, Zignol M, Van Deun A, Falzon D, Gerdes SR, Feldman K, Hoffner S, Drobniewski F, Barrera L, van Soolingen D, et al: Epidemiology of antituberculosis drug resistance 2002-07: an updated analysis of the Global Project on Anti-Tuberculosis Drug Resistance Surveillance. Lancet. 2009, 373 (4): 1861-1873.View ArticlePubMedGoogle Scholar
- Heifets L, Lindholm-Levy P: Pyrazinamide sterilizing activity in vitro against semidormant Mycobacterium tuberculosis bacterial populations. Am Rev Respir Dis. 1992, 145 (4): 1223-1225.View ArticlePubMedGoogle Scholar
- Chang KC, Yew WW, Zhang Y: Pyrazinamide susceptibility testing in Mycobacterium tuberculosis: a systematic review with meta-analyses. Antimicrob Agents Chemother. 2011, 55 (4): 4499-4505.PubMed CentralView ArticlePubMedGoogle Scholar
- Konno K, Feldmann FM, McDermott W: Pyrazinamide susceptibility and amidase activity of tubercle bacilli. Am Rev Respir Dis. 1967, 95 (4): 461-469.PubMedGoogle Scholar
- Zhang H, Deng JY, Bi LJ, Zhou YF, Zhang ZP, Zhang CG, Zhang Y, Zhang XE: Characterization of Mycobacterium tuberculosis nicotinamidase/ pyrazinamidase. FEBS J. 2008, 275 (4): 753-762. 10.1111/j.1742-4658.2007.06241.x.View ArticlePubMedGoogle Scholar
- Yeager RL, Munroe WG, Dessau FI: Pyrazinamide (aldinamide) in the treatment of pulmonary tuberculosis. Am Rev Tuberc. 1952, 65 (4): 523-546.PubMedGoogle Scholar
- Scorpio A, Lindholm-Levy P, Heifets L, Gilman R, Siddiqi S, Cynamon M, Zhang Y: Characterization of pncA mutations in pyrazinamide-resistant Mycobacterium tuberculosis. Antimicrob Agents Chemother. 1997, 41 (4): 540-543.PubMed CentralPubMedGoogle Scholar
- Morlock GP, Crawford JT, Butler WR, Brim SE, Sikes D, Mazurek GH, Woodley CL, Cooksey RC: Phenotypic characterization of pncA mutants of Mycobacterium tuberculosis. Antimicrob Agents Chemother. 2000, 44 (4): 2291-2295.PubMed CentralView ArticlePubMedGoogle Scholar
- Jureen P, Werngren J, Toro JC, Hoffner S: Pyrazinamide resistance and pncA gene mutations in Mycobacterium tuberculosis. Antimicrob Agents Chemother. 2008, 52 (4): 1852-1854.PubMed CentralView ArticlePubMedGoogle Scholar
- Lemaitre N, Sougakoff W, Truffot-Pernot C, Jarlier V: Characterization of new mutations in pyrazinamide-resistant strains of Mycobacterium tuberculosis and identification of conserved regions important for the catalytic activity of the pyrazinamidase PncA. Antimicrob Agents Chemother. 1999, 43 (4): 1761-1763.PubMed CentralPubMedGoogle Scholar
- Sreevatsan S, Pan X, Zhang Y, Kreiswirth BN, Musser JM: Mutations associated with pyrazinamide resistance in pncA of Mycobacterium tuberculosis complex organisms. Antimicrob Agents Chemother. 1997, 41 (4): 636-640.PubMed CentralPubMedGoogle Scholar
- Sheen P, Ferrer P, Gilman RH, Christiansen G, Moreno-Roman P, Gutierrez AH, Sotelo J, Evangelista W, Fuentes P, Rueda D, et al: Role of metal ions on the activity of Mycobacterium tuberculosis pyrazinamidase. Am J Trop Med Hyg. 2012, 87 (4): 153-161.PubMed CentralView ArticlePubMedGoogle Scholar
- Pandey S, Newton S, Upton A, Roberts S, Drinkovic D: Characterisation of pncA mutations in clinical Mycobacterium tuberculosis isolates in New Zealand. Pathology. 2009, 41 (4): 582-584.View ArticlePubMedGoogle Scholar
- Quiliano M, Gutierrez AH, Gilman RH, Lopez C, Evangelista W, Sotelo J, Sheen P, Zimic M: Structure-Activity relationship in mutated pyrazinamidases from Mycobacterium tuberculosis. Bioinformation. 2011, 6 (4): 335-339.PubMed CentralView ArticlePubMedGoogle Scholar
- Rajendran V, Sethumadhavan R: Drug resistance mechanism of PncA in Mycobacterium tuberculosis. Journal of biomolecular structure & dynamics. 2014, 32 (4): 209-221.View ArticleGoogle Scholar
- Petrella S, Gelus-Ziental N, Maudry A, Laurans C, Boudjelloul R, Sougakoff W: Crystal structure of the pyrazinamidase of Mycobacterium tuberculosis: insights into natural and acquired resistance to pyrazinamide. PLoS One. 2011, 6 (4): e15785-PubMed CentralView ArticlePubMedGoogle Scholar
- Du X, Wang W, Kim R, Yakota H, Nguyen H, Kim SH: Crystal structure and mechanism of catalysis of a pyrazinamidase from Pyrococcus horikoshii. Biochemistry. 2001, 40 (4): 14166-14172.View ArticlePubMedGoogle Scholar
- Lemaitre N, Callebaut I, Frenois F, Jarlier V, Sougakoff W: Study of the structure-activity relationships for the pyrazinamidase (PncA) from Mycobacterium tuberculosis. Biochem J. 2001, 353 (Pt 3): 453-458.PubMed CentralPubMedGoogle Scholar
- Ramaswamy S, Musser JM: Molecular genetic basis of antimicrobial agent resistance in Mycobacterium tuberculosis: 1998 update. Tubercle and lung disease: the official journal of the International Union against Tuberculosis and Lung Disease. 1998, 79 (4): 3-29.View ArticleGoogle Scholar
- Daum LT, Fourie PB, Bhattacharyya S, Ismail NA, Gradus S, Maningi NE, Omar SV, Fischer GW: Next-generation sequencing for identifying pyrazinamide resistance in Mycobacterium tuberculosis. Clinical infectious diseases: an official publication of the Infectious Diseases Society of America. 2014, 58 (4): 903-904.View ArticleGoogle Scholar
- Almeida Da Silva PE, Palomino JC: Molecular basis and mechanisms of drug resistance in Mycobacterium tuberculosis: classical and new drugs. J Antimicrob Chemother. 2011, 66 (4): 1417-1430.View ArticlePubMedGoogle Scholar
- Purohit R, Rajendran V, Sethumadhavan R: Relationship between mutation of serine residue at 315th position in M. tuberculosis catalase-peroxidase enzyme and Isoniazid susceptibility: an in silico analysis. J Mol Model. 2011, 17 (4): 869-877. 10.1007/s00894-010-0785-6.View ArticlePubMedGoogle Scholar
- Berman HM, Battistuz T, Bhat TN, Bluhm WF, Bourne PE, Burkhardt K, Feng Z, Gilliland GL, Iype L, Jain S, et al: The Protein Data Bank. Acta Crystallogr D Biol Crystallogr. 2002, 58 (Pt 6 No 1): 899-907.View ArticlePubMedGoogle Scholar
- Stahle L, Ljungberg T, Rodebjer A, Ogren SO, Ungerstedt U: Differential effects of the dopamine antagonist remoxipride on apomorphine induced behaviour in the rat. Pharmacol Toxicol. 1987, 60 (4): 227-232.View ArticlePubMedGoogle Scholar
- Wishart DS, Knox C, Guo AC, Cheng D, Shrivastava S, Tzur D, Gautam B, Hassanali M: DrugBank: a knowledgebase for drugs, drug actions and drug targets. Nucleic Acids Res. 2008, 36 (Database): D901-906.PubMed CentralView ArticlePubMedGoogle Scholar
- Joe Dundas ZO, Tseng Jeffery, Binkowski Andrew, Turpaz Yaron, Liang Jie: CASTp: computed atlas of surface topography of proteins with structural and topographical mapping of functionally annotated resiudes. Nucleic Acids Res. 2006, 34: W116-W118. 10.1093/nar/gkl282.PubMed CentralView ArticlePubMedGoogle Scholar
- Sehnal D, Vareková RS, Berka K, Pravda L, Navrátilová V, Banás P, Ionescu C-M, Otyepka M, Koca J: MOLE 2.0: advanced approach for analysis of biomacromolecular channels. J Cheminformatics. 2013, 5: 39-10.1186/1758-2946-5-39.View ArticleGoogle Scholar
- Friesner RA, Banks JL, Murphy RB, Halgren TA, Klicic JJ, Mainz DT, Repasky MP, Knoll EH, Shelley M, Perry JK, et al: Glide: a new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy. J Med Chem. 2004, 47 (4): 1739-1749.View ArticlePubMedGoogle Scholar
- Small-Molecule Drug Discovery Suite 2013-1: Glide v. 2013, Schrödinger, LLC, New York, NYGoogle Scholar
- Madhavi Sastry G, Adzhigirey M, Day T, Annabhimoju R, Sherman W: Protein and ligand preparation: parameters, protocols, and influence on virtual screening enrichments. J Comput Aided Mol Des. 2013, 27 (4): 221-234.View ArticlePubMedGoogle Scholar
- Schrödinger Release 2013-1: Schrödinger Suite 2013 Protein Preparation Wizard; Epik version 2.4. 2013, S, LLC, New York, NY, Impact version 5.9, Schrödinger, LLC, New York, NY, 2013; Prime version 3.2, Schrödinger, LLC, New York, NY, 2013Google Scholar
- Schrödinger Release 2013-1: LigPrep v. 2013, Schrödinger, LLC, New York, NYGoogle Scholar
- Jones S, Thornton JM: Principles of protein-protein interactions. Proc Natl Acad Sci USA. 1996, 93 (4): 13-20.PubMed CentralView ArticlePubMedGoogle Scholar
- Schneidman-Duhovny D, Inbar Y, Nussinov R, Wolfson HJ: PatchDock and SymmDock: servers for rigid and symmetric docking. Nucleic Acids Res. 2005, W363-367. 33 Web ServerGoogle Scholar
- Duhovny D, Nussinov R, Wolfson HJ: Efficient unbound docking of rigid molecules. Algorithms in bioinformatics. 2002, Springer, 185-200.View ArticleGoogle Scholar
- Grosdidier A, Zoete V, Michielin O: SwissDock, a protein-small molecule docking web service based on EADock DSS. Nucleic Acids Res. 2011, W270-277. 39 Web ServerGoogle Scholar
- Shivakumar D, Williams J, Wu Y, Damm W, Shelley J, Sherman W: Prediction of Absolute Solvation Free Energies using Molecular Dynamics Free Energy Perturbation and the OPLS Force Field. J Chem Theory Comput. 2010, 6 (4): 1509-1519.View ArticleGoogle Scholar
- Schrödinger Release 2013-1: Desmond Molecular Dynamics System v, D. E. 2013, Shaw Research, New York, NY, Maestro-Desmond Interoperability Tools, version 3.4, Schrödinger, New York, NY, 2013Google Scholar
- Wallace AC LR: LIGPLOT: a program to generate schematic diagrams of protein-ligand interactions. Protein Eng. 1995, 8 (4): 127-134.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.