Thermodynamics of internal C.T mismatches in DNA.

Thermodynamics of 23 oligonucleotides with internal single C.T mismatches were obtained by measuring UV absorbance as a function of temperature. Results from these 23 duplexes were combined with three measurements from the literature to derive nearest-neighbor thermodynamic parameters for seven linearly independent trimer sequences with internal C.T mismatches. The data show that the nearest-neighbor model is adequate for predicting thermodynamics of oligonucleotides with internal C.T with average deviations for Delta G degrees37, Delta H degrees, Delta S degrees and T m of 6.4%, 9.9%, 10.6%, and 1.9 degreesC respectively. C.T mismatches destabilize the duplex in all sequence contexts. The thermodynamic contribution of C. T mismatches to duplex stability varies weakly depending on the orientation of the mismatch and its context and ranges from +1.02 kcal/mol for GCG/CTC and CCG/GTC to +1.95 kcal/mol for TCC/ATG.     

The thermodynamics of 3'-terminal pyrene and guanosine for the design of isoenergetic 2'-O-methyl-RNA-LNA chimeric oligonucleotide probes of RNA structure

To facilitate design of short isoenergetic hybridization probes for RNA, we report the influence of adding 5'- or 3'-terminal 2'-O-methylguanosine (GM), LNA-guanosine (GL), or 3'-terminal pyrene pseudo-nucleotide (PPN) on the thermodynamic stability of 2'-O-methyl-RNA/RNA (2'-O-Me-RNA/RNA) duplexes with sequences 5'CMGMGMCMAM/3'AAXGCCGUXAA, where X is A, C, G, or U. A 3'-terminal GM or GL added to the 2'-O-Me-RNA strand to form a G-A, G-G or G-U mismatch enhances thermodynamic stability (DeltaDeltaG degrees 37) of the 2'-O-Me-RNA/RNA duplexes on average by 0.7 and 1.5 kcal/mol, respectively. A 3'-terminal GM or GL in a GM-C or GL-C pair stabilizes the 2'-O-Me-RNA/RNA duplex by 2.6 and 3.4 kcal/mol, respectively. A 5'-terminal GM or GL in a G-A or G-G mismatch provided less stabilization in comparison with a 3'-terminal G-A or G-G mismatch, but more stabilization in a G-C or G-U pair. In contrast to guanosine derivatives, pyrene residue (P) as PPN at the 3'-terminal position enhances thermodynamic stability of the 2'-O-Me-RNA/RNA duplexes on average by 2.3 +/- 0.1 kcal/mol, relatively independent of the type of ribonucleotide placed in the opposite strand. The thermodynamic data can be applied to design 2'-O-Me-RNA/RNA duplexes with enhanced thermodynamic stability that is also sequence independent. This is useful for design of hybridization probes to interrogate RNA structure and/or expression by microarray and other methods.     

Development of thermodynamic models for simulating probe dissociation profiles in fluorescence in situ hybridization

Stringency in ribosomal RNA (rRNA)-targeted fluorescence in situ hybridization (FISH) is typically adjusted with formamide, and the optimum formamide concentration at which the probe can hybridize with the target rRNA, but not with rRNAs with mismatches, is to be found experimentally. This is a difficult task when target or closest non-target organisms are not available in pure culture, or when there are numerous non-targets of concern. The objective of this work was to formulate mechanistic models capable of simulating the effect of formamide on probe dissociation. Using a previously described equilibrium model of FISH [Yilmaz and Noguera (2004) Applied and Environmental Microbiology 70(12):7126-7139] as the basis, the effect of formamide on free energy changes of probe-target duplex formation (DeltaG(1)(0)) and folding of target region (DeltaG(3)(0)) was simulated to be linear, and models differing in the definitions of the slopes of these relationships (m(1) and m(3)) were calibrated using experimental dissociation profiles for 27 probes targeting the 16S rRNA of Escherichia coli (E. coli). A good level of predictive power was obtained when m(1) was linearly related to probe length and when m(3) was made proportional to DeltaG(3)(0). The effect of single mismatches on probe dissociation with formamide was also studied, although at a preliminary level. The expected changes in DeltaG(1)(0) with the introduction of mismatches were not sufficient to capture the overall trends of mismatched dissociation profiles. In conclusion, this study offers the first theoretical method to calculate dissociation profiles for perfectly matched probes, and suggests a direction to systematically evaluate the effect of formamide on mismatched probes.     

Thermodynamics of sequence-specific binding of PNA to DNA

For further characterization of the hybridization properties of peptide nucleic acids (PNAs), the thermodynamics of hybridization of mixed sequence PNA-DNA duplexes have been studied. We have characterized the binding of PNA to DNA in terms of binding affinity (perfectly matched duplexes) and sequence specificity of binding (singly mismatched duplexes) using mainly absorption hypochromicity melting curves and isothermal titration calorimetry. For perfectly sequence-matched duplexes of varying lengths (6-20 bp), the average free energy of binding (DeltaG degrees ) was determined to be -6.5+/-0.3 kJ mol(-1) bp(-1), corresponding to a microscopic binding constant of about 14 M(-1) bp(-1). A variety of single mismatches were introduced in 9- and 12-mer PNA-DNA duplexes. Melting temperatures (T(m)) of 9- and 12-mer PNA-DNA duplexes with a single mismatch dropped typically 15-20 degrees C relative to that of the perfectly matched sequence with a corresponding free energy penalty of about 15 kJ mol(-1) bp(-1). The average cost of a single mismatch is therefore estimated to be on the order of or larger than the gain of two matched base pairs, resulting in an apparent binding constant of only 0.02 M(-1) per mismatch. The impact of a mismatch was found to be dependent on the neighboring base pairs. To a first approximation, increasing the stability of the surrounding region, i.e., the distribution of A.T and G.C base pairs, decreases the effect of the introduced mismatch.     

Effect of LNA- and OMeN-modified oligonucleotide probes on the stability and discrimination of mismatched base pairs of duplexes

Locked nucleic acid (LNA) and 2'-O-methyl nucleotide (OMeN) are the most extensively studied nucleotide analogues. Although both LNA and OMeN are characterized by the C3'-endo sugar pucker conformation, which is dominant in A-form DNA and RNA nucleotides, they demonstrate different binding behaviours. Previous studies have focused attention on their properties of duplex stabilities, hybridization kinetics and resistance against nuclease digestion; however, their ability to discriminate mismatched hybridizations has been explored much less. In this study, LNA- and OMeN-modified oligonucleotide probes have been prepared and their effects on the DNA duplex stability have been examined: LNA modifications can enhance the duplex stability, whereas OMeN modifications reduce the duplex stability. Next, we studied how the LNA:DNA and OMeN:DNA mismatches reduced the duplex stability. Melting temperature measurement showed that different LNA:DNA or OMeN:DNA mismatches indeed influence the duplex stability differently. LNA purines can discriminate LNA:DNA mismatches more effectively than LNA pyrimidines as well as DNA nucleotides. Furthermore, we designed five LNA- and five OMeN-modified oligonucleotide probes to simulate realistic situations where target-probe duplexes contain a complementary LNA:DNA or OMeN:DNA base pairs and a DNA:DNA mismatch simultaneously. The measured collective effect showed that the duplex stability was enhanced by the complementary LNA:DNA base pair but decreased by the DNA:DNA mismatch in a position-dependent manner regardless of the chemical identity and position of the complementary LNA:DNA base pair. On the other hand, the OMeN-modified probes also showed that the duplex stability was reduced by both the OMeN modification and the OMeN:DNA mismatch in a position-dependent manner.     

Mutation detection by stacking hybridization on genosensor arrays

A new strategy for analysis of point mutations using oligonucleotide array (genosensor) hybridization was investigated. In the new approach, a single-stranded target strand is preannealed with a labeled "stacking oligonucleotide," and then the partially duplex labeled target molecule is hybridized to an array of glass-tethered oligonucleotide probes, targeted to the region on the target immediately adjacent to the stacking oligomer. In this configuration, the base-stacking interactions between the "capture probe" and the contiguously stacking oligomer stabilize the binding of the target molecule to its complementary probe on the genosensor array. The temperature of hybridization can be adjusted so that the target molecule will bind to the glass-tethered probe only in the presence of the stacking oligomer, and a single mismatch at or near the terminal position ol the capture probe disrupts the stacking interactions and thereby eliminates or greatly reduces the hybridization. This stacking hybridization approach was investigated using a collection of synthetic targets, probes, and stacking oligonucleotides, which permitted identification of conditions for optimal base mismatch discrimination. The oligonucleotide probes were tethered to the glass using a simple, improved attachment chemistry in which a 3'-aminopropanol function introduced into the probe during chemical synthesis binds covalently to silanol groups on clean, underivatized glass. "Operating parameters" examined in the stacking hybridization system included length of capture probe, position, type and number of mismatches between the probe and the target, temperature of hybridization and length of washing, and the presence of terminal phosphate group in the probe, at its junction with the stacking oligomer. The results suggest that in the stacking hybridization configuration: 1. Optimal mismatch discrimination with 9-mer probes occurs at 45 degrees C, after which little or no improvement in mispair rejection occurred on lengthy continued washing at 45 degrees C. 2. At 25 degrees C optimal mismatch discrimination occurred with 7- or 8-mer probes, or with 9-mer probes containing an additional internal mismatch. 3. The presence of a phosphate group on the 5'-end of the glass-tethered probe had no general effect on mismatch discrimination, but influenced the relative stability of different mismatches in the sequence context studied. These results provide a motivation for continued development of the stacking hybridization technique for nucleic acid sequence analysis. This approach offers several advantages over the traditional allele-specific oligonucleotide hybridization technique, and is distinct from the contiguous stacking hybridization sitrategy that the Mirzabekov laboratory has introduced (Yershov et al. (1996) Proc. Natl. Acad. Sci. USA 93, 4913-4918; Parinov et al. (1996) Nucleic Acids Res. 24, 2998-3004).     

DNA hybridization in nanostructural molecular assemblies enables detection of gene mutations without a fluorescent probe

We have developed a simple single nucleotide polymorphisms (SNPs) analysis utilizing DNA hybridization in nanostructural molecular assemblies. The novel technique enables the detection of a single-base mismatch in a DNA sequence without a fluorescent probe. This report describes for the first time that DNA hybridization occurs in the nanostructural molecular assemblies (termed reverse micelles) formed in an organic medium. The restricted nanospace in the reverse micelles amplifies the differences in the hybridization rate between mismatched and perfectly matched DNA probes. For a model system, we hybridized a 20-mer based on the p53 gene sequence to 20-mer complementary oligonucleotides with various types of mismatches. Without any DNA labeling or electrochemical apparatus, we successfully detected the various oligonucleotide mismatches by simply measuring the UV absorbance at 260 nm.     

LNA guanine and 2,6-diaminopurine. Synthesis, characterization and hybridization properties of LNA 2,6-diaminopurine containing oligonucleotides

LNA guanine and 2,6-diaminopurine (D) phosphoramidites have been synthesized as building blocks for antisense oligonucleotides (ON). The effects of incorporating LNA D into ON were investigated. As expected, LNA D containing ON showed increased affinity towards complementary DNA (Delta Tm +1.6 to +3.0 degrees C) and RNA (Delta Tm +2.6 to +4.6 degrees C) ON. To evaluate if LNA D containing ON have an enhanced mismatch sensitivity compared to their complementary LNA A containing ON thermal denaturation experiments towards singly mismatched DNA and RNA ON were undertaken. Replacing one LNA A residue with LNA D, in fully LNA modified ON, resulted in higher mismatch sensitivity towards DNA ON (Delta Delta Tm -4 to >-17 degrees C). The same trend was observed towards singly mismatched RNA ON (Delta Delta Tm D-a = -8.7 degrees C and D-g = -4.5 degrees C) however, the effect was less clearcut and LNA A showed a better mismatch sensitivity than LNA D towards cytosine (Delta Tm +5.5 degrees C).     

Thermodynamic characterization of naturally occurring RNA single mismatches with G-U nearest neighbors

Due to their prevalence and roles in biological systems, single mismatches adjacent to G-U pairs are important RNA structural elements. Since there are only limited experimental values for the stability of single mismatches adjacent to G-U pairs, current algorithms using free energy minimization to predict RNA secondary structure from sequence assign predicted thermodynamic values to these types of single mismatches. Here, thermodynamic data are reported for frequently occurring single mismatches adjacent to at least one G-U pair. This experimental data can be used in place of predicted thermodynamic values in algorithms that predict secondary structure from sequence using free energy minimization. When predicting the thermodynamic contributions of previously unmeasured single mismatches, most algorithms apply the same thermodynamic penalty for an A-U pair adjacent to a single mismatch and a G-U pair adjacent to a single mismatch. A recent study, however, suggests that the penalty for a G-U pair adjacent to a tandem mismatch should be 1.2 +/- 0.1 kcal/mol, and the penalty for an A-U pair adjacent to a tandem mismatch should be 0.5 +/- 0.2 kcal/mol [Christiansen, M. E. and Znosko, B. M. (2008) Biochemistry 47, 4329-4336]. Therefore, the data reported here are combined with the existing thermodynamic dataset of single mismatches, and nearest neighbor parameters are derived for an A-U pair adjacent to a single mismatch (1.1 +/- 0.1 kcal/mol) and a G-U pair adjacent to a single mismatch (1.4 +/- 0.1 kcal/mol).     

Single nucleotide polymorphism discrimination assisted by improved base stacking hybridization using oligonucleotide microarrays

Efficiencies of mismatch discrimination using size-varied capture probes were examined at various hybridization temperatures. The probes were 17, 15, 13, 11, 9, and 7 nucleotides long and contained single-base mismatches at their 3' ends. The optimal signal intensity and efficiency of base stacking hybridization on mismatch discrimination were observed for capture probes with a melting temperature (Tm) value of 36 degrees C, in the detection of DNA sequence variations at 40 degrees C. We employed asymmetric PCR to prepare single-stranded target DNA labeled with a fluorescent dye, and the PCR product was hybridized on the DNA microarray with no further purification. Our efforts have enhanced the sensitivity and simplified the procedures of base stacking hybridization on mismatch discrimination. As a model experiment, this improved technology was used to identify plasmid templates of human leukocyte antigen (HLA)-A alleles 2601, 2902, and 0206 on oligonucleotide microarrays. It is now possible to apply this simple, rapid, sensitive, and reliable base stacking hybridization technology to detect DNA sequence variations on microarrays in clinical diagnosis and other applications.     

Recognition of base mismatches in DNA by 5,6-chrysenequinone diimine complexes of rhodium(III): a proposed mechanism for preferential binding in destabilized regions of the double helix

5,6-chrysenequinone diimine (chrysi) complexes of rhodium(III) have been shown to be versatile and specific recognition agents for mismatched base pairs in DNA. The design of these compounds was based on the hypothesis that the sterically expansive chrysi ligand, which should be too wide to readily intercalate into B-DNA, would bind preferentially in the destabilized regions of the DNA helix near base mismatches. In this work, this recognition hypothesis is comprehensively explored. Comparison of the recognition patterns of the complex [Rh(bpy)(2)(chrysi)](3+) with a nonsterically demanding analogue, [Rh(bpy)(2)(phi)](3+) (phi = 9,10-phenanthrenequinone diimine), demonstrates that the chrysi ligand does indeed disfavor binding to B-DNA and generate mismatch selectivity. Examination of mismatch recognition by [Rh(bpy)(2)(chrysi)](3+) in both constant and variable sequence contexts using photocleavage assays indicates that the recognition of base mismatches is influenced by the amount that a mismatch thermodynamically destabilizes the DNA helix. Thermodynamic binding constants for the rhodium complex at a range of mismatch sites have been determined by quantitative photocleavage titration and yield values which vary from 1 x 10(6) to 20 x 10(6) M(-)(1). These mismatch-specific binding affinities correlate with independent measurements of thermodynamic destabilization, supporting the hypothesis that helix destabilization is a factor determining the binding affinity of the metal complex for the mismatched site. Although not the only factor involved in the binding of [Rh(bpy)(2)(chrysi)](3+) to mismatch sites, a model is proposed where helix destabilization acts as the "door" which permits access of the sterically demanding intercalator to the base stack.     

Thermodynamics of single mismatches in RNA duplexes

The thermodynamic properties and structures of single mismatches in short RNA duplexes were studied in optical melting and imino proton NMR experiments. The free energy increments at 37 degrees C measured for non-GU single mismatches range from -2.6 to 1.7 kcal/mol. These increments depend on the identity of the mismatch, adjacent base pairs, and the position in the helix. UU and AA mismatches are more stable close to a helix end, but GG mismatch stability is essentially unaffected by the position in the helix. Approximations are suggested for predicting stabilities of single mismatches in short RNA duplexes.     

Influence of nearest neighbor sequence on the stability of base pair mismatches in long DNA; determination by temperature-gradient gel electrophoresis.

Temperature-gradient gel electrophoresis (TGGE) was employed to determine the thermal stabilities of 48 DNA fragments that differ by single base pair mismatches. The approach provides a rapid way for studying how specific base mismatches effect the stability of a long DNA fragment. Homologous 373 bp DNA fragments differing by single base pair substitutions in their first melting domain were employed. Heteroduplexes were formed by melting and reannealing pairs of DNAs, one of which was 32P-labeled on its 5'-end. Product DNAs were separated based on their thermal stability by parallel and perpendicular temperature-gradient gel electrophoresis. The order of stability was determined for all common base pairs and mismatched bases in four different nearest neighbor environments; d(GXT).d(AYC), d(GXG).d(CYC), d(CXA).d(TYG), and d(TXT).d(AYA) with X,Y = A, T, C, or G. DNA fragments containing a single mismatch were destabilized by 1 to 5 degrees C with respect to homologous DNAs with complete Watson-Crick base pairing. Both the bases at the mismatch site and neighboring stacking interactions influence the destabilization caused by a mismatch. G.T, G.G and G.A mismatches were always among the most stable mismatches for all nearest neighbor environments examined. Purine.purine mismatches were generally more stable than pyrimidine.pyrimidine mispairs. Our results are in very good agreement with data where available from solution studies of short DNA oligomers.     

Thermodynamics of DNA duplexes with adjacent G.A mismatches.

The sequence 5'-d(ATGAGCGAAT) forms a very stable self-complementary duplex with four G.A mismatch base pairs (underlined) out of ten total base pairs [Li et al. (1991) Proc. Natl. Acad. Sci. U.S.A. 88, 26-30]. The conformation is in the general B-family and is stabilized by base-pair hydrogen bonding of an unusual type, by favorable base dipole orientations, and by extensive purine-purine stacking at the mismatched sites. We have synthesized 13 decamers with systematic variations in the sequence above to determine how the flanking sequences, the number of G.A mismatches, and the mismatch sequence order (5'-GA-3' or 5'-AG-3') affect the duplex stability. Changing A.T to G.C base pairs in sequences flanking the mismatches stabilizes the duplexes, but only to the extent observed with B-form DNA. The sequence 5'-pyrimidine-GA-purine-3', however, is considerably more stable than 5'-purine-GA-pyrimidine-3'. The most stable sequences with two pairs of adjacent G.A mismatches have thermodynamic parameters for duplex formation that are comparable to those for fully Watson-Crick base-paired duplexes. Similar sequences with single G.A pairs are much less stable than sequences with adjacent G.A mismatches. Reversing the mismatch order from 5'-GA-3' to 5'-AG-3' results in an oligomer that does not form a duplex. These results agree with predictions from the model derived from NMR and molecular mechanics and indicate that the sequence 5'-pyrimidine-GA-purine-3' forms a stable conformational unit that fits quite well into a B-form double helix.(ABSTRACT TRUNCATED AT 250 WORDS)     

Hybridization of synthetic oligodeoxyribonucleotides to phi chi 174 DNA: the effect of single base pair mismatch.

Oligodeoxyribonucleotides complementary to the DNA of the wild type (wt) bacteriophage phi chi 174 have been synthesized by the phosphotriester method. The oligomers, 11, 14, and 17 bases long, are complementary to the region of the DNA which accounts for the am-3 point mutation. When hybridized to am-3 DNA, the oligonucleotides form duplexes with a single base pair mismatch. The thermal stability of the duplexes formed between wt and am-3 DNAs has been measured. The am-3 DNA:oligomer duplexes dissociate at a temperature about 10 degrees C lower than the corresponding wt DNA:oligomer duplexes. This dramatic decrease in thermal stability due to a single mismatch makes it possible to eliminate the formation of the mismatched duplexes by the appropriate choice of hybridization temperature. These results are discussed with respect to the use of oligonucleotides as probes for the isolation of specific cloned DNA sequences.     

Structure, dynamics, and thermodynamics of mismatched DNA oligonucleotide duplexes d(CCCAGGG)2 and d(CCCTGGG)2

The structures and hydrogen exchange properties of the mismatched DNA oligonucleotide duplexes d(CCCAGGG)2 and d(CCCTGGG)2 have been studied by high-resolution nuclear magnetic resonance. Both the adenine-adenine and thymine-thymine mismatches are intercalated in the duplexes. The structures of these self-complementary duplexes are symmetric, with the two strands in equivalent positions. The evidence indicates that these mismatches are not stably hydrogen bonded. The mismatched bases in both duplexes are in the anti conformation. The mismatched thymine nucleotide in d(CCCTGGG)2 is intercalated in the duplex with very little distortion of the bases or sugar-phosphate backbone. In contrast, the bases of the adenine-adenine mismatch in d(CCCAGGG)2 must tilt and push apart to reduce the overlap of the amino groups. The thermodynamic data show that the T-T mismatch is less destabilizing than the A-A mismatch when flanked by C-G base pairs in this sequence, in contrast to their approximately equal stabilities when flanked by A-T base pairs in the sequence d(CAAAXAAAG.CTTTYTTTG) where X and Y = A, C, G, and T [Aboul-ela, F., Koh, D., & Tinoco, I., Jr. (1985) Nucleic Acids Res. 13, 4811]. Although the mechanism cannot be determined conclusively from the limited data obtained, exchange of the imino protons with solvent in these destabilized heteroduplexes appears to occur by a cooperative mechanism in which half the helix dissociates.     

Conformational and thermodynamic characterization of the molten globule state occurring during unfolding of cytochromes-c by weak salt denaturants

The denaturation of bovine and horse cytochromes-c by weak salt denaturants (LiCl and CaCl(2)) was measured at 25 degrees C by observing changes in molar absorbance at 400 nm (Delta epsilon(400)) and circular dichroism (CD) at 222 and 409 nm. Measurements of Delta epsilon(400) and mean residue ellipticity at 409 nm ([theta](409)) gave a biphasic transition for both modes of denaturation of cytochromes-c. It has been observed that the first denaturation phase, N (native) conformation <--> X (intermediate) conformation and the second denaturation phase, X conformation <--> D (denatured) conformation are reversible. Conformational characterization of the X state by the far-UV CD, 8-anilino-1-naphthalene sulfonic acid (ANS) binding, and intrinsic viscosity measurements led us to conclude that the X state is a molten globule state. Analysis of denaturation transition curves for the stability of different states in terms of Gibbs energy change at pH 6.0 and 25 degrees C led us to conclude that the N state is more stable than the X state by 9.55 +/- 0.32 kcal mol(-1), whereas the X state is more stable than the D state by only 1.40 +/- 0.25 kcal mol(-1). We have also studied the effect of temperature on the equilibria, N conformation <--> X conformation and X conformation <--> D conformation in the presence of different denaturant concentrations using two different optical probes, namely, [theta](222) and Delta epsilon(400). These measurements yielded T(m), (midpoint of denaturation) and Delta H(m) (enthalpy change) at T(m) as a function of denaturant concentration. A plot of Delta H(m) versus corresponding T(m) was used to determine the constant-pressure heat capacity change, Delta C(p) (= ( partial differential Delta H(m)/ partial differential T(m))(p)). Values of Delta C(p) for N conformation <--> X conformation and X conformation <--> D conformation is 0.92 +/- 0.02 kcal mol(-1) K(-1) and 0.41 +/- 0.01 kcal mol(-1) K(-1), respectively. These measurements suggested that about 30% of the hydrophobic groups in the molten globule state are not accessible to the water.     

Thermodynamic basis of selectivity in guide-target-mismatched RNA interference

Silencing in RNAi is strongly affected by guide‐strand/target‐mRNA mismatches. Target nucleation is thought to occur at positions 2–8 of the guide (“seed region”); successful hybridization in this region is the primary determinant of target‐binding affinity and hence target cleavage. To define a molecular basis for the target sequence selectivity in RNAi, we studied all possible distinct single mismatches in seven positions of the seed region—a total of 21 substitutions. We report results from soft‐core thermodynamic integration simulations to determine changes in targeting binding‐free energies to Argonaute due to single mismatches in the guide strand, which arise during binding of an imperfectly matched target mRNA. In agreement with experiment, most mismatches impair target binding, consistent with a prominent role for binding affinity changes in RNAi sequence selectivity. Individual Argonaute residues located near the mismatched base pair are found to contribute significantly to binding affinity changes. We also use this methodology to analyze the mismatch‐dependent free energy changes for dissociation of a DNA?RNA hybrid from Argonaute, as a model for the escape of miRNAs from the silencing pathway. Several mismatched sequences of the miRNA have increased affinity to Argonaute, implying that some mismatches may reduce the probability for escape. Furthermore, calculations of base‐substitution‐dependent free energy changes for binding ssDNA reveal mild sequence sensitivity as expected for guide strand binding to Argonaute. Our findings give a thermodynamic basis for RNAi target sequence selectivity and suggest that miRNA mismatches may increase silencing effectiveness and thus could be evolutionarily advantageous. Proteins 2012; © 2011 Wiley Periodicals, Inc.     

Thermodynamic characterization of single mismatches found in naturally occurring RNA

Many naturally occurring RNA structures contain single mismatches. However, the algorithms currently used to predict RNA structure from sequence rely on a minimal set of data for single mismatches, most of which occur rather infrequently in nature. As a result, several approximations and assumptions are used to predict the stability of RNA duplexes containing the most common single mismatches. Therefore, the relative frequency of single mismatches was determined by compiling and searching a database of 955 RNA secondary structures. Thermodynamic parameters for duplex formation, derived from optical melting experiments, are reported for 28 oligoribonucleotides containing frequently occurring single mismatches. These data were then combined with previous data to construct a dataset of 64 single mismatches, including the 30 most common in the database. Because of this increase in experimental thermodynamic parameters for single mismatches that occur frequently in nature, more accurate free energy calculations have resulted. To improve the prediction of the thermodynamic parameters for duplexes containing single mismatches that have not been experimentally measured, single mismatch-specific nearest neighbor parameters were derived. The free energy of an RNA duplex containing a single mismatch that has not been thermodynamically characterized can be calculated by: DeltaG degrees 37,single mismatch = DeltaG degrees 37,mismatch nt + DeltaG degrees 37,mismatch-NN interaction + DeltaG degrees 37,AU/GU. Here, DeltaG degrees 37,mismatch is -0.4, -2.1, and -0.3 kcal/mol for A.G, G.G, and U.U mismatches, respectively; DeltaG degrees 37,mismatch-NN interaction is 0.7, -0.5, 0.4, -0.4, and -1.0 kcal/mol for 5'YRR3'/3'RRY5', 5'RYY3'/3'YYR5', 5'YYR3'/3'RYY5', 5'YRY3'/3'RYR5', and 5'RRY3'/3'YYR5' mismatch-nearest neighbor combinations, respectively, when A and G are categorized as purines (R) and C and U are categorized as pyrimidines (Y); and DeltaG degrees 37,AU/GU is a penalty of 1.2 kcal/mol for replacing a G-C base pair with either an A-U or G-U base pair. Similar predictive models were also derived for DeltaH degrees single mismatch and DeltaS degrees single mismatch. These new predictive models, in conjunction with the reported thermodynamics for frequently occurring single mismatches, should allow for more accurate calculations of the free energy of RNA duplexes containing single mismatches and, furthermore, allow for improved prediction of secondary structure from sequence.