Rare codons cluster

Wheat NAC transcription factor NAC5-1 is a positive regulator of senescence

Wheat (Triticum aestivum L.) is an important source of both calories and protein in global diets, but there is a trade-off between grain yield and protein content. The timing of leaf senescence could mediate this trade-off as it is associated with both declines in photosynthesis and nitrogen remobilization from leaves to grain. NAC transcription factors play key roles in regulating senescence timing. In rice, OsNAC5 expression is correlated with increased protein content and upregulated in senescing leaves, but the role of the wheat ortholog in senescence had not been characterized. We verified that NAC5-1 is the ortholog of OsNAC5 and that it is expressed in senescing flag leaves in wheat. To characterize NAC5-1, we combined missense mutations in NAC5-A1 and NAC5-B1 from a TILLING mutant population and overexpressed NAC5-A1 in wheat. Mutation in NAC5-1 was associated with delayed onset of flag leaf senescence, while overexpression of NAC5-A1 was associated with slightly earlier onset of leaf senescence. DAP-seq was performed to locate transcription factor binding sites of NAC5-1. Analysis of DAP-seq and comparison with other studies identified putative downstream target genes of NAC5-1 which could be associated with senescence. This work showed that NAC5-1 is a positive transcriptional regulator of leaf senescence in wheat. Further research is needed to test the effect of NAC5-1 on yield and protein content in field trials, to assess the potential to exploit this senescence regulator to develop high-yielding wheat while maintaining grain protein content.

Profiling expression strategies for a type III polyketide synthase in a lysate-based, cell-free system

Some of the most metabolically diverse species of bacteria (e.g., Actinobacteria) have higher GC content in their DNA, differ substantially in codon usage, and have distinct protein folding environments compared to tractable expression hosts like Escherichia coli. Consequentially, expressing biosynthetic gene clusters (BGCs) from these bacteria in E. coli often results in a myriad of unpredictable issues with regard to protein expression and folding, delaying the biochemical characterization of new natural products. Current strategies to achieve soluble, active expression of these enzymes in tractable hosts can be a lengthy trial-and-error process. Cell-free expression (CFE) has emerged as a valuable expression platform as a testbed for rapid prototyping expression parameters. Here, we use a type III polyketide synthase from Streptomyces griseus, RppA, which catalyzes the formation of the red pigment flaviolin, as a reporter to investigate BGC refactoring techniques. We applied a library of constructs with different combinations of promoters and rppA coding sequences to investigate the synergies between promoter and codon usage. Subsequently, we assess the utility of cell-free systems for prototyping these refactoring tactics prior to their implementation in cells. Overall, codon harmonization improves natural product synthesis more than traditional codon optimization across cell-free and cellular environments. More importantly, the choice of coding sequences and promoters impact protein expression synergistically, which should be considered for future efforts to use CFE for high-yield protein expression. The promoter strategy when applied to RppA was not completely correlated with that observed with GFP, indicating that different promoter strategies should be applied for different proteins. In vivo experiments suggest that there is correlation, but not complete alignment between expressing in cell free and in vivo. Refactoring promoters and/or coding sequences via CFE can be a valuable strategy to rapidly screen for catalytically functional production of enzymes from BCGs, which advances CFE as a tool for natural product research.

Metabolism-dependent secondary effect of anti-MAPK cancer therapy on DNA repair

Amino acid bioavailability impacts mRNA translation in a codon-dependent manner. Here, we report that the anti-cancer MAPK inhibitors (MAPKi) decrease the intracellular concentration of aspartate and glutamate in melanoma cells. This coincides with the accumulation of ribosomes on codons corresponding to these amino acids and triggers the translation-dependent degradation of mRNAs encoding aspartate- and glutamate-rich proteins, involved in DNA metabolism such as DNA replication and repair. Consequently, cells that survive MAPKi degrade aspartate and glutamate likely to generate energy, which simultaneously decreases their requirement for amino acids due to the downregulation of aspartate- and glutamate-rich proteins involved in cell proliferation. Concomitantly, the downregulation of aspartate- and glutamate-rich proteins involved in DNA repair increases DNA damage loads. Thus, DNA repair defects, and therefore mutations, are at least in part a secondary effect of the metabolic adaptation of cells exposed to MAPKi.

A small stretch of poor codon usage at the beginning of dengue virus open reading frame may act as a translational checkpoint

OBJECTIVE

Rare codons were previously shown to be enriched at the beginning of the dengue virus (DENV) open reading frame. However, the role of rare codons in regulating translation efficiency and replication of DENV remains unclear. The present study aims to clarify the significance of rare codon usage at the beginning of DENV transcripts using the codon adaptation index (CAI).

METHODOLOGY

CAIs of the whole starting regions of DENV transcripts as well as 18-codon sliding windows of the regions were analyzed.

RESULTS

One of the intriguing findings is that those rare codons do not typically result in uniformly low CAI in the starting region with rare codons. However, it shows a notable local drop in CAI around the 50th codon in all dengue serotypes. This suggests that there may be a translational checkpoint at this site and that the rare codon usage upstream to this checkpoint may not be related to translational control.

Profiling Expression Strategies for a Type III Polyketide Synthase in a Lysate-Based, Cell-free System

Some of the most metabolically diverse species of bacteria (e.g., Actinobacteria) have higher GC content in their DNA, differ substantially in codon usage, and have distinct protein folding environments compared to tractable expression hosts like Escherichia coli. Consequentially, expressing biosynthetic gene clusters (BGCs) from these bacteria in E. coli frequently results in a myriad of unpredictable issues with protein expression and folding, delaying the biochemical characterization of new natural products. Current strategies to achieve soluble, active expression of these enzymes in tractable hosts, such as BGC refactoring, can be a lengthy trial-and-error process. Cell-free expression (CFE) has emerged as 1) a valuable expression platform for enzymes that are challenging to synthesize in vivo, and as 2) a testbed for rapid prototyping that can improve cellular expression. Here, we use a type III polyketide synthase from Streptomyces griseus, RppA, which catalyzes the formation of the red pigment flaviolin, as a reporter to investigate BGC refactoring techniques. We synergistically tune promoter and codon usage to improve flaviolin production from cell-free expressed RppA. We then assess the utility of cell-free systems for prototyping these refactoring tactics prior to their implementation in cells. Overall, codon harmonization improves natural product synthesis more than traditional codon optimization across cell-free and cellular environments. Refactoring promoters and/or coding sequences via CFE can be a valuable strategy to rapidly screen for catalytically functional production of enzymes from BCGs. By showing the coordinators between CFE versus in vivo expression, this work advances CFE as a tool for natural product research.

An exploratory in silico comparison of open-source codon harmonization tools

BACKGROUND

Not changing the native constitution of genes prior to their expression by a heterologous host can affect the amount of proteins synthesized as well as their folding, hampering their activity and even cell viability. Over the past decades, several strategies have been developed to optimize the translation of heterologous genes by accommodating the difference in codon usage between species. While there have been a handful of studies assessing various codon optimization strategies, to the best of our knowledge, no research has been performed towards the evaluation and comparison of codon harmonization algorithms. To highlight their importance and encourage meaningful discussion, we compared different open-source codon harmonization tools pertaining to their in silico performance, and we investigated the influence of different gene-specific factors.

RESULTS

In total, 27 genes were harmonized with four tools toward two different heterologous hosts. The difference in %MinMax values between the harmonized and the original sequences was calculated (ΔMinMax), and statistical analysis of the obtained results was carried out. It became clear that not all tools perform similarly, and the choice of tool should depend on the intended application. Almost all biological factors under investigation (GC content, RNA secondary structures and choice of heterologous host) had a significant influence on the harmonization results and thus must be taken into account. These findings were substantiated using a validation dataset consisting of 8 strategically chosen genes.

CONCLUSIONS

Due to the size of the dataset, no complex models could be developed. However, this initial study showcases significant differences between the results of various codon harmonization tools. Although more elaborate investigation is needed, it is clear that biological factors such as GC content, RNA secondary structures and heterologous hosts must be taken into account when selecting the codon harmonization tool.

Effect of mRNA/tRNA mutations on translation speed: Implications for human diseases

Recent discoveries establish tRNAs as central regulators of mRNA translation dynamics, and therefore cotranslational folding and function of the encoded protein. The tRNA pool, whose composition and abundance change in a cell- and tissue-dependent manner, is the main factor which determines mRNA translation velocity. In this review, we discuss a group of pathogenic mutations, in the coding sequences of either protein-coding genes or in tRNA genes, that alter mRNA translation dynamics. We also summarize advances in tRNA biology that have uncovered how variations in tRNA levels on account of genetic mutations affect protein folding and function, and thereby contribute to phenotypic diversity in clinical manifestations.

Base Composition, Codon Usage, and Patterns of Gene Sequence Evolution in Butterflies

Coding sequence evolution is influenced by both natural selection and neutral evolutionary forces. In many species, the effects of mutation bias, codon usage, and GC-biased gene conversion (gBGC) on gene sequence evolution have not been detailed. Quantification of how these forces shape substitution patterns is therefore necessary to understand the strength and direction of natural selection. Here, we used comparative genomics to investigate the association between base composition and codon usage bias on gene sequence evolution in butterflies and moths (Lepidoptera), including an in-depth analysis of underlying patterns and processes in one species, Leptidea sinapis. The data revealed significant G/C to A/T substitution bias at third codon position with some variation in the strength among different butterfly lineages. However, the substitution bias was lower than expected from previously estimated mutation rate ratios, partly due to the influence of gBGC. We found that A/T-ending codons were overrepresented in most species, but there was a positive association between the magnitude of codon usage bias and GC-content in third codon positions. In addition, the tRNA-gene population in L. sinapis showed higher GC-content at third codon positions compared to coding sequences in general and less overrepresentation of A/T-ending codons. There was an inverse relationship between synonymous substitutions and codon usage bias indicating selection on synonymous sites. We conclude that the evolutionary rate in Lepidoptera is affected by a complex interaction between underlying G/C -> A/T mutation bias and partly counteracting fixation biases, predominantly conferred by overall purifying selection, gBGC, and selection on codon usage.

Multiplex genome engineering in Clostridium beijerinckii NCIMB 8052 using CRISPR-Cas12a

Clostridium species are re-emerging as biotechnological workhorses for industrial acetone-butanol-ethanol production. This re-emergence is largely due to advances in fermentation technologies but also due to advances in genome engineering and re-programming of the native metabolism. Several genome engineering techniques have been developed including the development of numerous CRISPR-Cas tools. Here, we expanded the CRISPR-Cas toolbox and developed a CRISPR-Cas12a genome engineering tool in Clostridium beijerinckii NCIMB 8052. By controlling the expression of FnCas12a with the xylose-inducible promoter, we achieved efficient (25-100%) single-gene knockout of five C. beijerinckii NCIMB 8052 genes (spo0A, upp, Cbei_1291, Cbei_3238, Cbei_3832). Moreover, we achieved multiplex genome engineering by simultaneously knocking out the spo0A and upp genes in a single step with an efficiency of 18%. Finally, we showed that the spacer sequence and position in the CRISPR array can affect the editing efficiency outcome.

Application of codon usage and context analysis in genes up- or down-regulated in neurodegeneration and cancer to combat comorbidities

Introduction

Neurodegeneration and cancer present in comorbidities with inverse effects due to the expression of genes and pathways acting in opposition. Identifying and studying the genes simultaneously up or downregulated during morbidities helps curb both ailments together.

Methods

This study examines four genes. Three of these (Amyloid Beta Precursor Protein (APP), Cyclin D1 (CCND1), and Cyclin E2 (CCNE2) are upregulated, and one protein phosphatase 2 phosphatase activator (PTPA) is simultaneously downregulated in both disorders. We investigated molecular patterns, codon usage, codon usage bias, nucleotide bias in the third codon position, preferred codons, preferred codon pairs, rare codons, and codon context.

Results

Parity analysis revealed that T is preferred over A, and G is preferred over C in the third codon position, suggesting composition plays no role in nucleotide bias in both the upregulated and downregulated gene sets and that mutational forces are stronger in upregulated gene sets than in downregulated ones. Transcript length influenced the overall %A composition and codon bias, and the codon AGG exerted the strongest influence on codon usage in both the upregulated and downregulated gene sets. Codons ending in G/C were preferred for 16 amino acids, and glutamic acid-, aspartic acid-, leucine-, valine-, and phenylalanine-initiated codon pairs were preferred in all genes. Codons CTA (Leu), GTA (Val), CAA (Gln), and CGT (Arg) were underrepresented in all examined genes.

Discussion

Using advanced gene editing tools such as CRISPR/Cas or any other gene augmentation technique, these recoded genes may be introduced into the human body to optimize gene expression levels to augment neurodegeneration and cancer therapeutic regimens simultaneously.

Differential production rates of cytosolic and transmembrane GFP reporters in C. elegans L3 larval uterine cells

Transgene driven protein expression is an important tool for investigating developmental mechanisms in C. elegans . Here, we have assessed protein production rates and levels in L3 larval uterine cells (UCs). Using ubiquitous promoter driven cytosolic and transmembrane tethered GFP, fluorescence recovery after photobleaching, and quantitative fluorescence analysis, we reveal that cytosolic GFP is produced at an ~two-fold higher rate than transmembrane tethered GFP and accumulates at ~five-fold higher levels in UCs. We also provide evidence that cytosolic GFP in the anchor cell, a specialized UC that mediates uterine-vulval connection, is more rapidly degraded through an autophagy-independent mechanism.

Effects of synonymous mutations on kinetic properties and structure of firefly luciferase: Molecular dynamics simulation, molecular docking, RNA folding, and experimental study

Although synonymous mutations have long been thought to lack striking results, a growing body of research shows these mutations have highly variable effects. In this study, the impact of synonymous mutations in the development of thermostable luciferase was investigated using a combination of experimental and theoretical approaches. Using bioinformatics analysis, the codon usage features in the Lampyridae family's luciferases were studied and four synonymous mutations of Arg in luciferase were created. An exciting result was that the analysis of kinetic parameters showed a slight increase in the thermal stability of the mutant luciferase. AutoDock Vina, %MinMax algorithm, and UNAFold Server were used to perform molecular docking, folding rate, and RNA folding, respectively. Here, it was assumed that in the region (Arg337) with a moderate propensity for coil, synonymous mutation altered the rate of translation, which in turn may lead to a slight change in the structure of the enzyme. According to the molecular dynamics simulation data, local minor global flexibility is observed in the context of the protein conformation. A plausible explanation is that this flexibility may strengthen hydrophobic interactions due to its sensitivity to a molecular collision. Accordingly, thermostability originated mainly from hydrophobic interaction.

Folding of heterologous proteins in bacterial cell factories: Cellular mechanisms and engineering strategies

The expression of correctly folded and functional heterologous proteins is important in many biotechnological production processes, whether it is enzymes, biopharmaceuticals or biosynthetic pathways for production of sustainable chemicals. For industrial applications, bacterial platform organisms, such as E. coli, are still broadly used due to the availability of tools and proven suitability at industrial scale. However, expression of heterologous proteins in these organisms can result in protein aggregation and low amounts of functional protein. This review provides an overview of the cellular mechanisms that can influence protein folding and expression, such as co-translational folding and assembly, chaperone binding, as well as protein quality control, across different model organisms. The knowledge of these mechanisms is then linked to different experimental methods that have been applied in order to improve functional heterologous protein folding, such as codon optimization, fusion tagging, chaperone co-production, as well as strain and protein engineering strategies.

Synonymous and Nonsynonymous Substitutions in Dictyostelium discoideum Ammonium Transporter amtA Are Necessary for Functional Complementation in Saccharomyces cerevisiae

Ammonium transporters are present in all three domains of life. They have undergone extensive horizontal gene transfer (HGT), gene duplication, and functional diversification and therefore offer an excellent paradigm to study protein evolution. We attempted to complement a mep1Δmep2Δmep3Δ strain of Saccharomyces cerevisiae (triple-deletion strain), which otherwise cannot grow on ammonium as a sole nitrogen source at concentrations of <3 mM, with amtA of Dictyostelium discoideum, an orthologue of S. cerevisiae MEP2. We observed that amtA did not complement the triple-deletion strain of S. cerevisiae for growth on low-ammonium medium. We isolated two mutant derivatives of amtA (amtA M1 and amtA M2) from a PCR-generated mutant plasmid library that complemented the triple-deletion strain of S. cerevisiae. amtA M1 bears three nonsynonymous and two synonymous substitutions, which are necessary for its functionality. amtA M2 bears two nonsynonymous substitutions and one synonymous substitution, all of which are necessary for functionality. Interestingly, AmtA M1 transports ammonium but does not confer methylamine toxicity, while AmtA M2 transports ammonium and confers methylamine toxicity, demonstrating functional diversification. Preliminary biochemical analyses indicated that the mutants differ in their conformations as well as their mechanisms of ammonium transport. These intriguing results clearly point out that protein evolution cannot be fathomed by studying nonsynonymous and synonymous substitutions in isolation. The above-described observations have significant implications for various facets of biological processes and are discussed in detail. IMPORTANCE Functional diversification following gene duplication is one of the major driving forces of protein evolution. While the role of nonsynonymous substitutions in the functional diversification of proteins is well recognized, knowledge of the role of synonymous substitutions in protein evolution is in its infancy. Using functional complementation, we isolated two functional alleles of the D. discoideum ammonium transporter gene (amtA), which otherwise does not function in S. cerevisiae as an ammonium transporters. One of them is an ammonium transporter, while the other is an ammonium transporter that also confers methylammonium (ammonium analogue) toxicity, suggesting functional diversification. Surprisingly, both alleles require a combination of synonymous and nonsynonymous substitutions for their functionality. These results bring out a hitherto-unknown pathway of protein evolution and pave the way for not only understanding protein evolution but also interpreting single nucleotide polymorphisms (SNPs).

Predictive Power of Long-Read Whole-Genome Sequencing for Rapid Diagnostics of Multidrug-Resistant Brachyspira hyodysenteriae Strains

Infections with Brachyspira hyodysenteriae, the etiological agent of swine dysentery, result in major economic losses in the pig industry worldwide. Even though microbial differentiation of various Brachyspira species can be obtained via PCR, no quick diagnostics for antimicrobial susceptibility testing are in place, which is mainly due to the time-consuming (4 to 7 days) anaerobic growth requirements of these organisms. Veterinarians often rely on a clinical diagnosis for initiating antimicrobial treatment. These treatments are not always effective, which may be due to high levels of acquired resistance in B. hyodysenteriae field isolates. By using long-read-only whole-genome sequencing and a custom-trained Bonito base-calling model, 81 complete B. hyodysenteriae genomes with median Q51 scores and 99% completeness were obtained from 86 field strains. This allowed the assessment of the predictive potential of genetic markers in relation to the observed acquired resistance phenotypes obtained via agar dilution susceptibility testing. Multidrug resistance was observed in 77% and 21% of the tested strains based on epidemiological cutoff and clinical breakpoint values, respectively. The predictive power of genetic hallmarks (genes and/or gene mutations) for antimicrobial susceptibility testing was promising. Sensitivity and specificity for tiamulin [tva(A) and 50SL3^N148S, 99% and 67%], valnemulin [tva(A), 97% and 92%), lincomycin (23S^A2153T/G and lnuC, 94% and 100%), tylvalosin (23S^A2153T/G, 99% and 93%), and doxycycline (16S^G1026C, 93% and 87%) were determined. The predictive power of these genetic hallmarks is promising for use in sequencing-based workflows to speed up swine dysentery diagnostics in veterinary medicine and determine proper antimicrobial use. IMPORTANCE Diagnostics for swine dysentery rely on the identification of Brachyspira species using molecular techniques. Nevertheless, no quick diagnostic tools are available for antimicrobial susceptibility testing due to extended growth requirements (7 to 14 days). To enable practitioners to tailor antimicrobial treatment to specific strains, long-read sequencing-based methods are expected to lead to rapid methods in the future. Nevertheless, their potential implementation should be validated extensively. This mainly implies assessing sequencing accuracy and the predictive power of genetic hallmarks in relation to their observed (multi)resistance phenotypes.

Codon Usage and Context Analysis of Genes Modulated during SARS-CoV-2 Infection and Dental Inflammation

The overexpression of SARS-CoV-2 primary receptors and co-receptors (ACE2, TMPRSS2, FURIN, and CD147) enhance the likeliness of SARS-CoV-2 infection. The genes for same receptors are overexpressed in the periodontal tissues of periodontitis patients. On the other hand, BMAL1 is recognized to play a crucial role in regulating pulmonary inflammation and enhancing susceptibility to viral infection. Silenced BMAL1 disrupts circadian transcriptional regulations, enhances vulnerability to SARS-CoV-2 infections, and may trigger the further production of TNF-α and other pro-inflammatory cytokines that propagate the cytokine storm and exacerbate periodontal inflammation. Therefore ACE2, TMPRSS2, FURIN, CD147, and BMAL1 are the crossroads between SARS-CoV-2 and Periodontitis genes. The enhanced expression of ACE2, TMPRSS2, FURIN, and CD147 and the diminished expression of BMAL1 may be a strategy to check both ailments simultaneously. In gene manipulation techniques, oligos are introduced, which contain all the necessary information to manipulate gene expression. The data are derived from the studies on genes' molecular patterns, including nucleotide composition, dinucleotide patterns, relative synonymous codon usage, codon usage bias, codon context, and rare and abundant codons. Such information may be used to manipulate the overexpression and underexpression of the genes at the time of SARS-CoV-2 infection and periodontitis to mitigate both ailments simultaneously; it can be explored to uncover possible future treatments.

Mechanisms governing codon usage bias and the implications for protein expression in the chloroplast of Chlamydomonas reinhardtii

Chloroplasts possess a considerably reduced genome that is decoded via an almost minimal set of tRNAs. These features make an excellent platform for gaining insights into fundamental mechanisms that govern protein expression. Here, we present a comprehensive and revised perspective of the mechanisms that drive codon selection in the chloroplast of Chlamydomonas reinhardtii and the functional consequences for protein expression. In order to extract this information, we applied several codon usage descriptors to genes with different expression levels. We show that highly expressed genes strongly favor translationally optimal codons, while genes with lower functional importance are rather affected by directional mutational bias. We demonstrate that codon optimality can be deduced from codon-anticodon pairing affinity and, for a small number of amino acids (leucine, arginine, serine, and isoleucine), tRNA concentrations. Finally, we review, analyze, and expand on the impact of codon usage on protein yield, secondary structures of mRNA, translation initiation and termination, and amino acid composition of proteins, as well as cotranslational protein folding. The comprehensive analysis of codon choice provides crucial insights into heterologous gene expression in the chloroplast of C. reinhardtii, which may also be applicable to other chloroplast-containing organisms and bacteria.

A unifying network modeling approach for codon optimization

MOTIVATION

Synthesizing genes to be expressed in other organisms is an essential tool in biotechnology. While the many-to-one mapping from codons to amino acids makes the genetic code degenerate, codon usage in a particular organism is not random either. This bias in codon use may have a remarkable effect on the level of gene expression. A number of measures have been developed to quantify a given codon sequence's strength to express a gene in a host organism. Codon optimization aims to find a codon sequence that will optimize one or more of these measures. Efficient computational approaches are needed since the possible number of codon sequences grows exponentially as the number of amino acids increases.

RESULTS

We develop a unifying modeling approach for codon optimization. With our mathematical formulations based on graph/network representations of amino acid sequences, any combination of measures can be optimized in the same framework by finding a path satisfying additional limitations in an acyclic layered network. We tested our approach on bi-objectives commonly used in the literature, namely, Codon Pair Bias versus Codon Adaptation Index and Relative Codon Pair Bias versus Relative Codon Bias. However, our framework is general enough to handle any number of objectives concurrently with certain restrictions or preferences on the use of specific nucleotide sequences. We implemented our models using Python's Gurobi interface and showed the efficacy of our approach even for the largest proteins available. We also provided experimentation showing that highly expressed genes have objective values close to the optimized values in the bi-objective codon design problem.

AVAILABILITY AND IMPLEMENTATION

http://alpersen.bilkent.edu.tr/NetworkCodon.zip.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Quantifying shifts in natural selection on codon usage between protein regions: a population genetics approach

Background

Codon usage bias (CUB), the non-uniform usage of synonymous codons, occurs across all domains of life. Adaptive CUB is hypothesized to result from various selective pressures, including selection for efficient ribosome elongation, accurate translation, mRNA secondary structure, and/or protein folding. Given the critical link between protein folding and protein function, numerous studies have analyzed the relationship between codon usage and protein structure. The results from these studies have often been contradictory, likely reflecting the differing methods used for measuring codon usage and the failure to appropriately control for confounding factors, such as differences in amino acid usage between protein structures and changes in the frequency of different structures with gene expression.

Results

Here we take an explicit population genetics approach to quantify codon-specific shifts in natural selection related to protein structure in S. cerevisiae and E. coli. Unlike other metrics of codon usage, our approach explicitly separates the effects of natural selection, scaled by gene expression, and mutation bias while naturally accounting for a region’s amino acid usage. Bayesian model comparisons suggest selection on codon usage varies only slightly between helix, sheet, and coil secondary structures and, similarly, between structured and intrinsically-disordered regions. Similarly, in contrast to prevous findings, we find selection on codon usage only varies slightly at the termini of helices in E. coli. Using simulated data, we show this previous work indicating “non-optimal” codons are enriched at the beginning of helices in S. cerevisiae was due to failure to control for various confounding factors (e.g. amino acid biases, gene expression, etc.), and rather than selection to modulate cotranslational folding.

Conclusions

Our results reveal a weak relationship between codon usage and protein structure, indicating that differences in selection on codon usage between structures are slight. In addition to the magnitude of differences in selection between protein structures being slight, the observed shifts appear to be idiosyncratic and largely codon-specific rather than systematic reversals in the nature of selection. Overall, our work demonstrates the statistical power and benefits of studying selective shifts on codon usage or other genomic features from an explicitly evolutionary approach. Limitations of this approach and future potential research avenues are discussed.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-022-08635-0).

Biosynthetic Protein Folding and Molecular Chaperons

The problem of linear polypeptide chain folding into a unique tertiary structure is one of the fundamental scientific challenges. The process of folding cannot be fully understood without its biological context, especially for big multidomain and multisubunit proteins. The principal features of biosynthetic folding are co-translational folding of growing nascent polypeptide chains and involvement of molecular chaperones in the process. The review summarizes available data on the early events of nascent chain folding, as well as on later advanced steps, including formation of elements of native structure. The relationship between the non-uniformity of translation rate and folding of the growing polypeptide is discussed. The results of studies on the effect of biosynthetic folding features on the parameters of folding as a physical process, its kinetics and mechanisms, are presented. Current understanding and hypotheses on the relationship of biosynthetic folding with the fundamental physical parameters and current views on polypeptide folding in the context of energy landscapes are discussed.

Rifampicin Increases Expression of Plant Codon-Optimized Bacillus thuringiensis δ-Endotoxin Genes in Escherichia coli

Transgenic crops expressing Cry δ-endotoxins of Bacillus thuringiensis for insect resistance have been commercialized worldwide with increased crop productivity and spectacular socioeconomic gains. To attain the enhanced level of protein expression, the cry genes have to be extensively modified for RNA stability and translation efficiency in the plant systems. However, such modifications in nucleotide sequences make it difficult to express the cry genes in Escherichia coli because of the presence of E. coli rare codons. Induction of gene expression through the T7 promoter/lac operator system results in high levels of transcription but limits the availability of activated tRNA corresponding to rare codons that leads to translation stalling at ribosomes. In the present study, an Isopropyl ß-D-1-thiogalactopyranoside (IPTG)/rifampicin combination-based approach was adopted to induce transcription of cry genes through T7 promoter/lac operator while simultaneously inhibiting the transcription of host genes through rifampicin. The results show that the IPTG/rifampicin combination leads to high-level expression of four plant codon-optimized cry genes (cry2Aa, cry1F, cry1Ac, and cry1AcF). Northern blot analysis of the cry gene expressing E. coli samples showed that the RNA expression level in the IPTG-induced samples was higher as compared to that in the IPTG/rifampicin-induced samples. Diet overlay insect bioassay of IPTG/rifampicin-induced Cry toxins with Helicoverpa armigera larvae showed bioactivity (measured as LC₅₀) similar to the previous studies. The experiment has proved that recombinant synthetic gene (plant codon-optimized gene) with the combination of Rifampicin which inhibits DNA-dependent bacterial RNA polymerase and reduces the excessive baggage of translational machinery of the bacterial cell triggers the production of synthetic protein. Purification of protein using high pH buffer increases the solubility of the protein. Further, LC₅₀ analysis shows no reduction of protein activity leads to protein stability. Further, purified cry toxin protein can be used for crop protection against pests and a purified form of the synthetic protein can be used for antibody production and perform the immunoassay for the identification of the transgenic plant. The crystallographic structure of synthetic protein could be used for interaction study with another insect to see insecticidal activity.

A switch from α-helical to β-strand conformation during co-translational protein folding

Cellular proteins begin to fold as they emerge from the ribosome. The folding landscape of nascent chains is not only shaped by their amino acid sequence but also by the interactions with the ribosome. Here, we combine biophysical methods with cryo‐EM structure determination to show that folding of a β‐barrel protein begins with formation of a dynamic α‐helix inside the ribosome. As the growing peptide reaches the end of the tunnel, the N‐terminal part of the nascent chain refolds to a β‐hairpin structure that remains dynamic until its release from the ribosome. Contacts with the ribosome and structure of the peptidyl transferase center depend on nascent chain conformation. These results indicate that proteins may start out as α‐helices inside the tunnel and switch into their native folds only as they emerge from the ribosome. Moreover, the correlation of nascent chain conformations with reorientation of key residues of the ribosomal peptidyl‐transferase center suggest that protein folding could modulate ribosome activity.

CHARMING: Harmonizing synonymous codon usage to replicate a desired codon usage pattern

There is a growing appreciation that synonymous codon usage, although historically regarded as phenotypically silent, can instead alter a wide range of mechanisms related to functional protein production, a term we use here to describe the net effect of transcription (mRNA synthesis), mRNA half-life, translation (protein synthesis) and the probability of a protein folding correctly to its active, functional structure. In particular, recent discoveries have highlighted the important role that sub-optimal codons can play in modifying co-translational protein folding. These results have drawn increased attention to the patterns of synonymous codon usage within coding sequences, particularly in light of the discovery that these patterns can be conserved across evolution for homologous proteins. Because synonymous codon usage differs between organisms, for heterologous gene expression it can be desirable to make synonymous codon substitutions to match the codon usage pattern from the original organism in the heterologous expression host. Here we present CHARMING (for Codon HARMonizING), a robust and versatile algorithm to design mRNA sequences for heterologous gene expression and other related codon harmonization tasks. CHARMING can be run as a downloadable Python script or via a web portal at http://www.codons.org.

Comprehensive Genome Engineering Toolbox for Microalgae Nannochloropsis oceanica Based on CRISPR-Cas Systems

Microalgae can produce industrially relevant metabolites using atmospheric CO₂ and sunlight as carbon and energy sources, respectively. Developing molecular tools for high-throughput genome engineering could accelerate the generation of tailored strains with improved traits. To this end, we developed a genome editing strategy based on Cas12a ribonucleoproteins (RNPs) and homology-directed repair (HDR) to generate scarless and markerless mutants of the microalga Nannochloropsis oceanica. We also developed an episomal plasmid-based Cas12a system for efficiently introducing indels at the target site. Additionally, we exploited the ability of Cas12a to process an associated CRISPR array to perform multiplexed genome engineering. We efficiently targeted three sites in the host genome in a single transformation, thereby making a major step toward high-throughput genome engineering in microalgae. Furthermore, a CRISPR interference (CRISPRi) tool based on Cas9 and Cas12a was developed for effective downregulation of target genes. We observed up to 85% reduction in the transcript levels upon performing CRISPRi with dCas9 in N. oceanica. Overall, these developments substantially accelerate genome engineering efforts in N. oceanica and potentially provide a general toolbox for improving other microalgal strains.

A synonymous mutation in IGF-1 impacts the transcription and translation process of gene expression

Insulin-like growth factor 1 (IGF-1) is considered to be a crucial gene in the animal development of bone and body size. In this study, a unique synonymous mutation (c.258 A > G) of the IGF-1 gene was modified with an adenine base editor to observe the growth and developmental situation of mutant mice. Significant expression differences and molecular mechanisms among vectors with different alanine synonymous codons were explored. Although modification of a single synonymous codon rarely interferes with animal phenotypes, we observed that the expression and secretion of IGF-1 were different between 8-week-old homozygous (Ho) and wild-type (WT) mice. In addition, the IGF-1 with optimal codon combinations showed a higher expression content than other codon combination modes at both transcription and translation levels and performed proliferation promotion. The gene stability and translation initiation efficiency also changed significantly. Our findings illustrated that the synonymous mutation altered the IGF-1 gene expression in individual mice and suggested that the synonymous mutation affected the IGF-1 expression and biological function through the transcription and translation processes.

Efficient expression of a cnidarian peptide-gated ion channel in mammalian cells

Hydra Na+ channels (HyNaCs) are peptide-gated ion channels of the DEG/ENaC gene family that are directly activated by neuropeptides of the Hydra nervous system. They have previously been successfully characterized in Xenopus oocytes. To establish their expression in mammalian cells, we transiently expressed heteromeric HyNaC2/3/5 in human HEK 293 and monkey COS-7 cells. We found that the expression of HyNaC2/3/5 using native cDNAs was inefficient and that codon optimization strongly increased protein expression and current amplitude in patch-clamp experiments. We used the improved expression of codon-optimized channel subunits to perform Ca2+ imaging and to demonstrate their glycosylation pattern. In summary, we established efficient expression of a cnidarian ion channel in mammalian cell lines.

Codon usage bias

Codon usage bias is the preferential or non-random use of synonymous codons, a ubiquitous phenomenon observed in bacteria, plants and animals. Different species have consistent and characteristic codon biases. Codon bias varies not only with species, family or group within kingdom, but also between the genes within an organism. Codon usage bias has evolved through mutation, natural selection, and genetic drift in various organisms. Genome composition, GC content, expression level and length of genes, position and context of codons in the genes, recombination rates, mRNA folding, and tRNA abundance and interactions are some factors influencing codon bias. The factors shaping codon bias may also be involved in evolution of the universal genetic code. Codon-usage bias is critical factor determining gene expression and cellular function by influencing diverse processes such as RNA processing, protein translation and protein folding. Codon usage bias reflects the origin, mutation patterns and evolution of the species or genes. Investigations of codon bias patterns in genomes can reveal phylogenetic relationships between organisms, horizontal gene transfers, molecular evolution of genes and identify selective forces that drive their evolution. Most important application of codon bias analysis is in the design of transgenes, to increase gene expression levels through codon optimization, for development of transgenic crops. The review gives an overview of deviations of genetic code, factors influencing codon usage or bias, codon usage bias of nuclear and organellar genes, computational methods to determine codon usage and the significance as well as applications of codon usage analysis in biological research, with emphasis on plants.

In Silico Analysis of Dengue Virus Serotype 2 Mutations Detected at the Intrahost Level in Patients with Different Clinical Outcomes

Intrahost genetic diversity is thought to facilitate arbovirus adaptation to changing environments and hosts, and it may also be linked to viral pathogenesis. Intending to shed light on the viral determinants for severe dengue pathogenesis, we previously analyzed the DENV-2 intrahost genetic diversity in 68 patients clinically classified as dengue fever (n = 31), dengue with warning signs (n = 19), and severe dengue (n = 18), performing viral whole-genome deep sequencing from clinical samples with an amplicon-free approach. From it, we identified a set of 141 relevant mutations distributed throughout the viral genome that deserved further attention. Therefore, we employed molecular modeling to recreate three-dimensional models of the viral proteins and secondary RNA structures to map the mutations and assess their potential effects. Results showed that, in general lines, disruptive variants were identified primarily among dengue fever cases. In contrast, potential immune-escape variants were associated mainly with warning signs and severe cases, in line with the latter's longer intrahost evolution times. Furthermore, several mutations were located on protein-surface regions, with no associated function. They could represent sites of further investigation, as the interaction of viral and host proteins is critical for both host immunomodulation and virus hijacking of the cellular machinery. The present analysis provides new information about the implications of the intrahost genetic diversity of DENV-2, contributing to the knowledge about the viral factors possibly involved in its pathogenesis within the human host. Strengthening our results with functional studies could allow many of these variants to be considered in the design of therapeutic or prophylactic compounds and the improvement of diagnostic assays. IMPORTANCE Previous evidence showed that intrahost genetic diversity in arboviruses may be linked to viral pathogenesis and that one or a few amino acid replacements within a single protein are enough to modify a biological feature of an RNA virus. To assess dengue virus serotype 2 determinants potentially involved in pathogenesis, we previously analyzed the intrahost genetic diversity of the virus in patients with different clinical outcomes and identified a set of 141 mutations that deserved further study. Thus, through a molecular modeling approach, we showed that disruptive variants were identified primarily among cases with mild dengue fever, while potential immune-escape variants were mainly associated with cases of greater severity. We believe that some of the variants pointed out in this study were attractive enough to be potentially considered in future intelligent designs of therapeutic or prophylactic compounds or the improvement of diagnostic tools. The present analysis provides new information about DENV-2 viral factors possibly involved in its pathogenesis within the human host.

Characterizing RNA Pseudouridylation by Convolutional Neural Networks

Pseudouridine (Ψ) is the most prevalent post-transcriptional RNA modification and is widespread in small cellular RNAs and mRNAs. However, the functions, mechanisms, and precise distribution of Ψs (especially in mRNAs) still remain largely unclear. The landscape of Ψs across the transcriptome has not yet been fully delineated. Here, we present a highly effective model based on a convolutional neural network (CNN), called PseudoUridyLation Site Estimator (PULSE), to analyze large-scale profiling data of Ψ sites and characterize the contextual sequence features of pseudouridylation. PULSE, consisting of two alternatively-stacked convolution and pooling layers followed by a fully-connected neural network, can automatically learn the hidden patterns of pseudouridylation from the local sequence information. Extensive validation tests demonstrated that PULSE can outperform other state-of-the-art prediction methods and achieve high prediction accuracy, thus enabling us to further characterize the transcriptome-wide landscape of Ψ sites. We further showed that the prediction results derived from PULSE can provide novel insights into understanding the functional roles of pseudouridylation, such as the regulations of RNA secondary structure, codon usage, translation, and RNA stability, and the connection to single nucleotide variants. The source code and final model for PULSE are available at https://github.com/mlcb-thu/PULSE.

A Code Within a Code: How Codons Fine-Tune Protein Folding in the Cell

The genetic code sets the correspondence between the sequence of a given nucleotide triplet in an mRNA molecule, called a codon, and the amino acid that is added to the growing polypeptide chain during protein synthesis. With four bases (A, G, U, and C), there are 64 possible triplet codons: 61 sense codons (encoding amino acids) and 3 nonsense codons (so-called, stop codons that define termination of translation). In most organisms, there are 20 common/standard amino acids used in protein synthesis; thus, the genetic code is redundant with most amino acids (with the exception of Met and Trp) are being encoded by more than one (synonymous) codon. Synonymous codons were initially presumed to have entirely equivalent functions, however, the finding that synonymous codons are not present at equal frequencies in mRNA suggested that the specific codon choice might have functional implications beyond coding for amino acid. Observation of nonequivalent use of codons in mRNAs implied a possibility of the existence of auxiliary information in the genetic code. Indeed, it has been found that genetic code contains several layers of such additional information and that synonymous codons are strategically placed within mRNAs to ensure a particular translation kinetics facilitating and fine-tuning co-translational protein folding in the cell via step-wise/sequential structuring of distinct regions of the polypeptide chain emerging from the ribosome at different points in time. This review summarizes key findings in the field that have identified the role of synonymous codons and their usage in protein folding in the cell.

The impact of codon choice on translation process in Saccharomyces cerevisiae: folding class, protein function and secondary structure

The genetic code consists in a set of rules used by living organisms to translate genomic information, contained in genes, into proteins; every amino acid is coded by a set of nucleotide triplets or codons. We refer to codon choice as the choice of a given codon, among the synonymous available ones, to code a given amino acid occurrence. The aim of this work is to shed light on the pivotal role that codon choice plays in regulating the timing of translation process, through patterns of low and high translation efficiency codons. A translation efficiency value, namely codon score, was associated to each codon through a formula based on the number of tRNAs gene copies able to translate the given codon. By using codon scores, those k-mers of the proteome of Saccharomyces cerevisiae, showing low and high average scores associated to the correspondent codons, were computed. The analysis of distribution of both low and high average score k-mers clearly showed that, in particular for higher k-mer size, they occur much more than expected, strongly suggesting a functional role. Moreover performed analysis highlighted that significant k-mers preferentially occur in some protein folding classes, such as those containing alpha helices, and in some functional classes mainly involved in transcription process while codon choice seems to have a very low impact in proteins associated to energy production and metabolism. The relationship between secondary structures and significant k-mers was investigated, revealing that low score k-mers tend to preferentially occur in coil or close to coil regions and almost never in beta sheets, while high score k-mers preferentially occur in alpha helices, avoiding beta sheets, and close to coil regions for high k-mer sizes. Finally the analysis of distribution of significant codon patterns along the proteins highlighted a relevant enrichment of low average score k-mers at the 5' end of protein-coding sequences in the region from 5th to 25th amino acid.

Mechanisms of Cotranslational Protein Maturation in Bacteria

Growing cells invest a significant part of their biosynthetic capacity into the production of proteins. To become functional, newly-synthesized proteins must be N-terminally processed, folded and often translocated to other cellular compartments. A general strategy is to integrate these protein maturation processes with translation, by cotranslationally engaging processing enzymes, chaperones and targeting factors with the nascent polypeptide. Precise coordination of all factors involved is critical for the efficiency and accuracy of protein synthesis and cellular homeostasis. This review provides an overview of the current knowledge on cotranslational protein maturation, with a focus on the production of cytosolic proteins in bacteria. We describe the role of the ribosome and the chaperone network in protein folding and how the dynamic interplay of all cotranslationally acting factors guides the sequence of cotranslational events. Finally, we discuss recent data demonstrating the coupling of protein synthesis with the assembly of protein complexes and end with a brief discussion of outstanding questions and emerging concepts in the field of cotranslational protein maturation.

Codon Bias Can Determine Sorting of a Potassium Channel Protein

Due to the redundancy of the genetic code most amino acids are encoded by multiple synonymous codons. It has been proposed that a biased frequency of synonymous codons can affect the function of proteins by modulating distinct steps in transcription, translation and folding. Here, we use two similar prototype K⁺ channels as model systems to examine whether codon choice has an impact on protein sorting. By monitoring transient expression of GFP-tagged channels in mammalian cells, we find that one of the two channels is sorted in a codon and cell cycle-dependent manner either to mitochondria or the secretory pathway. The data establish that a gene with either rare or frequent codons serves, together with a cell-state-dependent decoding mechanism, as a secondary code for sorting intracellular membrane proteins.

Assessing optimal: inequalities in codon optimization algorithms

BACKGROUND

Custom genes have become a common resource in recombinant biology over the last 20 years due to the plummeting cost of DNA synthesis. These genes are often "optimized" to non-native sequences for overexpression in a non-native host by substituting synonymous codons within the coding DNA sequence (CDS). A handful of studies have compared native and optimized CDSs, reporting different levels of soluble product due to the accumulation of misfolded aggregates, variable activity of enzymes, and (at least one report of) a change in substrate specificity. No study, to the best of our knowledge, has performed a practical comparison of CDSs generated from different codon optimization algorithms or reported the corresponding protein yields.

RESULTS

In our efforts to understand what factors constitute an optimized CDS, we identified that there is little consensus among codon-optimization algorithms, a roughly equivalent chance that an algorithm-optimized CDS will increase or diminish recombinant yields as compared to the native DNA, a near ubiquitous use of a codon database that was last updated in 2007, and a high variability of output CDSs by some algorithms. We present a case study, using KRas4B, to demonstrate that a median codon frequency may be a better predictor of soluble yields than the more commonly utilized CAI metric.

CONCLUSIONS

We present a method for visualizing, analyzing, and comparing algorithm-optimized DNA sequences for recombinant protein expression. We encourage researchers to consider if DNA optimization is right for their experiments, and work towards improving the reproducibility of published recombinant work by publishing non-native CDSs.

Potential role of the X circular code in the regulation of gene expression

The X circular code is a set of 20 trinucleotides (codons) that has been identified in the protein-coding genes of most organisms (bacteria, archaea, eukaryotes, plasmids, viruses). It has been shown previously that the X circular code has the important mathematical property of being an error-correcting code. Thus, motifs of the X circular code, i.e. a series of codons belonging to X and called X motifs, allow identification and maintenance of the reading frame in genes. X motifs are significantly enriched in protein-coding genes, but have also been identified in many transfer RNA (tRNA) genes and in important functional regions of the ribosomal RNA (rRNA), notably in the peptidyl transferase center and the decoding center. Here, we investigate the potential role of X motifs as functional elements of protein-coding genes. First, we identify the codons of the X circular code which are frequent or rare in each domain of life (archaea, bacteria, eukaryota) and show that, for the amino acids with the highest codon bias, the preferred codon is often an X codon. We also observe a correlation between the 20 X codons and the optimal codons/dicodons that have been shown to influence translation efficiency. Then, we examined recently published experimental results concerning gene expression levels in diverse organisms. The approach used is the analysis of X motifs according to their density d_s(X), i.e. the number of X motifs per kilobase in a gene sequence s. Surprisingly, this simple parameter identifies several unexpected relations between the X circular code and gene expression. For example, the X motifs are significantly enriched in the minimal gene set belonging to the three domains of life, and in codon-optimized genes. Furthermore, the density of X motifs generally correlates with experimental measures of translation efficiency and mRNA stability. Taken together, these results lead us to propose that the X motifs may represent a genetic signal contributing to the maintenance of the correct reading frame and the optimization and regulation of gene expression.

Translational Control by Ribosome Pausing in Bacteria: How a Non-uniform Pace of Translation Affects Protein Production and Folding

Protein homeostasis of bacterial cells is maintained by coordinated processes of protein production, folding, and degradation. Translational efficiency of a given mRNA depends on how often the ribosomes initiate synthesis of a new polypeptide and how quickly they read the coding sequence to produce a full-length protein. The pace of ribosomes along the mRNA is not uniform: periods of rapid synthesis are separated by pauses. Here, we summarize recent evidence on how ribosome pausing affects translational efficiency and protein folding. We discuss the factors that slow down translation elongation and affect the quality of the newly synthesized protein. Ribosome pausing emerges as important factor contributing to the regulatory programs that ensure the quality of the proteome and integrate the cellular and environmental cues into regulatory circuits of the cell.

Disome-seq reveals widespread ribosome collisions that promote cotranslational protein folding

BACKGROUND

The folding of proteins is challenging in the highly crowded and sticky environment of a cell. Regulation of translation elongation may play a crucial role in ensuring the correct folding of proteins. Much of our knowledge regarding translation elongation comes from the sequencing of mRNA fragments protected by single ribosomes by ribo-seq. However, larger protected mRNA fragments have been observed, suggesting the existence of an alternative and previously hidden layer of regulation.

RESULTS

In this study, we performed disome-seq to sequence mRNA fragments protected by two stacked ribosomes, a product of translational pauses during which the 5'-elongating ribosome collides with the 3'-paused one. We detected widespread ribosome collisions that are related to slow ribosome release when stop codons are at the A-site, slow peptide bond formation from proline, glycine, asparagine, and cysteine when they are at the P-site, and slow leaving of polylysine from the exit tunnel of ribosomes. The structure of disomes obtained by cryo-electron microscopy suggests a different conformation from the substrate of the ribosome-associated protein quality control pathway. Collisions occurred more frequently in the gap regions between α-helices, where a translational pause can prevent the folding interference from the downstream peptides. Paused or collided ribosomes are associated with specific chaperones, which can aid in the cotranslational folding of the nascent peptides.

CONCLUSIONS

Therefore, cells use regulated ribosome collisions to ensure protein homeostasis.

Stop codon read-through of mammalian MTCH2 leading to an unstable isoform regulates mitochondrial membrane potential

Stop codon read-through (SCR) is a process of continuation of translation beyond a stop codon. This phenomenon, which occurs only in certain mRNAs under specific conditions, leads to a longer isoform with properties different from that of the canonical isoform. MTCH2, which encodes a mitochondrial protein that regulates mitochondrial metabolism, was selected as a potential read-through candidate based on evolutionary conservation observed in the proximal region of its 3' UTR. Here, we demonstrate translational read-through across two evolutionarily conserved, in-frame stop codons of MTCH2 using luminescence- and fluorescence-based assays, and by analyzing ribosome-profiling and mass spectrometry (MS) data. This phenomenon generates two isoforms, MTCH2x and MTCH2xx (single- and double-SCR products, respectively), in addition to the canonical isoform MTCH2, from the same mRNA. Our experiments revealed that a cis-acting 12-nucleotide sequence in the proximal 3' UTR of MTCH2 is the necessary signal for SCR. Functional characterization showed that MTCH2 and MTCH2x were localized to mitochondria with a long t1/2 (>36 h). However, MTCH2xx was found predominantly in the cytoplasm. This mislocalization and its unique C terminus led to increased degradation, as shown by greatly reduced t1/2 (<1 h). MTCH2 read-through-deficient cells, generated using CRISPR-Cas9, showed increased MTCH2 expression and, consistent with this, decreased mitochondrial membrane potential. Thus, double-SCR of MTCH2 regulates its own expression levels contributing toward the maintenance of normal mitochondrial membrane potential.

Development of monoclonal antibodies against recombinant LipL21 protein of pathogenic Leptospira through phage display technology

Leptospirosis is a potentially fatal zoonosis that is caused by spirochete Leptospira. The signs and symptoms of leptospirosis are usually varied, allowing it to be mistaken for other causes of acute febrile syndromes. Thus, early diagnosis and identification of a specific agent in clinical samples is crucial for effective treatment. This study was aimed to develop specific monoclonal antibodies against LipL21 antigen for future use in leptospirosis rapid and accurate immunoassay. A recombinant LipL21 (rLipL21) antigen was optimized for expression and evaluated for immunogenicity. Then, a naïve phage antibody library was utilized to identify single chain fragment variable (scFv) clones against the rLipL21 antigen. A total of 47 clones were analysed through monoclonal phage ELISA. However, after taking into consideration the background OD₄₀₅ values, only 4 clones were sent for sequencing to determine human germline sequences. The sequence analysis showed that all 4 clones are identical. The in silico analysis of scFv-lip-1 complex indicated that the charged residues of scFv CDRs are responsible for the recognition with rLipL21 epitopes. The generated monoclonal antibody against rLipL21 will be evaluated as a detection reagent for the diagnosis of human leptospirosis in a future study.

Network analysis of synonymous codon usage

MOTIVATION

Most amino acids are encoded by multiple synonymous codons, some of which are used more rarely than others. Analyses of positions of such rare codons in protein sequences revealed that rare codons can impact co-translational protein folding and that positions of some rare codons are evolutionarily conserved. Analyses of their positions in protein 3-dimensional structures, which are richer in biochemical information than sequences alone, might further explain the role of rare codons in protein folding.

RESULTS

We model protein structures as networks and use network centrality to measure the structural position of an amino acid. We first validate that amino acids buried within the structural core are network-central, and those on the surface are not. Then, we study potential differences between network centralities and thus structural positions of amino acids encoded by conserved rare, non-conserved rare and commonly used codons. We find that in 84% of proteins, the three codon categories occupy significantly different structural positions. We examine protein groups showing different codon centrality trends, i.e. different relationships between structural positions of the three codon categories. We see several cases of all proteins from our data with some structural or functional property being in the same group. Also, we see a case of all proteins in some group having the same property. Our work shows that codon usage is linked to the final protein structure and thus possibly to co-translational protein folding.

AVAILABILITY AND IMPLEMENTATION

https://nd.edu/∼cone/CodonUsage/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Codon bias analyses on thyroid carcinoma genes

BACKGROUND

Thyroid carcinoma is one of the most common cancers in the world. Although the genetics of thyroid carcinoma was intensively studied, new mechanisms could be involved in its development as the codon bias. In this paper, we studied the codon bias of thyroid-cancer genes, considering not only the sequences but also the synonymous mutations.

METHODS

Different measures and statistical analyses were employed to characterize the thyroid-cancer genes. We considered classical measures as RSCU and ENC, the compositional and protein characteristics, but also the codon bias landscape via the %MinMax algorithm.

RESULTS

The compositional analyses highlighted two groups of thyroid cancer genes according to the GC% and GC3% content. The ENC did not show a clear codon bias in the genes. Differently, the RSCU analyses showed interesting codons that could play an important role in the development of thyroid cancer as the codon Ser-tcG. Furthermore, interesting synonymous mutations were detected that could affect the codon bias. The codon bias landscape detected genes enriched in rare codons as AKAP9 and KTN1. A cluster analysis based on %MinMax classified the thyroid cancer genes in four different groups according to the distribution of rare/frequent codons in the sequence.

CONCLUSIONS

This is the first study that analyzed the codon bias in thyroid cancer genes based also on synonymous mutations. This study provided different hints that should be further investigated by wet-lab validation and that it could open new scenarios in the understanding the molecular mechanisms involved in thyroid cancer development based on codon bias.

Synthesis runs counter to directional folding of a nascent protein domain

Folding of individual domains in large proteins during translation helps to avoid otherwise prevalent inter-domain misfolding. How folding intermediates observed in vitro for the majority of proteins relate to co-translational folding remains unclear. Combining in vivo and single-molecule experiments, we followed the co-translational folding of the G-domain, encompassing the first 293 amino acids of elongation factor G. Surprisingly, the domain remains unfolded until it is fully synthesized, without collapsing into molten globule-like states or forming stable intermediates. Upon fully emerging from the ribosome, the G-domain transitions to its stable native structure via folding intermediates. Our results suggest a strictly sequential folding pathway initiating from the C-terminus. Folding and synthesis thus proceed in opposite directions. The folding mechanism is likely imposed by the final structure and might have evolved to ensure efficient, timely folding of a highly abundant and essential protein.

In vivo experiments and optical tweezers force-spectroscopy measurements assessing the co-translational folding of the G-domain from bacterial elongation factor G reveal a sequential folding pathway initiating from the C-terminus. These results suggest that protein folding and synthesis proceed in opposite directions.

Analysis of computational codon usage models and their association with translationally slow codons

Improved computational modeling of protein translation rates, including better prediction of where translational slowdowns along an mRNA sequence may occur, is critical for understanding co-translational folding. Because codons within a synonymous codon group are translated at different rates, many computational translation models rely on analyzing synonymous codons. Some models rely on genome-wide codon usage bias (CUB), believing that globally rare and common codons are the most informative of slow and fast translation, respectively. Others use the CUB observed only in highly expressed genes, which should be under selective pressure to be translated efficiently (and whose CUB may therefore be more indicative of translation rates). No prior work has analyzed these models for their ability to predict translational slowdowns. Here, we evaluate five models for their association with slowly translated positions as denoted by two independent ribosome footprint (RFP) count experiments from S. cerevisiae, because RFP data is often considered as a “ground truth” for translation rates across mRNA sequences. We show that all five considered models strongly associate with the RFP data and therefore have potential for estimating translational slowdowns. However, we also show that there is a weak correlation between RFP counts for the same genes originating from independent experiments, even when their experimental conditions are similar. This raises concerns about the efficacy of using current RFP experimental data for estimating translation rates and highlights a potential advantage of using computational models to understand translation rates instead.

Codon Harmonization of a Kir3.1-KirBac1.3 Chimera for Structural Study Optimization

The expression of functional, folded, and isotopically enriched membrane proteins is an enduring bottleneck for nuclear magnetic resonance (NMR) studies. Indeed, historically, protein yield optimization has been insufficient to allow NMR analysis of many complex Eukaryotic membrane proteins. However, recent work has found that manipulation of plasmid codons improves the odds of successful NMR-friendly protein production. In the last decade, numerous studies showed that matching codon usage patterns in recombinant gene sequences to those in the native sequence is positively correlated with increased protein yield. This phenomenon, dubbed codon harmonization, may be a powerful tool in optimizing recombinant expression of difficult-to-produce membrane proteins for structural studies. Here, we apply this technique to an inward rectifier K⁺ Channel (Kir) 3.1-KirBac1.3 chimera. Kir3.1 falls within the G protein-coupled inward rectifier K⁺ (GIRK) channel family, thus NMR studies may inform on the nuances of GIRK gating action in the presence and absence of its G Protein, lipid, and small molecule ligands. In our hands, harmonized plasmids increase protein yield nearly two-fold compared to the traditional ‘fully codon optimized’ construct. We then employ a fluorescence-based functional assay and solid-state NMR correlation spectroscopy to show the final protein product is folded and functional.

Synonymous codon substitutions perturb cotranslational protein folding in vivo and impair cell fitness

In the cell, proteins are synthesized from N to C terminus and begin to fold during translation. Cotranslational folding mechanisms are therefore linked to elongation rate, which varies as a function of synonymous codon usage. However, synonymous codon substitutions can affect many distinct cellular processes, which has complicated attempts to deconvolve the extent to which synonymous codon usage can promote or frustrate proper protein folding in vivo. Although previous studies have shown that some synonymous changes can lead to different final structures, other substitutions will likely be more subtle, perturbing predominantly the protein folding pathway without radically altering the final structure. Here we show that synonymous codon substitutions encoding a single essential enzyme lead to dramatically slower cell growth. These mutations do not prevent active enzyme formation; instead, they predominantly alter the protein folding mechanism, leading to enhanced degradation in vivo. These results support a model in which synonymous codon substitutions can impair cell fitness by significantly perturbing cotranslational protein folding mechanisms, despite the chaperoning provided by the cellular protein homeostasis network.

[Synonymous Codon Usage-a Guide for Co-Translational Protein Folding in the Cell]

In the cell, protein folding begins during protein synthesis/translation and thus is a co-translational process. Co-translational protein folding is tightly linked to translation elongation, which is not a uniform process. While there are many reasons for translation non-uniformity, it is generally believed that non-uniform synonymous codon usage is one of the key factors modulating translation elongation rates. Frequent/optimal codons as a rule are translated more rapidly than infrequently used ones and vice versa. Over 30 years ago, it was hypothesized that changes in synonymous codon usage affecting translation elongation rates could impinge on co-translation protein folding and that many synonymous codons are strategically placed within mRNA to ensure a particular translation kinetics facilitating productive step-by-step co-translational folding of proteins. It was suggested that this particular translation kinetics (and, specifically, translation pause sites) may define the window of opportunity for the protein parts to fold locally, particularly at the critical points where folding is far from equilibrium. It was thus hypothesized that synonymous codons may provide a secondary code for protein folding in the cell. Although, mostly accepted now, this hypothesis appeared to be difficult to prove and many convincing results were obtained only relatively recently. Here, I review the progress in the field and explain, why this simple idea appeared to be so challenging to prove.

Conservation of location of several specific inhibitory codon pairs in the Saccharomyces sensu stricto yeasts reveals translational selection

Synonymous codons provide redundancy in the genetic code that influences translation rates in many organisms, in which overall codon use is driven by selection for optimal codons. It is unresolved if or to what extent translational selection drives use of suboptimal codons or codon pairs. In Saccharomyces cerevisiae, 17 specific inhibitory codon pairs, each comprised of adjacent suboptimal codons, inhibit translation efficiency in a manner distinct from their constituent codons, and many are translated slowly in native genes. We show here that selection operates within Saccharomyces sensu stricto yeasts to conserve nine of these codon pairs at defined positions in genes. Conservation of these inhibitory codon pairs is significantly greater than expected, relative to conservation of their constituent codons, with seven pairs more highly conserved than any other synonymous pair. Conservation is strongly correlated with slow translation of the pairs. Conservation of suboptimal codon pairs extends to two related Candida species, fungi that diverged from Saccharomyces ∼270 million years ago, with an enrichment for codons decoded by I•A and U•G wobble in both Candida and Saccharomyces. Thus, conservation of inhibitory codon pairs strongly implies selection for slow translation at particular gene locations, executed by suboptimal codon pairs.

Analysis and preliminary characterisation of the cytochrome P450 monooxygenases from Frankia sp. EuI1c (Frankia inefficax sp.)

Frankia bacteria are nitrogen fixing species from the Actinobacterium phylum which live on the root nodules of plants. They have been hypothesised to have significant potential for natural product biosynthesis. The cytochrome P450 monooxygenase complement of Frankia sp. EuI1c (Frankia inefficax sp.), which comprises 68 members, was analysed. Several members belonged to previously uncharacterised bacterial P450 families. There was an unusually high number of CYP189 family members (21) suggesting that this family has undergone gene duplication events which are classified as "blooms". The likely electron transfer partners for the P450 enzymes were also identified and analysed. These consisted of predominantly [3Fe-4S] cluster containing ferredoxins (eight), a single [2Fe-2S] ferredoxin and a couple of ferredoxin reductases. Three of these CYP family members were produced and purified, using Escherichia coli as a host, and their substrate range was characterised. CYP1027H1 and CYP150A20 bound a broad range of norisoprenoids and terpenoids. CYP1074A2 was highly specific for certain steroids including testosterone, progesterone, stanolone and 4-androstene-3,17-dione. It is likely that steroids are the physiological substrates of CYP1074A2. These results also give an indication that terpenoids are the likely substrates of CYP1027H1 and CYP150A2. The large number of P450s belonging to distinct families as well as the associated electron transfer partners found in different Frankia strains highlights the importance of this family of enzymes has in the secondary metabolism of these bacteria.