Email updates

Keep up to date with the latest news and content from BMC Bioinformatics and BioMed Central.

Open Access Software

EGNAS: an exhaustive DNA sequence design algorithm

Alfred Kick12, Martin Bönsch1 and Michael Mertig12*

Author affiliations

1 Professur für Physikalische Chemie, Mess- und Sensortechnik, Technische Universität Dresden, Dresden 01062, Germany

2 Kurt-Schwabe-Institut für Mess- und Sensortechnik e.V. Meinsberg, Kurt-Schwabe-Straße 4, Ziegra-Knobelsdorf 04720, Germany

For all author emails, please log on.

Citation and License

BMC Bioinformatics 2012, 13:138  doi:10.1186/1471-2105-13-138

The electronic version of this article is the complete one and can be found online at:

Received:3 February 2012
Accepted:7 June 2012
Published:20 June 2012

© 2012 Kick et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.



The molecular recognition based on the complementary base pairing of deoxyribonucleic acid (DNA) is the fundamental principle in the fields of genetics, DNA nanotechnology and DNA computing. We present an exhaustive DNA sequence design algorithm that allows to generate sets containing a maximum number of sequences with defined properties. EGNAS (Exhaustive Generation of Nucleic Acid Sequences) offers the possibility of controlling both interstrand and intrastrand properties. The guanine-cytosine content can be adjusted. Sequences can be forced to start and end with guanine or cytosine. This option reduces the risk of “fraying” of DNA strands. It is possible to limit cross hybridizations of a defined length, and to adjust the uniqueness of sequences. Self-complementarity and hairpin structures of certain length can be avoided. Sequences and subsequences can optionally be forbidden. Furthermore, sequences can be designed to have minimum interactions with predefined strands and neighboring sequences.


The algorithm is realized in a C++ program. TAG sequences can be generated and combined with primers for single-base extension reactions, which were described for multiplexed genotyping of single nucleotide polymorphisms. Thereby, possible foldback through intrastrand interaction of TAG-primer pairs can be limited. The design of sequences for specific attachment of molecular constructs to DNA origami is presented.


We developed a new software tool called EGNAS for the design of unique nucleic acid sequences. The presented exhaustive algorithm allows to generate greater sets of sequences than with previous software and equal constraints. EGNAS is freely available for noncommercial use at webcite.

DNA sequence design algorithm; Hairpin; stem-and-loop structure; Single nucleotide polymorphism (SNP); Single-base extension (SBE); Polymerase chain reaction (PCR); DNA computing; DNA origami


Deoxyribonucleic acid (DNA) has the remarkable ability of specific molecular recognition with regard to its sequence of the nucleic bases adenine (A), thymine (T), guanine (G) and cytosine (C). This sequence recognition is based on the WATSON-CRICK base pairing of complementary bases A-T and G-C [1].

However, in mixtures of many different DNA strands in solution or on surfaces a correct hybridization is crucial for most applications referring to genetics, DNA nanotechnology, DNA origami, and DNA computing [2-8]. Cross hybridizations have to be minimized by controlling the uniqueness of all possible subsequence motifs in the set of used sequences. Additionally, specific care has to be taken with regard to secondary structures that can occur by folding due to intrastrand interactions. Hairpin structures reduce the hybridization efficiency, the binding rates, and thus, the detection limits on DNA microarrays [9].

Here we describe the exhaustive DNA sequence design algorithm EGNAS (Exhaustive Generation of Nucleic Acid Sequences). This algorithm is realized in a C++ program and is used to generate sequences with controlled intra- and interstrand properties. EGNAS is compared with previous tools. Data from selected publications are reanalyzed by applying EGNAS to proof the viability of this new algorithm.


We realized the sequence design algorithm EGNAS in a program written in C++. It is currently a command line program that was compiled by GNU COMPILER COLLECTION for LINUX, MAC OS X and MICROSOFT WINDOWS operating systems. EGNAS is freely available for noncommercial use at webcite. A manual and the first version of EGNAS are attached as additional files [see Additional files 1, 2, 3 and 4].

Additional file 1. Manual for the use of EGNAS. This DOC file contains the manual for the use of EGNAS. The sequence design criteria and options are explained [27,28].

Format: DOC Size: 20KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Additional file 2. Executable file for LINUX operating systems. This ZIP file contains the executable file of EGNAS for LINUX operating systems.

Format: ZIP Size: 51KB Download fileOpen Data

Additional file 3. Executable file for MAC OS X operating systems. This ZIP file contains the executable file of EGNAS for MAC OS X operating systems.

Format: ZIP Size: 46KB Download fileOpen Data

Additional file 4. Executable file for MICROSOFT WINDOWS operating systems. This ZIP file contains the executable file of EGNAS for MICROSOFT WINDOWS operating systems.

Format: ZIP Size: 196KB Download fileOpen Data

For comparison of computing time, sequence generations were performed on one and the same computer system with INTELCORETM i5 CPU 3.20 GHz and 4 GB RAM. The operating system was a 32-bit UBUNTU 10.04 (LINUX). The intra- and interstrand properties were verified by the analysis option of the software SEED[10,11].

Results and discussion

Previous software and algorithms

Numerous strategies for the DNA sequence design are described in literature [2-4,10-25]. We do not intend to analyze all algorithms, but we evaluate our results directly by comparing with previously published data. BRENNEMAN and CODON gave a brief overview to the topic of strand design [26]. To our best knowledge, so far, there is no software tool available offering the possibility to design oligonucleotide sequences with adequate consideration of intra- and interstrand interactions. Furthermore, a maximum set size of generated sequences is desirable for microarrays or DNA strands used as TAGs for addressing a high number of different targets.

Sets of sequences with defined properties can be achieved with the program SEED developed by SEIFFERT et al. [10,11], as well as with the software tools DNASEQUENCEGENERATOR and CANADA by FELDKAMP et al. [15,20,22,25]. These programs are freely available, work efficiently, provide satisfactory set sizes and meet the criteria of uniqueness among all sequences of a set concerning interstrand properties. However, it is difficult to achieve a big set size with these tools and concurrently to control intrastrand properties. It is hardly possible to avoid at the same time self-complementary subsequences, stem-loop structures and repetitions within one and the same strand, as these intrastrand properties can not be defined separately from the interstrand properties.

Both SEIFFERT and FELDKAMP use the criton concept introduced by SEEMAN[2,3]. Critons are all Lc bases long overlapping parts of a sequence strand. Their complements are termed anti-critons. Based on the criton rules, strands are generated consisting of unique basic sequences. There are <a onClick="popup('','MathML',630,470);return false;" target="_blank" href="">View MathML</a> critons and anti-critons. This fact limits the maximum possible number of valid sequences Ns with the length of Ls bases. Every sequence consists of (Ls Lc + 1) overlapping segments. To obey the criton rules, every basic sequence is used only once in the whole set and its complement is not allowed. If Lc is even, one half of all possible sequences is complementary to the other half, and Ns can be estimated using

<a onClick="popup('','MathML',630,470);return false;" target="_blank" href="">View MathML</a>


If Lc is odd, the maximum set size is estimated applying Equation (2), because there are <a onClick="popup('','MathML',630,470);return false;" target="_blank" href="">View MathML</a> self-complementary basic sequences that are not allowed.

<a onClick="popup('','MathML',630,470);return false;" target="_blank" href="">View MathML</a>


Further limitations result from restrictions with respect to the guanine-cytosine content (GC content) and forbidden sequence motifs defined by the user.

thumbnailFigure 1. Sliding. Example for the sliding of a complementary sequence pair due to the repetition of 7 bases long subsequences.

A novel sequence design algorithm

Sequence design criteria and options

The sequence design algorithm EGNAS offers the user different options. Consequently, the generated sequences meet certain criteria:

1. Sequence length Ls.

2. Length of basic sequences (criton length) Lc.

3. Exact GC content or its range.

4. No terminal adenine or thymine in the strand./The demand on “GC ends”.

5. Forbidden sequences./Included sequences.

6. Length of forbidden self-complementary subsequences Lsc.

7. Forbidden stem length of hairpin structures Lhp.

8. Length of subsequences that are not allowed to be repeated within one and the same sequence (“sliding”) Lsl.

9. Forbidden length of subsequences that could interact with complementary neighboring sequences Lni.

The EGNAS software provides the option to calculate the molar free enthalpy of DNA duplex formation. This calculation is based on the nearest-neighbor model [27] with parameters taken from SANTALUCIA et al. [28]. MARKY et al. investigated the helix-to-coil transition and described the “fraying” of a DNA double strand at the terminal T-A base pairs [29]. In the nearest-neighbor model, SANTALUCIA et al. assigned a penalty of 0.4 kcal/mol for strands with a thymine base at the 5’ end [28]. Thus, with EGNAS it is possible to generate sequences that do not contain terminal adenine or thymine. The risk of fraying in the designed DNA strands can be reduced by only allowing guanine or cytosine at the terminal strand positions. We refer to this option as the demand on “GC ends”.

Guanine-rich motifs in DNA can form parallel four-stranded complexes [30]. Furthermore, it is known that telomeric ends of eukaryotic chromosomes contain guanine-rich overhangs and form intra- and intermolecular structures [31]. Therefore, the subsequences “GGG” or “CCC” are often forbidden in DNA strand design to circumvent the formation of guanine tetrads between hairpins.

thumbnailFigure 2. Flowchart of the exhaustive design algorithm. Simplified description of the exhaustive design algorithm of EGNAS. The input of the design criteria and options is used to generate basic sequences and the filter. The filter extracts valid basic sequences before the set generation. Additionally, strings are filtered within the sequence generation.

Homopolymeric runs of adenine or thymine are unwanted for certain applications requiring DNA polymerization. This process depends on the fidelity of DNA polymerases, which is influenced by the strand sequence that should be replicated. The repeat of one and the same base in the sequence can lead to error-prone replication through template-primer slippage [32]. Therefore, adenine or thymine runs, starting from four repeats, are preferentially forbidden for the sequence design. Forbidden sequences or subsequence motifs are specified by listing each of them in a set that is denoted further on in curly brackets. For instance, if both subsequences “GGG” and “CCC” are not allowed, the set containing the forbidden elements will be {GGG;CCC}.

Intrastrand properties

EGNAS offers novel options concerning the intrastrand properties. These are essential prerequisites to avoid secondary structures due to self-complementary sequences and hairpin formation. In contrast to the criton concept [2,3], as applied by SEIFFERT et al. [10,11] and FELDKAMP et al. [15,20,22], the novel algorithm treats intrastrand properties separately. Below, we define the used terms and describe these properties.

Hairpins and self-complementarity

Hairpin structures are also called stem-and-loop structures. They consist of two complementary arm sequences and the loop sequence. The arm sequences are able to form the double-stranded stem while being connected by the single-stranded loop sequence. Self-complementarity is treated as a special case of a hairpin when the loop size is zero. Therefore, a self-complementary sequence has always an odd number of bases. If hairpin structures with a Lhp bases long stem are forbidden, self-complementary subsequences that are equal to or longer than 2Lhp bases will consequentially be omitted.


We use the term “sliding” for an intrastrand property of a sequence. Sliding denotes that a subsequence can be found several times at different positions of one and the same strand. If a complementary strand hybridizes with such a strand, different positions will be possible. Thus, sliding between the hybridizing strands would take place. For example, the sequence pair in Figure 1 allows sliding due to the repetition of 7 bases long subsequences (Lsl = 8).

Interaction with the neighboring sequences

If primers are paired with TAGs, a special criterion for the TAG sequences will arise. In this case, primer foldback can become a problem. For example, the formation of hairpins will cause signals in single-base extension (SBE) reactions even lacking a template.

In the literature different approaches are given for the primer-TAG pairing. Those methods work with sets of previously found TAGs. HIRSCHHORN et al. suggested the calculation of an empirical foldback score [7]. Accordingly, if a foldback score is greater than a threshold value, the SBE primer will be paired with another TAG. KADERALI et al. used a free energy alignment algorithm for primer-TAG pairing [19].

We limit the interaction with neighboring strands already during the sequence generation. Therefore, all possible Lni bases long complementary subsequences of a neighboring strand are forbidden for the generation of the corresponding TAG. This is especially intended for designing strands where a molecular spacer is located between the neighboring sequence and the TAG. For example, such a spacer could be a hexaethylene glycol moiety [9]. Actually, overlapping subsequences, which would evolve through directly attaching one neighboring sequence to either the 3’ or 5’ end of a TAG, are not considered. Nevertheless, in Section “Combination of TAGs with primers” we show that even if TAGs are attached to neighboring sequences without an intended spacer, hairpin formation will still be diminished significantly.

Description of the exhaustive design algorithm

In the following section we describe, how EGNAS provides a set of valid sequences. A simplified flowchart of the underlying algorithm is given in Figure 2. Initially, all Lc bases long subsequences in the included and neighboring sequences are read out and saved together with their complements as forbidden sequences. The included and neighboring sequences are predefined by the user, meaning that Lc bases long cross-hybridizations are not allowed to occur between one of them and any of the generated sequences. Subsequently, all possible basic sequences with the length Lc are generated in such a way that they obey the criteria specified by the user. Thereafter, a string is formed stepwise from left to right by randomly combining allowed basic sequences. If the sequence length is not divisible by the criton length, additional basic sequences of the remainder length will be generated and used for the termination of the sequence design. After each step, the growing sequence is checked against the criteria chosen by the user.

thumbnailFigure 3. Influence of the GC content and GC ends on the set size. Dependence of the set sizes on the GC content. Sets of 10 bases long sequences with Lc = 4 for global criton rules were generated. Averages ± standard deviations were calculated from 1,000 sets for restricted and from 10,000 sets for unrestricted conditions [0; 100]. The maximum set size N(max) = 17 is shown.

After trying all of the valid basic sequences as suffixes, there can still be combinations of basic sequences that do not allow appending of any remaining basic sequence, because the forming strings do not meet at least one of the necessary criteria. In such a case, this combination is forbidden to occur at the beginning of a string. Thus, a new trail will skip the basic sequences leading to the forbidden ones.

A sequence will be saved, if it has the defined length and meets all criteria. Every Lc bases long subsequence and its complement within this generated strand are forbidden and removed from the list of basic sequences. Then, the next sequence generation starts by the analog combination of the remaining basic sequences.

If all basic sequences are either forbidden or used, the sequence generation of the current set will be finished. Furthermore, the user can force every set to be complete. Namely, no additional sequences exist with the same user-defined restrictions of the set. The generation of a complete set is at the expense of computing time. Hence, limiting the maximum number of generated sequences can reduce the computing time significantly.

Software comparison and performance tests

Global criton rules and intrastrand properties

In order to compare the novel algorithm with other available software tools, it was necessary to distinguish between two different settings of options for the sequence generation. On the one hand, the preassigned criton rules in Section “Intrastrand properties” are valid for the whole set of generated sequences. This holds for SEED[10,11], CANADA[25] and DNASEQUENCEGENERATOR[15,20]. In this cases, we refer to “global criton rules”. It means for the intrastrand properties (Lsc, Lsl and Lhpthat they are treated equally to interstrand criteria. In detail, if the criton length is Lc, Lsc bases long self-complementary subsequences will be forbidden, as well as sliding of Lsl bases long subsequences. Hairpin structures with a stem lengths of Lhp bases and longer will not be allowed. Thereby, the following relations hold:

1. Lc = Lsl = Lhp

2. If Lc is odd: Lsc = Lc

3. If Lc is even: Lsc = Lc + 1

On the other hand, EGNAS offers a control of intrastrand properties independent of the criton rules. In this case, the criton rules only control the cross hybridization among different sequences not taking into account self-complementarity, hairpins or sliding.

Influence of the GC content and GC ends on the set size

The influence of the GC content and GC ends were investigated by generating sets of sequences with global criton rules, Lc = 4. Those sets were generated for every possible GC content of 10 bases long sequences.

The restriction of the GC content can lead to a reduced set size as shown in Figure 3. While 50% GC content is not necessarily a restriction to the set size, extreme values of the GC content limit the set size significantly. At GC contents of 30% and 70%, the set size is almost half compared with the set sizes achieved without restrictions or with 50% GC content.

The demand on GC ends in combination with the restriction of the GC content to exactly 50% lowers the set size significantly. Otherwise, no evident reduction of the set size is observed. With Equation (2), the maximum number of sequences is 17. With 16 sequences a yield of 94% could be achieved where no restrictions were set or only the GC content was forced to be 50%. These are high set sizes in relation to other calculations performed by FELDKAMP[22] and those that are presented in Table 1.

Table 1. Influence of the criton and sequence length on the set size

Influence of the criton and sequence length on the set size

Table 1 shows the dependence of the set size from the combination of sequence and criton length. For each of the 124 combinations, 10 sets were calculated with Lc for global criton rules. In accordance with Equations (1) and (2), the set sizes grow with higher criton length and lower sequence length. Of course, the criton length has the highest influence, as it is in the exponent of those equations. FELDKAMP presented an analog table [22].

In average, FELDKAMP generated 84.4 ± 4.3% of the maximum possible numbers of sequences per set. Our new approach has a slightly higher average yield of 87.0 ± 4.6%. In 120 of the 124 different combinations of criton and sequence lengths, examined in Table 1, set sizes equal to or higher than those of FELDKAMP’s tool were calculated with EGNAS.

Variation of the set size

Another interesting issue of sequence sets is, how their sizes will be distributed, if the same options are chosen for multiple sequence generations. Figure 4 compares the results of FELDKAMP’s [22] and of our software. The combination of Lc = 6 and 20 bases long sequences was chosen for global criton rules. FELDKAMP generated 100 sets with these settings. However, we evaluated 10,000 sets to gain a higher confidence level. The most of FELDKAMP’s sets have 112 sequences, whereas the most of our completed sets contain 114 sequences. In both cases, these set sizes appear with a relative frequency of about 30%.

thumbnailFigure 4. Variation of the set size. Histogram of the relative frequency of the set sizes. The results of EGNAS and FELDKAMP[22] are compared with each other 10,000 sets of 20 bases long sequences with Lc = 6 for a global criton rule were generated.

Sliding and its influence on the set size

Within one and same strand, EGNAS can control the repetition of subsequences of a certain length. As an example, we generated sets of sequences with Lc = 6 for global criton rules but with varying Lsl. The results of this calculations are shown in Figure 5. For Lsl = 6, the global criton rules are valid, and in average, 113.8 ± 1.1 sequences are generated. Significant reduction of the set size to 100.4 ± 1.1 sequences is observable for Lsl = 3. This is because the decrease of Lsl reduces the number of usable basic sequences. A considerable increase of the set size to 142.7 ± 3.5 sequences is obtained by setting Lsl = 13. In this case, Lsl is large enough to allow the repetition of one basic sequence. Thus, for a complete strand, less basic sequences are used, and consequently, more usable basic sequences are left for the generation of further strands.

thumbnailFigure 5. Sliding and its influence on the set size. Dependence of the set size on the sliding length (Lsl). The sequence length is 20 and Lc = 6 was used for global criton rules but with varying Lsl. Averages ± standard deviation of 10 sets were calculated for each Lsl.

Controlling intrastrand properties

Self-complementarity and hairpins are the major reasons for the formation of secondary structures. Such structure formation competes against the hybridization with target DNA strands. Hence, a fast and sensitive detection on microarrays is hindered by secondary structures [9]. Accordingly, we generated sequences not only with minimum criton length Lc for a global criton rule but we also prevented self-complementarity and the risk of hairpin formation. The DINAMELT WEB SERVER was used to estimate the stability of secondary structures [33] ( webcite (February 1, 2012)). The stabilities were calculated as the molar free enthalpies (ΔG) of the most stable hairpins at 37°C and 1 mol/l sodium ions.

FELDKAMP et al. had generated 22 bases long sequences with the uniqueness Lc = 5 and chose sequences with the most unstable secondary structures afterwards [34]. In Table 2 the generated sequences of FELDKAMP et al. and our program are shown to compare the stabilities of possible secondary structures. Hairpins could be avoided by choosing Lhp = 2. Even so, enough sequences could be generated to compete with FELDKAMP’s sequence set. FELDKAMP’s software tool is not able to treat intrastrand properties separately. For instance, a stem length of 4 bases can be found in sequence No. 11 (5'‐CAAGgtctgCTTG.....‐3').

Table 2. Avoiding hairpin formation in 22 bases long sequences

The restriction of Lhp will improve sequence quality significantly, if no strand folding is desired. However, smaller values for Lhp result in smaller set sizes. Calculations of 10 sets were performed with varying Lhp and Lc = Lsc = Lsl = 6 for 20 bases long sequences. A clear reduction of the set size can be found from 111.6 ± 1.3 (Lhp = 3) to 61.8 ± 1.2 (Lhp = 2). For Lhp = 6, the set size is 113.8 ± 1.5 and not greater than the set sizes with Lhp = 4 or 5.

Further examples for preventing hairpins are given in Tables 3, 4 and 5. TANAKA et al. [17] and FELDKAMP et al. [20] published sets of 20 bases long sequences that we compare with our results with respect to the stabilities of hairpin structures (Table 3). Sequences with lower stabilities of hairpins were designed by applying EGNAS. Additionally, the set size could be increased from 14 to 16 sequences.

Table 3. Avoiding hairpin formation in 20 bases long sequences

Table 4. Calculations with included sequences to design additional 20 bases long sequences

Table 5. Calculations with included sequences to design additional 15 bases long sequences

Calculations with included sequences

Another useful feature of a sequence design algorithm is the possibility to include already existing sequences into the calculations. As a result, the cross hybridization of newly formed sequences with the included strands is controllable. The EGNAS user may include naturally occurring sequences of genomic DNA, plasmids, cloning vectors, primers or previously designed sets of sequences with certain properties. For instance, we performed stepwise generation of sequences shown in Table 4. In the first step nine 15 bases long sequences were found. These sequences were included in the second run to design four additional 20 bases long sequences. With EGNAS, one more 20 bases long sequence was generated in comparison to sequences suggested by ARITA et al. [13] and FELDKAMP et al. [20]. This was possible in spite of extra restrictions according to GC ends and secondary structure. Possible hairpin structures could be completely eliminated in these computations.

The optimization of sequences is also possible with a different stepwise strategy. Here we first start with strong restrictions on sequence properties and weaken them until the desired set size is achieved.

Accordingly, the maximum stem length of hairpins were first limited to one base pair by setting Lhp = 2. A set size of maximum 14 sequences was possible with this restriction. Finally, in the next sequence generation hairpins with stem length of 2 base pairs were allowed by setting Lhp = 3 and guanine was forbidden. Through higher Lhp, more variability of sequences is possible. Nevertheless, by forbidding guanine, there are only stem structures possible consisting of two consecutive adenine/thymine base pairs. These structures are less stable than those comprising guanine/cytosine base pairs. Therefore, other authors presented design algorithms that only use {ATC} for the sequence generation, and this aims to minimize the risk of hairpin formation [14,35]. In practice, this concept was realized for the TAG/anti-TAG system offered by LUMINEX[36]. For our calculations, the 14 previously generated sequences were included in order to expand the sequence set. Consequently, six further sequences could be generated. The result of this approach is shown in Table 5. It depicts a comparison between sets of 15 bases long sequences generated by FAULHAMMER et al. [14], FELDKAMP et al. [20] and EGNAS. EGNAS offers better sequences with respect to secondary structures, as the minimum molar free enthalpies of the most stable hairpin structures are -0.7, 0.1 and 1.8 kcal/mol for FELDKAMP et al., FAULHAMMER et al. and EGNAS, respectively. Comparing the criton length Lc = 5 and the GC content of 40%, there is no improvement versus FELDKAMP’s sequences. But a significant reduction of the secondary structure stabilities is evident, since the forbidden stem length was reduced from Lhp = 5 to Lhp = 3. Moreover, in contrast to FELDKAMP’s sequences, only GC ends were allowed. In spite of this further restriction, an equal set size was achieved.

Balancing cross hybridization, sliding and hairpins

We investigated further advantages of discriminating between global criton rules and intrastrand properties. The results are summarized in Table 6. The generated sequences are compared with those of SHIN et al. [16] and of FELDKAMP et al. [20]. There are only seven 20 bases long sequences in one set. We show that the user of EGNAS can balance intra- and interstrand properties to achieve a required set size.

Table 6. Balancing cross hybridization, sliding and hairpins

Because with Lc = 5, Lhp = 2 and Lsl = 5, the set sizes clearly extended seven sequences, sliding could be decreased to Lsl = 3 and even only GC ends were allowed for the first set. For the second set of sequences, we raised Lsl from 3 to 6, Lhp from 2 to 3 and lowered the criton length Lc from 5 to 4. Consequently, we have to accept sliding of 5 bases long subsequences. Additionally, adenine or thymine had to be allowed as terminal bases to generate seven sequences. Indeed, this is against the global criton rule with Lc = 5. However, no cross hybridization between different strands would take place with 4 bases long subsequences.

In both of our sets, the formation of hairpins is less probable then in the sets published by SHIN et al. [16] and by FELDKAMP et al. [20].

Computing time

In this section we demonstrate that the sequence generation lasts only a few seconds, if the user limits the set size reasonably. To this aim, 10 sets of 20 bases long sequences were generated for each set size with Lc = 6 for a global criton rule. We varied the maximum number of sequences to be generated starting from 100. The highest defined set size was 109 sequences, since we know from calculations concerning the set size distribution that this is the minimum set size that could be generated for a complete set.

After the defined set size had been achieved, the sequence generation stopped and the computing time was saved. Insisting on the completeness of a set leads to considerable increase of the computing time from a few seconds (4 to 9 s) to approximately one minute (56.2 s ± 12.9 s). Nevertheless, this is still in reasonable computing time for a complete set size. Actually, the computing time grows considerably with Lc. The generation of a complete set with Lc ≥ 8 can last hours, as the number of basic sequences also grows exponentially with Lc(Equations (1) and (2)).

Combination of TAGs with primers

Genotyping of single nucleotide polymorphisms (SNPs) is one example for a particular application of SBE reactions (Figure 6). This technique, as described by FAN et al., requires SBE-TAG primers [8]. Marker-specific polymerase chain reaction (PCR) primers are needed for the amplification of SNPcontaining regions. Thereafter, the PCR products are used as templates for the SBE reaction with SBE-TAG primers. The 3’ ends of these primers are complementary to the specifically flanking regions of the SNP loci and the 5’ ends are complementary to specific probes on an array. The generation of suitable SBE-TAGs was performed with the EGNAS algorithm. In this case, the TAG sequences were generated with respect to minimal interaction with the neighboring SNP flanking sequence. Consequently, no 3 bases long subsequences were complementary between one SNP flanking sequence and its corresponding TAG sequence (Lni = 3). The primers and SNP flanking sequences were taken from FAN et al. [8] to generate 148 TAGs that were 20 bases long with Lc = 8, Lsl = 8, Lhp = 2 and GC content 50%. The subsequences {GGG;CCC;AAAAA;TTTTT} were forbidden and GC ends were demanded. Forward and reverse primers as well as the flanking regions of the SNPs were chosen as included sequences. Thus, no cross hybridizations with 8 bases long subsequences occur among the TAGs and included sequences, secondary structures were avoided within the TAGs and foldback due to interactions with SNP flanking sequences were minimized.

thumbnailFigure 6. Principle of a TAG-SBE genotyping assay. Principle of a TAG-SBE genotyping assay to detect an A/G SNP. TAG sequences were designed and paired with primers by using EGNAS.

We generated TAGs without stable hairpins (Lhp = 2, ΔG > 0.4 kcal/mol) by using EGNAS. In contrast, the TAGs presented by FAN et al. show one very stable hairpin with a stem length of 3 base pairs (Lhp = 4, ΔG =−2.7 kcal/mol). Moreover, with regard to foldbacks, the TAG-primer pairing with EGNAS is better (Lni = 3) than presented by FAN et al. (Lni = 7) for 148 SBE-TAG primers. The advantages of EGNAS become also apparent by comparing the uniqueness of the TAG sequences. FAN et al. used TAGs with Lc = 11 and EGNAS provided TAGs with Lc = 8. Thus, the quality of the TAGs can be evaluated by the values of Lhp, Lni and Lc, when hairpins, TAG-primer pairing and cross hybridization are in the focus, respectively. Corresponding results are shown in Tables 7 and 8. All TAGs and primer sequences as well as the molar free enthalpies of the most stable hairpins are given in additional files [see Additional file 5 and Additional file 6].

Additional file 5. SBE-TAGs and primer sequences. This MICROSOFT OFFICE EXCEL 2003 sheet contains SBE-TAGs and primer sequences as well as the molar free enthalpies of the most stable hairpins.

Format: XLS Size: 63KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Additional file 6. SNPs and PCR primer sequences. This MICROSOFT OFFICE EXCEL 2003 sheet contains SNPs and PCR primer sequences as well as the molar free enthalpies of the most stable hairpins.

Format: XLS Size: 34KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Table 7. Combination of TAGs with Primers — comparison of TAGs

TAGs for staple strands of DNA origami

Self-assembled structures are fundamental in the field of nanotechnology. DNA allows the assembly of programmable building blocks. One promising method to control and organize functional materials is based on DNA origami [37]. The underlying concept requires an accurate design of DNA sequences. A spatially precise functionalization is accessible through the unique addressability of DNA origami on the nanometer scale. The DNA origami scaffold is a long single strand of DNA, which is folded by hundreds of short synthetic oligonucleotides called staple strands. The staple strands are designed to bind several desired parts of the scaffold, and thus, to connect otherwise distant sites of this long single-stranded DNA. ROTHEMUND presented a variety of different shapes based on this concept [6].

There is software available for the design of 3D DNA origami shapes, for instance the CADNANO software [38]. However, we show that EGNAS will be useful to find proper sequences, if staple strands have to be extended by anti-TAGs that stay single-stranded during origami formation. These anti-TAGs are used as sticky ends or capture probes to bind the TAGs of DNA-modified nanoparticles or peptide-DNA conjugates to a given DNA origami structure. In the following, we examine two examples with single-stranded circular M13mp18 viral DNA as scaffold. Its sequence [39] was included in the calculations below.

Sticky ends for triangular DNA origami

In the work of DING et al., triangular DNA origami is used for the assembly of six gold nanoparticles through DNA hybridization (Figure 7) [40]. The particles have sizes of either 5, 10 or 15 nm. Each particle is captured by three probes on the DNA origami with the M13mp18 DNA as a scaffold. To this aim, 18 staple strands are modified with 24 bases long sticky ends (anti-TAGs). DING et al. designed four different TAG sequences to label the 10 and 15 nm gold particles each with one sequences, and to modify two differently labeled 5 nm particles.

Table 8. Combination of TAGs with Primers — TAG-primer pairing

thumbnailFigure 7. Principle of the arrangement of gold nanoparticles on a triangular DNA origami. Principle of the arrangement of gold nanoparticles on a triangular DNA origami by extending staple strands with sticky ends. The corresponding sequences of the sticky ends were designed and paired by using EGNAS.

Table 9. Modification of staple strands for DNA origami

We analyzed possible cross hybridization of the sticky ends with the staple strands and with the scaffold. As a result, the 9 bases long sequence 5’-GAATCCTGA-3’ is identical in the staple strand “C28” and in the TAG “S5a” of the 5 nm gold nanoparticles. This could result in a cross hybridization of three modified staple strands with the unmodified “C28” strand. There are also numerous 8 and 9 bases long sequences that could cause cross hybridizations with the scaffold strand. Additionally, we found possible hairpin structures with a 4 base pairs long stem (5’-AGTC-3’) in the sequences of the TAG “S5b” and the corresponding anti-TAG. Hairpin structures lower the hybridization efficiency of the TAGs and anti-TAGs. This could be one explanation for missing 5 nm particles on the DNA origami. Actually, DING et al. mentioned this problem when only two sticky ends per particle were used.

We tried to improve the sequences of the sticky ends by applying EGNAS. To this aim, we included the scaffold and all staple strand sequences in the calculations. More precisely, we chose Lc = 8, and thus, avoided cross hybridizations with more than 7 consecutive base pairs of the sticky ends with the modified and unmodified staple strands, as well as with the scaffold. Secondary structures were evaded by setting Lhp = 2. We generated the sticky ends stepwise and included the previously designed sequences of the modified staple strands for the design of the next sticky end. Furthermore, all 4 bases long sequences that were partly complementary to the neighboring staple strands were forbidden for the current generation of the corresponding sticky end (Lni = 4). The risk of foldbacks is minimized by that procedure.

Consequently, the sequences used by DING et al. to modify staple strands exhibit the possibility of more stable cross hybridizations and secondary structures compared with sequences that were generated with EGNAS (Table 9). The sequences of staple strands before and after modification as well as their molar free enthalpies of the most stable hairpins are given in an additional file [see Additional file 7].

Additional file 7. Sequences for triangular DNA origami. This MICROSOFT OFFICE EXCEL 2003 sheet contains the sticky end sequences for triangular DNA origami. The sequences of staple strands are shown before and after modification. Their molar free enthalpies of the most stable hairpins are presented.

Format: XLS Size: 15KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Capture probes for six-helix bundle nanotubes

STEARNS et al. tried to organize a peptide-DNA conjugate on a six-helix bundle forming DNA origami (Figure 8) [41]. The peptide A3 was used for in situ nucleation and growth of gold nanoparticles, as it recognizes gold surfaces and reduces soluble gold ions. Only one sequence for capture probes was used to modify 10 staple strands to fold the M13mp18 DNA as a scaffold. Again, we analyzed the cross hybridization of the probe sequence with the staple strands and with the scaffold. Only one 8 base pairs long possible cross hybridization with the sequence 5’-GCCGTTGA-3’ of the staple strand “70” with the peptide-DNA conjugate was found. There are possible hairpin structures with a 5 base pairs long stem (5’-CGTTG-3’) in the probe sequences. Furthermore, foldback structures are possible through 5 bases long complementary subsequences (5’-AACGG-3’) of the probe and the staple strand “122”.

thumbnailFigure 8. Principle of in situ nucleation and growth of gold nanoparicles on a six-helix bundle nanotube. Principle of in situ nucleation and growth of gold nanoparicles on a six-helix bundle nanotube with extended staple strands. The corresponding sequence of a capture probe was designed and paired by using EGNAS.

Here, we optimized the probe sequence with equal parameters like presented above for the sticky ends for triangular DNA origami. As a result, to avoid secondary structures, the sequences that were generated by using EGNAS are more suitable to capture the peptide-DNA conjugate compared with those used by STEARNS et al. to modify the staple strands (Table 9). Sequences of staple strands before and after modification as well as their molar free enthalpies of the most stable hairpins are given in an additional file [see Additional file 8]. STEARNS et al. reported that the particles were not organized in the designed 27 nm intervals. Hence, the peptide was probably immobilized not efficiently enough. Hairpin formation of the probe and foldback with the staple strands might be the reasons for the low density of nanoparticles on the DNA origami.

Additional file 8. Sequences for DNA origami forming six-helix bundle nanotubes. This MICROSOFT OFFICE EXCEL 2003 sheet contains the capture probe sequences for DNA origami forming six-helix bundle nanotubes. The sequences of staple strands are shown before and after modification. Their molar free enthalpies of the most stable hairpins are presented.

Format: XLS Size: 11KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data


We developed a new software tool called EGNAS for the design of unique nucleic acid sequences. Sets of sequences with defined intra- and interstrand properties can be generated in reasonable computing time. A maximum set size with given constraints can be achieved. The presented exhaustive algorithm allows to generate greater sets of sequences than with previous software and equal constraints.

In the present stage, the developed program is suitable for generating sequences for DNA-TAGs avoiding secondary structures and cross hybridizations. Furthermore, predefined sequences can be taken into consideration. This is in principle useful for applications where no interaction of TAGs with template strands is desired. The offered option of TAG-primer pairing with regard to minimal foldback facilitates the generation of TAG sequences for multiplexed genotyping of SNPs. This kind of genotyping can be performed with PCR and SBE reactions on microarrays or bead surfaces. Additionally, EGNAS affords the computer aided design of sequences for specific attachment of molecular constructs to DNA origami. In further development the novel algorithm could be optimized to efficiently include large gene sequences. An extension of EGNAS to design branched structures like nucleic acid junctions is possible.

Availability and requirements

Project name: EGNAS

Project home page: webcite

Operating systems:LINUX, MAC OS X, and MICROSOFT WINDOWS

Programming language: C++

Other requirements: None

License: Free for noncommercial use

Any restrictions to use by nonacademics: License needed


EGNAS: Exhaustive generation of nucleic acid sequences; DNA: Deoxyribonucleic acid; A: Adenine; T: Thymine; G: Guanine; C: Cytosine; GC: Content, Guanine-cytosine content; SBE: Single-base extension; SNP: Single nucleotide polymorphism; PCR: Polymerase chain reaction; ddUTP: Dideoxyuridine triphosphate; ddCTP: Dideoxycytosine triphosphate.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

AK conceived the algorithm, performed the calculations and drafted the manuscript. MB participated in developing the algorithm and wrote the software program in C++. MM helped to draft the manuscript and revised it critically. All authors read and approved the final manuscript.


This work was financially supported by the Deutsche Forschungsgemeinschaft (Contract No.: ME 1256/10-1) and the Sächsische Aufbaubank (Contract No.: 14120/2447).


  1. Watson JD, Crick FH: Molecular structure of nucleic acids: A structure for deoxyribose nucleic acid.

    Nature 1953, 171(4356):737-738. PubMed Abstract | Publisher Full Text OpenURL

  2. Seeman NC: Nucleic acid junctions and lattices.

    J Theor Biol 1982, 99(2):237-247. PubMed Abstract | Publisher Full Text OpenURL

  3. Seeman NC, Kallenbach NR: Design of immobile nucleic acid junctions.

    Biophys J 1983, 44(2):201-209. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Seeman NC: De novo design of sequences for nucleic acid structural engineering.

    J Biomol Struct Dyn 1990, 8(3):573-581. PubMed Abstract | Publisher Full Text OpenURL

  5. Adleman LM: Molecular computation of solutions to combinatorial problems.

    Science 1994, 266(5187):1021-1024. PubMed Abstract | Publisher Full Text OpenURL

  6. Rothemund PWK: Folding DNA to create nanoscale shapes and patterns.

    Nature 2006, 440(7082):297-302. PubMed Abstract | Publisher Full Text OpenURL

  7. Hirschhorn JN, Sklar P, Lindblad-Toh K, Lim YM, Ruiz-Gutierrez M, Bolk S, Langhorst B, Schaffner S, Winchester E, Lander ES: SBE-TAGS: An array-based method for efficient single-nucleotide polymorphism genotyping.

    Proc Natl Acad Sci USA 2000, 97(22):12164-12169. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. Fan JB, Chen X, Halushka MK, Berno A, Huang X, Ryder T, Lipshutz RJ, Lockhart DJ, Chakravarti A: Parallel Genotyping of human SNPs using generic high-density oligonucleotide tag arrays.

    Genome Res 2000, 10(6):853-860. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Kick A, Bönsch M, Katzschner B, Voigt J, Herr A, Brabetz W, Jung M, Sonntag F, Klotzbach U, Danz N, Howitz S, Mertig M: DNA microarrays for hybridization detection by surface plasmon resonance spectroscopy.

    Biosens Bioelectron 2010, 26(4):1543-1547. PubMed Abstract | Publisher Full Text OpenURL

  10. Seiffert J: Ein Sequenzdesign-Algorithmus für verzweigte DNA-Strukturen.

    PhD thesis.

    Technischen Universität Dresden, Fakultät Informatik 2008


  11. Seiffert J, Huhle A: A full-automatic sequence design algorithm for branched DNA structures.

    J Biomol Struct Dyn 2008, 25(5):453-466. PubMed Abstract | Publisher Full Text OpenURL

  12. Deaton R, Garzon M, Murphy RC, Rose JA, Franceschetti DR, Stevens JrSE: Genetic search of reliable encodings for DNA-based computation. In Late Breaking Papers at the Genetic Programming 1996 Conference Stanford University July 28-31, 1996. Edited by Koza JR. Stanford University, CA, USA: Stanford Bookstore; 1996:9-15. OpenURL

  13. Arita M, Nishikawa A, Hagiya M, Komiya K, Gouzu H, Sakamoto K: Improving sequence design for DNA computing. In Proceedings of Genetic and Evolutionary Computation Conference 2000 (GECCO ’00) July 8-12, 2000. Edited by Whitley LD, Goldberg DE, Cantú-Paz E, Spector L, Parmee IC, Vegas BeyerH-G. Las. Nevada, USA: Morgan Kaufmann; 2000:875-882. OpenURL

  14. Faulhammer D, Cukras AR, Lipton RJ, Landweber LF: Molecular computation: RNA solutions to chess problems.

    Proc Natl Acad Sci USA 2000, 97(4):1385-1389. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Feldkamp U, Saghafi S, Banzhaf W, Rauhe H: DNA Sequence Generator: A program for the construction of DNA sequences DNA computing. In DNA Computing Volume 2340 of Lecture Notes in Computer Science. Edited by Jonoska N, Seeman N. Berlin, Heidelberg: Springer Berlin / Heidelberg; 2002:23-32. OpenURL

  16. Shin SY, Kim DM, Lee IH, Zhang BT: Evolutionary sequence generation for reliable DNA computing. In Proceedings of the 2002 Congress on Evolutionary Computation (CEC ’02) May 12-17 2002, Volume 1. Hilton, Hawaiian Village Hotel, Honolulu, Hawaii, USA; 2002:79-84. OpenURL

  17. Tanaka F, Nakatsugawa M, Yamamoto M, Shiba T, Ohuchi A: Towards a general-purpose sequence design system in DNA computing, Volume 1. In Proceedings of the 2002 Congress on Evolutionary Computation (CEC ’02) May 12-17, 2002. Hilton, Hawaiian Village Hotel, Honolulu, Hawaii, USA; 2002:73-78. OpenURL

  18. Tanaka F, Nakatsugawa M, Yamamoto M, Shiba T, Ohuchi1 A: Developing support system for sequence design in DNA computing. In DNA Computing Volume 2340 of Lecture Notes in Computer Science. Edited by Jonoska N, Seeman N. Berlin, Heidelberg: Springer Berlin / Heidelberg; 2002:129-137. OpenURL

  19. Kaderali L, Deshpande A, Nolan JP, White PS: Primer-design for multiplexed genotyping.

    Nucleic Acids Res 2003, 31(6):1796-1802. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Feldkamp U, Rauhe H, Banzhaf W: Software tools for DNA sequence design.

    Genet Programming Evolvable Machines 2003, 4(2):153-171. Publisher Full Text OpenURL

  21. Yin P, Guo B, Belmore C, Palmeri W, Winfree E, LaBean TH, Reif JH: TileSoft: Sequence optimization software for designing DNA secondary structures.


    (February 1, 2012). [ http://www.cs.duke. edu/reif/paper/peng/TileSoft/TileSoft.pdf webcite]

  22. Feldkamp U: Computer aided DNA sequence design.

    PhD thesis.

    Universität Dortmund, Fachbereich Informatik 2005


  23. Feldkamp U, Niemeyer CM: Rational design of DNA nanoarchitectures.

    Angew Chem Int Ed 2006, 45(12):1856-1876. Publisher Full Text OpenURL

  24. Zhu J, Wei B, Yuan Y, Mi Y: UNIQUIMER 3D, a software system for structural DNA nanotechnology design, analysis and evaluation.

    Nucleic Acids Res 2009, 37(7):2164-2175. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Feldkamp U: CANADA: Designing nucleic acid sequences for nanobiotechnology applications.

    J Comput Chem 2010, 31(3):660-663. PubMed Abstract | Publisher Full Text OpenURL

  26. Brenneman A, Condon A: Strand design for biomolecular computation.

    Theor Comput Sci 2002, 287:39-58. Publisher Full Text OpenURL

  27. Breslauer KJ, Frank R, Blöcker H, Marky LA: Predicting DNA duplex stability from the base sequence.

    Proc Nat Acad Sci 1986, 83(11):3746-3750. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. SantaLucia JJr, Allawi HT, Seneviratne PA: Improved nearest-neighbor parameters for predicting DNA duplex stability.

    Biochemistry 1996, 35(11):3555-3562. PubMed Abstract | Publisher Full Text OpenURL

  29. Marky LA, Canuel L, Jones RA, Breslauer KJ: Calorimetric and spectroscopic investigation of the helix-to-coil transition of the self-complementary deoxyribonucleotide ATGCAT.

    Biophys Chem 1981, 13(2):141-149. PubMed Abstract | Publisher Full Text OpenURL

  30. Sen D, Gilbert W: Formation of parallel four-stranded complexes by guanine-rich motifs in DNA and its implications for meiosis.

    Nature 1988, 334(6180):364-366. PubMed Abstract | Publisher Full Text OpenURL

  31. Sundquist WI, Klug A: Telomeric DNA dimerizes by formation of guanine tetrads between hairpin loops.

    Nature 1989, 342(6251):825-829. PubMed Abstract | Publisher Full Text OpenURL

  32. Kunkel TA: Misalignment-mediated DNA synthesis errors.

    Biochemistry 1990, 29(35):8003-8011. PubMed Abstract | Publisher Full Text OpenURL

  33. Markham NR, Zuker M: DINAMelt web server for nucleic acid melting prediction.

    Nucleic Acids Res 2005, 33(Web Server issue):W577—W581. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Feldkamp U, Wacker R, Schroeder H, Banzhaf W, Niemeyer CM: Microarray-based in vitro evaluation of DNA oligomer libraries designed in silico.

    ChemPhysChem 2004, 5(3):367-372. PubMed Abstract | Publisher Full Text OpenURL

  35. Penchovsky R, Ackermann J: DNA library design for molecular computation.

    J Comput Biol 2003, 10(2):215-229. PubMed Abstract | Publisher Full Text OpenURL

  36. Dunbar SA: Applications of LuminexⓇ xMAPTM technology for rapid, high-throughput multiplexed nucleic acid detection.

    Clin Chim Acta 2006, 363:71-82. PubMed Abstract | Publisher Full Text OpenURL

  37. Tørring T, Voigt NV, Nangreave J, Yan H, Gothelf KV: DNA origami: a quantum leap for self-assembly of complex structures.

    Chem Soc Rev 2011, 40(12):5636-5646. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  38. Douglas SM, Marblestone AH, Teerapittayanon S, Vazquez A, Church GM, Shih WM: Rapid prototyping of 3D DNA-origami shapes with caDNAno.

    Nucleic Acids Res 2009, 37(15):5001-5006. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  39. Yanisch-Perron C, Vieira J, Messing J: Improved M13 phage cloning vectors and host strains: nucleotide sequences of the M13mp18 and pUC19 vectors.

    Gene 1985, 33:103-119. PubMed Abstract | Publisher Full Text OpenURL

  40. Ding B, Deng Z, Yan H, Cabrini S, Zuckermann RN, Bokor J: Gold nanoparticle self-similar chain structure organized by DNA origami.

    J Am Chem Soc 2010, 132(10):3248-3249. PubMed Abstract | Publisher Full Text OpenURL

  41. Stearns LA, Chhabra R, Sharma J, Liu Y, Petuskey WT, Yan H, Chaput JC: Template-directed nucleation and growth of inorganic nanoparticles on DNA scaffolds.

    Angew Chem, Int Ed 2009, 48(45):8494-8496. Publisher Full Text OpenURL