Properties of the Nucleic-acid Bases in Free and Watson-Crick Hydrogen-bonded States: Computational Insights into the Sequence-dependent Features of Double-helical DNA

Base pair dynamics, electrostatics, and thermodynamics at the LTR-III quadruplex:duplex junction

G-quadruplexes (GQs) play key regulatory roles within the human genome and have also been identified to play similar roles in other eukaryotes, bacteria, archaea, and viruses. Human immunodeficiency virus 1, the etiological agent of acquired immunodeficiency syndrome, can form two GQs in its long terminal repeat (LTR) promoter region, each of which act to regulate viral gene expression in opposing manners. The major LTR GQ, called LTR-III, is a distinct hybrid GQ containing a 12-nucleotide duplex loop attached to the quadruplex motif. The resulting quadruplex:duplex junction (QDJ) has been hypothesized to serve as a selective drug targeting site. To better understand the dynamics of this QDJ, we performed conventional and enhanced-sampling molecular dynamics simulations using the Drude-2017 force field. We observed unbiased and reversible formation of additional base pairs in the QDJ, between Ade4:Thy14 and Gua3:Thy14. Both base pairs were electrostatically favored, but geometric constraints within the junction may drive the formation of, and preference for, the Ade4:Thy14 base pair. Finally, we demonstrated that the base pairs are separated only by small energy barriers that may enable transitions between both base-paired states. Together, these simulations provide new insights into the dynamics, electrostatics, and thermodynamics of the LTR-III QDJ.

Insights into DNA solvation found in protein-DNA structures

The proteins that bind double-helical DNA present various microenvironments that sense and/or induce signals in the genetic material. The high-resolution structures of protein-DNA complexes reveal the nature of both the microenvironments and the conformational responses in DNA and protein. Complex networks of interactions within the structures somehow tie the protein and DNA together and induce the observed spatial forms. Here we show how the cumulative buildup of amino acid atoms around the sugars, phosphates, and bases in different protein-DNA complexes produces a binding cloud around the double helix and how different types of atoms fill that cloud. Rather than focusing on the principles of molecular binding and recognition suggested by the arrangements of amino acids and nucleotides in the macromolecular complexes, we consider the proteins in contact with DNA as organized solvents. We describe differences in the mix of atoms that come in closest contact with DNA, subtle sequence-dependent features in the microenvironment of the sugar-phosphate backbone, a direct link between the localized buildup of ionic species and the electrostatic potential surfaces of the DNA bases, and sites of atomic buildup above and below the basepair planes that transmit the unique features of the base environments along the chain backbone. The inferences about solvation that can be drawn from the survey provide new stimuli for improvement of nucleic acid force fields and fresh ideas for exploration of the properties of DNA in solution.

DNA-Templated Modulation in the Photophysical Properties of a Fluorescent Molecular Rotor Auramine O by Varying the DNA Composition

This work delineates an integrative approach combining spectroscopic and computational studies to decipher the association-induced fluorescence properties of a fluorescent molecular rotor, viz., auramine O (AuO), after interacting with 20-mer duplex DNA having diverse well-matched base pairs. While exploring the scarcely explored sequence-dependent interaction mechanism of AuO and DNA, we observed that DNA could act as a conducive scaffold to the formation of AuO dimer through noncovalent interactions at lower molecular density. The photophysical properties of AuO depend on the nucleotide compositions as described from sequence-dependent shifting in the emission and absorption maxima. Furthermore, we explored such DNA base pair-dependent fluorescence spectral characteristics of AuO toward discriminating the thermodynamically most stable single nucleotide mismatch in a 20-mer sequence. Our results are interesting and could be useful in developing analogues with further enhanced emission properties toward mismatched DNA sequences.

Electron ratcheting in self-assembled soft matter

Ratcheted multi-step hopping electron transfer systems can plausibly produce directional charge transport over very large distances without requiring a source-drain voltage bias. We examine molecular strategies to realize ratcheted charge transport based on multi-step charge hopping, and we illustrate two ratcheting mechanisms with examples based on DNA structures. The charge transport times and currents that may be generated in these assemblies are also estimated using kinetic simulations. The first ratcheting mechanism described for nanoscale systems requires local electric fields on the 10⁹ V/m scale to realize nearly 100% population transport. The second ratcheting mechanism for even larger systems, based on electrochemical gating, is estimated to generate currents as large as 0.1 pA for DNA structures that are a few μm in length with a gate voltage of about 5 V, a magnitude comparable to currents measured in DNA wires at the nanoscale when a source-drain voltage bias of similar magnitude is applied, suggesting an approach to considerably extend the distance range over which DNA charge transport devices may operate.

Quantification of the mesh structure of bundled actin filaments

Biopolymer networks are essential for a wide variety of cellular functions. The biopolymer actin is known to self-assemble into a variety of spatial structures in response to physiological and physical mechanisms. So far, the mechanics of networks of single actin filaments and bundles has previously been described. However, the spatial structure of actin bundles remains poorly understood. Here, we investigate this question by bundling actin filaments with systematically varied concentrations of known physical bundling agents (MgCl2 and PEG) and physiological bundling agents (α-actinin and fascin). We image bundled actin networks with confocal microscopy and perform analysis to describe their mesh size and the nearest-distance distribution, which we call "mesh structure". We find that the mesh size ξ scales universally with actin concentration as ξ ∼ [actin]-1/2. However, the dependence of ξ on the concentration of the bundling agent depends on the agent used. Finally, we find that nearest-distance distributions are best fit by Weibull and Gamma distributions. A complete understanding of the mesh structure of biopolymer networks leads to a more mechanistic understanding of the structure of the cytoskeleton, and can be exploited to design filters with variable porosity for microfluidic devices.

Preferential Binding of Thioflavin T to AT-Rich DNA: White Light Emission through Intramolecular Förster Resonance Energy Transfer

Herein we report the effect of different nucleobase pair compositions on the association-induced fluorescence enhancement property of Thioflavin T (ThT), upon binding with 20 base pair long double-stranded DNA (dsDNA). Analysis of binding and decay constants along with the association (K_ass) and dissociation (K_diss) rate constants obtained from the fluctuation in the fluorescence intensity of ThT after binding with different DNA revealed selective affinity of ThT toward AT-rich dsDNA. Molecular docking also substantiates the experimental results. We also observed that addition of orange-emitting ethidium bromide (EtBr) to cyan-emitting ThT-DNA complexes leads to bright white light emission (WLE) through Förster resonance energy transfer. Additionally, the emission of white light is far greater in the case of intra-DNA strands. Besides endorsing the binding insights of ThT to AT-rich dsDNA, the present investigations open a new perspective for realizing promising WLE from two biomarkers without labeling the DNA.

Theoretical analyses and experimental validation of the effects caused by the fluorinated substituent modification of DNA

Halogen-modified nucleic acid molecules, such as trifluorothymidine (FTD) and 5-fluorouracil, are widely used in medical science and clinical site. These compounds have a very similar nucleobase structure. It is reported that both of these compounds could be incorporated into DNA. The incorporation of FTD produces highly anti-tumor effect. However, it is not known whether to occur a significant effect by the incorporation of 5-fluorouracil. Nobody knows why such a difference will occur. To understand the reason why there is large differences between trifluorothymidine and 5-fluorouracil, we have performed the molecular dynamics simulations and molecular orbital calculations. Although the active interaction energy between Halogen-modified nucleic acids or and complementary adenine was increased, in only FTD incorporated DNA, more strongly dispersion force interactions with an adjacent base were detected in many thermodynamic DNA conformations. As the results, the conformational changes occur even if it is in internal body temperature. Then the break of hydrogen bonding between FTD and complementary adenine base occur more frequently. The double helix structural destabilization of DNA with FTD is resulted from autoagglutination caused by the bonding via halogen orbitals such as halogen bonding and the general van der Waals interactions such as CH–\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\rm{\pi }}$$\end{document}π, lone pair (LP)–\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\rm{\pi }}$$\end{document}π, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\rm{\pi }}$$\end{document}π–\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\rm{\pi }}$$\end{document}π interactions. Therefore, it is strongly speculated that such structural changes caused by trifluoromethyl group is important for the anti-tumor effect of FTD alone.

SERS assessment of the cancer-specific methylation pattern of genomic DNA: towards the detection of acute myeloid leukemia in patients undergoing hematopoietic stem cell transplantation

In this label-free surface-enhanced Raman scattering (SERS) study of genomic DNA, we demonstrate that the cancer-specific DNA methylation pattern translates into specific spectral differences. Thus, DNA extracted from an acute myeloid leukemia (AML) cell line presented a decreased intensity of the 1005 cm^-1 band of 5-methylcytosine compared to normal DNA, in line with the well-described hypomethylation of cancer DNA. The unique methylation pattern of cancer DNA also influences the DNA adsorption geometry, resulting in higher adenine SERS intensities for cancer DNA. The possibility of detecting cancer DNA based on its SERS spectrum was validated on peripheral blood genomic DNA samples from n = 17 AML patients and n = 17 control samples, yielding an overall classification of 82% based on the 1005 cm^-1 band of 5-methylcytosine. By demonstrating the potential of SERS in assessing the methylation status in the case of real-life DNA samples, the study paves the way for novel methods of diagnosing cancer. Graphical abstract.

2'-O,4'-C-Methylene-Bridged Nucleic Acids Stabilize Metal-Mediated Base Pairing in a DNA Duplex

The 2'-O,4'-C-methylene-bridged or locked nucleic acid (2',4'-BNA/LNA), with an N-type sugar conformation, effectively improves duplex-forming ability. 2',4'-BNA/LNA is widely used to improve gene knockdown in nucleic acid based therapies and is used in gene diagnosis. Metal-mediated base pairs (MMBPs), such as thymine (T)-Hg^II -T and cytosine (C)-Ag^I -C have been developed and used as attractive tools in DNA nanotechnology studies. This study aimed to investigate the application of 2',4'-BNA/LNA in the field of MMBPs. 2',4'-BNA/LNA with 5-methylcytosine stabilized the MMBP of C with Ag^I ions. Moreover, the 2',4'-BNA/LNA sugar significantly improved the duplex-forming ability of the DNA/DNA complex, relative to that by the unmodified sugar. These results suggest that the sugar conformation is important for improving the stability of duplex-containing MMBPs. The results indicate that 2',4'-BNA/LNA can be applied not only to nucleic acid based therapies, but also to MMBP technologies.

Contributions of Sequence to the Higher-Order Structures of DNA

One of the critical unanswered questions in genome biophysics is how the primary sequence of DNA bases influences the global properties of very-long-chain molecules. The local sequence-dependent features of DNA found in high-resolution structures introduce irregularities in the disposition of adjacent residues that facilitate the specific binding of proteins and modulate the global folding and interactions of double helices with hundreds of basepairs. These features also determine the positions of nucleosomes on DNA and the lengths of the interspersed DNA linkers. Like the patterns of basepair association within DNA, the arrangements of nucleosomes in chromatin modulate the properties of longer polymers. The intrachromosomal loops detected in genomic studies contain hundreds of nucleosomes, and given that the simulated configurations of chromatin depend on the lengths of linker DNA, the formation of these loops may reflect sequence-dependent information encoded within the positioning of the nucleosomes. With knowledge of the positions of nucleosomes on a given genome, methods are now at hand to estimate the looping propensities of chromatin in terms of the spacing of nucleosomes and to make a direct connection between the DNA base sequence and larger-scale chromatin folding.

On the interpretation of electron microscopic maps of biological macromolecules

The images of flash-frozen biological macromolecules produced by cryo-electron microscopy (EM) can be used to generate accurate, three-dimensional, electric potential maps for these molecules that resemble X-ray-derived electron density maps. However, unlike electron density maps, electric potential maps can include negative features that might for example represent the negatively charged, backbone phosphate groups of nucleic acids or protein carboxylate side chains, which can complicate their interpretation. This study examines the images of groups that include charged atoms that appear in recently-published, high-resolution EM potential maps of the ribosome and β-galactosidase. Comparisons of simulated maps of these same groups with their experimental counterparts highlight the impact that charge has on the appearance of electric potential maps.

Coarse-grained modeling of RNA 3D structure

Functional RNA molecules depend on three-dimensional (3D) structures to carry out their tasks within the cell. Understanding how these molecules interact to carry out their biological roles requires a detailed knowledge of RNA 3D structure and dynamics as well as thermodynamics, which strongly governs the folding of RNA and RNA-RNA interactions as well as a host of other interactions within the cellular environment. Experimental determination of these properties is difficult, and various computational methods have been developed to model the folding of RNA 3D structures and their interactions with other molecules. However, computational methods also have their limitations, especially when the biological effects demand computation of the dynamics beyond a few hundred nanoseconds. For the researcher confronted with such challenges, a more amenable approach is to resort to coarse-grained modeling to reduce the number of data points and computational demand to a more tractable size, while sacrificing as little critical information as possible. This review presents an introduction to the topic of coarse-grained modeling of RNA 3D structures and dynamics, covering both high- and low-resolution strategies. We discuss how physics-based approaches compare with knowledge based methods that rely on databases of information. In the course of this review, we discuss important aspects in the reasoning process behind building different models and the goals and pitfalls that can result.

Homologous Pairing between Long DNA Double Helices

Molecular recognition between two double stranded (ds) DNA with homologous sequences may not seem compatible with the B-DNA structure because the sequence information is hidden when it is used for joining the two strands. Nevertheless, it has to be invoked to account for various biological data. Using quantum chemistry, molecular mechanics, and hints from recent genetics experiments, I show here that direct recognition between homologous dsDNA is possible through the formation of short quadruplexes due to direct complementary hydrogen bonding of major-groove surfaces in parallel alignment. The constraints imposed by the predicted structures of the recognition units determine the mechanism of complexation between long dsDNA. This mechanism and concomitant predictions agree with the available experimental data and shed light upon the sequence effects and the possible involvement of topoisomerase II in the recognition.

Triple helical DNA in a duplex context and base pair opening

It is fundamental to explore in atomic detail the behavior of DNA triple helices as a means to understand the role they might play in vivo and to better engineer their use in genetic technologies, such as antigene therapy. To this aim we have performed atomistic simulations of a purine-rich antiparallel triple helix stretch of 10 base triplets flanked by canonical Watson–Crick double helices. At the same time we have explored the thermodynamic behavior of a flipping Watson–Crick base pair in the context of the triple and double helix. The third strand can be accommodated in a B-like duplex conformation. Upon binding, the double helix changes shape, and becomes more rigid. The triple-helical region increases its major groove width mainly by oversliding in the negative direction. The resulting conformations are somewhere between the A and B conformations with base pairs remaining almost perpendicular to the helical axis. The neighboring duplex regions maintain a B DNA conformation. Base pair opening in the duplex regions is more probable than in the triplex and binding of the Hoogsteen strand does not influence base pair breathing in the neighboring duplex region.

B-DNA characteristics are preserved in double stranded d(A)3·d(T)3 and d(G)3·d(C)3 mini-helixes: conclusions from DFT/M06-2X study

We report the results of the first comprehensive DFT study on the d(A)3·d(T)3 and d(G)3·d(C)3 nucleic acid duplexes. The ability of mini-helixes to preserve the conformation of B-DNA in the gas phase and under the influence of such factors as: solvent, uncompensated charge, and counter-ions was evaluated using M06-2X functional with 6-31G(d,p) basis set. The accuracy of the models was ascertained based on their ability to reproduce key structural features of natural B-DNA. Analysis of the helicity suggests that the helical conformations adopt geometrical parameters which are close to those of the B-DNA form. The torsion angles fall somewhere between the values observed for BI/BII conformational classes. The comparative analysis of parameters of isolated Watson-Crick base pairs versus B-DNA-like conformations indicates the same tendency of base-pair polarization and hydration. Specifically, effects of polarization of nucleobases in continuum type dielectric medium mimicking water are stronger than those caused by the presence of backbone. Polar environment as well as the presence of counterions stabilizes duplexes, facilitating helix formation. Substantial conformational changes of nucleotides upon duplex formation decrease the binding energy. In spite of structural and energetic changes, the placement of a mini-helix into the gas phase does not lead to significant disruption of the structure. On the contrary, the duplex preserves its helicity and the strands remain bound.

In Darwinian evolution, feedback from natural selection leads to biased mutations

Natural selection provides feedback through which information about the environment and its recurring challenges is captured, inherited, and accumulated within genomes in the form of variations that contribute to survival. The variation upon which natural selection acts is generally described as "random." Yet evidence has been mounting for decades, from such phenomena as mutation hotspots, horizontal gene transfer, and highly mutable repetitive sequences, that variation is far from the simplifying idealization of random processes as white (uniform in space and time and independent of the environment or context).  This paper focuses on what is known about the generation and control of mutational variation, emphasizing that it is not uniform across the genome or in time, not unstructured with respect to survival, and is neither memoryless nor independent of the (also far from white) environment. We suggest that, as opposed to frequentist methods, Bayesian analysis could capture the evolution of nonuniform probabilities of distinct classes of mutation, and argue not only that the locations, styles, and timing of real mutations are not correctly modeled as generated by a white noise random process, but that such a process would be inconsistent with evolutionary theory.

DNA self-assembly: from chirality to evolution

Transient or long-term DNA self-assembly participates in essential genetic functions. The present review focuses on tight DNA-DNA interactions that have recently been found to play important roles in both controlling DNA higher-order structures and their topology. Due to their chirality, double helices are tightly packed into stable right-handed crossovers. Simple packing rules that are imposed by DNA geometry and sequence dictate the overall architecture of higher order DNA structures. Close DNA-DNA interactions also provide the missing link between local interactions and DNA topology, thus explaining how type II DNA topoisomerases may sense locally the global topology. Finally this paper proposes that through its influence on DNA self-assembled structures, DNA chirality played a critical role during the early steps of evolution.

Mg2+ in the major groove modulates B-DNA structure and dynamics

This study investigates the effect of Mg(2+) bound to the DNA major groove on DNA structure and dynamics. The analysis of a comprehensive dataset of B-DNA crystallographic structures shows that divalent cations are preferentially located in the DNA major groove where they interact with successive bases of (A/G)pG and the phosphate group of 5'-CpA or TpG. Based on this knowledge, molecular dynamics simulations were carried out on a DNA oligomer without or with Mg(2+) close to an ApG step. These simulations showed that the hydrated Mg(2+) forms a stable intra-strand cross-link between the two purines in solution. ApG generates an electrostatic potential in the major groove that is particularly attractive for cations; its intrinsic conformation is well-adapted to the formation of water-mediated hydrogen bonds with Mg(2+). The binding of Mg(2+) modulates the behavior of the 5'-neighboring step by increasing the BII (ε-ζ>0°) population of its phosphate group. Additional electrostatic interactions between the 5'-phosphate group and Mg(2+) strengthen both the DNA-cation binding and the BII character of the 5'-step. Cation binding in the major groove may therefore locally influence the DNA conformational landscape, suggesting a possible avenue for better understanding how strong DNA distortions can be stabilized in protein-DNA complexes.

Theoretical studies of the base pair fidelity of selenium-modified DNA

The introduction of selenium into DNA in the place of oxygen provides a unique opportunity for studying the fidelity of DNA replication, as well as providing advantages in the growth of DNA crystals and the greater resolution of their structures. However, the atomic mechanisms of the relative stability and base pair recognition of the selenium-modified DNA are poorly understood. In the present study, quantum mechanics calculations were performed on base pairings, base stacking, and base-water interactions for both unmodified thymine and thymine with the 2-exo-oxygen replaced with selenium, and the results were used to develop and validate CHARMM force field parameters for the 2-Se-thymine. Subsequently, molecular dynamics simulations and free-energy perturbation calculations were performed on 11-base DNA sequences containing native thymine and the 2-Se-thymine. The calculated relative free-energy values are in good agreement with experimentally determined relative stability, where the 2-Se-thymine offers similar stability to T-A in cognate DNA, while it dramatically destabilizes the DNA containing the T-G mismatch base pair when 2-Se-thymine is incorporated. Thus, 2-Se-thymine largely increases the native T-A base pair fidelity by discouraging the T-G wobble pair. Insights into the high pairing specificity and the relative stability of selenium-modified DNA were obtained based on detailed structural and energetic analysis of molecular dynamics trajectories. Our studies move one step further toward an understanding of the high base pair fidelity and thermodynamic properties of Se-DNA in solution and in protein-DNA complexes.

Nucleic-acid structural deformability deduced from anisotropic displacement parameters

The growing numbers of very well resolved nucleic-acid crystal structures with anisotropic displacement parameters provide an unprecedented opportunity to learn about the natural motions of DNA and RNA. Here we report a new Monte-Carlo approach that takes direct account of this information to extract the distortions of covalent structure, base pairing, and dinucleotide geometry intrinsic to regularly organized double-helical molecules. We present new methods to test the validity of the anisotropic parameters and examine the apparent deformability of a variety of structures, including several A, B, and Z DNA duplexes, an AB helical intermediate, an RNA, a ligand-DNA complex, and an enzyme-bound DNA. The rigid-body parameters characterizing the positions of the bases in the structures mirror the mean parameters found when atomic motion is taken into account. The base-pair fluctuations intrinsic to a single structure, however, differ from those extracted from collections of nucleic-acid structures, although selected base-pair steps undergo conformational excursions along routes suggested by the ensembles. The computations reveal surprising new molecular insights, such as the stiffening of DNA and concomitant separation of motions of contacted nucleotides on opposite strands by the binding of Escherichia coli endonuclease VIII, which suggest how the protein may direct enzymatic action.

Revisiting the association of cationic groove-binding drugs to DNA using a Poisson-Boltzmann approach

Proper modeling of nonspecific salt-mediated electrostatic interactions is essential to understanding the binding of charged ligands to nucleic acids. Because the linear Poisson-Boltzmann equation (PBE) and the more approximate generalized Born approach are applied routinely to nucleic acids and their interactions with charged ligands, the reliability of these methods is examined vis-à-vis an efficient nonlinear PBE method. For moderate salt concentrations, the negative derivative, SK(pred), of the electrostatic binding free energy, DeltaG(el), with respect to the logarithm of the 1:1 salt concentration, [M(+)], for 33 cationic minor groove drugs binding to AT-rich DNA sequences is shown to be consistently negative and virtually constant over the salt range considered (0.1-0.4 M NaCl). The magnitude of SK(pred) is approximately equal to the charge on the drug, as predicted by counterion condensation theory (CCT) and observed in thermodynamic binding studies. The linear PBE is shown to overestimate the magnitude of SK(pred), whereas the nonlinear PBE closely matches the experimental results. The PBE predictions of SK(pred) were not correlated with DeltaG(el) in the presence of a dielectric discontinuity, as would be expected from the CCT. Because this correlation does not hold, parameterizing the PBE predictions of DeltaG(el) against the reported experimental data is not possible. Moreover, the common practice of extracting the electrostatic and nonelectrostatic contributions to the binding of charged ligands to biopolyelectrolytes based on the simple relation between experimental SK values and the electrostatic binding free energy that is based on CCT is called into question by the results presented here. Although the rigid-docking nonlinear PB calculations provide reliable predictions of SK(pred), at least for the charged ligand-nucleic acid complexes studied here, accurate estimates of DeltaG(el) will require further development in theoretical and experimental approaches.

On the role of Hoogsteen:Hoogsteen interactions in RNA: ab initio investigations of structures and energies

We use a combination of database analysis and quantum chemical studies to investigate the role of cis and trans Hoogsteen:Hoogsteen (H:H) base pairs and associated higher-order structures in RNA. We add three new examples to the list of previously identified base-pair combinations belonging to these families and, in addition to contextual classification and characterization of their structural and energetic features, we compare their interbase interaction energies and propensities toward participation in triplets and quartets. We find that some base pairs, which are nonplanar in their isolated minimum energy geometries, attain planarity and stability upon triplet formation. A:A H:H trans is the most frequent H:H combination in RNA structures. This base pair occurs at many distinct positions in known rRNA structures, where it helps in the interaction of ribosomal domains in the 50S subunit. It is also present as a part of tertiary interaction in tRNA structures. Although quantum chemical studies suggest an intrinsically nonplanar geometry for this base pair in isolated form, it has the tendency to attain planar geometry in RNA crystal structures by forming higher-order tertiary interactions or in the presence of additional base-phosphate interactions. The tendency of this base pair to form such additional interactions may be helpful in bringing together different segments of RNA, thus making it suitable for the role of facilitator for RNA folding. This also explains the high occurrence frequency of this base pair among all H:H interactions.

Pairing geometry of the hydrophobic thymine analogue 2,4-difluorotoluene in duplex DNA as analyzed by X-ray crystallography

Certain DNA polymerases (pols) were found to efficiently insert A opposite the hydrophobic T isostere 2,4-difluorotoluene (F) and vice versa, resulting in the widely held belief that some pols rely on shape rather than H-bonding for accurate replication. Using X-ray crystallography we have analyzed the geometry of F:A pairs in duplex DNA and observed a distance between fluorine and the exocyclic amino group of A that is consistent with a H-bond, thus challenging the assumption that the F analogue is unable to engage in H-bonding as well as the steric hypothesis of DNA replication. Therefore, shape and H-bonding are inherently related, and steric constraints at a pol active site, or conferred by stacking or the DNA backbone conformation, may enable H-bonding by F.