A Supersymmetry and Quantum Cryptosystem with Path Integral Approach in Biology

: The source of cancerous mutations and the relationship to telomeres is explained in an alternative way. We deﬁne the smallest subunit in the genetic code as a loop braid group element. The loop braid group is suitable to be deﬁned as a conﬁguration space in the process of converting the information written in the DNA into the structure of a folded protein. This smallest subunit, or a ﬂying ring in our deﬁnition, is a representation of 8-spinor ﬁeld in the supermanifold of the genetic code. The image of spectral analysis from the tensor correlation of mutation genes as our biological system is produced. We apply the loop braid group for biology and authentication in quantum cryptography to understand the cell cocycle and division mechanism of telomerase aging. A quantum biological cryptosystem is used to detect cancer signatures in 36 genotypes of the bone ALX1 cancer gene. The loop braid group with the RSA algorithm is applied for the calculation of public and private keys as cancer signatures in genes. The key role of this approach is the use of the Chern–Simons current and then the ﬁber bundle representation of the genetic code that allows a quantization procedure.


Introduction
For a long time, scientists have been trying to understand the source of self and nonself authentication keys [1] for cancer cells and the immune system's recognition cells. It might be related to information encryption [2] inside 64 codons and a (CG) n repeated sequence in the non-protein-coding DNA as key pairs [3] of a biological cipher in the cancer quantum public-key cryptosystem-some information from the genus of an ancient organism. This information is a biological Chern-Simons current [4,5] induced from an adaptive behavior field in the genetic code. Some questions arise from a new theory of cell biology. What is the smallest subunit used for a representations of the structures of the folded proteins in a living organism? Is it possible to quantize such a structure by adopting standard procedures of theoretical physics? We try to show that it can be done by a new way of unifying triplet states of DNA, RNA, and protein. The new approach of the Chern-Simons current in the genotype with a mathematical model of the hidden state in a genome [6] is trying to answer some questions. The suppersymetry approach [7] was also used for the model of graphene wormhole and for the computation of the Chern-Simons current in a Josephson junction of superconductor states in the graphene. Additionally, the applications of neural networks for density estimation [8] outperformed parametric and non-parametric approaches.
The paper is organized as follows. In Section 2 we define the concept of using a loop braid group for the genetic code. We discuss the source of translating the DNA sequence into a folded protein by using loop braid group operation. In Section 3 we define new asymmetric cryptography for biology, as a source of the public and private keys in the system of protein-protein interactions. We apply the moduli state space model to RSA algorithm to calculate public and private keys in asymmetric cryptography for biology for samples of 36 genes as a new signature. We use tensor correlation for gene expression to generate the image of a geneonic spectral sequence for a mutation gene signature in Section 4 and an empirical analysis with some concrete cases is developed. For brevity and easier understanding of the issue, we are using the term cancer gene for all genetic mutations and the term cancer signature defines all patterns for known types of damage and abnormalities. In the last section we conclude and discuss all results, giving future perspectives of this approach. Some technical aspects of our model are discussed in Appendixes A and B. Amino acid sequences are given in Appendixes C and D; the cipher text for the encryption of a gene is reported in Appendix E.

Algebraic Construction of G-Quadruplex in Telomere
In each cell division, the telomeres (see Figure 1a,b) are shortened [25] and the total length of DNA becomes shorter. As the result of the shorter biological clock from cell division, living things die. In order to understand the cell cocycle and the division mechanism of telomerase aging, we are looking on the source of cancer as the source of age acceleration and we relate it to the telomere shortening mechanism [26,27] (also see Figure 2 for detail of the biological clock), which is a source of braid group operation [28]-so-called self-diffeomorphism-in the genetic code. In the next part we are also trying to explain the real source and state of cancer as damage to the biological clock which can induce age acceleration [29]. The age acceleration is a relative measurement between the chronological clock and the biological clock in the telomere [30]. Up to now, scientists understood telomeres and telomerase as remnants of ancient viruses that rely on DNA in the chromosomes of living organisms. The telomere is composed of the repeated sequence of (TTAGGG) dt * where 1000 ≤ dt * ≤ 2000. The size of the duplicate sequence at the end of this open chromosome is amplified by six braids caused by six superspaces in time series data of organisms. It is naturally selected by the evolution of transgenic reproduction in all possible directions of the transcription process. The G alphabet might be suitable for the hidden time scale in the biological clock. Definition 1. Let d ∈ D be a superspace of DNA, r ∈ R be a superspace of RNA, and p ∈ P be a superspace of protein. We have a Wigner ray [31] of biological time series data of living organism x t ∈ X t defined by p : D → R, (d, r) → p(d, r) W X t d, r := λd, λr = λ d, r = d, r . (1) The protein structure can be considered a transition state. It is a partition function with d ∈ [d], an equivalent class of homotopic paths of gene expression in DNA, [d] ∈ [D, S 2 ]. Definition 2. Let t = (t 1 , t 2 , t 3 , t 4 ) be a biological clock time scale in four corners of G-quadruplex. We denote t G = t 1 , t G * = t 2 , t NG = t 3 , t NG * = t 4 to be a possible pair of times in four corners of G-quadraplex.
There are two types of cell cocycles (Figure 1c) in the loop structure of G-quadruplex of telomere and in telomerase enzymes; we call them β-and α-cocycles. They are associated with T-loop and D-loop in G-quadruplex ( Figure 2). We visualize the gene expression as an anyon transition state with behavior field β i . It is the cocycles of the flying ring in geneon wave function Ψ(β 1 , β 2 , · · · , β n ) = e i ∑ n i=1 β i = e iβ 1 e iβ 2 · · · e iβ n (see Figure 3). The sequence of a gene in a protein is separable and can be processed by permutation as the mutation in the new order of the gene with transposon Ψ transposon (β i ) and retrotransposon states Ψ retrotransposon (β i ). The dark line in the middle of the plot is DNA with the transcription process. We define a DNA flying ring as an orbit of geneonic transition in superspace of DNA d ∈ D := S 2 , for the RNA flying ring is r ∈ R := S 2 . The trajectory of intersection between these two rings induces a loop braid group in genetic codes in (3 + 3) extra dimensions. The fiber bundle along the trajectory is a folding structure of the secondary protein in a superspace p ∈ P. The upper layer is a starting point for flying; the lower layer is the passive layer of protein folding in the feedback loop. The yRNA and intron are just a result of the Reidemeister move for loop braid group in three trajectories of (d, r, p) (Figure 4b). This interaction of D-brane in a living organism is a coupling in loop space in biological time series data with parasitism state between a host cell of the living organism. One state is used for shortening the telomere, and the other reversed transcription state is used for an adenovirus to prolong the length of telomere during cell division. It is a real function of ncRNA embedded in telomerase. They both have G-quadruplex stabilizers. We can represent the T-loop and D-loop in the end-point of the chromosome by using the loop braid group for biology. The loop braid group for biology is not just a representation for DNA, RNA, and protein transitions ( Figure 2). It can be used also to define the intrinsic source of the energy-momentum tensor with undivided properties in genotypes of underlying behavior fields as a flying ring (Figure 1d) in the genetic code.  . We define gene expression as a cocycle with behavior field β i in geneonic wave function (a). It is equivalent tothe exon state in mRNA. In the right picture (b) is a loop braid group diagram for extended central dogma. There exist six superspaces in Kolmogorov space of biological time series data. The trajectory of d ∈ D means DNA, r ∈ R means RNA, and p ∈ P means a protein. The upper layer is an active layer of gene expression with a flying ring along the trajectory as a loop braid element. The Reidemeister move switches state of gene expression between two loop braid group elements in genetic code.
(a) (b) Figure 4. (a) The picture shows a loop braid group structure in a telomere. The passive layer of a telomerase protein at the end of the chromosome is an element of superspace p * ∈ P * . Inside the layer of telomerase there exists a flying ring of RNA r ∈ R and DNA d ∈ D. When the cell divides the loop braid group element, the genetic code will fly to the active side of a histone protein p ∈ P. The four circular rings inside the loop braid group are T-loop and D-loop structures in a telomere. (b) The picture shows gene expression from DNA to mRNA with a stem-loop structure. It is an example of a flying ring in loop braid group of a superspace of DNA and RNA.

Anyon in Biology and Configuration Phase Space for Protein Folding
The permutation of an alphabet code from an evolutional process induces a mutation and it appears as an element of a symmetric group Ω BIO n . → Ω BIO n . The permutation of alphabet code in mutation induced a symmetric group as group operation. Let G × X t → X t be a loop braid group operation over Kolmogorov space in time series data X t , as configuration space in the process of converting the DNA sequence into a folded protein Definition 4. Let d = (d 1 , d 2 , · · · , d n ) ∈ D and r = (r 1 , r 2 , r 3 , r 4 , · · · , r n ) ∈ R. We have an anyon in biology defined by the wave function of gene expression Ψ d i := Ψ i (d 1 , d 2 , · · · , d n ) ∈ X t with behavior field (α i , β i ) as the cocycle in adaptive behavior in the immune system.
It is used to define the configuration space in the wave function of the transition state for gene expression by Ψ i (d 1 , d 2 , · · · , d n ) ∈ X t , where anyon fields are induced from the behavior of immune system. In DNA there is a pair of genotypes in the form of a pair of alleles (α i , β i ); it is an adaptive behavior of cocycle of gene expression for undocking behavior in the protein-protein docking system.

Circular Artin Braid Group Representation for Spinor Field in Genetic Code
In this section, we assume that all genetic code cannot be completely separated. The biological clock in telomere length is parametrized by a hidden state of the number of dt * := [G] alphabet in a repeated pattern of the telomere (TTAGGG) n=dt * .
The element is Grothendieck topology over an adjoint cofunctor. It is a self-diffeomorphism The loop braid generator for B c 3 is a quaternionic field in genetic code. We define their explicit forms and their permutations over the symmetric group by a chosen basis in Clifford algebra.

Definition 5. Let
therefore, one can write eight bases for spinor field in the genetic code in braid form as follows.
One can also use θ = 2πs with spin quantum number s to be an integer number for the retrotransposon, so that e iθ = e 2iπs = (−1) 2s | ψ 2 ψ 1 .
The red color line represents the curvature from the physiology of biological time series data and blue color represents the active and passive behavior field layers. (b) The picture shows a member of loop braid group. We have three circles S 1 . They are a source of closed 3-balls B c 3 , the structure group of the affine transform of 3 behavioral fields in genetic code σ [A µ ] i ∈ B 2 action on an affine fiber bundle of the behavioral field in the genetic code σ ω ± , σ im ± , σ m ± ∈ B c 3 . The affine group is as a loop braid group in genetic code by B 3 2 B c 3 .

Classification of the Loop Braid Group in Genetic Code
One can classify three types of loop braid group operations as representations of an anyon for protein folding. For a two dimensional representation of D-brane in the loop braid group for the genetic code, one can define the abelian anyon for biology in (2 + 1) dimensions. The extra dimensions are used to represent the homotopy path of protein folding. Definition 6. The loop braid group, LB BIO n , for a genetic code has three types of generators, σ = (σ D , σ R , σ P ), ρ = (σ [A] , σ [U] , σ [C] , σ [G] ), and τ = (σ ω , σ im , σ m ). The generators σ i , ρ i ∈ B 3 2 | i = 1, · · · , n − 1 and τ i ∈ B c 3 | i = 1, · · · , n fulfill following relations Definition 7. Let G be a group operation of the genotype. Let ρ be a representation for gene translation as an anyon. We have ρ : G → U(1) = S 2 as a representation in Ω BIO n and in LB BIO n . In order to visualize 3D folding structure of the protein, we define three types of loop braid group operations in biology. All loop elements of the representation of amino acids arose from group operations over the superspace of time series data. These tree types are a translation, reflection, and rotation; ρ T−I , ρ T−I I , and ρ T−I I I . For the translation as a string of amino acids, we have the anyon T-I for biology: Considering an analogy with the action of the symmetric group by permutations, one can find a natural action of the braid group on n-tuples of objects or on the n-folded tensor product that involves some twistors. For this purpose let us consider an arbitrary group G and let X be the set of all n-tuples of elements of G whose product is the identity element of G. The kernel of the homomorphism LB BIO n → Ω BIO n is a subgroup of LB BIO n called pure loop braid group for biology on n strands and it is denoted as LP BIO n . In the pure braid, the beginning and the end of each strand are in the same position. Pure braid groups fit into a short exact sequence As is obvious, the sequence splits and therefore pure braid groups are realized as iterated semi-direct products of free groups. The braid group B 3 is the universal central extension of the modular group PSL(2, Z).
In that way, a braid group operation gives Ψ P which implies that Ψ P is in the center of B 3 . It is a wave function of protein transition anyon state. If G = LB BIO n acts on X t If G = LB BIO n acts on Y t , one gets From above, it follows that the elements d i and d i+1 exchange places in DNA strand by an analogy with genetic variation. One can check that the braid group relations are satisfied and this formula indeed defines a group action of LB BIO n on X t .

Affine Loop Braid Group in Public Key Cryptosystem of Protein Folding
Topological cryptography for biology releases a self versus non-self protein-protein authentication as a public key, a non-protein-coding piece of DNA. An example of such an area of DNA is mitochondrial DNA (miDNA). The genetic codes of tRNA in miDNA and in the chromosome are the same codes but in a different places. The methyl transfer to AdoMet [32] of tRNA represents the role of knot protein as the RSA algorithm in public key biological cryptosystem. The encryption of the moduli state of a protein comes from 20 modulus states in tRNA. The real function of tRNA is a hash function for encrypting the digital signature of biological information.

Public Key Biological Cryptosystem in Protein-Protein Docking
The system of protein-protein interaction in an immune system can be realized as a public key biological cryptosystem. It is the main application of the loop braid group relations in self and non-self protein-protein authentication in a cryptosystem of antigen-antibody recognition.
Secret information is encoded in the genetic code and immune system with some help from y-RNA.
There exists a normal form for the elements of LB BIO n in the terms of generators σ n . It induces the application of loop braid groups in the genetic code to pubic key biological cryptosystem. In order to explain the authentication mechanism of the digital signature in the system of protein-protein interaction, we have to recall the natural phenomena in the immune system. It is a classification between self and non-self cell death program in apoptosis. Unlike the use of a two-factor authentication certificate, the user needs to know the digital signature to authenticate the self and non-self certificate when a cancer cell or virus attacks the immune system. There must be a private key associated with the noncoding y-RNA to transfer the certificate through a nuclear pore, so the immune system will be able to access the defense authentication system.
We are generated by the generators σ l−1 and the generators σ l+r−1 respectively. We have commutative properties between left and right supersymmetry of loop braid group as an asymmetric cryptosystem for the authentication public and private keys in the RSA algorithm for biology (see Figure 6c), for any Ψ a ∈ LB BIO,L l and Figure 6. (a) The picture shows eight anyon smallest subunits in biology. It is a representation of eight spinor fields in the genetic code of living organisms. We define them as a flying ring in the unification between three elements from DNA, RNA, and protein. (b) The picture shows three layers of DNA, RNA, and protein in the asymmetric cryptosystem model for biology. (c) The picture demonstrates the encryption and decryption modules in the public key cryptosystem. The non-protein-coding DNA part is a source of the private key. The telomere is a biological clock subunit; it is equivalent to the modulus part in a moduli state space model. Let a one-way function of the gene expression be defined by It is a one-way function to carry public-private keys given by a pair (Ψ a , σ x ). Now we define key agreement in the protein docking system between the proteins A and B with the help of the braid group version of the Diffie-Hellman key agreement in a quantum cryptosystem for biology. We have the following steps.

1.
Preparatory step: A suitable pair of integers (l, r) are chosen, and a sufficiently complicated (l + r)-loop braid σ x ∈ LB BIO l+r is selected and published through the gene expression.

2.
Key agreement: A key is shared through the protein folding by performing the following steps each time a shared key is required: • Protein A chooses a random secret loop braid Ψ a ∈ LB BIO,L l and sends σ y 1 = Ψ a σ x Ψ a −1 to B.
• Protein B chooses a random secret loop braid Ψ b ∈ LB BIO,R r and sends σ Protein A receives σ y 2 and computes the shared key Ψ K = Ψ a σ y 2 Ψ a −1 . • Protein B receives σ y 1 and computes the shared key Thus both proteins A and B obtain the same loop braid group element with the same curvature. Therefore they can dock to each other; otherwise the system will not be in equilibrium and the moduli state space for the gene expression with control equation will recursively loop back to send the key again.

Moduli State Space Model for Protein-Protein Interaction
The protein-protein interaction can be realized as a share key agreement in the biological cryptosystem, if we consider the central dogma as the short exact sequence The moduli state space model for protein-protein interaction with G-quadruplex, t G = t NG , in the D-brane plane, comes from this short exact sequence: If we approximate this equation with the gene expression at first time period t 1 = t 1 G with a short exact sequence, then β 1 is the first genotype at time period t 1 = t 1 G in the equilibrium moduli state of the gene expression: We have a recursive loop of reversed transcription of jumping gene of retrotransposon transition state y t with repeated geneon state x t G by so we have For no genetic variation, µ t = 0, we have a recursive loop of repeated jumping gene of retrotransposon in the non-protein-coding DNA as The initial value equation of moduli state space model for the non-protein-coding DNA with starting state of retrotransposon y 0 is Therefore we have We have moduli state space model for DNA, The definitions of public and private keys in a biological cryptosystem are the following.
Encryption module: Given a biological message authentication code by y-RNA, tRNA, and ncRNA as a hash function with spinor state m ∈ (Z/2) k and the public key (σ x , σ y ), protein B chooses a loop braid group Ψ b ∈ LB r at a rate of random mutation t . Biological cipher state is Decryption module: Given the biological cipher state (Ψ c , Ψ d ) and the private key Ψ a , compute the Because Ψ a and Ψ b commute to each other we can recover the information that transmission from DNA to protein by encryption process in RNA by

Empirical Analysis of Cancer Gene Signature
Our empirical analysis is based on a novel algorithm based on (ITD − IMF)chain 1 (n) in time series prediction introduced in [10]. The topological approach to computation on DNA time series data was successfully applied in the empirical mode decomposition of the Chern-Simons current related to genetic variations of viral glycoprotein gene and host T-cell receptor gene [4].
We also calculated the Chern-Simons current in the genotype. We used the genetic code of telomerase RNA component TERC from GenBank [33] (see also Appendix C) and the genetic code of ncRNA in telomerase (Appendix D). The results of the plots for the Chern-Simons current in the genetic code of PD1 and ncRNA are presented in Figures 7 and 8. The plots are the Chern-Simons currents in genetic code, for different values of n; they represent the wave functions of retrotransposon states in ncRNA. Higher values of n yield smoother curves; they have parabolic shapes with concave or convex curvature, as one can see from the ncRNA or protein PD1 plots. Similar plot results for the Chern-Simons current can be used to detect a cancer signature, as can be seen in Figure 9.
For that purpose we downloaded data of 36 samples of cancer genes (Table 1) from the online database Catalogue of Somatic Mutations in Cancer (COSMIC) with sample name PD7301a and COSMIC sample ID COSS1540693 [34]; the tumour location is in the bone; the pelvis (chondrosarcoma; central); see also Appendix E for the AX1 gene. In Figure 9, one can see the difference in the calculation of the Chern-Simons current in the ALX1 gene for cancer (red) and normal (blue) gene versions; the change is better visible in the computation of (ITD − IMF)chain 1 (8). In Figures 10 and 11 are presented results for the tensor correlation of the Chern-Simons current from (ITD − IMF)chain 1 (1) and (ITD − IMF)chain 1 (6) between normal (blue) and cancer (red) genes in selected 28 samples of bone cancer. In the region of amino acids 700 and 1000, the mutation area is noticeable. In Figures 12  and 13 are plots of the (ITD − IMF)chain 1 (1) and (ITD − IMF)chain 1 (8) of the Chern-Simons current in the genetic code for mutation for samples 33-36 presented in Table 1. Next are presented the images of the spectrum of tensor correlation of 28 cancer genes, computed for the amino acid number 531 ( Figure 14) and the amino acid number 9980 ( Figure 15). The computation uses canonical tensor correlation structure from (ITD − IMF)chain 1 (1)) to (ITD − IMF)chain 1 (6) from 28 genes listed in Table 1 with matrix size 6× 19,244. The number 19,244 comes from total sum of all amino acids in 28 samples. One can see the differences in the color patterns. The visualization of comparative difference between the cancer and normal genes can help in building the databases and the classification of different types of genes in the future.  Figure 11. (a) The picture shows the tensor correlation of the Chern-Simons current from (ITD − IMF)chain 1 (6), between normal (blue) and cancer (red) genes in 28 selected samples of bone cancer for the first 1000 amino acids. We can see that the signature of mutation appears between the amino acids 400 and 600. On the plot (b), if we zoom the region between amino acids 700 and 1000 we can also notice some signs of the mutation in that area.  Table 1.  Table 1. The same image at the same position of amino acid which we computed from selected 28 normal reference genes. The computation uses canonical tensor correlation structure from (ITD − IMF)chain 1 (1)) to (ITD − IMF)chain 1 (6) from 28 genes listed in Table 1. We used slide window with the size 6 because we had total samples of amino acids with matrix size 6× 19,244. The number 19,244 comes from total sum of all amino acids in 28 samples.

Discussion and Conclusions
We use loop braid group to define the flying ring in the (3 + 3) extra dimensions of triplet state of transition between the DNA, RNA, and protein. These three molecules are interconnected by combining them into a four-dimensional coordinate system. We can unify them as a loop braid group in genetic code as the smallest subunit in the translating of the DNA sequence into a folded protein. The flying ring is just an 8-spinor field in the homogeneous coordinate of supermanifold in the genetic code.
The interaction between the D-brane of a telomere in the non-protein-coding DNA and the anti-D-brane of the telomerase enzyme with embedded noncoding RNA is the source of the biological clock and cell aging in a living organism. This process of the biological clock generator is the rotation of a homogeneous coordinate in the gene expression. The rotation of the clock produces a spinor field in time series data as a signal for cell transduction. It is the information inside codon which generates the public and private keys in quantum cryptography along the geneonic state. The signal of protein expression is expressed as immune receptors and signal transduction with the public and private keys. The protein state belongs to a superspace P. The gene expression in the eukaryotic cell with mRNA will have approximately eight exon states with β i cocycle geneonic states.
The other part of noncoding RNA in mRNA is an intron state, denoted as α i cocycle of retrotransposon state in ncRNA space R * . The function of intron state is to control the gene expression as public and private keys to deliver the sequence of mRNA out of the nuclear membrane with the help of y-RNA. The keys will match to each other outside the nuclear pore in ribosomes in the rough endoplasmic reticulum for the process of t-RNA to match the public and private keys inside ribosome for the translation process of producing the amino acid. It might be a real function of the non-protein-coding DNA. The acceleration of the life clock is caused by the imbalance of collisions between the telomere and telomerase by passing the energy-momentum tensor in protein-protein interactions.
In this work, a new theory of quantum cryptography for biology has been developed. In particular, we used the loop braid group for biology to visualize the G-quadruplex with six repeated sequences in a telomere. We plotted the Chern-Simons current in the genotypes of all samples to detect the physiology in the biological time series data. We found that the graph of cancer gene has a degree of nonstationary time series data more than the normal gene. The tensor correlation network for 36 cancer genes has been generated and the image of spectra of geneon state in a cancer gene has been produced. We perform RSA image encryption to see the private and public keys over the image of geneon. In our approach, quantum cryptography is applied in biology to find the pattern of mutation genes by using an image encryption method.
We think that this new methodology could evolve into the more practical visualization of genes in the future, instead of showing the sequence of amino acids. DNA (and genes) is permanently changing, it has a dynamic structure. The mutations on the gene level affect the genes, which subsequently can lead to production malfunctions of various proteins. Gene-level mutations are caused by different types of processes. It may be a swapping of one nucleotide for another one, a deleting or omitting of a specific nucleotide from the strand, or an inserting of an extra nucleotide into the replication strand. Besides the single codon substitutions, on the level of codons there exist other categories of mutations. The incorrect amino acids can be inserted into a protein molecule, the synthesis of a protein can be prematurely terminated, and the mutated codon can code the same amino acid as the non-mutated codon. As a consequence of DNA constant change, the difference between the mutated and normal versions of the gene changes over time. We can train our algorithm on the different versions of the genes, on normal genes with non-modified sequences of nucleotides or on their mutated versions as cancer genes, which code different structures than the normal ones.
With our model we are looking for the effective change of topology of proteins, i.e., the stable structures of noncoding DNA binding to genes; this procedure selects only certain structures that can be detected. We hope for the application in cancer genes, but at first the genes on which the measurements can be made must be identified by genetic experts. However, we provide the tool for the terminal measurements incorporated into more complex biophysical experiments. In the future, the application of the methodology can be developed into more practical automatic detection of mutated genes, either for medical purposes or for imaging production of the cancer genome project.
In conclusion, we applied a direct application of algebraic topological quantum cryptography and quantum field theory with loop braid group theory directly to perform the analysis of the ALX1 cancer gene for the detection of a cancer signature. We proposed and used a new theory of quantum public key encryption in a cancer cryptosystem. The result of the computation of the RSA algorithm for public and private keys in normal and cancer genes appeared to be significant for the biological cryptosystem. Finally, we used the Chern-Simons current to plot the tensor correlation and compute the image of spectra for cancer signature.
The relevant point of our approach is the fact that quantum field theory developed for physics, in particular, the Chern-Simons theory, is extremely useful for describing biological systems [35,36], despite the huge difference in orders of magnitude between fundamental particles and genetic code. However, our mathematical method can also be used not only for mapping and determination of the genomes of viruses and subsequent prediction of their mutations, but also in astrophysical and cosmological systems. Following symmetry in biological nature gives us a chance to also compute the structures of biological systems in our description of the mutation genes.
Author Contributions: Conceptualization: S.C. and R.P.; investigation: S.C., R.P. and E.B.; writing-original draft preparation: S.C. and R.P.; writing-review and editing: S.C. and E.B.; visualization: E.B.; project administration: S.C. and R.P. All authors have read and agreed to the published version of the manuscript.

Acknowledgments:
This article is based upon work from CANTATA COST (European Cooperation in Science and Techonology) action CA15117, EU Framework Programme Horizon 2020. S.C. acknowledges INFN for partial support. R.P. would like to thank the TH division in CERN and BLTP in JINR, Dubna for hospitality.

Conflicts of Interest:
All authors declare that they have no conflict of interests.

Appendix A. Path Integral and Quantum Biology Theory
New definition of a Dirac operator for the gene expression can be written with the help of Hamiltonian H(p, q) as where E d n is an eigenvalue of transition energy with a transition state d n , n = 1, 2, 3, · · · in DNA. The operator D has two parts. The kinetic part is the source of energy and the momentum tensor in the genetic code. The potential part is the source of an RNA state or a copy state of DNA to RNA. The first part is the momentum square of protein-protein interaction, , where β is an behavior field in the protein folding structure; p is a momentum in protein p = −i ∂ ∂p 3 . The imaginary number means a hidden state in DNA structure associated with the protein. The second part is the potential part of RNA molecule. The interaction of RNA to DNA is expressed in this part analogically with potential energy in the Dirac operator.
The path integral [37] in loop braid group can reproduce the Schrödinger equation for gene expression. We can derive following well known result in quantum mechanics; it is true also in quantum biology because the quantum concept in this work proves the validation of a biological system of gene expression in the form of a path integral over a loop braid group.
Let r 1 , r 2 , and r 2 be three alphabet letters of the anti-codons of t-RNA. They are associated with the smallest elements of proteins as the momentum component p = (p 1 , p 2 , p 3 ) for the component of angular momentum in the gene expression:d We have d = (d 1 , d 2 , d 3 ), r = (r 1 , r 2 , r 3 ) and p = (p 1 , p 2 , We have energy and momentum tensor T ik defined by Lagrangian L(d, r, p) = L(q, p), p = q , and q DNA state q = d or RNA state q = r. Since r = 1 d in an equivalent relation, q is the generalized coordination of the configuration space of a Lagrangian. It maps from the tangent of a super manifold of living organism to a scalar field as a source of energy and momentum tensor in the genetic code. If the path integral along the configuration space is parametrized by a path as a functional f , We generalize three angular momenta p µν , of three states of (d, r, p) as the unification of central dogma from d → r → p with the flow of energy and momentum tensor along the gene expression in every possible direction of path integral. Let x = (x 1 , x 2 , x 3 , x 4 ) be states of four genetic codes in DNA; it is an element of SO(3, 1). Let p i be a protein state as the momentum vector and let t be a rate of an evolutional field. The angular momentum tensor is analogous with a source of gene expression in DNA state. That blends the molecular structure of DNA into curve d µν = d i dp k − x k dp i = 1 If we rewrite the relation r = p(d) in a scalar field to a vector field with a cross product analogous to a vector field in Lie algebra, we will get the induced field of transition state in DNA as the interaction between RNA and the protein (the behavior of interaction between two docking proteins) with d = r × p for the induced ncRNA in ancient virus t-RNA with their associated amino acids in proteins for the codon alphabet in the DNA template.
If t = ∆d is an evolutional field in the behavior field of host cell DNA, we can formulate the so-called uncertainty principle in the gene expression Let p be momentum of protein in active layer of docking system and let p * be a protein in passive layer; we have the interaction as transition state in DNA as kinetic energy in gene expression The explicit form of a wave function in DNA contains the momentum for three states of momentum p = (p 1 , p 2 , p 3 ) in one amino acid as a linear combination where ϕ i = e 2πip i d .

Largangian of Path Integral Formalism for System of DNA-RNA Transcription
We start from the formulation of Lagrangian of the system of DNA-RNA transcription with extended central dogma. Let Lagragian L = L(p, r, d) = L(p(d), r(d)) be a map from the tangent of the supermanifold of a living organism with eight spinor states T p M =: S 7 → HP 1 with a fiber of S 3 . It is a source of Hopf fibration as pa ath integal line and it is exactly a line of the loop braid group in genetic code parametrized by d ∈ D. We choose the coordinates d, r, p ∈ S 3 = SU(2) = H satisfying the equation ∂r ∂d Definition A1. We have the Euler-Lagrange equation for gene expression in the system of extended central dogma with L = L(p, r, d) = L(p(d), r(d)) given by The action of path integral to minimize the geodesic curvature of the path in DNA is α(d) = It is a flow of genetic code as a line in braid group in genetic code. In this work we choose r = p and d is a time variable.
The path in quantum biology theory means the string of protein folding which can deform continuously from primary protein structure to secondary protein structure by using homotopy equivalent relation in loop braid group (see Figure A1). If the sequence of protein folding starts from discretized time scale in the path integral, one can choose t 0 := d 0 as a starting point of gene expression in DNA and t n = d n as an ending point. It is a chosen localized base point in braid group as line in our parameter of localized path in Lagrangian. When these lines intersect each other we have path ordering of Wilson loop [38] in genetic code. The homotopic path is the homotopy equivalent of a line in the path integral from an equivalent class of curve map from the Lagrangian of the tangent of the manifold to the based point of the intersection line between fibration. We define an element of loop braid group in the path integral over loop braid group in genetic code, σ i ∼ σ j , if and only if there exists a homotopic path as a map in [α(d)] ∈ [S 1 , * ] from α : L(d, r, p)) → * ∈ R.
The localized path in loop braid group d a = d 0 < d 1 < · · · d i · · · < d n−1 < d n = d b can be divided into n smaller subintervals d i − d i−1 , where i = 1, · · · , n. We can permute this path ordering by using Wilson loop over loop braid group. The source of mutation is and it is called a source of uncertainty principle in so-called quantum biology theory, and a source of evolutional field in the form of three types of behavior fields in the genetic code. An approximation for the path integral can be computed to be proportional to number of gene expression in tRNA. The codes of DNA (codon state) and tRNA (in anticodon state) are assumed to be the same length (without exon and intron involved) in this simple system of DNA-RNA transcription. We have We define a probability amplitudes of gene expression as a source of protein folding by r b , d b |r a , d a . It is a spectrum of gene expression with the probability to find the gene particle (geneon) at d a in the initial state of RNA r a and at d b in the final state r b . In the simple case of Lagrangian we have L(r, p, d) The wave propagation obtained from the path integral so-called geneon wave function in gene expression with eigenvalues as transition states in protein folding structure is where Dr denotes integration over all paths r, with r(0) = r and Z is a normalization factor. Here S is the action, given by S[r, p] = dd L(r(d), p(d)). We obtain the exponential of the action by and we can use the approximation Then, after the rearrangement of the terms properly, we get which is the Schrödinger equation for the gene expression. Let cell cocycle be a function of behavior in the genetic code of living organism denoted by g 00 g 01 g 02 g 03 g 10 g 11 g 12 g 13 g 20 g 21 g 22 g 23 g 30 g 31 g 32 g 33 g 00 g 01 g 02 g 03 g 10 g 11 g 12 g 13 g 20 g 21 g 22 g 23 g 30 g 31 g 32 g 33 There exists a question about the smallest subunit of genotypes as a pair of alleles in the pair of chromosomes. The smallest subunit should exist in two DNAs, one from a mother's chromosome and other from a father chromosome. It is the appearance of four states A, U, C, G from one DNA inactive part of allele AA, Aa, aa and other four hidden states from a passive part of other one from father DNA, denoted as A * , U * , C * , G * . The induced four states from homozygous and heterozygous to other four states in the proof is a Hodge star operator in type-III Yang-Mills field * F µν = −F µν . It is a source of a supersymmetry and it answers why the base in DNA molecule has only right-hand chiral symmetry. The left-hand one is in the hidden state of passive layer N A, NU, NC, NG. Let the Yang-Mills field induced from the interaction of hidden evolutional field be F µν = t | * t . This field comes from an adaptive behavior field in the genetic code in the form of a mutation.
Definition A2. For immature behavior field, we use finite state for physiology of time series to accepted protein folding pattern defined by where σ y = σ 2 is a Pauli spin matrix along y-axis.
Definition A3. For a maturing behavior field, we use a finite state for the physiology of a time series for the accepted folding pattern defined by where σ z = σ 3 is a Pauli spin matrix along z-axis.

Definition A4.
Let ω be a instinct behavior of "junk" DNA. From the well known induced relationship between spinor fields [σ y , We introduce three evolutional behavior fields of adaptive behavior of changing curvature in the protein folding where + means a positive adaptive behavior of a protein, − is negative adaptive behavior, m is a maturing behavior, im is an immature behavior, and ω is an instinct behavior. A 1 is a connection behavior of a maturing protein, A 2 is a connection behavior field of an immature protein, A 3 is a connection behavior field of the instinct-adaptive behavior field of the protein.
Using the definitions above we have It is a vierbien of a transition of cell cocycles in two rounds with the active and passive side of the DNA-RNA transcription (see Figure A2). In order to compute the interaction of the behavior field in the genetic code in terms of evolution, we simplify by analogy a multiplication of a predefined number with vierbein: We use the well-known formula for the connection A µ := Γ µ ij for the enegy-momentum tensor in the genotype. It is the transition between the ghost field and anti-ghost field in the active and passive layer of the DNA-RNA transcription. With the existence of a knot protein, the underlying knot of quantum observation over a superspace of genetic code is a modified Wilson loop. It is a hyperbolic knotted-over loop space in biological time series data. The allele in the genotype is a partition function of the active site in DNA. The partition function of a lattice gauge theory for biology is If we compare this expression with an effective action of Josephson junction in a biological system with critical Chern-Simons current in biology given by where J µ c is a supercurrent, ϕ = ϕ R − ϕ L is a phase difference between D-brane and anti-D-brane sheet in the DNA-RNA transcription with underlying two-child manifolds of knot state with partition function Z(β) = e −βF ; the free energy of our system is measured by a holomony of Wilson loop inside the Chern-Simons current. We have expansion of the Wilson loop over Green function by µ 1 ,µ 2 ,··· ,µ n (x 1 , x 2 , · · · , x n ).
The invariance of Wilson loop can be obtained from the Laurant polynomial in Witten invariance with an explicit form of the Chern-Simons current in which the closed form of sine is an analogy with the phase shift in a current of energy-momentum tensor, so we have W 4 1 ,ρ,SU(2) (C L ) = β 4 with a strong coupling constant β = β ω ± ,im ± ,m ± . We have a cell cocycle over behavior field by The path integral in a superspace of genetic code is a Green function of the propagator of a behavior field in an immune system A µ along a principle fiber space in secondary protein structure. Let x t 0 be a time series of gene expression at time t 0 and let et x τ be a biological time series at time t = τ. The path integral defines a gene expression equation with the coupling state between geneonic eight states between active and passive layers of a behavior field with parallel transport A µ .