ISPUB.com / IJMB/6/1/8877
  • Author/Editor Login
  • Registration
  • Facebook
  • Google Plus

ISPUB.com

Internet
Scientific
Publications

  • Home
  • Journals
  • Latest Articles
  • Disclaimers
  • Article Submissions
  • Contact
  • Help
  • The Internet Journal of Microbiology
  • Volume 6
  • Number 1

Original Article

New Approaches to design anti-tuberculosis drugs: theoretical modeling, docking and dynamics studies on Mycobacterium tuberculosis serine hydroxymethyltransferase in complex with PLP and FFO

B Babajan, M Chaitanya, C Anuradha, C Rajasekhar, S Chitta

Citation

B Babajan, M Chaitanya, C Anuradha, C Rajasekhar, S Chitta. New Approaches to design anti-tuberculosis drugs: theoretical modeling, docking and dynamics studies on Mycobacterium tuberculosis serine hydroxymethyltransferase in complex with PLP and FFO. The Internet Journal of Microbiology. 2008 Volume 6 Number 1.

Abstract

Tuberculosis (TB) is a common and deadly infectious disease that is caused by Mycobacterium
tuberculosis (Mtb). The reemergence of TB as a public health threat has created a need to develop improved anti-mycobacterial agents. The folate pathway was found to be an attractive target for the development of anti-microbial agents. Due to indispensable role in folate metabolism Serine Hydroxymethyltransferase (SHMT) enzyme has been chosen in our present study. In this scenario we investigated the molecular model of Mtb-SHMT, by comparing with human SHMT (h-SHMT), together with enzymatic docking studies, to reveal key differences that could be useful for development of new anti-tuberculosis drugs.

 

Introduction

Tuberculosis (TB) continues to be a major cause of morbidity and mortality through out the world. It has been estimated a prevalence of 1/3 the world's population, an incidence of 9 million cases each year, and 5% of the cases are bacteria resistant to anti-TB drugs (1). The treatment of TB remains unsatisfactory. The resurgence of this disease has primarily been due to the emergence of drug resistant tubercle, especially to the most effective drug, isoniazid and rifampicin (2). Mycobacterium tuberculosis (Mtb), the main causative organism of TB, is a successful pathogen that overcomes the numerous challenges presented by the immune systems of the host. Mtb is difficult to kill for a number of reasons such as its slow growth, dormancy, complex cell envelope, intracellular pathogenesis and genetic homogeneity (3). In addition, a rising number of people in the developed world are contracting tuberculosis because their immune systems are compromised by immunosuppressive drugs, substance abuse, or AIDS. (4). There are some pharmaceuticals that can be employed in TB treatment; they explore metabolic differences between the parasite and the host, using certain enzymes as targets. However, there are several evidences indicating that these drugs have been losing their efficiency for the last decades, in most cases due to the high mutation rate in Mtb (5, 6). For this, need to develop a precise drug for TB, we focus on insilico characterization and analysis of the drug targeted enzyme in Mtb which is applicable to development of a new anti-tuberculosis drugs. It is well known that inhibitors of folate metabolism are quite important drugs, not only in the chemotherapy of TB, but also of bacterial infections and cancer (7). The effectiveness of antifolates is based on the perturbations they cause in the folate pathways, which rapidly lead to nucleotide imbalances and cell death (8, 9). This turns the enzymes thymidylate synthase (TS), dihydrofolate reductase (DHFR) and serine hydroxymethyltransferase (SHMT) involved in the folate cycle are good targets to anti-tuberculosis chemotherapy.

Serine hydroxymethyltransferase (SHMT), L-serine:tetrahydrofolate 5,10-hydroxymethyltransferase is a pyridoxyl 5`- phosphate (PLP)-dependent enzyme. Mtb-SHMT reaction plays a major role in cell physiology as it is considered to be a key enzyme in the pathway for interconversion of folate coenzymes that provide almost exclusively one-carbon fragments for the biosynthesis of a variety of end products such as DNA, RNA, ubiquinone, methionine, etc. (11, 12, 13). The physiological role of Mtb-SHMT is the reversible interconversion of serine to glycine and irreversible hydrolysis of N5, N10-methylenetetrahydrofolate to N5-methyltetrahydrofolate. In addition to these physiological reactions, Mtb-SHMT also catalyzes, in the absence of H4PteGlu, the retroaldol cleavage of several 3-hydroxyamino acids, such as allo-threonine, and the transamination and slow racemization of D- and L-alanine (14, 15). The indispensable role of Mtb-SHMT in DNA biosynthesis is valid target to develop new anti-tuberculosis drug. But the sequence similarity at active sites in the Mtb-SHMT and hSHMT is quite comparable. This fact makes the discovery of selective inhibitors for the enzyme quite difficult. On the other hand, if those inhibitors are somehow developed, it is likely that the natural selection for resistant strains would be much slower (16). And another disadvantage is very little known about Mtb-SHMT, even its crystallographic structure is yet unknown, only the gene sequence is revealed. The solution for this problem is to study the key structural difference in between hSHMT and Mtb-SHMT using computational biology; if structural difference is there it is very need full to development of new anti-tuberculosis drugs. In this paper, we propose a 3D model, molecular dynamics (MD) and docking of Mtb-SHMT with PLP and Tetrahydrofolate (FFO). The model is built based on templates available in the Protein Data Bank (PDB). The quality of the refined Mtb-SHMT structure thus obtained checked with PROCHECK program and then docking with ligands can be performed using AUTODOCK 4.0 and molecular dynamics stimulations are carried out through GROMACS.

Computitonal Methods

3-D model building

The amino acid sequence (gi: 57116827) of Mtb-SHMT was obtained from National Center for Biotechnology Information (NCBI, http:// www.ncbi.nlm.nih.gov/). Homologous sequence identity was carrier out through the BLASTp server (21), to search for short nearly exact matches and corrected the dataset to remove redundant sequences. The SHMT from Thermus thermophilus SHMT (tSHMT) (2DKJ) (17), Geobacillus Steorothermophillus SHMT (gSHMT) (1KL2) (18), B. Stearo thermophilus SHMT (bSHMT) (1YJS) (19), Human mitochondrial serine SHMT (hSHMT) (2A7V) (20) has sequence aligned with sequences of Mtb-SHMT using the ClustalW server (21. The alignments from BLAST were manually refined by comparison with the alignment of Mtb-SHMT of other species SHMT. (Fig.1) The MODELLER 9v3 program uses the spatial constrains determined from the crystal structures, to build a 3-D model of target protein with unknown tertiary structures. 100 runs were carried out to obtain the most reasonable model (22). Based on sequence alignment as explained in Figure1 and by taking default parameters one hundred models were generated and their quality assessed using the Modeller Objective Function parameter. The built Mtb-SHMT secondary structure characterization were carried out by submitting Mtb-SHMT model to PDBSUM online server (23)

Structural refinement and dynamics

Energy minimization of the model is necessary in order to relieve short contacts and correct bad geometry that may be present in the model. The best initial model obtained from the homology modeling was solvated with solvent water molecule and was roughly energy-minimized in order to make it suitable for performing molecular dynamics (MD) simulation to relax the loop and side chain (24). The MD stimulation of Mtb-SHMT was performed using GROMACS and in particular “ffG43a1” (GROMOS96) force field (25). All bonds were constrained using LINCS (26) First, 100 ps of molecular dynamics at 300 K in the water molecules inside the box to allow for the equilibration of the solvent around the protein residues. In this dynamics, all protein atoms had their positions restrained. Then, a full molecular dynamics simulation of 5000 ps at 300 K with no restrictions using 1 fs of integration time and a cut-off of 14 Å, for long-range interactions was carried out. Because these proteins have only a residual net charge, which was balanced by contra-ions for these simulations, a twin-range cutoff was used for long-range interactions: 1.5 nm for electrostatic interactions and 1.2 nm for vanderwaal interactions. The Shake algorithm was used to constrain hydrogen bond lengths (27). The effects of the electrostatic potential truncation were minor at 14 Å, and the cut-off procedure made the calculations faster than other methods. Considering that the protein and water in the simulation box can achieve about 130,000 atoms, this method can save significant computational time.

Structural validation of the Mtb-SHMT model

Stereochemical quality and structure analysis (backbone ?and dihedral angle values) for the refined Mtb-SHMT model was done with the PROCHECK ()28 WHAT IF (29) and ProSA-Web (30) programs, respectively, through the server: http://biotech.emblebi.ac.uk (Biotech for Validation Suite Protein Structures). The results were compared with the corresponding analysis of the crystallographic structures (1YJS, 2A7V, 2DKJ and 1KL2), whose coordinates were taken directly from the PDB.

Docking with PLG and FFO

After building the structure of Mtb-SHMT apoenzyme, the next step was to build the corresponding holoenzyme by docking PLG and FFO into the respective active sites. The atomic partial charges of PLG and FFO were added by online PRODRG server (31). AutoDock 4.0 (32) was used for the docking study of holoenzyme complex of Mtb-SHMT combined exploits the Lamarckian genetic algorithm.. The docking grid size was prepared with the autogrid utility of autodock setting to 60X60X60 points with a grid spacing of 0.370 A?. The grid center was placed in the active site pocket center. The grid boxes included the entire binding site of the enzyme and provided enough space for the ligand translational and rotational walk. . The consistencies of the maps were ascertained by checking the minimum and maximum values of the vanderwaal energies and electrostatic potentials for each calculated grid map. Flexible ligand docking was performed for all synthesized compounds. Docking was carried out using the empirical free energy function and the Lamarckian genetic algorithm, applying a standard protocol, the energy evaluations were 250,000, the maximum number of iterations 27,000 for an initial population of 150 randomly placed individuals with a mutation rate of 0.02, a crossover rate of 0.80, and an elitism value of 1. The number of docking runs was 100 and, after docking, the 100 solutions were clustered into groups with the RMS deviations lower than 0.5 A?. The clusters were ranked by the lowest-energy representative of each binding mode

In order to perform high frequency computational analysis such as molecular modeling, molecular dynamics and protein ligand interaction studies are performed using, a Hi-end server (Pentium IV 3.4 MHzs, AMD Athlon 64 bit, Dual processor with 1 GB video graphics card). Structural diagrams were prepared using the programs VMD (33), and Other analyses were performed using scripts included with the Gromacs distribution.

Result and Discussion

Homology modeling

Fig.1 explains final alignment, which was modeled as Mtb-SHMT in Modeller 9v3. This alignment was obtained after manual adjustments of the initial alignment from the BLAST server. Four reference proteins, hSHMT (2A7V), tSHMT (2DKJ), sSHMT (1YKJ) and rSHMT (1KL2) were used to model the structure of the Mtb-SHMT and homology scores comparing to target proteins were 60%, 57%, 57% and 45% respectively. High levels of sequence identity could guarantee more accurate alignment between the target sequence and template structure. In order to define structurally conserved regions (SCRs) of the protein family, the multi-dimensional alignment based on the structural identity was used to superimpose four reference structures, and 126 SCRs were determined (Fig.1).

Figure 1
Figure 1: Sequence alignment of Mtb-SHMT with the experimental structures gSHMT, bSHMT, tSHMT, hSHMT show the conserved region in stars (*) and deleted regions with dashes (-)

The Modeller program uses the spatial constraints, determined from the crystal structure of a template protein, to build a 3D model of the target protein with unknown tertiary structures, on the basis of amino acid sequence homology to the sequence alignment. The Ramachandran plots reveal more than 95% of the amino acid residues in the favorable regions of the plot for the whole enzyme. The main structural elements of the optimized Mtb-SHMT homology model are appearing in (Fig.2A) In secondary structure build model consists of three domains, the N-terminal domain, the second N-terminal domain and small domain, The N-terminal domain mediates inter subunit contacts and folds into two ?-helices and one ?-stand. The second N-terminal domain or large domain binds PLP, has most of the active site residues and folds into an ?-?-? sandwich containing nine confirmations, The confirmation clearly appears four beta sheets, three beta-alpha-beta motifs, four beta hairpins, one beta bulge, fifteen stands, twenty two helices, thirty five helix-helix interactions, nineteen beta turns, three gamma turns (Fig.2B).

Figure 2
Figure 2: A) Shows the refined model of Mtb-SHMT, it revels the secondary structure elements alpha helices, beta pleated sheets, loops. B) Secondary structure of Mtb-SHMT

Molecular Dynamics

The modeled Mtb-SHMT energy minimized using gromacs force fields. Later the energy minimization, the 3D Mtb-SHMT model and the decent identity crystallographic hSHMT holoenzymes were submitted to molecular dynamics simulations using the GROMOS 96 force field. This procedure had the goal of putting both enzymes in the same physiologic conditions for further superposition and also of providing a way for additional validation and refinement of the structure of Mtb-SHMT. In fact, an overall average structure conserved in time during a molecular dynamics simulation is essential to consider a model acceptable. After the dynamics simulations, Averages and fluctuations of several energies were calculated for the 5000 ps trajectory to examine stability of the Mtb-SHMT the total energy fluctuation, root mean square fluctuation (RMSF), was as small as 0.2%. The kinetic and potential energies showed fluctuations 3.5 and 0.49%, respectively (Fig.3A, B and C). The time evolution of the root mean square deviation (RMSD) of the MD trajectory from the hSHMT X-ray structure was computed to analyze a structure stability of Mtb-SHMT (not shown in Figure).

Figure 3
Figure 3: A) Variation in potential energy during the 5000 ps of MD on the Mtb-SHMT. The potential energy is averaged over 1000 ps interval. Variation in Total energy during 5000ps C) RMS fluctuation peaks in Mtb-SHMT side chain residues (D1, D2 and

The RMSD increased rapidly to 0.15 Å for the main and 0.31 Å for the side chain atoms, respectively, within the first 10 psec of the simulation (Fig.4)

Figure 4
Figure 4: Shows from the MD simulation of Mtb-SHMT, the time evolution of RMS deviation between the back bone, side chain at the 300 K initial temperature.

then, the RMSD values gradually decreased for the main and increased for side chain atoms. The total RMSD values were less than 0.31 Å indicating high structural stability of Mtb-SHMT. From the MD simulations the side chain of Lys227, His201 and Thr224 in Mtb-SHMT is oriented about at 3 to 4 Ao when it compared with the active site amino acids of hSHMT (Fig.5). This orientation make possible for efficient binding of FFO monoglutamate tail and PLP. On the other hand, the corresponding amino acid in the human enzyme, Asp146, is forced by the protein backbone to orient its carboxylate group in opposition to the monoglutamate tail of FFO, thus avoiding any possibility of interaction even if it were a positively charged side chain. Another significant difference between both enzymes is that Thr145 in hSHMT sets its side chain at about 4 to 5 Ao from FFO, with its OH group pointing towards FFO, while Asp198 in Mtb-S HMT have its carboxylate group in an opposite position, i.e., about 8 to 11 Ao from FFO. Finally, Ser193 in hSHMT and Thr224 in Mtb-SHMT set their side chains between PLP and FFO, and, in both cases, they establish hydrogen bonds with FFO, thus behaving practically in the same way. These significant orientation differences lead to development of new drugs for deadly Mycobacterium.

Figure 5
Figure 5: A) shows the active site amino acids of Mtb-SHMT (orange) and hSHMt (cyan) before refinement revels, structural orientation of active site amino acids of Mtb-SHMT after the refinement

Validation of the models

In general, the evaluation parameters of the structural model obtained for Mtb-SHMT by WHATIF PROCHECK and ProsA-Web are within the interval of values derived for their homologs and for highly refined structures RMSD Z-score values for bonds and angle parameters for the model are within values typical of highly refined structures. The fact that the RMS Z-score values of bonding distances and angles for the crystal structures are small might indicate that too-strong constrains have been used in the original refinement of 2A7V, 2DKJ, 1YJS and 1KL2. From the analysis of backbone conformations, the Mtb-SHMT interface presents only two residues located in a generous region and the remaining interface residues are in the favorable region of the Ramachandran Plot. One indication that our model is a well-refined structure is the fact that its evaluation criterion of stereo chemical and structural parameters is better than those for 2DKJ, which is the worst of the 2A7V crystal structures. The low overall RMS values for backbone superposition reflect the high structural conservation of this complex through evolution, making it a good system for homology modeling. From ProSA-Web analysis of a Mtb-SHMT protein structure shows the energy graphs having negative values correspond to stable parts of the structure (Fig.6A, B). The ProSA-web analysis of high identical template structure 2DKJand Mtb-SHMT model showed the energy graphs nearly similar to each other and in correlation with the energy pattern between the X-ray structures. The output graphs showed the Z scores,–9.69 for 2DKJ and -9.12 for Mtb-SHMT model, from this data we can consider the structure of the Mtb-SHMT interface as a good representation of the actual system.

Figure 6: (A) ProSA-web Z-score of Mtb-SHMT (black spot) revels the good quality of the model when it compared with crystal structure 2DKJ amino acid sequence position showing negative values like crystal structures2DKJ. (B) Prosa energy plot shows local model quality by plotting energies as a function of amino acid sequence position; positive values correspond to problematic or erroneous parts of the input structure. Where as the negative values corresponds to good quality of the structure. Mtb-SHMT showing 98% amino acids in negative region similar as in 2DKJ.

Active site determination

The active site determination in our model was accomplished based on its alignment to the templates. There were found degrees of identity above 95% for the active site in all the cases when the high identity template model is submitted to the PDBsum server it deliver a catalytical residues or active site amino acids, taking that it compared to our model which exhibit much more similar identities at active site amino acids region (Lys227, His 126, Thr224, Asp198and His 201). At first hand, it seems that it would be very difficult to design selective inhibitor for Mtb-SHMT, however, if this is possible, we expect that those compounds would be more durable as efficient anti-tuberculosis drugs, once the close similarity between the SHMT active sites of the different species is an indication of the low potential for mutation of the enzyme. Furthermore, Fu et al., (34) propose that the canonical SHMT structure supports the inference that the folate interacting site evolved from a type I PLP precursor enzyme through sequence insertions and not by domain swapping from other folate requiring enzymes. This may be the reason why antifolate compounds developed as chemotherapeutic agents are ineffective as inhibitors of SHMT and suggests that effective antifolate inhibitor of SHMT may not inhibit other folate enzymes. Thus, the challenge of developing selective inhibitor for SHMT would bring the additional advantage of achieving selective to the other enzyme of folate cycle.

Molecular Docking

The computer simulated automated docking studies were performed using the widely distributed molecular docking software, AutoDock 4.0. Energy minimized FFO and PLP from Dundee PRODRG2 server, were docked with Mtb-SHMT. The PLP and FFO specifically bind at active site amino acids Lys227, His 124, His127, Thr224, Asp198 and His201, give different docked energies AutoDock binding affinities of the PLP and FFO, showing binding affinity evaluated by

Table 1: Mtb-SHMT docking results with PLP and FFO

the binding free energies (?Gb, kcal/ mol), inhibition constants (Ki), hydrogen bonds, and RMSD values. The obtained success rates of AutoDock is highly excellent shown in (Table. 1), where the docked PLP binding free energies -28.43 kcal/mol and for FFO binding free enrgies -18.08 Kcal/mol (Fig.7) From the results it will revel that PLP give lowest docked energy when compared with FFO this analysis reveal that the PLP able to bind tightly with Mtb-SHMT then the FFO and showing greater binding energies, on other hand, experimental hSHMT also showing less similar interactions but due to the structural orientation in binding of hSHMT FFO and PLP is distinctive with Mtb-SHMT FFO and PLP holoenzyme complex. This unique feature may be help full to design a novel drug to against the Mycobacterium tuberculosis.

Fig.7 Active groove of Mtb-SHMT. FFO and PLP involve in binding of active site amino acids of Mtb-SHMT and the binding pocket shown in transparent solid surface with labeled amino acids (yellow) and the ligands is shown as a pink.

Conclusion

In spite of the availability of effective chemotherapy like isoniazid, rifampicin drugs and Bacille-Calmette–Guerin (BCG) vaccine, tuberculosis remains a leading infectious killer world-wide. Many factors such as, human immunodeficiency virus (HIV) co-infection, drug resistance, lack of patient compliance with chemotherapy, delay in diagnosis, variable efficacy of BCG vaccine and drugs, various other factors contribute to the mortality due to tuberculosis, for this there is need to development of a new anti tuberculosis drug. In this work we choose Mtb-SHMT, which plays a role in biosynthesis of nucleic acids, the enzyme deemed necessary for survival of Mycobacterium tuberculosis. It seems to be good target to develop a new anti-chemotherapy against mycobacterium. In this work the 3D structure Mtb-SHMT model built using homology modeling; and comparison of Mtb-SHMT with crystal structures reveals the structure similarities and conserved regions, We have also used molecular dynamics and docking studies to explore the opportunities opened by the differences found for the interactions of a monoglutamate tailed substrate with the active site of the Mtb-SHMT model and that of the crystallographic structure of hSHMT. The structural orientations of the FFO clearly indicates distinctive affinities of hSHMT FFO and Mtb-SHMT FFO, this distinctive feature helps to development of a new drugs against Mtb-SHMT and the supporting experimental studies on this data have been conducting in our lab.

References

1. Raviglione MC, Smith IM (2007) XDR tuberculosis--implications for global public health. N Engl J Med. 356: 656-659
2. Purushothaman S, Gupta G, Srivastava R, Ramu VG, Surolia A (2008) Ligand Specificity of Group I Biotin Protein Ligase of Mycobacterium tuberculosis. PLoS ONE 3: e2320. doi:10.1371/journal.pone.0002320
3. Besra G.S. , D. Chatterjee (1994) Lipids and carbohydrates of Mycobacterium tuberculosis, in: B.R. Blood (Ed.), Tuberculosis. Pathogenesis, protection, and Cantrol, ASM Press, Washington. DC: 285-306.
4. Corbett E.L, C.J. Watt, N. Walker, D. Maher, B.G. Williams, M.C. Raviglione, C. Dye (2003) The growing burden of tuberculosis: global trends and interactions with the HIV epidemic. Arch. Int. Med. 163: 1009-1021.
5. Cole, S. T. (1998) Deciphering the biology of Mtb from the complete genome sequence. Nature 393: 537-543.
6. Winstanley, P. A. (2000) Chemotherapy for falciparum malaria: the armoury, the problems and the prospects Parasitol. Today 16: 140- 146.
7. Huennekens F.M., T.H. Duffy, K.S. Vitols (1987) Folic acid metabolism and its disruption by pharmacological agents. NCI Monogr. 5: 1-8.
8. Ingraham H.A., L. Dickey, M. Goulian (1986) DNA fragmentation and cytotoxicity from increased cellular deoxyuridylate. Biochemistry 25: 3225- 3230.
9. Yoshida A., S. Tanaka, O. Hiraoka, Y. Koyama, Y. Hirota, D. Ayusawa, T. Seno, C. Garrett, Y. Wataya (1987) Deoxyribonuceloside triphosphate imbalance. 5-fluorodeoxyuridine-induced DNA double strand breaks in mouse FM3A cells and the mechanism of cell death. J. Biol. Chem. 262: 8235- 8241.
10. Chaturvedi S., V. Bhakuni (2002) Unusual Structural, Functional, and Stability Properties of Serine Hydroxymethyltransferase from Mycobacterium tuberculosis. J. Biol. Chem 278: 40793-40805.
11. Schirch, V. (1982) Serine hydroxymethyltransferase Adv. Enzymol. Relat. Areas Mol. Biol. 53: 83-112
12. Mc Neil, J., Bognar, A., and Perlman, R. (1996) In-vivo analysis of folate co-enzymes and their compartmentation in Saccharomyces cerevisiae. Genetics 142: 371-381
13. Shane, B. (1989) Folylpolyglutamate synthesis and role in the regulation of one-carbon metabolism. Vit. Horm. 45: 263-335
14. Shostak, K., and Schirch, V. (1988) Serine hydroxymethyltransferase: mechanism of the racemization and transamination of D- and L-alanine. Biochemistry 27: 8007-8017
15. Contestabile, R., Paiardini, A., Pascarella, S., di Salvo, M. L., D Aguanno, S., and Bossa, F. (2001) L-Threonine aldolase, serine hydroxymethyltransferase and fungal alanine racemase. A subgroup of strictly related enzymes specialized for different functions. Eur. J. Biochem. 268: 6508-6525
16. Tanos C., Franca C., Perdo G., Pascutti, Tedorico C., Ramalho, Jose D., Figureueoroa, A (2005) three-dimensional structure of Plasmodium falciparum serine hydroxymethyltransferase in complex with glycine and 5-formyl-tetrahydrofolate Homology modeling and molecular dynamics. Bio. Phy. Chem 115: 1-10
17. Goto, M., Kai, K., Miyahara, I., Hirotsu, K. Crystal Structure of T.th.HB8 Serine Hydroxymethyltransferase To be Published (http://www.rcsb.org/pdb/exp- lore/explore.do?structureId=2DKJ)
18. Trivedi, V., Gupta, A., Jala, V.R., Saravanan, P., Rao, G.S., Rao, N.A., Savithri, H.S., Subramanya (2002) Crystal structure of binary and ternary complexes of serine hydroxymethyltransferase from Bacillus stearothermophilus: insights into the catalytic mechanism. J.Biol.Chem. 277: 17161-17169.
19. Bhavani, S., Trivedi, V., Jala, V.R., Subramanya, H.S., Kaul, P., Purnima, K., Prakash, V., Appaji, R.N., Savithri, H.S (2005) Role of Lys-226 in the Catalytic Mechanism of Bacillus Stearothermophilus Serine Hydroxymethyltransferase-Crystal Structure and Kinetic Studies Biochemistry 44: 6929-6937.
20. Dong, A., Wu, H., Zeng, H., Loppnau, P., Sundstrom, M., Arrowsmith, C.H., Edwards, A.M., Bochkarev, A., Plotnikov, A.N. Crystal structure of human mitochondrial serine hydroxymethyltransferase 2 To be Published (http://www.rcsb.org /pdb/explore/explore.do?structureId=2A7V)
21. Altschul S. F., T.L. Madden, A.A. Schaffer, J. Zhang, Z. Zhang, W. Mile, F. Lipmsn, (1997) Gapped BLAST and PSIBLST a new generation of protein database search programes, J. Nucleic Acid Res. 25: 2714-2723.
22. Eswar Narayanan, Eramian David, Webb Ben, Shen Min-Yi, Sali Andrej (2008) Protein structure modeling with MODELLER. Methods in molecular biology. 426:145-59.
23. Laskowski R A, Chistyakov V V (2005) Thornton J M. PDBsum more: new summaries and analyses of the known 3D structures proteins and nucleic acids. Nucleic Acids Res., 33: D266-D268
24. Hamza, A., Hoon, C., Hsin-Hsiung, T., Chang-Guo, Z. (2005) Understanding human 15-hydroxyprostaglandin dehydrogenase binding with NAD+ and PGE2 by homology modeling, docking and molecular dynamics simulation Bioorganic & Medicinal Chemistry. 13: 4544-455.
25. van der Spoel D., E. Lindahl, B. Hess, A. R. van Buuren, E. Apol, P. J. Meulenhoff, D. P. Tieleman, A. L. T. M. Sijbers, K. A. Feenstra, R. van Drunen and H. J. C. Berendsen (2004) Gromacs User Manual version 3.2, www.gromacs.org
26. Hess B., H. Bekker, H.J.C. Berendsen and J.G.E.M. Fraaije(1997) LINCS: a linear constraint solver for molecular simulations. J. Comput. Chem. 18: 1463-1472.
27. Ryckaert, J.P., Ciccotti, G. and Berendsen, H.J.C. (1977) Numerical integration of the cartesian equations of motion of a system with constraints: Molecular dynamics of n-alkanes. J. Comput. Phys. 23: 327-341.
28. Laskowski R.A., M.W. MacArthur, D.S. Moss, J.M. Thornton (1993) PROCHECK: a program to check the stereochemical quality of protein structures. J. Appl. Crystallogr. 26: 283-291.
29. Vriend G (1990) WHAT IF: a molecular modeling and drug design program. J.Mol.Graph; 8: 52-6
30. Sippl M.J. (1993) Recognition of errors in three-dimensional structures in proteins. Proteins 17: 355-362
31. Schuettelkopf A. W. and D. M. F. van Aalten (2004) PRODRG - a tool for high-throughput crystallography of protein-ligand complexes. Acta Crystallographica D60: 1355--1363
32. Goodsell, D.S., Morris, G.M (1998) Automated docking using a Lamarckian genetic algorithm and an empirical binding free energy function. J. Comput Chem. 19: 1639-1662
33. Humphrey, W., Dalke, A., Schulten, K. (1996) VMD -Visual Molecular Dynamics. J. Molec. Graphics, 14: 33-38
34. Fu T.F., J.N. Scarsdale, G. Kazanina, V. Schirch, H.T. Wright (2003) Location of the pteroylpolyglutamate-binding site on rabbit cyto- solic serine hydroxymethyltransferase. J. Biol. Chem. 24: 2645- 2653.

Author Information

B. Babajan
UGC-Major Resesarch fellow,Research Scholar, Dept. of Biochemistry, Sri Krishnadevaraya university

M. Chaitanya
Research Scholar, Dept. of Biochemistry, Sri Krishnadevaraya university

C.M. Anuradha
Lecturer, Dept. of Biotechnology, Sri Krishnadevaraya university

C. Rajasekhar
Research Scholar, Dept. of Biochemistry, Sri Krishnadevaraya university

Suresh Kumar Chitta
Co-ordinator, Bioinformatics Facility, Dept. of Biochemistry, Sri Krishnadevaraya university

Download PDF

Your free access to ISPUB is funded by the following advertisements:

 

BACK TO TOP
  • Facebook
  • Google Plus

© 2013 Internet Scientific Publications, LLC. All rights reserved.    UBM Medica Network Privacy Policy