Part I of this thesis describes the development of bioactive small molecules of relevance to the study of the apicomlexan parasite Toxoplasma gondii into useful chemical tools. The research includes the target identification and validation studies, using both chemical and biological methods. Chapter 1 provides an overview of chemical genetics with a particular emphasis on methods for the identification of the protein targets of bioactive small molecules. The concept of biochemical protein target identification techniques was introduced with a detailed discussion of interesting applications from the literature. Chapter 2 focuses on the development of a tetrahydro-β-carboline based lead molecule into a chemical tool through target identification studies. The structure activity relationship (SAR) data associated with this core structure, the design of a chemical inducer of dimerisation (CID) and the synthesis of this CID are discussed in detail. Chapter 3 described work done to identify the potential protein target(s) of Conoidin A. Experiments to assess whether Conoidin A can inhibit a proposed target in vitro are also included. Further optimisation of this structural class to develop more potent inhibitors is discussed in the second part of this chapter. Part II of this thesis describes the development of methods for the synthesis of medium-sized ring containing compounds using oxidative fragmentation and rearrangement strategies. Chapter 5 provides an overview of the existing oxidative fragmentation methodology, with an emphasis on the use of oxidative fragmentation reactions for the synthesis of medium-sized ring systems (8-11 ring atoms). Chapter 6 focuses on using the established oxidative fragmentation method in the oxizino carbazolone system to investigate the diasteroselectivity of this reaction. Possible mechanisms for this transformation are investigated and discussed using both chemical and computational methods. An interesting rearrangement reaction has also been observed during this study. Chapter 7 focuses on developing an asymmetric oxidative fragmentation method, for use in the diazabenz[e]aceathrylenes system. Asymmetric oxidative fragmentation reactions using [Ru(pybox)(pydic)] catalysts are discussed. Attempts to optimise the enantiomeric excesses of the reaction by varying reaction conditions and substituents in the substrate are also included.

CoxsackievirusB3 (CVB3), a member of the Picornaviridae family is the causative agent of Virus-induced Myocarditis and Dilated Cardiomyopathy. The 5’UTR contains an Internal Ribosome Entry Site or IRES element that recruits ribosomes in a cap-independent manner. The ribosomes are recruited upstream of the AUG triplet at 591 (AUG591), also called as the cryptic AUG, after which they scan downstream for about 150 nucleotide, before initiating at the initiator AUG or AUG741. The 3’UTR of CVB3 is 99 nts long, highly structured RNA containing conserved domains, and is followed by a poly (A) tail of variable lengths. We have investigated possible involvement of host proteins which may interact with CVB3 IRES and influence its activity. We have demonstrated the role of Poly-pyrimidine tract binding protein (PTB) and established PTB as a bona-fide ITAF for CVB3, by characterizing the effect of partial silencing of PTB ex-vivo in HeLa cells. The IRES activity in BSC-1 cells, reported to have very low level of endogenous PTB, is found to be significantly low compared to that in HeLa cells. PTB is observed to interact with both the 5’ and 3’ UTR of CVB3, although with different affinities. Finer mapping of the interaction between PTB and the UTRs showed that the protein interacts with multiple regions of both UTRs. We have also shown the cis-acting effect of the CVB3-3’UTR on IRES mediated translation. The PTB contact points on the 3’UTRwas found to map to conserved regions, the deletion of which abrogates the 3’UTR mediated enhancement of the IRES activity. The possible role played by PTB in enhancing IRES activity by CVB3 3’UTR suggests that PTB protein might help in circularization of the CVB3 RNA by bridging the ends necessary for efficient translation of the viral RNA. In the second part, we have investigated possible role of some of the cis-acting element present in the CVB3 5’UTR RNA particularly the cryptic AUG. We have shown that mutation in cryptic AUG reduces the efficiency of translation mediated by the CVB3 IRES. Mutation in cryptic AUG moiety also reduces the interaction of mutant RNA with La protein. We have demonstrated that binding of 48S ribosomal complex with mutant IRES RNA was weaker compared to wt IRES RNA. We have investigated the possible alteration in secondary structure in the mutant RNA by chemical and enzymatic modification, which suggests that there is marginal alteration in the local structure due to mutation. It appears that integrity of cryptic AUG is important for efficient translation initiation by the CVB3 IRES. Results suggest that cryptic AUG plays a significant role in mediating internal initiation of translation of CVB3 RNA by mediating precise La binding and correct positioning of the 48S ribosomal complex. Finally, we have investigated the importance of a conserved hexa-nucleotide stretch in the apical loop within stem loop C (SLC, nt104-180), upstream of the ribosome landing site, on CVB3 IRES function. It has been already shown from our laboratory that the deletion at this apical loop resulted in significant decrease in IRES activity. This deletion mutant was shown to alter the secondary structure of the CVB3 5’UTR RNA. Here we have investigated the effect of point mutation in the apical loop SLC/c on CVB3 IRES activity by generating substitution mutation in the apical loop SLC/c in order to avoid possible alteration in secondary structure. Both the deletion or substitution mutation at this apical loop resulted in significant decrease in IRES activity. Both the mutant IRES RNAs (deletion and substitution mutant) failed to interact with certain trans-acting factors. Furthermore, expression of CVB3 2A protease significantly enhanced IRES activity of the wild type, but the effect was not so pronounced on the mutant IRESs. It is possible that the mutant RNAs were unable to interact with some trans-acting factors critical for enhanced IRES function. We have short-listed three proteins of approximate molecular mass of 56, 64 and 90 kDa, which showed reduced binding with mutant IRESs. By using RNA affinity column with biotinylated UTP labeled RNA we have purified couple of proteins and identified p64 as Cyto Keratin 1 protein by performing in-gel trypsin digestion followed by MALDI analysis. Overall, the results characterize the CVB3 IRES structurally and functionally, which could be useful in targeting critical RNA-protein interactions to develop candidate antiviral agent against Coxsackievirus infection.

No description available.

No description available.

Promoter region is the key regulatory region, which enables the gene to be transcribed or repressed by anchoring RNA polymerase and other transcription factors, but it is difficult to determine experimentally. Hence an in silico identification of promoters is crucial in order to guide experimental work and to pin point the key region that controls the transcription initiation of a gene. Analysis of various genome sequences in the vicinity of experimentally identified transcription start sites (TSSs) in prokaryotic as well as eukaryotic genomes had earlier indicated that they have several structural features in common, such as lower stability, higher curvature and less bendability, when compared with their neighboring regions. In this thesis work, the variation observed for these DNA sequence dependent structural properties have been used to identify and delineate promoter regions from other genomic regions. Since the number of bacterial genomes being sequenced is increasing very rapidly, it is crucial to have procedures for rapid and reliable annotation of their functional elements such as promoter regions, which control the expression of each gene or each transcription unit of the genome. The thesis work addresses this requirement and presents step by step protocols followed to get a generic method for promoter prediction that can be applicable across organisms. The each paragraph below gives an overall idea about the thesis organization into chapters. An overview of prokaryotic transcriptional regulation, structural polymorphism adapted by DNA molecule and its impact on transcriptional regulation has been discussed in introduction chapter of this thesis (chapter 1). Standardization of promoter prediction methodology - Part I Based on the difference in stability between neighboring upstream and downstream regions in the vicinity of experimentally determined transcription start sites, a promoter prediction algorithm has been developed to identify prokaryotic promoter sequences in whole genomes. The average free energy (E) over known promoter sequences and the difference (D) between E and the average free energy over the random sequence generated using the downstream region of known TSS (REav) are used to search for promoters in the genomic sequences. Using these cutoff values to predict promoter regions across entire E. coli genome, a reliability of 70% has been achieved, when the predicted promoters were cross verified against the 960 transcription start sites (TSSs) listed in the Ecocyc database. Reliable promoter prediction is obtained when these genome specific threshold values were used to search for promoters in the whole E. coli genome sequence. Annotation of the whole E. coli genome for promoter region has been carried out with 49% accuracy. Reference Rangannan, V. and Bansal, M. (2007) Identification and annotation of promoter regions inmicrobial genome sequences on the basis of DNA stability. J Biosci, 32, 851-862. Standardization of promoter prediction methodology - Part II In this chapter, it has been demonstrated that while the promoter regions are in general less stable than the flanking regions, their average free energy varies depending on the GC composition of the flanking genomic sequence. Therefore, a set of free energy threshold values (TSS based threshold values), from the genomic DNA with varying GC content in the vicinity of experimentally identified TSSs have been obtained. These threshold values have been used as generic criteria for predicting promoter regions in E. coli and B. subtilis and M. tuberculosis genomes, using an in-house developed tool ‘PromPredict’. On applying it to predict promoter regions corresponding to the 1144 and 612 experimentally validated TSSs in E. coli (genome %GC : 50.8) and B. subtilis (genome %GC : 43.5) sensitivity of 99% and 95% and precision values of 58% and 60%, respectively, were achieved. For the limited data set of 81 TSSs available for M. tuberculosis (65.6% GC) a sensitivity of 100% and precision of 49% was obtained. Reference Rangannan, V. and Bansal, M. (2009) Relative stability of DNA as a generic criterion for promoter prediction: whole genome annotation of microbial genomes with varying nucleotide base composition. Mol Biosyst, 5, 1758 - 1769. Standardization of promoter prediction methodology - Part III In this chapter, the promoter prediction algorithm and the threshold values have been improved to predict promoter regions on a large scale over 913 microbial genome sequences. The average free energy (AFE) values for the promoter regions as well as their downstream regions are found to differ, depending on their GC content even with respect to translation start sites (TLSs) from 913 microbial genomes. The TSS based cut-off values derived in chapter 3 do not have cut-off values for both extremes of GC-bins at 5% interval. Hence, threshold values have been derived from a subset of translation start sites (TLSs) from all microbial genomes which were categorized based on their GC-content. Interestingly the cut-off values derived with respect to TSS data set (chapter 3) and TLS data set are very similar for the in-between GC-bins. Therefore, TSS based cut-off values derived in chapter 2 with the TLS based cut-off values have been combined (denoted as TSS-TLS based cutoff values) to predict promoters over the complete genome sequences. An average recall value of 72% (which indicates the percentage of protein and RNA coding genes with predicted promoter regions assigned to them) and precision of 56% is achieved over the 913 microbial genome dataset. These predicted promoter regions have been given a reliability level (low, medium, high, very high and highest) based on the difference in its relative average free energy, which can help the users design their experiments with more confidence by using the predictions with higher reliability levels. Reference Rangannan, V. and Bansal, M. (2010) High Quality Annotation of Promoter Regions for 913 Bacterial Genomes. Bioinformatics, 26, 3043-3050. Web applications PromBase : The predicted promoter regions for 913 microbial genomes were deposited into a public domain database called, PromBase which can serve as a valuable resource for comparative genomics study for their general genomic features and also help the experimentalist to rapidly access the annotation of the promoter regions in any given genome. This database is freely accessible for the users via the World Wide Web http://nucleix.mbu.iisc.ernet.in/prombase/. EcoProm : EcoProm is a database that can identify and display the potential promoter regions corresponding to EcoCyc annotated TSS and genes. Also displays predictions for whole genomic sequence of E. coli and EcoProm is available at http://nucleix.mbu.iisc.ernet.in/ecoprom/index.htm. PromPredict : The generic promoter prediction methodology described in previous chapters has been implemented in to an algorithm ‘PromPredict’ and available at http://nucleix.mbu.iisc.ernet.in/prompredict/prompredict.html. Analysing the DNA structural characteristic of prokaryotic promoter sequences for their predominance Sequence dependent structural properties and their variation in genomic DNA are important in controlling several crucial processes such as transcription, replication, recombination and chromatin compaction. In this chapter 6, quantitative analysis of sequences motifs as well as sequence dependent structural properties, such as curvature, bendability and stability in the upstream region of TSS and TLS from E. coli, B. subtilis and M. tuberculosis has been carried out in order to assess their predictive power for promoter regions. Also the correlation between these structural properties and GC-content has been investigated. Our results have shown that AFE values (stability) gives finer discrimination rather than %GC in identifying promoter regions and stability have shown to be the better structural property in delineating promoter regions from non-promoter regions. Analysis of these DNA structural properties has been carried out in human promoter sequences and observed to be correlating with the inactivation status of the X-linked genes in human genome. Since, it is deviating from the theme of main thesis; this chapter has been included as appendix A to the main thesis. General conclusion Stability is the ubiquitous DNA structural property seen in promoter regions. Stability shows finer discrimination for promoter prediction rather than directly using %GC-content. Based on relative stability of DNA, a generic promoter prediction algorithm has been developed and implemented to predict promoter regions on a large scale over 913 microbial genome sequences. The analysis of the predicted regions across organisms showed highly reliable predictive performance of the algorithm.

This thesis work deals with different aspects of the chemistry of La(III) and Gd(III) complexes, their interaction with DNA and proteins, photo-induced cleavage of double-stranded DNA, photocytotoxic effect on cancer cells, cell death mechanism and cellular localization behaviour. Chapter I gives an introduction to the metal-based anticancer agents with special emphasis on clinically used drugs and the growing field of lanthanide therapeutics. An overview of the current strategies of cancer treatment, especially photodynamic therapy (PDT), is presented. Mode of small molecule-DNA interactions and the mechanistic aspects associated with DNA photodamage reactions and PDT effect are discussed with selected examples of compounds that are known to photocleave DNA on exposure to light of different wavelengths. A brief discussion on the various therapeutic applications of the lanthanide compounds is also made. Chapter II presents the synthesis, characterization, DNA binding, BSA binding, photo-induced DNA cleavage activity and photocytotoxicity of La(III) and Gd(III) complexes of phenanthroline bases to explore the UV-A light-induced DNA cleavage activity and photocytotoxicity of the complexes. Chapter III describes the synthesis, characterization, DNA binding, photo-induced DNA cleavage activity and photocytotoxicity of La(III) and Gd(III) complexes of phenanthroline bases with an aim to improve the design of the complexes to achieve better solution stability and DNA binding of the complexes. Chapter IV presents the synthesis, characterization, DNA binding, and UV-A light-induced DNA photocleavage activity and photocytotoxicity of La(III) and Gd(III) complexes of pyridyl phenanthroline bases with an objective to improve the photoactivity of the complexes by introducing an additional pyridyl group. Cell death mechanism and confocal microscopic studies are also carried out to gain more insight into the PDT effect caused by light in the presence of the complex. Chapter V describes the synthesis and characterization of La(III) and Gd(III) complexes of terpyridine bases and acetylacetonate to study the complexes as a new class of photosensitizers to explore their DNA photocleavage activity and photocytotoxicity in HeLa cells. Effect of attaching a glucose moiety to the acetyl acetone (Hacac) ligand has been studied. The cellular uptake behaviour of the La(III) pyrenyl-terpyridine complexes has also been investigated. Finally, Chapter VI presents the synthesis and characterization of curcumin and glycosylated curcumin La(III) and Gd(III) complexes having terpyridine base with an objective to study the photoactivated anticancer activity of the complexes in visible light. This chapter describes the visible light-induced DNA cleavage activity and photocytotoxicity of the complexes by exploiting curcumin and glycosylated curcumin as the photosensitizer ligands. Study on the cellular uptake behavior of curcumin La(III) complexes having pyrenyl terpyridine ligand is also presented. The references have been assembled at the end of each chapter and indicated as superscript numbers in the text. The complexes presented in this thesis are represented by bold-faced numbers. Crystallographic data of the complexes which are characterized structurally by single crystal X-ray crystallography are provided in CIF format in the enclosed CD (Appendix-I). Due acknowledgements have been made wherever the work described is based on the findings of other investigators. Any unintentional omission that might have happened due to oversight or mistake is sincerely regretted.

The flow of genetic information from genes to proteins is mediated through proteins which interact with the nucleic acids at several stages to successfully transmit the information from the nucleus to the cell cytoplasm. Unlike in the case of protein-protein interactions, the principles behind protein-nucleic acid interactions are still not very (Pabo and Nekludova, 2000) and efforts are still underway to arrive at the basic principles behind the specific recognition of nucleic acids by proteins (Prabakaran et al., 2006). This is mainly due to the innate complexity involved in recognition of nucleotides by proteins, where, even within a given family of DNA binding proteins, different modes of binding and recognition strategies are employed to suit their function (Luscomb et al., 2000). Such difficulties have also not made possible, a thorough classification of DNA/RNA binding proteins based on the mode of interaction as well as the specificity of recognition of the nucleotides. The availability of a large number of structures of protein-nucleic acids complexes (albeit lesser than the number of protein structures present in the PDB) in the past few decades has provided the knowledge-base for understanding the details behind their molecular mechanisms (Berman et al., 1992). Previously, studies have been carried out to characterize these interactions by analyzing specific non-covalent interactions such as hydrogen bonds, van der Walls, and hydrophobic interactions between a given amino acid and the nucleic acid (DNA, RNA) in a pair-wise manner, or through the analysis of interface areas of the protein-nucleic acid complexes (Nadassy et al., 1998; Jones et al., 1999). Though the studies have deciphered the common pairing preferences of a particular amino acid with a given nucleotide of DNA or RNA, there is little room for understanding these specificities in the context of spatial interactions at a global level from the protein-nucleic acid complexes. The representation of the amino acids and the nucleotides as components of graphs, and trying to explore the nature of the interactions at a level higher than exploring the individual pair-wise interactions, could provide greater details about the nature of these interactions and their specificity. This thesis reports the study of protein-nucleic interactions using graph and network based approaches. The evaluation of the parameters for characterizing protein-nucleic acid graphs have been carried out for the first time and these parameters have been successfully employed to capture biologically important non-covalent interactions as clusters of interacting amino acids and nucleotides from different protein-DNA and protein-RNA complexes. Graph and network based approaches are well established in the field of protein structure analysis for analyzing protein structure, stability and function (Kannan and Vishveshwara, 1999; Brinda and Vishveshwara, 2005). However, the use of graph and network principles for analyzing structures of protein-nucleic acid complexes is so far not accomplished and is being reported the first time in this thesis. The matter embodied in the thesis is presented as ten chapters. Chapter 1 lays the foundation for the study, surveying relevant literature from the field. Chapter 2 describes in detail the methods used in constructing graphs and networks from protein-nucleic acid complexes. Initially, only protein structure graphs and networks are constructed from proteins known to interact with specific DNA or RNA, and inferences with regard to nucleic acid binding and recognition were indirectly obtained . Subsequently, parameters were evaluated for representing both the interacting amino acids and the nucleotides as components of graphs and a direct evaluation of protein-DNA and Protein-RNA interactions as graphs has been carried out. Chapter 3 and 4 discuss the graph and network approaches applied to proteins from a dataset of DNA binding proteins complexed with DNA. In chapter 3, the protein structure graphs were constructed on the basis of the non-covalent interactions existing between the side chains of amino acids. Clusters of interacting side chains from the graphs were obtained using the graph spectral method. The clusters from the protein-DNA interface were analyzed in detail for the interaction geometry and biological importance (Sathyapriya and Vishveshwara, 2004). Chapter 4 also uses the same dataset of DNA binding proteins, but a network-based approach is presented. From the analysis of the protein structure networks from these DNA binding proteins, interesting observations relating the presence of highly connected nodes(or hubs) of the network to functionally important amino acids in the structure, emerged. Also, the comparison between the hubs identified from the protein-protein and the protein-DNA interfaces in terms of their amino acid composition and their connectivity are also presented (Sathyapriya and Vishveshwara, 2006) Chapter 5 and 6 deal with the graph and network applications to a specific system of protein-RNA complex (aminoacyl-tRNA synthetases) to gain insights into their interface biology based on amino acid connectivity. Chapter 5 deals with a dataset of aminoacyl-tRNA synthetase (aaRS) complexes obtained with various ligands like ATP, tRNA and L-amino acids. A graph based identification of side chain clusters from these ligand-bound aaRS structures has highlighted important features of ligand-binding at the catalytic sites of the two structurally different classes of aaRS (Class I and Class II). Side chain clusters from other regions of aaRS such as the anticodon binding region and the ligand-activation sites are discussed. A network approach is used in a specific system of aaRS(E.coli Glutaminyl-tRNA synthetase (GlnRS) complexed with its ligands, to specifically understand the effects of different ligand binding., in chapter 6. The structure networks of E.coli GlnRS in the ligand-free and different ligand-bound states are constructed. The ligand-free and the ligand-bound complexes are compared by analyzing their network properties and the presence of hubs to understand the effect of ligand-binding. These properties have elegantly captured the effects of ligand-binding to the GlnRS structure and have also provided an alternate method for comparing three dimensional structures of proteins in different ligand-bound states (Sathyapriya and Vishveshwara, 2007). In contrast to protein structure graphs (PSG), both the interacting amino acids and nucleotides (DNA/RNA) form the components of the protein-nucleic acid graphs (PNG) from protein-nucleic acid complexes. These graphs are constructed based on the non-covalent interactions existing between the side chains of the amino acids and nucleotides. After representing the interacting nucleotides and amino acids as graphs, clusters of the interacting components are identified. These clusters are the strongly interacting amino acids and nucleotides from the protein-nucleic acid complexes. These clusters can be generated at different strengths of interaction between the amino acid side chain and the nucleotide (measured in terms of its atomic connectivity) and can be used for detecting clusters of non-specific as well as specific interactions of amino acids and nucleotides. Though the methodology of graph construction and cluster identification are given in chapter 2, the details of the parameters evaluated for constructing PNG are given in chapter 7. Unlike in the previous chapters, the succeeding chapters deal exclusively with results that are obtained from the analyses of PNG. Two examples of obtaining clusters from a PNG are given, one each for a protein-DNA and a protein-RNA complex. In the first example, a nucleosome core particle is subjected to the graph based analysis and different clusters of amino acids with different regions of the DNA chain such as phosphate, deoxyribose sugar and the base are identified. Another example of aminoacyl-tRNA synthetase complexed with its cognate tRNA is used to illustrate the method with a protein-RNA complex. Further, the method of constructing and analyzing protein-nucleic acid graphs has been applied to the macromolecular machinery of the pre-translocation complex of the T. thermophilus 70S ribosome. Chapter 8 deals exclusively with the results identified from the analysis of this magnificent macromolecular ensemble. The availability of the method that can handle interactions between both amino acids and the nucleotides of the protein-nucleic acid complexes has given us the basis fro evaluating these interactions in a level higher than that of analyzing pair-wise interactions. A study on the evaluation of short hydrogen bonds(SHB) in proteins, which does not fall under the realm of the main objective of the thesis, is discussed in the Chapter 9. The short hydrogen bonds, defined by the geometrical distance and angle parameters, are identified from a non-redundant dataset of proteins. The insights into their occurrence, amino acid composition and secondary structural preferences are discussed. The SHB are present in distinct regions of protein three-dimensional structures, such that they mediate specific geometrical constraints that are necessary for stability of the structure (Sathyapriya and Vishveshwara, 2005). The significant conclusions of various studies carried out are summarized in the last chapter (Chapter 10). In conclusion, this thesis reports the analyses performed with protein-nucleic acid complexes using graph and network based methods. The parameters necessary for representing both amino acids and the nucleotides as components of a graph, are evaluated for the first time and can be used subsequently for other analyses. More importantly, the use of graph-based methods has resulted in considering the interaction between the amino acids and the nucleotides at a global level with respect to their topology of the protein-nucleic acid complexes. Such studies performed on a wide variety of protein-nucleic acid complexes could provide more insights into the details of protein-nucleic acid recognition mechanisms. The results of these studies can be used for rational design of experimental mutations that ascertain the structure-function relationships in proteins and protein-nucleic acid complexes.

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.

La plupart des conditions détectées par le dépistage néonatal sont reliées à l'une des enzymes qui dégradent les acyls-CoA mitochondriaux. Le rôle physiopathologique des acyls-CoA dans ces maladies est peu connue, en partie parce que les esters liés au CoA sont intracellulaires et les échantillons tissulaires de patients humains ne sont généralement pas disponibles. Nous avons créé une modèle animal murin de l'une de ces maladies, la déficience en 3-hydroxy-3-methylglutaryl-CoA lyase (HL), dans le foie (souris HLLKO). HL est la dernière enzyme de la cétogenèse et de la dégradation de la leucine. Une déficience chronique en HL et les crises métaboliques aigües, produisent chacune un portrait anormal et distinct d'acyls-CoA hépatiques. Ces profils ne sont pas prévisibles à partir des niveaux d'acides organiques urinaires et d'acylcarnitines plasmatiques. La cétogenèse est indétectable dans les hépatocytes HLLKO. Dans les mitochondries HLLKO isolées, le dégagement de 14CO2 à partir du [2-14C]pyruvate a diminué en présence de 2-ketoisocaproate (KIC), un métabolite de la leucine. Au test de tolérance au pyruvate, une mesure de la gluconéogenèse, les souris HLLKO ne présentent pas la réponse hyperglycémique normale. L'hyperammoniémie et l'hypoglycémie, des signes classiques de plusieurs erreurs innées du métabolisme (EIM) des acyls-CoA, surviennent de façon spontanée chez des souris HLLKO et sont inductibles par l'administration de KIC. Une charge en KIC augmente le niveau d'acyls-CoA reliés à la leucine et diminue le niveau d'acétyl-CoA. Les mitochondries des hépatocytes des souris HLLKO traitées avec KIC présentent un gonflement marqué. L'hyperammoniémie des souris HLLKO répond au traitement par l'acide N-carbamyl-L-glutamique. Ce composé permet de contourner une enzyme acétyl-CoA-dépendante essentielle pour l’uréogenèse, le N-acétylglutamate synthase. Ceci démontre un mécanisme d’hyperammoniémie lié aux acyls-CoA. Dans une deuxième EIM des acyls-CoA, la souris SCADD, déficiente en déshydrogénase des acyls-CoA à chaînes courtes. Le profil des acyls-CoA hépatiques montre un niveau élevé du butyryl-CoA particulièrement après un jeûne et après une charge en triglycérides à chaîne moyenne précurseurs du butyryl-CoA. / Most conditions detected by expanded newborn screening result from deficiency of one of the enzymes that degrade acyl-CoA esters in mitochondria. The role of acyl-CoAs in the pathophysiology of these disorders is poorly understood, in part because CoA esters are intracellular and samples are not generally available from human patients. We created a mouse model of one such condition, deficiency of 3-hydroxy-3-methylglutaryl-CoA lyase (HL), in liver (HLLKO mice). HL catalyses a reaction of ketone body synthesis and of leucine degradation. Chronic HL deficiency and acute crises each produced distinct abnormal liver acyl-CoA patterns, which would not be predictable from levels of urine organic acids and plasma acylcarnitines. In HLLKO hepatocytes, ketogenesis was undetectable. Measures of Krebs cycle flux diminished following incubation of HLLKO mitochondria with the leucine metabolite 2-ketoisocaproate (KIC). HLLKO mice also had suppression of the normal hyperglycemic response to a systemic pyruvate load, a measure of gluconeogenesis. Hyperammonemia and hypoglycemia, cardinal features of many inborn errors of acyl-CoA metabolism, occurred spontaneously in some HLLKO mice and were inducible by administering KIC. KIC loading also increased levels of several leucine-related acyl-CoAs and reduced acetyl-CoA levels. Ultrastructurally, hepatocyte mitochondria of KIC-treated HLLKO mice show marked swelling. KIC-induced hyperammonemia improved following administration of carglumate (N-carbamyl-L-glutamic acid), which bypasses an acetyl-CoA-dependent reaction essential for urea cycle function, thus demonstrating an acyl-CoA-related mechanism for this complication. In a second animal model of an inborn error of acyl-CoA metabolism, short chain acyl-CoA dehydrogenase (SCAD)-deficient mice, the main finding in liver acyl-CoAs is increased butyryl-CoA, particularly during fasting or after enteral loading with medium chain triglyceride precursor of butyryl-CoA.

La plupart des conditions détectées par le dépistage néonatal sont reliées à l'une des enzymes qui dégradent les acyls-CoA mitochondriaux. Le rôle physiopathologique des acyls-CoA dans ces maladies est peu connue, en partie parce que les esters liés au CoA sont intracellulaires et les échantillons tissulaires de patients humains ne sont généralement pas disponibles. Nous avons créé une modèle animal murin de l'une de ces maladies, la déficience en 3-hydroxy-3-methylglutaryl-CoA lyase (HL), dans le foie (souris HLLKO). HL est la dernière enzyme de la cétogenèse et de la dégradation de la leucine. Une déficience chronique en HL et les crises métaboliques aigües, produisent chacune un portrait anormal et distinct d'acyls-CoA hépatiques. Ces profils ne sont pas prévisibles à partir des niveaux d'acides organiques urinaires et d'acylcarnitines plasmatiques. La cétogenèse est indétectable dans les hépatocytes HLLKO. Dans les mitochondries HLLKO isolées, le dégagement de 14CO2 à partir du [2-14C]pyruvate a diminué en présence de 2-ketoisocaproate (KIC), un métabolite de la leucine. Au test de tolérance au pyruvate, une mesure de la gluconéogenèse, les souris HLLKO ne présentent pas la réponse hyperglycémique normale. L'hyperammoniémie et l'hypoglycémie, des signes classiques de plusieurs erreurs innées du métabolisme (EIM) des acyls-CoA, surviennent de façon spontanée chez des souris HLLKO et sont inductibles par l'administration de KIC. Une charge en KIC augmente le niveau d'acyls-CoA reliés à la leucine et diminue le niveau d'acétyl-CoA. Les mitochondries des hépatocytes des souris HLLKO traitées avec KIC présentent un gonflement marqué. L'hyperammoniémie des souris HLLKO répond au traitement par l'acide N-carbamyl-L-glutamique. Ce composé permet de contourner une enzyme acétyl-CoA-dépendante essentielle pour l’uréogenèse, le N-acétylglutamate synthase. Ceci démontre un mécanisme d’hyperammoniémie lié aux acyls-CoA. Dans une deuxième EIM des acyls-CoA, la souris SCADD, déficiente en déshydrogénase des acyls-CoA à chaînes courtes. Le profil des acyls-CoA hépatiques montre un niveau élevé du butyryl-CoA particulièrement après un jeûne et après une charge en triglycérides à chaîne moyenne précurseurs du butyryl-CoA. / Most conditions detected by expanded newborn screening result from deficiency of one of the enzymes that degrade acyl-CoA esters in mitochondria. The role of acyl-CoAs in the pathophysiology of these disorders is poorly understood, in part because CoA esters are intracellular and samples are not generally available from human patients. We created a mouse model of one such condition, deficiency of 3-hydroxy-3-methylglutaryl-CoA lyase (HL), in liver (HLLKO mice). HL catalyses a reaction of ketone body synthesis and of leucine degradation. Chronic HL deficiency and acute crises each produced distinct abnormal liver acyl-CoA patterns, which would not be predictable from levels of urine organic acids and plasma acylcarnitines. In HLLKO hepatocytes, ketogenesis was undetectable. Measures of Krebs cycle flux diminished following incubation of HLLKO mitochondria with the leucine metabolite 2-ketoisocaproate (KIC). HLLKO mice also had suppression of the normal hyperglycemic response to a systemic pyruvate load, a measure of gluconeogenesis. Hyperammonemia and hypoglycemia, cardinal features of many inborn errors of acyl-CoA metabolism, occurred spontaneously in some HLLKO mice and were inducible by administering KIC. KIC loading also increased levels of several leucine-related acyl-CoAs and reduced acetyl-CoA levels. Ultrastructurally, hepatocyte mitochondria of KIC-treated HLLKO mice show marked swelling. KIC-induced hyperammonemia improved following administration of carglumate (N-carbamyl-L-glutamic acid), which bypasses an acetyl-CoA-dependent reaction essential for urea cycle function, thus demonstrating an acyl-CoA-related mechanism for this complication. In a second animal model of an inborn error of acyl-CoA metabolism, short chain acyl-CoA dehydrogenase (SCAD)-deficient mice, the main finding in liver acyl-CoAs is increased butyryl-CoA, particularly during fasting or after enteral loading with medium chain triglyceride precursor of butyryl-CoA.