Download Presentation

Loading in 3 Seconds

This presentation is the property of its rightful owner.

X

Sponsored Links

- 67 Views
- Uploaded on
- Presentation posted in: General

C. Karanikas, N. Atreas, P. Polychronidou, A. Bakalakos Department of Informatics

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

Problems for Effective Analysis of Biological Data:Discrete Transforms on Symbolic Sequences for String-Matching, Pattern-Recognition and Grammar Detection

C. Karanikas, N. Atreas,

P. Polychronidou, A. Bakalakos

Department of Informatics

Aristotle University of Thessaloniki

- We draw our inspiration from the basic operations that nature performs in biological sequences (Replication, Dilation, Translation, Splicing, etc).
- We introduce a variety of new discrete (linear or non linear) invertible transforms on symbolic sequences:
- The Cyclic Class Transform
- The Stern Brocot Transform
- The Generalisation of Haar Transform
- The Haar-Riesz Product

- Our main target with these transforms is
- to encode-decode localinformation on strings
- to make fast non-exact string matching and pattern recognition.
- to identify some of the grammatical rules of the string-collections.

- Thus we deal with the notion of similarity and distances (such as the edit distance), i.e. distances measuring the number of the operations delete, insert and substitute required to identify two strings.

- Global Security depends on Information Security
- Mobility, Heterogeneity, Size, Complexity of Information Systems make it difficult to detect, prevent, respond to, overcome security threats.
- We deal with a fundamental problem: Detect, recognize, interpret and ultimately assign a meaning to symbolic information such as:
- Biological/Biometric data
- Intelligence Information
- Any other form of information

- We formulate the problem in mathematical/information theoretic terms: Given a pattern and a string find parts of the string similar to the pattern
- Or, find strings that are identical.
- Or, find strings that are similar.
- Also, make the algorithms work in a “noisy” environment (I.e. “forgive” a few accidental errors on a big string)
- Also, find the underlying grammar of the string.
****We need new mathematical tools for effective analysis of data.

- CGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGATCGAT

- Each protein can be written in 3D as a string in an alphabet of 20 letters.

- Let {x1 ,…,xn} be a symbolic sequence in an alphabet {ε1,…, εk} (w.l.g. consider epsilons primes)
- Let p1,…,pn,… , be an increasing sequence of primes.
- The map {x1 ,…,xn} Σxi (pi/pi+1) provides an invertible coding for the collection of all symbolic sequences (as above)
- Similar maps: {x1 ,…,xn} Σi (1-xi/pi) or
- {x1 ,…,xn} Πi (1-xi/pi)

Given a collection of strings (genes/proteins written in an alphabet of four/twenty letters), one can assign a number to pairs of this collection, which tell one how distant or how similar or how close they are.

The Edit (Levenshtein) distance counts the minimum number of the operations: Delete, Insert and Substitute to make to strings identical. E.g. the Edit distance of {1,0,1,1,0,1} and {1,1,1,0,1,0}is 2. (delete the second digit, insert 0 in 5 th position)

- The problem given two strings : ACCGHAAGGHCC , ACGHHAGGHACC
- Replace A 0, C1,H2 and G 3 and consider the corresponding numbers.
- Find a new (fast) transform on symbolic sequences and a measure on them, which is “invariant” under a small number of changes as insert, delete and replace. Next we provide a new distance suitable for biodata (biological or biometric)

- Consider the matrices T(0) and T(1) respectively
- Each string {ε1,…,εn}, where εi = 0 or 1, corresponds to the (dot) product of matrices:
- T(ε1).T(ε2)…T(εn ). For example the string {1,0,1,0,1}T(1).T(0).T(1).T(0).T(1) =
- The pair {8,13} (the sum of each row) is unique and so {1,0,1,0,1} {8,13}

- Now T(0),T(1) and T(2) are
- So we have
- {1,2,1,0,2}
- = {5, 31 ,17}
- By a simple algorithm we get {5, 31 ,17}
- {1,2,1,0,2}

- The second string below, has two differences with the first one, the triples of corresponding numbers are similar. We used successfully this idea to find similar parts of genes.
- {B,A,A,B,A,A,B,G,A,B,G,G,G,B,B}
- {B,A,B,B,A,A,G,A,B,G,G,G,B,B}
- {0.740721, 0.0476564, 0.211623}
- {0.735667,0.0487892,0.215544}

- Definition 1: Let Mn,m the space of all n x m matrices, we define the p-adic dilation operation Dp: Mn,m --> Mn,pm , p=2, 3, …. such that:
- Dp(M)={ Mi,[j/p], i=,1,…,n, j = 1, …,m p}
- where [x] is the largest integer greater than or equal to the number x.

- A mathematical transform mimicking antigen processing must code local information (as. the Haar wavelet transform) cutting data in pieces (peptides). Note that a protein (symbolic sequence in an alphabet of 20 letters) is a union of peptides.
- Introduce transforms/ computational tools for string marching, for pattern recognition of grammars, for detecting dynamical systems as hidden Markov process and Riesz products.
- Results and experience from reserch projects:
- European Project: IST-2000-26016, Immunocomputing
- GSRT 2005-6: Mathematics for Bioinformatical applications – Multiresolution methods for the study of biodata.
- Greek-Bulgarian project (2006-7), on Application of Wavelet Theory in Bioinformatics

- Definition 2 Let 0m= {0,0,…,0} is in M1,mThe p-adic vector translations Tp: M1,m -> M1,pm, p=2, 3, …. such that: Tp(v)= Join[0km,v, 0jm], such that k+j+1=p, where Join means splicing vectors together.

- The matrices above called SPMM, are iteratively generated by dilation and translation of block sub-matrices ( work by N. Atreas, C.K. and P. Polychronidou).
- The determinant of SPMM is ± 1
- The Inverse of SPMM is iteratively generated too.

- Theorem Let rj ,j = 1,…, pn be the j row of the pn x pn Haar matrix, and R(p,n) = ∏ ( 1+ x(k) rk ) its Haar Riesz product with coefficients {x(k), k = 1,…pn} , If t={t(k), k = 1,…pn} is any non-negative data the system: R(p,n) . rj = t . rj has a unique solution w.r.t. {x(k), k = 1,…pn} .

- Input 729 (3^6) samples of a Cantor collection. Output the Riesz Haar co-efficients. Observe that the collection is orthogonal on certain rows of the Haar matrix.