Probing the Allosteric Mechanism in Pyrrolysyl-tRNA Synthetase Using Energy-Weighted Network Formalism

Pyrrolysyl-tRNA synthetase (PyIRS) is an atypical enzyme responsible for charging tRNA(Pyl) with pyrrolysine, despite lacking precise tRNA anticodon recognition. This dimeric protein exhibits allosteric regulation of function, like any other tRNA synthetases. In this study we examine the paths of allosteric communication at the atomic level, through energy-weighted networks of Desulfitobacterium hafniense PyIRS (DhPyIRS) and its complexes with tRNA(Pyl) and activated pyrrolysine. We performed molecular dynamics simulations of the structures of these complexes to obtain an ensemble conformation-population perspective. Weighted graph parameters relevant to identifying key players and ties in the context of social networks such as edge/node betweenness, closeness index, and the concept of funneling are explored in identifying key residues and interactions leading to shortest paths of communication in the structure networks of DhPylRS. Further, the changes in the status of important residues and connections and the costs of communication due to ligand induced perturbations are evaluated. The optimal, suboptimal, and preexisting paths are also investigated. Many of these parameters have exhibited an enhanced asymmetry between the two subunits of the dimeric protein, especially in the pretransfer complex, leading us to conclude that encoding of function goes beyond the sequence/structure of proteins. The local and global perturbations mediated by appropriate ligands and their influence on the equilibrium ensemble of conformations also have a significant role to play in the functioning of proteins. Taking a comprehensive view of these observations, we propose that the origin of many functional aspects (allostery rand half-sites reactivity in the case of DhPyIRS) lies in subtle rearrangements of interactions and dynamics at a global level.

Allosteric Communication in Cysteinyl tRNA Synthetase A NETWORK OF DIRECT AND INDIRECT READOUT

Protein structure networks are constructed for the identification of long-range signaling pathways in cysteinyl tRNA synthetase (CysRS). Molecular dynamics simulation trajectory of CysRS-ligand complexes were used to determine conformational ensembles in order to gain insight into the allosteric signaling paths. Communication paths between the anticodon binding region and the aminoacylation region have been identified. Extensive interaction between the helix bundle domain and the anticodon binding domain, resulting in structural rigidity in the presence of tRNA, has been detected. Based on the predicted model, six residues along the communication paths have been examined by mutations (single and double) and shown to mediate a coordinated coupling between anticodon recognition and activation of amino acid at the active site. This study on CysRS clearly shows that specific key residues, which are involved in communication between distal sites in allosteric proteins but may be elusive in direct structure analysis, can be identified from dynamics of protein structure networks.

Structures of new crystal forms of Mycobacterium tuberculosis peptidyl-tRNA hydrolase and functionally important plasticity of the molecule

The X-ray structures of new crystal forms of peptidyl-tRNA hydrolase from M.similar to tuberculosis reported here and the results of previous X-ray studies of the enzyme from different sources provide a picture of the functionally relevant plasticity of the protein molecule. The new X-ray results confirm the connection deduced previously between the closure of the lid at the peptide-binding site and the opening of the gate that separates the peptide-binding and tRNA-binding sites. The plasticity of the molecule indicated by X-ray structures is in general agreement with that deduced from the available solution NMR results. The correlation between the lid and the gate movements is not, however, observed in the NMR structure.

Catalysis of tRNA Aminoacylation: Single Turnover to Steady-State Kinetics of tRNA Synthetases

Aminoacyl-tRNA synthetases (aaRS) catalyze the bimolecular association reaction between amino acid and tRNA by specifically and unerringly choosing the cognate amino acid and tRNA. There are two classes of such synthetases that perform tRNA-aminoacylation reaction. Interestingly, these two classes of aminoacyl-tRNA synthetases differ not only in their structures but they also exhibit remarkably distinct kinetics under pre-steady-state condition. The class I synthetases show initial burst of product formation followed by a slower steady-state rate. This has been argued to represent the influence of slow product release. In contrast, there is no burst in the case of class H enzymes. The tight binding of product with enzyme for class I enzymes is correlated with the enhancement of rate in presence of elongation factor. EF-TU. In spite of extensive experimental studies, there is no detailed theoretical analysis that can provide a quantitative understanding of this important problem. In this article, we present a theoretical investigation of enzyme kinetics for both classes of aminoacyl-tRNA synthetases. We present an augmented kinetic scheme and then employ the methods of time-dependent probability statistics to obtain expressions for the first passage time distribution that gives both the time-dependent and the steady-state rates. The present study quantitatively explains all the above experimental observations. We propose an alternative path way in the case of class II enzymes showing the tRNA-dependent amino acid activation and the discrepancy between the single-turnover and steady-state rate.

Kinetic proofreading at single molecular level: aminoacylation of tRNA(Ile) and the role of water as an editor

Proofreading/editing in protein synthesis is essential for accurate translation of information from the genetic code. In this article we present a theoretical investigation of efficiency of a kinetic proofreading mechanism that employs hydrolysis of the wrong substrate as the discriminatory step in enzyme catalytic reactions. We consider aminoacylation of tRNA(Ile) which is a crucial step in protein synthesis and for which experimental results are now available. We present an augmented kinetic scheme and then employ methods of stochastic simulation algorithm to obtain time dependent concentrations of different substances involved in the reaction and their rates of formation. We obtain the rates of product formation and ATP hydrolysis for both correct and wrong substrates (isoleucine and valine in our case, respectively), in single molecular enzyme as well as ensemble enzyme kinetics. The present theoretical scheme correctly reproduces (i) the amplitude of the discrimination factor in the overall rates between isoleucine and valine which is obtained as (1.8x10(2)).(4.33x10(2)) = 7.8x10(4), (ii) the rates of ATP hydrolysis for both Ile and Val at different substrate concentrations in the aminoacylation of tRNA(Ile). The present study shows a non-michaelis type dependence of rate of reaction on tRNA(Ile) concentration in case of valine. The overall editing in steady state is found to be independent of amino acid concentration. Interestingly, the computed ATP hydrolysis rate for valine at high substrate concentration is same as the rate of formation of Ile-tRNA(Ile) whereas at intermediate substrate concentration the ATP hydrolysis rate is relatively low. We find that the presence of additional editing domain in class I editing enzyme makes the kinetic proofreading more efficient through enhanced hydrolysis of wrong product at the editing CP1 domain.

How Many Initiator tRNA Genes Does Escherichia coli Need?

Multiple copies of a gene require enhanced investment on the part of the cell and, as such, call for an explanation. The observation that Escherichia coli has four copies of initiator tRNA (tRNA(i)) genes, encoding a special tRNA (tRNA(fMet)) required to start protein synthesis, is puzzling particularly because the cell appears to be unaffected by the removal of one copy. However, the fitness of an organism has both absolute and relative connotations. Thus, we carried out growth competition experiments between E. coli strains that differ in the number of tRNA(i) genes they contain. This has enabled us to uncover an unexpected link between the number of tRNA(i) genes and protein synthesis, nutritional status, and fitness. Wild-type strains with the canonical four tRNA(i) genes are favored in nutrient-rich environments, and those carrying fewer are favored in nutrient-poor environments. Auxotrophs behave as if they have a nutritionally poor internal environment. A heuristic model that links tRNA(i) gene copy number, genetic stress, and growth rate accounts for the findings. Our observations provide strong evidence that natural selection can work through seemingly minor quantitative variations in gene copy number and thereby impact organismal fitness.

Insights into the Structural Dynamics of Nucleocytoplasmic Transport of tRNA by Exportin-t

Exportin-t (Xpot) transports mature 5'- and 3'-end processed tRNA from the nucleus to the cytoplasm by associating with a small G-protein Ran (RAs-related nuclear protein), in the nucleus. The release of tRNA in cytoplasm involves RanGTP hydrolysis. Despite the availability of crystal structures of nuclear and cytosolic forms of Xpot, the molecular details regarding the sequential events leading to tRNA release and subsequent conformational changes occurring in Xpot remain unknown. We have performed a combination of classical all-atom and accelerated molecular dynamics simulations on a set of complexes involving Xpot to study a range of features including conformational flexibility of free and cargo-bound Xpot and functionally critical contacts between Xpot and its cargo. The systems investigated include free Xpot and its different complexes, bound either to Ran (GTP/GDP) or tRNA or both. This approach provided a statistically reliable estimate of structural dynamics of Xpot after cargo release. The mechanistic basis for Xpot opening after cargo release has been explained in terms of dynamic structural hinges, about which neighboring region could be displaced to facilitate the nuclear to cytosolic state transition. Post-RanGTP hydrolysis, a cascade of events including local conformational change in RanGTP and loss of critical contacts at Xpot/tRNA interface suggest factors responsible for eventual release of tRNA. The level of flexibility in different Xpot complexes varied depending on the arrangement of individual HEAT repeats. Current study provides one of the most comprehensive and robust analysis carried out on this protein using molecular dynamics schemes.

An in vivo approach to tRNA identity

A leucine-inserting tRNA has been transformed into a serine-inserting tRNA by changing 12 nucleotides. Only 8 of the 12 changes are required to effect the conversion of the leucine tRNA to serine tRNA identity. The 8 essential changes reside in basepair 11-24 in the D stem, basepairs 3-70, 2-71 and nucleotides 72 and 73, all of the acceptor stem.

Functional amber suppressor tRNA genes were generated for 14 species of tRNA in E. coli, and their amino acid specificities determined. The suppressors can be classified into three groups, based upon their specificities. Class I suppressors, tRNA^(Ala2)_(CUA), tRNA^(GlyU)_(CUA), tRNA^(HisA)_(CUA), tRNA^(Lys)_(CUA), and tRNA^(ProH)_(CUA), inserted the predicted amino acid. The Class II suppressors, tRNA^(GluA)_(CUA) , tRNA^(GlyT)_(CUA), and tRNA^(Ile1)_(CUA) were either partially or predominantly mischarged by the glutamine aminoacyl tRNA synthetase (AAS). The Class III suppressors, tRNA^(Arg)_(CUA), tRNA^(AspM)_(CUA), tRNA^(Ile2)_(CUA), tRNA^(Thr2)_(CUA), tRNA^(Met(m))_(CUA) and tRNA^(Val)_(CUA) inserted predominantly lysine.

Studies of the serotonin type 3A receptor and the chemical preparation of tRNA

This thesis describes studies surrounding a ligand-gated ion channel (LGIC): the serotonin type 3A receptor (5-HT₃AR). Structure-function experiments using unnatural amino acid mutagenesis are described, as well as experiments on the methodology of unnatural amino acid mutagenesis. Chapter 1 introduces LGICs, experimental methods, and an overview of the unnatural amino acid mutagenesis.

In Chapter 2, the binding orientation of the clinically available drugs ondansetron and granisetron within 5-HT₃A is determined through a combination of unnatural amino acid mutagenesis and an inhibition based assay. A cation-π interaction is found for both ondansetron and granisetron with a specific tryptophan residue (Trp183, TrpB) of the mouse 5-HT₃AR, which establishes a binding orientation for these drugs.

In Chapter 3, further studies were performed with ondansetron and granisetron with 5-HT₃A. The primary determinant of binding for these drugs was determined to not include interactions with a specific tyrosine residue (Tyr234, TyrC2). In completing these studies, evidence supporting a cation-π interaction of a synthetic agonist, meta-chlorophenylbiguanide, was found with TyrC2.

In Chapter 4, a direct chemical acylation strategy was implemented to prepare full-length suppressor tRNA mediated by lanthanum(III) and amino acid phosphate esters. The derived aminoacyl-tRNA is shown to be translationally competent in Xenopus oocytes.

Appendix A.1 gives details of a pharmacological method for determining the equilibrium dissociation constant, K_B, of a competitive antagonist with a receptor, known as Schild analysis. Appendix A.2 describes an examination of the inhibitory activity of new chemical analogs of the 5-HT₃A antagonist ondansetron. Appendix A.3 reports an organic synthesis of an intermediate for a new unnatural amino acid. Appendix A.4 covers an additional methodological examination for the preparation of amino-acyl tRNA.

Mitochondrial gene sequences useful for species identification of western North Atlantic Ocean sharks

Molecular-based approaches for shark species identification have been driven largely by issues specific to the fishery. In an effort to establish a more comprehensive identification data set, we investigated DNA sequence variation of a 1.4-kb region from the mitochondrial genome covering partial sequences from the 12S rDNA, 16S rDNA, and the complete valine tRNA from 35 shark species from the Atlantic fishery. Generally, within-species variability was low in relation to interspecific divergence because species haloptypes formed monophyletic groups. Phylogenetic analyses resolved ordinal relationships among Carcharhiniformes and Lamniformes, and revealed support for the families Sphyrnidae and Triakidae (within Carcharhiniformes) and Lamnidae and Alopidae (within Lamniformes). The combination of limited intraspecific variability and sufficient between-species divergence indicates that this locus is suitable for species identification.

大银鱼和小齿日本银鱼线粒体细胞色素b 和16S rRNA基因部分序列分析

　对大银鱼和小齿日本银鱼的线粒体细胞色素b 和16S rRNA 基因片段进行了扩增和序列测定,分析比较了2 种 间的序列差异。大银鱼2 个个体间细胞色素b 基因序列无差异,小齿日本银鱼3 个个体的序列出现了2 种单倍型;大银鱼 同小齿日本银鱼细胞色素b 序列(单倍型SB21) 之间存在86 个碱基差异(差异率21 %) ,变异较大。大银鱼、小齿日本银鱼 的16S rRNA 基因片段种内个体间无序列差异,两种间存在21 个碱基的差异(差异率5 %) ,有1 个碱基的插入/ 缺失。由 此可见,2 种银鱼间线粒体细胞色素b 基因的进化速率较快,约为16S rRNA 基因片段的4 倍。根据细胞色素b 序列数据, 推算出2 种银鱼大概在中新世晚期发生分化。