Dissertations / Theses on the topic 'TrnK intron'

The recent call to reconstruct a detailed picture of the tree of life for all organisms has forever changed the field of molecular phylogenetics. Sequencing technology has improved to the point that scientics can now routinely sequence complete plastid/mitochondrial genomes and thus, vast amounts of data can be used to reconstruct phylogenies. These data are accumulating in DNA sequence repositories, such as GenBank, where everyone can benefit from the vast growth of information. The trend of generating genomic-region rich datasets has far outpaced the expasion of datasets by sampling a broader array of taxa. We show here that expanding a dataset both by increasing genomic regions and species sampled using GenBank data, despite the inherent missing DNA that comes with GenBank data, can provide a robust phylogeny for the plant order Caryophyllales (angiosperms). We also investigate the utility of trnK intron in phylogeny reconstruction at relativley deep evolutionary history (the caryophyllid order) by comparing it with rapidly evolving matK. We show that trnK intron is comparable to matK in terms of the proportion of variable sites, parsimony informative sites, the distribution of those sites among rate classes, and phylogenetic informativness across the history of the order. This is especailly useful since trnK intron is often sequenced concurrently with matK which saves on time and resources by increasing the phylogenetic utility of a single genomic region (rapidly evolving matK/trnK). Finally, we show that the inclusion of RNA edited sites in datasets for phylogeny reconstruction did not appear to impact resolution or support in the Gnetales indicating that edited sites in such low proportions do not need to be a consideration when building datasets. We also propose an alternate start codon for matK in Ephedra based on the presense of a 38 base pair indel in several species that otherwise result in pre-mature stop codons, and present 20 RNA edited sites in two Zamiaceae and three Pinaceae species.
Ph. D.

Evolution of Piperales – matK gene and trnK intron sequence data reveal lineage specific resolution contrast. Piperales are one of the largest basal angiosperm orders with a nearly worldwide distribution. The order includes three species rich genera, Piper (ca. 1,000 species), Peperomia (ca. 1,500-1,700 species), and Aristolochia s. l. (ca. 500 species). Sequences of the matK gene and the non-coding trnK group II intron are analysed for a dense set of 105 taxa representing all families (except Hydnoraceae) and all generic segregates (except Euglypha within Aristolochiaceae) of Piperales. A large number of highly informative indels are found in the Piperales trnK/matK dataset. Within a narrow region approximately 500 nt downstream in the matK coding region (CDS), a length variable simple sequence repeat (SSR) expansion segment occurs, in which insertions and deletions have led to short frame-shifts. These are corrected shortly afterwards, resulting in a maximum of 6 amino acids being affected. Furthermore, additional non-functional matK copies were found in Zippelia begoniifolia, which can easily be discriminated from the functional open reading frame (ORF). The trnK/matK sequence data fully resolve relationships within Peperomia, whereas they are not effective within Piper. The resolution contrast is correlated with the rate heterogenity between those lineages. Parsimony, Bayesian and likelihood analyses result in virtually the same topology, and converge on the monophyly of Piperaceae and Saururaceae. Lactoris gains high support as sister to Aristolochiaceae subf. Aristolochioideae, but the different tree inference methods yield conflicting results with respect to the relationships of subfam. Asaroideae. In Piperaceae, a clade formed by the monotypic genus Zippelia and the small genus Manekia (=Sarcorhachis) is sister to the two large genera Piper and Peperomia. Systematics of pipevines – Combining morphological and fast-evolving molecular characters to investigate the relationships within subfamily Aristolochioideae (Aristolochiaceae) A combined phylogenetic analysis of the Aristolochioideae was conducted based on 72 morphological characters and molecular datasets (matK gene, trnK intron, trnL intron, trnL-trnF spacer). The analysis sampled 33 species as the ingroup, including two species of Thottea and 30 species of Aristolochia and the monotypic genus Euglypha, which represent all the infrageneric taxa formally described; Saruma henryi and Asarum caudatum were used as the outgroup. The results corroborate a sister-group relationship between Thottea and Aristolochia, and the paraphyly of Aristolochia with respect to Euglypha that consequently should be included into Aristolochia. Two of the three subgenera within Aristolochia (Isotrema and Pararistolochia) are shown to be monophyletic, whereas the signal obtained from the different datasets about the relationships within subg. Aristolochia is low and conflicting, resulting in collapsed or unsupported branches. The relationship between the New World and the Old World species of subgenus Aristolochia is conflictive because morphological data support these two groups as monophyletic, whereas molecular data show the monophyletic Old World species of Aristolochia nested within the New World species. A sister group relationship is proposed between A. lindneri and pentandrous species, which suggests that a group of five species from central and southern South America (including A. lindneri) could be monophyletic and sister to Aristolochia subsection Pentandrae, a monophyletic taxon consisting of about 35 species from southern USA, Mesoamerica, and the West Indies. Colonisation, phylogeography and evolution of endemism in Mediterranean Aristolochia (Aristolochiaceae). This study provides evidence for a multiple colonisation of the western Old World from Asian ancestors within Aristolochia section Diplolobus (subsection Aristolochia and Podanthemum). Within subsection Podanthemum it is assumed, that the colonisation of the African continent happened at least two times independently. In contrast, for subsection Aristolochia, a rapid morphological radiation in the Near East (or close to this area) with subsequent star like colonisation of the different current distribution areas, which is not paralleled on the molecular level, appears to be more likely. Phylogenetic tree reconstruction is unsupported for these clades, but most clades are highly supported as monophyletic. Interestingly the Mediterranean and temperate Eurasian species, which are morphologically distinct (A. pistolochia, A. clematitis) are not clustering within the main clades, but are independent lineages. Analogue, A. rigida a species from Somalia is well-supported sister to the subsection Aristolochia. Within subsection Podanthemum the colonisation event from an Asian ancestor is clearly traceable, whereas in subsection Aristolochia the path is not traceable, since the ancestors are extinct or not present in the connecting areas. Within the Mediterranean, Near East and Caucasian species of subsection Aristolochia two morphologically and biogeographically well supported groups can be identified: the Near East/Caucasian species and the West Mediterranean species. The previous groupings for the latter, based on morphological characters, could be substantiated only partly by our results. This study provides the first phylogeny of all West Mediterranean species. In addition an independent complex is established including some micro endemic species. The phylogenetic results are discussed with respect to biogeography, and morphology, to give a first insight into the radiation and colonisation of the genus Aristolochia in the Mediterranean region. Universal primers for a large cryptically simple cpDNA microsatellite region in Aristolochia. We provide a new and valuable marker to study species relationships and population genetics in order to trace evolutionary, ecological, and conservational aspects in the genus Aristolochia. Universal primers for amplification and subsequent sequencing of a chloroplast microsatellite locus inside the trnK intron are presented. Utility of the primers has been tested in 32 species representing all clades of Aristolochia, including population studies within the A. pallida complex, A. clusii and A. rotunda. The microsatellite region is characterized as a (AnTm)k repeat of 22–438 bp containing a combination of different repeats arranged as ‘cryptically simple’. Trapped! Pollination of Aristolochia pallida Willd. in the Mediterranean A first study of the pollination biology of a Mediterranean Aristolochia species in its natural habitat is presented. 183 flowers of Aristolochia pallida were investigated, which in total contained 73 arthropods, dominated by two groups of Diptera, Sciaridae (37%) and Phoridae (19%). However, only Phoridae are regarded as potential pollinators, since pollen has been found exclusively on the body surfaces of these insects. All Phoridae belong to the genus Megaselia and are recognised as four undescribed species. The measurements of flower and insect dimensions suggest that size is an important constrain for successful pollination: 1) the insects must have a definitive size for being able to enter the flower and 2) must be able to get in touch with the pollen. Only very few insect groups found in Aristolochia pallida fulfil these size requirements. However, size alone is not a sufficient constrain as too many fly species of the same size might be trapped but not function as pollinators. Instead, specific attraction is required as otherwise pollen is lost. Since all trapped Phoridae are males, a chemical attraction (pheromones) is proposed as an additional constrain. Since A. pallida flowers are protogynous, the record of Megaselia loaded with pollen found in a flower during its female stage proves that this insect must have been visited at least one different flower during its male stage before. Further on, this observation provides strong evidence that the flowers are cross-pollinated. All these factors indicate a highly specialised pollination of Aristolochia pallida by Megaselia species.

Evolution of Piperales – matK gene and trnK intron sequence data reveal lineage specific resolution contrast. Piperales are one of the largest basal angiosperm orders with a nearly worldwide distribution. The order includes three species rich genera, Piper (ca. 1,000 species), Peperomia (ca. 1,500-1,700 species), and Aristolochia s. l. (ca. 500 species). Sequences of the matK gene and the non-coding trnK group II intron are analysed for a dense set of 105 taxa representing all families (except Hydnoraceae) and all generic segregates (except Euglypha within Aristolochiaceae) of Piperales. A large number of highly informative indels are found in the Piperales trnK/matK dataset. Within a narrow region approximately 500 nt downstream in the matK coding region (CDS), a length variable simple sequence repeat (SSR) expansion segment occurs, in which insertions and deletions have led to short frame-shifts. These are corrected shortly afterwards, resulting in a maximum of 6 amino acids being affected. Furthermore, additional non-functional matK copies were found in Zippelia begoniifolia, which can easily be discriminated from the functional open reading frame (ORF). The trnK/matK sequence data fully resolve relationships within Peperomia, whereas they are not effective within Piper. The resolution contrast is correlated with the rate heterogenity between those lineages. Parsimony, Bayesian and likelihood analyses result in virtually the same topology, and converge on the monophyly of Piperaceae and Saururaceae. Lactoris gains high support as sister to Aristolochiaceae subf. Aristolochioideae, but the different tree inference methods yield conflicting results with respect to the relationships of subfam. Asaroideae. In Piperaceae, a clade formed by the monotypic genus Zippelia and the small genus Manekia (=Sarcorhachis) is sister to the two large genera Piper and Peperomia. Systematics of pipevines – Combining morphological and fast-evolving molecular characters to investigate the relationships within subfamily Aristolochioideae (Aristolochiaceae) A combined phylogenetic analysis of the Aristolochioideae was conducted based on 72 morphological characters and molecular datasets (matK gene, trnK intron, trnL intron, trnL-trnF spacer). The analysis sampled 33 species as the ingroup, including two species of Thottea and 30 species of Aristolochia and the monotypic genus Euglypha, which represent all the infrageneric taxa formally described; Saruma henryi and Asarum caudatum were used as the outgroup. The results corroborate a sister-group relationship between Thottea and Aristolochia, and the paraphyly of Aristolochia with respect to Euglypha that consequently should be included into Aristolochia. Two of the three subgenera within Aristolochia (Isotrema and Pararistolochia) are shown to be monophyletic, whereas the signal obtained from the different datasets about the relationships within subg. Aristolochia is low and conflicting, resulting in collapsed or unsupported branches. The relationship between the New World and the Old World species of subgenus Aristolochia is conflictive because morphological data support these two groups as monophyletic, whereas molecular data show the monophyletic Old World species of Aristolochia nested within the New World species. A sister group relationship is proposed between A. lindneri and pentandrous species, which suggests that a group of five species from central and southern South America (including A. lindneri) could be monophyletic and sister to Aristolochia subsection Pentandrae, a monophyletic taxon consisting of about 35 species from southern USA, Mesoamerica, and the West Indies. Colonisation, phylogeography and evolution of endemism in Mediterranean Aristolochia (Aristolochiaceae). This study provides evidence for a multiple colonisation of the western Old World from Asian ancestors within Aristolochia section Diplolobus (subsection Aristolochia and Podanthemum). Within subsection Podanthemum it is assumed, that the colonisation of the African continent happened at least two times independently. In contrast, for subsection Aristolochia, a rapid morphological radiation in the Near East (or close to this area) with subsequent star like colonisation of the different current distribution areas, which is not paralleled on the molecular level, appears to be more likely. Phylogenetic tree reconstruction is unsupported for these clades, but most clades are highly supported as monophyletic. Interestingly the Mediterranean and temperate Eurasian species, which are morphologically distinct (A. pistolochia, A. clematitis) are not clustering within the main clades, but are independent lineages. Analogue, A. rigida a species from Somalia is well-supported sister to the subsection Aristolochia. Within subsection Podanthemum the colonisation event from an Asian ancestor is clearly traceable, whereas in subsection Aristolochia the path is not traceable, since the ancestors are extinct or not present in the connecting areas. Within the Mediterranean, Near East and Caucasian species of subsection Aristolochia two morphologically and biogeographically well supported groups can be identified: the Near East/Caucasian species and the West Mediterranean species. The previous groupings for the latter, based on morphological characters, could be substantiated only partly by our results. This study provides the first phylogeny of all West Mediterranean species. In addition an independent complex is established including some micro endemic species. The phylogenetic results are discussed with respect to biogeography, and morphology, to give a first insight into the radiation and colonisation of the genus Aristolochia in the Mediterranean region. Universal primers for a large cryptically simple cpDNA microsatellite region in Aristolochia. We provide a new and valuable marker to study species relationships and population genetics in order to trace evolutionary, ecological, and conservational aspects in the genus Aristolochia. Universal primers for amplification and subsequent sequencing of a chloroplast microsatellite locus inside the trnK intron are presented. Utility of the primers has been tested in 32 species representing all clades of Aristolochia, including population studies within the A. pallida complex, A. clusii and A. rotunda. The microsatellite region is characterized as a (AnTm)k repeat of 22–438 bp containing a combination of different repeats arranged as ‘cryptically simple’. Trapped! Pollination of Aristolochia pallida Willd. in the Mediterranean A first study of the pollination biology of a Mediterranean Aristolochia species in its natural habitat is presented. 183 flowers of Aristolochia pallida were investigated, which in total contained 73 arthropods, dominated by two groups of Diptera, Sciaridae (37%) and Phoridae (19%). However, only Phoridae are regarded as potential pollinators, since pollen has been found exclusively on the body surfaces of these insects. All Phoridae belong to the genus Megaselia and are recognised as four undescribed species. The measurements of flower and insect dimensions suggest that size is an important constrain for successful pollination: 1) the insects must have a definitive size for being able to enter the flower and 2) must be able to get in touch with the pollen. Only very few insect groups found in Aristolochia pallida fulfil these size requirements. However, size alone is not a sufficient constrain as too many fly species of the same size might be trapped but not function as pollinators. Instead, specific attraction is required as otherwise pollen is lost. Since all trapped Phoridae are males, a chemical attraction (pheromones) is proposed as an additional constrain. Since A. pallida flowers are protogynous, the record of Megaselia loaded with pollen found in a flower during its female stage proves that this insect must have been visited at least one different flower during its male stage before. Further on, this observation provides strong evidence that the flowers are cross-pollinated. All these factors indicate a highly specialised pollination of Aristolochia pallida by Megaselia species.

I studied evolutionary history in the angiosperm order Apiales, with a special emphasis on interactions between form, time, and space. Four broad categories of problems were addressed: interfamilial relationships in Apiales, the assignment of genera traditionally assigned to the Apiaceae subfamily Hydrocotyloideae, the estimation of divergence times of the major clades, and the reconstruction of the biogeographic history of Apiales. We used molecular markers with different evolutionary properties and rates derived from the plastid (trnD-trnT and rpl16), nuclear (RPB2), and mitochondrial (nad1 intron 2) genomes, from more than 250 species representing all major clades in the order. The nuclear RPB2 region exhibited evidence of at least six duplication events in Apiales and provided a rich source of information for understanding the origins of polyploid lineages, especially in Araliaceae. Sequence comparisons among the copies show that exon regions are highly conserved. All copies appear to be functional but may have undergone subfunctionalization. Phylogenetic analyses of the three genomes suggest that Hydrocotyloideae should be divided into as many as six evolutionary lineages, but that most taxa should be included in subfamilies Azorelloideae and Mackinlayoideae. Relationships among and within the major clades of Azorelloideae need further analyses since many genera appeared non-monophyletic (e.g., Azorella, Schizeilema, and Eremocharis). Mackinlayoideae appeared as the earliest diverging lineage of Apiaceae, but the plastid and nuclear trees were incongruent in the placement of the Platysace clade relative to Mackinlayoideae and the rest of Apiaceae. Among the remaining clades of suborder Apiineae, Myodocapaceae appeared sister to Apiaceae in both plastid and nuclear trees, preceded by the divergence of Araliaceae and then Pittosporaceae. At the base of the gene trees in Apiales, Griseliniaceae and Torricelliaceae formed successive sisters to Apiineae. The placement of Pennantiaceae as sister to the rest of Apiales was confirmed by plastid data, but was not found in the nuclear trees. The order appears to have originated in the Cretaceous, with Apiineae having an age of c. 100 Mya. Australasia appears to be the most likely center of origin for Apiineae and most of its major clades, except Azorelloideae (South America) and Apioideae-Saniculoideae (sub-Saharan Africa).

Evolutionary relationships within and among three Astragalus sections (Incani DC., Hypoglottidei DC., and Dissitiflori DC.) that were native to Turkey were inferred from variations of nucleotide sequences of both chloroplast and nuclear genome regions. In the current study, Fifty-six species included in the three Astragalus sections were utilized to figure out phylogenetic relationships and estimate evolutionary divergence time based on DNA sequence of trnL intron (trnL5&rsquo
-L3&rsquo
) , trnL3&rsquo
-F(GAA) (trnL-F intergenic spacer), trnV intron, matK (maturase kinase) cpDNA (chloroplast) and ITS (internal transcribed spacer) nDNA (nuclear) regions. Fifty-six Astragalus species with their replicas and one Cicer species as outgroup were analyzed by polymerase chain reaction amplification and DNA sequencing methods. Eleven unknown samples were also used in the current study to understand their section and species name. The results of the study indicated that unknown A35 and A52 samples could be named as A. dasycarpus, while unknown A65 and A66 samples as A. ovatus and lastly unknown A2 sample as A. nitens or A. aucheri. Section of unknown A3, A16, A20, A108, A109 and A110 samples were determined as Incani, but the exact species identification of these samples were not possible because of their close phylogenetic associations with more than one species. Highest genetic diversity was observed when the DNA sequences of ITS nrDNA (nuclear ribosomal) region comprising three subregions as ITS1, 5.8S and ITS2 was used, while the lowest one was calculated when DNA sequence of trnL-F cpDNA region was analyzed. The genetic divergence between Incani and Dissitiflori sections was highest whereas between Hypoglottidei and Dissitiflori was lowest based on all used regions. To figure out phylogenetic relationships among Astragalus species distributed in Turkey and in other regions of the World, DNA sequences of studied regions of foreign samples were collected from the NCBI database and were evaluated with DNA sequence of Turkish species used in the curent study. The Iranian samples either scattered in the phylogenetic tree or attached to our samples externally. South and North American samples (New World Astragalus or Neo Astragalus group) were nested within a different subcluster, which was located in the main cluster produced by samples of Old World Astragalus group (Turkish samples). With these results, we can say that New World Astragalus group is monophyletic and diverged from Old World Astragalus group. Evolutionary divergence time for Astragalus genus was estimated as about 12.5 - 14.5 million years (Ma), and that of New World Astragalus group as 5.0 - 4.0 Ma when rates of nucleotide substitutions of trnL intron and matK cpDNA regions were analyzed. In addition to evolutionary divergence time estimation for Astragalus and New World Astragalus group, divergence times among used three sections of the genus were also calculated by using DNA sequences of trnL, trnV intron and matK cpDNA regions and results indicated that Hypoglottidei and Dissitiflori sections diverged about 5.0-7.0 million years later than Incani section.

RNA molecules are diverse in form and function. They include messenger RNAs (mRNAs) that are templates for proteins, splice products such as introns that can generate functional noncoding RNAs, and a slew of smaller RNAs such as transfer RNAs (tRNAs) that help decode mRNAs into proteins. RNAs can show distinct patterns of subcellular localization that play an important role in protein localization. However, RNA distribution in cells is incompletely understood, with prior studies focusing primarily on RNAs that are long (>200 nucleotides), fully processed, and polyadenylated. We examined the distribution of RNAs in neurons. Neuronal compartments can be separated by long distances and play distinct roles, raising the possibility that RNA localization is especially overt and functionally meaningful in these cells. In our exploration, we physically dissected projections from cell bodies of neurons from the rat brain and sequenced total RNA. We describe two main findings. First, we identified excised introns that are enriched in neuronal projections and confirmed their localization by single- molecule fluorescence in situ hybridization. These are a previously unknown set of circular RNAs in neuronal projections: tailless lariats that possess a non- canonical C branchpoint. Second, we observed a highly abundant population of small (20-150 nucleotide) RNAs in neuronal projections, most of which are tRNAs. For both circular introns and tRNAs, we did not observe known RNA localization signals. Thus, many types of RNA, if sufficiently stable, appear free to diffuse to distant locations, their localization perhaps aided by the movement of large organelles in the confines of neuronal projections. Our survey of RNA molecules across subcellular compartments provides a foundation for investigating the function of these molecules and the mechanisms that localize them.

La reconnaissance specifique du trna#m#e#t par la metrs de levure fait intervenir des determinants de structure primaire, secondaire et tertiaire, a savoir: la boucle anticodon et la base discriminatrice, le bras d et la conformation en l respectivement. Dans la deuxieme partie de ce travail, nous avons aborde l'etude de la reconnaissance entre un trna et son aminoacyl-trna synthetase sous un aspect evolutif. La phers mitochondriale de levure se distingue de son homologue cytoplasmique par un modele d'interaction enzyme-trna simplifie faisant appel a deux points de contacts majeurs: l'anticodon et le bras accepteur. Nous avons egalement analyse l'identite d'un trna#i#l#e ambre de levure dans un contexte bacterien in vivo demontrant que la paire u#3#0-g#4#0 dans le bras anticodon est capable de moduler l'identite et l'efficacite de suppression du trna. Enfin, nous avons etudie le role des bases modifiees i et psi en premiere position de l'anticodon des trna#i#l#e majeur et mineur de levure. Les effets de substitution de bases en position 34 de l'anticodon du trna#i#l#e par ordre decroissant donnent la hierarchie suivante: i psi > g u > a >> c suggerant que l'ilers reconnait une structure c=o n-h de la base 34 dependante d'une conformation particuliere de la boucle anticodon liee a la presence d'une autre base modifiee. D'autre part, les proprietes traductionnelles liees a la presence de psi#3#4, nous ont amene a etudier la biosynthese de ces psis dans l'anticodon du trna#i#l#e mineur. Nous avons montre qu'elle est: i) est intron dependante, ii) n'est pas sensible aux mutations ponctuelles dans la region de l'anticodon (positions 32 a 37), iii) ne requiert pas la structure tridimensionnelle du trna, iv) les deux psis en positions 34 et 36 sont probablement synthetisees par un meme enzyme

碩士
中央警察大學
鑑識科學研究所
90
The intron of trnL(UAA) and intergenic spacer (IGS) of trnL(UAA) - trnF(GAA) in the chloroplast genome were used as phylogenetic makers in many researches. The sequence variances represent the relationship of evolution among each plant. This research collected 335 samples from common plants including 221 species of 72 families in north Taiwan. The trnL intron and trnL-trnF IGS of each sample was sequenced to assess phylogenetic analysis for forensic applications. The 55 of 65 species for trnL intron and the 62 of 66 species for trnL-trnF IGS conformed to the species which have been known in the phylogenetic analysis. The remaining 4 species could not refer to the known groups by combining analyses of trnL intron and trnL-trnF IGS. Also 2 of 25 genuses could not entirely refer to the identified genuses in the phylogenetic analyses. The genetic distance and bootstrap values in phylogenetic analysis could refer to the species and genus of unknown sample. This study indicated that the trnL intron and trnL-trnF IGS could be used as efficient phylogenetic markers for forensic sciences.

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.