DNA base pairs are known to open more easily at the helix terminal, a process usually called end fraying, the details of which are still poorly understood. Here, we present a mesoscopic model calculation based on available experimental data where we consider separately the terminal base pairs of a DNA duplex. Our results show an important reduction of hydrogen bond strength for terminal cytosine-guanine (CG) base pairs which is uniform over the whole range of salt concentrations, while for AT base pairs, we obtain a nearly 1/3 reduction but only at low salt concentrations. At higher salt concentrations, terminal adenine-thymine (AT) pair has almost the same hydrogen bond strength than interior bases. The calculated terminal stacking interaction parameters display some peculiarly contrasting behavior. While there is mostly no perceptible difference to internal stacking, for some cases, we observe an unusually strong dependence with salt concentration which does not appear follow any pattern or trend.

We calculate the nearest-neighbour enthalpies and entropies at 5 salt concentrations of 18 RNA sequences, each for at least 9 different species concentrations, totalling 757 melting temperatures, using a melting temperature optimization method. These new parameters do not need to be saltcorrected and are shown to provide overall improved melting temperature predictions. They show a marked quadratic dependence with salt concentrations which are compensated to form linear Gibbs free energies. Two different parameter schemes were tested, with fixed or variable initial parameters. We have found that using variable initial parameters provides better predictive results than using salt correction factors and that the prediction uncertainty is considerably reduced for a validation set of independent sequences. An interpolation scheme is introduced to generate model parameters for arbitrary salt concentrations which performs better against a validation set than predictions using salt corrections.

No abstract

The use of mesoscopic models to describe the thermodynamic properties of locked nucleic acid (LNA)-modified nucleotides can provide useful insights into their properties, such as hydrogen-bonding and stacking interactions. In addition, the mesoscopic parameters can be used to optimize LNA insertion in probes, to achieve accurate melting temperature predictions, and to obtain duplex opening profiles at the base-pair level. Here, we applied this type of model to parameterize a large set of melting temperatures for LNA-modified sequences, from published sources, covering all possible nearest-neighbor configurations. We have found a very large increase in Morse potentials, which indicates very strong hydrogen bonding as the main cause of improved LNA thermodynamic stability. LNA-modified adenine−thymine (AT) was found to have similar hydrogen bonding to unmodified cytosine−guanine (CG) base pairs, while for LNA CG, we found exceptionally large hydrogen bonding. In contrast, stacking interactions, which were thought to be behind the stability of LNA, were similar to unmodified DNA in most cases. We applied the new LNA parameters to the design of BRAF, KRAS, and EGFR oncogene variants by testing all possible LNA modifications. Selected sequences were then synthesized and had their hybridization temperatures measured, achieving a prediction accuracy within 1 °C. We performed a detailed base-pair opening analysis to discuss specific aspects of these probe hybridizations that may be relevant for probe design.