TaqMan probes based on oligonucleotide-hairpin minor groove binder conjugates

Hybridization of TaqMan probes derived from oligonucleotides containing fluorophores (fluorescein, FAM, or tetramethylrhodamine, (Tamra)), fluorescence quenchers (BHQ1 or BHQ2), and a conjugated hairpin binder (MGB) composed of two tripyrrolcarboxamide residues connected through an aminobutyric acid residue were proposed for discrimination of single base mismatch using the real time PCR technique. Identification of A/C mismatch was shown to be highly specific for hepatitis C virus subtypes 1a and 1b with two variants of the probe (5′-3′): Tamra-ATTGAGCGGGTTTAp-BHQ2-MGB for subtype 1a and FAMATTGAGCGGGTTGAp-BHQ1-MGB for subtype 1b. Perfect duplexes (A·T-and G·C pairs) increase fluorescence in the process of amplification, whereas imperfect duplexes (A·G-and T·C pairs) induce no fluorescence changes. This phenomenon enables simultaneous genotyping of hepatitis C virus subtypes 1a and 1b.     

Association of hepatitis C virus NS5B variants with resistance to new antiviral drugs among untreated patients

Mutations located in the 109-amino acid fragment of NS5B are typically associated with resistance to interferon (IFN) and ribavirin (RIB) and to new antiviral drugs. The prevalence of these mutations was examined in 69 drug-na?ve individuals with hepatitis C virus (HCV) infections in Rio de Janeiro, Brazil. Mutations related to non-response to IFN/RIB were observed in all subtypes studied (1a, 1b, 2b, 3a and 4). The most common mutation was Q309R, present in all subtypes, except subtype 2b with frequency above 20%. D244N was detected only in subtype 3a and A333E was detected only in subtype 2b. We did not detect the S282T, S326G or T329I mutations in any of the samples analysed. Of note, the C316N mutation, previously related to a new non-nucleoside compound (HCV796 and AG-021541), was observed in only eight of 33 (24%) samples from subtype 1b. Site 316 was under positive selection in this HCV variant. Our data highlight the presence of previously described resistance mutations in HCV genotypes from drug-na?ve patients.     

3'-minor groove binder-DNA probes increase sequence specificity at PCR extension temperatures

DNA probes with conjugated minor groove binder (MGB) groups form extremely stable duplexes with single-stranded DNA targets, allowing shorter probes to be used for hybridization based assays. In this paper, sequence specificity of 3′-MGB probes was explored. In comparison with unmodified DNA, MGB probes had higher melting temperature (T_m) and increased specificity, especially when a mismatch was in the MGB region of the duplex. To exploit these properties, fluorogenic MGB probes were prepared and investigated in the 5′-nuclease PCR assay (real-time PCR assay, TaqMan assay). A 12mer MGB probe had the same T_m (65°C) as a no-MGB 27mer probe. The fluorogenic MGB probes were more specific for single base mismatches and fluorescence quenching was more efficient, giving increased sensitivity. A/T rich duplexes were stabilized more than G/C rich duplexes, thereby leveling probe T_m and simplifying design. In summary, MGB probes were more sequence specific than standard DNA probes, especially for single base mismatches at elevated hybridization temperatures.     

Point mutation in the S gene of hepatitis B virus for a d/y or w/r subtypic change in two blood donors carrying a surface antigen of compound subtype adyr or adwr

Genomes of hepatitis B virus (HBV) were cloned from the plasma of a blood donor who carried subviral particles of three distinct subtypes in the following proportions: adr, 25%; ayr, 63%; and adyr, 12%. HBV DNA clones were classified into two groups based on a difference at only one nucleotide in the S gene. Two clones had A as nucleotide 365 that formed part of the codon for lysine as amino acid residue 122 and produced a surface antigen of subtype adr in transfected NIH 3T3 cells. The remaining four clones had G determining the codon for arginine and produced a surface antigen of subtype ayr in transfected cells. Similarly, HBV genomes were cloned from the plasma of an individual who carried subviral particles of subtypes adr (71%) and adwr (29%). Two clones had T and A as nucleotides 476 and 479, respectively. The other seven clones had C and G as the respective nucleotides. Based on a comparison with previously reported HBV genomes of various subtypes, the mutation of nucleotide 479, forming part of the codon for lysine or arginine as amino acid residue 160, was deduced to determine the w or r subtype, respectively. When NIH 3T3 cells were transfected separately with the genome of subtype adw or adr, derived from plasma containing a surface antigen of subtype adwr, and then cocultured, they produced subviral particles of either subtype adw or adr. When cells were transfected with the genomes of subtypes adw and adr simultaneously, however, subviral particles were produced that possessed w and r determinants on the selfsame particles. These results attributed the d/y or w/r subtypic change to a point mutation in the S gene and favored coinfection of hepatocytes with an HBV genome and its mutant as the mechanism of compound subtypes.     

Structural characterisation of the bromouracil.guanine base pair mismatch in a Z-DNA fragment.

The deoxyoligonucleotide d(BrU-G-C-G-C-G) was crystallised at pH 8.2 and its structure analysed by X-ray diffraction. The unit cell, of dimensions a = 17.94, b = 30.85, c = 49.94A contains four DNA duplexes in space group P2(1)2(1)2(1). The duplexes are in the Z conformation, with four Watson-Crick G.C base pairs and two BrU.G base pairs. The structure was refined to an R factor of 0.16 at a resolution of 2.2A with 64 solvent molecules located. The BrU.G base pair mismatch is of the wobble type, with both bases in the major tautomer form and hydrogen bonds linking 0-2 of BrU with N-1 of G and N3 of BrU with 0-6 of G. There is no indication of the presence of ionised base pairs, in spite of the high pH of crystallisation. The results are discussed in terms of the mutagenic properties of 5- bromouracil.     

Determination of affinity,stoichiometry and sequence selectivity of minor groove binder complexes with double-stranded oligodeoxynucleotides by electrospray ionization mass spectrometry

Electrospray mass spectrometry was evaluated regarding the reliability of the determination of the stoichiometries and equilibrium association constants from single spectra. Complexes between minor groove binders (Hoechst 33258, Hoechst 33342, DAPI, netropsin and berenil) and 12mer oligonucleotide duplexes with a central sequence (A/T)₄ flanked by G/C base pairs were chosen as model systems. To validate the electrospray ionization mass spectrometry (ESI-MS) method, comparisons were made with circular dichroism and fluorescence spectroscopy measurements. ESI-MS allowed the detection of minor (2 drug + DNA) species for Hoechst 33258, Hoechst 33342, DAPI and berenil with duplex d(GGGG(A/T)₄GGGG)· d(CCCC(A/T)₄CCCC), which were undetectable with the other techniques. Assuming that the duplexes and the complexes have the same electrospray response factors, the equilbrium association constants of the 1:1 and 2:1 complexes were determined by ESI-MS, and the values show a good quantitative agreement with fluorescence determined constants for Hoechst 33258 and Hoechst 33342. It is also shown that ESI-MS can quickly give reliable information on the A/T sequence selectivity of a drug: the signal of a complex is directly related to the affinity of the drug for that particular duplex. The potential of ESI-MS as a qualitative and quantitative affinity screening method is emphasized.     

Unusual duplex formation in purine rich oligodeoxyribonucleotides

The purine rich oligodeoxyribonucleotides 1C, d(ATGACGGAATA) and 2C, d(ATGAGCGAATA) alone exhibit highly cooperative melting transitions. Analysis of the concentration dependence of melting, and electrophoretic studies indicate that these oligomers can form an unusual purine rich offset double helix. The unusual duplex is predicted to contain four A.T, two G.C, and four G.A mismatch base pairs as well as a single A base stacked on the 3' end of each chain of the helix. Other possible models for the duplex are unlikely because they are predicted to contain many base pairs of low stability. Changing the central sequence to CGG or GGG should destabilize the duplex and this is observed. The unusual duplex of 2C is more stable than the duplex of 1C indicating that the stability of G.A base pairs is quite sensitive to the surrounding sequence. Addition of 1C and 2C to their complementary pyrimidine strands results in normal duplexes of similar stability. We feel that the unusual duplexes are significantly stabilized by the intrinsic stacking tendency of purine bases.     

Parallel-stranded duplex DNA containing dA · dU base pairs

DNA oligonucleotides with dA and dU residues can form duplexes with trans d(A · U) base pairing and the sugar-phosphate backbone in a parallel-stranded orientation, as previously established for oligonucleotides with d(A · T) base pairs. The properties of such parallel-stranded DNA (ps-DNA) 25-mer duplexes have been characterized by absorption (uv), CD, ir, and fluorescence spectroscopy, as well as by nuclease sensitivity. Comparisons were made with duplex molecules containing (a) dT in both strands, (b) dU in one strand and dT in the second, and (c) the same base combinations in reference antiparallel-stranded (aps) structures. Thermodynamic analysis revealed that total replacement of deoxythymine by deoxyuridine was accompanied by destabilization of the ps-helix (reduction in T_m by −13°C in 2 mM MgGl₂, 10 mM Na-cacodylate). The U-containing ps-helix (U1 · U2) also melted 14°C lower than the corresponding aps-helix under the same ionic conditions; this difference was very close to that observed between ps and aps duplexes with d(A · T) base pairs. Force field minimized structures of the various ps and aps duplexes with either d(A · T) or d(A · U) base pairs ps/aps and dT/dU combinations are presented. The energy-minimized helical parameters did not differ significantly between the DNAs containing dT and dU. © 1996 John Wiley & Sons, Inc.     

Subtypes of human immunodeficiency virus type 1 and disease stage among women in Nairobi, Kenya

In sub-Saharan Africa, where the effects of human immunodeficiency virus type 1 (HIV-1) have been most devastating, there are multiple subtypes of this virus. The distribution of different subtypes within African populations is generally not linked to particular risk behaviors. Thus, Africa is an ideal setting in which to examine the diversity and mixing of viruses from different subtypes on a population basis. In this setting, it is also possible to address whether infection with a particular subtype is associated with differences in disease stage. To address these questions, we analyzed the HIV-1 subtype, plasma viral loads, and CD4 lymphocyte levels in 320 women from Nairobi, Kenya. Subtype was determined by a combination of heteroduplex mobility assays and sequence analyses of envelope genes, using geographically diverse subtype reference sequences as well as envelope sequences of known subtype from Kenya. The distribution of subtypes in this population was as follows: subtype A, 225 (70.3%); subtype D, 65 (20.5%); subtype C, 22 (6.9%); and subtype G, 1 (0.3%). Intersubtype recombinant envelope genes were detected in 2.2% of the sequences analyzed. Given that the sequences analyzed represented only a small fraction of the proviral genome, this suggests that intersubtype recombinant viral genomes may be very common in Kenya and in other parts of Africa where there are multiple subtypes. The plasma viral RNA levels were highest in women infected with subtype C virus, and women infected with subtype C virus had significantly lower CD4 lymphocyte levels than women infected with the other subtypes. Together, these data suggest that women in Kenya who are infected with subtype C viruses are at more advanced stages of immunosuppression than women infected with subtype A or D. There are at least two models to explain the data from this cross-sectional study; one is that infection with subtype C is associated with a more rapid disease progression, and the second is that subtype C represents an older epidemic in Kenya. Discriminating between these possibilities in a longitudinal study will be important for increasing our understanding of the role of specific subtypes in the transmission and pathogenesis of HIV-1.     

Thermodynamics of RNA duplexes with tandem mismatches containing a uracil-uracil pair flanked by C.G/G.C or G.C/A.U closing base pairs

The thermodynamics governing the denaturation of RNA duplexes containing 8 bp and a central tandem mismatch or 10 bp were evaluated using UV absorbance melting curves. Each of the eight tandem mismatches that were examined had one U-U pair adjacent to another noncanonical base pair. They were examined in two different RNA duplex environments, one with the tandem mismatch closed by G.C base pairs and the other with G.C and A.U closing base pairs. The free energy increments (Delta Gdegrees(loop)) of the 2 x 2 loops were positive, and showed relatively small differences between the two closing base pair environments. Assuming temperature-independent enthalpy changes for the transitions, (Delta Gdegrees(loop)) for the 2 x 2 loops varied from 0.9 to 1.9 kcal/mol in 1 M Na(+) at 37 degrees C. Most values were within 0.8 kcal/mol of previously estimated values; however, a few sequences differed by 1.2-2.0 kcal/mol. Single strands employed to form the RNA duplexes exhibited small noncooperative absorbance increases with temperature or transitions indicative of partial self-complementary duplexes. One strand formed a partial self-complementary duplex that was more stable than the tandem mismatch duplexes it formed. Transitions of the RNA duplexes were analyzed using equations that included the coupled equilibrium of self-complementary duplex and non-self-complementary duplex denaturation. The average heat capacity change (DeltaC(p)) associated with the transitions of two RNA duplexes was estimated by plotting DeltaH degrees and DeltaS degrees evaluated at different strand concentrations as a function of T(m) and ln T(m), respectively. The average DeltaC(p) was 70 +/- 5 cal K(-)(1) (mol of base pairs)(-)(1). Consideration of this heat capacity change reduced the free energy of formation at 37 degrees C of the 10 bp control RNA duplexes by 0.3-0.6 kcal/mol, which may increase Delta Gdegrees(loop) values by similar amounts.     

Frequent Detection of Hepatitis C Virus Subtype 3a (HCV-3a) Isolates in Thailand by PCR Using Subtype-Specific Primers

By means of a polymerase chain reaction (PCR) method using subtype-specific primers for hepatitis C virus (HCV) subtypes 1a, 1b, 2a, 2b and 3a, the prevalence of each subtype among HCV isolates in Chiang Mai, Thailand, was determined. HCV-3a appeared to be the most common subtype in blood donors, and was also frequently found in patients with liver disease. HCV-1b, but not HCV-2a or ?2b, was also commonly found in this area, while a considerable percentage of the total HCV isolates still remained unclassifiable by the above methods. Serotype analysis of the HCV isolates using C14-1 and C14-2 recombinant peptides revealed that HCV-3a was likely to carry an antigenic determinant(s) different from those of the major types 1 (HCV-1a and ?1b) and 2 (HCV-2a and ?2b).     

Mismatches in DNA double strands: thermodynamic parameters and their correlation to repair efficiencies.

The helix-coil transitions of the 16 octadecameric DNA duplexes dCGTCGTTTXACAACGTCG X dCGACGTTGTX1AAACGACG with A, T, G, and C for X and X1 were measured by UV-absorption. This sequence was taken from former studies of in vivo determination of efficiencies of mismatch repair (Kramer, Kramer, and Fritz (1984) Cell 38, 879-887). The thermodynamic parameters for double strand and mismatch formation have been obtained by evaluating the partition function of a stack model which allowed for loop formation. As a result the mismatches could be classified into wobble base pairs (T/G, G/G, C/A, A/A, A/G), open base pairs, i.e. permanent loops (T/T, C/T, T/C, C/C), and intermediate or weak base pairs (G/T, A/C, G/A). There is no correlation between Tm and the biological repair efficiency of X/X1. The structure classes, however, as described above show a close correlation: Open base pairs show the lowest repair efficiencies, whereas mismatches with high repair efficiency always belong to the structural class of wobble base pairs. Because of the palindromic nearest neighbors of the variation site X/X1, the influence of next-nearest neighbor interactions could be detected and be estimated to about 1 kJ/mol for one stack.     

Base pairing involving deoxyinosine: implications for probe design.

The thermal stability of oligodeoxyribonucleotide duplexes containing deoxyinosine (I) residues matched with each of the four normal DNA bases were determined by optical melting techniques. The duplexes containing at least one I were obtained by mixing equimolar amounts of an oligonucleotide of sequence dCA3XA3G with one of sequence dCT3YT3G where X and Y were A, C, G, T, or I. Comparison of optical melting curves yielded relative stabilities for the I-containing standard base pairs in an otherwise identical base-pair sequence. I:C pairs were found to be less stable than A:T pairs in these duplexes. Large neighboring-base effects upon stability were observed. For example, when (X,Y) = (I,A), the duplex is eight-fold more stable than when (X,Y) = (A,I). Independent of sequence effects the order of stabilities is: I:C greater than I:A greater than I:T congruent to I:G. This order differs from that of deoxyguanosine which pairs less strongly with dA; otherwise each deoxyinosine base pair is less stable than its deoxyguanosine counterpart in the same sequence environment. Implications of these results for design of DNA oligonucleotide probes are discussed.     

Consideration of Viral Resistance for Optimization of Direct Antiviral Therapy of Hepatitis C Virus Genotype 1-Infected Patients

Different highly effective interferon-free treatment options for chronic hepatitis C virus (HCV) infection are currently available. Pre-existence of resistance associated variants (RAVs) to direct antiviral agents (DAAs) reduces sustained virologic response (SVR) rates by 3–53% in hepatitis C virus (HCV) genotype 1 infected patients depending on different predictors and the DAA regimen used. Frequencies of single and combined resistance to NS3, NS5A and NS5B inhibitors and consequences for the applicability of different treatment regimens are unknown. Parallel population based sequencing of HCV NS3, NS5A and NS5B genes in 312 treatment-naïve Caucasian HCV genotype 1 infected patients showed the presence of major resistant variants in 20.5% (NS3), 11.9% (NS5A), and 22.1% (NS5B) with important differences for HCV subtypes. In NS3, Q80K was observed in 34.7% and 2.1% of subtype 1a and 1b patients, respectively while other RAVs to second generation protease inhibitors were detected rarely (1.4%). Within NS5A RAVs were observed in 7.1% of subtype 1a and 17.6% in subtype 1b infected patients. RAVs to non-nucleoside NS5B inhibitors were observed in 3.5% and 44.4% of subtype 1a and 1b patients, respectively. Considering all three DAA targets all subtype 1a and 98.6% of subtype 1b infected patients were wildtype for at least one interferon free DAA regimen currently available. In conclusion, baseline resistance testing allows the selection of at least one RAVs-free treatment option for nearly all patients enabling a potentially cost- and efficacy-optimized treatment of chronic hepatitis C.     

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.     

Identification of hepatitis C virus subtype 2c by sequencing analysis in patients from Córdoba, Argentina

In Argentina, most information on hepatitis C virus (HCV) genotype distribution comes from studies carried out in Buenos Aires (east province). In order to identify HCV subtypes in central Argentina, nucleotide sequencing of core region was performed in samples from 36 patients living in Córdoba, the second most populated province of Argentina. The sequence analysis identified subtype 2c as the most prevalent (50%), followed by subtype 1b (33%) and to a lesser extent by subtypes 1a (11%), 3a (3%) and 4a (3%). This is the first report of circulation of HCV subtype 2c in this region of Argentina and also such high prevalence has never been found before in the genotype distribution of South America.     

Nearest-neighbor thermodynamics of deoxyinosine pairs in DNA duplexes

Nearest-neighbor thermodynamic parameters of the ‘universal pairing base’ deoxyinosine were determined for the pairs I·C, I·A, I·T, I·G and I·I adjacent to G·C and A·T pairs. Ultraviolet absorbance melting curves were measured and non-linear regression performed on 84 oligonucleotide duplexes with 9 or 12 bp lengths. These data were combined with data for 13 inosine containing duplexes from the literature. Multiple linear regression was used to solve for the 32 nearest-neighbor unknowns. The parameters predict the T_m for all sequences within 1.2°C on average. The general trend in decreasing stability is I·C > I·A > I·T ≈ I· G > I·I. The stability trend for the base pair 5′ of the I·X pair is G·C > C·G > A·T > T·A. The stability trend for the base pair 3′ of I·X is the same. These trends indicate a complex interplay between H-bonding, nearest-neighbor stacking, and mismatch geometry. A survey of 14 tandem inosine pairs and 8 tandem self-complementary inosine pairs is also provided. These results may be used in the design of degenerate PCR primers and for degenerate microarray probes.     

Early evolution of the human immunodeficiency virus type 1 subtype C epidemic in rural Malawi

We have tracked the early years of the evolution of the human immunodeficiency virus type 1 (HIV-1) epidemic in a rural district of central east Africa from the first documented introductions of subtypes A, D, and C to the present predominance of subtype C. The earliest subtype C sequences ever reported are described. Blood samples were collected on filter papers from 1981 to 1984 and from 1987 to 1989 from more than 44,000 individuals living in two areas of Karonga District, Malawi. These samples included HIV-1-positive samples from 200 people. In 1982 to 1984, HIV-1 subtypes A, C, and D were all present, though in small numbers. By 1987 to 1989, 152 (90%) of a total of 168 sequences were subtype C and AC, AD, and DC recombinants had emerged. Four of the subtype C sequences from 1983 to 1984 were closely related and were found at the base of a large cluster of low diversity that by the late 1980s accounted for 40% of C sequences. The other two early C sequences fell into a separate and more diverse cluster. Three other clusters containing sequences from the late 1980s were identified. Each cluster contained at least one sample from a person who had recently arrived in the district. From 18 HIV-1-positive spouse pairs, 12 very closely related pairs of sequences were identified. We conclude that there were multiple introductions of HIV-1 with limited spread, followed by explosive growth of a subtype C cluster, probably arising from a single introduction in or before 1983.     

Genetic diversity of NS3 protease from Brazilian HCV isolates and possible implications for therapy with direct-acting antiviral drugs

The hepatitis C virus (HCV) NS3 protease has been one of the molecular targets of new therapeutic approaches. Its genomic sequence variability in Brazilian HCV isolates is poorly documented. To obtain more information on the magnitude of its genetic diversity, 114 Brazilian HCV samples were sequenced and analysed together with global reference sequences. Genetic distance (d) analyses revealed that subtype 1b had a higher degree of heterogeneity (d = 0.098) than subtypes 1a (d = 0.060) and 3a (d = 0.062). Brazilian isolates of subtype 1b were distributed in the phylogenetic tree among sequences from other countries, whereas most subtype 1a and 3a sequences clustered into a single branch. Additional characterisation of subtype 1a in clades 1 and 2 revealed that all but two Brazilian subtype 1a sequences formed a distinct and strongly supported (approximate likelihood-ratio test = 93) group of sequences inside clade 1. Moreover, this subcluster inside clade 1 presented an unusual phenotypic characteristic in relation to the presence of resistance mutations for macrocyclic inhibitors. In particular, the mutation Q80K was found in the majority of clade 1 sequences, but not in the Brazilian isolates. These data demonstrate that Brazilian HCV subtypes display a distinct pattern of genetic diversity and reinforce the importance of sequence information in future therapeutic approaches.     

A search for base analogs to enhance third-strand binding to 'inverted' target base pairs of triplexes in the pyrimidine/parallel motif.

Eight base analogs were tested as third strand residues in otherwise homopyrimidine strands opposite each of the 'direct' (A.T and G.C) and 'inverted' (T.A and C.G) Watson-Crick base pairs, using UV melting profiles to assess triplex stability. The target duplexes contained 20 A.T base pairs and a central test base pair X.Y, while the third strand contained 20 T residues and a central Z test base. Z included 5-bromo-uracil, 5-propynyluracil, 5-propynylcytosine, 5-methyl-cytosine, 5-bromocytosine, hypoxanthine, 2-amino-purine and 2,6-diaminopurine. Some of the base analogs enhanced third strand binding to the target duplex with one or other 'inverted' central base pair relative to the binding afforded by any of the canonical bases. Other analogs did the same for the duplexes with the 'direct' target pairs. The increasing order of triplex stabilization by these base analogs is: opposite the 'inverted' base pairs, for T.A, A < C < 5-pC < 5-pU < T < 5-BrC < 5-meC < 5-BrU < 2-AP < 2,6-DAP < Hy < G, for C.G, 2-AP < A < Hy < G < 5-pC < 5-BrC < 5-meC < C < 2,6-DAP < T < 5-BrU < 5-pU; opposite the 'direct' base pairs, for A.T, 2-AP < A < 5-meC < C < G < Hy < 2,6-DAP < 5-pU < T = 5-BrU < 5-BrC < 5-pC, for G.C, G < 2,6-DAP < 2-AP < A < Hy < T < 5-BrU < 5-pU < 5-pC < 5-BrC < C < 5-meC.