Communication pathways from the anti-codon region to the aminoacylation site in Methionyl tRNA Synthetase: Molecular dynamics simulations and the structure network analysis

The enzymes of the family of tRNA synthetases perform their functions with high precision by synchronously recognizing the anticodon region and the aminoacylation region, which are separated by ?70 in space. This precision in function is brought about by establishing good communication paths between the two regions. We have modeled the structure of the complex consisting of Escherichia coli methionyl-tRNA synthetase (MetRS), tRNA, and the 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 have been evaluated. Furthermore, the network analysis on these simulated structures has been carried out to elucidate the paths of communication between the activation site and the anticodon recognition site. This study has provided the detailed paths of communication, which are consistent with experimental results. Similar studies also have been carried out on the complexes (MetRS + activated methonine) and (MetRS + tRNA) along with ligand-free native enzyme. A comparison of the paths derived from the four simulations clearly has 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 details of the method of our investigation and the biological implications of the results are presented in this article. The method developed here also could be used to investigate any protein system where the function takes place through long-distance communication.

WebGeSTer DB-a transcription terminator database

We present WebGeSTer DB, the largest database of intrinsic transcription terminators (http://pallab.serc.iisc.ernet.in/gester). The database comprises of a million terminators identified in 1060 bacterial genome sequences and 798 plasmids. Users can obtain both graphic and tabular results on putative terminators based on default or user-defined parameters. The results are arranged in different tiers to facilitate retrieval, as per the specific requirements. An interactive map has been incorporated to visualize the distribution of terminators across the whole genome. Analysis of the results, both at the whole-genome level and with respect to terminators downstream of specific genes, offers insight into the prevalence of canonical and non-canonical terminators across different phyla. The data in the database reinforce the paradigm that intrinsic termination is a conserved and efficient regulatory mechanism in bacteria. Our database is freely accessible.

Functional Significance of an Evolutionarily Conserved Alanine (GCA) Resume Codon in tmRNA in Escherichia coli

Occasionally, ribosomes stall on mRNAs prior to the completion of the polypeptide chain. In Escherichia coli and other eubacteria, tmRNA-mediated trans-translation is a major mechanism that recycles the stalled ribosomes. The tmRNA possesses a tRNA-like domain and a short mRNA region encoding a short peptide (ANDENYALAA in E. coli) followed by a termination codon. The first amino acid (Ala) of this peptide encoded by the resume codon (GCN) is highly conserved in tmRNAs in different species. However, reasons for the high evolutionary conservation of the resume codon identity have remained unclear. In this study, we show that changing the E. coli tmRNA resume codon to other efficiently translatable codons retains efficient functioning of the tmRNA. However, when the resume codon was replaced with the low-usage codons, its function was adversely affected. Interestingly, expression of tRNAs decoding the low-usage codon from plasmid-borne gene copies restored efficient utilization of tmRNA. We discuss why in E. coli, the GCA (Ala) is one of the best codons and why all codons in the short mRNA of the tmRNA are decoded by the abundant tRNAs.

Use of a coupled transcriptional system for consistent overexpression and purification of UDG-Ugi complex and Ugi from Escherichia coli

We have designed a novel coupled transcriptional construct wherein Escherichia coil uracil DNA glycosylase (UDC:) and Bacillus subtilis phage PBS-2 encoded uracil DNA glycosylase inhibitor protein (Ugi) genes were cloned in tandem, downstream of an inducible promoter (P-trc). Use of this bicistronic operon has allowed purification of large amounts of UDG-Ugi complex formed in vivo. The system has also been exploited for purification of large amounts of Ugi. While establishing the expression system, one of the constructs showed detectable suppression of UAG termination codon and resulted in accumulation of a minor population of a putative readthrough polypeptide cor responding to UDG. We discuss the likely occurrence of such a phenomenon in overproduction of other recombinant proteins. Finally, the usefulness of the operon construct in convenient mutational analysis to study the mechanism of UDG-Ugi interaction is also discussed.

Allostery in tRNA Synthetases Elucidated from MD Simulations and Protein Structure Networks

tRNA synthetases (aaRS) are enzymes crucial in the translation of genetic code. The enzyme accylates the acceptor stem of tRNA by the congnate amino acid bound at the active site, when the anti-codon is recognized by the anti-codon site of aaRS. In a typical aaRS, the distance between the anti-codon region and the amino accylation site is approximately 70 Å. We have investigated this allosteric phenomenon at molecular level by MD simulations followed by the analysis of protein structure networks (PSN) of non-covalent interactions. Specifically, we have generated conformational ensembles by performing MD simulations on different liganded states of methionyl tRNA synthetase (MetRS) from Escherichia coli and tryptophenyl tRNA synthetase (TrpRS) from Human. The correlated residues during the MD simulations are identified by cross correlation maps. We have identified the amino acids connecting the correlated residues by the shortest path between the two selected members of the PSN. The frequencies of paths have been evaluated from the MD snapshots[1]. The conformational populations in different liganded states of the protein have been beautifully captured in terms of network parameters such as hubs, cliques and communities[2]. These parameters have been associated with the rigidity and plasticity of the protein conformations and can be associated with free energy landscape. A comparison of allosteric communication in MetRS and TrpRS [3] elucidated in this study highlights diverse means adopted by different enzymes to perform a similar function. The computational method described for these two enzymes can be applied to the investigation of allostery in other systems.

p53 and little brother p53/47: linking IRES activities with protein functions

The tumor suppressor p53 represents a paradigm for gene regulation. Its rapid induction in response to DNA damage conditions has been attributed to both increased half-life of p53 protein and also increased translation of p53 mRNA. Recent advances in our understanding of the post-transcriptional regulation of p53 include the discovery of internal ribosome entry sites (IRESs) within the p53 mRNA. These IRES elements regulate the translation of the full length as well as the N-terminally truncated isoform, p53/47. The p53/47 isoform is generated by alternative initiation at an internal AUG codon present within the p53 ORF. The aim of this review is to summarize the role of translational control mechanisms in regulating p53 functions. We discuss here in detail how diverse cellular stress pathways trigger alterations in the cap-dependent and cap-independent translation of p53 mRNA and how changes in the relative expression levels of p53 isoforms result in more differentiated p53 activity.

Occurrence, divergence and evolution of intrinsic terminators across Eubacteria

In Escherichia coli, the canonical intrinsic terminator of transcription includes a palindrome followed by a U-trail on the transcript. The apparent underrepresentation of such terminators in eubacterial genomes led us to develop a rapid and accurate algorithm, GeSTer, to predict putative intrinsic terminators. Now, we have analyzed 378 genome sequences with an improved version of GeSTer. Our results indicate that the canonical E. coli type terminators are not overwhelmingly abundant in eubacteria. The atypical structures, having stem-loop structures but lacking ‘U’ trail, occur downstream of genes in all the analyzed genomes but different phyla show conserved preference for different types of terminators. This propensity correlates with genomic GC content and presence of the factor, Rho. 60–70% of identified terminators in all the genomes show “optimized” stem-length and ΔG. These results provide evidence that eubacteria extensively rely on the mechanism of intrinsic termination, with a considerable divergence in their structure, positioning and prevalence. The software and detailed results for individual genomes are freely available on request

Implications of high level pseudogene transcription in Mycobacterium leprae

Background: The Mycobacterium leprae genome has less than 50% coding capacity and 1,133 pseudogenes. Preliminary evidence suggests that some pseudogenes are expressed. Therefore, defining pseudogene transcriptional and translational potentials of this genome should increase our understanding of their impact on M. leprae physiology. Results: Gene expression analysis identified transcripts from 49% of all M. leprae genes including 57% of all ORFs and 43% of all pseudogenes in the genome. Transcribed pseudogenes were randomly distributed throughout the chromosome. Factors resulting in pseudogene transcription included: 1) co-orientation of transcribed pseudogenes with transcribed ORFs within or exclusive of operon-like structures; 2) the paucity of intrinsic stem-loop transcriptional terminators between transcribed ORFs and downstream pseudogenes; and 3) predicted pseudogene promoters. Mechanisms for translational ``silencing'' of pseudogene transcripts included the lack of both translational start codons and strong Shine-Dalgarno (SD) sequences. Transcribed pseudogenes also contained multiple ``in-frame'' stop codons and high Ka/Ks ratios, compared to that of homologs in M. tuberculosis and ORFs in M. leprae. A pseudogene transcript containing an active promoter, strong SD site, a start codon, but containing two in frame stop codons yielded a protein product when expressed in E. coli. Conclusion: Approximately half of M. leprae's transcriptome consists of inactive gene products consuming energy and resources without potential benefit to M. leprae. Presently it is unclear what additional detrimental affect(s) this large number of inactive mRNAs has on the functional capability of this organism. Translation of these pseudogenes may play an important role in overall energy consumption and resultant pathophysiological characteristics of M. leprae. However, this study also demonstrated that multiple translational ``silencing'' mechanisms are present, reducing additional energy and resource expenditure required for protein production from the vast majority of these transcripts.

Cloning, overexpression, folding and purification of a biosynthetically derived three disulfide scorpion toxin (BTK-2) from Mesobuthus tamulus

BTK-2, a 32 residue scorpion toxin initially identified in the venom of red Indian scorpion Mesobuthus tamulus was cloned, overexpressed and purified using Cytochrome 155 fusion protein system developed in our laboratory. The synthetic gene coding for the peptide was designed taking into account optimal codon usage by Escherichia coli. High expression levels of the fusion protein enabled facile purification of this peptide. The presence of disulfide bonded isomers, occurring as distinctly populated states even in the fusion protein, were separated by gel filtration chromatography. The target peptide was liberated from the host protein by Tev protease cleavage and subsequent purification was achieved using RP-HPLC methods. Reverse phase HPLC clearly showed the presence of at least two isomeric forms of the peptide that were significantly populated. The oxidative folding of BTK-2 was achieved under ambient conditions during the course of purification. Structural characterization of the two forms, by solution homonuclear and heteronuclear NMR methods, has shown that these two forms exhibit significantly different structural properties, and represent the natively folded and a "misfolded" form of the peptide. The formation of properly folded BTK-2 as a major fraction without the use of in vitro oxidative refolding methods clearly indicate the versatility of the Cytochrome b(5) fusion protein system for the efficient production of peptides for high resolution NMR studies.

Predictive testing for familial adenomatous polyposis in a rural South Indian community

During the course of genome studies in a rural community in the South Indian state of Karnataka, DNA-based investigations and counselling for familial adenomatous polyposis (FAP) were requested via the community physician. The proposita died in 1940 and FAP had been clinically diagnosed in 2 of her 5 children, both deceased. DNA samples from 2 affected individuals in the third generation were screened for mutations in the APC gene, and a frame-shift mutation was identified in exon 15 with a common deletion at codon 1061. Predictive testing for the mutation was then organized on a voluntary basis. There were 11 positive tests, including confirmatory positives on 2 persons diagnosed by colonoscopy, and to date surgery has been successfully undertaken on 3 previously undiagnosed adults. The ongoing success of the study indicates that, with appropriate access to the facilities offered by collaborating centres, predictive testing is feasible for diseases such as FAP and could be of significant benefit to communities in economically less developed countries.

Analysis of the beta-glucoside utilization (bgl) genes of Shigella sonnei: evolutionary implications for their maintenance in a cryptic state

The pattern of expression of the genes involved in the utilization of aryl beta-glucosides such as arbutin and salicin is different in the genus Shigella compared to Escherichia coli. The results presented here indicate that the homologue of the cryptic bgl operon of E. coli is conserved in Shigella sonnei and is the primary system involved in beta-glucoside utilization in the organism. The organization of the bgl genes in 5. sonnei is similar to that of E. coli; however there are three major differences in terms of their pattern of expression. (i) The bglB gene, encoding phospho-beta-glucosidase B, is insertionally inactivated in 5. sonnei. As a result, mutational activation of the silent bgl promoter confers an Arbutin-positive (Arb(+)) phenotype to the cells in a single step; however, acquiring a Salicin-positive (Sal(+)) phenotype requires the reversion or suppression of the bglB mutation in addition. (ii) Unlike in E. coli, a majority of the activating mutations (conferring the Arb(+) phenotype) map within the unlinked hns locus, whereas activation of the E. coli bgl operon under the same conditions is predominantly due to insertions within the bglR locus. (iii) Although the bgl promoter is silent in the wild-type strain of 5. sonnei (as in the case of E. coli), transcriptional and functional analyses indicated a higher basal level of transcription of the downstream genes. This was correlated with a 1 bp deletion within the putative Rho-independent terminator present in the leader sequence preceding the homologue of the bglG gene. The possible evolutionary implications of these differences for the maintenance of the genes in the cryptic state are discussed.

The fate of the initiator tRNAs is sensitive to the critical balance between interacting proteins

Formylation of the initiator tRNA is essential for normal growth of Escherichia coil, The initiator tRNA containing the U35A36 mutation (CUA anticodon) initiates from UAG codon, However, an additional mutation at position 72 (72A --> G) renders the tRNA (G72/U35A36) inactive in initiation because it is defective in formylation, In this study, we isolated U1G72/U35A36 tRNA containing a wobble base pair at 1-72 positions as an intragenic suppressor of the G72 mutation. The U1G72/U35A36 tRNA is formylated and participates in initiation. More importantly, we show that the mismatch at 1-72 positions of the initiator tRNA, which was thus far thought to be the hallmark of the resistance of this tRNA against peptidyl-tRNA hydrolase (PTH), is not sufficient, The amino acid attached to the initiator tRNA is also important in conferring protection against PTH. Further, we show that the relative levels of PTH and IF2 influence the path adopted by the initiator tRNAs in protein synthesis. These findings provide an important clue to understand the dual function of the single tRNA(Met) in initiation and elongation, in the mitochondria of various organisms.

sRNAscanner: A Computational Tool for Intergenic Small RNA Detection in Bacterial Genomes

Background:Bacterial non-coding small RNAs (sRNAs) have attracted considerable attention due to their ubiquitous nature and contribution to numerous cellular processes including survival, adaptation and pathogenesis. Existing computational approaches for identifying bacterial sRNAs demonstrate varying levels of success and there remains considerable room for improvement. Methodology/Principal Findings: Here we have proposed a transcriptional signal-based computational method to identify intergenic sRNA transcriptional units (TUs) in completely sequenced bacterial genomes. Our sRNAscanner tool uses position weight matrices derived from experimentally defined E. coli K-12 MG1655 sRNA promoter and rho-independent terminator signals to identify intergenic sRNA TUs through sliding window based genome scans. Analysis of genomes representative of twelve species suggested that sRNAscanner demonstrated equivalent sensitivity to sRNAPredict2, the best performing bioinformatics tool available presently. However, each algorithm yielded substantial numbers of known and uncharacterized hits that were unique to one or the other tool only. sRNAscanner identified 118 novel putative intergenic sRNA genes in Salmonella enterica Typhimurium LT2, none of which were flagged by sRNAPredict2. Candidate sRNA locations were compared with available deep sequencing libraries derived from Hfq-co-immunoprecipitated RNA purified from a second Typhimurium strain (Sittka et al. (2008) PLoS Genetics 4: e1000163). Sixteen potential novel sRNAs computationally predicted and detected in deep sequencing libraries were selected for experimental validation by Northern analysis using total RNA isolated from bacteria grown under eleven different growth conditions. RNA bands of expected sizes were detected in Northern blots for six of the examined candidates. Furthermore, the 5'-ends of these six Northern-supported sRNA candidates were successfully mapped using 5'-RACE analysis. Conclusions/Significance: We have developed, computationally examined and experimentally validated the sRNAscanner algorithm. Data derived from this study has successfully identified six novel S. Typhimurium sRNA genes. In addition, the computational specificity analysis we have undertaken suggests that similar to 40% of sRNAscanner hits with high cumulative sum of scores represent genuine, undiscovered sRNA genes. Collectively, these data strongly support the utility of sRNAscanner and offer a glimpse of its potential to reveal large numbers of sRNA genes that have to date defied identification. sRNAscanner is available from: http://bicmku.in:8081/sRNAscanner or http://cluster.physics.iisc.ernet.in/sRNAscanner/.

A novel approach to design of cis-acting DNA structural element for regulation of gene expression in vivo

Taking advantage of the degeneracy of the genetic code we have developed a novel approach to introduce, within a gene, DNA sequences capable of adopting unusual structures and to investigate the role of such sequences in regulation of gene expression in vivo. We used a computer program that generates alternative codon sequences for the same amino-acid sequence to convert a stretch of nucleotides into an inverted-repeat sequence with the potential to adopt cruciform structure. This approach was used to replace a 51-base-pair EcoRI-HindIII segment in the N-terminal region of the beta-galactosidase gene in plasmid pUC19 with a 51-bp synthetic oligonucleotide sequence with the potential to adopt a cruciform structure with 18 bp in the stem region. In selecting the 51-bp sequence, care was taken to include those codons that are preferred in E. coli. E. coli DH5-alpha cells harbouring the plasmid containing the redesigned sequence showed drastic reduction in expression of the beta-galactosidase gene compared to cells harbouring the plasmid with the native sequence. This approach demonstrates the possibility of introducing DNA secondary-structure elements to alter regulation of gene expression in vivo.

Identification and characterization of a new plasmid carrying genes for degradation of 2,4-dichlorophenoxyacetate from Pseudomonas cepacia CSV90

Pseudomonas cepacia CSV90 is able to utilize 2,4-dichlorophenoxyacetate (2,4-D) and 2-methyl-4-chlorophenoxyacetate as sole sources of carbon and energy. Mutants of the strain CSV90 which had lost this ability appeared spontaneously on a nonselective medium. The wild-type strain harbored a 90-kb plasmid, pMAB1, whereas 2,4-D-negative mutants either lost the plasmid or had a 70-kb plasmid, pMAB2. The plasmid pMAB2 was found to have undergone a deletion Of a 20-kb fragment of pMAB1. The plasmid-free mutants regained the ability to degrade 2,4-D after introduction of purified pMAB1 by electroporation. Cloning in Escherichia coli of a 10-kb BamHI fragment from pMAB1, the region absent in pMAB2, resulted in the expression of the gene tfdC encoding 3,5-dichlorocatechol 1,2-dioxygenase. After subcloning, the tfdC gene was located in a 1.6-kb HindIII fragment. The nucleotide sequence of the tfdC gene and the restriction map of its contiguous region are identical to those of the well-characterized 2,4-D-degradative plasmid pJP4 of Alcaligenes eutrophus, whereas the overall restriction maps of the two plasmids are different. The N-terminal 44-amino-acid sequence of the enzyme purified from the strain CSV90 confirmed the reading frame in the DNA sequence for tfdC and indicated that the initiation codon GUG is read as methionine instead of valine.