Structural Studies on Thymidylate Kinase : Evolution, Specificity and Catalysis

Abstract

Thymidylate kinase (TMK) is a key enzyme for DNA synthesis. It occurs at the junction of the de novo and salvage pathways for the synthesis of deoxythymidine triphosphate (dTTP). Its inhibition affects cell viability, thereby making it an important target for the development of anticancer, antibacterial and antiparasitic drugs. This thesis describes the analyses of the sequence, structure and dynamics of thymidylate kinase to obtain insights into its function. Two thermophilic variants of the enzyme were chosen for our studies. The studies provide valuable insights about the active site residues and the mechanism of catalysis, which have implications in protein engineering and design of specific inhibitors. Following is a chapter-wise description of the overall layout of the thesis. Chapter 1 | Introduction: This chapter provides a brief survey of the literature on TMKs and the scope of the work presented in the thesis. TMK belongs to the nucleoside monophosphate kinase (NMPK) family of enzymes, which includes adenylate kinase (AMK), guanylate kinase (GMK), uridylate kinase (UMK) and cytidylate kinase (CMK). The NMPK family of enzymes is associated with the reversible transfer of the terminal phosphoryl group from a nucleoside triphosphate (NTP) (usually adenosine triphosphate, i.e., ATP) to a nucleoside monophosphate (NMP). The identity of the NMP substrate varies among different enzymes. NMPKs share a common Rossmann fold and are comprised of a conserved P-loop, Lid region, CORE and NMP domains. The enzymes in the NMPK family also contain structurally similar active site architecture. Besides the three signature motifs, there are other conserved residues at the active site of TMK which are involved in interactions with the substrates ATP and dTMP. Despite the overall similarity, TMKs exhibit significant variations in sequence, residue conformation, substrate specificity and oligomerization mode. However, the residues responsible for these differences have not been studied. This thesis describes a comprehensive analysis of the sequence space of TMKs to detect the residues involved in such diversity. Subsequently, TMKs from a thermophilic archaeon (Sulfolobus tokodaii) and a hyperthermophilic bacterium (Aquifex aeolicus) were chosen for biochemical characterization and structural studies. Of these, the Sulfolobus tokodaii TMK (StTMK) has low sequence identity to the other TMKs of known three dimensional (3D) structures. Crystal structure analyses depicted the presence of some novel structural features and provided insights into the role of a conserved Arginine residue in function, which was verified through computational studies and mutagenesis experiments. Finally, the study on Aquifex aeolicus TMK resulted in multiple crystal structures of the apo form and different holo forms. These helped us to understand the mechanistic details of TMK-mediated catalysis, namely, the order of substrate binding and the reaction mechanism for phosphate transfer. Chapter 2 | Materials and Methods: This chapter provides a brief description of the procedures used to carry out the thesis work. The protein samples were purified to a high degree using column chromatography, and the purity was assessed using SDS-PAGE (sodium dodecyl sulfate polyacrylamide gel electrophoresis) analysis. Circular dichroism (CD) spectroscopy was employed to assess if the purified protein was well-folded. The pure and properly folded protein samples were used in further experiments. Differential scanning fluorimetry (DSF) was performed to determine the melting temperature of the thermophilic protein. MicroScale Thermophoresis (MST) and Surface Plasmon Resonance (SPR) were carried out to detect protein-substrate interactions. The protein samples were crystallized using the hanging drop vapour diffusion and microbatch under-oil techniques using commercially available crystallization screens, and the conditions which gave crystals were further optimized. Diffraction data, collected at either the home source or the synchrotron, were processed and scaled. Subsequently, phase information was obtained using the molecular replacement (MR) calculations. The MR solution was refined till convergence and its geometry was validated using different softwares. Finally, molecular dynamics (MD) simulations were performed to study the functionally important motions in the protein. Chapter 3 Insights into substrate specificity and oligomerization mode of Thymidylate Kinases from sequence evolution and conformational dynamics: Thymidylate kinase homologs exhibit significant variations in sequence, residue conformation, substrate specificity and oligomerization mode. However, the influence of sequence evolution and conformational dynamics on its quaternary structure and function has not been studied before. Based on extensive sequence and structure analyses, our study detected several non-conserved residues which are linked by co-evolution and are implicated in the observed variations in flexibility, oligomeric assembly and substrate specificity among the homologs. These lead to differences in the pattern of interactions at the active site in TMKs of different specificity. The method was further tested on thymidylate kinase from Sulfolobus tokodaii (StTMK) which has substantial differences in sequence and structure compared to other TMKs. Our sequence analyses pointed to a more flexible dTMP-binding site in StTMK compared to the other homologs, which was also indicated in MD simulations on the protein 3D structure. Binding assays proved that the protein can accommodate both purine and pyrimidine nucleotides at the dTMP binding site with comparable affinity. Additionally, the residues responsible for the narrow specificity of Brugia malayi TMK, whose three dimensional structure is unavailable, were detected. Our study provides a residue-level understanding of the differences observed among TMK homologs in previous experiments. It also illustrates the correlation among sequence evolution, conformational dynamics, oligomerization mode and substrate recognition in thymidylate kinases and detects co-evolving residues that affect binding, which should be taken into account while designing novel inhibitors. Chapter 4 | Biochemical and Structural characterization of a thermophilic variant of thymidylate kinase: This chapter reports the biochemical characterization and crystal structure determination of thymidylate kinase from the hyperthermophilic organism Sulfolobus tokodaii (StTMK) in its apo and ADP-bound forms. Our study describes the first three-dimensional structure of an archaeal TMK. The different structures had resolution ranging between 1.60 Å and 2.40 Å. StTMK is a thermostable enzyme with a melting temperature of 85.3 °C, as observed from thermal unfolding studies. The protein exists as a dimer in solution. A coupled enzyme assay, performed using thermo-stable lactate dehydrogenase (TtLdh) and pyruvate kinase (TtPk) from Thermus thermophilus, showed that StTMK has optimum activity at 80 °C. Despite the overall similarity to homologous TMKs, StTMK structures revealed several residue substitutions at the active site. However, enzyme assays demonstrated specificity to its natural substrates ATP and dTMP. A novel insertion (9 residues long) is observed in the C-terminal stretch of the Lid region. However, it is relatively rigid, which may be attributed to the presence of two proline residues and a hydrogen bond with an arginine residue in the α4/α5 loop. The C-terminus of the α2 helix points away from the Lid region in StTMK to avoid steric clashes with the Lid insertion. The main chain dihedral angles of the conserved Arg in the DRX motif are in the disallowed region of the Ramachandran plot in all holo TMK structures, wherein it forms several conserved hydrogen bonds with residues in the P-loop, α4 helix and α7 helix, as well as with the phosphate groups of both the substrates. A similar feature is observed in some of the StTMK structures. However, torsion angles in the allowed region of the Ramachandran plot are observed in one chain each in two of the apo structures. Further, conformational rearrangements of this Arg and its neighboring residues at the binding site of the second substrate are observed. The functional implication of this variation is described in the next chapter (chapter 5). Chapter 5 | Role of a conserved active site Arginine residue in Thymidylate kinase: Analysis of the structures of StTMK revealed multiple conformational states of Arg93 which is located at the reaction centre and is a part of the highly conserved DRX motif. Conformational heterogeneity of Arg can also be observed in some structures of Staphylococcus aureus and human TMK. However, the functional implication of this feature has not been probed before. The rearrangements of Arg93 are accompanied by related changes in the conformations of its neighbouring residues at the active site. This leads to three distinct conformational states in the dTMP-binding site, namely ‘Arg in’, ‘Arg intermediate’ and ‘Arg out’. Only the ‘Arg in’ state was found to be suitable for the proper positioning of the α-phosphate group of dTMP at the active site. This is hindered in the ‘Arg out’ and ‘Arg intermediate’ states. MD simulations showed that the torsion angles of the DRX Arg can sample between allowed and disallowed values in the apo-protein, with a preference for the catalytically suitable disallowed conformation in the holo-protein. Computational alanine scanning and MM/PBSA binding energy calculation further revealed the importance of Arg93 side chain in substrate binding. Subsequent site directed mutagenesis at this position to an Ala resulted in the loss of activity. Our work provides the first experimental evidence for the functional importance of Arg93 and gives insight into its regulatory role in the catalytically competent placement of dTMP. Our study also has implications for the development of potent inhibitors to lock the enzyme in the catalytically non-productive state. Chapter 6 | Characterizing active site dynamics from structural studies on the Intermediates along the reaction coordinate of a hyperthermophilic Thymidylate Kinase: TMK belongs to the family of nucleoside monophosphate kinases (NMPKs), several of which undergo structure-encoded conformational changes to perform their function. However, the absence of three dimensional structures for all the different reaction intermediates of a single TMK homolog hinders a clear understanding of its functional mechanism. We herein report the different conformational states along the reaction coordinate of a hyperthermophilic TMK from Aquifex aeolicus, determined via X-ray diffraction and further validated through normal mode studies. The analyses implicate an arginine residue in the Lid region in catalysis, which was confirmed through site-directed mutagenesis and subsequent enzyme assays on the wild type protein and mutants. Further, the enzyme was found to exhibit broad specificity towards phosphate group acceptor nucleotides. Our comprehensive analyses of the conformational landscape of TMK, together with associated biochemical experiments, provide insights into the mechanistic details of TMK-driven catalysis, for example, the order of substrate binding and the reaction mechanism for phosphate transfer. Such a study has utility in the design of potent inhibitors for these enzymes. Finally, the implications of the work described in this thesis and its future applications have been discussed in the section titled ‘Future prospects’. The work described in chapters 3 – 6 have been published in peer reviewed journals. Additionally, the author was involved in several collaborative projects which also resulted in publications (reprints attached in appendix).