|
|
||||||||
Articles |
1
Department of Clinical Chemistry, Georg-August University, Robert Koch Strasse 40, 37075 Goettingen, Germany.
2
Department of Pathology, University of Utah Medical School, Salt Lake City, UT 84132.
3
Chronic Illness Research Foundation, San Francisco, CA 94107.
aAuthor for correspondence. Fax 49-551-39-12504; e-mail nahsen{at}gwdg.de
| Abstract |
|---|
|
|
|---|
Methods: The experimental Tms of 475 matched or mismatched target/probe duplexes were obtained in our laboratories or were compiled from the literature based on studies using the same real-time PCR platform. This data set was used to evaluate the contributions of [Mg2+], [dNTPs], and [DMSO] in N-N calculations. In addition, best-fit coefficients for common empirical formulas based on GC content, length, and the equivalent sodium ion concentration of cations [Na+eq] were obtained by multiple regression.
Results: When we used [Na+eq] = [Monovalent cations] + 120(
) (the concentrations in this formula are mmol/L) to correct
S0 and a DMSO term of 0.75 °C (%DMSO), the SE of the N-N Tm estimate was 1.76 °C for perfectly matched duplexes (n = 217). Alternatively, the empirical formula Tm (°C) = 77.1 °C + 11.7 x log[Na+eq] + 0.41(%GC) - 528/bp - 0.75 °C(%DMSO) gave a slightly higher SE of 1.87 °C. When all duplexes (matched and mismatched; n = 475) were included in N-N calculations, the SE was 2.06 °C.
Conclusions: This robust model, accounting for the effects of Mg2+, DMSO, and dNTPs on oligonucleotide Tm in PCR, gives reliable Tm predictions using thermodynamic N-N calculations or empirical formulas.
| Introduction |
|---|
|
|
|---|
The most accurate prediction of Tm for oligonucleotide DNA uses the thermodynamic nearest-neighbor (N-N) model [see Ref. (9) for review and parameters]. N-N calculations for Tm prediction are useful on microarrays (10) and for the selection of PCR primers and hybridization probes (11). Empirical data from probe melting curve analysis during real-time PCR correlate well with theoretical predictions (12)(13). The N-N model is based on the assumption that probe hybridization energy can be calculated from the enthalpy and entropy of all N-N pairs, including a contribution from each dangling end. Thermodynamic values for entropy and enthalpy of each possible matched N-N and dangling ends have been determined (9)(14). Dangling-end effects account for the stacking energy of a shorter probe on a more lengthy target (15)(16). The entropy (
S) is salt dependent, and
S0 must be corrected if the ionic environment is different from 1 mol/L NaCl, the salt concentration at which most thermodynamic values have been derived. However, Mg2+ is present in PCR as an important cofactor for Taq DNA polymerase and strongly influences
S (17). Deoxynucleotide triphosphates (dNTPs) are also essential and chelate some of the available Mg2+ (18). In addition, dimethyl sulfoxide (DMSO) is commonly used as a cosolvent (19) to facilitate amplification from difficult templates. Addition of DMSO decreases the Tm (20)(21)(22), which must be taken into account when primer Tm is calculated (23).
To get a deeper insight into probe Tm under common PCR buffer conditions, we have compiled Tm data from 475 different DNA duplexes. These assays were performed with different concentrations of Mg2+ and DMSO, reflecting current PCR laboratory practice. We provide an empirical
S compensation for the Mg2+ and dNTP influence on ionic strength. In addition, best-fit coefficients for simpler formulas based on GC content, length, and the equivalent sodium ion concentration are determined for convenient bench-side use that offer accurate Tm predictions.
| Materials and Methods |
|---|
|
|
|---|
data collection
All assays were performed using the LightCycler real-time PCR instrument. Most Tms were measured in our laboratories (n = 388) during the course of genotyping experiments [for examples, see Refs (12)(28)]. Some additional Tms were extracted from the literature when complete experimental conditions were published (n = 87). In total, 162 different probes were used with various templates and conditions, including 221 completely matched hybridizations, 237 single mismatches, and 17 two-point mismatches. Forty assays were based on melting curves detected with SYBR Green I, whereas the remainder were based on the melting of fluorescent oligonucleotide probes. DMSO was used in 206 assays in concentrations ranging from 2.5% to 10%. These data are available as an online supplement at Clinical Chemistry Online (http://www.clinchem.org/content/vol47/issue11).
statistical analysis
N-N calculations.
The entropy and enthalpy were calculated from probe sequences at standard conditions (1 mol/L NaCl) as described in more detail elsewhere (9)(29)(30). In addition we considered the published thermodynamic data for dangling-end contributions (14). Mismatches were accounted for by the thermodynamic data reported by Allawi and SantaLucia (31)(31)(32)(33)(34) and Peyret et al. (35). The PCR DNA target concentration was set to 50 nM. SYBR Green I, if present, was assumed to increase the Tm by 1 °C at 1:20 000 dilution, based on own preliminary data. Calculations were performed with ExcelTM for Windows (Microsoft), using the built-in statistical functions. The Pearson r2 was used for correlations, and standard linear regression was used for relating observed to measured Tm. Thermodynamic N-N stability calculations were performed using MeltCalc, a spreadsheet software for Excel (36).
The observed Tm was used as the dependent variable in a multiple variable fit to determine the DMSO coefficient and the best formula for Mg2+ influence. Prior evidence suggested that the influence of Mg2+ on Tm is stoichiometrically reduced by dNTPs. Initial calculations based on r2 indicated that the relationship of DMSO to Tm was linear, whereas [Mg2+] was nonlinear. Therefore, our model was: Tm (observed) = Tm (predicted) - a x DMSO (%), with [Na+eq] = [monovalent cations] + b x ([Mg2+] - [dNTP])c. The parameters a, b, and c were optimized to minimize the prediction error by stepwise incremental iterations. With a = 0.75, b = 120, and c = 0.5, only 18% of the predicted values fell outside a 5% error limit. The nonlinear effect of [Mg2+] on [Na+eq] was best approximated by the square-root function, which is in agreement with a previous report (37).
Alternative formulas.
Using our empirical data set, we evaluated several simpler formulas for their ability to predict Tm. These formulas cannot properly account for the presence of single mismatches. Therefore, only data for matched probe/template duplexes were used.
The WallaceIkatura rule is often used as a rule of thumb when primer Tm is to be estimated at the bench (1)(38). However, the formula was originally applied to the hybridization of probes in 1 mol/L NaCl (1) and is an estimate of the denaturation temperature (Td):
![]() | (1) |
![]() | (2) |
Marmur and Doty (2) originally established a formula to correlate GC content (%GC) to the Tm of long duplexes at a given ionic strength. Chester and Marshak (23) added a term to account for DNA strand length (n in base pairs) to estimate primer Tm:
![]() | (3) |
![]() | (4) |
Another modification is that of Wetmur (1):
![]() | (5) |
![]() | (6) |
| Results |
|---|
|
|
|---|
statistical analysis
N-N calculations.
The best fit values for [Mg2+] and [dNTP] coefficients revealed that the Na+ equivalents (Na+eq) were approximated by:
![]() | (7) |
S0 was corrected for the salt concentration as follows (9):
![]() | (8) |
Each percentage of DMSO (by volume) decreased the Tm by 0.75 °C. GC content had no obvious influence on the DMSO factor. Using these equations, we found a good regression of predicted vs observed Tm (Fig. 1
). The mean prediction error was 0.2 ± 2.18 °C, which is within the error range for N-N calculations.
|
Alternative formulas.
The WallaceIkatura rule (Eq. 1
) overestimates the Tm of long duplexes and gives reasonable results only in the range of 1420 bp (1)(38). Therefore, only duplexes shorter than 21 bp were included in the analysis (Table 1
). When the analysis was extended to include duplexes of up to 24 bp, r2 decreased to 0.64.
|
The equation for the effective priming temperature by Wu et al. (25) (Eq. 2
) is similar to the WallaceIkatura rule. Only oligonucleotides with Ln <39 (see Eq. 2
) were included in the analysis, as suggested by the authors. The oligonucleotide Tm is overestimated by a mean of 2.5 °C by this formula (Table 1
).
Including the length dependence of Chester and Marshak [Eq. 3
; (23)] improved r2, but the intercept and difference between observed and predicted Tms were poor (data not shown). However, a better fit was obtained with new constants obtained from stepwise iterations (fit variable: length dependency term):
![]() | (9) |
![]() | (10) |
![]() | (11) |
Eq. 6
has been suggested as the best predictor of polymer DNA Tm (41)(42). The best fit of constants in this formula based on our oligonucleotide data produced (fit variables: Tm offset, oligomer [Na+] dependency term, length dependency term added):
![]() | (12) |
| Discussion |
|---|
|
|
|---|
). The N-N model accurately predicts probe Tm in the LightCycler analysis system (12). We compiled and analyzed measured Tm data for matched and mismatched hybridization probes from different laboratories. Therefore, some interlaboratory variation is expected from different reagents and protocols. The temperature transition rates during melting curve acquisition (usually 0.10.2 °C/s) are too fast to achieve equilibrium conditions and cause a slight overestimation of the true Tm (43).
Probe sequence choices reflect the demands of mutation detection and are not designed particularly for two-state behavior. Even with these limitations, the predictive accuracy we have achieved underscores the robustness of the parameterization. Systematic errors introduced by fluorescent dyes were negligible in another study (10), and fluorescence resonance energy transfer probes themselves have successfully been used to derive thermodynamic parameters (44). Our findings should be equally applicable to probe and primer oligonucleotides, thereby allowing in silico optimization of primers and probes, saving both on time required for optimization and costs for probe resynthesis (12)(29)(45)(46)(47).
An extensive survey of alternative formulas used for the prediction of perfectly matched oligonucleotide DNA was also performed. The WallaceIkatura rule (Eq. 1
) is often used as a rough predictor of primer Tm but has limited accuracy, especially for longer oligonucleotides. This rule assumes a salt concentration of 1 mol/L NaCl, which is typical for dot blots and other hybridizations but not PCR. The fact that it is used for PCR applications is more a testament of the robustness of PCR toward different annealing temperatures than evidence for accurate Tm estimates. The same is true for the formula of Wu et al. (25). Many formulas were originally designed to relate measured Tm and GC content of polymer DNA (Eqs. 36
). The inclusion of additional terms for ionic strength (6), length dependency (5), and GC dependency of ionic strength (39) has led to more accurate estimates for polymer DNA. Eqs. 912
have been specifically optimized for oligonucleotide Tm estimation by best-fit estimates of our data set. Eq. 10
is recommended as a tradeoff between accuracy and ease of use. Table 2
gives primer Tms for common primer compositions and PCR conditions. These estimates may not be accurate for certain sequences with a biased N-N composition (41). Furthermore, mismatches are strongly dependent on their N-N bases and require more laborious N-N calculations.
|
The effect of DMSO on thermal stability of DNA has been investigated before. Our factor of 0.75 °C decrease in Tm per 1% DMSO is similar to previous findings of 0.6 °C per 1% DMSO (20), 0.675 °C per 1% DMSO (22), and 0.5 °C per 1% DMSO (21). These prior studies were performed on polymer DNA, suggesting that DMSO may have a slightly greater effect on oligomer DNA.
Because template priming during PCR is a kinetic process, efficient, specific priming should occur at the primer Tm (25), suggesting that the Tm can be used as the annealing temperature in PCR. Because efficient amplification is dependent on hybridization of both primers, it is rational to use the Tm of the least stable primer. The temperature used for annealing in PCR also depends on the annealing time. For example, allele-specific amplification can be achieved with rapid cycling ("0" s annealing) at a lower annealing temperature than conventional cycling with a longer annealing time (48). Finally, Tm is not just a property of an oligonucleotide, but a property of an oligonucleotide under specific conditions and at a given concentration.
In conclusion, we have developed a robust model for the effects of Mg2+, DMSO, and dNTPs on oligonucleotide Tm under common PCR buffer conditions. This enables reliable Tm predictions and in silico primer and probe optimization using thermodynamic N-N calculations. We performed an extensive evaluation of different equations advocated for PCR primer Tm prediction. These formulas have been parameterized to accommodate for standard PCR conditions and are now useful for rapid calculation of the Tm of perfectly matched oligonucleotides. It is rational to use the Tm of the least stable PCR primer as the annealing temperature in PCR.
| Acknowledgments |
|---|
| Footnotes |
|---|
| References |
|---|
|
|
|---|
F508 locus. Clin Chem 1993;39:804-809.The following articles in journals at HighWire Press have cited this article:
![]() |
H. Ye, M. C. Huang, M.-H. Li, and J. Y. Ying Experimental analysis of gene assembly with TopDown one-step real-time gene synthesis Nucleic Acids Res., April 1, 2009; 37(7): e51 - e51. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. Qu, Z. Shen, D. Zhao, Y. Yang, and C. Zhang MFEprimer: multiple factor evaluation of the specificity of PCR primers Bioinformatics, January 15, 2009; 25(2): 276 - 278. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Zhou, R. J. Errigo, H. Lu, M. A. Poritz, M. T. Seipp, and C. T. Wittwer Snapback Primer Genotyping with Saturating DNA Dye and Melting Analysis Clin. Chem., October 1, 2008; 54(10): 1648 - 1656. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Pont-Kingdon, R. L. Margraf, K. Sumner, A. Millson, E. Lyon, and E. Schutz Design and Application of Noncontinuously Binding Probes Used for Haplotyping and Genotyping Clin. Chem., June 1, 2008; 54(6): 990 - 999. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y.-F. Chen, R.-C. Chen, L.-Y. Tseng, E. Lin, Y.-K. Chan, and R.-H. Pan NTMG (N-terminal Truncated Mutants Generator for cDNA): an automatic multiplex PCR assays design for generating various N-terminal truncated cDNA mutants Nucleic Acids Res., July 13, 2007; 35(suppl_2): W66 - W70. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Koressaar and M. Remm Enhancements and modifications of primer design program Primer3 Bioinformatics, May 15, 2007; 23(10): 1289 - 1291. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. D. Levin, D. Fiala, M. F. Samala, J. D. Kahn, and R. J. Peterson Position-dependent effects of locked nucleic acid (LNA) on DNA sequencing and PCR primers Nucleic Acids Res., November 6, 2006; 34(20): e142 - e142. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. S. Robinson, P. T. Monis, and P. J. Dobson Rapid, Sensitive, and Discriminating Identification of Naegleria spp. by Real-Time PCR and Melting-Curve Analysis Appl. Envir. Microbiol., September 1, 2006; 72(9): 5857 - 5863. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Schutz, M. Scharfenstein, and B. Brenig Genotyping of Ovine Prion Protein Gene (PRNP) Variants by PCR with Melting Curve Analysis, Clin. Chem., July 1, 2006; 52(7): 1426 - 1429. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Chavali, A. Mahajan, R. Tabassum, S. Maiti, and D. Bharadwaj Oligonucleotide properties determination and primer designing: a critical examination of predictions Bioinformatics, October 15, 2005; 21(20): 3918 - 3925. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Zhou, A. N. Myers, J. G. Vandersteen, L. Wang, and C. T. Wittwer Closed-Tube Genotyping with Unlabeled Oligonucleotide Probes and a Saturating DNA Dye Clin. Chem., August 1, 2004; 50(8): 1328 - 1335. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Liew, R. Pryor, R. Palais, C. Meadows, M. Erali, E. Lyon, and C. Wittwer Genotyping of Single-Nucleotide Polymorphisms by High-Resolution Melting of Small Amplicons Clin. Chem., July 1, 2004; 50(7): 1156 - 1164. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Wang and B. Seed A PCR primer bank for quantitative gene expression analysis Nucleic Acids Res., December 15, 2003; 31(24): e154 - e154. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Muller, K. Zopf, J. Bachofer, and W. Steimer Optimized Strategy for Rapid Cytochrome P450 2D6 Genotyping by Real-Time Long PCR Clin. Chem., October 1, 2003; 49(10): 1624 - 1631. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Maheshri and D. V. Schaffer Computational and experimental analysis of DNA shuffling PNAS, March 18, 2003; 100(6): 3071 - 3076. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. N. Gundry, J. G. Vandersteen, G. H. Reed, R. J. Pryor, J. Chen, and C. T. Wittwer Amplicon Melting Analysis with Labeled Primers: A Closed-Tube Method for Differentiating Homozygotes and Heterozygotes Clin. Chem., March 1, 2003; 49(3): 396 - 406. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |