The degeneracy of the genetic code and Hadamard matrices. Sergey V. Petoukhov

Similar documents
Genetic Code, Attributive Mappings and Stochastic Matrices

Aoife McLysaght Dept. of Genetics Trinity College Dublin

The genetic code, 8-dimensional hypercomplex numbers and dyadic shifts. Sergey V. Petoukhov

Objective: You will be able to justify the claim that organisms share many conserved core processes and features.

ATTRIBUTIVE CONCEPTION OF GENETIC CODE, ITS BI-PERIODIC TABLES AND PROBLEM OF UNIFICATION BASES OF BIOLOGICAL LANGUAGES *

Genetic code on the dyadic plane

Mathematics of Bioinformatics ---Theory, Practice, and Applications (Part II)

Using an Artificial Regulatory Network to Investigate Neural Computation

A p-adic Model of DNA Sequence and Genetic Code 1

UNIVERSITY OF CAMBRIDGE INTERNATIONAL EXAMINATIONS General Certifi cate of Education Advanced Subsidiary Level and Advanced Level

Edinburgh Research Explorer

Lecture IV A. Shannon s theory of noisy channels and molecular codes

In previous lecture. Shannon s information measure x. Intuitive notion: H = number of required yes/no questions.

Reducing Redundancy of Codons through Total Graph

A Minimum Principle in Codon-Anticodon Interaction

Biology 155 Practice FINAL EXAM

CHEMISTRY 9701/42 Paper 4 Structured Questions May/June hours Candidates answer on the Question Paper. Additional Materials: Data Booklet

A modular Fibonacci sequence in proteins

Abstract Following Petoukhov and his collaborators we use two length n zero-one sequences, α and β,

HADAMARD MATRICES AND QUINT MATRICES IN MATRIX PRESENTATIONS OF MOLECULAR GENETIC SYSTEMS

SEQUENCE ALIGNMENT BACKGROUND: BIOINFORMATICS. Prokaryotes and Eukaryotes. DNA and RNA

A Mathematical Model of the Genetic Code, the Origin of Protein Coding, and the Ribosome as a Dynamical Molecular Machine

Crystal Basis Model of the Genetic Code: Structure and Consequences

Natural Selection. Nothing in Biology makes sense, except in the light of evolution. T. Dobzhansky

Lect. 19. Natural Selection I. 4 April 2017 EEB 2245, C. Simon

Three-Dimensional Algebraic Models of the trna Code and 12 Graphs for Representing the Amino Acids

THE GENETIC CODE INVARIANCE: WHEN EULER AND FIBONACCI MEET

The Genetic Code Degeneracy and the Amino Acids Chemical Composition are Connected

Analysis of Codon Usage Bias of Delta 6 Fatty Acid Elongase Gene in Pyramimonas cordata isolate CS-140

C CH 3 N C COOH. Write the structural formulas of all of the dipeptides that they could form with each other.

Natural Selection. Nothing in Biology makes sense, except in the light of evolution. T. Dobzhansky

The genetic system and algebras of projection operators Sergey Petoukhov Head of Lab of Biomechanical systems, Russian Academy of Sciences, Moscow

Ribosome kinetics and aa-trna competition determine rate and fidelity of peptide synthesis

Foundations of biomaterials: Models of protein solvation

Get started on your Cornell notes right away

CODING A LIFE FULL OF ERRORS

In Silico Modelling and Analysis of Ribosome Kinetics and aa-trna Competition

PROTEIN SYNTHESIS INTRO

SYMMETRIES, GENERALIZED NUMBERS AND HARMONIC LAWS IN MATRIX GENETICS

In Silico Modelling and Analysis of Ribosome Kinetics and aa-trna Competition

Fundamental mathematical structures applied to physics and biology. Peter Rowlands and Vanessa Hill

Lesson Overview. Ribosomes and Protein Synthesis 13.2

Lecture 27. Phylogeny methods, part 4 (Models of DNA and protein change) p.1/26

The Trigram and other Fundamental Philosophies

Molecular Evolution and Phylogenetic Analysis

Interpolated Markov Models for Gene Finding. BMI/CS 776 Spring 2015 Colin Dewey

2013 Japan Student Services Origanization

GENETIC CODE AND THE ANCIENT CHINESE BOOK OF CHANGES

Slide 1 / 54. Gene Expression in Eukaryotic cells

The translation machinery of the cell works with triples of types of RNA bases. Any triple of RNA bases is known as a codon. The set of codons is

Six Fractal Codes of Biological Life: Perspectives in Astrobiology and Emergence of Binary Logics

Practical Bioinformatics

Recent Evidence for Evolution of the Genetic Code

Week 6: Protein sequence models, likelihood, hidden Markov models

ANALOGIES OF GENETIC AND CHEMICAL CODE. Miloje M. Rakočević

Protein Threading. Combinatorial optimization approach. Stefan Balev.

Transcription Attenuation

In this article, we investigate the possible existence of errordetection/correction

An algebraic hypothesis about the primeval genetic code

Modelling and Analysis in Bioinformatics. Lecture 1: Genomic k-mer Statistics

Structure and Function. Overview of BCOR 11. Various forms of Carbon molecules and functional groups

arxiv: v2 [physics.bio-ph] 8 Mar 2018

High throughput near infrared screening discovers DNA-templated silver clusters with peak fluorescence beyond 950 nm

Nature Structural & Molecular Biology: doi: /nsmb Supplementary Figure 1

NSCI Basic Properties of Life and The Biochemistry of Life on Earth

3. Evolution makes sense of homologies. 3. Evolution makes sense of homologies. 3. Evolution makes sense of homologies

Videos. Bozeman, transcription and translation: Crashcourse: Transcription and Translation -

Chemistry Chapter 26

Protein Synthesis. Unit 6 Goal: Students will be able to describe the processes of transcription and translation.

1. Contains the sugar ribose instead of deoxyribose. 2. Single-stranded instead of double stranded. 3. Contains uracil in place of thymine.

Chapter 4: Amino Acids

It is the author's version of the article accepted for publication in the journal "Biosystems" on 03/10/2015.

Lecture 15: Realities of Genome Assembly Protein Sequencing

Expression pattern and, surprisingly, gene length shape codon usage in Caenorhabditis, Drosophila, and Arabidopsis

THE MATHEMATICAL STRUCTURE OF THE GENETIC CODE: A TOOL FOR INQUIRING ON THE ORIGIN OF LIFE

Crick s early Hypothesis Revisited

From Gene to Protein

On Eigenvalues of Row-Inverted Sylvester Hadamard Matrices

DO NOT OPEN THE EXAMINATION PAPER UNTIL YOU ARE TOLD BY THE SUPERVISOR TO BEGIN

Protein Synthesis. Unit 6 Goal: Students will be able to describe the processes of transcription and translation.

NIH Public Access Author Manuscript J Theor Biol. Author manuscript; available in PMC 2009 April 21.

Proteins: Characteristics and Properties of Amino Acids

Advanced topics in bioinformatics

From DNA to protein, i.e. the central dogma

Codon Distribution in Error-Detecting Circular Codes

Genome and language two scripts of heredity

Gene Finding Using Rt-pcr Tests

Journal of Biometrics & Biostatistics

Supplemental data. Pommerrenig et al. (2011). Plant Cell /tpc

Practice Problems 6. a) Why is there such a big difference between the length of the HMG CoA gene found on chromosome 5 and the length of the mrna?

On the optimality of the standard genetic code: the role of stop codons

Answers to Chapter 6 (in-text & asterisked problems)

Protein Struktur. Biologen und Chemiker dürfen mit Handys spielen (leise) go home, go to sleep. wake up at slide 39

Packing of Secondary Structures

A Brief History of Life on Earth

SUPPORTING INFORMATION FOR. SEquence-Enabled Reassembly of β-lactamase (SEER-LAC): a Sensitive Method for the Detection of Double-Stranded DNA

Advanced Topics in RNA and DNA. DNA Microarrays Aptamers

Ranjit P. Bahadur Assistant Professor Department of Biotechnology Indian Institute of Technology Kharagpur, India. 1 st November, 2013

Could Genetic Code Be Understood Number Theoretically?

Characterization of Pathogenic Genes through Condensed Matrix Method, Case Study through Bacterial Zeta Toxin

Transcription:

The degeneracy of the genetic code and Hadamard matrices Sergey V. Petoukhov Department of Biomechanics, Mechanical Engineering Research Institute of the Russian Academy of Sciences petoukhov@hotmail.com, petoukhov@imash.ru, http://symmetry.hu/isabm/petoukhov.html Abstract: The matrix form of the presentation of the genetic code is described as a cognitive form to analyze structures of the genetic code. A similar matrix form is utilized in the theory of signal processing. The Kronecker family of the genetic matrices is investigated, which is based on the genetic matrix [C A; U G], where C, A, U, G are the letters of the genetic alphabet. This matrix in the third Kronecker power is the (8*8)-matrix, which contains 64 triplets. Peculiarities of the degeneracy of the vertebrate mitochondrial genetic code are reflected in the symmetrical black-and-white mosaic of this genetic (8*8)-matrix. This mosaic matrix is connected algorithmically with Hadamard matrices unexpectedly, which are famous in the theory of signal processing, spectral analysis, quantum mechanics and quantum computers. KEYWORDS: genetic code, degeneracy, Hadamard matrix, symmetry, permutation 1 Introduction Genetic information is transferred by means of discrete elements: 4 letters of the genetic alphabet, 64 triplets, 20 amino acids, etc. General theory of signal processing utilizes the encoding of discrete signals by means of special mathematical matrices and spectral representations of signals to increase reliability and efficiency of information transfer [Sklar, 2001; Ahmed, Rao, 1975, etc]. A typical example of such matrices is the family of Hadamard matrices. Rows of Hadamard matrices form an orthogonal system of functions Hadamard-Walsh, which is used for the spectral presentation and transfer of discrete signals [Ahmed, Rao, 1975; Geramita, 1979; Yarlagadda, Hershey, 1997]. An investigation of structural analogies between computer informatics and genetic informatics is one of the important tasks of modern science in a connection with a creation of DNA-computers and with a development of bioinformatics. The author investigates molecular structures of the genetic code from the viewpoint of matrix methods of encoding discrete signals. 2 Method This section describes details of the method of the matrix presentation of genetic multiplets to study symmetries and other structural peculiarities of genetic code systems. The utility of the matrix approach for investigations of the genetic code systems was demonstrated in the works [Petoukhov, 2001-2008; He et al., 2004; He, Petoukhov, 2007]. Hadamard matrices. By a definition a Hadamard matrix of dimension n is the (n*n)-matrix H(n) with elements +1 and -1. It satisfies the condition H(n)*H(n) T = n*i n, where H(n) T is the transposed matrix and I n is the (n*n)-identity matrix. The Hadamard matrices of dimension 2 k are formed, for example, by the recursive formula H(2 k ) = H(2) (k) = H(2) H(2 k-1 ) for 2 k N, where denotes the Kronecker (or tensor) product, (k) means the Kronecker exponentiation, k and N are integers, H(2) is demonstrated in Figure 1. Rows of a Hadamard matrix are mutually orthogonal. It means that every two different rows in a Hadamard matrix represent two perpendicular vectors, a scalar product of which is equal to 0. The element -1 can be disposed in any of four positions in the Hadamard matrix H(2).

1 1 1 1 1 1-1 1-1 1 H(2 k-1 ) H(2 k-1 ) H(2) = -1 1 ; H(4) = -1-1 1 1 ; H(2 k ) = -H(2 k-1 ) H(2 k-1 ) 1-1 -1 1 Figure 1: The family of Hadamard matrices H(2 k ) based on the Kronecker product. Such matrices are used in many fields due to their advantageous properties: in error-correcting codes such as the Reed-Muller code; in spectral analysis and multi-channel spectrometers with Hadamard transformations; in quantum computers with Hadamard gates, etc. The author discovers unexpectedly that Hadamard matrices reflect essential peculiarities of molecular genetic systems. Normalized Hadamard (2x2)-matrices are matrices of rotation on 45 0 or 135 0 depending on an arrangement of signs of its individual elements. A Kronecker product of two Hadamard matrices is a Hadamard matrix as well. A permutation of any columns or rows of a Hadamard matrix leads to a new Hadamard matrix. Hadamard matrices and their Kronecker powers are used widely in spectral methods of analysis and processing of discrete signals and in quantum computers. A transform of a vector ā by means of a Hadamard matrix H gives the vector ū = Н*ā, which is named Hadamard spectrum. A greater analogy between Hadamard transforms and Fourier transforms exists [Ahmed & Rao, 1975]. In particular the fast Hadamard transform exists in parallel with the fast Fourier transform. The whole class of multichannel "spectrometers with Hadamard transforms is known [Tolmachev, 1976], where the principle of tape masks (or chain masks) is used, and it reminds one of the principles of a chain construction of genetic texts in DNA. Hadamard matrices are used widely in the theory of coding (for example, they are connected with Reed Muller error correcting codes and with Hadamard codes [Peterson & Weldon, 1972; Solovieva, 2006], the theory of compression of signals and images, a realization of Boolean functions by means of spectral methods, the theory of planning of multiple-factor experiments and in many other branches of mathematics. Rows of Hadamard matrices are named Walsh functions or Hadamard functions. Walsh functions can be represented in terms of product of Rademacher functions r n (t) = sign(sin2 n πt), n = 1,2,3,, which accept the two values "+1" and "-1" only (here sign is the function of a sign on argument). Sets of numerated Walsh functions (or Hadamard functions), when they are united in square matrices, form systems depending on features of such union. Figure 2 shows two examples of systems of such functions, which are used widely in the theory of digital signals processing. They are connected with (8x8)-matrices by Hadamard and with the Walsh-Hadamard transform, which is the most famous among non-sinusoidal orthogonal transforms and which can be calculated by means of mathematical operations of addition and subtraction only (see more detail in [Ahmed & Rao, 1975; Trahtman & Trahtman, 1975; Yarlagadda, & Hershey, 1997]. Hereinafter we will use the simplified designations of matrix elements on illustrations of Hadamard matrices: the symbol + or the black color of a matrix cell means the element +1 ; the symbol - or the white color of a matrix cell means the element -1. The theory of discrete signals pays special attention to quantities of changes of signs + and - along each row and each column in Hadamard matrices. These quantities are connected with an important notion of sequency as a generalization of notion of frequency [Ahmed & Rao, 1975, p.85]. Figure 2

shows these quantities for each row and each column in presented matrix examples. Each of these two Hadamard matrices is symmetrical relative to its main diagonal. 0 0 7 1 3 3 4 2 1 7 6 6 2 4 5 5 0 7 3 4 1 6 2 5 0 1 3 2 7 6 4 5 Figure 2. Examples of two systems of Walsh functions (or Hadamard functions), which are used frequently in the theory of digital signals processing [Trahtman & Trahtman, 1975]. On the left side: the Walsh-Hadamard system. On the right side: the Walsh-Paley system. Quantities of changes of signs + and - are shown for each row and each column. Normalized Hadamard matrices are unitary operators. They serve as one of the important instruments to create quantum computers, which utilize so called Hadamard gates (as evolution of the closed quantum system is unitary) [Nielsen & Chuang, 2001]. Let us demonstrate now unexpected connections of Hadamard matrices with Kronecker families of genetic matrices. The Kronecker family of genetic matrices. Genetic multiplets are one of main peculiarities of the genetic code. Really, the alphabet of the genetic code is a set of 4 monoplets (nitrogenous bases): A (adenine), C (cytosine), G (guanine), U/T (uracil in RNA or thymine in DNA); 64 triplets encode amino acids; each protein is encoded by more or less long multiplets. On the base of the idea about analogies between computer informatics and genetic informatics, the author has represented all sets of genetic multiplets as proper parts of the Kronecker family of the square matrices (genomatrices) P (n) = [C A; U G] (n), where A, C, G, U are the letters of the genetic alphabet, (n) means the Kronecker exponentiation. Each genomatrix [C A; U G] (n) contains a complete set of n-plets as its matrix elements. For example, the (8x8)-genomatrix [C A; U G] (3) contains all 64 triplets which encode 20 amino acids and stop-signals. This family of matrices gives complete sets of n-plets in the universal mathematical form, which is based on the square matrix of the genetic alphabet. The black-and-white matrix mosaic of the degeneracy of the vertebrate mitochondrial genetic code. Modern science knows many dialects of the genetic code (http://www.ncbi.nlm.nih.gov/taxonomy/utils/wprintgc.cgi). In this article we shall investigate the case of the vertebrate mitochondria genetic code, which is considered in molecular genetics as the most ancient and perfect dialect of the genetic code [Frank-Kamenetskiy, 1988]. Figure 4 demonstrates the disposition of 20 amino acids and stop-signals in the genomatrix [C A; U G] (3) for this code. Each sub-quadrant (2*2) of the matrix [C A; U G] (3) contains a subfamily of those four triplets, which are identical to each other by two first letters. Such quadruple of triplets will be named a subfamily of NN-triplets. For the case of the vertebrate mitochondria genetic code, a complete set of 16 subfamilies of NN-triplets is divided by the nature into two subsets with 8 subfamilies in each. The first subset with 32 triplets, disposed in dark cells on Figure 4, contains those subfamilies of NN-triplets, coded meanings of which are independent of their third letter. In this reason, all four triplets of such subfamily encode the

same amino acids. We name these 32 triplets "black triplets" conditionally. On the contrary, the second subset with other 32 triplets, disposed in white matrix cells, contains those subfamilies of NN-triplets, coded meanings of which are dependent of their third letter. In this reason, each such subfamily has triplets, which encode different amino acids or stop-signals. We name these 32 triplets "white triplets" conditionally. Р = CC CA AC AA С А ; Р (2) = Р Р = CU CG AU AG U G UC UA GC GA UU UG GU GG CCC CCA CAC CAA ACC ACA AAC AAA CCU CCG CAU CAG ACU ACG AAU AAG CUC CUA CGC CGA AUC AUA AGC AGA P (3) = CUU CUG CGU CGG AUU AUG AGU AGG UCC UCA UAC UAA GCC GCA GAC GAA UCU UCG UAU UAG GCU GCG GAU GAG UUC UUA UGC UGA GUC GUA GGC GGA UUU UUG UGU UGG GUU GUG GGU GGG Figure 3: The beginning of the Kronecker family of symbolic genomatrices Р (n) = [C A; U G] (3) for n = 1, 2, 3. The disposition of these black and white triplets forms the very symmetric black-and-white mosaic in the genomatrix [C A; U G] (3) (Figure 4). For instance, left and right matrix halves are mirror-anti-symmetric to each other in its colors: any pair of cells, disposed by the mirrorsymmetrical manner in these halves, has opposite colors. Diagonal quadrants of the matrix are CCC PRO CCU PRO CUC LEU CUU LEU CCA PRO CCG PRO CUA LEU CUG LEU CAC HIS CAU HIS CGC ARG CGU ARG CAA GLN CAG GLN CGA ARG CGG ARG ACC THR ACU THR AUC ILE AUU ILE ACA THR ACG THR AUA MET AUG MET AAC ASN AAU ASN AGC SER AGU SER AAA LYS AAG LYS AGA STOP AGG STOP UCC SER UCU SER UUC PHE UUU PHE UCA SER UCG SER UUA LEU UUG LEU UAC TYR UAU TYR UGC CYS UGU CYS UAA STOP UAG STOP UGA TRP UGG TRP GCC ALA GCU ALA GUC VAL GUU VAL GCA ALA GCG ALA GUA VAL GUG VAL GAC ASP GAU ASP GGC GLY GGU GLY GAA GLU GAG GLU GGA GLY GGG GLY Figure 4: The representation of the genomatrix [C A; U G] (3) of 64 triplets for the case of the vertebrate mitochondrial genetic code. The matrix contains 20 amino acids with their traditional abbreviations. Stop-codons are marked as stop.

identical to each other from the viewpoint of their mosaic. The rows 1-2, 3-4, 5-6, 7-8 are identical to each other from the viewpoint of the mosaic and of the disposition of the same amino acids in their proper cells, etc. A sequence of black and white cells in each row corresponds to one of Rademacher functions. It should be mentioned that the quantity of variants of possible dispositions of 64 genetic triplets in 64 cells of a (8*8)-matrix is equal to the huge number 64!~10 89. It is obvious that the most of these variants have not symmetries in a mosaic disposition of these black and white triplets. Why the nature has chosen (from the huge number of possible variants) this specific code degeneracy, which is reflected in such symmetric mosaic? Below we shall demonstrate that this mosaic of the genomatrix [C A; U G] (3) is connected with the disposition mosaic of the signs +1 and -1 in a certain Hadamard matrix by a simple algorithm based on molecular peculiarities of the genetic alphabet. Perhaps, this fact will lead to the desired answer. For this demonstration we mark black and white cells in the genomatrix on Figure 4 by the signs +1 and -1 correspondingly. In the result we receive the numeric matrix G on Figure 5 (on the left). The genomatrix G is transformed into the Hadamard matrix G H on Figure 5 (on the right) by the simple algorithm described in the next section. +1 +1-1 -1 +1 +1-1 -1 +1 +1-1 -1 +1 +1-1 -1 3 +1 +1-1 -1 +1 +1-1 -1-1 +1 +1-1 -1 +1 +1-1 4 +1 +1 +1 +1-1 -1-1 -1 +1 +1 +1 +1-1 -1-1 -1 1 +1 +1 +1 +1-1 -1-1 -1-1 +1-1 +1 +1-1 +1 +1 6 +1 +1-1 -1 +1 +1-1 -1-1 -1 +1 +1 +1 +1-1 -1 2 +1 +1-1 -1 +1 +1-1 -1 +1-1 -1 +1-1 +1 +1-1 5-1 -1-1 -1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 +1 0-1 -1-1 -1 +1 +1 +1 +1-1 +1-1 +1-1 +1-1 +1 7 5 2 6 1 5 2 6 1 Figure 5: The genomatrix G (on the left) and the Hadamard matrices G H (on the right) Amino-groups NH 2 and the U-algorithm. It is known that amino-groups NH 2 play an important role in molecular genetics. For instance, the amino-group in amino acids provides a function of recognition of the amino acids by ferments [Shapeville, Haenni, 1974]. A detachment of amino-groups in nitrogenous bases A and C in RNA under action of nitrous acid HNO 2 determines a property of amino-mutating of these bases, which was used to divide the set of 64 triplets in eight natural subsets with 8 triplets in each [Wittmann, 1961]. But how the amino-groups are represented in the genetic alphabet? One can note that each of three nitrogenous bases A, C, G has one amino-group, but the fourth basis U/T has not it Figure 6). From the viewpoint of existence of the amino-group, the letters A, C, G are identical to each other and the letter U is opposite to them.

Figure 6. The complementary pairs of the four nitrogenous bases in DNA. A-T (adenine and thymine), C-G (cytosine and guanine). Hydrogen bonds in these pairs are shown by dotted lines. Black circles are atoms of carbon; small white circles are atoms of hydrogen; squares with the letter N are atoms of nitrogen; triangles with the letter O are atoms of oxygen. Amides (or amino-groups) NH 2 are marked by big circles. This fact of existence or absence of the amino-group in certain genetic letters can be reflected in the alphabetic genomatrix P=[C A; U G] by symbols +1 and -1 instead of the letters A, C, G and U/T correspondingly. In this case the Hadamard genomatrix P H(2) = H(2) is appeared (Figure 7). C A +1 +1 P = U G H(2) = -1 +1 Figure 7: The transformation of the matrix P=[C A; U G] of the genetic alphabet into the Hadamard matrix. Black cells of this Hadamard matrix contain elements +1, and the white cell contains the element -1. Hence the letter U in RNA (and the letter T in DNA) is the peculiar letter in the genetic alphabet in this sense. The letter U has also another unique property among 4 letters of the genetic alphabet: the letter U exists in RNA, but in DNA it is replaced by the letter T. These molecular properties of the letter U can be utilized in genetic computers of organisms. Taking into account this unique status of the letter U, the author revealed the existence of the following formal U-algorithm, which demonstrates the close connection between Hadamard matrices and the matrix mosaic of the degeneracy of the genetic code. The definition of the U-algorithm: each triplet in the black-and-white genomatrix [C A; U G] (3) on Figure 4 should change its own color into opposite color each time when the letter U stands in

an odd position (in the first or in the third position) inside the triplet. For example, the white triplet UUA (see Figure 4) should become the black triplet (and its matrix cell should be marked by black color) because of the letter U in its first position. Or the white triplet UUU should not change its color because of the letter U in its first and third positions (the color of this triplet is changed twice according to the described algorithm). The triplet ACG does not change its color because the letter U is absent in this triplet at all. By means of the U-algorithm, the genomatrix P (3) = [C A; U G] (3) from Figure 4 is transformed into the genomatrix P H (3) on Figure 7. CCC CCA CAC CAA ACC ACA AAC AAA CCU CCG CAU CAG ACU ACG AAU AAG CUC CUA CGC CGA AUC AUA AGC AGA CUU CUG CGU CGG AUU AUG AGU AGG UCC UCA UAC UAA GCC GCA GAC GAA UCU UCG UAU UAG GCU GCG GAU GAG UUC UUA UGC UGA GUC GUA GGC GGA UUU UUG UGU UGG GUU GUG GGU GGG Figure 8: The mosaic genomatrix P H (3), which is received from [C A; U G] (3) by the U-algorithm. 3 Results The genomatrix P H (3) (Figure 8) possesses the black-and-white mosaic which is identical to the mosaic of the matrix G H on Figure 5 (on the right). If we mark each black (white) cell of the P H (3) by the number +1 ( -1 ), the matrix P H (3) is transformed into the matrix G H. But the matrix G H is one of Hadamard matrices because it satisfies the general condition of Hadamard matrices H(n)*H(n) T = n*i n. The Hadamard matrix G H has the interesting property of a fractal character. Each (4*4)-quadrant of the (8*8)-matrix G H is a Hadamard matrix also. Furthermore, each (2*2)-sub-quadrant of the (8*8)-matrix G H is a Hadamard matrix as well (this situation can be named Hadamard fractals conditionally). One can mention additionally that each black-and-white row in the matrix G (Figure 4) corresponds to a proper Rademacker function, which is connected with Hadamard- Walsh functions of Hadamard matrices [Ahmed, Rao, 1975]. CCC Pro CUC Leu UCC Ser UUC Phe CCU Pro CUU Leu UCU Ser UUU Phe CAC His ACC Thr AAC Asn CCA Pro CAA Gln ACA Thr CGC AUC AGC CUA CGA AUA Arg Ile Ser Leu Arg Met UAC GCC GAC UCA UAA GCA Tyr Ala Asp Ser Stop Ala UGC GUC GGC UUA UGA GUA Cys Val Gly Leu Trp Val CAU ACU AAU CCG CAG ACG His Thr Asn Pro Gln Thr CGU AUU AGU CUG CGG AUG Arg Ile Ser Leu Arg Met UAU GCU GAU UCG UAG GCG Tyr Ala Asp Ser Stop Ala UGU GUU GGU UUG UGG GUG Cys Val Gly Leu Trp Val Figure 9: The genomatrix P (3) 231. AAA Lys AGA Stop GAA Glu GGA Gly AAG Lys AGG Stop GAG Glu GGG Gly

The theory of signal processing pays a special attention to permutations of components. Let us analyze properties of the genomatrix P (3) = [C A; U G] (3) relative to positional permutations in all triplets. Any triplet has the six permutation variants of a sequence of its three positions: 1-2-3, 2-3-1, 3-1-2, 1-3-2, 2-1-3, 3-2-1. Let us take for example the permutation of the cyclic shift: 1-2-3 2-3-1. In the result of this permutation the triplet CAG is replaced in its matrix cell by the triplet AGC, etc. And the whole genomatrix P (3) 123 = [C A; U G] (3) is reconstructed cardinally (3) into the new mosaic matrix P 231 (Figure 9), which has considerable symmetries also (all quadrants have the identical mosaics; upper and bottom halves have identical contents, etc). The other four variants of the positional permutations in the triplets produce the four genomatrices P 312 (3), P 132 (3), P 213 (3), P 321 (3), each of which has considerable symmetries also (the bottom indexes show the positional permutations). It means that the degeneracy of the genetic code is connected with the positional permutations in the set of 64 triplets. The same U-algorithm transforms the mosaics of the matrices P 231 (3), P 312 (3), P 132 (3), P 213 (3), P 321 (3) into new mosaics, which coincide with the mosaics of the proper Hadamard matrices on Figure 10. One can note additionally a special feature connected with quantities of changes of the signs + and - in rows and columns of all genetic Hadamard matrices on Figures 5 and 10. Each sum of such quantities in the first four rows and in second four rows of such matrix is equal to 14. And each sum of such quantities in the first four columns and in second four columns of such matrix is equal to 14. It can be named a phenomenon of a balanced sequency in these genetic Hadamard matrices (the notion of sequency is a generalization of notion of frequency [Ahmed & Rao, 1975, p.85]). It should be noted that Hadamard matrices, which usually are used in the theory of signal processing (see Figure 2), do not possess such feature. In addition all genetic Hadamard matrices on Figures 5 and 10 are asymmetric relative to both diagonals in contrast to the Hadamard matrices on Figure 2. Besides the described cyclic permutations of positions in the triplets (123 2 231 312 and 321 213 132), a few other types of cyclic permutations of the genetic elements exist, which reveal new genetic Hadamard matrices. The speech is about the fact that each of such permutations transform the initial genetic matrix [C A; U G] (3) into a new genetic mosaic matrix which can be transformed into a relevant Hadamard matrix by the same U-algorithm unexpectedly. For example, one of such types of permutations is a cyclic alphabetical permutation of the genetic letters C G U A C. The application of this alphabetical permutation to the initial genomatrix P CAUG 123 = [C A; U G] (3) transforms it into a new genomatrix P GCAU 123 = [G C; A U] (3). The permutations of positions in the triplets in this genomatrix P GCAU 123 leads to new genomatrices P GCAU 231, P GCAU 312, P GCAU 321, P GCAU 213, P GCAU 132 by analogy with the case considered above (see Figure 10). The same U-algorithm transforms all these genomatrices into genetic Hadamard matrices presented on Figure 11. The second part of Figure 11 demonstrates additionally a case of Hadamard matrices which are received by means of the U-algorithm from genomatrices P CAGU 123 = [C A; G U] (3) and P CAGU 231, P CAGU 312, P CAGU 321, P CAGU 213, P CAGU 132. All kinds of genetic Hadamard matrices, which were met by the author in this study, possess the following analogical features with the Hadamard matrices on Figures 5 and 10. Firstly, each of such genetic Hadamard matrices is asymmetrical relative to both diagonals. Secondly, each of them demonstrates the described phenomenon of a balanced sequency as well. More details on this theme can be found in the book which is in print now: S.Petoukhov & M. He Symmetrical Analysis Techniques for Genetic Systems and Bioinformatics: Advanced Patterns and Applications, Hershey, IGI Global.

FOR Р CAUG 123: FOR Р CAUG 231: 3 7 4 3 1 4 6 0 2 6 5 2 0 5 7 1 5 2 6 1 5 2 6 1 5 2 5 2 4 3 4 3 FOR Р CAUG 312: FOR Р CAUG 132: 1 7 6 1 2 4 5 2 7 6 0 0 4 5 3 3 3 3 4 4 4 4 3 3 3 4 4 3 3 4 4 3 FOR Р CAUG 213: FOR Р CAUG 321: 1 3 6 4 2 7 5 0 3 2 4 5 0 6 7 1 5 2 5 2 6 1 6 1 5 5 2 2 4 4 3 3 Figure 10. The six Hadamard matrices, which are produced from the six mentioned genomatrices by means of the U-algorithm. The black cells correspond to elements +1 and the white cells correspond to elements -1. Quantities of changes of signs + and - (or changes of colors) are shown for each row and each column. 4 Discussion The described phenomenological facts reveal a close connection of the genetic code (including, firstly, the phenomenon of the degeneracy of the vertebrate mitochondrial genetic code and, secondly, described permutation properties of the genetic code) with a special set of Hadamard matrices. This connection between Hadamard matrices and the genetic code is interesting in a few aspects. Genetic molecules are objects of quantum mechanics, where normalized Hadamard matrices play important role as unitary operators (it is known that an evolution of closed quantum system is described by unitary transformation). In particular, quantum computers use these matrices as Hadamard gates [Nielsen, Chuang, 2001]. In view of this, new theoretical possibilities are appeared to transfer achievements of quantum computer conceptions into the field of molecular genetics and to consider genetic system as a quantum computer. From the viewpoint of quantum mechanics and its unitary operators, first of all, Hadamard operators, a possible answer on the fundamental question about reasons for the nature to choose the fourletters genetic alphabet is the following one: the reason is that simplest unitary matrices in

FOR P GCAU 123 : FOR Р GCAU 231 : FOR Р GCAU 312 : FOR Р GCAU 132 : FOR Р GCAU 213 : FOR Р GCAU 321 : FOR Р CAGU 123 : FOR Р CAGU 231 : FOR Р CAGU 312 : FOR Р CAGU 312 : FOR Р CAGU 213 : FOR Р CAGU 321 : Figure 11. The 12 balanced Hadamard matrices, which are produced from the indicated 12 genomatrices of triplets by means of the U-algorithm. Black cells correspond to elements +1 and white cells correspond to elements -1. two-dimensional space, first of all, Hadamard matrices (and also Pauli matrices, etc.) consist of four elements exactly. It seems very probably that principles of quantum mechanics and quantum computers underlie structural peculiarities of the genetic code. One can suppose that Hadamard genomatrices can be used in genetic systems by analogy with applications of Hadamard matrices in different fields of science and technology: signal processing, error-correcting and other codes, spectral analysis, multi-channel spectrometers, etc. Rows of Hadamard matrices represent orthogonal systems of Hadamard-Walsh functions. Such orthogonal system can be a natural base to organize storage and transfer of genetic information with noise immunity properties by means of the decomposition of genetic sequences on the base of these orthogonal systems and by means of using proper codes (orthogonal, bi-orthogonal,

etc.). In particular, advantages of Hadamard matrices can be exploited in genetic system for spectral analysis of genetic sequences and for utilizing emission spectrums of genetic elements. It is essential that orthogonal systems of Hadamard-Walsh functions were revealed by a few authors in macro-physiological systems, which should be agreed structurally with genetic system for their transferring along a chain of generations [Shiozaki, 1980; Carl, 1974; Ginsburg et all, 1974]. Acknowledgments: Described researches were made by the author in the frame of a long-term cooperation between Russian and Hungarian Academies of Sciences and in the frames of programs of International Society of Symmetry in Bioinformatics (USA, http://polaris.nova.edu/mst/issb) and of International Symmetry Association (Hungary, http://symmetry.hu/). The author is grateful to Frolov K.V., Darvas G., Ne eman Y., He M., Smolianinov V.V., Vladimirov Y.S. for their support. REFERENCES Ahmed N.U., Rao K.R. (1975). Orthogonal transforms for digital signal processing, Springer- Verlag New York, Inc. Carl J.V. (1974). On the use of Walsh functions in mane-made and biological pattern recognition systems. Applications of Walsh function and sequence theory: Proc. 1974 Symposium, Washington, 9-25. Geramita, A. V. (1979). Orthogonal designs: quadratic forms and Hadamard matrices. Dekker. Ginsburg A.P., Carl J.V., Kabrisky M., Hall C.F., Gill R.A. (1974). Physiological aspects of a model for the classification of visual images. Adv. Cybern. and Syst., London e.a., 3, 1289-1305. Frank-Kamenetskiy M.D. (1988) The most principal molecule. Moscow, Nauka. He M., Petoukhov S.V. (2007) Harmony of living nature, symmetries of genetic systems and matrix genetics. International journal of integrative medicine, 1, 1, 41-43. He M., Petoukhov S.V., Ricci P. (2004) Genetic code, hamming distance and stochastic matrices. Bulletin for Mathematical Biology, 66, 1405-1421. Nielsen M.A., Chuang I.L. (2001) Quantum computation and quantum information, Cambridge University Press. Peterson, W.W., & Weldon, E.J. (1972). Error-correcting codes. Cambridge: MIT Press. Petoukhov S.V. (2001) Genetic codes I: binary sub-alphabets, bi-symmetric matrices and golden section; Genetic codes II: numeric rules of degeneracy and the chronocyclic theory. Symmetry: Culture and Science, 12, 3-4, 255-306. Petoukhov S.V. (2003-2004) Attributive conception of genetic code, its bi-periodic tables and problem of unification bases of biological languages. Symmetry: Culture and Science, 14-15, part 1, 281-307. Petoukhov S.V. (2005) The rules of degeneracy and segregations in genetic codes. The chronocyclic conception and parallels with Mendel s laws. - In: Advances in Bioinformatics and its Applications, Series in Mathematical Biology and Medicine, 8, 512-532, World Scientific. Petoukhov S.V. (2008a) Matrix genetics, algebras of the genetic code, noise-immunity. Moscow: Regular and Chaotic Dynamics, 316 p. (in Russian; summary in English is on the http://www.geocities.com/symmetrion/matrix_genetics/matrix_genetics.html) Petoukhov S.V. (2008b) Matrix genetics, part 1: permutations of positions in triplets and symmetries of genetic matrices. 1-12. Retrieved March 06, 2008, from http://arxiv:0803.0888.

Petoukhov, S.V. (2008c). Matrix genetics, part 2: the degeneracy of the genetic code and the octave algebra with two quasi-real units (the Yin-Yang octave algebra ). 1-23. Retrieved March 23, 2008, from http://arxiv:0803.3330. Petoukhov, S.V. (2008d). Matrix genetics, part 3: the evolution of the genetic code from the viewpoint of the genetic octave Yin-Yang-algebra. 1-22. Retrieved May 30, 2008, from http:// arxiv:0805.4692. Petoukhov, S.V. (2008e). Matrix genetics, part 4: cyclic changes of the genetic 8-dimensional Yin-Yang-algebras and the algebraic models of physiological cycles. 1-22.Retrieved September 17, 2008, from http://arxiv:0809.2714. Shapeville F., Haenni A.-L. (1974) Biosynthese des proteins. Hermann. Shiozaki A. (1980) A model of distributed type associated memory with quantized Hadamard transform., Biol.Cybern., 38, 1, 19-22. Sklar B. (2001) Digital communication. Fundamentals and applications. Prentice Hall. Solovieva, F. I. (2006). Inroduction to the theory of coding. Novosibirsk: Novosibirsk University (in Russian). Tolmachev, Yu. A. (1976). New optic spectrometers. Leningrad: Leningrad University (in Russian). Trahtman, A. M., Trahtman, V. A. (1975). The foundations of the theory of discrete signals on finite intervals. Moscow: Sovetskoie Radio (in Russian). Wittmann, H.G. (1961) Ansatze zur Entschlusselung des genetischen Codes. Die Naturwissenschaften, 48, 24, 55. Yarlagadda R., Hershey J. (1997) Hadamard matrix analysis and synthesis with applications to communications and signal/image processing. Kluwer Academic Publ.