Stability Of Base Pairs Research Articles

Design of the Crosslinking Reactions for Nucleic Acids-Binding Protein and Evaluation of the Reactivity

Selective chemical reactions of biomolecules are some of the important tools for investigations by biological studies. We have developed the selective crosslinking reactions to form covalent bonds to DNA or RNA using crosslinking oligonucleotides (CFO) bearing reactive bases. In this study, we designed the cross-linkable 4-amino-6-oxo-2-vinyltriazine derivative with an acyclic linker (acyAOVT) to react with the nucleic acids-binding protein based on our previous results. We hypothesized that the acyAOVT base would form a stable base pair with guanine by three hydrogen bonds at the positions of the vinyl group in the duplex DNA major groove, and the vinyl group can react with the nucleophilic species in the proximity, for example, the cysteine or lysine residue in the nucleic acids-binding protein. The synthesized oligonucleotides bearing the acyAOVT derivative showed a higher reactivity than that of the corresponding pyrimidine derivative without one nitrogen. The duplex containing acyAOVT-guanine (G) formed complexes with Hha1 DNMT even in the presence of 2-mercaptoethanol. We expect that our system will provide a useful tool for the molecular study of nucleic acids-binding proteins.

Complex RNA-DNA hybrid nanoshapes from iterative mix-and-match screening

Hybrid nucleic acid nanostructures partition architectural and functional roles between ribonucleic acid (RNA) joints and deoxyribonucleic acid (DNA) connectors. Nanoshapes self-assemble from nucleic acid modules through synergistic stabilization of marginally stable base pairing interactions within circularly closed polygons. Herein, we report the development of hybrid nanoshapes that include multiple different RNA modules such as internal loop and three-way junction (3WJ) motifs. An iterative mix-and-match screening approach was used to identify suitable DNA connectors that furnished stable nanoshapes for combinations of different RNA modules. The resulting complex multicomponent RNA-DNA hybrid nanoshapes were characterized by atomic force microscopy (AFM) imaging. Our research provides proof of concept for modular design, assembly and screening of RNA-DNA hybrid nanoshapes as building blocks for complex extended nucleic acid materials with features at the sub-10 nm scale.

Base pairing, structural and functional insights into N4-methylcytidine (m4C) and N4,N4-dimethylcytidine (m42C) modified RNA.

The N4-methylation of cytidine (m4C and m42C) in RNA plays important roles in both bacterial and eukaryotic cells. In this work, we synthesized a series of m4C and m42C modified RNA oligonucleotides, conducted their base pairing and bioactivity studies, and solved three new crystal structures of the RNA duplexes containing these two modifications. Our thermostability and X-ray crystallography studies, together with the molecular dynamic simulation studies, demonstrated that m4C retains a regular C:G base pairing pattern in RNA duplex and has a relatively small effect on its base pairing stability and specificity. By contrast, the m42C modification disrupts the C:G pair and significantly decreases the duplex stability through a conformational shift of native Watson-Crick pair to a wobble-like pattern with the formation of two hydrogen bonds. This double-methylated m42C also results in the loss of base pairing discrimination between C:G and other mismatched pairs like C:A, C:T and C:C. The biochemical investigation of these two modified residues in the reverse transcription model shows that both mono- or di-methylated cytosine bases could specify the C:T pair and induce the G to T mutation using HIV-1 RT. In the presence of other reverse transcriptases with higher fidelity like AMV-RT, the methylation could either retain the normal nucleotide incorporation or completely inhibit the DNA synthesis. These results indicate the methylation at N4-position of cytidine is a molecular mechanism to fine tune base pairing specificity and affect the coding efficiency and fidelity during gene replication.

Control of Transcription Initiation by Biased Thermal Fluctuations on Repetitive Genomic Sequences.

In the process of transcription initiation by RNA polymerase, promoter DNA sequences affect multiple reaction pathways determining the productivity of transcription. However, the question of how the molecular mechanism of transcription initiation depends on the sequence properties of promoter DNA remains poorly understood. Here, combining the statistical mechanical approach with high-throughput sequencing results, we characterize abortive transcription and pausing during transcription initiation by Escherichia coli RNA polymerase at a genome-wide level. Our results suggest that initially transcribed sequences, when enriched with thymine bases, contain the signal for inducing abortive transcription, whereas certain repetitive sequence elements embedded in promoter regions constitute the signal for inducing pausing. Both signals decrease the productivity of transcription initiation. Based on solution NMR and in vitro transcription measurements, we suggest that repetitive sequence elements within the promoter DNA modulate the nonlocal base pair stability of its double-stranded form. This stability profoundly influences the reaction coordinates of the productive initiation via pausing.

7-Iodo-5-aza-7-deazaguanine ribonucleoside: crystal structure, physical properties, base-pair stability and functionalization.

The positional change of nitrogen-7 of the RNA constituent guanosine to the bridgehead position-5 leads to the base-modified nucleoside 5-aza-7-deazaguanosine. Contrary to guanosine, this molecule cannot form Hoogsteen base pairs and the Watson-Crick proton donor site N3-H becomes a proton-acceptor site. This causes changes in nucleobase recognition in nucleic acids and has been used to construct stable `all-purine' DNA and DNA with silver-mediated base pairs. The present work reports the single-crystal X-ray structure of 7-iodo-5-aza-7-deazaguanosine, C10H12IN5O5 (1). The iodinated nucleoside shows an anti conformation at the glycosylic bond and an N conformation (O4'-endo) for the ribose moiety, with an antiperiplanar orientation of the 5'-hydroxy group. Crystal packing is controlled by interactions between nucleobase and sugar moieties. The 7-iodo substituent forms a contact to oxygen-2' of the ribose moiety. Self-pairing of the nucleobases does not take place. A Hirshfeld surface analysis of 1 highlights the contacts of the nucleobase and sugar moiety (O-H...O and N-H...O). The concept of pK-value differences to evaluate base-pair stability was applied to purine-purine base pairing and stable base pairs were predicted for the construction of `all-purine' RNA. Furthermore, the 7-iodo substituent of 1 was functionalized with benzofuran to detect motional constraints by fluorescence spectroscopy.

2-Ethynylnaphthalene-modified 8-aza-3,7-dideaza-2′-deoxyadenosine derivative discriminates thymine in target DNA via changes in the fluorescence wavelength

Novel environmentally sensitive fluorescent purine nucleosides, 2n37zA (1a), 2c37zA (1b), and 2a37zA (1c), were developed. Among them, 2c37zA (1b) and 2a37zA (1c), which contain electron-withdrawing substituents, exhibited highly solvatochromic properties and emitted at longer wavelengths when compared with the previously reported 1n37zA (2) (ca. 50–120 nm). The newly developed 2-ethynylnaphthalene derivatives, 2n37zA (1a), 2c37zA (1b), and 2a37zA (1c), switched emission modes between the ICT and the locally excited states by changing their molecular geometry. With their electron-withdrawing substituents, 2c37zA (1b) and 2a37zA (1c) had a wide range of fluorescence wavelength that was due to the combination of these two emission modes (Δλ = 142 and 96 nm, respectively). In the ODN duplexes, 2c37zA (1b) formed a stable base pair with thymine and was able to differentiate it through the respective fluorescence wavelength shift. Therefore, 2c37zA (1b) can be used for the detection of single-base alterations in the target nucleic acid sequences.

Adsorption of DNA/RNA nucleobases and base pairs on penta-graphene from first principles

The density functional theory with Grimme correction is performed to explore the interaction and adsorption mechanism of single nucleobases and hydrogen-bond base pairs in the deoxyribonucleic/ribonucleic acid (DNA/RNA) onto a penta-graphene (PG) unilaminar quasi two-dimensional (2D) sheet for the application of sensing of biomolecules. Adenine (A), guanine (G), cytosine (C), thymine (T), and uracil (U) are physisorbed onto PG as the adsorbates, reaching binding energy of −[0.50–1.10] eV and ranking the binding strength as G>A>C>T>U. Similarly, the hydrogen-bond base pairs are adsorbed onto PG with the binding energy up to −[1.40–1.70] eV and sorting the adsorption stability as GC>AT>AU. Above results confirm that the stability of nucleobases and base pairs on PG is improved reaching up to 8–44% compared with some traditional 2D materials, such as graphene, phosphorene and molybdenum disulfide. The remarkable electronic property changes of adsorbed PG make it be a potential candidate for rapid sensing, sequencing and identification of DNA/RNA elements. The new insights in this work reveal some physical and chemical mechanisms for the interaction between biomolecules and penta-graphene, suggesting that PG presents the capacity as a promising template for nanobiological devices and could be used for sensing of biomolecules.

Expanding the analytical applications of nucleic acid hybridization using junction probes.

Nucleic acid hybridization is crucial in target recognition with respect to in vitro and in vivo nucleic acid biosensing. Conventional linear probes and molecular beacons encounter challenges in multiplexing and specific recognition of intractable nucleic acids. Advances in nucleic acid nanotechnologies have resulted in a set of novel structural probes: junction probes (JPs), which make full use of the advantages of specificity, stability, programmability and predictability of Watson-Crick base pairing. In recent years, junction probes have been regularly implemented in constructing systems related to biosensing, synthetic biology and gene regulation. Herein, we summarize the latest advances in JP designs as potential nucleic acid biosensing systems and their expansive applications, and provide some general guidelines for developing JP based sensing strategies for implementation of such systems.

Melting temperature measurement and mesoscopic evaluation of single, double and triple DNA mismatches.

Unlike the canonical base pairs AT and GC, the molecular properties of mismatches such as hydrogen bonding and stacking interactions are strongly dependent on the identity of the neighbouring base pairs. As a result, due to the sheer number of possible combinations of mismatches and flanking base pairs, only a fraction of these have been studied in varying experiments or theoretical models. Here, we report on the melting temperature measurement and mesoscopic analysis of contiguous DNA mismatches in nearest-neighbours and next-nearest neighbour contexts. A total of 4032 different mismatch combinations, including single, double and triple mismatches were covered. These were compared with 64 sequences containing all combinations of canonical base pairs in the same location under the same conditions. For a substantial number of single mismatch configurations, 15%, the measured melting temperatures were higher than the least stable AT base pair. The mesoscopic calculation, using the Peyrard–Bishop model, was performed on the set of 4096 sequences, and resulted in estimates of on-site and nearest-neighbour interactions that can be correlated to hydrogen bonding and base stacking. Our results confirm many of the known properties of mismatches, including the peculiar sheared stacking of tandem GA mismatches. More intriguingly, it also reveals that a number of mismatches present strong hydrogen bonding when flanked on both sites by other mismatches. To highlight the applicability of our results, we discuss a number of practical situations such as enzyme binding affinities, thymine DNA glycosylase repair activity, and trinucleotide repeat expansions.

A hitherto unknown stability of DNA basepairs.

The resonance capture of ubiquitous very low energy electrons (vLEEs) into the π* orbitals of nucleobases is a potential doorway to DNA damage. Our ab initio quantum chemical calculations reveal a possible protection function, which is specific to the complementary basepairing, against such vLEE induced DNA damage.

Role of different tautomers in the base-pairing abilities of some of the vital antiviral drugs used against COVID-19.

Repurposed drugs are now considered as attractive therapeutics against COVID-19. It is shown that Remdesivir, a nucleoside drug that was originally invented for the Ebola virus, is effective in suppressing the replication of SARS-CoV-2 that causes COVID-19. Similarly, Galidesivir, Favipiravir, Ribavirin, N4-hydroxycytidine (EIDD-1931), and EIDD-2801 (a prodrug of EIDD-1931) were also found to be effective against COVID-19. However, the mechanisms of action of these drugs are not yet fully understood. For example, in some experimental studies, these drugs were proposed to act as a RNA-chain terminator, while in other studies, these were proposed to induce base-pair mutations above the error catastrophe limit to stall the replication of the viral RNA. To understand the mutagenic effects of these drugs, the role of different tautomers in their base-pairing abilities is studied here in detail by employing a reliable dispersion-corrected density functional theoretic method. It is found that Remdesivir and Galidesivir can adopt both amino and imino tautomeric conformations to base-pair with RNA bases. While the insertions of G and U are preferred against the amino tautomers of these drugs, the insertion of C is mainly possible against the imino tautomers. However, although Favipiravir and Ribavirin can make stable base pair interactions by using their keto and enol tautomers, the formation of the latter pairs would be less probable due to the endothermic nature of the products. Interestingly, the insertions of all of the RNA bases are found to be possible against the keto tautomer of Favipiravir, while the keto tautomer of Ribavirin has a clear preference for G. Remarkably, due to the negligible difference in the stability of EIDD-2801 and EIDD-1931, these tautomers would coexist in the biological environment. The insertion of G is found to be preferred against EIDD-1931 and the incorporations of U, A, and G are preferred opposite EIDD-2801. These findings suggest that base-pair mutations are the main causes of the antiviral properties of these drugs.

Accurate Base Pair Energies of Artificially Expanded Genetic Information Systems (AEGIS): Clues for Their Mutagenic Characteristics.

Recently, several artificial nucleobases, such as B, S, J, V, X, K, P, and Z, have been proposed to help in the expansion of the genetic information system and diagnosis of diseases. Among these bases, P and Z were identified to form stable DNA and to participate in the replication. However, the stabilities of P:Z and other artificial base pairs are not fully understood. The abilities of these unnatural nucleobases in mispairing with themselves and with natural bases are also not known. Here, the ωB97X-D dispersion-corrected density functional theoretical and complete basis set (CBS-QB3) methods are used to obtain accurate structural and energetic data related to base pair interactions involving these unnatural nucleobases. The roles of protonation and deprotonation of certain artificial bases in inducing mutations are also studied. It is found that each artificial purine has a complementary artificial pyrimidine, the base pair interactions between which are similar to those of the natural Watson-Crick base pairs. Hence, these base pairs will function naturally and would not impart mutagenicity. Among these base pairs, the J:V complex is found to be the most stable and promising artificial base pair. Remarkably, the noncomplementary artificial nucleobases are found to form stable mispairs, which may generate mutagenic products in DNA. Similarly, the misinsertions of natural bases opposite artificial bases are also found to be mutagenic. The mechanisms of these mutations are explained in detail. These results are in agreement with earlier biochemical studies. It is thus expected that this study would aid in the advancement of the synthetic biology to design more robust artificial nucleotides.

Anomeric 5-Aza-7-deaza-2'-deoxyguanosines in Silver-Ion-Mediated Homo and Hybrid Base Pairs: Impact of Mismatch Structure, Helical Environment, and Nucleobase Substituents on DNA Stability.

Nucleoside configuration (α-d vs. β-d), nucleobase substituents, and the helical DNA environment of silver-mediated 5-aza-7-deazaguanine-cytosine base pairs have a strong impact on DNA stability. This has been demonstrated by investigations on oligonucleotide duplexes with silver-mediated base pairs of α-d and β-d anomeric 5-aza-7-deaza-2'-deoxyguanosines and anomeric 2'-deoxycytidines incorporated in 12-mer duplexes. To this end, a new synthetic protocol has been developed to access the pure anomers of 5-aza-7-deaza-2'-deoxyguanosine by glycosylation of either the protected nucleobase or its salt followed by separation of the glycosylation products by crystallization and chromatography. Thermal stability measurements were performed on duplexes with α-d/α-d and β-d/β-d homo base pairs or α-d/β-d and β-d/α-d hybrid pairs within two sequence environments, positions 6 or 7, of oligonucleotide duplexes. The respective Tm stability increases observed after silver ion addition differ significantly. Homo base pairs with β-d/β-d or α-d/α-d nucleoside combinations are more stable than α-d/β-d hybrid base pairs. The positional switch of silver-ion-mediated base pairs has a significant impact on stability. Nucleobase substituents introduced at the 5-position of the dC site of silver-mediated base pairs affect base pair stability to a minor extent. Our investigation might lead to applications in the construction of bioinspired nanodevices, in DNA diagnostics, or metal-DNA hybrid materials.

Estimating Strengths of Individual Hydrogen Bonds in RNA Base Pairs: Toward a Consensus between Different Computational Approaches.

Noncoding RNA molecules are composed of a large variety of noncanonical base pairs that shape up their functionally competent folded structures. Each base pair is composed of at least two interbase hydrogen bonds (H-bonds). It is expected that the characteristic geometry and stability of different noncanonical base pairs are determined collectively by the properties of these interbase H-bonds. We have studied the ground-state electronic properties [using density functional theory (DFT) and DFT-D3-based methods] of all the 118 normal base pairs and 36 modified base pairs, belonging to 12 different geometric families (cis and trans of WW, WH, HH, WS, HS, and SS) that occur in a nonredundant set of high-resolution RNA crystal structures. Having addressed some of the limitations of the earlier approaches, we provide here a comprehensive compilation of the average energies of different types of interbase H-bonds (EHB). We have also characterized each interbase H-bond using 13 different parameters that describe its geometry, charge distribution at its bond critical point (BCP), and n → σ*-type charge transfer from filled π orbitals of the H-bond acceptor to the empty antibonding orbital of the H-bond donor. On the basis of the extent of their linear correlation with the H-bonding energy, we have shortlisted five parameters to model linear equations for predicting EHB values. They are (i) electron density at the BCP: ρ, (ii) its Laplacian: ∇2ρ, (iii) stabilization energy due to n → σ*-type charge transfer: E(2), (iv) donor–hydrogen distance, and (v) hydrogen–acceptor distance. We have performed single variable and multivariable linear regression analysis over the normal base pairs and have modeled sets of linear relationships between these five parameters and EHB. Performance testing of our model over the set of modified base pairs shows promising results, at least for the moderately strong H-bonds.

Functionalized Silver-Ion-Mediated α-dC/β-dC Hybrid Base Pairs with Exceptional Stability: α-d-5-Iodo-2'-Deoxycytidine and Its Octadiynyl Derivative in Metal DNA.

Silver-mediated α-dC-Ag+ -β-dC hybrid base pairs decorated with 5-iodo- or 5-octadiynyl residues are well accommodated in duplex DNA. A strong Tm increase and favorable thermodynamic data for duplex DNA were observed after addition of silver ions. The phenomenon is particularly obvious when both nucleobases of the base pairs are functionalized. Neither the position of the base pair, nor the type of 5-substituent had a negative influence. On the contrary, functionalization of conventional silver-mediated β-dC-Ag+ -β-dC homo base pairs showed a negative impact induced by the bulky substituents. To this end, cytosine modified 12-mer oligodeoxynucleotides were prepared by solid-phase synthesis employing new α-anomeric 2'-deoxycytidine phosphoramidites. A multigram scale synthesis was developed for 5-iodo-α-d-2'-deoxycytidine (1) employing the direct glycosylation of cytosine with Hoffer's α-d-halogenose followed by separation of anomeric DMT nucleosides. Regarding base-pair stability and functionalization silver-mediated α/β-dC hybrid base pairs were found to be superior to β/β-dC homo pairs. According to their extraordinary properties, they might find applications in DNA diagnostics, material science, or nanotechnology.

Evidence for Hidden Involvement of N3-Protonated Guanine in RNA Structure and Function.

Charged nucleobases have been found to occur in several known RNA molecules and are considered essential for their structure and function. The mechanism of their involvement is however not yet fully understood. Revelation of the role of N7-protonated guanine, in modulating the geometry and stability of noncanonical base pairs formed through its unprotonated edges [Watson–Crick (WC) and sugar], has triggered the need to evaluate the feasibility of similar roles of other protonated nucleobases [Halder et al., Phys Chem Chem Phys, 2015, 17, 26249]. In this context, N3 protonation of guanine makes an interesting case as its influence on the charge distribution of the WC edge is similar to that of N7 protonation, though its thermodynamic cost of protonation is significantly higher. In this work, we have carried out structural bioinformatics analyses and quantum mechanics-based calculations to show that N3 protonation of guanine may take place in a cellular environment, at least in the G:C W:W Trans and G:G W:H Cis base pairs. Our results provide a reasonable starting point for future investigations in order to address the larger mechanistic question.

Cyano Modification on Uridine Decreases Base-Pairing Stability and Specificity through Neighboring Disruption in RNA Duplex.

5-Cyanomethyluridine (cnm5 U) and 5-cyanouridine (cn5 U), the two uridine analogues, were synthesized and incorporated into RNA oligonucleotides. Base-pairing stability and specificity studies in RNA duplexes indicated that cnm5 U slightly decreased the stability of the duplex but retained the base-pairing preference. In contrast, cn5 U dramatically decreased both base-pairing stability and specificity between U:A and other noncanonical U:G, U:U, and U:C pairs. In addition, the cn5 U:G pair was found to be stronger than the cn5 U:A pair and the other mismatched pairs in the context of a RNA duplex; this implied that cn5 U might slightly prefer to recognize G over A. Our mechanistic studies by molecular simulations showed that the cn5 U modification did not directly affect the base pairing of the parent nucleotide; instead, it weakened the neighboring base pair in the 5' side of the modification in the RNA duplexes. Consistent with the simulation data, replacing the Watson-Crick A:U pair to a mismatched C:U pair in the 5'-neighboring site did not affect the overall stability of the duplex. Our work reveals the significance of the electron-withdrawing cyano group in natural tRNA systems and provides two novel building blocks for constructing RNA-based therapeutics.

Optimal Destabilization of DNA Double Strands by Single-Nucleobase Caging.

Photolabile protecting groups are widely used to trigger oligonucleotide activity. The ON/OFF-amplitude is a critical parameter. An experimental setup has been developed to identify protecting group derivatives with superior caging properties. Bulky rests are attached to the cage moiety via Cu-catalyzed azide-alkyne cycloaddition post-synthetically on DNA. Interestingly, the decrease in melting temperature upon introducing o-nitrobenzyl-caged (NPBY-) and diethylaminocoumarin-cages (DEACM-) in DNA duplexes reaches a limiting value. NMR spectroscopy was used to characterize individual base-pair stabilities and determine experimental structures of a selected number of photocaged DNA molecules. The experimental structures agree well with structures predicted by MD simulations. Combined, the structural data indicate that once a sterically demanding group is added to generate a tri-substituted carbon, the sterically less demanding cage moiety points towards the neighboring nucleoside and the bulkier substituents remain in the major groove.

A dinuclear Cu(i)-mediated complex: Theoretical studies of the G2Cu2 4+ cluster ion.

Recently, the T-Hg(ii)2-A base pair containing two equivalents of Hg(ii) has been prepared and characterized experimentally, which implies that there might exist considerable stable metal-mediated base pairs holding two neighbouring metal centers. Here we report a quantum chemical study on geometries, electronic structures, and bonding of various G2Cu2 4+ (G = guanine) isomers including one di-copper(i) unit. Different density functional methods [Becke 3-parameter-Lee-Yang-Parr, Perdew-Becke-Ernzerhof, Becke-Perdew, Density Functional Theory with Dispersion Corrections (DFT-D)] assign ambiguous relative energies to these isomers with the singlet and triplet states. High-level ab initio [domain-based local pair natural orbital (DLPNO) coupled-cluster with single and double excitations and DLPNO-coupled-cluster with single, double, and perturbative triple excitations] calculations confirm that the lowest-lying isomer of the G2Cu2 4+ ion has C 2h symmetry with the singlet state and is comparable to the singly and doubly charged homologues (G2Cu2 + and G2Cu2 2+). The extended transition state (ETS)-natural orbitals for the chemical valence (ETS-NOCV) calculations point out that it has larger instantaneous interaction energy and bond dissociation energy than the corresponding singly and doubly charged complexes due to its relatively stronger attractive energies and weaker Pauli repulsion. The orbital interactions in the quadruply charged cluster chiefly come from Cu2 4+ ← G⋯G π donations. The results may help the understanding of the bonding properties of other potential metal-base pair complexes with the electron transfer.

Structural Dynamics of Lateral and Diagonal Loops of Human Telomeric G-Quadruplexes in Extended MD Simulations.

The NMR solution structures of human telomeric (Htel) G-quadruplexes (GQs) are characterized by the presence of two lateral loops complemented by either diagonal or propeller loops. Bases of a given loop can establish interactions within the loop as well as with other loops and the flanking bases. This can lead to a formation of base alignments above and below the GQ stems. These base alignments are known to affect the loop structures and relative stabilities of different Htel GQ folds. We have carried out a total of 217 μs of classical (unbiased) molecular dynamics (MD) simulations starting from the available solution structures of Htel GQs to characterize structural dynamics of the lateral and diagonal loops, using several recent AMBER DNA force-field variants. As the loops are involved in diverse stacking and H-bonding interactions, their dynamics is slow, and extended sampling is required to capture different conformations. Nevertheless, although the simulations are far from being quantitatively converged, the data suggest that multiple 10 μs-scale simulations can provide a quite good assessment of the loop conformational space as described by the force field. The simulations indicate that the lateral loops may sample multiple coexisting conformations, which should be considered when comparing simulations with the NMR models as the latter include ensemble averaging. The adenine-thymine Watson-Crick arrangement was the most stable base pairing in the simulations. Adenine-adenine and thymine-thymine base pairs were also sampled but were less stable. The data suggest that the description of lateral and diagonal GQ loops in contemporary MD simulations is considerably more realistic than the description of propeller loops, though definitely not flawless.