Rare codons cluster

The characteristics of rare codon clusters in the genome and proteins of hepatitis C virus; a bioinformatics look

BACKGROUND Recent studies suggest that rare codon clusters are functionally important for protein activity. METHODS Here, for the first time we analyzed and reported rare codon clusters in Hepatitis C Virus (HCV) genome and then identified the location of these rare codon clusters in the structure of HCV protein. This analysis was performed using the Sherlocc program that detects statistically relevant conserved rare codon clusters. RESULTS By this program, we identified the rare codon cluster in three regions of HCV genome; NS2, NS3, and NS5A coding sequence of HCV genome. For further understanding of the role of these rare codon clusters, we studied the location of these rare codon clusters and critical residues in the structure of NS2, NS3 and NS5A proteins. We identified some critical residues near or within rare codon clusters. It should be mentioned that characteristics of these critical residues such as location and situation of side chains are important in assurance of the HCV life cycle. CONCLUSION The characteristics of these residues and their relative status showed that these rare codon clusters play an important role in proper folding of these proteins. Thus, it is likely that these rare codon clusters may have an important role in the function of HCV proteins. This information is helpful in development of new avenues for vaccine and treatment protocols.

The effects of the context-dependent codon usage bias on the structure of the nsp1α of porcine reproductive and respiratory syndrome virus

The information about the crystal structure of porcine reproductive and respiratory syndrome virus (PRRSV) leader protease nsp1α is available to analyze the roles of tRNA abundance of pigs and codon usage of the nsp1 α gene in the formation of this protease. The effects of tRNA abundance of the pigs and the synonymous codon usage and the context-dependent codon bias (CDCB) of the nsp1 α on shaping the specific folding units (α-helix, β-strand, and the coil) in the nsp1α were analyzed based on the structural information about this protease from protein data bank (PDB: 3IFU) and the nsp1 α of the 191 PRRSV strains. By mapping the overall tRNA abundance along the nsp1 α, we found that there is no link between the fluctuation of the overall tRNA abundance and the specific folding units in the nsp1α, and the low translation speed of ribosome caused by the tRNA abundance exists in the nsp1 α. The strong correlation between some synonymous codon usage and the specific folding units in the nsp1α was found, and the phenomenon of CDCB exists in the specific folding units of the nsp1α. These findings provide an insight into the roles of the synonymous codon usage and CDCB in the formation of PRRSV nsp1α structure.

Guanine nucleotide-binding protein (Gα) endocytosis by a cascade of ubiquitin binding domain proteins is required for sustained morphogenesis and proper mating in yeast

Heterotrimeric G proteins are well known to transmit signals from cell surface receptors to intracellular effector proteins. There is growing appreciation that G proteins are also present at endomembrane compartments, where they can potentially interact with a distinct set of signaling proteins. Here, we examine the cellular trafficking function of the G protein α subunit in yeast, Gpa1. Gpa1 contains a unique 109-amino acid insert within the α-helical domain that undergoes a variety of posttranslational modifications. Among these is monoubiquitination, catalyzed by the NEDD4 family ubiquitin ligase Rsp5. Using a newly optimized method for G protein purification together with biophysical measures of structure and function, we show that the ubiquitination domain does not influence enzyme activity. By screening a panel of 39 gene deletion mutants, each lacking a different ubiquitin binding domain protein, we identify seven that are necessary to deliver Gpa1 to the vacuole compartment including four proteins (Ede1, Bul1, Ddi1, and Rup1) previously not known to be involved in this process. Finally, we show that proper endocytosis of the G protein is needed for sustained cellular morphogenesis and mating in response to pheromone stimulation. We conclude that a cascade of ubiquitin-binding proteins serves to deliver the G protein to its final destination within the cell. In this instance and in contrast to the previously characterized visual system, endocytosis from the plasma membrane is needed for proper signal transduction rather than for signal desensitization.

Understanding the influence of codon translation rates on cotranslational protein folding

Protein domains can fold into stable tertiary structures while they are synthesized by the ribosome in a process known as cotranslational folding. If a protein does not fold cotranslationally, however, it has the opportunity to do so post-translationally, that is, after the nascent chain has been fully synthesized and released from the ribosome. The rate at which a ribosome adds an amino acid encoded by a particular codon to the elongating nascent chain can vary significantly and is called the codon translation rate. Recent experiments have illustrated the profound impact that codon translation rates can have on the cotranslational folding process and the acquisition of function by nascent proteins. Synonymous codon mutations in an mRNA molecule change the chemical identity of a codon and its translation rate without changing the sequence of the synthesized protein. This change in codon translation rate can, however, cause a nascent protein to malfunction as a result of cotranslational misfolding. In some situations, such dysfunction can have profound implications; for example, it can alter the substrate specificity of an ABC transporter protein, resulting in patients who are nonresponsive to chemotherapy treatment. Thus, codon translation rates are crucial in coordinating protein folding in a cellular environment and can affect downstream cellular processes that depend on the proper functioning of newly synthesized proteins. As the importance of codon translation rates makes clear, a necessary aspect of fully understanding cotranslational folding lies in considering the kinetics of the process in addition to its thermodynamics. In this Account, we examine the contributions that have been made to elucidating the mechanisms of cotranslational folding by using the theoretical and computational tools of chemical kinetics, molecular simulations, and systems biology. These efforts have extended our ability to understand, model, and predict the influence of codon translation rates on cotranslational protein folding and misfolding. The application of such approaches to this important problem is creating a framework for making quantitative predictions of the impact of synonymous codon substitutions on cotranslational folding that has led to a novel hypothesis regarding the role of fast-translating codons in coordinating cotranslational folding. In addition, it is providing new insights into proteome-wide cotranslational folding behavior and making it possible to identify potential molecular mechanisms by which molecular chaperones can influence such behavior during protein synthesis. As we discuss in this Account, bringing together these theoretical developments with experimental approaches is increasingly helping answer fundamental questions about the nature of nascent protein folding on the ribosome.

An optimized expression vector for improving the yield of dengue virus-like particles from transfected insect cells

Recombinant virus-like particles (rVLPs) of flaviviruses are non-infectious particles released from cells expressing the envelope glycoproteins prM and E. Dengue virus rVLPs are recognized as a potential vaccine candidate, but large scale production of these particles is hindered by low yields and the occurrence of cytopathic effects. In an approach to improve the yield of rVLPs from transfected insect cells, several components of a dengue serotype 2 virus prM+E expression cassette were modified and the effect of these modifications was assessed during transient expression. Enhancement of extracellular rVLP levels by simultaneous substitutions of the prM signal peptide and the stem-anchor region of E with homologous cellular and viral counterparts, respectively, was further augmented by codon optimization. Extensive formation of multinucleated cells following transfection with the codon-optimized expression cassette was abrogated by introducing an E fusion loop mutation. This mutation also helped restore the extracellular E levels affected negatively by alteration of a charged residue at the pr-M junction, which was intended to promote maturation of rVLPs during export. Optimized expression cassettes generated in this multiple add-on modification approach should be useful in the generation of stably expressing clones and production of dengue virus rVLPs for immunogenicity studies.

Expanding Anfinsen's principle: contributions of synonymous codon selection to rational protein design

Anfinsen’s principle asserts that all information required to specify the structure of a protein is encoded in its amino acid sequence. However, during protein synthesis by the ribosome, the N-terminus of the nascent chain can begin to fold before the C-terminus is available. We tested whether this cotranslational folding can alter the folded structure of an encoded protein in vivo, versus the structure formed when refolded in vitro. We designed a fluorescent protein consisting of three half-domains, where the N- and C-terminal half-domains compete with each other to interact with the central half-domain. The outcome of this competition determines the fluorescence properties of the resulting folded structure. Upon refolding after chemical denaturation, this protein produced equimolar amounts of the N- and C-terminal folded structures, respectively. In contrast, translation in Escherichia coli resulted in a 2-fold enhancement in the formation of the N-terminal folded structure. Rare synonymous codon substitutions at the 5′ end of the C-terminal half-domain further increased selection for the N-terminal folded structure. These results demonstrate that the rate at which a nascent protein emerges from the ribosome can specify the folded structure of a protein.

Co-translational mechanisms of protein maturation

Protein biogenesis integrates multiple finely regulated mechanisms, ensuring nascent polypeptide chains are correctly enzymatically processed, targeted to membranes and folded to native structure. Recent studies show that the cellular translation machinery serves as hub that coordinates the maturation events in space and time at various levels. The ribosome itself serves as docking site for a multitude of nascent chain-interacting factors. The movement of ribosomes along open reading frames is non-uniformous and includes pausing sites, which facilitates nascent chain folding and perhaps factor engagement. Here we summarize current knowledge and discuss emerging concepts underlying the critical interplay between translation and protein maturation in E. coli.

Molecular chaperone functions in protein folding and proteostasis

The biological functions of proteins are governed by their three-dimensional fold. Protein folding, maintenance of proteome integrity, and protein homeostasis (proteostasis) critically depend on a complex network of molecular chaperones. Disruption of proteostasis is implicated in aging and the pathogenesis of numerous degenerative diseases. In the cytosol, different classes of molecular chaperones cooperate in evolutionarily conserved folding pathways. Nascent polypeptides interact cotranslationally with a first set of chaperones, including trigger factor and the Hsp70 system, which prevent premature (mis)folding. Folding occurs upon controlled release of newly synthesized proteins from these factors or after transfer to downstream chaperones such as the chaperonins. Chaperonins are large, cylindrical complexes that provide a central compartment for a single protein chain to fold unimpaired by aggregation. This review focuses on recent advances in understanding the mechanisms of chaperone action in promoting and regulating protein folding and on the pathological consequences of protein misfolding and aggregation.

The distribution of synonymous codon choice in the translation initiation region of dengue virus

Dengue is the most common arthropod-borne viral (Arboviral) illness in humans. The genetic features concerning the codon usage of dengue virus (DENV) were analyzed by the relative synonymous codon usage, the effective number of codons and the codon adaptation index. The evolutionary distance between DENV and the natural hosts (Homo sapiens, Pan troglodytes, Aedes albopictus and Aedes aegypti) was estimated by a novel formula. Finally, the synonymous codon usage preference for the translation initiation region of this virus was also analyzed. The result indicates that the general trend of the 59 synonymous codon usage of the four genotypes of DENV are similar to each other, and this pattern has no link with the geographic distribution of the virus. The effect of codon usage pattern of Aedes albopictus and Aedes aegypti on the formation of codon usage of DENV is stronger than that of the two primates. Turning to the codon usage preference of the translation initiation region of this virus, some codons pairing to low tRNA copy numbers in the two primates have a stronger tendency to exist in the translation initiation region than those in the open reading frame of DENV. Although DENV, like other RNA viruses, has a high mutation to adapt its hosts, the regulatory features about the synonymous codon usage have been 'branded' on the translation initiation region of this virus in order to hijack the translational mechanisms of the hosts.

Quantifying elongation rhythm during full-length protein synthesis

Pauses regulate the rhythm of ribosomal protein synthesis. Mutations disrupting even minor pauses can give rise to improperly formed proteins and human disease. Such minor pauses are difficult to characterize by ensemble methods, but can be readily examined by single-molecule (sm) approaches. Here we use smFRET to carry out real-time monitoring of the expression of a full-length protein, the green fluorescent protein variant Emerald GFP. We demonstrate significant correlations between measured elongation rates and codon and isoacceptor tRNA usage, and provide a quantitative estimate of the effect on elongation rate of replacing a codon recognizing an abundant tRNA with a synonymous codon cognate to a rarer tRNA. Our results suggest that tRNA selection plays an important general role in modulating the rates and rhythms of protein synthesis, potentially influencing simultaneous co-translational processes such as folding and chemical modification.

The effects of the synonymous codon usage and tRNA abundance on protein folding of the 3C protease of foot-and-mouth disease virus

The 3C protease of foot-and-mouth disease virus (FMDV) has a conserved amino acid sequence and is responsible for most cleavage in the viral polyprotein. The effects of the synonymous codon usage of FMDV 3C gene and tRNA abundance of the hosts on shaping different folding units (α-helix, β-strand and the coil) in the 3C protease were analyzed based on the structural information of the FMDV 3C protease from Protein Data Bank (PDB: 2BHG) and 210 genes of 3C for all serotypes of FMDV. The strong correlation between some codons usage and the specific folding unit in the FMDV 3C protease is found. As for the effect of translation speed caused by tRNA abundance on the formation of folding units, the codon positions with lowly abundant tRNA scatter in the 3C gene and there is the obvious fluctuation of tRNA abundance locating in the transition boundaries from the β-strand to the α-helix and the coil, respectively. However, codon positions with lowly abundant tRNA clustering into these boundaries are not found, suggesting that the adjustment of translation speed is likely also achieved by the single codon position with low tRNA abundance rather than a cluster. The observations can provide the information for insight into the role of the synonymous codon usage in the formation of 3C protease of FMDV and effect of the tRNA abundance of the hosts on this formation of protease.

The ribosome as a hub for protein quality control

Cells face a constant challenge as they produce new proteins. The newly synthesized polypeptides must be folded properly to avoid aggregation. If proteins do misfold, they must be cleared to maintain a functional and healthy proteome. Recent work is revealing the complex mechanisms that work cotranslationally to ensure protein quality control during biogenesis at the ribosome. Indeed, the ribosome is emerging as a central hub in coordinating these processes, particularly in sensing the nature of the nascent protein chain, recruiting protein folding and translocation components, and integrating mRNA and nascent chain quality control. The tiered and complementary nature of these decision-making processes confers robustness and fidelity to protein homeostasis during protein synthesis.

CBDB: the codon bias database

Background

In many genomes, a clear preference in the usage of particular codons exists. The mechanisms that induce codon biases remain an open question; studies have attributed codon usage to translational selection, mutational bias and drift. Furthermore, correlations between codon usage within host genomes and their viral pathogens have been observed for a myriad of host-virus systems. As such, numerous studies have investigated codon usage and codon bias in an effort to better understand how species evolve. Numerous metrics have been developed to identify biases in codon usage. In addition, a few data repositories of codon bias data are available, differing in the metrics reported as well as the number and taxonomy of strains examined.

Description

We have created a new web resource called the Codon Bias Database (CBDB) which provides information regarding the codon bias within the set of highly expressed genes for 300+ bacterial genomes. CBDB was developed to provide a resource for researchers investigating codon bias in bacteria, facilitating comparisons between strains and species. Furthermore, the site was created to serve those studying adaptation in phage; the genera selected for this first release of CBDB all have sequenced, annotated bacteriophages. The annotations and sequences for the highly expressed gene set are available for each strain in addition to the strain’s codon bias measurements.

Conclusions

Comparing species and strains provides a comprehensive look at how codon usage has been shaped over evolutionary time and can elucidate the putative mechanisms behind it. The Codon Bias Database provides a centralized repository of look-up tables and codon usage bias measures for a wide variety of genera, species and strains. Through our analysis of the variation in codon usage within the strains presently available, we find that most members of a genus have a codon composition most similar to other members of its genus, although not necessarily other members of its species.

Effects of rare codon clusters on the expression of a high-turnover chloroplast protein in Chlamydomonas reinhardtii

Expression of foreign proteins in chloroplasts has become an important field of plant genetic engineering. Optimized codon usage is generally thought to increase translational efficiency, but high speed translation of codon bias-adjusted mRNAs can also result in protein misfolding due to a lack of rare codons. In order to analyze the effect of rare codons on a native chloroplast protein in vivo, we modified the D1 subunit of photosystem II by fusing small peptides with different codons into a loop region which tolerates insertions without loss of function. Because of its high-turnover properties, the D1 protein represents an excellent test object to investigate the impact of rare codons on its translation. We choose codons for amino acids Arg, Leu, Ser, Ala and Gly which are rarely used and compared translation of the modified D1 proteins with the respective mutant proteins containing insertions with frequently used codons. Our data indicate that only rare Arg codons drastically affect synthesis of the D1 protein and cluster of rare Ser-codon can induce strategic ribosomal pausing sites.

Large-scale analysis of conserved rare codon clusters suggests an involvement in co-translational molecular recognition events

MOTIVATION

An increasing amount of evidence from experimental and computational analysis suggests that rare codon clusters are functionally important for protein activity. Most of the studies on rare codon clusters were performed on a limited number of proteins or protein families. In the present study, we present the Sherlocc program and how it can be used for large scale protein family analysis of evolutionarily conserved rare codon clusters and their relation to protein function and structure. This large-scale analysis was performed using the whole Pfam database covering over 70% of the known protein sequence universe. Our program Sherlocc, detects statistically relevant conserved rare codon clusters and produces a user-friendly HTML output.

RESULTS

Statistically significant rare codon clusters were detected in a multitude of Pfam protein families. The most statistically significant rare codon clusters were predominantly identified in N-terminal Pfam families. Many of the longest rare codon clusters are found in membrane-related proteins which are required to interact with other proteins as part of their function, for example in targeting or insertion. We identified some cases where rare codon clusters can play a regulating role in the folding of catalytically important domains. Our results support the existence of a widespread functional role for rare codon clusters across species. Finally, we developed an online filter-based search interface that provides access to Sherlocc results for all Pfam families.

AVAILABILITY

The Sherlocc program and search interface are open access and are available at http://bcb.med.usherbrooke.ca

Genetic code redundancy and its influence on the encoded polypeptides

The genetic code is said to be redundant in that the same amino acid residue can be encoded by multiple, so-called synonymous, codons. If all properties of synonymous codons were entirely equivalent, one would expect that they would be equally distributed along protein coding sequences. However, many studies over the last three decades have demonstrated that their distribution is not entirely random. It has been postulated that certain codons may be translated by the ribosome faster than others and thus their non-random distribution dictates how fast the ribosome moves along particular segments of the mRNA. The reasons behind such segmental variability in the rates of protein synthesis, and thus polypeptide emergence from the ribosome, have been explored by theoretical and experimental approaches. Predictions of the relative rates at which particular codons are translated and their impact on the nascent chain have not arrived at unequivocal conclusions. This is probably due, at least in part, to variation in the basis for classification of codons as “fast” or “slow”, as well as variability in the number and types of genes and proteins analyzed. Recent methodological advances have allowed nucleotide-resolution studies of ribosome residency times in entire transcriptomes, which confirm the non-uniform movement of ribosomes along mRNAs and shed light on the actual determinants of rate control. Moreover, experiments have begun to emerge that systematically examine the influence of variations in ribosomal movement and the fate of the emerging polypeptide chain.

A new helper phage for improved monovalent display of Fab molecules

Phage display technology is a powerful tool for the identification of novel antibodies for drug discovery. Phage display libraries have been constructed with massive diversity, but their use may be hindered by limited antibody display levels when rescued with the M13KO7 helper phage. Variants of M13KO7 have been constructed previously that increase the levels of display of rescued phage, but all produce phage that display multiple copies of the antibody fragment on their surface and have reduced titer and infectivity. In this study, we describe a new helper phage, XP5, which increased the display level of Fab molecules more than two-fold compared to phage rescued with M13KO7. XP5 uses a combination of ribosome binding site spacing alterations and rare codon clusters to reduce the expression of pIII from the helper phage. This reduction in pIII expression leads to an increase in the incorporation of pIII-Fab fusions during phage rescue. The rescued phage displayed a single copy of the Fab molecule, preventing any avidity effects during the selection process. This also suggests that the percentage of the population of phage displaying a Fab molecule is increased when rescued with XP5. Additionally, the phage titers and infectivity are comparable to libraries rescued with M13KO7. After two rounds of panning we observed a nearly 5-fold increase in the number of antigen binding Fab molecules compared to panning conducted with the same library rescued with M13KO7. The nature of the mutations in XP5 makes it a universal substitute for M13KO7 in pIII-based phage display, compatible with most phagemids and bacterial strains.

LaTcOm: a web server for visualizing rare codon clusters in coding sequences

We present LaTcOm, a new web tool, which offers several alternative methods for 'rare codon cluster' (RCC) identification from a single and simple graphical user interface. In the current version, three RCC detection schemes are implemented: the recently described %MinMax algorithm and a simplified sliding window approach, along with a novel modification of a linear-time algorithm for the detection of maximally scoring subsequences tailored to the RCC detection problem. Among a number of user tunable parameters, several codon-based scales relevant for RCC detection are available, including tRNA abundance values from Escherichia coli and several codon usage tables from a selection of genomes. Furthermore, useful scale transformations may be performed upon user request (e.g. linear, sigmoid). Users may choose to visualize RCC positions within the submitted sequences either with graphical representations or in textual form for further processing.

AVAILABILITY

LaTcOm is freely available online at the URL http://troodos.biol.ucy.ac.cy/latcom.html.

Manipulating the genetic code for membrane protein production: what have we learnt so far?

With synthetic gene services, molecular cloning is as easy as ordering a pizza. However choosing the right RNA code for efficient protein production is less straightforward, more akin to deciding on the pizza toppings. The possibility to choose synonymous codons in the gene sequence has ignited a discussion that dates back 50 years: Does synonymous codon use matter? Recent studies indicate that replacement of particular codons for synonymous codons can improve expression in homologous or heterologous hosts, however it is not always successful. Furthermore it is increasingly apparent that membrane protein biogenesis can be codon-sensitive. Single synonymous codon substitutions can influence mRNA stability, mRNA structure, translational initiation, translational elongation and even protein folding. Synonymous codon substitutions therefore need to be carefully evaluated when membrane proteins are engineered for higher production levels and further studies are needed to fully understand how to select the codons that are optimal for higher production. This article is part of a Special Issue entitled: Protein Folding in Membranes.

The codon usage model of the context flanking each cleavage site in the polyprotein of foot-and-mouth disease virus

To investigate the codon usage pattern of the contexts flanking 11 cleavage sites of foot-and-mouth disease virus (FMDV) polyprotein, the codon usage model of the corresponding codon position and the synonymous codon usage in the target contexts of 66 strains were characterized by two simple methods based on the relative synonymous codon usage value. The synonymous codons usage pattern was also compared between this virus and two species of hosts (cattle and domestic pig). It is indicated that FMDV bore a general resemblance to the hosts in terms of the synonymous codon usage pattern. This feature may help FMDV to utilize translational resources of host efficiently. The two amino acid residues constituting each cleavage site contain at least one conserved residue. It was noticed that the codon usage model with the strong bias appeared in some specific positions in the target contexts, and the under-represented synonymous codons, AUA for Ile, CUA for Leu, UUA for Leu and GUA for Val, are preferentially used in these positions. These under-represented synonymous codons likely play role in regulating the translation rate and influencing the secondary structure of the contexts flanking the cleavage sites.

Birth, life and death of nascent polypeptide chains

The journey of nascent polypeptides from synthesis at the peptidyl transferase center of the ribosome (“birth”) to full function (“maturity”) involves multiple interactions, constraints, modifications and folding events. Each step of this journey impacts the ultimate expression level and functional capacity of the translated protein. It has become clear that the kinetics of protein translation is predominantly modulated by synonymous codon usage along the mRNA, and that this provides an active mechanism for coordinating the synthesis, maturation and folding of nascent polypeptides. Multiple quality control systems ensure that proteins achieve their native, functional form. Unproductive co-translational folding intermediates that arise during protein synthesis may undergo enhanced interaction with components of these systems, such as chaperones, and/or be subjects of co-translational degradation (“death”). This review provides an overview of our current understanding of the complex co-translational events that accompany the synthesis, maturation, folding and degradation of nascent polypeptide chains.

Protein folding at the exit tunnel

Over five decades of research have yielded a large body of information on how purified proteins attain their native state when refolded in the test tube, starting from a chemically or thermally denatured state. Nevertheless, we still know little about how proteins fold and unfold in their natural biological habitat: the living cell. Indeed, a variety of cellular components, including molecular chaperones, the ribosome, and crowding of the intracellular medium, modulate folding mechanisms in physiologically relevant environments. This review focuses on the current state of knowledge in protein folding in the cell with emphasis on the early stage of a protein's life, as the nascent polypeptide traverses and emerges from the ribosomal tunnel. Given the vectorial nature of ribosome-assisted translation, the transient degree of chain elongation becomes a relevant variable expected to affect nascent protein foldability, aggregation propensity and extent of interaction with chaperones and the ribosome.

Folding at the birth of the nascent chain: coordinating translation with co-translational folding

In the living cells, the folding of many proteins is largely believed to begin co-translationally, during their biosynthesis at the ribosomes. In the ribosomal tunnel, the nascent peptide may establish local interactions and stabilize α-helical structures. Long-range contacts are more likely outside the ribosomes after release of larger segments of the nascent chain. Examples suggest that domains can attain native-like structure on the ribosome with and without population of folding intermediates. The co-translational folding is limited by the speed of the gradual extrusion of the nascent peptide which imposes conformational restraints on its folding landscape. Recent experimental and in silico modeling studies indicate that translation kinetics fine-tunes co-translational folding by providing a time delay for sequential folding of distinct portions of the nascent chain.

Identification of three single nucleotide polymorphisms in Anopheles gambiae immune signaling genes that are associated with natural Plasmodium falciparum infection

Background

Laboratory studies have demonstrated that a variety of immune signaling pathways regulate malaria parasite infection in Anopheles gambiae, the primary vector species in Africa.

Methods

To begin to understand the importance of these associations under natural conditions, an association mapping approach was adopted to determine whether single nucleotide polymorphisms (SNPs) in selected immune signaling genes in A. gambiae collected in Mali were associated with the phenotype of Plasmodium falciparum infection.

Results

Three SNPs were identified in field-collected mosquitoes that were associated with parasite infection in molecular form-dependent patterns: two were detected in the Toll5B gene and one was detected in the gene encoding insulin-like peptide 3 precursor. In addition, one infection-associated Toll5B SNP was in linkage disequilibrium with a SNP in sequence encoding a mitogen-activated protein kinase that has been associated with Toll signaling in mammalian cells. Both Toll5B SNPs showed divergence from Hardy-Weinberg equilibrium, suggesting that selection pressure(s) are acting on these loci.

Conclusions

Seven of these eight infection-associated and linked SNPs alter codon frequency or introduce non-synonymous changes that would be predicted to alter protein structure and, hence, function, suggesting that these SNPs could alter immune signaling and responsiveness to parasite infection.

Synonymous codon usage influences the local protein structure observed

Translation of mRNA into protein is a unidirectional information flow process. Analysing the input (mRNA) and output (protein) of translation, we find that local protein structure information is encoded in the mRNA nucleotide sequence. The Coding Sequence and Structure (CSandS) database developed in this work provides a detailed mapping between over 4000 solved protein structures and their mRNA. CSandS facilitates a comprehensive analysis of codon usage over many organisms. In assigning translation speed, we find that relative codon usage is less informative than tRNA concentration. For all speed measures, no evidence was found that domain boundaries are enriched with slow codons. In fact, genes seemingly avoid slow codons around structurally defined domain boundaries. Translation speed, however, does decrease at the transition into secondary structure. Codons are identified that have structural preferences significantly different from the amino acid they encode. However, each organism has its own set of ‘significant codons’. Our results support the premise that codons encode more information than merely amino acids and give insight into the role of translation in protein folding.

Directionality in protein fold prediction

BACKGROUND

Ever since the ground-breaking work of Anfinsen et al. in which a denatured protein was found to refold to its native state, it has been frequently stated by the protein fold prediction community that all the information required for protein folding lies in the amino acid sequence. Recent in vitro experiments and in silico computational studies, however, have shown that cotranslation may affect the folding pathway of some proteins, especially those of ancient folds. In this paper aspects of cotranslational folding have been incorporated into a protein structure prediction algorithm by adapting the Rosetta program to fold proteins as the nascent chain elongates. This makes it possible to conduct a pairwise comparison of folding accuracy, by comparing folds created sequentially from each end of the protein.

RESULTS

A single main result emerged: in 94% of proteins analyzed, following the sense of translation, from N-terminus to C-terminus, produced better predictions than following the reverse sense of translation, from the C-terminus to N-terminus. Two secondary results emerged. First, this superiority of N-terminus to C-terminus folding was more marked for proteins showing stronger evidence of cotranslation and second, an algorithm following the sense of translation produced predictions comparable to, and occasionally better than, Rosetta.

CONCLUSIONS

There is a directionality effect in protein fold prediction. At present, prediction methods appear to be too noisy to take advantage of this effect; as techniques refine, it may be possible to draw benefit from a sequential approach to protein fold prediction.

Detecting positive and purifying selection at synonymous sites in yeast and worm

We present a new computational method to identify positive and purifying selection at synonymous sites in yeast and worm. We define synonymous substitutions that change codons from preferred to unpreferred or vice versa as nonconservative synonymous substitutions and all other substitutions as conservative. Using a maximum-likelihood framework, we then test whether conservative and nonconservative synonymous substitutions occur at equal rates. Our approach replaces the standard rate of synonymous substitutions per synonymous site, dS, with two new rates, the conservative synonymous substitution rate (dS(C)) and the nonconservative synonymous substitution rate (dS(N)). Based on the ratio dS(N)/dS(C), we find that 0.05% of all yeast genes and none of worm genes show evidence of positive selection at synonymous sites (dS(N)/dS(C) > 1). On the other hand, 9.44% of all yeast genes and 5.12% of all worm genes show evidence of significant purifying selection on synonymous sites (dS(N)/dS(C) < 1). We also find that dS(N) correlates strongly with gene expression level, whereas the correlation between expression level and dS(C) is very weak. Thus, dS(N) captures most of the signal of selection for translational accuracy and speed, whereas dS(C) is not strongly influenced by this selection pressure. We suggest that the ratio dN/dS(C) may be more appropriate than the ratio dN/dS to identify positive or purifying selection on amino acids.

Increased incidence of rare codon clusters at 5' and 3' gene termini: implications for function

BACKGROUND

The process of translation can be affected by the use of rare versus common codons within the mRNA transcript.

RESULTS

Here, we show that rare codons are enriched at the 5' and 3' termini of genes from E. coli and other prokaryotes. Genes predicted to be secreted show significant enrichment in 5' rare codon clusters, but not 3' rare codon clusters. Surprisingly, no correlation between 5' mRNA structure and rare codon usage was observed.

CONCLUSIONS

Potential functional roles for the enrichment of rare codons at terminal positions are explored.

Structure and function of the molecular chaperone Trigger Factor

Newly synthesized proteins often require the assistance of molecular chaperones to efficiently fold into functional three-dimensional structures. At first, ribosome-associated chaperones guide the initial folding steps and protect growing polypeptide chains from misfolding and aggregation. After that folding into the native structure may occur spontaneously or require support by additional chaperones which do not bind to the ribosome such as DnaK and GroEL. Here we review the current knowledge on the best-characterized ribosome-associated chaperone at present, the Escherichia coli Trigger Factor. We describe recent progress on structural and dynamic aspects of Trigger Factor's interactions with the ribosome and substrates and discuss how these interactions affect co-translational protein folding. In addition, we discuss the newly proposed ribosome-independent function of Trigger Factor as assembly factor of multi-subunit protein complexes. Finally, we cover the functional cooperation between Trigger Factor, DnaK and GroEL in folding of cytosolic proteins and the interplay between Trigger Factor and other ribosome-associated factors acting in enzymatic processing and translocation of nascent polypeptide chains.

Low Complexity Regions behave as tRNA sponges to help co-translational folding of plasmodial proteins

In most organisms, the information necessary to specify the native 3D-structures of proteins is encoded in the corresponding mRNA sequences. Translational accuracy and efficiency are coupled and sequences that are slowly translated play an essential role in the concomitant folding of protein domains. Here, we suggest that the well-known mechanisms for the regulation of translational efficiency, which involves mRNA structure and/or asymmetric tRNA abundance, do not apply to all organisms. We propose that Plasmodium, the parasite responsible for malaria, uses an alternative strategy to slow down ribosomal speed and avoid multidomain protein misfolding during translation. In our model, the abundant Low Complexity Regions present in Plasmodium proteins replace the codon preferences, which influence the assembly of protein secondary structures.

Measuring cotranslational folding of nascent polypeptide chains on ribosomes

Protein folding has been studied extensively in vitro, but much less is known about how folding proceeds in vivo. A particular distinction of folding in vivo is that folding begins while the nascent polypeptide chain is still undergoing synthesis by the ribosome. Studies of cotranslational protein folding are inherently much more complex than classical in vitro protein folding studies, and historically there have been few methods available to produce the quantities of pure material required for biophysical studies of the nascent chain, or assays to specifically interrogate its conformation. However, the past few years have produced dramatic methodological advances, which now place cotranslational folding studies within reach of more biochemists, enabling a detailed comparison of the earliest stages of protein folding on the ribosome to the wealth of information available for the refolding of full-length polypeptide chains in vitro.

The comprehensive microbial resource

The Comprehensive Microbial Resource or CMR (http://cmr.jcvi.org) provides a web-based central resource for the display, search and analysis of the sequence and annotation for complete and publicly available bacterial and archaeal genomes. In addition to displaying the original annotation from GenBank, the CMR makes available secondary automated structural and functional annotation across all genomes to provide consistent data types necessary for effective mining of genomic data. Precomputed homology searches are stored to allow meaningful genome comparisons. The CMR supplies users with over 50 different tools to utilize the sequence and annotation data across one or more of the 571 currently available genomes. At the gene level users can view the gene annotation and underlying evidence. Genome level information includes whole genome graphical displays, biochemical pathway maps and genome summary data. Comparative tools display analysis between genomes with homology and genome alignment tools, and searches across the accessions, annotation, and evidence assigned to all genes/genomes are available. The data and tools on the CMR aid genomic research and analysis, and the CMR is included in over 200 scientific publications. The code underlying the CMR website and the CMR database are freely available for download with no license restrictions.

On finding poorly translated codons based on their usage frequency

Long stretches of “rare” codons are known to severely inhibit the efficiency of translation. Understanding the distribution of such rare codons is of critical importance in improving the efficiency of heterologous gene expression systems. Accurate estimates of codon usage take the abundance of each protein into consideration. In this paper, we analyze the correlation between approximate measures of codon usage and the availability of tRNA at various growth rates in E coli. We show that the computationally derived estimates of tRNA isoacceptor concentration enable the finding of poorly translated codons.

The ribosome as a platform for co-translational processing, folding and targeting of newly synthesized proteins

The early events in the life of newly synthesized proteins in the cellular environment are remarkably complex. Concurrently with their synthesis by the ribosome, nascent polypeptides are subjected to enzymatic processing, chaperone-assisted folding or targeting to translocation pores at membranes. The ribosome itself has a key role in these different tasks and governs the interplay between the various factors involved. Indeed, the ribosome serves as a platform for the spatially and temporally regulated association of enzymes, targeting factors and chaperones that act upon the nascent polypeptides emerging from the exit tunnel. Furthermore, the ribosome provides opportunities to coordinate the protein-synthesis activity of its peptidyl transferase center with the protein targeting and folding processes. Here we review the early co-translational events involving the ribosome that guide cytosolic proteins to their native state.

Clustering of codons with rare cognate tRNAs in human genes suggests an extra level of expression regulation

In species with large effective population sizes, highly expressed genes tend to be encoded by codons with highly abundant cognate tRNAs to maximize translation rate. However, there has been little evidence for a similar bias of synonymous codons in highly expressed human genes. Here, we ask instead whether there is evidence for the selection for codons associated with low abundance tRNAs. Rather than averaging the codon usage of complete genes, we scan the genes for windows with deviating codon usage. We show that there is a significant over representation of human genes that contain clusters of codons with low abundance cognate tRNAs. We name these regions, which on average have a 50% reduction in the amount of cognate tRNA available compared to the remainder of the gene, RTS (rare tRNA score) clusters. We observed a significant reduction in the substitution rate between the human RTS clusters and their orthologous chimp sequence, when compared to non-RTS cluster sequences. Overall, the genes with an RTS cluster have higher tissue specificity than the non-RTS cluster genes. Furthermore, these genes are functionally enriched for transcription regulation. As genes that regulate transcription in lower eukaryotes are known to be involved in translation on demand, this suggests that the mechanism of translation level expression regulation also exists within the human genome.

Generic algorithm to predict the speed of translational elongation: implications for protein biogenesis

Synonymous codon usage and variations in the level of isoaccepting tRNAs exert a powerful selective force on translation fidelity. We have developed an algorithm to evaluate the relative rate of translation which allows large-scale comparisons of the non-uniform translation rate on the protein biogenesis. Using the complete genomes of Escherichia coli and Bacillus subtilis we show that stretches of codons pairing to minor tRNAs form putative sites to locally attenuate translation; thereby the tendency is to cluster in near proximity whereas long contiguous stretches of slow-translating triplets are avoided. The presence of slow-translating segments positively correlates with the protein length irrespective of the protein abundance. The slow-translating clusters are predominantly located down-stream of the domain boundaries presumably to fine-tune translational accuracy with the folding fidelity of multidomain proteins. Translation attenuation patterns at highly structurally and functionally conserved domains are preserved across the species suggesting a concerted selective pressure on the codon selection and species-specific tRNA abundance in these regions.