Auswahl der wissenschaftlichen Literatur zum Thema „Cysteine tRNA“

The sulfur-containing nucleosides in transfer RNA (tRNAs) are present in all three domains of life; they have critical functions for accurate and efficient translation, such as tRNA structure stabilization and proper codon recognition. The tRNA modification enzymes ThiI (in bacteria and archaea) and Ncs6 (in archaea and eukaryotic cytosols) catalyze the formation of 4-thiouridine (s4U) and 2-thiouridine (s2U), respectively. The ThiI homologs were proposed to transfer sulfur via cysteine persulfide enzyme adducts, whereas the reaction mechanism of Ncs6 remains unknown. Here we show that ThiI from the archaeon Methanococcus maripaludis contains a [3Fe-4S] cluster that is essential for its tRNA thiolation activity. Furthermore, the archaeal and eukaryotic Ncs6 homologs as well as phosphoseryl-tRNA (Sep-tRNA):Cys-tRNA synthase (SepCysS), which catalyzes the Sep-tRNA to Cys-tRNA conversion in methanogens, also possess a [3Fe-4S] cluster similar to the methanogenic archaeal ThiI. These results suggest that the diverse tRNA thiolation processes in archaea and eukaryotic cytosols share a common mechanism dependent on a [3Fe-4S] cluster for sulfur transfer.

ABSTRACT The complete genomic sequencing of Methanococcus jannaschii cannot identify the gene for the cysteine-specific member of aminoacyl-tRNA synthetases. However, we show here that enzyme activity is present in the cell lysate of M. jannaschii. The demonstration of this activity suggests a direct pathway for the synthesis of cysteinyl-tRNACys during protein synthesis.

The thiolation of biomolecules is a complex process that involves the activation of sulfur. The L-cysteine desulfurase IscS is the main sulfur mobilizing protein in Escherichia coli that provides the sulfur from L-cysteine to several important biomolecules in the cell such as iron sulfur (FeS) clusters, molybdopterin (MPT), thiamine, and thionucleosides of tRNA. Various proteins mediate the transfer of sulfur from IscS to various biomolecules using different interaction partners. A direct connection between the sulfur-containing molecules FeS clusters, thiolated tRNA, and the molybdenum cofactor (Moco) has been identified. The first step of Moco biosynthesis involves the conversion of 5′GTP to cyclic pyranopterin monophosphate (cPMP), a reaction catalyzed by a FeS cluster containing protein. Formed cPMP is further converted to MPT by insertion of two sulfur atoms. The sulfur for this reaction is provided by the L-cysteine desulfurase IscS in addition to the involvement of the TusA protein. TusA is also involved in the sulfur transfer for the thiolation of tRNA. This review will describe the biosynthesis of Moco in E. coli in detail and dissects the sulfur transfer pathways for Moco and tRNA and their connection to FeS cluster biosynthesis.

ABSTRACT Methanococcus maripaludis and Methanocaldococcus jannaschii produce cysteine for protein synthesis using a tRNA-dependent pathway. These methanogens charge tRNACys with l-phosphoserine, which is also an intermediate in the predicted pathways for serine and cystathionine biosynthesis. To establish the mode of phosphoserine production in Methanococcales, cell extracts of M. maripaludis were shown to have phosphoglycerate dehydrogenase and phosphoserine aminotransferase activities. The heterologously expressed and purified phosphoglycerate dehydrogenase from M. maripaludis had enzymological properties similar to those of its bacterial homologs but was poorly inhibited by serine. While bacterial enzymes are inhibited by micromolar concentrations of serine bound to an allosteric site, the low sensitivity of the archaeal protein to serine is consistent with phosphoserine's position as a branch point in several pathways. A broad-specificity class V aspartate aminotransferase from M. jannaschii converted the phosphohydroxypyruvate product to phosphoserine. This enzyme catalyzed the transamination of aspartate, glutamate, phosphoserine, alanine, and cysteate. The M. maripaludis homolog complemented a serC mutation in the Escherichia coli phosphoserine aminotransferase. All methanogenic archaea apparently share this pathway, providing sufficient phosphoserine for the tRNA-dependent cysteine biosynthetic pathway.

Cysteine persulfide (CysSSH) and cysteine polysulfides (CysSSnH, n > 1) are cysteine derivatives that have sulfane sulfur atoms bound to cysteine thiol. Advances in analytical methods that detect and quantify persulfides and polysulfides have shown that CysSSH and related species such as glutathione persulfide occur physiologically and are prevalent in prokaryotes, eukaryotes, and mammals in vivo. The chemical properties and abundance of these compounds suggest a central role for reactive persulfides in cell-regulatory processes. CysSSH and related species have been suggested to act as powerful antioxidants and cellular protectants and may serve as redox signaling intermediates. It was recently shown that cysteinyl-tRNA synthetase (CARS) is a new cysteine persulfide synthase. In addition, we discovered that CARS is involved in protein polysulfidation that is coupled with translation. Mitochondrial activity in biogenesis and bioenergetics is supported and upregulated by CysSSH derived from mitochondrial CARS. In this review article, we discuss the mechanisms of the biosynthesis of CysSSH and related persulfide species, with a particular focus on the roles of CARS. We also review the antioxidative and anti-inflammatory actions of persulfides.

As aminoacil tRNA sintetases (AaRSs) são enzimas essenciais na síntese de proteínas assegurando a correta relação entre os aminoácidos e seus tRNA cognatos. O genoma mitocondrial dos tripanossomatídeos perdeu os genes codificantes dos tRNAs, assim os tRNA mitocondriais são codificados no núcleo e importados do citoplasma. O código genético do kinetoplasto desvia do código genético pela utilização do códon de terminação UGA para a decodificação do códon do triptofano. Um único gene codificando o tRNATrp(CCA) observado no genoma de Leismania é responsável pela incorporação do aminoácido triptofano durante a síntese proteíca na mitocôndria. Para decodificar os dois códons do Trp (UGA e UGG) a base na posição 34 do tRNATrp(CCA) passa por um evento de editoração, convertendo o ribunuclotídeo C34 em U34, produzindo o tRNATrp(UCA) capaz de decodificar o códon UGA. Nesse trabalho foram caracterizadas duas triptofanil tRNA sintetases de Leishmania major. De acordo com experimentos de ?western blotting? e análises ?in silico? das seqüências de aminoácidos, uma enzima tem localização citoplasmática (LmTrpRS1) enquanto a outra mitocondrial (LmTrpRS2). Os mRNAs dos dois genes foram definidos por experimentos de 5? e 3? RT-PCR. As duas enzimas foram clonadas em diversos vetores de expressão procariotos e eucariotos. A LmTrpRS1 foi obtida somente na fração insolúvel, já a LmTrpRS2 foi obtida na fração solúvel quando clonada no vetor de expressão pET28a. Esta porém mostrou-se instável precipitando rapidamente após sua purificação. Os ensaios enzimáticos realizados com a mesma mostraram que ela é capaz de reconhecer os tRNAsTrp editado e não editado. Modelagem molecular por homologia com as duas proteínas foi realizada usando a proteína citoplasmática humana como molde, para estudar a interação entre a proteína e o tRNATrp. Xylella fastidiosa é um bactéria gram negativa limitada ao xilema, responsável por um grande número de doenças economicamente importantes, como a doença de Pierces em videiras, Clorose variegata do Citrus (CVC) e a doença da requeima das folhas em outras plantas incluindo, amendoeira, ameixeira, louro, amoreira e café. Em todos os casos a X. fastidiosa afeta o xylema da planta causando redução na produção de frutos. Nesse trabalho nós mostramos a estrutura da Xylellaína, uma cisteíno protease desse patógeno. A estrutura foi resolvida por dispersão anômala a um único comprimento de onda, utilizando cristais de xylellaína selenometionina substituídos. A estrutura da Xylellaína foi refinada até 1,65 Å de resolução, mostrando enovelamento similar às proteínas da família da papaína, porém algumas características interessantes como uma região N-terminal composta por 38 aminoácidos cobrindo o sulco ativo da enzima, um intrigante ribonucleotídeo encontrado fora do sítio ativo da enzima e um ?loop? semelhante ao ?loop? de oclusão presente na catepsina B.
The aminoacyl tRNA synthetases (aaRSs) are essential enzymes in protein synthesis that ensure the correct match between amino acids and their cognate tRNAs. The mitochondrial (kinetoplast) genome of trypanossomatids lacks tRNA genes, and therefore nucleus-encoded tRNAs are imported from the cytoplasm, the kinetoplast genetic code deviates from the universal code in that UGA instead of UGG encodes for tryptophan. A single nucleus-encoded tRNATrp(CCA) is responsible for Trp insertion during organellar protein synthesis. To decode both Trp codons (UGA and UGG), tRNATrp(CCA) undergoes a single C to U editing event at position 34 of the anticodon yielding to versions of the tRNA in the mitochondria with anticodon CCA and UCA, permitting UGA decoding. This work have characterized two Leishmania major tryptophanyl-tRNA synthetase, acording western blotting experiments and ?in silico? sequence analisis one of cytoplasmatic localization (LmTrpRS1) and another from mitochondria localization (LmTrpRS2). The mature mRNA transcripts for both genes were defined by 5? and 3? RT-PCR. Both enzymes were cloned into several expressions vectors. LmTrpRs1 was obtained as an insoluble protein and LmTrpRs2 expressed into the soluble fraction in pET28a expression system. LmTrpRS2 protein, however, is unstable precipitating shortly after purification. The enzymatic assay showed that this enzyme is able to recognize both tRNATrp. Molecular modeling for LmTrpRS1 and LmTrpRS2 were constructed using the cytoplasmatic human tryptophanyl tRNA synthetase as a model, to study the interaction between proteins and tRNATrp. Xylella fastidiosa is a xylem-limited, gram-negative bacteria responsible for a large number of economically important plant diseases, such as Pierces disease in grapevines, citrus variegated chlorosis (CVC) in sweet oranges and leaf scorch diseases in other plants, including almond, plum, oleander, mulberry and coffee. In all cases, X. fastidiosa infects the plant xylem and impairs fruit production. Here, we report the crystal structure of xylellain, a cystein protease from X. fastidiosa. The structure was solved by single-wavelength anomalous dispersion (SAD) using seleno-methionine containing xylellain crystals. The final structure of Xylellaína was refined against the best native data set (1.65 Å) showing R/Rfree= 17/21. Xylellain shares fold similar to Papain like Family, but contains some interesting features, like a 38 N-terminal tail covering the active site cleft; one intriguing ribonucleotide found outside the active site and one loop that resemble the ocluding loop from cathepsin B.

During protein synthesis, elongation factor Tu (EF-Tu) delivers aminoacyl-tRNA (aa-tRNA) to the A-site of mRNA-programmed ribosomes in a GTP-dependent manner. To enable future studies on the functional and structural requirement of EF-Tu’s function, a Cysteine-free variant of EF-Tu was constructed suitable for subsequent labelling of the protein and use in kinetic studies. Here, the kinetic properties of three Cysteine-less EF-Tu variants are reported, demonstrating that only the variant with the Alanine substitution in position 81 retains wild-type activity with respect to the interaction with guanine nucleotides, aa-tRNA and the ribosome. To explore a possible tRNA independent pathway for the GTPase activation signal, three residues in domain II of EF-Tu (Arginine 279, Glycine 280, Arginine 283) were mutated; the activity of EF-Tu variants were analyzed. Results suggest that these residues are indeed required for efficient ribosome-dependent stimulation of the GTPase activity of EF-Tu.
x, 85 leaves : ill. (some col.) ; 29 cm

碩士
國立中央大學
生命科學系
106
Aminoacyl-tRNA synthetases (aaRSs) belongs to a group of enzymes necessary for protein synthesis. Their main function is to attach an amino acid to its corresponding tRNA to form aminoacyl-tRNA, which is then brought to the ribosome for protein synthesis. One aaRSs corresponds to each amino acid. Previous studies have shown that although some bacteria lack CysRS, they can synthesize Cys-tRNACys through an indirect pathway using serine as a substrate, whereas most bacteria have a CysRS to synthesize Cys-tRNACys. Herein, we present the evidence that that Mycobacterium abscessus possesses two cysteinyl-tRNA synthetase (CysRS) homologues genes that CysS1 and CysS2 (which encode MaCysRS1 and MaCysRS2, respectively). Two homologous CysRSs in M. abscessus have a 37% identity, 80% similarity and 37-42% identity with E. coli CysRS, it is totally different with CysRS in Escherichia coli. Further sequence and phylogenetic analyses showed that MaCysRS2 is actually a mycothiol cysteine-ligase (MshC) which is involved in Mycothiol (MSH) synthesis as a protective thiol. It is not only different protein involved in protein synthesis but also lacks anticodon-binding domain. The result of complementation assay showed that both MaCysS1 and MaCysS2 were moderately expressed in the yeast but failed to complement the cytoplasmic function of the knockout strain, i.e., these two genes cannot provide the required CysRS activity to support the growth of the null allele on 5’-FOA medium. However, if a mitochondrial targeting signal (MTS) was attached to the N-terminal of the MaCysRS1, the fusion protein successfully rescued the growth defect of the knockout strain on YPG, suggesting that this fusion protein can substitute the mitochondrial activity of yeast CysRS. In contrast, MaMshC, even fused to an MTS, could not do so, probably because MaMshC lacks an anticodon-binding domain (ABD). Most surprisingly, fusion of a tRNA-binding domain of Arc1p to MTS-MaMshC, yielding an MTS-MaMshC-Arc1p (M+C), enabled the enzyme to restore the growth of the yeast knockout strain on YPG. This result shows that MaMshC, a bacterial protective thiol-producing enzyme, can be converted to a functional cysteinyl-tRNA synthetase through fusion of a non-specific tRNA-binding domain. Keywords: CysRS, MaCysRS, Arc1p

A significant number of known human genetic diseases is associated with nonsense mutations leading to the introduction of a premature termination codon into the coding sequence. A termination codon can be read through by its near-cognate tRNA (tRNA with two anticodon nucleotides base-pairing with a stop codon); potentially generating C-terminally extended protein variants. In yeast, UGA stop codon was described to be read through by tRNA-Trp and tRNA-Cys. Similar was observed for tRNA-Trp in human HEK293T cell line. The aim of this thesis was to investigate if human tRNA-Cys can act as a near-cognate tRNA in human HEK293T cell line. There are two isoacceptors which constitute the tRNA-Cys family, with ACA and GCA anticodon. There are 1 and 23 isodecoders to the ACA and GCA anticodons, respectively. Here, altogether as many as nine tRNA-Cys isodecoders (distinct in their sequence and with varying levels of expression) were tested for their ability to increase UGA readthrough in HEK293T using p2luci and pSGDluc dual-luciferase reporter vectors. In both p2luci and pSGDluc, we observed that at least one tRNA-Cys isodecoder, tRNA-Cys-GCA-4-1, is capable of significantly elevating the UGA readthrough levels when overexpressed in HEK293T. This indicates that similarly to yeast, tRNA-Cys is capable of...

A study of the recognition of tRNA^(Cys) by E. coli cysteinyl-tRNA synthetase using in vivo and in vitro methods was performed. All three anticodon nucleotides, the discriminator U73, and some element(s) within the tertiary domain (the D stem/loop, the TΨC stem/loop and extra loop) are important for recognition; the anticodon stem and acceptor stem appear to contain no essential elements. A T7 RNA polymerase transcribed tRNA^(Cys) is a 5.5-fold worse substrate than native tRNA^(Cys)(in terms of the selectivity constant, k_cat/K_m) mainly due to an increase in K_m. This may reflect recognition of modified nucleotides or subtle effects on the folding of the tRNA. The greatest loss of specificity caused by mutation of a single nucleotide occurs when the discriminator U73 is changed; k_cat/K_m declines 3 to 4 orders of magnitude depending on the substitution. Mutations in the wobble nucleotide of the anticodon also cause reductions in the selectivity constant of 3 orders of magnitude, while mutations in the other anticodon nucleotides caused lesser effects. Interestingly, a C35A mutation had no effect on aminoacylation by the cysteinyl-tRNA synthetase. Several amber suppressor tRNAs were constructed whose in vivo identity did not correlate with their in vitro specificity, indicating the need for both types of experiments to understand the factor(s) which maintain tRNA specificity. Future in vitro experiments will attempt to explain the in vivo discrimination between the glycine, phenylalanine, and cysteine tRNAs by the cysteinyl-tRNA synthetase. Finally, these results suggest that the notion that a small set of isoacceptor specific elements define tRNA identity (the socalled "second genetic code") is incorrect. A better model is based on competition between synthetases for tRNA substrates which contain differing amounts of partially overlapping identity determinants.

博士
國立臺灣大學
化學研究所
103
Ensuring labeling specificity for fluorophore labeling is a challenge in fluorescent spectroscopy. This biophysical technique is one of the primary research tools for studying cotranslational protein folding which studies protein conformation before it has been released from the ribosome. Due to the nature of ribosome-bound nascent chains (RNCs), fluorescent labeling must be coupled with translation during which tRNA acts as the carrier of fluorescent amino acid. In this work, a novel overexpressed suppressor tRNAcysAmber is developed for the production of BODIPY FL-labeled RNCs. In vitro transcription as well as overexpression is tested as the methods of suppressor tRNA production. In order to simplify the purification procedures, Bacillus subtilis tRNAcysAmber has been selected for its distinctive sequence from any endogenous E. coli RNA. In a single purification step, ample amounts of tRNAcysAmber have been obtained. As cysteinyl-tRNA synthetase was previously shown to aminoacylate tRNAcysAmber with low efficiency, several point mutations were introduced into the C-terminus of cysteinyl-tRNA synthetase to compensate for the Amber mutation in the tRNA anticodon loop. Out of the cysteinyl-tRNA synthetase mutants, D436S mutant is shown with improved aminoacylation efficiency and specificity towards tRNAcysAmber. In addition, overexpressed BODIPY FL-cysteinyl-tRNAcysAmber indicates improved stability of this tRNA compared to the in vitro transcribed tRNA. Applying this tRNA, the dynamics of single-labeled RNCs by time-resolved anisotropy was studied to reveal information about the protein folding on the ribosome. The natively unfolded phosphorylated insulin receptor domain (PIR) and the zinc-induced folding in zinc-finger RNC help to correlate the fluorescence correlation time with different nascent chain movements. To further study the impact of chaperones on the RNC dynamics at different stages of translation, Entner-Douderoff aldolase (Eda) RNCs with four predetermined chain length are generated in either the wild-type or chaperone-depleted cell-free system. By applying BODIPY FL-cysteinyl-tRNAcysAmber, our results indicate that Eda may start folding without chaperones after approximately half of the protein emerges from the ribosome. In addition, chaperones increase the nascent chain confinement in the full length Eda RNC, which may cause the decrease the binding of the trigger factor, a co-translational chaperone, due to growing hindrances between RNCs and chaperones. Overall, the facile preparation of suppressor tRNA for labeling with fluorophores is demonstrated together with the application of single-residue labeled nascent chains in studying the effect of chaperones on the RNC dynamics.

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.