Knot Theory for Proteins: Gauss Codes, Quandles and Bondles

Proteins are linear molecular chains that often fold to function. The topology of folding is widely believed to define its properties and function, and knot theory has been applied to study protein structure and its implications. More that 97% of proteins are, however, classified as unknots when intra-chain interactions are ignored. This raises the question as to whether knot theory can be extended to include intra-chain interactions and thus be able to categorize topology of the proteins that are otherwise classified as unknotted. Here, we develop knot theory for folded linear molecular chains and apply it to proteins. For this purpose, proteins will be thought of as an embedding of a linear segment into three dimensions, with additional structure coming from self-bonding. We then project to a two-dimensional diagram and consider the basic rules of equivalence between two diagrams. We further consider the representation of projections of proteins using Gauss codes, or strings of numbers and letters, and how we can equate these codes with changes allowed in the diagrams. Finally, we explore the possibility of applying the algebraic structure of quandles to distinguish the topologies of proteins. Because of the presence of bonds, we extend the theory to define bondles, a type of quandle particularly adapted to distinguishing the topological types of proteins.


INTRODUCTION
Folded linear molecular chains are ubiquitous in biology. Proteins and nucleic acids are linear polymers responsible for most cellular functions, for the inheritance of biological information, and are subject to changes during evolution and pathologies [Dob02,RC18]. These chains often fold to function and their 3D structure contains information about their dynamics, evolution, and inter-molecular interactions and can be used for designing drugs [KLN + 19, AKB + 17]. Geometric and chemical properties of folded proteins and genomic DNA have been widely studied using various methods including NMR spectroscopy, X-ray crystallography, chromosome conformation capture, and mass spectrometry among others. Topological properties of these molecules have remained relatively unexplored due to lack of a relevant conceptual framework. Knot theory was successfully applied to study proteins and nucleic acids and protein and DNA knots were studied using various experimental techniques including nanopore technology and probe microscopy [TL03,SSSC08,SF13,MRJ10]. Despite being interesting and innovative, these studies had a limited impact on protein science as the majority of identified proteins fall into one topology class, i.e., the unknot [SRM + 12]. Thus, standard knot theory cannot be effectively used to classify proteins [MvWT14]. Another major shortcoming of the standard knot theoretic approach is that intra-molecular interactions or contacts are ignored. These interactions drive the folding of the molecular chains and are functionally important [MR15, MTM14, HST + 17, SHM + 17, HSM19]. Thus, there is a need for a new topology framework that includes intra-chain interactions and is able to classify fold topology of biomolecular chains and in particular the proteins. This paper presents a new knot theory for folded linear molecular chains and looks to classify the topological structure of proteins through the application of certain aspects of knot theory. More specifically, we will apply a modified singular knot theory, Gauss codes to keep track of the structure, and an associated singular quandle called a bondle to distinguish structures.
Proteins are continuous linear molecules with the ends unbonded, so we will look at them as linear segments embedded in three-space. Protein structure for a single protein is formed on three different levels, shown in Figure  1. The Primary Structure is defined by the ordering of the amino acids, or building blocks of the protein. These amino acids are bonded in a sequential chain from the N-terminus to the C-terminus, with each amino acid presenting an exposed R-group that can interact chemically or through the electrostatic effect with other R-groups and other molecules. The Secondary Structure is defined by the coiling or local structuring of the amino acids. This is where structural patterns appear such as β-pleated sheets and α-helices, both held together by hydrogen bonds. Tertiary Structure is defined by the interactions between different R-groups or backbone interactions, forming a structure for the entire protein.
A knot is an embedding of S 1 , the circle, into three dimensions, called a conformation, given by a function f : [0, 2π] → R 3 where f (0) = f (2π). We define an equivalence on conformations of knots, considering two conformations equivalent if there is an ambient isotopy from one to the other. This means that we can deform the one through space to the other without passing the knot through itself. As is common, we will use the word "knot" to represent both a given conformation of a knot and the equivalence class of conformations corresponding to the given conformation, clarifying which is which when necessary. Knot theory can be extended to singular knot theory, where we allow a finite number of singular points where two points on the circle are sent to the same point in 3-space by f . (See [Dan18] for more on singular knots.) Allowing singularities will provide us with the ability to model intra-molecular bonds.
While it is often preferable to describe a knot in three dimensions, it is not always tractable. Therefore, we project the knot in a particular direction onto a plane, obtaining a projection. We consider only regular projections where a finite number of pairs of points on the knot are identified with each other and result in what we call crossings. We keep track of which of the two points in the pair is the top one. A projection is essentially a shadow that retains information at the crossings of the strands. Using a defined set of rules called Reidemeister moves, we can change one projection of a knot to any other projection of that same knot. In a two-dimensional projection, we define a classical crossing as a place where one strand goes over another, and if we are allowing singularities, a singular crossing where two or more strands intersect each other at a point in the three-dimensional conformation.
We can similarly think of proteins as a conformation of a line segment [0,1] in three dimensions, which we will call a protein model, with an associated function f : [0, 1] → R 3 . If we utilize the same equivalence of ambient isotopy, then every protein can be disentangled and they are all equivalent. But once we include singularities reflecting bonding of points along the conformation, this will no longer be true.
In a protein, we define singular crossings to exist where a protein has intra-chain interactions (also called contacts). These contacts take one of two forms. The first is covalent bonds, defined as two atoms sharing electrons. Of special interest to proteins are disulfide bridges, in which two thiol R-groups, made of one sulfur and one hydrogen, bond and release their hydrogen atoms. The second type of bond is formed through noncovalent interactions. A major form of non-covalent interactions are the electrostatic ones, particularly in the form of hydrogen bonds, which are partially electrostatic, but often come in multiples, making their strength significant enough to control the protein's structure. Hydrogen bonds mediate the interactions between beta strands and the formation of alpha-helical structures.
Proteins differ from knots in that they have two endpoints that are not connected. One approach to this difference has been work on modelling proteins as knotoids, which are the projection of the conformations of a linear segment [DTRG + 19, GGL + 17]. Reidemeister moves (which we will look at in the next section) apply in the theory of knotoids, but in this theory, the endpoints cannot be moved across other strands in the projection. Although this may be useful when dealing with proteins that are somewhat rigid and have so-called lassos, here we desire a general theory that allows the movement of endpoints across strands in a projection. Thus we do not consider proteins as they relate to knotoids.
In finding a notation for protein structure, it is true that proteins have variable flexibility and restricted length, providing more limitations than the ones we place on curves in 3-space when defining knots. But the goal of this paper is topological in nature. Thus, we do not capture the full sense of rigidity or steric hindrance in a protein. Since we are allowing for the deformation of a protein's strands in the following sections, we will at times allow the same for the ends of the strand. The notation defined in this paper looks to balance simplicity and mathematical utility with chemical precision.
In Section 2, we discuss Reidemeister moves, which are moves one can do on a projection of a knot to obtain a new diagram of the same knot. We extend them to allow features present in proteins, including bonds, endpoint β-pleated sheets and α-helices.
In Section 3, we introduce Gauss codes, which can be used to describe in symbols a projection of a knot. We extend them to proteins.
In Section 4, we introduce quandles, which are algebraic objects that can be used to distinguish between different knots. In Section 5, we introduce singquandles, which are an extension of quandles that have been used to distinguish knots with singularities. We further extend this idea to the idea of a bondle, which is a quandle that can be applied in the presence of bonds.
In Section 6, we introduce the oriented bondle, which seems particularly suited to distinguishing between the topological types of proteins. We then identify several families of oriented bondles. In Section 7, we provide several examples of pairs of proteins that can be distinguished using oriented bondles.

REIDEMEISTER MOVES
Reidemeister moves are a set of changes to the combinatorial pattern that is a projection of a knot. Planar isotopy is a deformation of the projection that does not change the combinatorial pattern. The critical result from [AB26] or [Rei27] says that two knots K 1 and K 2 are equivalent if and only if there exists a sequence of Reidemeister moves and planar isotopy in the plane that transforms a projection of K 1 to a projection of K 2 . The three Reidemeister moves are depicted in Figure 2. We refer to the monogonal face found in a Type I Reidemeister move as a kink. While these three Reidemeister moves are able to completely describe equivalence in classical knots, for our purposes we need to add in singular crossings to represent the covalent and hydrogen bonds that proteins form with themselves ( [Dan18]). The use of singular crossings requires an additional two Reidemeister moves as shown in Figure 3   As appearing in this illustration, we denote a singularity by a small rectangle with two parallel edges on two strands. By drawing singular crossings in this fashion, we show that the strands do not cross over each other, but are instead bound together to more closely replicate the structure we see in proteins. By doing so, we do not lose any of the structure or properties we would hope to retain.
In a protein, a β-pleated sheet consists of multiple segments of a protein that run parallel to each other, roughly in a plane, with hydrogen bonds connecting each segment to its adjacent segment in multiple places. If we collapse the sheet to a point, this functionally looks like a singular crossing with more than two strands. Therefore, in order to represent β-pleated sheets, we must extend singular knot theory to contain multi-singularities. We define a multi-singularity as a place where two or more strands intersect each other at a single point, as shown in Figure 4. We can then incorporate this into the moves shown in Figure 3 and define the set of singular Reidemeister moves shown in Figure 5 to describe topological isotopy.
Note that when applying a Type V Reidemeister move to a multi-singular crossing, we allow any number of segments to be included, even though we have only depicted three segments to simplify the illustration. When the singularity is flipped, this causes a half-twist of all of the segments above and below the multi-singular crossing in the process.
Whereas in a Type IV move we only explicitly allow a segment to pass through a multi-singularity from left to right, we can use the given Reidemeister moves to show that we can slide a horizontal segment past a multisingularity from top to bottom a well, as shown in Figure 6.
Proteins also contain α-helices, which we can represent in multiple ways. An α-helix appears where a segment of a protein coils, with hydrogen bonds holding the coils together, as shown in Figure 1. For now, we simply say that they do not add topological structure, but are important in protein function. In a projection, we allow α-helices to slide along a segment, freely FIGURE 5. Types IV and V Reidemeister moves for multisingular crossings.
Type IV FIGURE 6. A vertical Type IV Reidemeister move produces a horizontal Type IV Reidemeister move. traversing a classical crossing, but not to pass through a singular crossing. Thus, we create a Reidemeister move Type VI to allow α-helices, depicted as the small set of jagged lines, to pass over or under other segments as in Figure 7.
Since a protein is a continuous strand, we are able to equate its projection to a segment of a knot, allowing us the aforementioned Reidemeister moves. But because the protein has two unbonded endpoints, denoted N and C, that are free to move around, we need to be able to slide these endpoints past another strand. Therefore, we define a final Reidemeister Type VII move for this action as well as in Figure 8. We would like to prove that this set of seven moves captures all possible equivalences between projections of protein conformations. To achieve this, as is done in the knot case as well, we replace smooth conformations with piece-wise linear conformations, meaning that our conformation can be represented by a finite number of line segments glued end-to-end.
Theorem 2.1. Given two protein conformations, they are equivalent if and only if there is a sequence of planar isotopies and Reidemeister moves from this set of seven moves to get from a projection of the first to a projection of the second.
Proof. As previously mentioned, it has already been proved that for a classical knot, planar isotopy and Type I, Type II, and Type III moves suffice to represent ambient isotopy in a knot (see [Kau89] for a readable proof). The proof uses triangle moves to realize ambient isotopy between two piecewise linear knots. A triangle move is realized by taking a solid triangle that intersects the knot in one or two edges and replacing those line segments on the knot by the non-intersecting edges on the triangle, as shown in Figure  9. Any isotopy we attain from deforming the original conformation can be represented by triangle moves. When we project down to a projection, one can show the triangle moves appear as planar isotopy and Reidemeister moves.

FIGURE 9. Example of Triangle Move
A spatial graph is a conformation of a graph (consisting of a collection of edges sharing a collection of vertices as their endpoints) in 3-space. A rigid vertex graph further posits that adjacent edges coming into a vertex cannot twist about one another. The vertex can be flipped, intertwining the edges appropriately as in our Type V move. In [Kau89], it is shown that two conformations of a rigid vertex spatial graph with all vertices of valency 4 (called an RV4 graph) are equivalent if and only if there is a sequence of planar isotopies and Reidemesister moves of Types I, II, III, IV, and V from a projection of one to a projection of the other. This proof also considers how triangle moves impact the projection.
Our situation for a protein conformation has three differences from this one. First, we allow our vertices to have an even number of edges that can be four or greater. However, this does not impact the proof as in [Kau89] and it goes through exactly as before.
Second, we have α-helices. These can be treated as vertices of valency two, and then the same arguments as in [Kau89] go through to generate the Type VI Reidemeister move.
Third, we have the N and C endpoints of the protein. But it is straightforward to see that a triangle move that projects to overlap an endpoint will simply generate the Type VII Reidemeister move.
Thus, all isotopies allowed in deforming our structure can be represented by triangle moves, which in the projections can be represented by the defined Reidemeister moves. Hence, the seven moves are sufficient to convert one projection of a protein conformation to any projection of a topologically equivalent protein conformation.
With these seven Reidemeister moves, we are able to transform one projection of a protein to another. However, it is inconvenient to convey these transformations through diagrams. Encoding diagrams as strings of text allows us to both interpret and transmit the information. Therefore, we turn to Gauss codes to represent protein structure in text form.

GAUSS CODES
Gauss codes are a means to represent knot projections using only symbols instead of diagrams. Using a Gauss code, we are able to easily go back and forth between the code and a projection while also being able to change entries in the code when Reidemeister moves are applied.
Since proteins are constructed and written from the N-terminus to the Cterminus, there is a definitive start and end point to proteins. Therefore, we start the Gauss code of a protein projection with N and end the code with C. We also have a natural orientation to the protein. Crossings are oriented as in Figure 10. As with traditional Gauss codes, when following a strand, if the strand crosses over another strand, we denote this with an O for over. Similarly, when traversing beneath another strand, we denote this with a U for under. We assign an orientation to the crossing denoted by a superscript of ±.
If an α-helix appears, we denote it as α with a superscript of + if the helix is right-handed (coils clockwise) and − if the helix is left-handed (coils counterclockwise). Bonds are written as B, and β-pleated sheets are written as β. Strands in bonds and β-pleated sheets are labeled with a superscript of ±, using + if the strand runs parallel to the strand that first occurs in the bond or sheet (so the first strand always receives a +), and − if the strand runs anti-parallel to the strand that first occurs.
With β-pleated sheets, strands are numbered with a subscript. The zero strand is defined as the first strand that appears in the protein's sequence. Numbers are assigned as sequential integers to the left and right of the initial strand, with positive integers appearing on the side to which the second strand appears in the sequence, and all strands on the opposite side of the initial strand are defined as negative, as in Figure 11.
Finally, each crossing, α-helix, bond, or β-pleated sheet is denoted with a sequential numbering based on its first appearance in the protein. Figure  12 gives an example of a Gauss code for a complete protein projection.
When we apply Reidemeister moves to a given projection, the move will be reflected in certain changes to the Gauss code. For instance, a Type I move inserts On ± U n ± or U n ± On ± into the Gauss code at the relevant point. Similar operations hold for all the Reidemeister moves.
For example, for the protein projection appearing in Figure 12, we could apply a Type VI Reidemeister move to slide the α-helix out of the bigon bounded by crossings 5 and 6, and then remove the bigon by a Type II Reidemeister move to result in the Gauss code N β1 But caution should be exercised. The corresponding operations on the Gauss codes do not always correspond to actual Reidemeister moves. For example, to do a Type II Reidemeister move, we must have two strands of the projection that are on the same complementary face of the projection. This is not visible from the Gauss code.

QUANDLES
A knot invariant is a map I : K → S from all knot diagrams K to some set S such that for any two projections K 1 and K 2 of the same knot type, I(K 1 ) = I(K 2 ). The set S can be a collection of integers, groups, polynomials or other mathematical objects. An invariant is said to be a complete invariant if the converse is true, which is to say I(K 1 ) = I(K 2 ) implies K 1 and K 2 represent the same knot type.
A quandle is an algebraic object that was introduced as an invariant for knots in 1982 independently in [Joy82] and [Mat84]. It has turned out to be a particularly effective means to distinguish knots. For more details on quandles, see for example [EN15].
Definition 4.1. A quandle is a set X with an operation £ : X × X → X such that the following three conditions are satisfied. (1) For all x ∈ X, x £ x = x.
A slightly more restrictive algebraic structure than a quandle is a kei, also called an involutory quandle.
Definition 4.2. A kei, or involutory quandle is a set X and operation £ : X × X → X that satisfy the following three conditions. (1) For all x ∈ X, x £ x = x.
Note that the only difference is that for an involutory quandle, £ is equivalent to £ −1 .
Depending on the situation, as we will discuss, one or the other of these algebraic structures may be the more appropriate to apply.
A coloring of an oriented knot projection by a quandle is an assignment of a value from X to each arc, where an arc is defined as part of a strand in a projection that both starts and ends at an under crossing, but going over zero or as many crossings as we like. We require that the labels assigned to the arcs be related through the quandle operation as in Figure 13. The relevance of quandles to knots becomes apparent when we consider how the Reidemeister moves affect our labelled diagram as in Figure 14. We see that the quandle axioms satisfied by the labels ensure that the quandle coloring is still valid after the Reidemeister moves. This means that the validity of the quandle coloring does not depend on the particular projection. It just depends on the knot type.
Thus, given a particular quandle, we can generate an invariant for knots by seeing how many distinct colorings by that quandle a particular knot has. Two knots with different numbers of colorings by that quandle must then be distinct knots.
We can also drop the orientation on the knots, in which case £ and £ −1 become identical, the arrows disappear in Figure 14, and we color with involutory quandles instead. This simplifies things as we only have one operation to consider instead of two.

Type I Type II
Type III

QUANDLES AND SINGULARITIES
In order to allow for singularities in knots, the authors of [CEHN17] introduced the singquandle. We first consider the singquandle for an unoriented knot, which will be an involutory quandle that satisfies additional conditions.
An arc in a singular knot projection is a strand that begins and ends at either an under-crossing or a singularity. Given an involutory quandle coloring of the arcs of a projection, we require the labels to satisfy conditions at the singular crossings as in Figure 15, where R 1 (x, y) and R 2 (x, y) are maps from X × X to X yet to be specified.
Since the diagram in Figure 15 can be rotated by 90 degrees, 180 degrees and 270 degrees clockwise and the relation between the top pair of labels and the bottom pair of labels must be maintained, we immediately obtain certain relations that must be satisfied: x y R 1 (x,y) R 2 (x,y) FIGURE 15. Labels at a singularity.
In [CEHN17], the authors show that in the presence of singularities, the only additional Reidemeister moves necessary are those coming from sliding a separate vertical strand on the left to the right behind or in front of a singularity, or flipping a singularity as in Figure 16. These moves generate the additional relations: Definition 5.1. A singquandle is an involutory quandle, with a choice of R 1 (x, y) and R 2 (x, y) that satisfy all of the additional relations (1)-(11).
But the singularities we wish to consider for proteins are not of this type. In our case, we have bonds across two parallel strands, as in Figure 17. Such a bond does not have a four-fold rotational symmetry, but only a two-fold rotational symmetry. Thus, we have the following definition: Definition 5.2. An involutory bondle is an involutory quandle that satisfies the relations (1), (2), (7), (8), (9), (10) and (11).
Although this choice allows us to incorporate bonds into our quandle, we do not yet have a way to represent β-pleated sheets. To deal with them, we replace a β-pleated sheet by a sequence of independent adjacent singular crossings as follows.
We have already assigned positive and negative integer values to the strands in a β-pleated sheet from the subscripts of the Gauss code. Therefore, if we define the direction of the zero strand as downwards, we can define the relative heights of the individual singularities replacing a β-pleated sheet to be strictly increasing as the numbering increases, as shown in Figure 18. The bonds appear as a set of stairs, either rising to the right or left, depending on which is the positive side of the labels on the β-pleated sheet.
This transformation of a β-pleated sheet into adjacent singularities is called a segmentation of the β-pleated sheet. With this, we can still perform the Type IV and Type V moves on multisingular crossing utilizing a sequence of Reidemeister moves on order two singularites, as shown in Figure 19. Thus, no multi-singularity Reidemeister moves are needed. However, this choice for how to represent a β-pleated sheet does mean that we cannot distinguish between a protein with a βpleated sheet and an identical one that has the corresponding sequence of bonds in place of the β-pleated sheet. Bondle invariants will be equivalent for the diagrams in Figure 20. The next issue we need to consider is the endpoints of the protein model. When considering proteins, we can view them as knot segments, with the ends free to move. Although in the physical realization of a protein, ends FIGURE 20. Three indistinguishable protein models.
are sometimes tucked inside the protein or are subject to constraints and are therefore not free to move, in our model, we allow them to slide past strands in any given projection. Even for a fixed rigid conformation, as we change our projection direction, the endpoints in the projections can slide past strands, eliminating or creating crossings. Therefore, we treat the end strands as insignificant until they reach the first bond. We think of the ends as only being relevant in defining the first and last bonds, and ignoring them otherwise, as in Figure 21.

FIGURE 21. Reducing End Arcs
The final structure we need to consider is the α-helix. We view it as a sequence of n kinks, where n is the number of full rotations that the helix contains, all having either + or − crossings depending on whether it is a clockwise or counterclockwise α-helix, as in Figure 22. These kinks are referred to as residues. Just as we are unable to distinguish a segmented β-pleated sheet from a sequence of adjacent bonds, we cannot distinguish an α-helix from a sequence of kinks.
When coloring a protein with a quandle, the α-helix becomes invisible because the Reidemeister Type I move in Figure 14 allows for the removal of kinks. However, there is a generalization of a quandle called a rack that does not allow for the removal of kinks, and therefore does see the existence of an α-helix. A rack is simply a set that satisfies the second and
third axioms of a quandle but not the first. We will not pursue racks further here.

THE ORIENTED BONDLE
Since proteins do have a natural orientation, we should also consider the oriented version of the bondle. The oriented singquandle was defined in [BEHY18]. The authors showed that in addition to the four traditional Reidemeister moves on oriented diagrams that were shown in [Pol10] to suffice for oriented links (appearing as the first four in Figure 23), the 14 possible Reidemeister moves involving singularities for oriented links can be reduced to the three depicted in Figure 23. Thus seven Reidemeister moves suffice for equivalency of singular diagrams.

Ia
Ib II

III IV
Va Vb FIGURE 23. A generating set of Reidemeister moves for oriented singular knots.
Inserting our labels as in Figure 15 (but with both strands pointing downward) into these three possibilities, we obtain the a set of axioms to go with our three traditional quandle axioms coming from the non-singular moves.
Definition 6.1. Let (X, ) be a quandle. Then if R 1 and R 2 are two maps from X × X to X satisfying the following relations, we say that (X, ) is an oriented singquandle.
Note that for the oriented singquandle, there are no axioms coming from successive rotations by 90 degrees of Figure 15. The authors of [BEHY18] give the following two examples of singquandles.
Example 6.2. Let n be a positive integer, and let a be an invertible element in Z n and b any element in Z n . Then the binary operations x y = ax + (1 − a)y, x −1 y = a −1 x + (1 − a −1 )y, R 1 (x, y) = bx + (1 − b)y and a)]y make the triple (Z n , , R 1 , R 2 ) satisfy the conditions to be an oriented singquandle. Example 6.3. Let X = G be a non-abelian group with the binary operation x y = y −1 xy. Then, for n ≥ 1, the following families of maps R 1 and R 2 make (X, , R 1 , R 2 ) into an oriented singquandle: (1) R 1 (x, y) = x(xy −1 ) n and R 2 (x, y) = y(x −1 y) n , (2) R 1 (x, y) = (xy −1 ) n x and R 2 (x, y) = (x −1 y) n y, In the case of proteins, we would like to consider bonds rather than singularities. There are two distinct types of oriented bonds, one where the orientations on the two strands are parallel and one where they are antiparallel, as in Figure 24.
Each of the fourteen moves involving singularities from [BEHY18] yields two possibilities corresponding to whether the singularity is replaced with a vertical or horizontal bond. However, it is still true that for each of a vertical or horizontal bond, the fourteen moves reduce to three. So in addition to the four non-singular Reidemeister moves, we have six more moves to consider.
x y R 1 (x,y) R 2 (x,y) x y R 3 (x,y) FIGURE 24. Labels at bonds with parallel and anti-parallel strands.
The first three correspond to bond diagram (A) in Figure 24, and we inherit the same set of relations as for the singquandle, namely (12)-(16).
Considering bond diagram (B) from Figure 24, we pick up two more functions R 3 (x, y) and R 4 (x, y). But note that rotation by 180 degrees switches the roles of x and y and the roles of R 3 and R 4 . Thus, it is always the case that R 4 (x, y) = R 3 (y, x). We will use this to eliminate R 4 (x, y) from all subsequent relations.
From Figure 25, we obtain four additional relations.
Definition 6.4. An oriented bondle is a quandle with operation £ and choices for functions R 1 (x, y), R 2 (x, y) and R 3 (x, y) such that they satisfy relations (12)-(16) and the additional relations: Note that the maps R 3 (x, y) = x and R 3 (x, y) = y do always satisfy the relations (17), (18), (19) and (20) for any quandle (X, ). We call these trivial solutions as they do not recognize the existence of the bond.
Since we already have examples of the desired maps R 1 and R 2 for both Example 6.2 and Example 6.3, we would like to find some solutions for the map R 3 satisfying relations (17), (18), (19) and (20) .
Lemma 6.5. Let n be a positive odd integer greater than or equal to 3 and let a be an invertible element of Z n . Consider the quandle (Z n , ) with x y = ax+(1−a)y and inverse operation x −1 y = a −1 x+(1−a −1 )y. Let m be an element in Z n and let R 3 be given by R 3 (x, y) = mx + (1 − m)y.
The following is a list of some (n, m) satisfying Corollary 6.6.
(1) If n = 15 then m = 6 or m = 10. Now we consider the case when the quandle is a group G with conjugation. First recall that the commutator of two elements x and y in a group G is given by [x, y] := xyx −1 y −1 . We have the following Lemma.
Lemma 6.7. Let X = G be a non-abelian group and let the quandle operation on G be given by x y = y −1 xy, so that x −1 y = yxy −1 . Assume that R 3 is given by R 3 (x, y) = x p y q , where p and q are integers, then (1) The map R 3 satisfies both equation (17) and equation (18) for any integers p and q.
(2) If for all x, y ∈ G, x p−1 y q = x −q y 1−p then R 3 satisfies equation Proof. Assume that R 3 has the form R 3 (x, y) = x p y q , then (1) One can see that equation (17) is satisfied for all integers p and q from the following.
Similarly, one has showing that equation (18) is satisfied also for all integers p and q.
(2) Now we check equation (19). Assume that the equation x p−1 y q = x −q y 1−p holds in G. Now we compute both the left hand side (LHS) and the right hand side (RHS) of equation (19).
Since x p−1 y q = x −q y 1−p , then LHS = RHS giving the result.
(3) We finish by checking equation (20). Here also we compute separately the LHS and the RHS. We thus have Now since the commutator [x p , y 1−p ] = 1, then we have x p y 1−p = y 1−p x p . Multiplying this equation by y p from the left and by y q−1 from the right gives the equation y p x p y q−p = yx p y q−1 , thus we have RHS = LHS giving equation (20).
In order to give a more explicit example of a non-abelian group with a map R 3 satisfying equations (19) and (20), we use the group of symmetries of a square.
Definition 6.8. Given a square with vertices labeled by 1, 2, 3 and 4, let G be the set of all rigid motions of the square that send vertices to vertices. Under composition, this set forms a non-abelian group called the dihedral group of order 8 and denoted D 4 . Precisely, D 4 = {1, r, r 2 , r 3 , s, sr, sr 2 , sr 3 }, where the permutation r = (1 2 3 4) is the clockwise rotation of 90 degrees and s is the reflection s = (1 2)(3 4).
Recall that in D 4 , the elements r and s satisfy the relations r 4 = 1 = s 2 and srs = r −1 . By iterating this last identity, we obtain sr i s = r −i , for 0 ≤ i ≤ 3. The square of any element of D 4 is either the identity element 1 or r 2 . Then we see that r 2 commutes with any other element of D 4 , since sr i r 2 = (sr i+2 s) s = r −i−2 s = r 4−i−2 s = r 2−i s = r 2 (sr i s)s = r 2 sr i . Corollary 6.9. In the dihedral group D 4 , the maps R 3 (x, y) = x 2 y −1 and R 3 (x, y) = x −1 y 2 both satisfy equations (19) and (20).
We thus obtain the following family of bondles.
Note that this example still holds if we change R 3 (x, y) = x 2 y −1 to R 3 (x, y) = x −1 y 2 .

EXAMPLES
Given two projections of proteins and a choice of bondle, we can count the number of distinct colorings of each projection by that bondle, and if those numbers are distinct, we know the two proteins are not topologically equivalent. This provides an opportunity for the categorization of proteins into distinct topological types. In the following we give two examples demonstrating the use of oriented bondles to topologically distinguish proteins with bonds.
Example 7.1. In this example, we use the oriented bondle (Z 15 , , R 1 , R 2 , R 3 ) from Corollary 6.6 with b = 2 to distinguish the topological type of the following two two proteins P 1 and P 2 .
Precisely, x y = 8(x + y), x −1 y = 2x − y, R 1 (x, y) = 2x − y and R 2 (x, y) = 7x − 6y. Note that we do not need to define R 3 (x, y) because there are no anti-parallel bonds in the diagrams.
A coloring of P 1 gives the following equation R 2 (y, x) = y which simplifies to 6(y − x) = 0 So if 5 divides y−x, we obtain a nontrivial coloring. Thius, the total number of colorings, including the trivial colorings is 45.
On the other hand, a coloring of P 2 gives the following equation R 2 (y, x) −1 R 1 (y, x) = y which simplifies to 11(y − x) = 0. Since 11 is invertible in Z 15 , we see that x = y, implying that P 2 has only trivial colorings, of which there are 15. Thus P 1 and P 2 are distinct.
But at the crossing with a dot in P 2 , we obtain the relation: This yields 0 = 7(y − x), and as 7 is invertible, we only obtain the 15 trivial colorings corresponding to y = x. Thus, the two proteins must be topologically distinct.

ACKNOWLEDGEMENT
Thanks to Jack Roche for suggesting the term "bondle".