Crick Pairs Research Articles

The stability and number of nucleating interactions determine DNA hybridization rates in the absence of secondary structure.

The kinetics of DNA hybridization are fundamental to biological processes and DNA-based technologies. However, the precise physical mechanisms that determine why different DNA sequences hybridize at different rates are not well understood. Secondary structure is one predictable factor that influences hybridization rates but is not sufficient on its own to fully explain the observed sequence-dependent variance. In this context, we measured hybridization rates of 43 different DNA sequences that are not predicted to form secondary structure and present a parsimonious physically justified model to quantify our observations. Accounting only for the combinatorics of complementary nucleating interactions and their sequence-dependent stability, the model achieves good correlation with experiment with only two free parameters. Our results indicate that greater repetition of Watson–Crick pairs increases the number of initial states able to proceed to full hybridization, with the stability of those pairings dictating the likelihood of such progression, thus providing new insight into the physical factors underpinning DNA hybridization rates.

A Novel AgI‐DNA Rod Comprising a One‐Dimensional Array of 11 Silver Ions within a Double Helical Structure

AbstractDNA/RNA duplexes containing metal‐ion‐mediated base pairs (metallo‐base pairs) have potential applications in developing nucleic acid‐based nanodevices and genetic code expansion. Many metallo‐base pairs are formed within duplexes stabilized by Watson–Crick base pairs. Recently, the crystal structure of an AgI‐DNA nanowire with an uninterrupted one‐dimensional silver array was determined. Here, we present a new DNA helical wire, the “AgI‐DNA rod”, containing an uninterrupted array of 11 AgI ions. The AgI‐DNA rod consisted of only C‐AgI‐C, G‐AgI‐G, G‐AgI‐5‐bromouracil (BrU), and BrU‐AgI‐BrU metallo base pairs, with no Watson–Crick pairs. The AgI‐DNA rods were connected by non‐canonical G–G pairs in crystals. Notably, data from our absorbance, circular dichroism, nuclear magnetic resonance, and mass spectrometry analyses suggested that the AgI‐DNA rods formed in solution, as well as within crystals.

How proton transfer affects the helical parameters in DNA:DNA microhelices

Proton transfer reactions are a widespread phenomenon in many areas of the life sciences and it is one of the origins of the spontaneous point mutations during DNA replication. Because of its importance, many studies have been reported on these reactions. However, the present work is the first one focused on the structural geometrical changes by double proton transfer (DPT). Thus, different Watson–Crick (WC) pairs were optimized first in a simple model with one nucleoside base pair, and in a microhelix form with three nucleoside base pairs. The canonical and few tautomeric forms were considered in DNA:DNA microhelices with A-type and B-type helical forms. The stability of these structures and how the DPT process affects the main geometrical parameters was analyzed, in particular the deformation of the helical parameters. The M06-2X DFT method was used for this purpose. The purine/pyrimidine ring in the keto form appears easier to be deformed than when it is in the enol form. The weaker WC base pair formed with mixed microhelices than with nucleobases alone and the significant deformation of the helical and backbone parameters with the DPT appears to complicate this process in microhelices. Communicated by Ramaswamy H. Sarma

Chemically Stabilized DNA Barcodes for DNA‐Encoded Chemistry

DNA‐encoded compound libraries are a widely used small molecule screening technology. One important aim in library design is the coverage of chemical space through structurally diverse molecules. Yet, the chemical reactivity of native DNA barcodes limits the toolbox of reactions for library design. Substituting the chemically vulnerable purines by 7‐deazaadenine, which exhibits tautomerization stability similar to natural adenine with respect to the formation of stable Watson–Crick pairs, yielded ligation‐competent, amplifiable, and readable DNA barcodes for encoded chemistry with enhanced stability against protic acid‐ and metal ion‐promoted depurination. The barcode stability allowed for straightforward translation of 16 exemplary reactions that included isocyanide multicomponent reactions, acid‐promoted Pictet–Spengler and Biginelli reactions, and metal‐promoted pyrazole syntheses on controlled pore glass‐coupled barcodes for diverse DEL design. The Boc protective group of reaction products offered a convenient handle for encoded compound purification.

U5 snRNA Interactions With Exons Ensure Splicing Precision.

Imperfect conservation of human pre-mRNA splice sites is necessary to produce alternative isoforms. This flexibility is combined with the precision of the message reading frame. Apart from intron-termini GU_AG and the branchpoint A, the most conserved are the exon-end guanine and +5G of the intron start. Association between these guanines cannot be explained solely by base-pairing with U1 snRNA in the early spliceosome complex. U6 succeeds U1 and pairs +5G in the pre-catalytic spliceosome, while U5 binds the exon end. Current U5 snRNA reconstructions by CryoEM cannot explain the conservation of the exon-end G. Conversely, human mutation analyses show that guanines of both exon termini can suppress splicing mutations. Our U5 hypothesis explains the mechanism of splicing precision and the role of these conserved guanines in the pre-catalytic spliceosome. We propose: (1) optimal binding register for human exons and U5—the exon junction positioned at U5Loop1 C39|C38; (2) common mechanism for base-pairing of human U5 snRNA with diverse exons and bacterial Ll.LtrB intron with new loci in retrotransposition—guided by base pair geometry; and (3) U5 plays a significant role in specific exon recognition in the pre-catalytic spliceosome. Statistical analyses showed increased U5 Watson–Crick pairs with the 5′exon in the absence of +5G at the intron start. In 5′exon positions −3 and −5, this effect is specific to U5 snRNA rather than U1 snRNA of the early spliceosome. Increased U5 Watson–Crick pairs with 3′exon position +1 coincide with substitutions of the conserved −3C at the intron 3′end. Based on mutation and X-ray evidence, we propose that −3C pairs with U2 G31 juxtaposing the branchpoint and the 3′intron end. The intron-termini pair, formed in the pre-catalytic spliceosome to be ready for transition after branching, and the early involvement of the 3′intron end ensure that the 3′exon contacts U5 in the pre-catalytic complex. We suggest that splicing precision is safeguarded cooperatively by U5, U6, and U2 snRNAs that stabilize the pre-catalytic complex by Watson–Crick base pairing. In addition, our new U5 model explains the splicing effect of exon-start +1G mutations: U5 Watson–Crick pairs with exon +2C/+3G strongly promote exon inclusion. We discuss potential applications for snRNA therapeutics and gene repair by reverse splicing.

Structure of an RNA helix with pyrimidine mismatches and cross-strand stacking.

The structure of a 22-base-pair RNA helix with mismatched pyrimidine base pairs is reported. The helix contains two symmetry-related CUG sequences: a triplet-repeat motif implicated in myotonic dystrophy type 1. The CUG repeat contains a U-U mismatch sandwiched between Watson-Crick pairs. Additionally, the center of the helix contains a dimerized UUCG motif with tandem pyrimidine (U-C/C-U) mismatches flanked by U-G wobble pairs. This region of the structure is significantly different from previously observed structures that share the same sequence and neighboring base pairs. The tandem pyrimidine mismatches are unusual and display sheared, cross-strand stacking geometries that locally constrict the helical width, a type of stacking previously associated with purines in internal loops. Thus, pyrimidine-rich regions of RNA have a high degree of structural diversity.

On the Cooperativity Effect in Watson and Crick and Wobble Pairs for a Halouracil Series and Its Potential Quantitative Application Studied through Surface-Enhanced Raman Spectroscopy.

The nature of the cooperativity effect of hydrogen bonds in Watson and Crick and wobble base pairs formed with thymine, uracil, and its 5-halogenated derivatives (5-fluoro, -chloro, and -bromouracil) has been studied through SERS and by using chemometric tools to process data and extract relevant information. Remarkable differences between the two kinds of pairs were clearly observed, and the behavior correlated to the withdrawing character of different substituents at the 5-position of uracil was verified. Multivariate analyses have also unveiled information about the pair's stability, and a stronger cooperativity effect seems to rule the Watson and Crick pairs when compared to wobble pairs. Defined patterns in the behavior of Watson and Crick pairs allowed the design of an indirect methodology for quantifying 5-bromouracil using a partial least squares (PLS) method with variable selection. Limit of detection (LOD) values of 0.037 and 0.112 mmol L-1 in the absence and presence of structurally similar interferences were reached, while its direct surface-enhanced Raman spectroscopy (SERS) quantification is only possible at ∼45 mmol L-1.

Improving RNA nearest neighbor parameters for helices by going beyond the two-state model

RNA folding free energy change nearest neighbor parameters are widely used to predict folding stabilities of secondary structures. They were determined by linear regression to datasets of optical melting experiments on small model systems. Traditionally, the optical melting experiments are analyzed assuming a two-state model, i.e. a structure is either complete or denatured. Experimental evidence, however, shows that structures exist in an ensemble of conformations. Partition functions calculated with existing nearest neighbor parameters predict that secondary structures can be partially denatured, which also directly conflicts with the two-state model. Here, a new approach for determining RNA nearest neighbor parameters is presented. Available optical melting data for 34 Watson–Crick helices were fit directly to a partition function model that allows an ensemble of conformations. Fitting parameters were the enthalpy and entropy changes for helix initiation, terminal AU pairs, stacks of Watson–Crick pairs and disordered internal loops. The resulting set of nearest neighbor parameters shows a 38.5% improvement in the sum of residuals in fitting the experimental melting curves compared to the current literature set.

Nuclear Magnetic Resonance Structure of an 8 × 8 Nucleotide RNA Internal Loop Flanked on Each Side by Three Watson-Crick Pairs and Comparison to Three-Dimensional Predictions.

The prediction of RNA three-dimensional structure from sequence alone has been a long-standing goal. High-resolution, experimentally determined structures of simple noncanonical pairings and motifs are critical to the development of prediction programs. Here, we present the nuclear magnetic resonance structure of the (5'CCAGAAACGGAUGGA)2 duplex, which contains an 8 × 8 nucleotide internal loop flanked by three Watson-Crick pairs on each side. The loop is comprised of a central 5'AC/3'CA nearest neighbor flanked by two 3RRs motifs, a known stable motif consisting of three consecutive sheared GA pairs. Hydrogen bonding patterns between base pairs in the loop, the all-atom root-mean-square deviation for the loop, and the deformation index were used to compare the structure to automated predictions by MC-sym, RNA FARFAR, and RNAComposer.

Structures of RNA repeats associated with neurological diseases.

All RNA molecules possess a 'propensity' to fold into complex secondary and tertiary structures. Although they are composed of only four types of nucleotides, they show an enormous structural richness which reflects their diverse functions in the cell. However, in some cases the folding of RNA can have deleterious consequences. Aberrantly expanded, repeated RNA sequences can exhibit gain-of-function abnormalities and become pathogenic, giving rise to many incurable neurological diseases. Most RNA repeats form long hairpin structures whose stem consists of noncanonical base pairs interspersed among Watson-Crick pairs. The expanded hairpins have an ability to sequester important proteins and form insoluble nuclear foci. The RNA pathology, common to many repeat disorders, has drawn attention to the structures of the RNA repeats. In this review, we summarize secondary structure probing and crystallographic studies of disease-related RNA repeat sequences. We discuss the unique structural features which can contribute to the pathogenic properties of the repeated runs. In addition, we present the newest reports concerning structural data linked to therapeutic approaches. WIREs RNA 2017, 8:e1412. doi: 10.1002/wrna.1412 For further resources related to this article, please visit the WIREs website.

Valid Plane Trees: Combinatorial Models for RNA Secondary Structures with Watson--Crick Base Pairs

The combinatorics of RNA plays a central role in biology. Mathematical biologists have several commonly used models for RNA: words in a fixed alphabet (representing the primary sequence of nucleotides) and plane trees (representing the secondary structure, or folding of the RNA sequence). This paper considers an augmented version of the standard model of plane trees, specifically one that incorporates some observed constraints on how the folding can occur. In particular we assume the alphabet consists of complementary pairs, for instance the Watson--Crick pairs A-U and C-G of RNA. Given a word in the alphabet, a valid plane tree is a tree for which, when the word is folded around the tree, each edge matches two complementary letters. Consider the graph whose vertices are valid plane trees for a fixed word and whose edges are given by Condon, Heitsch, and Hoos's local moves (see [C. E. Heitsch, Combinatorics on Plane Trees, Motivated by RNA Secondary Structure Configurations, preprint, Genome Center of Wisconsin and Department of Mathematics, University of Wisconsin--Madison]). We prove that this graph is connected. We give an explicit algorithm to construct a valid plane tree from a primary sequence, assuming that at least one valid plane tree exists. The tree produced by our algorithm has other useful characterizations, including a uniqueness condition defined by local moves. We also study enumerative properties of valid plane trees, analyzing how the number of valid plane trees depends on the choice of sequence length and alphabet size. Finally, we show that the proportion of words with at least one valid plane tree goes to zero as the word size increases. We also present some open questions.

The loss of a hydrogen bond: Thermodynamic contributions of a non-standard nucleotide.

Non-standard nucleotides are ubiquitous in RNA. Thermodynamic studies with RNA duplexes containing non-standard nucleotides, whether incorporated naturally or chemically, can provide insight into the stability of Watson–Crick pairs and the role of specific functional groups in stabilizing a Watson–Crick pair. For example, an A-U, inosine•U and pseudouridine•A pair each form two hydrogen bonds. However, an RNA duplex containing a central I•U pair or central Ψ•A pair is 2.4 kcal/mol less stable or 1.7 kcal/mol more stable, respectively, than the corresponding duplex containing an A-U pair. In the non-standard nucleotide purine, hydrogen replaces the exocyclic amino group of A. This replacement results in a P•U pair containing only one hydrogen bond. Optical melting studies were performed with RNA duplexes containing P•U pairs adjacent to different nearest neighbors. The resulting thermodynamic parameters were compared to RNA duplexes containing A-U pairs in order to determine the contribution of the hydrogen bond involving the exocyclic amino group. Results indicate a loss of 1.78 kcal/mol, on average, when an internal P•U replaces A-U in an RNA duplex. This value is compared to the thermodynamics of a hydrogen bond determined by similar methods. Nearest neighbor parameters were derived for use in free energy and secondary structure prediction software.

Tautomerism in some pyrimidine nucleoside analogues used in the treatment of cancer: an ab initio study

The tautomerism of pyrimidine \(2^\prime\)-deoxynucleoside, 2DN, has fundamental importance to understanding the mechanism of action in DNA replication and consequently to development of new drugs. In this work, the tautomerism of pyrimidine 2DN and some of its analogues with anticancer activity was analyzed along their potential energy surfaces by means of the MP2/cc-aug-pVDZ and QCISD/cc-aug-pVDZ theoretical methods of quantum chemistry. In gas phase and in solvent implicit (water), the energy barriers for hydrogen migration are above of 38.0 kcal/mol while tautomerization energies are between −4.9 and 10.2 kcal/mol. When hydrogen migration is catalyzed by one water molecule, the energy barriers for hydrogen migration are above of 14.2 kcal/mol, and tautomerization energies are found between 3.1 and 7.8 kcal/mol. This finding, supported by equilibrium constant and kinetic data, suggests that both tautomers, canonical and non-canonical, of each 2DN and analogues coexist in monohydrated medium forming mainly Watson–Crick pairs (WC) and some non-WC pairs. However, the calculated interaction energies of the pairs formed by guanine (G) and non-canonical pyrimidine (Pyr*) are higher than those WC pairs formed by G and Pyr, and therefore it is expected that the formation of these non-WC pairs plays an important role in the action of 2DN analogues in DNA replication.

Unravelling Protein-DNA Interactions at Molecular Level: A DFT and NCI Study.

Histone-DNA interactions were probed computationally at a molecular level, by characterizing the bimolecular clusters constituted by selected amino acid derivatives with polar (asparagine and glutamine), nonpolar (alanine, valine, and isoleucine), and charged (arginine) side chains and methylated pyrimidinic (1-methylcytosine and 1-methylthymine) and puric (9-methyladenine and 9-methylguanine) DNA bases. The computational approach combined different methodologies: a molecular mechanics (MMFFs forced field) conformational search and structural and vibrational density-functional calculations (M06-2X with double and triple-ζ Pople's basis sets). To dissect the interactions, intermolecular forces were analyzed with the Non-Covalent Interactions (NCI) analysis. The results for the 24 different clusters studied show a noticeable correlation between the calculated binding energies and the propensities for protein-DNA base interactions found in the literature. Such correlation holded even for the interaction of the selected amino acid derivatives with Watson and Crick pairs. Therefore, the balance between hydrogen bonds and van der Waals interactions (specially stacking) in the control of the final shape of the investigated amino acid-DNA base pairs seems to be well reproduced in dispersion-corrected DFT molecular models, reinforcing the idea that the specificity between the amino acids and the DNA bases play an important role in the regulation of DNA.

Ultraviolet Absorption Induces Hydrogen‐Atom Transfer in G⋅C Watson–Crick DNA Base Pairs in Solution

AbstractUltrafast deactivation pathways bestow photostability on nucleobases and hence preserve the structural integrity of DNA following absorption of ultraviolet (UV) radiation. One controversial recovery mechanism proposed to account for this photostability involves electron‐driven proton transfer (EDPT) in Watson–Crick base pairs. The first direct observation is reported of the EDPT process after UV excitation of individual guanine–cytosine (G⋅C) Watson–Crick base pairs by ultrafast time‐resolved UV/visible and mid‐infrared spectroscopy. The formation of an intermediate biradical species (G[−H]⋅C[+H]) with a lifetime of 2.9 ps was tracked. The majority of these biradicals return to the original G⋅C Watson–Crick pairs, but up to 10 % of the initially excited molecules instead form a stable photoproduct G*⋅C* that has undergone double hydrogen‐atom transfer. The observation of these sequential EDPT mechanisms across intermolecular hydrogen bonds confirms an important and long debated pathway for the deactivation of photoexcited base pairs, with possible implications for the UV photochemistry of DNA.

Ultraviolet Absorption Induces Hydrogen-Atom Transfer in G⋅C Watson-Crick DNA Base Pairs in Solution.

Ultrafast deactivation pathways bestow photostability on nucleobases and hence preserve the structural integrity of DNA following absorption of ultraviolet (UV) radiation. One controversial recovery mechanism proposed to account for this photostability involves electron-driven proton transfer (EDPT) in Watson-Crick base pairs. The first direct observation is reported of the EDPT process after UV excitation of individual guanine-cytosine (G⋅C) Watson-Crick base pairs by ultrafast time-resolved UV/visible and mid-infrared spectroscopy. The formation of an intermediate biradical species (G[-H]⋅C[+H]) with a lifetime of 2.9 ps was tracked. The majority of these biradicals return to the original G⋅C Watson-Crick pairs, but up to 10% of the initially excited molecules instead form a stable photoproduct G*⋅C* that has undergone double hydrogen-atom transfer. The observation of these sequential EDPT mechanisms across intermolecular hydrogen bonds confirms an important and long debated pathway for the deactivation of photoexcited base pairs, with possible implications for the UV photochemistry of DNA.

Polymerase Bypass of N(6)-Deoxyadenosine Adducts Derived from Epoxide Metabolites of 1,3-Butadiene.

N(6)-(2-Hydroxy-3-buten-1-yl)-2'-deoxyadenosine (N(6)-HB-dA I) and N(6),N(6)-(2,3-dihydroxybutan-1,4-diyl)-2'-deoxyadenosine (N(6),N(6)-DHB-dA) are exocyclic DNA adducts formed upon alkylation of the N(6) position of adenine in DNA by epoxide metabolites of 1,3-butadiene (BD), a common industrial and environmental chemical classified as a human and animal carcinogen. Since the N(6)-H atom of adenine is required for Watson-Crick hydrogen bonding with thymine, N(6)-alkylation can prevent adenine from normal pairing with thymine, potentially compromising the accuracy of DNA replication. To evaluate the ability of BD-derived N(6)-alkyladenine lesions to induce mutations, synthetic oligodeoxynucleotides containing site-specific (S)-N(6)-HB-dA I and (R,R)-N(6),N(6)-DHB-dA adducts were subjected to in vitro translesion synthesis in the presence of human DNA polymerases β, η, ι, and κ. While (S)-N(6)-HB-dA I was readily bypassed by all four enzymes, only polymerases η and κ were able to carry out DNA synthesis past (R,R)-N(6),N(6)-DHB-dA. Steady-state kinetic analyses indicated that all four DNA polymerases preferentially incorporated the correct base (T) opposite (S)-N(6)-HB-dA I. In contrast, hPol β was completely blocked by (R,R)-N(6),N(6)-DHB-dA, while hPol η and κ inserted A, G, C, or T opposite the adduct with similar frequency. HPLC-ESI-MS/MS analysis of primer extension products confirmed that while translesion synthesis past (S)-N(6)-HB-dA I was mostly error-free, replication of DNA containing (R,R)-N(6),N(6)-DHB-dA induced significant numbers of A, C, and G insertions and small deletions. These results indicate that singly substituted (S)-N(6)-HB-dA I lesions are not miscoding, but that exocyclic (R,R)-N(6),N(6)-DHB-dA adducts are strongly mispairing, probably due to their inability to form stable Watson-Crick pairs with dT.

The Importance of Short- and Long-Range Exchange on Various Excited State Properties of DNA Monomers, Stacked Complexes, and Watson-Crick Pairs.

We present a detailed analysis of several time-dependent DFT (TD-DFT) methods, including conventional hybrid functionals and two types of nonempirically tuned range-separated functionals, for predicting a diverse set of electronic excitations in DNA nucleobase monomers and dimers. This large and extensive set of excitations comprises a total of 50 different transitions (for each tested DFT functional) that includes several n → π and π → π* valence excitations, long-range charge-transfer excitations, and extended Rydberg transitions (complete with benchmark calculations from high-level EOM-CCSD(T) methods). The presence of localized valence excitations as well as extreme long-range charge-transfer excitations in these systems poses a serious challenge for TD-DFT methods that allows us to assess the importance of both short- and long-range exchange contributions for simultaneously predicting all of these various transitions. In particular, we find that functionals that do not have both short- and full long-range exchange components are unable to predict the different types of nucleobase excitations with the same accuracy. Most importantly, the current study highlights the importance of both short-range exchange and a nonempirically tuned contribution of long-range exchange for accurately predicting the diverse excitations in these challenging nucleobase systems.

Carcinogenesis of urethane: simulation versus experiment.

The carcinogenesis of urethane (ethyl carbamate), a byproduct of fermentation that is consistently found in various food products, was investigated with a combination of kinetic experiments and quantum chemical calculations. The main objective of the study was to find ΔG(⧧), the activation free energy for the rate-limiting step of the SN2 reaction among the ultimate carcinogen of urethane, vinyl carbamate epoxide (VCE), and different nucleobases of the DNA. In the experimental part, the second-order reaction rate constants for the formation of the main 7-(2-oxoethyl)guanine adduct in aqueous solutions of deoxyguanosine and in DNA were determined. A series of ab initio, density functional theory (DFT), and semiempirical molecular orbital (MO) calculations was then performed to determine the activation barriers for the reaction between VCE and nucleobases methylguanine, methyladenine, and methylcytosine. Effects of hydration were incorporated with the use of the solvent reaction field method of Tomasi and co-workers and the Langevine dipoles model of Florian and Warshel. The computational results for the main adduct were found to be in good agreement with the experiment, thus presenting strong evidence for the validity of the proposed SN2 mechanism. This allowed us to predict the activation barriers of reactions leading to side products for which kinetic experiments have not yet been performed. Our calculations have shown that the main 7-(2-oxoethyl)deoxyguanosine adduct indeed forms preferentially because the emergence of other adducts either proceeds across a significantly higher activation barrier or the geometry of the reaction requires the Watson-Crick pairs of the DNA to be broken. The computational study also considered the questions of stereoselectivity, the ease of the elimination of the leaving group, and the relative contributions of the two possible reaction paths for the formation of the 1,N(2)-ethenodeoxyguanosine adduct.

The contribution of pseudouridine to stabilities and structure of RNAs

Thermodynamic data are reported revealing that pseudouridine (Ψ) can stabilize RNA duplexes when replacing U and forming Ψ-A, Ψ-G, Ψ-U and Ψ-C pairs. Stabilization is dependent on type of base pair, position of Ψ within the RNA duplex, and type and orientation of adjacent Watson–Crick pairs. NMR spectra demonstrate that for internal Ψ-A, Ψ-G and Ψ-U pairs, the N3 imino proton is hydrogen bonded to the opposite strand nucleotide and the N1 imino proton may also be hydrogen bonded. CD spectra show that general A-helix structure is preserved, but there is some shifting of peaks and changing of intensities. Ψ has two hydrogen donors (N1 and N3 imino protons) and two hydrogen bond acceptors because the glycosidic bond is C-C rather than C-N as in uridine. This greater structural potential may allow Ψ to behave as a kind of structurally driven universal base because it can enhance stability relative to U when paired with A, G, U or C inside a double helix. These structural and thermodynamic properties may contribute to the biological functions of Ψ.