Solvent accessibility, residue charge and residue volume, the three ingredients of a robust amino acid substitution matrix

Fail-safe genetic codes designed to intrinsically contain engineered organisms

One challenge in engineering organisms is taking responsibility for their behavior over many generations. Spontaneous mutations arising before or during use can impact heterologous genetic functions, disrupt system integration, or change organism phenotype. Here, we propose restructuring the genetic code itself such that point mutations in protein-coding sequences are selected against. Synthetic genetic systems so-encoded should fail more safely in response to most spontaneous mutations. We designed fail-safe codes and simulated their expected effects on the evolution of so-encoded proteins. We predict fail-safe codes supporting expression of 20 or 15 amino acids could slow protein evolution to ∼30% or 0% the rate of standard-encoded proteins, respectively. We also designed quadruplet-codon codes that should ensure all single point mutations in protein-coding sequences are selected against while maintaining expression of 20 or more amino acids. We demonstrate experimentally that a reduced set of 21 tRNAs is capable of expressing a protein encoded by only 20 sense codons, whereas a standard 64-codon encoding is not expressed. Our work suggests that biological systems using rationally depleted but otherwise natural translation systems should evolve more slowly and that such hypoevolvable organisms may be less likely to invade new niches or outcompete native populations.

Global importance of RNA secondary structures in protein-coding sequences

Motivation

The protein-coding sequences of messenger RNAs are the linear template for translation of the gene sequence into protein. Nevertheless, the RNA can also form secondary structures by intramolecular base-pairing.

Results

We show that the nucleotide distribution within codons is biased in all taxa of life on a global scale. Thereby, RNA secondary structures that require base-pairing between the position 1 of a codon with the position 1 of an opposing codon (here named RNA secondary structure class c₁) are under-represented. We conclude that this bias may result from the co-evolution of codon sequence and mRNA secondary structure, suggesting that RNA secondary structures are generally important in protein-coding regions of mRNAs. The above result also implies that codon position 2 has a smaller influence on the amino acid choice than codon position 1.

Supplementary information

Supplementary data are available at Bioinformatics online.

A binary representation of the genetic code

This article introduces a novel binary representation of the canonical genetic code based on both the structural similarities of the nucleotides, as well as the physicochemical properties of the encoded amino acids. Each of the four mRNA bases is assigned a unique 2-bit identifier, so that the 64 triplet codons are each indexed by a 6-bit label. The ordering of the bits reflects the hierarchical organization manifested by the DNA replication/repair and tRNA translation systems. In this system, transition and transversion mutations are naturally expressed as binary operations, and the severities of the different point mutations can be analyzed. Using a principal component analysis, it is shown that the physicochemical properties of amino acids related to protein folding also correlate with certain bit positions of their respective labels. Thus, the likelihood for a point mutation to be conservative, and less likely to cause a change in protein functionality, can be estimated.

Principal components analysis of protein sequence clusters

Sequence analysis of large protein families can produce sub-clusters even within the same family. In some cases, it is of interest to know precisely which amino acid position variations are most responsible for driving separation into sub-clusters. In large protein families composed of large proteins, it can be quite challenging to assign the relative importance to specific amino acid positions. Principal components analysis (PCA) is ideal for such a task, since the problem is posed in a large variable space, i.e. the number of amino acids that make up the protein sequence, and PCA is powerful at reducing the dimensionality of complex problems by projecting the data into an eigenspace that represents the directions of greatest variation. However, PCA of aligned protein sequence families is complicated by the fact that protein sequences are traditionally represented by single letter alphabetic codes, whereas PCA of protein sequence families requires conversion of sequence information into a numerical representation. Here, we introduce a new amino acid sequence conversion algorithm optimized for PCA data input. The method is demonstrated using a small artificial dataset to illustrate the characteristics and performance of the algorithm, as well as a small protein sequence family consisting of nine members, COG2263, and finally with a large protein sequence family, Pfam04237, which contains more than 1,800 sequences that group into two sub-clusters.