BioSimplex

Representing DNA sequences as tetahedrals (Simplex)

Overview

Representing DNA sequences as regular tetrahedrals (Simplex)

This packages has a single public function biosimplex that takes a BioSequence and returns a Simplex representation of a BioSequence. The Simplex representation is a 3D representation of the BioSequence where each base can be represented as unit vectors pointing into a regular tetrahedron (Silverman et al., 1986; Coward, 1997).

Installation

You can install BioSimplex from the julia REPL. Press ] to enter pkg mode, and enter the following:

pkg> add BioSimplex

Usage

using BioSequences, BioSimplex

# Create a BioSequence
seq = dna"ATCG"

# Convert the BioSequence to a Simplex representation 
biosimplex(seq)

3×4 Matrix{Float64}:
 0.0   0.942809  -0.471405  -0.471405
 0.0   0.0        0.816497  -0.816497
 1.0  -0.333333  -0.333333  -0.333333

Applications

The Simplex representation is useful for to generate a numerical representation of the sequences so that it can be used in machine learning models.

References

Coward, E. (1997). Equivalence of two Fourier methods for biological sequences. Journal of Mathematical Biology, 36(1), 64–70. https://doi.org/10.1007/s002850050090

Silverman, B. D., & Linsker, R. (1986). A measure of DNA periodicity. Journal of Theoretical Biology, 118(3), 295–300. https://doi.org/10.1016/S0022-5193(86)80060-1