<p>This thesis addresses different aspects of the question about accuracy of protein synthesis: i) the mechanism of tRNA selection during translation ii) study of ribosomal mutations that affect accuracy and iii) the choice of aminoacyl-tRNA isoacceptors on synonymous codons.</p><p>By measuring the codon reading efficiencies of cognate and near-cognate ternary complexes we demonstrate that in optimal physiological conditions accuracy of substrate selection is much higher than previously reported; that during translation the ribosomal A site is not blocked by unspecific binding of the non-cognate tRNAs which would inhibit the speed of protein synthesis. Our results suggest that there is an asymmetry between initial selection and proofreading step concerning the wobble position, and that binding of non-cognate substrate does not induce GTP hydrolysis on the ribosome.</p><p>The knowledge obtained from the ribosomal mutant strains can be used to explain the general relation between the structure of the ribosome and the mechanism of codon recognition, as well as the streptomycin resistance or dependence phenomenon.</p><p>Our work showed experimentally that the probability for binding certain tRNA to the A site of the ribosome is not based on the simple codon-anticodon base pair matching. In the living cell the availability of cognate tRNAs versus the demand for them (the frequency of codon usage) is finely balanced to ensure critical protein synthesis in stress conditions. We have also discovered a new codon assignment for a specific tRNALeu isoacceptor and detected a base modification in its anticodon, which has not been previously observed. The motivation for the later findings comes from a system biology modeling and the results are an example of an interdisciplinary collaboration.</p>

<p>RNA is the only biological molecule that can function both as a repository of information and as a catalyst. This, together with the ability to self-replicate, led to recognition of RNA as ‘prelude to life’.</p><p>My work highlights some of the important features of RNA as a catalyst, exemplified by RNase P. It addresses questions of evolutionary preservations of residues and structure, involvement of metal ions and finally structure evolution towards minimal catalytically competent RNA motifs.</p><p>RNase P is the only enzyme involved in 5’ end processing of all pre-tRNAs. Until recently, it was believed that the RNA moiety of RNase P is responsible for mediating catalysis only in Bacteria. However, my recent study conclusively demonstrated that eukaryotic RNase P RNA is catalytically competent in vitro in absence of proteins. These findings evidenced evolutionary preservation of RNA-mediated catalysis in RNase P.</p><p>RNase P RNA is a metalloeznyme. In my studies I analyzed the contributions of individual chemical groups at the cleavage site to catalysis. My findings suggested that the 2’OH of N_-1 and the exocyclic amine of G₊₁ are involved in positioning of functionally important metal ions. Additionally, data appointed the function of Pb²⁺ as both structural metal ion and important in generating the nucleophile. My studies further indicate a conformational change upon RNase P RNA -substrate complex formation in keeping with an induced fit mechanism. </p><p>Studying the effects of reducing the ribozyme size upon dissection of bacterial RNase P RNAs, we defined the smallest catalytically competent domain i.e. P15-loop. Derivatives of this autonomous metal ion binding domain, (the smallest being 31nt-s), are able to cleave both whole-length pre-tRNAs as well as hairpin substrates, though with severely reduced rates relative to their parent ribozymes. The study has inferred that partite ES interactions at the cleavage site prove sufficient for catalysis.</p>

The stringent response plays a significant role in the survival of bacteria during different environmental conditions. It is activated by the binding of stringent factor (SF) to stalled ribosomes that have an unacylated tRNA in the ribosomal A-site which leads to the synthesis of (p)ppGpp. ppGpp binds to the RNA polymerase, resulting in a rapid down-regulation of rRNA and tRNA transcription and up-regulation of mRNAs coding for enzymes involved in amino acid biosynthesis. The importance of the A-site and unacylated tRNA in the activation of SF was confirmed by chemical modification and subsequent primer extension experiments (footprinting experiments) which showed that binding of SF to ribosomes resulted in the protection of regions in 23S rRNA, the A-loop and helix 89 that are involved in the binding of the A-site tRNA. An in vitro assay showed that the ribosomal protein L11 and its flexible N-terminal part was important in the activation of SF. Interestingly the N-terminal part of L11 was shown to activate SF on its own and this activation was dependent on both ribosomes and an unacylated tRNA in the A-site. The N-terminal part of L11 was suggested to mediate an interaction between ribosome-bound SF and the unacylated tRNA in the A-site or interact with SF and the unacylated tRNA independently of each other. Footprinting experiments showed that SF bound to the ribosome protected bases in the L11 binding domain of the ribosome that were not involved in an interaction with ribosomal protein L11. The sarcin/ricin loop, in close contact with the L11 binding domain on the ribosome and essential for the binding and activation of translation elongation factors was also found to be protected by the binding of SF. Altogether the presented results suggest that SF binds to the factor-binding stalk of the ribosome and that activation of SF is dependent on the flexible N-terminal domain of L11 and an interaction of SF with the unacylated tRNA in the A-site of the 50S subunit.

RNA is the only biological molecule that can function both as a repository of information and as a catalyst. This, together with the ability to self-replicate, led to recognition of RNA as ‘prelude to life’. My work highlights some of the important features of RNA as a catalyst, exemplified by RNase P. It addresses questions of evolutionary preservations of residues and structure, involvement of metal ions and finally structure evolution towards minimal catalytically competent RNA motifs. RNase P is the only enzyme involved in 5’ end processing of all pre-tRNAs. Until recently, it was believed that the RNA moiety of RNase P is responsible for mediating catalysis only in Bacteria. However, my recent study conclusively demonstrated that eukaryotic RNase P RNA is catalytically competent in vitro in absence of proteins. These findings evidenced evolutionary preservation of RNA-mediated catalysis in RNase P. RNase P RNA is a metalloeznyme. In my studies I analyzed the contributions of individual chemical groups at the cleavage site to catalysis. My findings suggested that the 2’OH of N-1 and the exocyclic amine of G+1 are involved in positioning of functionally important metal ions. Additionally, data appointed the function of Pb2+ as both structural metal ion and important in generating the nucleophile. My studies further indicate a conformational change upon RNase P RNA -substrate complex formation in keeping with an induced fit mechanism. Studying the effects of reducing the ribozyme size upon dissection of bacterial RNase P RNAs, we defined the smallest catalytically competent domain i.e. P15-loop. Derivatives of this autonomous metal ion binding domain, (the smallest being 31nt-s), are able to cleave both whole-length pre-tRNAs as well as hairpin substrates, though with severely reduced rates relative to their parent ribozymes. The study has inferred that partite ES interactions at the cleavage site prove sufficient for catalysis.

Transfer RNA from all organisms has modified nucleosides and position 34 (the wobble position) is one of the most extensively modified positions. Some wobble nucleoside modifications restrict codon choice (e.g. 5-methylaminomethyl-2-thiouridine, mnm5s2U) while some extend the decoding capacity (e.g. uridine-5-oxyacetic acid, cmo5U). In this thesis the influence of wobble nucleoside modification on cell physiology and translation efficiency and accuracy is described. A mutant proL tRNA (proL207) was isolated that had an unmodified adenosine in the wobble position. Surprisingly, the proL207 mutant grows normally and is efficiently selected at the non-complementary CCC codon. The explanation of how an A34 containing tRNA can read CCC codon could be that a protonated A can form a base pair with C. cmo5U (uridine-5-oxyacetic acid) is present in the wobble position of five tRNA species in S.enterica. Two genes (cmoA and cmoB) have been identified that are involved in the synthetic pathway of cmo5U. Mutants were constructed in alanine, valine, proline, and threonine codon boxes which left only a cmo5U containing tRNA present in the cell. The influence of cmo5U on growth or on A site selection rates of the ternary complex was found to be tRNA dependent. During the study of the frameshift suppressor sufY of the hisC3737 frameshift mutation, a dominant mutation was found in YbbB protein, a selenouridine synthetase. The frameshifting occurs at CCC-CAA codon contexts and is specific for CAA codons, which are read by tRNAGlncmnm5s2UUG . The sufY204 mutation is a dominant mutation resulting in a change from Gly67 to Glu67 in the YbbB protein, and mediates the synthesis of several novel modified nucleosides/nucleotides (UKs) with unknown structure. The synthesis of these UKs is connected to the synthesis of cmnm5s2U34. The presence of UK on tRNAGlnU*UG reduced aminoacylation and therefore might account for the slow entry at CAA codons which could result in +1 frameshifting by P site tRNA. The selenourdine synthetase activity is not required for the synthesis of UKs. We hypothesize that an intrinsic activity that is low in the wild type protein has been elevated by the single amino acid substitution and results in the synthesis of UKs.

This thesis addresses different aspects of the question about accuracy of protein synthesis: i) the mechanism of tRNA selection during translation ii) study of ribosomal mutations that affect accuracy and iii) the choice of aminoacyl-tRNA isoacceptors on synonymous codons. By measuring the codon reading efficiencies of cognate and near-cognate ternary complexes we demonstrate that in optimal physiological conditions accuracy of substrate selection is much higher than previously reported; that during translation the ribosomal A site is not blocked by unspecific binding of the non-cognate tRNAs which would inhibit the speed of protein synthesis. Our results suggest that there is an asymmetry between initial selection and proofreading step concerning the wobble position, and that binding of non-cognate substrate does not induce GTP hydrolysis on the ribosome. The knowledge obtained from the ribosomal mutant strains can be used to explain the general relation between the structure of the ribosome and the mechanism of codon recognition, as well as the streptomycin resistance or dependence phenomenon. Our work showed experimentally that the probability for binding certain tRNA to the A site of the ribosome is not based on the simple codon-anticodon base pair matching. In the living cell the availability of cognate tRNAs versus the demand for them (the frequency of codon usage) is finely balanced to ensure critical protein synthesis in stress conditions. We have also discovered a new codon assignment for a specific tRNALeu isoacceptor and detected a base modification in its anticodon, which has not been previously observed. The motivation for the later findings comes from a system biology modeling and the results are an example of an interdisciplinary collaboration.

tRNA molecules contain more than 80 chemically unique nucleotide base modifications that contribute to the chemical and physical diversity of RNAs as well as add to the overall fitness of the cell. For instance, base modifications have been shown to play a critical role in tRNA molecules by improving the fidelity and efficiency of translation. Most of this work has been carried out extensively in Gram-negative bacteria, however, the role of modified bases in tRNAs as they relate to thermostability, structure, and transcriptional regulation in Gram-positive bacteria, such as Bacillus subtilis and Bacillus anthracis, are not well characterized. Infections by Gram-positive bacteria that have become more resistant to established drug regiments are on the rise, making Gram-positive bacteria a serious threat to public safety. My thesis work examined what role partial base modification of the tyrosyl-anticodon stem-loops (ASLTyr ) of B. subtilis and B. anthracis have on thermostability, structure, and transcriptional regulation. The ASLTyr molecules have three modified residues which include Queuine (Q34), 2-thiomethyl-N6-dimethylallyl (ms2i6A37), and pseudouridine (Y39). Differential Scanning Calorimetry (DSC) and UV melting were employed to examine the thermodynamic effects of partial modification on ASLTyr stability. The DSC and UV data indicated that the Y39 and i6A37 modifications improved the molecular stability of the ASL. To examine the effects of partial base modification on ASLTyr structure, NMR spectroscopy was employed. The NMR data indicated that the unmodified and [Y39]-ASLTyr form a protonated C-A+ Watson-Crick-like base pair instead of the canonical bifurcated C-A+ interaction. Additionally, the loop regions of the unmodified and [Y39]-ASLTyr molecules were well ordered. Interestingly, the [i6A37]- and [i6A37; Y39]- ASLTyr molecules did not form a protonated C-A+ base pair and the bases of the loop region were not well ordered. The NMR data also suggested that the unmodified and partially modified molecules do not adopt the canonical U-turn structure. The structures of the unmodified, [Y39]-, and [i6A37;Y39]-ASLTyr molecules did not depend on the presence of Mg2+, but the structure of the [i6A37]-ASLTyr molecule did depend on the presence of multivalent cations. Finally, to determine the repercussions that partial modification has on physiology and tRNA mediated transcriptional regulation in B. anthracis, antibiotic sensitivity tests, growth curves, and quantitative real-time polymerase chain reaction (qRT-PCR) were employed. Strains deficient in ms2 showed comparable growth to the parent strain when cultured in defined media, but Q deficient strains did not. The loss of ms2i6A37 conferred resistance to spectinomycin and ciprofloxacin, whereas the loss of Q34 resulted in sensitivity to erythromycin. Changes in the ratio full-length to truncated transcripts of the tyrS1 and tyrS2 genes were used to monitor tRNA mediated transcriptional regulation. The qRT-PCR data suggested that tyrS1 and tyrS2 are T-box regulated and that the loss of ms2i6A37 and Q34 might affect the interaction of the tRNATyr molecule with the specifier sequence, which is located in the 5’-untranscribed region (UTR) of the messenger RNA (mRNA).

The studies reported in this thesis address the aspects of initiator tRNA selection in Escherichia Coli. A summary of the relevant literature discussing the process of ptotein biosynthesis in general and initiator tRNA selection, in particular is presented in chapter 1. The next chapter (Chapter2) describes the ‘Materials and Methods’ used throughout the experimental work carried out in this thesis. It is followed by two chapters(Chapter 3 and Chapter 4) which describe the isolation and characterization of an E. coli mutant, to understand the mechanism of initiator tRNA selection. Chapter 5 comprises of some experimental work and future perspectives on the utility of the E.coli mutant. The last chapter (Chapter 6) summarizes the published work where I have contributed to besides the work described in Chapters 3 to 5. The summary of chapters 3-5 is as described below:- (i)Isolation and genetic mapping of extragenic suppressors of mutant initiator tRNA lacking the three consecutive G, C base pairs in the anticodon stem Initiator tRNA selection on the ribosomes is a result of several steps, some of which are unique to the prokaryotic world. Structure-function analyses of E.Coli tRNAfMet have revealed that the most important features of tRNAfMet, pertinent to its in vivo function as an initiator, are located in the acceptor stem and the anticodon arm regions. The three consecutive G-C base pairs in the anticodon stem of the tRNAfMet, conserved across all kingdoms of life, have been implicated in preferential binding to 30S ribosomal P-site. How the 3G-C base pairs are exploited by ribosomes in selecting the initiator tRNA, has been a long standing question. In the present work, a genetic screen was developed to isolate second site compensatory mutations of the mutant tRNAfMet, inactive in initiation because the 3G-C base pairs in it were changed to those found in the elongator tRNAMet(‘3G-C mutant’). Two extragenic suppressors were mapped to defined regions in the 12 min and 85 min locations in the E. Coli genome and three others were classified in these two broad groups. A super suppressor strain exhibiting synergistic suppression was generated. Further genetic mapping identified a G122D mutation in the folD gene encoding 5, 10 methylene tetrahydrofolate dehydrogenase/cyclohydrolase in one of the suppressor strains E. Coli A48. Complementation analysis using over expression of fold confirmed the results obtained by genetic mapping. (ii) Role of the intracellular S-adenosylmethionine flux in initiation with an initiator versus elongator tRNAs in Escherichia Coli How a defect in folD gene product (in E. Coli A48) leads to initiation with the ‘3G-C mutant’ initiator tRNA, has been addressed in this work. The FolD enzyme plays a key role in the one-carbon metabolism. The mutation in folD resulted in a lethal phenotype in minimal medium. The end-products of the pathway, 10 formyl-THF, methionine and S-adenosylmethionine(SAM) were analyzed for their possible role in initiation with the ‘3G-C mutant’ tRNAfMet, which revealed that lowering of the steady-state abundance of methionine and SAM had a direct role in initiation with the ‘3G-C mutant” tRNAfMet. Analysis of the 16S tRNA revealed that the methylations, as a result of reduced levels of SAM, were undetectable in the E.Coli A48. This prompted us to generate targeted mutations in the methyltransferase genes, which have highlighted the importance of methylations in initiator tRNA selection. Consistent with the growth retardation phenotype of methylase deficient strains at higher temperatures, the E. Coli A48 also displays temperature sensitivity. Further analysis of mycoplasma genomes, which do not follow the strong conservation of three G-C base pairs in the anicodon stem of initiator tRNA has uncovered an hitherto unknown evolutionary connection between methylations of 16S rRNA and initiator tRNA selection. We observed genetic interaction between infC(encoding IF3) and fold (encoding FolD). We also demonstrate initiation with tRNAfMet containing mutations in one, two or all the three G-C base pairs, as also with the elongator tRNA (tRNAGln). (iii) Utility of E. Coli A48 in investigation of biological processes: Some Preliminary studies and future perspectives. The availability of the E. Coli A48 strain is a valuable addition to the field of initiator tRNA selection and opens up further opportunities for its application. In this study, we have analyzed some of the properties of the E. Coli A48 strain viz. sensitivity to UV light and formylation independent initiation. E. Coli possess multiple copies of initiator tRNA, encoded by the metZVW operon and the metY gene. We reasoned that the abundance of cellular initiator tRNA might be a contributing factor in maintenance of specificity of initiation. Consistent with our prediction, we observed initiation with the ‘3G-C mutant’ tRNAfMet in E. Coli strains deficient in initiator tRNA genes. The various aspects of SAM limitation, biological functions of post-transcriptional modifications, incorporation of non-methionine amino acids in then-terminus of proteins and genetic approaches to system biology for the understanding of one-carbon metabolism are discussed.

Aminoacyl-tRNA synthetases (aaRSs) are at the center of the question of the origin of life and are essential proteins found in all living organisms. AARSs arose early in evolution to interpret genetic code and are believed to be a group of ancient proteins. They constitute a family of enzymes integrating the two levels of cellular organization: nucleic acids and proteins. These enzymes ensure the fidelity of transfer of genetic information from the DNA to the protein. They are responsible for attaching amino acid residues to their cognate tRNA molecules by virtue of matching the nucleotide triplet, which is the first step in the protein synthesis. The translation of genetic code into protein sequence is mediated by tRNA, which accurately picks up the cognate amino acids. The attachment of the cognate amino acid to tRNA is catalyzed by aaRSs, which have binding sites for the anticodon region of tRNA and for the amino acid to be attached. The two binding sites are separated by ≈ 76 Å and experiments have shown that the communication does not go through tRNA (Gale et al., 1996). The problem addressed here is how the information of binding of tRNA anticodon near the anticodon binding site is communicated to the active site through the protein structure. These enzymes are modular with distinct domains on which extensive kinetic and mutational experiments and supported by structural data are available, highlighting the role of inter-domain communication (Alexander and Schimmel, 2001). Hence these proteins present themselves as excellent systems for in-silico studies. Various methods involved for the construction of protein structure networks are well established and analyzed in a variety of ways to gain insights into different aspects of protein structure, stability and function (Kannan and Vishveshwara, 1999; Brinda and Vishveshwara, 2005). In the present study, we have incorporated network parameters for the analysis of molecular dynamics (MD) simulation data, representing the global dynamic behavior of protein in a more elegant way. MD simulations have been performed on the available (and modeled) structures of aaRSs bound to a variety of ligands, and the protein structure networks (PSN) of non-covalent interactions have been characterized in dynamical equilibrium. The changes in the structure networks are used to understand the mode of communication, and the paths between the two sites of interest identified by the analysis of the shortest path. The allosteric concept has played a key role in understanding the biological functions of aaRSs. The rigidity/plasticity and the conformational population are the two important ideas invoked in explaining the allosteric effect. We have explored the conformational changes in the complexes of aaRSs through novel parameters such as cliques and communities (Palla et al., 2005), which identify the rigid regions in the protein structure networks (PSNs) constructed from the non-covalent interactions of amino acid side chains. The thesis consists of 7 chapters. The first chapter constitutes the survey of the literature and also provides suitable background for this study. The aims of the thesis are presented in this chapter. Chapter 2 describes various techniques employed and the new techniques developed for the analysis of PSNs. It includes a brief description of well -known methods of molecular dynamics simulations, essential dynamics, and cross correlation maps. The method used for the construction of graphs and networks is also described in detail. The incorporation of network parameters for the analysis of MD simulation data are done for the first time and has been applied on a well studied protein lysozyme, as described in chapter 3. Chapter 3 focuses on the dynamical behavior of protein structure networks, examined by considering the example of T4-lysozyme. The equilibrium dynamics and the process of unfolding are followed by simulating the protein with explicit water molecules at 300K and at higher temperatures (400K, 500K) respectively. Three simulations of 10ns duration have been performed at 500K to ensure the validity of the results. The snapshots of the protein structure from the simulations are represented as Protein Structure Networks (PSN) of non-covalent interactions. The strength of the non-covalent interaction is evaluated and used as an important criterion in the construction of edges. The profiles of the network parameters such as the degree distribution and the size of the largest cluster (giant component) have been examined as a function of interaction strength (Ghosh et al., 2007). We observe a critical strength of interaction (Icritical) at which there is a transition in the size of the largest cluster. Although the transition profiles at all temperatures show behavior similar to those found in the crystal structures, the 500K simulations show that the non-native structures have lower Icritical values. Based on the interactions evaluated at Icritical value, the folding/unfolding transition region has been identified from the 500K simulation trajectories. Furthermore, the residues in the largest cluster obtained at interaction strength higher than Icritical have been identified to be important for folding. Thus, the compositions of the top largest clusters in the 500K simulations have been monitored to understand the dynamical processes such as folding/unfolding and domain formation/disruption. The results correlate well with experimental findings. In addition, the highly connected residues in the network have been identified from the 300K and 400K simulations and have been correlated with the protein stability as determined from mutation experiments. Based on these analyses, certain residues, on which experimental data is not available, have been predicted to be important for the folding and the stability of the protein. The method can also be employed as a valuable tool in the analysis of MD simulation data, since it captures the details at a global level, which may elude conventional pair-wise interaction analysis. After standardizing the concept of dynamical network analysis using Lysozyme, it was applied to our system of interest, the aaRSs. The investigations carried out on Methionyl-tRNA synthetases (MetRS) are presented in chapter 4. This chapter is divided into three parts: Chapter 4A deals with the introduction to aminoacyl tRNA synthetases (aaRS). Classification and functional insights of aaRSs obtained through various studies are presented. Chapter 4B is again divided into parts: BI and BII. Chapter 4BI elucidates a new technique developed for finding communication pathways essential for proper functioning of aaRS. The enzymes of the family of tRNA synthetases perform their functions with high precision, by synchronously recognizing the anticodon region and the amino acylation region, which is separated by about 70Å in space. This precision in function is brought about by establishing good communication paths between the two regions. We have modelled the structure of E.coli Methionyl tRNA synthetase, which is complexed with tRNA and activated methionine. Molecular dynamics simulations have been performed on the modeled structure to obtain the equilibrated structure of the complex and the cross correlations between the residues in MetRS. Furthermore, the network analysis on these structures has been carried out to elucidate the paths of communication between the aminoacyl activation site and the anticodon recognition site (Ghosh and Vishveshwara, 2007). This study has provided the detailed paths of communication, which are consistent with experimental results. A similar study on the (MetRS + activated methionine) and (MetRS+tRNA) complexes along with ligand free-native enzyme has also been carried out. A comparison of the paths derived from the four simulations has clearly shown that the communication path is strongly correlated and unique to the enzyme complex, which is bound to both the tRNA and the activated methionine. The method developed here could also be utilized to investigate any protein system where the function takes place through long distance communication. The details of the method of our investigation and the biological implications of the results are presented in this chapter. In chapter 4BII, we have explored the conformational changes in the complexes of E.coli Methionyl tRNA synthetase (MetRS) through novel parameters such as cliques and communities, which identify the rigid regions in the protein structure networks (PSNs). The rigidity/plasticity and the conformational population are the two important ideas invoked in explaining the allosteric effect. MetRS belongs to the aminoacyl tRNA Synthetases (aaRSs) family that play a crucial role in initiating the protein synthesis process. The network parameters evaluated here on the conformational ensembles of MetRS complexes, generated from molecular dynamics simulations, have enabled us to understand the inter-domain communication in detail. Additionally, the characterization of conformational changes in terms of cliques/communities has also become possible, which had eluded conventional analyses. Furthermore, we find that most of the residues participating in clique/communities are strikingly different from those that take part in long-range communication. The cliques/communities evaluated here for the first time on PSNs have beautifully captured the local geometries in their detail within the framework of global topology. Here the allosteric effect is revealed at the residue level by identifying the important residues specific for structural rigidity and functional flexibility in MetRS. Chapter 4C focuses on MD simulations of Methionyl tRNA synthetase (AmetRS) from a thermophilic bacterium, Aquifex aeolicus. As describe in Chapter 4B, we have explored the communication pathways between the anticodon binding region and the aminoacylation site, and the conformational changes in the complexes through cliques and communities. The two MetRSs from E.coli and Aquifex aeolicus are structurally and sequentially very close to each other. But the communication pathways between anticodon binding region and the aminoacylation site from A. aeolicus have differed significantly with the communication paths obtained from E.coli. The residue composition and cliques/communities structure participating in communication are not similar in the MetRSs of both these organisms. Furthermore the formation of cliques/communities and hubs in the communication paths are more in A. aeolicus compared to E.coli. The participation of structurally homologous linker peptide, essential for orienting the two domains for efficient communication is same in both the organisms although, the residues composition near domain interface regions including the linker peptide is different. Thus, the diversity in the functioning of two different MetRS has been brought out, by comparing the E.coli and Aquifex aeolicus systems. Protein Structure network analysis of MD simulated trajectories of various ligand bound complexes of Escherichia coli Cysteinyl-tRNA synthetase (CysRS) have been discussed in Chapter 5. The modeling of the complex is done by docking the ligand CysAMP into the tRNA bound structure of E.coli Cysteinyl tRNA synthetase. Molecular dynamics simulations have been performed on the modeled structure and the paths of communications were evaluated using a similar method as used in finding communication paths for MetRS enzymes. Compared to MetRS the evaluation of communication paths in CysRS is complicated due to presence of both direct and indirect readouts. The direct and indirect readouts (DR/IR) involve interaction of protein residues with base-specific functional group and sugar-phosphate backbone of nucleic acids respectively. Two paths of communication between the anticodon region and the activation site has been identified by combining the cross correlation information with the protein structure network constructed on the basis of non-covalent interaction. The complete paths include DR/IR interactions with tRNA. Cliques/communities of non-covalently interacting residues imparting structural rigidity are present along the paths. The reduction of cooperative fluctuation due to the presence of community is compensated by IR/DR interaction and thus plays a crucial role in communication of CysRS. Chapter 6 focuses on free energy calculations of aminoacyl tRNA synthetases with various ligands. The free energy contributions to the binding of the substrates are calculated using a method called MM-PBSA (Massova and Kollman, 2000). The binding free energies were calculated as the difference between the free energy of the enzyme-ligand complex, and the free ligand and protein. The ligand unbinding energy values obtained from the umbrella sampling MD correlates well with the ligand binding energies obtained from MM-PBSA method. Furthermore the essential dynamics was captured from MD simulations trajectories performed on E.coli MetRS, A. aeolius MetRS and E.coli CysRS in terms of the eigenvalues. The top two modes account for more than 50% of the motion in essential space for systems E.coli MetRS, A. aeolius MetRS and E.coli CysRS. Population distribution of protein conformation states are looked at the essential plane defined by the two principal components with highest eigenvalues. This shows how aaRSs existed as a population of conformational states and the variation with the addition of ligands. The population of conformational states is converted into Free energy contour surface. From free energy surfaces, it is evident that the E.coli tRNAMet bound MetRS conformational fluctuations are more, which attributes to less rigidity in the complex. Whereas E.coli tRNACys bound CysRS conformational fluctuations are less and this is reflected in the increase in rigidity of the complex as confirmed by its entropic contribution. Future directions have been discussed in the final chapter (Chapter 7). Specifically, it deals with the ab-initio QM/MM study of the enzymatic reaction involved in the active site of E.coli Methionyl tRNA synthetase. To achieve this, two softwares are integrated: the Quantum Mechanics (QM) part includes small ligands and the Molecular Mechanics (MM) part as protein MetRS are handled using CPMD and Gromacs respectively. The inputs for two reactions pathways are prepared. First reaction involves cyclization reaction of homocysteine in the active site of MetRS and the second reaction deals with charging of methionine in the presence of ATP and magnesium ion. These simulations require very high power computing systems and also time of computation is also very large. With the available computational power we could simulate up to 10ps and it is insufficient for analysis. The future direction will involve the simulations of these systems for longer time, followed by the analysis for reaction pathways.

Thèse de doctorat : Aspects moléculaires et cellulaires de la biologie : Strasbourg 1 : 2007. / Thèse soutenue sur un ensemble de travaux. Titre provenant de l'écran-titre. Bibliogr. p. 152-162.