No description available.

Haematological cancers like leukemia and lymphoma are characterized by genetic abnormalities, specifically chromosomal translocations. Analyses of the translocation breakpoint regions in patients have shown that some loci in the genome are more susceptible to breakage than others. However, very little is known about the mechanism of generation of many such chromosomal translocations. In the present study, we have attempted to understand the mechanism of fragility of three regions, which are prone to breaks during translocations in follicular lymphoma (FL) and T-cell leukemia. The t(14;18) translocation in FL is one of the most common chromosomal translocations. Most breaks on chromosome 18 are located at the 3’ UTR of the BCL2 gene and are broadly classified into three clusters, namely major breakpoint region (mbr), minor breakpoint cluster region (mcr) and the intermediate cluster region (icr). The RAG complex has been shown to cleave BCL2 mbr by recognizing an altered DNA structure. In the present study, by using a gel based assay, nature of the non-B DNA structure at BCL2 mbr was identified as parallel intramolecular G-quadruplex. Various studies including circular dichroism (CD), mutagenesis, DMS modification assay and 1H NMR showed the presence of three guanine tetrads in the structure. Further, evidence was also found for the formation of such a G-quadruplex structure within mammalian cells. In an effort to characterize the mechanism of fragility of mcr, a unique pattern of RAG cleavage was observed in a sequence dependent manner. Three independent nicks of equal efficiency were generated by RAGs at the cryptic sequence, “CCACCTCT”, at mcr and at a cytosine upstream of it, unlike a single specific nick at the 5’ of heptamer during V(D)J rearrangement. Interestingly, RAG nicking at mcr occured in the presence of both Mg2+ and Mn2+. Using recombination assay, followed by sequencing of the junctions, we find that mcr can recombine with standard RSS in vivo, albeit at a very low frequency. Mutations to this novel motif abolish recombination at the mcr within the cells. In order to determine the prevalence of t(14;18) translocation in the healthy Indian population, nested PCR approach followed by Southern hybridization was used. Results showed 34% prevalence of t(14;18) translocation in the Indian population. Although, no gender based difference was observed, an age dependent increase was found in adults. Further, presence of the t(14;18) transcripts was also detected. The mechanism underlying the fragility of the t(10;14) translocation involving HOX11 gene in T-cell leukemia is not known. Using primer extension assays on a plasmid DNA containing HOX11 breakpoint region, presence of consistent pause sites corresponding to two G-quadruplex forming regions, flanking the patient breakpoints, were detected. These replication blocks were dependent on K+ ions. Native gel shift assays, mutation analysis, S1 nuclease and CD, further revealed formation of intermolecular G-quadruplexes, unlike the BCL2 mbr. Further, sodium bisulfite modification assay indicated the presence of such structures in the genomic DNA within cells. Hence, we propose that two independent G-quadruplex structures formed in the HOX11 gene could interact with each other, thereby resulting in fragility of the intervening sequences, where majority of the patient breakpoints are mapped. Overall, this study has attempted to understand the role of both sequence and structure of DNA, in generating chromosomal fragility during t(14;18) translocation in FL and t(10;14) translocation in T-cell leukemia. These results may facilitate future studies in unraveling the mechanism leading to genomic instability in other lymphoid cancers.

Introduction Nucleotide excision repair (NER) is a highly conserved pathway involved in repair of a wide variety of structurally unrelated DNA lesions. One of the well characterized NER systems is from E. coli which involves UvrABC nucleases. NER consists of two related sub-pathways: global genomic repair (GGR), which removes lesions from the overall genome, and transcription coupled repair (TCR), which removes lesions from the transcribed strand of active genes. Bulky DNA lesions such as cyclobutane pyrimidine photodimers (CPD) induced by UV irradiation block RNA polymerase (RNAP) during transcription. In bacteria, a gene product of mfd called transcription repair coupling factor (TRCF) or Mfd is required for TCR. Bacterial Mfd interacts with the stalled RNAP, displaces it from the DNA and recruits NER proteins at the site of damage. Mfd, thus contributes to the faster repair of the transcribed strand compared to the non-transcribed strand for similar kind of lesions. Intracellular pathogens like M. tuberculosis are constantly exposed to a variety of stress conditions inside the host, mainly due to host defense systems and antibiotic treatments. It is therefore, extremely important for bacteria to have DNA damage repair and reversal mechanisms that can efficiently counteract these effects. However, very little is known about DNA repair systems in M. tuberculosis compared to other bacteria. Sequencing of M. tuberculosis genome revealed the presence of NER associated genes including a putative mfd. Additionally, due to the high GC content of genome as well as the DNA damage prone host environment, the transcription in M. tuberculosis may encounter the problems, which are not apparent in other bacteria. Therefore, the gene like mfd may play very important role in physiology of M. tuberculosis. In the present study, we describe the biochemical and functional characterization of Mfd from M. tuberculosis (MtbMfd) and discuss its unusual properties. Biochemical characterization of MtbMfd Genome analysis of M. tuberculosis as well as the sequence alignment studies revealed that MtbMfd is 1234 amino acids long multifunctional protein having various domains specialized for different functions. Cloning of Mtbmfd was carried out by reconstructing the full length gene from three PCR amplified fragments using genomic DNA as a template. Complementation study using Mtbmfd suggested that the gene of interest complements E. coli counterpart and increases survival of UV irradiated cells. To further characterize the function of Mtbmfd, a road block reporter assay was performed, which indicates that the MtbMfd interacts with stalled E. coli RNAP and displaces it from the site of transcription resulting in low reporter gene activity. The MtbMfd protein was expressed and purified by using various chromatographic techniques, and confirmed by mass spectrometry. In addition to full length protein, a number of truncated MtbMfd constructs were generated and purified to homogeneity. Mfd is a motor protein and requires ATP hydrolysis in order to translocate along DNA. The signature motifs of superfamily 2 helicases / ATPases are present at the C-terminal of Mfd along with translocase motif which is highly homologous to motif present in RecG helicase. To analyze the kinetics of ATP hydrolysis of MtbMfd and its truncated proteins, ATPase reactions were carried out using γ32P-ATP as a tracer. Wild-type MtbMfd exhibited ATPase activity, which was stimulated ~1.5 fold in presence of dsDNA. The mutant MtbMfd (D778A), which harbors mutation in one of the key residues of Walker B motif of the ATPase domain showed negligible ATPase activity indicating the importance of residue D778 for ATP hydrolysis. While the C-terminal domain (CTD) comprising amino acids 600 to 1234 showed elevated ATPase activity, the N-terminal domain (NTD) containing the first 500 amino acid residues was able to bind ATP but deficient in hydrolysis. Deletion of 184 amino acids from the C-terminal end of MtbMfd (MfdΔC) increased the ATPase activity by ~10-fold compared to full-length MtbMfd. The translocase activity of MtbMfd was measured by an oligonucleotide displacement assay and it was found that full length MtbMfd and CTD have a very weak translocase activity whereas, MfdΔC exhibited efficient translocation along DNA in ATP dependent manner. These results provide a direct correlation between translocase and ATPase activity of MtbMfd, and suggest possibly an auto-regulatory function for the extreme C-terminus of MtbMfd. Oligomeric status of MtbMfd was determined using various techniques including gel filtration chromatography and it was found that MtbMfd exists as monomer and hexamer in solution. The monomer showed increased ATPase activity and susceptibility to proteases compared to the hexameric form. MfdΔC, on the other hand, was predominantly monomer in solution implicating importance of the extreme C-terminal region in oligomerization of protein. Taken together, the biochemical evidence suggests that monomeric MtbMfd is an active form and oligomerization provides stability to the protein. One important finding of the present study is the binding of ATP to NTD of MtbMfd. All Mfd NTDs resemble UvrB and possesses the degenerate ATPase motifs. Indeed, on the basis of sequence and structural similarities, it has been suggested that Mfds have evolved from UvrB incorporating an additional translocase activity. UvrB has a cryptic ATPase activity while the NTD of Mfd may have lost the activity as it possesses degenerate Walker motifs. In contrast, NTD of MtbMfd binds ATP but is hydrolysis deficient. A closer comparison of the amino acid sequences in the Walker A motif reveal that conserved K 45 of UvrB has been replaced by R in case of NTD of MtbMfd. It has been shown previously that mutation of K 45 to A, D and R led to a loss of ATPase activity of UvrB. Thus, MtbMfd seems to be a natural mutant of UvrB. Since NTD harbors an intact UvrA interacting domain, when it is expressed it may sequester the cellular pool of UvrA leading to dominant negative phenotype. When UV survival assays were carried out, cells expressing NTD showed hyper-sensitivity to UV light – a typical characteristic of NER deficiency. In addition, in vitro NER assay clearly suggested that NTD sequesters pool of UvrA inside the cell and blocks both GGR and TCR which further affects the mutation frequency of bacterial cells. Influence of MtbMfd on elongation state of RNAP The movement of RNAP along the template during transcription elongation is not uniform and is interrupted due to various factors. To overcome transcription elongation interruptions, a number of proteins viz. Mfd, Gre and Nus act on RNAP and modify its activity. RNAP displacement and transcript release experiments showed that MtbMfd influenced the elongating RNAP by more than one way. MtbMfd displaced stalled RNAP, which was blocked by NTP starvation on T7A1 promoter based template in a concentration and time-dependent manner. RNAP displacement activity of MtbMfd was shown to depend on ATP or dATP hydrolysis. On the other hand nucleotides like ADP, GTP, CTP and ATPγS did not support the RNAP displacement activity. However, in presence of ATPγS, MtbMfd was able to bind stalled complex but unable to displace RNAP suggesting that ATP or dATP hydrolysis is important for MtbMfd function. On the other hand, MtbMfd did not affect initiating RNAP when σ factor was still bound suggesting that upstream DNA is necessary for Mfd function. To assay RNA or transcript release activity of MtbMfd after transcription complex disruption, immobilized transcription complex assay was carried out. Immobilized stalled complex was generated by UTP and CTP starvation on biotinylated T7A1 promoter based template which can be affixed to temporary pellet in presence of streptavidin beads. It was found that MtbMfd released RNA into a supernatant fraction in a concentration-dependent manner suggesting that MtbMfd releases transcript after ternary complex disruption. MtbMfd released transcript in an energy-dependent manner and both ATP and dATP supported the activity, which allows the complete separation of RNA release from RNA synthesis inside the cell. An ATPase mutant of MtbMfd (MfdD778A) failed to release transcript, which further supported that ATP hydrolysis is important for MtbMfd function. Since both Mfds and RNAPs are evolutionary conserved proteins, to analyze the effect of MtbMfd on other bacterial RNAPs, displacement and release assays were carried out. Stalled complexes were generated using EcoRNAP (E. coli), MsRNAP (M. smegmatis) and MtbRNAP (M. tuberculosis) on T7A1 promoter based template. It was observed that MtbMfd was able to displace all the three RNAPs from stalled elongation complex as well as released transcript with varying efficiency. MtbMfd showed optimal displacement and release activity in presence of mycobacterial RNAPs. Transcription elongation complexes adopt various conformations and exist as different isomerized states during elongation. In an active elongation complex the 3'-OH polymerizing end of transcript aligns with an active centre of the RNAP. However, one of the most common and intrinsic properties of RNAP is backtracking or reverse translocation, which leads to misalignment of 3'-OH polymerizing end from an active centre of the polymerase. It is of interest to know if backtracking affects MtbMfd function. It is likely that complexes blocked by lesions inside the cell might tend to backtrack, and different translocational isomers possibly have different sensitivities to MtbMfd action which may illuminate the overall mechanism of MtbMfd. Backtracking of RNAP was induced on +20 and +39 stalled complexes and the effect of MtbMfd was analyzed in presence of NTPs in the reaction. It was found that arrested or backtracked complexes were restored to the forward position by the activity of MtbMfd in presence of NTP resulting into productive elongation. These results suggest that arrested RNAP again resumes transcription if conditions are favorable; otherwise, MtbMfd further assists RNAP to dissociate which leads to release of transcript. Anti-backtracking activity of MtbMfd might have important function in cellular metabolism and it has been speculated that Mfd could play more general role during transcription apart from repair. To explore the role of MtbMfd as a transcription factor and effect of MtbMtb on transcription processes in the mycobacteria, a variety of T7A1 promoter based templates were generated. These templates were derived from genes of M. tuberculosis and E. coli having varying GC content (39-81 %). The rationale behind this experiment is that the high GC content of mycobacteria and the template derived from mycobacterial genes may pose as sequence dependent structural constraints and hence block the RNAP during transcription. By anti-backtracking activity of MtbMfd these paused complexes may get relieved, leading to efficient transcription by RNAP which may lead to the formation of more full length transcript. To analyze the effect of MtbMfd, purified templates of different GC content were incubated with RNAP and MtbMfd to carry out in vitro transcription. Although, in case of multiple rounds of transcription, multiple pauses were observed even in presence of MtbMfd. However, in presence MtbMfd around 1.5 - 2 fold increased full-length transcripts were observed suggesting that MtbMfd assisted RNAP during elongation to overcome sequence dependent pause. To avoid multiple pauses that are likely to occur due to the initiation of multiple round of transcription, and trailing effect of RNAP itself, single round of transcriptions were carried out in presence of heparin. Sequence specific pauses were observed with increasing GC percentage in template suggesting that indeed high GC content contributes to transcription pause. At the same time, MtbMfd in the reaction increased the amount of full length transcript by 1.5 - 2.0 fold probably by pushing paused RNAP forward to resume elongation. Taken together, this study investigates the biochemical properties of MtbMfd and its mechanism of action. In addition, it explores the importance of the coupling of transcription to repair in M. tuberculosis as well as the overall proof reading mechanism of transcription elongation in the GC rich genome of mycobacteria.

Tuberculosis continues to be one of the major causes of morbidity and mortality worldwide. Several mycobacterial species such as M. tuberculosis and M. africanum are responsible for causing this disease in humans. Reports of high cAMP levels in mycobacterial species (as compared to other bacteria such as E. coli) suggested that this second messenger may play an important role in the biology of mycobacteria. Further, it was reported that infection with mycobacteria led to an increase in the cAMP levels within the host macrophage. More recent studies have shown that this cAMP increase may be due to bacterially derived cAMP, hinting at a role for cAMP in mycobacterial pathogenesis. Given this background, the study of cAMP in mycobacteria proves to be an interesting field of research. Signalling through cAMP involves an interaction of this cyclic nucleotide with a cAMP-binding protein. These proteins typically contain a cyclic nucleotide-binding domain (CNB domain) linked to another (effector) domain. The CNB domain is thought to allosterically control the activity of the effector domain, thus mediating cellular responses to altered cAMP levels. For example, in the case of eukaryotic protein kinase A (PKA), binding of cAMP to the CNB domain results in relieving the inhibitory effects of the regulatory subunit on the catalytic subunit. The catalytic subunit then phosphorylates its target substrates, eliciting a variety of cellular responses. This work involves the characterisation of novel cAMP-binding proteins from mycobacteria, in an attempt to better understand cAMP signalling mechanisms in these organisms. The genome of M .tuberculosis H37Rv is predicted to code for ten CNB domain-containing proteins. One of these genes is Rv0998 (KATmt). KATmt was found to contain a GCN5 related N-acetyltransferase (GNAT) domain linked to a CNB domain. KATmt finds orthologues throughout the genus Mycobacterium, thereby suggesting its role in the basic physiology of these organisms. In addition, such a domain fusion is unique to mycobacteria and hence promises to deliver insights into the biology of this medically important genus. Presented here are the biochemical and functional characterisation of KATmt and its orthologue from M. smegmatis, MSMEG_5458 (KATms). Recombinant KATms bound cAMP with high affinity, validating the functionality of its CNB domain. Mutational and analogue-binding studies showed that the biochemical properties of the CNB domain were similar to mammalian protein kinase A and G-like CNB domains. The substrate for the GNAT acetyltransferase domain was identified to be a universal stress protein from M. smegmatis (MSMEG_4207). MSMEG_4207 was acetylated at a single lysine residue (Lys 104) by KATms in vitro. Further, cAMP binding to KATms increased the initial rate of acetylation of MSMEG_4207 by 2.5-fold, suggesting allosteric control of acetyltransferase activity by the CNB domain. To ascertain that KATms acetylated MEMEG_4207 in vivo, an in-frame deletion of the KATms gene was generated in M. smegmatis (ΔKATms). MSMEG_4207 was immunoprecipitated from wild-type M. smegmatis and the ΔKATms strains, followed by mass spectrometric analysis. Acetylated MSMEG_4207 was only present in the wild-type strain, confirming that KATms and MSMEG_4207 is an in vivo enzyme-substrate pair. Key biochemical differences were observed between KATms and KATmt. KATmt had an affinity for cAMP in the micromolar range, close to three log orders lower than that of KATms. In addition, KATmt showed strictly cAMP-dependent acetylation of MSMEG_4207. This demonstrates that orthologous proteins often evolve under varied selective pressures, resulting in divergent properties. Using a combination of bioluminescence resonance energy transfer (BRET) and amide hydrogen/deuterium exchange mass spectrometry (HDXMS), the conformational changes that occur upon cAMP binding to the CNB domain of KATms were monitored. A BRET-based conformation sensor was constructed for KATms by inserting KATms between GFP2 (green fluorescent protein) and Rluc (Renilla luciferase). An increase in BRET upon cAMP binding to the sensor was observed. HDXMS analysis revealed that besides the CNB domain, the only other region that showed conformational changes in KATms upon cAMP-binding was the linker region. To confirm that the linker region was important in propagating the effects of cAMP-binding to the acetyltransferase domain, an additional construct for BRET analysis encompassing the CNB domain and the linker region was generated. The magnitude of the increase in BRET was similar to the full length BRET-based sensor, validating the crucial role of the linker region in propagating cAMP-mediated conformational changes. A ‘PXXP’ motif found in the linker region, showed maximum exchange in HDXMS analysis. Mutation of both these proline residues to alanine in KATms, as well as KATmt, resulted in decoupling of cAMP-binding and allosteric potentiation of acetyltransferase activity. In contrast to the intricate parallel allosteric relays observed in other CNB domain-containing proteins, the CNB domain in KATms functions as a simpler cyclic nucleotide binding-induced switch involving stabilization of the CNB and linker domain alone. Therefore, KATms is an example of a primordial CNB domain where conformational changes are a consequence of binding-induced ordering alone. Using a computational approach, putative substrate proteins of KATmt from M. tuberculosis were identified. The substrate specificity of lysine acetyltransferases is determined loosely by a consensus sequence around the lysine residue which is acetylated. Using this property of protein acetyltransferases, the genome of M. tuberculosis H37Rv was mined for proteins harboring lysine residues in a similar sequence context as seen in MSMEG_4207. In vitro biochemical analysis of some of the predicted substrates helped confirm a subset of enzymes belonging to the fatty acyl CoA synthetase (FadD) class as substrates of KATmt. The acetylation of FadDs by KATmt was cAMP-dependent. In each of the four proteins tested, acetylation was found to occur at a single conserved lysine residue. To confirm that FadDs were acetylated by KATmt in vivo, BCG_1055, the orthologue of KATmt in M. bovis BCG, was deleted using the specialised transduction method. FadD13, one of the FadDs acetylated by KATmt in vitro, was immunoprecipitated from wild-type M. bovis and the ΔBCG_1055 strains using a FadD13-specific polyclonal antibody. Acetylated FadD13 was almost completely absent in ΔBCG_1055 but substantial amounts of acetylated FadD13 were present in the wild-type strain, indicating that FadD13 was indeed an in vivo substrate of KATmt. The functional consequences of acetylation of FadDs were analysed using an in vitro fatty acyl CoA synthetase assay. The activities of FadD2 and FadD13 were inhibited on acetylation with KATmt, while acetylation of FadD5 resulted in the formation of a novel product. Therefore, modification of the highly conserved lysine residue in these enzymes by acetylation led to loss or alteration of their enzymatic activity, suggesting that acetylation may be used as a regulatory mechanism to modulate the activities of some of the FadDs by KATmt in a cAMP-dependent manner. Given the extensive role of FadDs in cell wall biosynthesis and lipid degradation in mycobacteria, it seems possible that post-translational control by KATmt in a cAMP-dependent manner constitutes a novel mechanism utilised by these bacteria to regulate these pathways. This direct regulation of protein lysine acetylation by cAMP appears to be unique to mycobacteria, as orthologues of KATmt are not found outside this genus. In addition, the biochemical differences between KATmt and its orthologue from M. smegmatis KATms, indicate species specific variation, on a common theme. This study is the first report of protein lysine acetylation in mycobacteria. In addition to the identification of several proteins subject to this post-translational modification, the effect of acetylation on the enzymatic activities of some of them has been elucidated.

Studies on stationary phase cultures of microorganisms in the laboratory have helped us understand the genetic and physiological basis of adaptation in their natural habitats. Stasis or decline in bacterial populations due to nutrient depletion during stationary phase has been shown to lead to the selection of mutants that are able to survive better than their parent. This phenomenon, where a mutant exhibits relatively better growth than its immediate parent, under the conditions that prevailed during its appearance, is termed as Growth Advantage in Stationary Phase (GASP) and constitutes a general strategy to survive prolonged stationary phase. A mutation conferring growth advantage in one environment, however, can result in the loss of fitness in another environment, typifying the environment-specificity of the adaptation. Several GASP loci have been isolated over the last two decades and they have enhanced our understanding of the rapid evolution of microorganisms under nutrient starvation and the forces driving it. The results presented in this thesis detail the isolation and characterization of one such survivor from the culture of an E. coli strain ZK819 that was grown for 28 days without any additional nutrients. The strain ZK819 carries the rpoS819 allele that was isolated as the first GASP mutation resulting in attenuated activity of the stationary phase - specific sigma factor RpoS. Among the several possible mutations accumulating in the survivors over this long incubation period, one was identified to be within the hns gene encoding the global regulator H-NS. Expression of a majority of genes regulated by H-NS is environmentally modulated and their products are required under various stress conditions. The presence of the hns66 mutation was identified by the derepression of the bgl operon, one of the well-studied targets of H-NS mediated repression. The primary question being addressed in this study is whether this naturally isolated hns66 mutation has any role in the survival of its bearer during prolonged stationary phase. The hns66 mutation results in a partially active longer polypeptide due to a single nucleotide alteration within the stop codon of the H-NS ORF. Strains carrying the hns66 allele showed differential fitness under two different environments relative to the parent. In early stationary phase there is a strong disadvantage associated with the hns66 allele in the rpoS819 background, when competed against the parent, highlighted by the observation that there is a sharp fall in the number of mutants and a concomitant appearance of revertants that lead to the rapid displacement of the original mutant. The reversion is predominantly due to suppressor mutations within the ssrA locus encoding a tmRNA involved in the release of stalled ribosomes and degradation of aberrant proteins. Analysis of the suppression phenomenon revealed that the altered H-NS66 protein is a substrate for SsrA-mediated tagging for proteolysis. The H-NS66 protein is sub-optimally functional and can repress the bgl operon when its level is increased either by over-expression from a plasmid or a mutation in ssrA. Mutation in ssrA suppressed the strong growth disadvantage displayed by the original survivor in early stationary phase. The growth disadvantage of the hns66 survivor was also lost upon introduction of the wild type rpoS allele, indicating that the combined presence of the hns66 and rpoS819 alleles is detrimental in early stationary phase. In the presence of the rpoS+ allele, the hns66 mutant showed a modet GASP phenotype relative to the parent in early stationary phase. In contrast to the crash observed during early stationary phase, strains carryingthe hns66 allele fared better in late stationary phase conditions compared to the parent. When present in its original genetic context, the hns66 allele conferred growth advantage to the original survivor in late stationary phase conditions when competed against the parent grown for 21 days. The hns66 allele, when introduced in the parent strain ZK819 by transduction could still confer a strong GASP phenotype when competed against one-day old parent cells in medium derived from a 21 day old culture. These observations suggest that the hns66 allele enables the cells to scavenge for limited quantities of specific nutrients available in the aged medium. This is consistent with the observation that the advantage is seen when the mutant is in minority and is lost when the mutant is in majority. These studies highlight the fact that the conditions under stationary phase are constantly changing and to adapt under these rapidly changing conditions, modulation of an already existing function would be a preferred strategy than selection for a new function or the complete loss of a function. In this context, genetic modulation of a global regulator offers a better option due to the pleiotropic effects it may generate. This study adds hns to the list of genes, which upon mutation confer a GASP phenotype to cells. It also highlights that the age of the cells as well as the medium has an impact on the GASP phenotype. Elucidation of the mechanism of the H-NS mediated growth advantage and its relation to the status of the rpoS locus can contribute to our understanding of the dynamics of the stationary phase.

RNA polymerase II (Pol II) is the enzyme responsible for the synthesis of all mRNAs in eukaryotic cells. As the central component of the eukaryotic transcription machinery, Pol II is the final target of transcription regulatory pathways. While the role for different Pol II associated proteins, co-activators and general transcription factors (GTFs) in regulation of transcription in response to different stimuli is well studied, a similar role for some subunits of the core Pol II is only now being recognized. The studies reported in this thesis address the role of the fourth largest subunit of Pol II, Rpb4, in transcription and stress response using Saccharomyces cerevisiae as the model system. Rpb4 is closely associated with another smaller subunit, Rpb7 and forms a dissociable complex (Edwards et al. 1991). The rpb4 null mutant is viable but is unable to survive at extreme temperatures (>34ºC and <12ºC) (Woychik and Young, 1989). This mutant has also been shown to be defective in activated transcription and unable to respond adequately to several stress conditions (Pillai et al. 2001; Sampath and Sadhale, 2005). In spite of wealth of available information, the exact role of Rpb4 in transcription process remains poorly understood. In the present work, we have used genetic, molecular and biochemical approaches to understand the role of Rpb4 as described in three different parts below: I. Role of Rpb4 in various pathways related to Transcription Elongation The genome-wide recruitment study of RNA pol II in presence and absence of Rpb4 has indicated role of Rpb4 in transcription elongation (Verma-Gaur et al. 2008). However, a recent proteomics based report has argued against it (Mosley et al. 2013). To address this conflict and understand Rpb4 functions, we monitored recruitment of RNA pol II on a few individual long genes in wild type and rpb4∆ cells. It was observed that RNA pol II recruitment on genes with longer coding regions is not significantly affected in rpb4∆ as compared to wild type thus ruling out role of Rpb4 in transcription elongation of these genes. However, our genetic interaction studies have shown a strong interaction (synthetic lethality) between RPB4 and the PAF1 and SPT4 genes, the products of which code for well-known transcription elongation factors. The studies based on Rpb4 overexpression in mutants for elongation factors, 6-Azauracil sensitivity of cells, effect of Dst1 overexpression in rpb4∆ cells and mitotic recombination rate in rpb4∆ cells have indicated functional interactions of Rpb4 with many of the transcription elongation factors. II. Studies on Genetic and Functional Interactions of Rpb4 with SAGA Complex in Promoter- Specific Transcription Initiation To carry out transcription, RNA pol II depends on several general transcription factors, mediators, activators, co-activators and chromatin remodeling complexes. In the present study, we explored the genetic and functional relationships between Rpb4 and the SAGA complex of transcription machinery, to gain some insight on the role of Rpb4 during transcription. Our chromatin immunoprecipitation data suggest that RNA pol II does not associate with promoters of heat shock genes during transcription activation of these heat stress induced genes in absence of Rpb4. SAGA coactivator complex is required for RNA pol II recruitment and transcription activation of these genes (Zanton and Pugh, 2004). However, recruitment of the SAGA complex at promoters of these heat shock genes was not affected in rpb4∆ cells after heat stress. Our genetic interaction analysis between RPB4 and components of SAGA complex (spt20∆) showed synthetic lethality indicating that fully functional Rpb4 and SAGA complex are required for cellular functions in the absence of heat stress and the simultaneous deletion of factors in the two complexes leads to cell death. III. Role of Rpb4 in phosphorylation cycles of Rpb1-CTD The C-Terminal Domain (CTD) of Rpb1 protein of RNA pol II undergoes several rounds of phosphorylation cycles at Ser-2 and Ser-5 residues on its heptad repeats during transcription. These phosphorylation marks are to be erased before the start of next round of transcription. Using protein pull down assay, we observed that hyperphosphorylated form of Rpb1 is reduced in rpb4∆ as compared to that seen in wild type cells among the free RNA pol II molecules. The level of Rpb2 protein was unaffected in both wild type and rpb4∆. These preliminary data hints at role of Rpb4 in the regulation of Rpb1 phosphorylation.

Biological processes are so complicated that to understand the mechanisms underlying the functioning of biomolecules it is inevitable to study them from various perspectives and with a wide range of tools. Understanding the function at the molecular level obviously requires the knowledge of the three dimensional structure of the biomolecules. Experimentally this can be obtained by techniques such as X‐ray crystallography and NMR studies. Computational biology has also played an important role in elucidating the structure function relationship in biomolecules. Computationally one can obtain the temporal as well as ensemble behavior of biomolecules at atomic level under conditions that are experimentally not accessible. Molecular dynamics(MD) study is a technique that can be used to obtain information of the dynamic behavior of the biomolecules. Dynamics of large systems like proteins can be investigated by classical force fields. However, the changes at the level of covalent bond involve the reorganization of electron density distribution which can be addressed only at Quantum mechanical level. In the present thesis, some of the biological systems have been characterized both at the classical and quantum mechanical level. The systems investigated by MD simulations and the insights brought from these studies are presented in Chapters 3 and 4. The unusual bonds such as pyrophosphate linkage in ATP and short strong hydrogen bonds in proteins, investigated through high level quantum chemical methods, are presented in Chapters 5, 6 and 7. Part of this thesis is aimed to address some important issues related to the dynamics of Tryptophanyl tRNA synthetase (TrpRS) which belongs to classic of aminoacyl‐tRNA synthetases (aaRS). aaRSs are extremely important class of enzymes involved in the translation of genetic code. These enzymes catalyze the aminoacylation of tRNAs to relate the cognate amino acids to the anticodon trinucleotide sequences. aaRSs are modular enzymes with distinct domains on which extensive kinetic and mutational experiments as well as structural analyses have been carried out, highlighting the role of inter‐domain communication (Alexander and Schimmel, 2001). The overall architecture of tRNA synthetases consists of primarily two domains. The active site domain is responsible for the activation of an amino acid with ATP in synthesizing an enzyme‐bound aminoacyl‐adenylate, and transfer of the aminoacyl‐adenylate intermediate to the 3’end of tRNA. The second domain is responsible for selection and binding of the cognate tRNA. aaRSs are allosteric proteins in which the binding of tRNA at the anticodon domain influences the activity at the catalytic region. These two binding sites are separated by a large distance. One of the aims of this thesis is to characterize such long distance communication (allosteric communication) at atomic level in Tryptophanyl tRNA synthetase. This is achieved by generating ensembles of conformations by MD simulations and analyzing the trajectories by novel graph theoretic approach. Graph and network based approaches are well established in the field of protein structure analysis for analyzing protein structure, stability and function (Kannan and Vishveshwara, 1999; Brinda and Vishveshwara, 2005). The parameters such as clusters, hubs and shortest paths provide valuable information on the structure and dynamics of the proteins. In this thesis, network parameters are used for the analysis of molecular dynamics MD) simulation data, to represent the global dynamic behavior of protein in a more elegant way. MD simulations are performed on some available (and modeled) structures of TrpRS bound to a variety of ligands, and the protein structure networks( PSN) of non‐covalent interactions are characterized in dynamical equilibrium. The ligand induced conformational changes are investigated through structure networks. These networks are used to understand the mode of communication between the anticodon domain and the active site. The interface dynamics is crucial for the function of TrpRS (since it is a functional dimer) and it is investigated through interface clusters. The matter embodied in the thesis is presented as 9 chapters. Chapter 1 lays the suitable background and foundation for the study, surveying relevant literature from different fields .Chapter 2 describes in detail the various materials, methods and techniques employed in the different analyses and studies presented in this thesis. A brief description of well‐known methods of molecular dynamics simulations, essential dynamics calculations, cross correlation maps, conformational clustering etc.is presented. The methods for constructing protein structure graphs and networks, developed in our lab, are described in detail. The use of network parameters for the analysis of MD simulation data to address the problem of communication between the two distal sites is also presented. Some descriptions of the ab initio quantum mechanical methods, which are used to investigate the unusual bonds in biomolecules, are also presented in this chapter. Chapter 3 is devoted in discussing the results from several normal as well as high temperature MD simulations of ligand‐free and ligand bound Bacillus stearothermophilus Tryptophanyl‐tRNA synthetase (bsTrpRS). The essential modes of the protein in the presence of different ligands are captured by essential dynamics calculations. Different conformations of the protein associated with the catalysis process of TrpRS, as captured through experiments, are discussed in the context of conformational sampling. High temperature simulations are carried out to explore the larger conformational space. Chapter 4 is focused on the results obtained from the MD simulation of human Tryptophanyl‐tRNA synthetase (hTrpRS). The structure of human TrpRS bound to the activated ligand (TrpAMP) and the cognate tRNA(tRNATRP) is modeled since no structure in the presence of both TrpAMP and tRNATRP is available. MD simulations on these modeled as well as other complexes of hTrpRS are performed to capture the dynamical process of ligand induced conformational changes (Hansiaetal., communicated). Both the local and the global changes in the protein conformation from the protein structure network (PSN) of MD snapshots are analyzed. Several important information such as the ligand induced correlation between different residues of the protein, asymmetric binding of the ligands to the two subunits of the protein, and the path of communication between the anticodon region and the aminoacylation site are obtained. Also, the role of the dimmer interface, from a dynamic perspective, is obtained for the first time. The interface dynamics which stabilize different quaternary structures of lectins (with high sequence and structure similarity) were investigated in a collaborative work (Hansiaetal.,2007). The lectin peanut agglutinin (PNA) is a tetramer with three different types of interfaces. The interface dynamics of this protein in the presence and in the absence of metal ions was investigated and the paper reporting the results from this study is included as appendix in this thesis. Chapter 5 deals with high level ab initio quantum chemical calculations on tri‐ and diphosphate fragments of adenosine triphosphate (ATP). Pyrophosphate prototypes such as methyl triphosphate and methyl diphosphate molecules in their different protonation states have been investigated at high levels of calculations (Hansiaetal., 2006a). The optimized geometries, the thermochemistry of the hydrolysis and the molecular orbitals contributing to the high energy of these compounds have been analyzed. These investigations provide insights into the‘‘highenergy’’character of ATP molecule. Further, the dependence of vibrational frequencies on the number of phosphate groups and the charged states has also been presented. These results aid in the interpretation of spectra obtained by experiments on complexes containing pyrophosphate prototypes. Hydrogen bonding is fundamental in understanding the structure and properties of molecules of biological interest including proteins. A recent analysis carried out in our lab showed that a significant number of short hydrogen bonds (SHB) are present in proteins (Rajagopal and Vishveshwara, 2005). Chapters 6 and 7 elucidate the results obtained from ab initio quantum chemical calculations on some of these SHBs to get aquantitative estimation of their geometry and strength. In chapter 6, asystematic analysis of the geometries and the energetics of possible SHB systems, which are frequently encountered in proteins, are presented at different levels of theory (HF,DFTandMP2). It is found that the SHBs involving both charged residues in the proteins are intrinsic in nature. However, two neutral residues form a SHB in the protein crystal structures either due to geometric constraints or due to the environment of these residues. This analysis enables one to distinguish SHBs which are formed because of geometric constraints from those which are formed because of the inherent property of the chemical groups involved in the hydrogen bonding. These results are useful in refining protein structures determined by crystallographic or NMR methods. In addition, sulfur atom of methionine and cysteinein proteins also participate in SHBs, which are not so well characterized. Chapter 7 presents the similar analysis carried out on short hydrogen bonds in proteins involving sulfur atom. A detailed analysis of SHBs of sulfur containing groups in a data set of proteins has been carried out. Some of the residue pairs from this analysis were considered for ab initio calculations. However, the optimization of these examples resulted in breaking of the hydrogen bonds involving sulfur atoms and formation of new hydrogen bonds with oxygen and/or nitrogen atoms. Hence model systems, which mimic the real examples, were designed to carry out ab initio studies and to investigate the short hydrogen bonds involving sulfur atoms. Another study on the protein‐water interaction, which does not fall under the realm of the main objective of the thesis, is discussed in Chapter 8. Protein–water interaction is crucial for accomplishing many biological functions of proteins. In the recent past, natural probe tryptophan, located at the protein surfaces, has been extensively investigated using femtosecond spectroscopy experiments to understand salvation dynamics (Peonetal.,2002). In this chapter a method is described to follow up the molecular events of the protein–water interactions in detail. Tryptophan–water interaction in the protein Monellin is investigated in order to get the atomic level insights into the hydration dynamics, by carrying out MD simulations on Monellin (Hansiaetal.,2006b). The results are compared with those obtained from femtosecond resolved fluorescence spectroscopy. The time constants of the survival correlation function match well with the reported experimental values.This validates the procedure, adapted here for Monellin, to investigate the hydration dynamics in general. The last chapter (Chapter9) summarizes the results obtained from various studies and discusses the future directions. First part of this thesis aims to present the analysis by carrying out MD simulations on monomeric and dimeric TrpRS protein in order to understand the two steps of the aminoacylation reaction: activation of the aminoacid Trp in the first step and the transfer of the activated amino acid in the next step. In the second part, quantitative estimation of the geometry and the strength of pyrophosphate bond and short hydrogen bonds in proteins are reported in detail by subjecting the systems to high levels of quantum mechanical calculations(QM). The use of ab initio QM/MM calculations by combining the quantum mechanics(QM) with the molecular mechanics(MM) in order to study the enzymatic reactions is discussed as the future

In bacteria, the RNA polymerase (RNAP) consists of the following subunits: α2, β, β’, ω and σ. The core RNAP (α2ββ’ω) possesses the polymerising activity and it associates with one of the sigma factors to initiate transcription from a promoter region on the DNA template. All bacteria carry an essential housekeeping sigma factor and a number of extra cytoplasmic function (ECF) sigma factors. During alternate physiological states, a major part of transcriptional regulation is carried out by sigma factors, which act as transcriptional switches, thus, making it possible for bacteria to adapt to varied environmental signals by transcribing the necessary set of genes. Bacteriophages utilise various mechanisms for subverting the bacterial biochemical machinery for their advantage. One such example in E. coli is AsiA protein encoded by an early gene of T4 bacteriophage. Because of its property of binding to σ70, AsiA can inhibit transcription from E. coli promoters bearing –10 and –35 DNA sequences leading to inhibition of growth. σ70 of E. coli is also regulated by a stationary phase specific protein, Rsd, whose major function seems to be helping the cell in switching the transcription in favour of stationary phase genes. In this study we have investigated the mechanism of interaction of T4 AsiA and E. coli Rsd to σ70 of E. coli and also tried to determine the basis of differential inhibition of E. coli growth by AsiA and Rsd. In chapter one we have reviewed the published literature on regulation of transcription in bacteria. Some of the well known mechanisms of regulating gene expression are: DNA supercoiling, two component signal transduction system (TCS), regulation by alarmone ppGpp and 6S RNA, and sigma-antisigma interactions. Most bacteria carry a number of sigma factors and each of them is dedicated to transcribing genes in response to environmental signals. Intracellular levels of sigma factors and their binding affinity to core RNAP are deciding factors for initiating transcription from specific subsets of genes. In addition, sigma factor activity is also controlled by specific proteins, which bind to sigma factors (anti-sigma factors) under certain environmental conditions. A number of anti-sigma factors have been isolated from a variety of bacteria and the mechanisms of action of binding to cognate sigma factors have been worked out by using genetic, biochemical and structural tools. In chapter two, using yeast two hybrid assay (YTH), we have identified the regions of σ70 which interact with AsiA, and it was observed that amino acid residues from 547-603, encompassing region 4.1 and 4.2 are involved in binding to σ70. Interestingly, we found that truncated σ70 fragments lacking the N-terminal regions, apparently bound to AsiA with higher affinity compared to full length σ70. As AsiA expression, because of its transcription inhibitory activity, is inhibitory to E.coli growth, co-expression of the truncated C-terminal σ70 fragments (e.g. residues 493-613, σ70C121), which bind to σ70 with high affinity, could relieve growth inhibition. The complex of GST:AsiA-σ70C121 could be purified from E. coli cells. GST:AsiA purified from E .coli cells was found to be associated with RNAP subunits. Since further studies on this interaction required GST:AsiA preparation devoid of RNAP subunits, we decided to express this protein in S. cerevisiae. Bioinformatics analysis indicated the absence of a σ70 homologue in S.cerevisiae. As expected, GST:AsiA purified from the yeast was found to be free from any RNAP like proteins. The protein purified from yeast was used for in-vitro binding experiments. Our YTH analysis had indicated that deletion a part of region 4.1 or 4.2 of σ70 leads to loss of binding to AsiA. However, the published NMR structure of AsiA in complex with peptides corresponding to region 4 of σ70, showed that either region 4.1 or 4.2 alone can bind to AsiA indicating at the possible existence of two binding sites for AsiA. In order to confirm the physiological significance of this finding, we studied the interaction of truncated σ70 fragments lacking either region 4.1 or 4.2 with AsiA in-vivo in E. coli and in-vitro by affinity pull down assays. It was observed that σ70 fragments lacking either region 4.1 (σ70∆4.1) or 4.2 (σ70∆4.2), did not neutralize the GST:AsiA toxicity, indicating lack of interaction. The affinity purified GST:AsiA from these E. coli cells did not have σ70∆4.1 or σ70∆4.2 associated with it. Similar results were obtained from pull down assays in-vitro, where we found that σ70∆4.1 or σ70∆4.2 do not show any observable interaction with AsiA. This clearly established that the minimum region of σ70 required for physiologically relevant interaction with AsiA consists of both the regions 4.1 and 4.2. Chapter 3 of this thesis has been devoted to this aspect of AsiA-σ70 interaction. Having defined the minimum region of σ70 interacting with AsiA, we sought to identify the regions and amino acid residues of AsiA, which are critical for interaction with σ70. The approach for identification of mutants and their characterisation has been discussed in chapter 4. For this purpose, we made systematic deletions in the N and C-terminal regions of the protein and also isolated random mutants of AsiA, which lack binding to σ70 and thus are non-inhibitory to E. coli growth. It was found that deletion of 5 amino acids from N-terminus and 17 amino acids from C-terminus did not alter the inhibitory activity of AsiA. In contrast, deletion of N-terminal 10 amino residues led to complete loss of activity, while in the C-terminus, a gradual loss of activity was observed when amino acid residues beyond 17 amino acids were deleted. A 34 amino acids C-terminal deletion mutant was found to be completely inactive. E10K mutant was found to be inactive, but changes of E to other amino acids such as S, Y, L, A and Q were tolerated, indicating that negative charge at E10 is not a crucial element for interaction with σ70. Inactive mutants could be overexpressed in E. coli and showed reduced binding in YTH assay and were also poor inhibitors of in-vivo transcription in E. coli. We concluded that the primary σ70 binding site of AsiA is present in the N-terminus, yet C-terminal 64-73 amino acid residues are required for effective binding in-vivo. These studies also correlate the inhibitory potential of AsiA with its σ70 binding proficiency. In chapter 5, we have made a comparative analysis of mechanism of interaction of AsiA and Rsd to E. coli RNAP. Overexpression of Rsd was found to be less inhibitory to E. coli cell growth than that of AsiA. The affinity purified GST-AsiA from E. coli was found to have all the RNAP subunits associated with it, whereas, only σ70 was found to be associated with similarly purified GST:Rsd, pointing towards differences in binding to RNAP. In affinity pull down assays, in-vitro, it was found that both AsiA and Rsd do not show any observable binding to core RNAP. Binding of AsiA to σ70 in holo RNAP led to the formation of a ternary complex, whereas no ternary complex was observed when Rsd was made to interact with holo RNAP. Analysis of protein-protein interaction by YTH showed that region 4.1 and 4.2 are critical for binding of both AsiA and Rsd to σ70. However, in the case of Rsd, the surface of interaction is not limited to this region only and other regions of σ70 make significant contribution to this binding. Possibly, the interaction of Rsd with the core binding regions of σ70 prevents its association with core RNAP. Kinetic analysis of binding by surface plasmon resonance (SPR) showed that binding affinities (Kd) of AsiA and Rsd to σ70 are in similar range. Therefore, we concluded that the ability of AsiA to trap the holo RNAP is, probably, responsible for higher inhibitory activity of this protein compared to that of Rsd. Thus, T4 AsiA and E. coli Rsd, which share regions of interaction on σ70, have evolved differences in their mechanism of binding to RNAP such that T4 AsiA, by trapping the holo RNAP subverts the complete bacterial transcription machinery to transcribe its own genes. Rsd, on the other hand, has evolved to interact primarily with σ70, which favours the utilisation of core RNAP by other sigma factors.

A signal transduction process refers to chain of highly regulated biochemical steps which results in the transfer of signal in response to a stimulus in the extracellular environment to the intracellular compartments such as nucleus. Variety of biomolecules such as proteins and lipids participate in such processes. One of the superfamilies of proteins which actively participate in signaling processes is protein kinase which transfers γ-phosphate from Adenosine Triphosphate (ATP) to the specific hydroxyl group(s) in the protein substrates. Phosphorylation and dephosphorylation events are critical in many signal transduction pathways affecting biological system as a whole. Protein phosphorylation carried out by protein kinases has emerged as pre-eminent mechanism for the regulation of variety of cellular processes such as cell growth, development, differentiation, homeostasis, apoptosis, metabolism, transcription and translation. The current thesis encompasses the investigations carried out by the author, using various bioinformatics tools and methods, to comprehend the structural and functional roles of diverse set of protein kinase subfamilies in various eukaryotic and prokaryotic organisms. The present thesis has been divided into various chapters. Chapter 1 of the thesis provides introduction to the superfamily of protein kinases and covers the relevant literature. The database of Kinases in Genomes (KinG) set-up in the author’s group a few years ago (Krupa et al, 2004a), comprises of a collection of Serine/Threonine/Tyrosine protein kinases recognized using bioinformatics approaches, from the genomic data of various eukaryotes, prokaryotes and viruses (Krupa et al, 2004a). KinG database also provides classification of protein kinases into various groups and subfamilies (Hanks et al, 1988). Information on non-kinase domains which are tethered to the catalytic kinase domains is also available for every kinase in the KinG database. KinG is periodically (annually) updated with rise in the number of genome sequence datasets of various organisms, increase in the number of known protein domain families and refinement or reannotation of genomic datasets (Anamika et al, 2008c). Author describes the work on annual update of KinG database in Chapter 2 of the thesis. Availability of an improved version of the human genomic data has provided an opportunity to re-investigate protein kinase complement of the human genome and enabled an analysis of the splice variants. This analysis is also described in Chapter 2. Chapter 3, Chapter 4 and Chapter 5 report recognition and analysis of repertoire of protein kinases in Chimpanzee, two Plasmodium species (Plasmodium falciparum and Plasmodium yoelii yoelli) and Entamoeba histolytica respectively. A detailed analysis of the non-kinase domains which are tethered to catalytic protein kinase domains in eukaryotic organisms is presented in Chapter 6. Chapter 7 discusses a systematic classification framework developed by the author to classify Serine/Threonine protein kinases in prokaryotic organisms. Investigation carried out on 3-D structural aspects of protein kinase-substrate interactions is described in Chapter 8. While identifying protein kinases from genomic data occurrence of protein kinase-like non-kinases (PKLNK), which lack aspartate in a specific position in the amino acid sequence (and hence are unlikely to function as a kinase), has also been observed. Chapter 9 presents an analysis of PKLNKs with an objective of obtaining clues to their functions. Chapter 10 summarizes the main conclusions of the thesis and provides an outlook of the current study. Chapter 1: Chapter 1 provides an introduction to cell signaling and the involvement of protein kinases in various signaling pathways compiled from author’s literature survey. This chapter provides a description of molecular events in cell signaling in prokaryotic and eukaryotic organisms. The diversity, specificity and cellular roles of protein kinases are discussed in detail. Chapter 2: Chapter 2 describes KinG (Kinases in Genomes) database which was first established by Krupa et al (2004a). The KinG database is an on-line compilation of the putative Serine/Threonine/Tyrosine protein kinases encoded in the completely sequenced genomes of archaea, eubacteria, viruses and eukaryotes. Surge in the datasets of genomes, improvements in the quality of the genomic data for various organisms and growing number of protein domain families necessitates periodic update of KinG database. The updated version of KinG holds information on protein kinases for 483 organisms (Anamika et al, 2008c). Availability of draft version of the human genome data in 2001 enabled recognition of repertoire of human protein kinases (Krupa and Srinivasan, 2002a; Manning et al, 2002; Kostich et al, 2002). Over the last 7 years human genomic data is being refined and at present the quality of the human genomic data available is much superior to the one available in 2001. By gleaning the latest version of human genome data, 46 new human protein kinase splice variants have been identified which were not recognized in the earlier studies on human kinome. Improper regulation or mutant forms of many of these newly identified protein kinase splice variants are directly involved in various diseases such as different kinds of cancer, Severe Combined Immunodeficiency Disease (SCID) and Huntington disease. In addition, abnormal forms of mouse orthologues of some of the newly identified human kinase splice variants are known to cause various diseases in mice. This raises the possibility of the human orthologues playing similar roles in the disease processes. Such observations and detailed analysis of these protein kinase splice variants would have a profound influence on drug design and development against various diseases. Chapter 3: Investigations on the identification and analysis of protein kinases encoded in the genome of chimpanzee (chimp) has been discussed in Chapter 3. Further, the kinome complement has been compared between chimp and its evolutionary close relative, human (Anamika et al, 2008b). The shared core biology between chimp and human is characterized by many orthologous protein kinases which are involved in conserved pathways. Domain architectures specific to chimp/human kinases have been observed. Chimp kinases with unique domain architectures are characterized by deletion of one or more non-kinase domains present in the human kinases. Interestingly, counterparts of some of the multi-domain human kinases in chimp are characterized by identical domain architectures but with kinase-like non-kinase domain (PKLNK). Remarkably, for 160 out of 587 chimpanzee kinases no human orthologue with sequence identity greater than 95% could be identified. Variations in chimpanzee kinases compared to human kinases are brought about also by differences in functions of domains tethered to the catalytic kinase domain. For example, the heterodimer forming PB1 domain related to the fold of ubiquitin / Ras-binding domain is seen uniquely tethered to PKC-like chimpanzee kinase. Though chimpanzee and human have close evolutionary relationship, there are chimpanzee kinases with no close counterpart in the human suggesting differences in their functions. This chapter provides a direction for experimental analysis of human and chimpanzee protein kinases in order to enhance our understanding on their specific biological roles. Chapter 4: Chapter 4 describes genome-wide comparative analysis for protein kinases encoded in the two apicomplexa namely Plasmodium falciparum (P. falciparum) (3D7 strain) and Plasmodium yoelii yoelii (P. yoelii yoelii) (17XNL strain) genomes which are causative agents of malaria in human and rodent respectively (Anamika and Srinivasan, 2007). Sensitive bioinformatics techniques enable identification of 82 and 60 putative protein kinases in P. falciparum and P. yoelii yoelii respectively. These protein kinases have been classified further into subfamilies based on the extent of sequence similarity of their catalytic domains (Hanks et al, 1988). The most populated kinase subfamilies in both the Plasmodium species correspond to CAMK and CMGC groups. Analysis of domain architectures enables detection of uncommon domain organisation in kinases of both the organisms such as kinase domain tethered to EF hands as well as pleckstrin homology domain. Components of MAPK signaling pathway are not well conserved in Plasmodium species. Such observations suggest that Plasmodium protein kinases are highly divergent from other eukaryotes. A trans-membrane kinase with 6 membrane spanning segments in P. falciparum seems to have no orthologue in P. yoelii yoelii. 19 P. falciparum kinases (Anamika et al, 2005; Anamika and Srinivasan, 2007) have been found to cluster separately from P. yoelii yoelii kinases and hence these kinases are unique to P. falciparum genome. Only 28 orthologous pairs of kinases could be identified between these two Plasmodium species. Comparative kinome analysis of the two Plasmodium species has thus provided clues to the function of many protein kinases based upon their classification and domain organisation and also implicate marked differences even between two Plasmodium species. Chapter 5: Identification and analysis of the repertoire of protein kinases in the intracellular parasite Entamoeba histolytica (E. histolytica) using sensitive sequence and profile search methods forms the basis of Chapter 5. A systematic analysis of a set of 307 protein kinases in E. histolytica genome has been carried out by classifying them into different subfamilies originally defined by Hanks and Hunter (Hanks et al, 1988) and by examining the functional domains which are tethered to the catalytic kinase domains (Anamika et al, 2008a). Compared to other eukaryotic organisms, protein kinases from E. histolytica vary in terms of their domain organisation and displays features that may have a bearing in the unusual biology of this organism. Some of the parasitic kinases show high sequence similarity in the catalytic domain region with calmodulin/calcium dependent protein kinase subfamily. However, they are unlikely to act like calcium/calmodulin dependent kinases as they lack non-catalytic domains characteristic of such kinases in other organisms. Such kinases form the largest subfamily of protein kinases in E. histolytica. Interestingly a Protein Kinase A/Protein Kinase G-like hybrid kinase subfamily member is tethered to pleckstrin homology domain. Although potential cyclins and cyclin-dependent kinases could be identified in the genome the likely absence of other cell cycle proteins suggests unusual nature of cell cycle in E. histolytica. Some of the unusual kinases recognized in the analysis include the absence of Mitogen activated protein kinase kinase (MEK) as a part of the Mitogen Activated Kinase signaling pathway and identification of trans-membraneous kinases with catalytic kinase region showing a good sequence similarity to Src kinases which are usually cytosolic. Sequences which could not be classified into known subfamilies of protein kinases have unusual domain architectures. Many such unclassified protein kinases are tethered to domains which are cysteine-rich and to domains known to be involved in protein-protein interactions. The current chapter on kinome analysis of E. histolytica suggests that the organism possesses a complex protein phosphorylation network that involves many unusual protein kinases. Chapter 6: Protein kinases phosphorylating Serine/Threonine/Tyrosine residues in several cellular proteins exert tight control over their biological functions. They constitute the largest protein family in most eukaryotic species. Classification based on sequence similarity of their catalytic domains, results in clustering of protein kinases sharing gross functional properties into various subfamilies. Many protein kinases are associated or tethered covalently to domains that serve as adapter or regulatory modules, aiding substrate recruitment, specificity, and also serve as scaffolds. Hence the modular organisation of the protein kinases serves as guidelines to their molecular interaction which has been discussed in Chapter 6. Recent studies on repertoires of protein kinases in eukaryotes have revealed wide spectrum of domain organisation in model organisms across various subfamilies. Occurrence of organism-specific novel domain combinations suggests functional diversity achieved by the protein kinase in order to regulate variety of biological processes. In addition, domain architectures of protein kinases revealed existence of hybrid protein kinase subfamilies and their emerging roles in the signaling of eukaryotic organisms. The repertoire of non-kinase domains tethered to multi-domain kinases in the higher eukaryotes is discussed in Chapter 6. Similarities and differences in the domain architectures of protein kinases in these organisms indicate conserved and unique features that are critical to functional specialization. Chapter 7: Chapter 7 describes systematic classification of Serine/Threonine protein kinases encoded in archaeal and eubacterial genomes. Majority of the Serine/Threonine protein kinases which have been identified in archaeal and eubacterial genomes could not be classified into any of the well known subfamilies (Hanks et al, 1988) of protein kinases suggesting their diversity from kinases in eukaryotes. The extensive prokaryotic Serine/Threonine protein kinase dataset obtained from KinG (Krupa et al, 2004a, Anamika et al, 2008c) has given an opportunity to classify these prokaryotic Serine/Threonine protein kinases mainly into three categories based upon sequence identity based clustering: 1) Species/Genus-specific clusters: Species/Genus-specific Serine/Threonine protein kinases contain members from a particular species or genus of the eubacteria or archaea suggesting requirement of these Serine/Threonine protein kinases for certain lineage specific function. 2) Organism-specific clusters: Organism specific clusters has members from certain specific types of organisms which suggests role of these Serine/Threonine protein kinases in some specific function being carried out by limited sets of prokaryotes. 3) Organism-diverse clusters: Organism diverse clusters suggest common function performed by such kinases in wide variety of organisms. Interestingly, occurrence of several species/genus or organism specific subfamilies of prokaryotic kinases contrasts with classification of eukaryotic protein kinases in which most of the popular subfamilies of eukaryotic protein kinases occur diversely in several eukaryotes. Function-based classification has also been proposed which shows that members of each cluster has specific function to perform. In this analysis, almost 50% of the “clusters” obtained have only one member suggesting their sequence and probably functional divergence. Many prokaryotic Serine/Threonine protein kinases exhibit a wide variety of modular organisation which indicates a degree of complexity and protein-protein interactions in the signaling pathways in these microbes. Chapter 8: A wide spectrum of protein kinases belonging to different Hanks and Hunter groups of kinases and subfamilies has been identified in various eukaryotes. However, specific biological targets (substrates) of many protein kinase subfamilies are still unknown and this is one of the active areas of research. In the current analysis reported in Chapter 8, an attempt has been made to understand protein kinase-substrate interaction and substrate consensus prediction by analyzing known 3-D structures of complexes of kinases and peptide substrates/pseudosubstrates. Considering protein kinase ternary complex structures in their active states, it has been observed that protein kinase residues which are interacting with the substrate residues having constraint are at topologically equivalent positions despite belonging to different Hanks and Hunter protein kinase subfamilies. In this analysis, it has also been observed that the residues in a given kinase subfamily interacting with consensus substrate residues are usually conserved across homologues. Interestingly, in Protein Kinase B and Phosphorylase Kinase subfamily homologues, residues interacting with substrate residue/s having no constraint are not well conserved even within the kinase subfamily suggesting different evolutionary rate of substrate interacting residues. This result is anticipated to be helpful in furthering our understanding of protein kinase-substrate relationship which is likely to be helpful in drug design. Chapter 9: Protein Kinase-Like Non-kinases (PKLNKs) are closely related to protein kinases but they lack the crucial catalytic aspartate in the catalytic loop and hence cannot function as a protein kinase. PKLNKs have been analyzed (Chapter 9) with an objective of obtaining clues about their functions. Using various sensitive sequence analysis methods, 82 PKLNKs from four higher eukaryotic organisms namely, Homo sapiens, Mus Musculus, Rattus norvegicus and Drosophila melanogaster have been recognized. On the basis of their domain combinations and functions of tethered domains, PKLNKs have been classified mainly into four categories: 1) Ligand binding PKLNKs 2) PKLNKs having extracellular protein-protein interaction domain 3) PKLNKs involved in dimerization 4) PKLNKs with cytoplasmic protein-protein interaction module. While members of the first two classes of PKLNKs have transmembrane domain tethered to the PKLNK domain, members of the other two classes of PKLNKs are entirely cytoplasmic in nature. The current classification scheme hopes to provide a convenient framework to classify the PKLNKs from other eukaryotes and it should be helpful in deciphering their roles in cellular processes. Chapter 10: This is a chapter on conclusions of the entire thesis work. Summary of the major outcomes of this thesis work is provided and implications of the work in the area of signal transduction are discussed. In addition to above mentioned work, studies on repertoire of protein kinases from two plant organisms have been carried out and the kinomes have been comparatively analyzed (Krupa et al, 2006) (Appendix 1). Comparison of plant protein kinases with other eukaryotes revealed remarkable differences. Trans-genomic comparison of the protein kinase repertoires of Arabidopsis thaliana and Oryza sativa has enabled identification of members that are uniquely conserved within the two species. Analysis on the domain organisation of plant protein kinases has also been carried out. Appendix 2 presents the work done on Entamoeba histolytica (E. histolytica) ornithine decarboxylase (ODC)-like protein which regulates the polyamine biosynthesis. DFMO (Difluoromethylornithine) is unable to inhibit the E. histolytica ODC-like protein while it inhibits the homologues of ODC in other organisms. Modelling study has suggested substitution of three amino acids in the E. histolytica ODC-like protein because of which DFMO is unable to inhibit the activity of ODC-like protein (Jhingran et al, 2008). All the computational modeling work reported in Appendix 2 was performed by the author while all the laboratory experiments were performed in the laboratory of the collaborator Prof. Madhubala of JNU, New Delhi. The supplementary data pertaining to this thesis is presented in an accompanying CD. The supplementary data in this CD is organized into different folders corresponding to various chapters.

Tumor suppressor p53, the guardian of the genome, has been intensely studied molecule owing to its central role in maintaining cellular integrity. While the level of p53 protein is maintained low in unstressed conditions, there is a rapid increase in the functional p53 protein levels during stress conditions. It is now well documented in literature that p53 protein accumulates in the cells following DNA damage by posttranslational modifications leading to increased stability and half life of protein. Additionally, recent studies have also highlighted the significance of increased p53 translation during stress conditions. Interestingly, an alternative initiation codon has been shown to be present within the coding region of p53 mRNA. Translation initiation from this internal AUG results in an N-terminally truncated p53 isoform, described as ΔN-p53. However, the mechanisms underlying co-translational regulation of p53 and ΔN-p53 are still poorly understood. Studies have suggested that synthesis of both p53 and its ΔN-p53 isoform is regulated during cell cycle and also stress and cell-type specific manner. Interestingly, reports also demonstrate continued synthesis of both p53 isoforms during stress conditions. In contrast, global rates of cap-dependent translation initiation are shown to be reduced during stress conditions. This translation attenuation is observed mainly due to restricted availability of critical initiation factors. Interestingly, preferential synthesis of a vital pool of survival factors persists even during these circumstances. Studies have suggested that this selective translation is mediated via alternative mechanisms of translation initiation. One of the important mechanisms used for protein synthesis during these conditions is internal initiation. In this mechanism, the ribosomes are recruited to a complex RNA structural element known as ‘Internal Ribosome Entry Site (IRES)’, generally present in the 5’ untranslated region (UTR) of mRNA. Therefore, it is possible that the translation of p53 and ΔN-p53 could also be regulated by IRES mediated translation, especially during stress conditions. In this thesis the role of internal initiation in translational control of p53 and ΔN-p53 has been investigated. Additionally, the putative secondary structure of p53 IRES RNA has been determined. Further, it has been shown that polypyrimidine tract binding (PTB) protein acts as an important regulator of p53 IRES activities. The probable mechanism of action of PTB protein has also been investigated. The results suggest that interaction with PTB alters the p53 IRES conformation which could facilitate translation initiation. Finally, the possible physiological significance of existence of p53 IRES elements has been addressed. In the first part of the thesis, the presence of internal ribosome entry site within p53 mRNA has been investigated. As a first step, the 5’UTRs mediating the translation of both p53 and ΔN-p53 were cloned in the intercistronic regions of bicistronic constructs. Results of in vivo transfection of these bicistronic constructs suggested the presence of two IRES elements within p53 mRNA, with activities comparable to known viral and cellular IRESs. The IRES directing the translation of p53 is in the 5'-untranslated region of the mRNA, whereas the IRES mediating the translation of ΔN-p53 extends further into the protein-coding region. To further validate, stringent assays were performed to rule out the possibility of any cryptic promoter activity, re-initiation/scanning or alternative splicing in the p53 mRNA. Transfection of in vitro synthesized bicistronic RNAs confirmed the presence of IRES elements within p53 mRNA. Incidentally, this constitutes the first report on translational control of p53 by internal initiation. In the second part of the thesis, the secondary structure of p53 IRES RNA has been investigated. Structural analysis of p53 RNA was performed using structure-specific nucleases and modifying chemicals. The results obtained from chemical modification and nuclease probing experiments were used to constrain Mfold predicted structures. Based on this, a putative secondary structure model for p53 IRES RNA has been derived. Sequence alignment suggested that the p53 IRES RNA showed significant sequence conservation across mammalian species. To study the effect of mutations on the IRES structure, mutant p53 IRESs were used that harbor silent mutations at critical locations within the p53 IRES element. Incidentally, one of the mutant constructs used in the study was observed to be a naturally occurring mutation in a chronic lymphocyte leukemia patient. RNA structure analyses of these two mutant p53 IRES RNAs were performed. The nuclease mapping data suggested conformational alteration in these mutant RNAs with respect to wild type. Consistently, a comparative Circular-Dichroism spectroscopy of the Wt and mutant RNAs also validated the conformational alteration of the mutant RNAs. This also suggested that the presence of mutations in p53 IRES might result in decreased induction of p53 protein following DNA damage due to altered RNA structure. This might constitute as one of the mechanisms leading to tumor development in some types of cancers. In the third part of the thesis, the role of important cellular proteins that might modulate p53 IRES mediated translation has been studied. These cellular proteins act as IRES interacting trans-acting factors (ITAFs). Polypyrimidine tract binding (PTB) protein is an important ITAF implicated in regulating IRES mediated gene expression during apoptosis. It was observed that PTB protein specifically interacts with both the IRES elements within p53 mRNA. Interestingly, the affinity of interaction of PTB protein with both p53 IRES RNAs was observed to be significantly different. In order to determine the contact points of PTB on p53 IRES, a foot-printing assay using structure specific nuclease and recombinant-PTB protein was performed on p53 RNA. The data from foot-printing as well as primer extension inhibition assay (toe-printing analysis) suggested the presence of multiple PTB binding sites on p53 IRES RNA. Based on these results, a deletion mutant was generated that showed reduced PTB binding and also reduced IRES activity as compared to wild type. Further, to study the role of PTB in mediating p53 translation, the expression of PTB gene was partially silenced by using PTB specific siRNA. Partial depletion of endogenous PTB protein showed a significant decrease in the p53 IRES activities. These results suggest that PTB protein is essential for the p53 IRES activities. To understand the probable mechanism by which PTB regulates p53 IRES mediated translation, CD spectroscopy analysis of p53 IRES RNA was performed in the absence and presence of PTB protein. Interestingly, CD spectra analysis of the p53 RNA in the presence of PTB suggested a specific conformational change in p53 IRES, which might probably facilitate ribosome loading during internal initiation. This also suggests that abnormal expression of p53 ITAFs might lead to reduced p53 induction following DNA damage conditions. It could also be another event leading to malignant transformation of cells bearing wild type p53. It is highly tempting to speculate that the levels of p53 ITAFs could also be used as tumor biomarkers. In the fourth part of the thesis, the physiological relevance of existence of IRES elements within p53 mRNA has been investigated. The levels of p53 and ΔN-p53 proteins are known to be regulated in a cell cycle phase-dependent manner. The IRES activities of both p53 IRES elements were investigated at different phases of cell cycle. The activity of the IRES responsible for translation of p53 protein was found to be highest at G2-M transition and the maximum IRES activity corresponding to ΔN-p53 synthesis was observed at G1-S transition. These results suggested that the p53 IRES activities are regulated in a cell-cycle phase-dependent manner. Next, the regulation of p53 IRES mediated translation during stress conditions was studied. Human lung carcinoma cell line, A549 cells (that endogenously express both the p53 isoforms), were exposed to DNA damaging drug, doxorubicin. The level of p53 protein was observed to increase in a time-dependent manner. Interestingly, PTB protein, which is predominantly nuclear, was found to translocate to the cytoplasm during stress condition in a time-dependent manner. Under similar conditions, p53 protein was observed to reverse translocate from the cytoplasm to nucleus, probably to function as a transcription factor. Next, the influence of partial PTB silencing on p53 isoforms in the presence of cell stress (mediated by doxorubicin) was investigated. The data indicated reduced levels of both p53 and ΔN-p53 when PTB gene expression was partially silenced. These observations constitute “the proof of concept” that relative abundance of an ITAF, such as PTB protein, might contribute to regulating the coordinated expression of the p53 isoforms. The thesis reveals the presence as well as the physiological relevance of existence of IRES elements within p53 mRNA. The novel discovery of p53 IRES elements may provide new insights into the underlying mechanism of translational regulation. The modulation of the p53 IRES activities by PTB protein suggests that the regulated expression of p53 isoforms depends on the integrity of IRES elements and availability of cellular proteins that can serve as p53 ITAFs. Thus, studies pertaining to the identification of mutations within p53 IRES region as well as abnormal expression of p53 ITAFs such as PTB in cancer cells may have far reaching implications. These studies might lead to further advances in the field of cancer detection, prognosis and design of novel therapeutic strategies.