Stability of double-stranded oligonucleotide DNA with a bulged loop: a microarray study
© Trapp et al; licensee BioMed Central Ltd. 2011
Received: 2 August 2011
Accepted: 13 December 2011
Published: 13 December 2011
Skip to main content
© Trapp et al; licensee BioMed Central Ltd. 2011
Received: 2 August 2011
Accepted: 13 December 2011
Published: 13 December 2011
DNA is a carrier of biological information. The hybridization process, the formation of the DNA double-helix from single-strands with complementary sequences, is important for all living cells. DNA microarrays, among other biotechnologies such as PCR, rely on DNA hybridization. However, to date the thermodynamics of hybridization is only partly understood. Here we address, experimentally and theoretically, the hybridization of oligonucleotide strands of unequal lengths, which form a bulged loop upon hybridization. For our study we use in-house synthesized DNA microarrays.
We synthesize a microarray with additional thymine bases in the probe sequence motifs so that bulged loops occur upon target hybridization. We observe a monotonic decrease of the fluorescence signal of the hybridized strands with increasing length of the bulged loop. This corresponds to a decrease in duplex binding affinity within the considered loop lengths of one to thirteen bases. By varying the position of the bulged loop along the DNA duplex, we observe a symmetric signal variation with respect to the center of the strand. We reproduce the experimental results well using a molecular zipper model at thermal equilibrium. However, binding states between both strands, which emerge through duplex opening at the position of the bulged loop, need to be taken into account.
We show that stable DNA duplexes with a bulged loop can form from short strands of unequal length and they contribute substantially to the fluorescence intensity from the hybridized strands on a microarray. In order to reproduce the result with the help of equilibrium thermodynamics, it is essential (and to a good approximation sufficient) to consider duplex opening not only at the ends but also at the position of the bulged loop. Although the thermodynamic parameters used in this study are taken from hybridization experiments in solution, these parameters fit our DNA microarray data well.
The hybridization process - the formation of the well-known double-helix structure from two complementary nucleic acid strands (such that A · T and C · G base pairs are formed) - is pivotal to the living organism. Among other important biotechnological methods, PCR or DNA microarray technology rely on it.
DNA microarrays consist of regular spaced domains of surface-attached probe sequences, which act as binding sites for their complementary fluorescently-labeled target sequences in solution. The probe sequence and position of each domain on the surface is known and the amount of bound target DNA can be determined quantitatively. Microarrays are important in many biotechnological methods such as gene expression profiling, where complex target oligonucleotides mixtures need to be analyzed in a highly parallel manner [1–3].
Due to the very sensitive molecular recognition process of DNA, one is in principle able to detect even small sequence deviations with the help of DNA microarrays. However, DNA targets that are not perfectly complementary can also form duplexes with the surface bound probes, albeit are less stable than the perfectly matching correspondent (PM).
Although DNA microarrays are widely used in biological and biotechnological applications, the underlying physical mechanisms of the hybridization process are poorly understood. Data analysis is mostly based on empirical, statistical methods [4–6]. To fully exploit the potential of the DNA microarray technology, it is desirable to pursue a more fundamental approach to the stability of hybridized or partly hybridized strands. Molecular simulations have greatly increased our understanding of DNA dynamics, thermal fluctuations and hybridization. DNA hybridization and mechanical properties of DNA e.g. the persistence length in the presence of surfaces were investigated on the molecular level [7–9]. While molecular simulations give a very detailed view of the molecular dynamics, here we are interested in a simple scheme to assess the stability of bulged loops on a DNA microarray. Systematic experiments on short bulged loops have hardly been performed.
The standard model for hybridization in solution is the so called two-state-nearest-neighbor model (NN-Model), which treats the formation of the DNA duplex as a two-state process where the duplex is either fully hybridized or fully denatured [10, 11]. The model calculates the binding free energy of a perfectly complementary double-stranded duplex by summing the nearest-neighbor interaction parameters (10 experimentally determined free energy parameters [12–14]). These parameters take into account, that DNA stability arises from hydrogen bonding and base stacking interactions. Furthermore it is possible to extend the model and include single base mismatch (MM) defect parameters . This model proved very successful for the prediction of duplex melting temperature T m in solution.
In several experiments, incorporated MMs have had a position dependent influence on the fluorescent signal. Zhang et al. suggested the position-dependent-nearest-neighbor-model (PDNN)  where the binding free energy of the duplex is calculated as a weighted sum of the nearest-neighbor parameters. The weight parameters are determined empirically.
In the past we experimentally and theoretically investigated the effect of single MMs on the duplex stability of a DNA microarray in the case where the lengths of probe and target match [16–18]. We have shown that a two state NN-model could not predict the MM binding affinities precisely. Therefore we developed a different theoretical approach, based on a double-ended molecular zipper [19–21]. The double-ended molecular zipper considers, that the duplex can only open from the ends. This simplification is justified because base pairs, which are located away from the duplex ends are less stable. This holds even if a single MM is incorporated into the duplex. Taking into account the heterogeneity of the binding affinities due to synthesis defects, the DNA microarray data could be reproduced with the model. We have shown that the double-ended zipper model maps to the PDNN model, while the former is derived from first principles . The purpose of this study is to investigate the case where probe and target have unequal lengths and bulged loops form upon hybridization. Bulged loops are referred to as loops in the following. With our DNA microarray setup, loops of different lengths and at different positions can be obtained in a controlled manner by inserting additional bases into the perfectly matching probe sequence. The formation of loops increases the complexity of the hybridized state: new binding states between probe and target strands may emerge. We show that a good reproduction of the experimental data remains possible with the molecular zipper, but only if duplex opening can also occur at the loop position.
We use in-house synthesized DNA Microarrays. All employed protocols including the preparation of dendrimer-functionalized microarray substrates, the light-directed synthesis (a "maskless" photolithographic technique based on NPPOC-phosphoramidites), as well as the data analysis methods are provided in Naiser et al. . The only difference to the previously published experimental setup is a more homogeneous illumination of the microarray surface as well as an increased resolution due to the improved optics.
Target sequence Cy3-labeled perfect matching target sequence in solution.
Images of the hybridized DNA microarray are taken for data analysis after thermal equilibrium is reached.
In order to test if the microarray surface has a significant influence on the hybridization, we repeat our experiments with the reversed probe sequence. The 5'-end of the sequence employed throughout this work corresponds to the 3'-end of the test sequence. No influence of the microarray surface on the hybridization could be detected (see additional file 1: Influence of the microarray surface on the hybridization signal).
We also synthesized probes with other loop sequences than the discussed poly-T. We investigated the influence of poly-C-sequences and random sequences on duplex stability as a function of loop length. The results are provided in additional file 2: Duplex stability of DNA duplexes with bulged loops of different sequences as a function of loop length. We didn't observe a significant change in the dependence of the fluorescent signal as a function of loop length as compared to the poly-T-sequences.
Figure 4b shows the measured hybridization signals as a function of loop position in 3' to 5' direction after averaging over all loop lengths (PM signal set to 1). The resulting "loop position defect profile" is symmetric with respect to the center of the duplex. The signal is strongest for loops at the end, as well as in the middle of the duplex, it is weakest for loops at a distance of about 3-4 bases from the center. The difference between maximum and minimum is about 10% only. This is a weak variation compared to the hybridization signal as a function of loop length.
Loops positioned close to either end of a duplex have less potential binding sites towards that end, and they can open to form a dangling end. However, in this case a large part of the duplex to the opposite side of the loop remains strongly bound (Figure 5a).
Loops located at a center position have many possible binding partners to the left and to the right resulting in a closed loop and higher duplex stability (Figure 5b).
In between both of the extremes above, the hybridization signal drops to a minimum. This is because on one side these loops have less binding partners than loops in the middle of the duplex. On the other, the large hybridized part is shorter than for loops occupying end positions (Figure 5c).
where K is the equilibrium binding constant of the probe-target duplex. Since the fluorescent signal of the array is proportional to the fraction of hybridized probes θ, we think of θ as the "hybridization signal" in the following.
To minimize computation time, synthesis defects are only considered up to a certain maximum number per strand . The bases of the loop are treated synthesis defect free. Since the bases in the loop are, most of the time, only weakly or not at all bound (there are almost no complementary bases in the target strand), the consideration of synthesis defects in the loop is not necessary.
In our case, N ' = 33, we took up to synthesis defects into account. This generates 6018 different probe sequences. Strands with more than 3 synthesis defects can be neglected (see additional file 3: Influence of the number of MMs on the fluorescent signal).
In the following, we calculate the binding constants K i as a function of loop position P and loop length L in thermodynamic equilibrium.
Δg init = -4.5 kcal/mol is the duplex initialization free energy . For the binding constant K i (P, L), K i (P, L) = Z D (P, L) taking the totally denatured state, S 0, as the reference state.
σ = 1, 75 · 10-4 is the so-called cooperativity parameter, μ is the connectivity constant and c = 2, 15 is the loop closure exponent. σ and c are universal constants whereas μ (μ = 4, 684 used here) depends on the considered geometry.
where L ranges from 1 to 13.
So far duplex opening was only possible from the ends of the duplex. States, in which the duplex zips at the loop position, are essential for the correct reproduction of our experimental results.
The partition function of a duplex Z on (P, L) as a function of loop position P and loop length L can be decomposed as a sum of five elements:
Z zipper (P, L): The hybridized strands zip from both ends. Partition function as presented in the section above.
Z extended,right (P, L): Probe and target strands right of the loop position can undergo every possible binding configuration among each other (not limited to a zipper). Thus, loops of different size in probe and target strand can appear. The duplex part left of the loop position zips to and from the loop position. The free energy of this part is considered in ΔG left . Figure 8 illustrates Z extended,right . The red and green dashed lines represent hybridized duplex parts. The black dashed lines are denatured (at the end) or they form a loop between probe and target (middle). Here in this particular case, the probe and target strand form a loop of 16 and 11 bases respectively. The two strands reunite after base 28 of the probe strand and base 23 of the target strand, the following 7 bases are hybridized. This results in the free energy ΔG7,28,23.
Z extended,left (P, L): analogous to
Z extended,right (P, L) but opposite side.
Z double zipper (P, L): Both parts, left and right from the loop position behave like an independent zipper. To avoid double count of states from adding Z extended,right (P, L) and Z extended,left (P, L), this partition function needs to be subtracted.
i and j mark the positions of the zipper forks in probe and target respectively, n bases of probe and target strand (n - 1 NN pairs) of the region to the right of i and j are hybridized. Thus, ΔG n,i,j is the NN energy of (n - 1) base pairs which are hybridized from zipper fork position i in the probe and zipper fork position j in the target. is the NN energy of a single hybridized base pair in this region.
Constants γ and μ are defined as before. ν = 0, 588 ± 1, 5 · 10-3 is the (universal) metric exponent. # total (M i ) is the total number of SAWs of M i steps (as defined in equation (9)).
In the case where probe and target length match, duplex zipping can only occur if the two strands are perfectly aligned. We consider the initiation energy, the entropic barrier to meet this constraint, as constant. We simply write K i = Z D and include the initiation energy in a prefactor.
L0 is the characteristic length of the problem, which is the persistence length (in our experimental conditions this corresponds to a single base). In the case of a short, loop-forming sequence located in the center of the strand, however, there are two positions, where parallel but shifted probe and target strands can initiate duplex formation. These positions correspond to the matching sequence left and the right from the loop implying a correction of ΔL/L0 by 1/2. However, if the loop forms towards the ends, we are close to the situation of a single strand above. In the following we neglect this dependence on loop position and use a factor 1/2 throughout. Either factor (1 or 1/2) does not drastically modify our result, if the factor C is adjusted accordingly.
Our approximation for ΔS init tends to overestimate the corresponding initiation energy penalty as ΔL increases. This is because for large ΔL the situation differs: in this case the separated matching sequences are almost independent and the initiation energy tends to its asymptotic value of two independent hybridization events. As a conclusion for large ΔL a weaker dependence of the initiation energy on ΔL can be expected.
The calculation of the hybridization signal is then straight forward.
We note, that the choice of the denominator of equation (26) following from (25) has an impact on the calculated hybridization signals. Our theory could possibly be improved by choosing a different denominator which, however, may be a subtle problem by itself, not the scope of this paper.
Figure 10 and 11 show that the model reproduces our experimental findings well. Parameters used here were: simulation temperature T sym = 317 K, synthesis error rate p = 0.084, energy penalty for synthesis related defects Δg def,syn = -1 kcal/mol (consideration up to three errors per probe during synthesis, Δg def,syn was determined in ). We use the temperature adjusted NN and MM defect parameters from [12, 13] and the references therein. Since MM defect parameters are only available for isolated MMs, we include another parameter MM def = -2 kcal/mol for the case of two adjacent MMs (we approximate two adjacent MMs as two independent synthesis defects next to each other, therefore: MM defect = 2.g def,syn . Furthermore, we use the (universal) parameters for a SAW . The only free parameters are the factor C = 1.5 · 10-3 that links the theoretical binding constants K i to fit our experimental data of fluorescent signals and the probability for synthesis related defects p = 0.084 (the latter is not completely free since it is used to check if our theory is consistent with the coupling and deprotection efficiency of the used oligonucleotides).
We note that the partition function Z double zipper (P, L) alone already reproduces the approximate shape of the symmetric loop defect profile as shown in Figure 11b (dependence of the hybridization signal on loop position). However, the resulting binding constants are smaller than the ones calculated with Z extended,left and Z extended,right respectively. Z extended,left and Z extended,right help in reproducing the shape and moreover, the absolute values of the experimental signals (see additional file 4: Hybridization signals resulting from Z double zipper and comparison to Z extended,right + Z extended,left ).
Small differences between theoretical and experimental results regarding the signal dependence on loop position can be explained by the particularities of the duplex sequence under study. Here we look at two differences:
region ranging from loop position 14 to 18: this duplex region has many A/T bases and the distance between two C bases is the largest for the whole sequence. The duplex destabilization of an A/T rich region may be underestimated.
loop position 21: the region has many C bases and the loop bases are inserted after two existing C bases. It has been shown [22, 29, 30], that degenerated base pairs may reinforce binding considerably. Stabilization by degenerated base pairs is not included in our theory.
Although there are differences between experiment and theory, the deviations are small (see Figure 10 and 11). An even better agreement could be obtained by choosing a different dependence of the duplex initiation energy on ΔL. Our approximation for it (see above) only holds for short ΔLand we suppose the systematic deviation visible in Figure 11a from theory and experiment to originate from our approximation. As expected, at longer ΔL, we tend to underestimate the binding constant. To our knowledge, although an often encountered problem, no simple scheme to assess the initiation energy is known. Working out the dependence of the initiation energy between the two regimes discussed above (short and very long ΔL) is beyond the scope of this paper. Molecular simulations could help to provide better understanding of the nucleation process .
In literature, internal DNA loops or bubbles of total length l = l 1 + l 2 e.g. occurring in DNA denaturation experiments are often treated as SAWs of the same length returning to their origin (l 1 : unbound bases in probe; l 2 : unbound bases in target) . Reproduction of our experimental data could not be achieved when the calculation is done in this way, because the calculated loop energy penalties were much too large. Treating a DNA loop as a SAW of length l = l 1 + l 2 returning to the origin is different from calculating the probability that two SAWs of given lengths l 1 and l 2 start at the same point and meet again at some distance. In the first case, the number of possible conformations is much higher because the constraint is weakened to any pair with , not just the given l 1, l 2. The first case could give the same same results if the calculation is done under the constraint that the loop of length l = l 1 + l 2 reaches the position where the two loops reunite after l 1 steps similar to the way described in .
This may not always matter so much: the length of the probe sequences used throughout this study is much shorter than the length of DNA strands used in DNA denaturation experiments. Since the free energy of a short DNA strand is small, the size of the loop energy penalties is more crucial.
In this paper we investigated the stability of DNA with a bulged loop. We inserted additional thymine bases into the surface-bound PM motif at a given position. By hybridizing DNA oligonucleotide targets onto the DNA microarray, bulged loops of different length and at different positions along the DNA duplex are formed.
We find that duplex stability decreases monotonically with the length of the bulged loop. Moreover, if the position of the bulged loop on the probe strand is varied, duplex stability exhibits a symmetric variation with respect to the center of the duplex. Duplex stability is highest for end- and middle-positions of the inserted bulged loop. For theoretical prediction we have shown that it is necessary and sufficient to consider strand opening at the position of the bulged loop. We have elaborated a successful approximation for the partition function of these new binding states. The signal dependence on loop length and on loop position could be reproduced with a limited amount of computing time (see Figure 11).
The employed NN free energy parameters from  are based on solution hybridization experiments. However, as we show in this study and in a previous paper , these parameters can be used to describe microarray hybridization well. The corresponding loop energy penalties can be obtained by considering the bulged loops as a self-avoiding walk on a lattice.
In our simulation, we use just two free parameters:
C = 1.5 · 10-3: scaling factor, which fits the calculated binding constants to the fluorescent light intensities. This parameter cannot be avoided.
p = 0.084: probability of a synthesis related defect. In a previous work, the value of p was determined to p = 0, 1. In our improved experimental setup, we have less stray light and a better resolution which result in a better synthesis quality (see Methods). Therefore, we chose p to be a free parameter. p is obtained as 0.084 in good agreement with the coupling and deprotection efficiency of the employed oligonucleotides and the achievable contrast of the optical setup [32–34]. Given this knowledge, p is not completely free and the resulting value is used to check the consistency of our theory.
The formation of bulged loops is an important aspect that needs to be considered when analyzing DNA microarray data or DNA hybridization of complex mixtures in general. Partly non-complementary sequences can form stable complementary duplexes through formation of a bulged loop resulting in false positive signals. The investigation of these bulged loop structures is therefore necessary to gain a deeper understanding of DNA hybridization and to make DNA microarrays and other, nucleic acid based high throughput technology based on DNA hybridization more reliable and accurate.
The authors thank Jona Kayser for many helpful discussions on this work. Our research was supported by the University of Saarland.
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.