A common periodic table of codons and amino acids

The Canonical Table of the Genetic Code as a periodic system of triplets

The Canonical Table of the Genetic Code (CTGC) is constructed theoretically on the basis of the similarity of PFs (PF) of proteins with the conformation of 4-arc chain graphs (Karasev, 2019). Of the 64 conformations of the graph, specified by the position of the connectivity edges, and the matrices of 6 variables (x₁ … x₆), x_i = (0, 1), 4 blocks of 16 elements each were formed. Then they were coded in the form of triplets based on the correspondence of pairs of variables to four letters of the code: 00 = C, 01 = U, 10 = G, 11 = A, and supplemented based on the known triplet-amino acid assignment. The resulting table is compared with the Periodic Table of Chemical Elements (PTCE). As in the PTCE, this CTGC has an initial element - a triplet that encodes graphs with zero number of connected edges. Within each block, vacancies are filled with connectivity edges in two alternative ways, both in rows and in the columns. As we move from the initial block 00 to the final block 11, there is a sequential filling of vacancies for variables x₃x₄: 00, 01, 10, 11. In general, the CTGC can be considered as a periodic system of triplets. Comparison with the previously described variety of tables of the genetic code made it possible to conclude that the CTGC more adequately reflects the properties of the genetic code. Prospects for the possible application of this table are being discussed.

Interplay between coding and exonic splicing regulatory sequences

The inclusion of exons during the splicing process depends on the binding of splicing factors to short low-complexity regulatory sequences. The relationship between exonic splicing regulatory sequences and coding sequences is still poorly understood. We demonstrate that exons that are coregulated by any given splicing factor share a similar nucleotide composition bias and preferentially code for amino acids with similar physicochemical properties because of the nonrandomness of the genetic code. Indeed, amino acids sharing similar physicochemical properties correspond to codons that have the same nucleotide composition bias. In particular, we uncover that the TRA2A and TRA2B splicing factors that bind to adenine-rich motifs promote the inclusion of adenine-rich exons coding preferentially for hydrophilic amino acids that correspond to adenine-rich codons. SRSF2 that binds guanine/cytosine-rich motifs promotes the inclusion of GC-rich exons coding preferentially for small amino acids, whereas SRSF3 that binds cytosine-rich motifs promotes the inclusion of exons coding preferentially for uncharged amino acids, like serine and threonine that can be phosphorylated. Finally, coregulated exons encoding amino acids with similar physicochemical properties correspond to specific protein features. In conclusion, the regulation of an exon by a splicing factor that relies on the affinity of this factor for specific nucleotide(s) is tightly interconnected with the exon-encoded physicochemical properties. We therefore uncover an unanticipated bidirectional interplay between the splicing regulatory process and its biological functional outcome.

Reciprocal Nucleopeptides as the Ancestral Darwinian Self-Replicator

Even the simplest organisms are too complex to have spontaneously arisen fully formed, yet precursors to first life must have emerged ab initio from their environment. A watershed event was the appearance of the first entity capable of evolution: the Initial Darwinian Ancestor. Here, we suggest that nucleopeptide reciprocal replicators could have carried out this important role and contend that this is the simplest way to explain extant replication systems in a mathematically consistent way. We propose short nucleic acid templates on which amino-acylated adapters assembled. Spatial localization drives peptide ligation from activated precursors to generate phosphodiester-bond-catalytic peptides. Comprising autocatalytic protein and nucleic acid sequences, this dynamical system links and unifies several previous hypotheses and provides a plausible model for the emergence of DNA and the operational code.

In Silico Molecular Modeling and Docking Studies on Novel Mutants (E229V, H225P and D230C) of the Nucleotide-Binding Domain of Homo sapiens Hsp70

In this study, we explored the possibility of determining the synergistic interactions between nucleotide-binding domain (NBD) of Homo sapiens heat-shock 70 kDa protein (Hsp70) and E1A 32 kDa of adenovirus serotype 5 motif (PNLVP) in the efficiency of killing of tumor cells in cancer treatment. At present, the protein interaction between NBD and PNLVP motif is still unknown, but believed to enhance the rate of virus replication in tumor cells. Three mutant models (E229V, H225P and D230C) were built and simulated, and their interactions with PNLVP motif were studied. The PNLVP motif showed the binding energy and intermolecular energy values with the novel E229V mutant at -7.32 and -11.2 kcal/mol. The E229V mutant had the highest number of hydrogen bonds (7). Based on the root mean square deviation, root mean square fluctuation, hydrogen bonds, salt bridge, secondary structure, surface-accessible solvent area, potential energy and distance matrices analyses, it was proved that the E229V had the strongest and most stable interaction with the PNLVP motif among all the four protein-ligand complex structures. The knowledge of this protein-ligand complex model would help in designing Hsp70 structure-based drug for cancer therapy.

Cross Talk between KGF and KITLG Proteins Implicated with Ovarian Folliculogenesis in Buffalo Bubalus bubalis

Molecular interactions between mesenchymal-derived Keratinocyte growth factor (KGF) and Kit ligand (KITLG) are essential for follicular development. These factors are expressed by theca and granulosa cells. We determined full length coding sequence of buffalo KGF and KITLG proteins having 194 and 274 amino acids, respectively. The recombinant KGF and KITLG proteins were solubilized in 10 mM Tris, pH 7.5 and 50 mM Tris, pH 7.4 and purified using Ni-NTA column and GST affinity chromatography, respectively. The purity and molecular weight of His-KGF (~23 kDa) and GST-KITLG (~57 kDa) proteins were confirmed by SDS-PAGE and western blotting. The co-immunoprecipitation assay accompanied with computational analysis demonstrated the interaction between KGF and KITLG proteins. We deduced 3D structures of the candidate proteins and assessed their binding based on protein docking. In the process, KGF specific residues, Lys123, Glu135, Lys140, Lys155 and Trp156 and KITLG specific ones, Ser226, Phe233, Gly234, Ala235, Phe236, Trp238 and Lys239 involved in the formation of KGF-KITLG complex were detected. The hydrophobic interactions surrounding KGF-KITLG complex affirmed their binding affinity and stability to the interacting interface. Additionally, in-silico site directed mutagenesis enabled the assessment of changes that occurred in the binding energies of mutated KGF-KITLG protein complex. Our results demonstrate that in the presence of KITLG, KGF mimics its native binding mode suggesting all the KGF residues are specific to their binding complex. This study provides an insight on the critical amino acid residues participating in buffalo ovarian folliculogenesis.

Modeling and docking studies on novel mutants (K71L and T204V) of the ATPase domain of human heat shock 70 kDa protein 1

The purpose of exploring protein interactions between human adenovirus and heat shock protein 70 is to exploit a potentially synergistic interaction to enhance anti-tumoral efficacy and decrease toxicity in cancer treatment. However, the protein interaction of Hsp70 with E1A32 kDa of human adenovirus serotype 5 remains to be elucidated. In this study, two residues of ATPase domain of human heat shock 70 kDa protein 1 (PDB: 1 HJO) were mutated. 3D mutant models (K71L and T204V) using PyMol software were then constructed. The structures were evaluated by PROCHECK, ProQ, ERRAT, Verify 3D and ProSA modules. All evidence suggests that all protein models are acceptable and of good quality. The E1A32 kDa motif was retrieved from UniProt (P03255), as well as subjected to docking interaction with NBD, K71L and T204V, using the Autodock 4.2 program. The best lowest binding energy value of −9.09 kcal/mol was selected for novel T204V. Moreover, the protein-ligand complex structures were validated by RMSD, RMSF, hydrogen bonds and salt bridge analysis. This revealed that the T204V-E1A32 kDa motif complex was the most stable among all three complex structures. This study provides information about the interaction between Hsp70 and the E1A32 kDa motif, which emphasizes future perspectives to design rational drugs and vaccines in cancer therapy.

Principles of self-organization in biological pathways: a hypothesis on the autogenous association of alpha-synuclein

Previous evidence indicates that a number of proteins are able to interact with cognate mRNAs. These autogenous associations represent important regulatory mechanisms that control gene expression at the translational level. Using the catRAPID approach to predict the propensity of proteins to bind to RNA, we investigated the occurrence of autogenous associations in the human proteome. Our algorithm correctly identified binding sites in well-known cases such as thymidylate synthase, tumor suppressor P53, synaptotagmin-1, serine/ariginine-rich splicing factor 2, heat shock 70 kDa, ribonucleic particle-specific U1A and ribosomal protein S13. In addition, we found that several other proteins are able to bind to their own mRNAs. A large-scale analysis of biological pathways revealed that aggregation-prone and structurally disordered proteins have the highest propensity to interact with cognate RNAs. These findings are substantiated by experimental evidence on amyloidogenic proteins such as TAR DNA-binding protein 43 and fragile X mental retardation protein. Among the amyloidogenic proteins, we predicted that Parkinson’s disease-related α-synuclein is highly prone to interact with cognate transcripts, which suggests the existence of RNA-dependent factors in its function and dysfunction. Indeed, as aggregation is intrinsically concentration dependent, it is possible that autogenous interactions play a crucial role in controlling protein homeostasis.

Coding nucleic acids are chaperons for protein folding: a novel theory of protein folding

The arguments for nucleic acid chaperons are reviewed and three new lines of evidence are added. (1) It was found that amino acids encoded by codons in short nucleic acid loops frequently form turns and helices in the corresponding protein structures. (2) The amino acids encoded by partially complementary (1st and 3rd nucleotides) codons are more frequently co-located in the encoded proteins than expected by chance. (3) There are significant correlations between thermodynamic changes (ddG) caused by codon mutations in nucleic acids and the thermodynamic changes caused by the corresponding amino acid mutations in the encoded proteins. We conclude that the concept of the Proteomic Code and nucleic acid chaperons seems correct from the bioinformatics point of view, and we expect to see direct biochemical experiments and evidence in the near future.

The concept of RNA-assisted protein folding: the role of tRNA

We suggest that tRNA actively participates in the transfer of 3D information from mRNA to peptides--in addition to its well-known, "classical" role of translating the 3-letter RNA codes into the one letter protein code. The tRNA molecule displays a series of thermodynamically favored configurations during translation, a movement which places the codon and coded amino acids in proximity to each other and make physical contact between some amino acids and their codons possible. This specific codon-amino acid interaction of some selected amino acids is necessary for the transfer of spatial information from mRNA to coded proteins, and is known as RNA-assisted protein folding.

Codon Deviation Coefficient: a novel measure for estimating codon usage bias and its statistical significance

Background

Genetic mutation, selective pressure for translational efficiency and accuracy, level of gene expression, and protein function through natural selection are all believed to lead to codon usage bias (CUB). Therefore, informative measurement of CUB is of fundamental importance to making inferences regarding gene function and genome evolution. However, extant measures of CUB have not fully accounted for the quantitative effect of background nucleotide composition and have not statistically evaluated the significance of CUB in sequence analysis.

Results

Here we propose a novel measure--Codon Deviation Coefficient (CDC)--that provides an informative measurement of CUB and its statistical significance without requiring any prior knowledge. Unlike previous measures, CDC estimates CUB by accounting for background nucleotide compositions tailored to codon positions and adopts the bootstrapping to assess the statistical significance of CUB for any given sequence. We evaluate CDC by examining its effectiveness on simulated sequences and empirical data and show that CDC outperforms extant measures by achieving a more informative estimation of CUB and its statistical significance.

Conclusions

As validated by both simulated and empirical data, CDC provides a highly informative quantification of CUB and its statistical significance, useful for determining comparative magnitudes and patterns of biased codon usage for genes or genomes with diverse sequence compositions.

On the organizational dynamics of the genetic code

The organization of the canonical genetic code needs to be thoroughly illuminated. Here we reorder the four nucleotides—adenine, thymine, guanine and cytosine—according to their emergence in evolution, and apply the organizational rules to devising an algebraic representation for the canonical genetic code. Under a framework of the devised code, we quantify codon and amino acid usages from a large collection of 917 prokaryotic genome sequences, and associate the usages with its intrinsic structure and classification schemes as well as amino acid physicochemical properties. Our results show that the algebraic representation of the code is structurally equivalent to a content-centric organization of the code and that codon and amino acid usages under different classification schemes were correlated closely with GC content, implying a set of rules governing composition dynamics across a wide variety of prokaryotic genome sequences. These results also indicate that codons and amino acids are not randomly allocated in the code, where the six-fold degenerate codons and their amino acids have important balancing roles for error minimization. Therefore, the content-centric code is of great usefulness in deciphering its hitherto unknown regularities as well as the dynamics of nucleotide, codon, and amino acid compositions.

Modeling compositional dynamics based on GC and purine contents of protein-coding sequences

Background

Understanding the compositional dynamics of genomes and their coding sequences is of great significance in gaining clues into molecular evolution and a large number of publically-available genome sequences have allowed us to quantitatively predict deviations of empirical data from their theoretical counterparts. However, the quantification of theoretical compositional variations for a wide diversity of genomes remains a major challenge.

Results

To model the compositional dynamics of protein-coding sequences, we propose two simple models that take into account both mutation and selection effects, which act differently at the three codon positions, and use both GC and purine contents as compositional parameters. The two models concern the theoretical composition of nucleotides, codons, and amino acids, with no prerequisite of homologous sequences or their alignments. We evaluated the two models by quantifying theoretical compositions of a large collection of protein-coding sequences (including 46 of Archaea, 686 of Bacteria, and 826 of Eukarya), yielding consistent theoretical compositions across all the collected sequences.

Conclusions

We show that the compositions of nucleotides, codons, and amino acids are largely determined by both GC and purine contents and suggest that deviations of the observed from the expected compositions may reflect compositional signatures that arise from a complex interplay between mutation and selection via DNA replication and repair mechanisms.

Reviewers

This article was reviewed by Zhaolei Zhang (nominated by Mark Gerstein), Guruprasad Ananda (nominated by Kateryna Makova), and Daniel Haft.

The rules of variation: amino acid exchange according to the rotating circular genetic code

General guidelines for the molecular basis of functional variation are presented while focused on the rotating circular genetic code and allowable exchanges that make it resistant to genetic diseases under normal conditions. The rules of variation, bioinformatics aids for preventative medicine, are: (1) same position in the four quadrants for hydrophobic codons, (2) same or contiguous position in two quadrants for synonymous or related codons, and (3) same quadrant for equivalent codons. To preserve protein function, amino acid exchange according to the first rule takes into account the positional homology of essential hydrophobic amino acids with every codon with a central uracil in the four quadrants, the second rule includes codons for identical, acidic, or their amidic amino acids present in two quadrants, and the third rule, the smaller, aromatic, stop codons, and basic amino acids, each in proximity within a 90 degree angle. I also define codifying genes and palindromati, CTCGTGCCGAATTCGGCACGAG.

RNA editing of 10 Didymium iridis mitochondrial genes and comparison with the homologous genes in Physarum polycephalum

Regions of the Didymium iridis mitochondrial genome were identified with similarity to typical mitochondrial genes; however, these regions contained numerous stop codons. We used RT-PCR and DNA sequencing to determine whether, through RNA editing, these regions were transcribed into mRNAs that could encode functional proteins. Ten putative gene regions were examined: atp1, atp6, atp8, atp9, cox1, cox2, cytb, nad4L, nad6, and nad7. The cDNA sequences of each gene could encode a functional mitochondrial protein that was highly conserved compared with homologous genes. The type of editing events and editing sequence features were very similar to those observed in the homologous genes of Physarum polycephalum, though the actual editing locations showed a variable degree of conservation. Edited sites were compared with encoded sites in D. iridis and P. polycephalum for all 10 genes. Edited sequence for a portion of the cox1 gene was available for six myxomycetes, which, when compared, showed a high degree of conservation at the protein level. Different types of editing events showed varying degrees of site conservation with C-to-U base changes being the least conserved. Several aspects of single C insertion editing events led to the preferential creation of hydrophobic amino acid codons that may help to minimize adverse effects on the resulting protein structure.

Genome wide exploration of the origin and evolution of amino acids

Background

Even after years of exploration, the terrestrial origin of bio-molecules remains unsolved and controversial. Today, observation of amino acid composition in proteins has become an alternative way for a global understanding of the mystery encoded in whole genomes and seeking clues for the origin of amino acids.

Results

In this study, we statistically monitored the frequencies of 20 alpha-amino acids in 549 taxa from three kingdoms of life: archaebacteria, eubacteria, and eukaryotes. We found that the amino acids evolved independently in these three kingdoms; but, conserved linkages were observed in two groups of amino acids, (A, G, H, L, P, Q, R, and W) and (F, I, K, N, S, and Y). Moreover, the amino acids encoded by GC-poor codons (F, Y, N, K, I, and M) were found to "lose" their usage in the development from single cell eukaryotic organisms like S. cerevisiae to H. sapiens, while the amino acids encoded by GC-rich codons (P, A, G, and W) were found to gain usage. These findings further support the co-evolution hypothesis of amino acids and genetic codes.

Conclusion

We proposed a new chronological order of the appearance of amino acids (L, A, V/E/G, S, I, K, T, R/D, P, N, F, Q, Y, M, H, W, C). Two conserved evolutionary paths of amino acids were also suggested: A→G→R→P and K→Y.

RNA editing in six mitochondrial ribosomal protein genes of Didymium iridis

Similarity searches with Didymium iridis mitochondrial genomic DNA identified six possible ribosomal protein-coding regions, however, each region contained stop codons that would need to be removed by RNA editing to produce functional transcripts. RT-PCR was used to amplify these regions from total RNA for cloning and sequencing. Six functional transcripts were verified for the following ribosomal protein genes: rpS12, rpS7, rpL2, rpS19, rpS3, and rpL16. The editing events observed, such as single C and U nucleotide insertions and a dinucleotide insertion, were consistent with previously observed editing patterns seen in D. iridis. Additionally, a new form of insertional editing, a single A insertion, was observed in a conserved region of the rpL16 gene. While the majority of codons created by editing specify hydrophobic amino acids, a greater proportion of the codons created in these hydrophilic ribosomal proteins called for positively charged amino acids in comparison to the previously characterized hydrophobic respiratory protein genes. This first report of edited soluble mitochondrial ribosomal proteins in myxomycetes expands upon the RNA editing patterns previously seen; there was: a greater proportion of created codons specifying positively charged amino acids, a shift in the codon position edited, and the insertion of single A nucleotides.

Editorial: hypotheses about protein folding--the proteomic code and wonderfolds

Theoretical biology journals can contribute in many ways to the progress of knowledge. They are particularly well-placed to encourage dialogue and debate about hypotheses addressing problematical areas of research. An online journal provides an especially useful forum for such debate because of the option of posting comments within days of the publication of a contentious article.

Discovery of proteomic code with mRNA assisted protein folding

The 3x redundancy of the Genetic Code is usually explained as a necessity to increase the mutation-resistance of the genetic information. However recent bioinformatical observations indicate that the redundant Genetic Code contains more biological information than previously known and which is additional to the 64/20 definition of amino acids. It might define the physico-chemical and structural properties of amino acids, the codon boundaries, the amino acid co-locations (interactions) in the coded proteins and the free folding energy of mRNAs. This additional information, which seems to be necessary to determine the 3D structure of coding nucleic acids as well as the coded proteins, is known as the Proteomic Code and mRNA Assisted Protein Folding.

Does codon bias have an evolutionary origin?

BACKGROUND

There is a 3-fold redundancy in the Genetic Code; most amino acids are encoded by more than one codon. These synonymous codons are not used equally; there is a Codon Usage Bias (CUB). This article will provide novel information about the origin and evolution of this bias.

RESULTS

Codon Usage Bias (CUB, defined here as deviation from equal usage of synonymous codons) was studied in 113 species. The average CUB was 29.3 +/- 1.1% (S.E.M, n = 113) of the theoretical maximum and declined progressively with evolution and increasing genome complexity. A Pan-Genomic Codon Usage Frequency (CUF) Table was constructed to describe genome-wide relationships among codons. Significant correlations were found between the number of synonymous codons and (i) the frequency of the respective amino acids (ii) the size of CUB. Numerous, statistically highly significant, internal correlations were found among codons and the nucleic acids they comprise. These strong correlations made it possible to predict missing synonymous codons (wobble bases) reliably from the remaining codons or codon residues.

CONCLUSION

The results put the concept of "codon bias" into a novel perspective. The internal connectivity of codons indicates that all synonymous codons might be integrated parts of the Genetic Code with equal importance in maintaining its functional integrity.

The Proteomic Code: a molecular recognition code for proteins

Background

The Proteomic Code is a set of rules by which information in genetic material is transferred into the physico-chemical properties of amino acids. It determines how individual amino acids interact with each other during folding and in specific protein-protein interactions. The Proteomic Code is part of the redundant Genetic Code.

Review

The 25-year-old history of this concept is reviewed from the first independent suggestions by Biro and Mekler, through the works of Blalock, Root-Bernstein, Siemion, Miller and others, followed by the discovery of a Common Periodic Table of Codons and Nucleic Acids in 2003 and culminating in the recent conceptualization of partial complementary coding of interacting amino acids as well as the theory of the nucleic acid-assisted protein folding.

Methods and conclusions

A novel cloning method for the design and production of specific, high-affinity-reacting proteins (SHARP) is presented. This method is based on the concept of proteomic codes and is suitable for large-scale, industrial production of specifically interacting peptides.

Codes in the codons: construction of a codon/amino acid periodic table and a study of the nature of specific nucleic acid-protein interactions

The theory of "codon-amino acid coevolution" was first proposed by Woese in 1967. It suggests that there is a stereochemical matching - that is, affinity - between amino acids and certain of the base triplet sequences that code for those amino acids. We have constructed a common periodic table of codons and amino acids, where the nucleic acid table showed perfect axial symmetry for codons and the corresponding amino acid table also displayed periodicity regarding the biochemical properties (charge and hydrophobicity) of the 20 amino acids and the position of the stop signals. The table indicates that the middle (2/sup nd/) amino acid in the codon has a prominent role in determining some of the structural features of the amino acids. The possibility that physical contact between codons and amino acids might exist was tested on restriction enzymes. Many recognition site-like sequences were found in the coding sequences of these enzymes and as many as 73 examples of codon-amino acid co-location were observed in the 7 known 3D structures (December 2003) of endonuclease-nucleic acid complexes. These results indicate that the smallest possible units of specific nucleic acid-protein interaction are indeed the stereochemically compatible codons and amino acids.

Protein folding information in nucleic acids which is not present in the genetic code

Nucleic acid subsequences comprising the 1st and/or 3rd codon residues in mRNAs express significantly higher free folding energy (FFE) than the subsequence containing only the 2nd residues (P < 0.0001, n = 81). This periodic FFE difference is not present in introns. The FFE in the 1st and 3rd residues is additive, which suggests that these residues contain a significant number of complementary bases and contribute to selection for local mRNA secondary structures. This periodic, codon-related structure forming of mRNAs indicates a connection between the structure of exons and the corresponding (translated) proteins. The folding energy dot plots of RNAs and the residue contact maps of the coded proteins are indeed similar. Residue contact statistics using 81 different protein structures confirmed that amino acids that are coded by partially reverse and complementary codons (Watson-Crick base pairs at the 1st and 3rd codon positions and translated in reverse orientation) are preferentially co-located in protein structures.

Neuronal Control of Skin Function: The Skin as a Neuroimmunoendocrine Organ

This review focuses on the role of the peripheral nervous system in cutaneous biology and disease. During the last few years, a modern concept of an interactive network between cutaneous nerves, the neuroendocrine axis, and the immune system has been established. We learned that neurocutaneous interactions influence a variety of physiological and pathophysiological functions, including cell growth, immunity, inflammation, pruritus, and wound healing. This interaction is mediated by primary afferent as well as autonomic nerves, which release neuromediators and activate specific receptors on many target cells in the skin. A dense network of sensory nerves releases neuropeptides, thereby modulating inflammation, cell growth, and the immune responses in the skin. Neurotrophic factors, in addition to regulating nerve growth, participate in many properties of skin function. The skin expresses a variety of neurohormone receptors coupled to heterotrimeric G proteins that are tightly involved in skin homeostasis and inflammation. This neurohormone-receptor interaction is modulated by endopeptidases, which are able to terminate neuropeptide-induced inflammatory or immune responses. Neuronal proteinase-activated receptors or transient receptor potential ion channels are recently described receptors that may have been important in regulating neurogenic inflammation, pain, and pruritus. Together, a close multidirectional interaction between neuromediators, high-affinity receptors, and regulatory proteases is critically involved to maintain tissue integrity and regulate inflammatory responses in the skin. A deeper understanding of cutaneous neuroimmunoendocrinology may help to develop new strategies for the treatment of several skin diseases.

Indications that "codon boundaries" are physico-chemically defined and that protein-folding information is contained in the redundant exon bases

Background

All the information necessary for protein folding is supposed to be present in the amino acid sequence. It is still not possible to provide specific ab initio structure predictions by bioinformatical methods. It is suspected that additional folding information is present in protein coding nucleic acid sequences, but this is not represented by the known genetic code.

Results

Nucleic acid subsequences comprising the 1st and/or 3rd codon residues in mRNAs express significantly higher free folding energy (FFE) than the subsequence containing only the 2nd residues (p < 0.0001, n = 81). This periodic FFE difference is not present in introns. It is therefore a specific physico-chemical characteristic of coding sequences and might contribute to unambiguous definition of codon boundaries during translation. The FFEs of the 1st and 3rd residues are additive, which suggests that these residues contain a significant number of complementary bases and that may contribute to selection for local RNA secondary structures in coding regions. This periodic, codon-related structure-formation of mRNAs indicates a connection between the structures of exons and the corresponding (translated) proteins. The folding energy dot plots of RNAs and the residue contact maps of the coded proteins are indeed similar. Residue contact statistics using 81 different protein structures confirmed that amino acids that are coded by partially reverse and complementary codons (Watson-Crick (WC) base pairs at the 1st and 3rd codon positions and translated in reverse orientation) are preferentially co-located in protein structures.

Conclusion

Exons are distinguished from introns, and codon boundaries are physico-chemically defined, by periodically distributed FFE differences between codon positions. There is a selection for local RNA secondary structures in coding regions and this nucleic acid structure resembles the folding profiles of the coded proteins. The preferentially (specifically) interacting amino acids are coded by partially complementary codons, which strongly supports the connection between mRNA and the corresponding protein structures and indicates that there is protein folding information in nucleic acids that is not present in the genetic code. This might suggest an additional explanation of codon redundancy.

Amino acid size, charge, hydropathy indices and matrices for protein structure analysis

BACKGROUND

Prediction of protein folding and specific interactions from only the sequence (ab initio) is a major challenge in bioinformatics. It is believed that such prediction will prove possible if Anfinsen's thermodynamic principle is correct for all kinds of proteins, and all the information necessary to form a concrete 3D structure is indeed present in the sequence.

RESULTS

We indexed the 200 possible amino acid pairs for their compatibility regarding the three major physicochemical properties--size, charge and hydrophobicity--and constructed Size, Charge and Hydropathy Compatibility Indices and Matrices (SCI & SCM, CCI & CCM, and HCI & HCM). Each index characterized the expected strength of interaction (compatibility) of two amino acids by numbers from 1 (not compatible) to 20 (highly compatible). We found statistically significant positive correlations between these indices and the propensity for amino acid co-locations in real protein structures (a sample containing total 34630 co-locations in 80 different protein structures): for HCI: p < 0.01, n = 400 in 10 subgroups; for SCI p < 1.3E-08, n = 400 in 10 subgroups; for CCI: p < 0.01, n = 175). Size compatibility between residues (well known to exist in nucleic acids) is a novel observation for proteins. Regression analyzes indicated at least 7 well distinguished clusters regarding size compatibility and 5 clusters of charge compatibility. We tried to predict or reconstruct simple 2D representations of 3D structures from the sequence using these matrices by applying a dot plot-like method. The location and pattern of the most compatible subsequences was very similar or identical when the three fundamentally different matrices were used, which indicates the consistency of physicochemical compatibility. However, it was not sufficient to choose one preferred configuration between the many possible predicted options.

CONCLUSION

Indexing of amino acids for major physico-chemical properties is a powerful approach to understanding and assisting protein design. However, it is probably insufficient itself for complete ab initio structure prediction.

A novel intra-molecular protein–protein interaction code based on partial complementary coding of co-locating amino acids

Proteins are assumed to contain all the information necessary for unambiguous folding and specific interaction with each other. However, ab initio structure prediction is often not successful because the amino acid sequence itself is simply not sufficient to guide between endless folding possibilities. It seems to be logical to try to find the "missing" information in nucleic acids, in the redundant codon. Statistical analyses of approximately 35K amino acid co-locations in 80 different protein structures indicate the existence of a weak intra-molecular protein-protein interaction code. Co-locating amino acids are preferentially coded by codons which are complementary in reverse orientation to each other at the 1st and 3rd codon positions, but not necessarily at the 2nd. This code, called D-1 X 3/RC-3 X 1, limits the number of preferred amino acid pairs from 20 to 10.3+/-0.8 (SEM, n=20) and emphasizes the importance of "strictly" defined amino acids (those having less synonymous codons). The existence of this code does not by any means violate the known physicochemical rules of protein folding or interaction. It is suggested that the biological source of preferential (specific) amino acid co-locations is the partial complementarity of their codons. This special coding of co-locating amino acids is important to better understanding of some fundamental biochemical processes and observations such as: (a) protein folding; (b) specific and high affinity protein-protein interactions; (c) the role of the wobble bases; (d) the significance of the redundant genetic code; (e) the origin of specific protein-protein interactions. Furthermore it might be useful even in protein design.

Mitochondrial-encoded membrane protein transcripts are pyrimidine-rich while soluble protein transcripts and ribosomal RNA are purine-rich

BACKGROUND

Eukaryotic organisms contain mitochondria, organelles capable of producing large amounts of ATP by oxidative phosphorylation. Each cell contains many mitochondria with many copies of mitochondrial DNA in each organelle. The mitochondrial DNA encodes a small but functionally critical portion of the oxidative phosphorylation machinery, a few other species-specific proteins, and the rRNA and tRNA used for the translation of these transcripts. Because the microenvironment of the mitochondrion is unique, mitochondrial genes may be subject to different selectional pressures than those affecting nuclear genes.

RESULTS

From an analysis of the mitochondrial genomes of a wide range of eukaryotic species we show that there are three simple rules for the pyrimidine and purine abundances in mitochondrial DNA transcripts. Mitochondrial membrane protein transcripts are pyrimidine rich, rRNA transcripts are purine-rich and the soluble protein transcripts are purine-rich. The transitions between pyrimidine and purine-rich regions of the genomes are rapid and are easily visible on a pyrimidine-purine walk graph. These rules are followed, with few exceptions, independent of which strand encodes the gene. Despite the robustness of these rules across a diverse set of species, the magnitude of the differences between the pyrimidine and purine content is fairly small. Typically, the mitochondrial membrane protein transcripts have a pyrimidine richness of 56%, the rRNA transcripts are 55% purine, and the soluble protein transcripts are only 53% purine.

CONCLUSION

The pyrimidine richness of mitochondrial-encoded membrane protein transcripts is partly driven by U nucleotides in the second codon position in all species, which yields hydrophobic amino acids. The purine-richness of soluble protein transcripts is mainly driven by A nucleotides in the first codon position. The purine-richness of rRNA is also due to an abundance of A nucleotides. Possible mechanisms as to how these trends are maintained in mtDNA genomes of such diverse ancestry, size and variability of A-T richness are discussed.

Nucleic acid chaperons: a theory of an RNA-assisted protein folding

Background

Proteins are assumed to contain all the information necessary for unambiguous folding (Anfinsen's principle). However, ab initio structure prediction is often not successful because the amino acid sequence itself is not sufficient to guide between endless folding possibilities. It seems to be a logical to try to find the "missing" information in nucleic acids, in the redundant codon base.

Results

mRNA energy dot plots and protein residue contact maps were found to be rather similar. The structure of mRNA is also conserved if the protein structure is conserved, even if the sequence similarity is low. These observations led me to suppose that some similarity might exist between nucleic acid and protein folding. I found that amino acid pairs, which are co-located in the protein structure, are preferentially coded by complementary codons. This codon complementarity is not perfect; it is suboptimal where the 1st and 3rd codon residues are complementary to each other in reverse orientation, while the 2nd codon letters may be, but are not necessarily, complementary.

Conclusion

Partial complementary coding of co-locating amino acids in protein structures suggests that mRNA assists in protein folding and functions not only as a template but even as a chaperon during translation. This function explains the role of wobble bases and answers the mystery of why we have a redundant codon base.

In search of the nature of specific nucleic acid-protein interactions

The theory of "codon-amino acid coevolution" was first proposed by Woese in 1967. It suggests that there is a stereochemical matching - that is, affinity - between amino acids and certain of the base triplet sequences that code for those amino acids. We have constructed a Common Periodic Table of Codons and Amino Acids, where the Nucleic Acid Table showed perfect axial symmetry for codons and the corresponding Amino Acid Table also displayed periodicity regarding the biochemical properties (charge and hydrophobicity) of the 20 amino acids and the position of the stop signals. The Table indicates that the middle (2nd) amino acid in the codon has a prominent role in determining some of the structural features of the amino acids. The possibility that physical contact between codons and amino acids might exist was tested on restriction enzymes. Many recognition site-like sequences were found in the coding sequences of these enzymes and as many as 73 examples of codon-amino acid co-location were observed in the 7 known 3D structures (December 2003) of endonuclease-nucleic acid complexes. These results indicate that the smallest possible units of specific nucleic acid-protein interaction are indeed the stereochemically compatible codons and amino acids.

Seven fundamental, unsolved questions in molecular biology. Cooperative storage and bi-directional transfer of biological information by nucleic acids and proteins: an alternative to "central dogma"

The Human Genome Mapping Project provided us a large amount of sequence data. However our understanding of these data did not grow proportionally, because old dogmas still set the limits of our thinking. The gene-centric, reductionistical side of molecular biology is reviewed and seven problems are formulated, each indicating the insufficiency of the "central dogma". The following is concluded and suggested: 1. Genes are located and expressed on both DNA strands; 2. Introns are the source of important biological regulation and diversity; 3. Repeats are the frame of the chromatin structure and participate in the chromatin regulation; 4. The molecular accessibility of the canonical dsDNA structure is poor; 5. The genetic code is co-evolved with the amino acids and there is a stereochemical matching between the codes andamino acids; 6. The flow of information between nucleic acids and proteins is bi-directional and reverse translation might exist; 7. Complex genetic information is always carried and stored by nucleic acids and proteins together.

Frequent occurrence of recognition site-like sequences in the restriction endonucleases

Background

There are two different theories about the development of the genetic code. Woese suggested that it was developed in connection with the amino acid repertoire, while Crick argued that any connection between codons and amino acids is only the result of an "accident". This question is fundamental to understand the nature of specific protein-nucleic acid interactions.

Results

The nature of specific protein-nucleic acid interaction between restriction endonucleases (RE) and their recognition sequences (RS) was studied by bioinformatics methods. It was found that the frequency of 5–6 residue long RS-like oligonucleotides is unexpectedly high in the nucleic acid sequence of the corresponding RE (p < 0.05 and p < 0.001 respectively, n = 7). There is an extensive conservation of these RS-like sequences in RE isoschizomers. A review of the seven available crystallographic studies showed that the amino acids coded by codons that are subsets of recognition sequences were often closely located to the RS itself and they were in many cases directly adjacent to the codon-like triplets in the RS.

Fifty-five examples of this codon-amino acid co-localization are found and analyzed, which represents 41.5% of total 132 amino acids which are localized within 8 Å distance to the C1' atoms in the DNA. The average distance between the closest atoms in the codons and amino acids is 5.5 +/- 0.2 Å (mean +/- S.E.M, n = 55), while the distance between the nitrogen and oxygen atoms of the co-localized molecules is significantly shorter, (3.4 +/- 0.2 Å, p < 0.001, n = 15), when positively charged amino acids are involved. This is indicating that an interaction between the nucleic- and amino acids might occur.

Conclusion

We interpret these results in favor of Woese and suggest that the genetic code is "rational" and there is a stereospecific relationship between the codes and the amino acids.