Hydrogen bonding of nucleotide base pairs: Application of the PM3 method

The ability of the pm3 semiempirical quantum mechanical method to reproduce hydrogen bonding in nucleotide base pairs was assessed. Results of pm3 calculations on the nucleotides 2′-deoxyadenosine 5′-monophosphate (pdA), 2′-deoxyguanosine 5′-monophosphate (pdG), 2′-deoxycytidine 5′-monophosphate (pdC), and 2′-deoxythymidine 5′-monophosphate (pdT) and the base pairs pdA–pdT, pdG–pdC, and pdG(syn)–pdC are presented and discussed. The pm3 method is the first of the parameterized nddo quantum mechanical models with any ability to reproduce hydrogen bonding between nucleotide base pairs. Intermolecular hydrogen bond lengths between nucleotides displaying Watson–Crick base pairing are 0.1–0.2 Å less than experimental results. Nucleotide bond distances, bond angles, and torsion angles about the glycosyl bond (χ), the C4′C5′ bond (γ), and the C5′O5′ bond (β) agree with experimental results. There are many possible conformations of nucleotides. pm3 calculations reveal that many of the most stable conformations are stabilized by intramolecular CHO hydrogen bonds. These interactions disrupt the usual sugar puckering. The stacking interactions of a dT–pdA duplex are examined at different levels of gradient optimization. The intramolecular hydrogen bonds found in the nucleotide base pairs disappear in the duplex, as a result of the additional constraints on the phosphate group when part of a DNA backbone. Sugar puckering is reproduced by the pm3 method for the four bases in the dT–pdA duplex. pm3 underestimates the attractive stacking interactions of base pairs in a B-DNA helical conformation. The performance of the pm3 method implemented in SPARTAN is contrasted with that implemented in MOPAC. At present, accurate ab initio calculations are too timeconsuming to be of practical use, and molecular mechanics methods cannot be used to determine quantum mechanical properties such as reaction-path calculations, transition-state structures, and activation energies. The pm3 method should be used with extreme caution for examination of small DNA systems. Future parameterizations of semiempirical methods should incorporate base stacking interactions into the parameterization data set to enhance the ability of these methods.

The structural basis of DNA binding by the single-stranded DNA-binding protein from Sulfolobus solfataricus

Canonical single-stranded DNA-binding proteins (SSBs) from the oligosaccharide/oligonucleotide-binding (OB) domain family are present in all known organisms and are critical for DNA replication, recombination and repair. The SSB from the hyperthermophilic crenarchaeote Sulfolobus solfataricus (SsoSSB) has a ‘simple’ domain organization consisting of a single DNA-binding OB fold coupled to a flexible C-terminal tail, in contrast with other SSBs in this family that incorporate up to four OB domains. Despite the large differences in the domain organization within the SSB family, the structure of the OB domain is remarkably similar all cellular life forms. However, there are significant differences in the molecular mechanism of ssDNA binding. We have determined the structure of the SsoSSB OB domain bound to ssDNA by NMR spectroscopy. We reveal that ssDNA recognition is modulated by base-stacking of three key aromatic residues, in contrast with the OB domains of human RPA and the recently discovered human homologue of SsoSSB, hSSB1. We also demonstrate that SsoSSB binds ssDNA with a footprint of five bases and with a defined binding polarity. These data elucidate the structural basis of DNA binding and shed light on the molecular mechanism by which these ‘simple’ SSBs interact with ssDNA.

Sequence effects in structures of the dinucleotides d-pApT AND d-pTpA

The H1',H2' and H2″ regions of the 270-MHz PMR spectra of two deoxydinucleotides, d-pTpA and d-pApT, have been analyzed. The coupling constants in the sugar ring indicate that both A and T sugars have a tendency to acquire 2E conformations. There is also a marginal difference in the 2E populations of the T sugar in the two dinucleotides. The trends in the chemical shifts of base protons indicate different stacking of the bases in d-pApT and d-pTpA. The sequence effects on base stacking and pentose conformation are discussed.

Interaction between the Z-type DNA duplex and 1,3-propanediamine:crystal structure of d(CACGTG)(2) at 1.2 angstrom resolution

The crystal structure of a hexamer duplex d(CACGTG)(2) has been determined and refined to an R-factor of 18.3% using X-ray data up to 1.2 angstrom resolution. The sequence crystallizes as a left-handed Z-form double helix with Watson-Crick base pairing. There is one hexamer duplex, a spermine molecule, 71 water molecules, and an unexpected diamine (Z-5, 1,3-propanediamine, C3H10N2)) in the asymmetric unit. This is the high-resolution non-disordered structure of a Z-DNA hexamer containing two AT base pairs in the interior of a duplex with no modifications such as bromination or methylation on cytosine bases. This structure does not possess multivalent cations such as cobalt hexaammine that are known to stabilize Z-DNA. The overall duplex structure and its crystal interactions are similar to those of the pure-spermine form of the d(CGCGCG)(2) structure. The spine of hydration in the minor groove is intact except in the vicinity of the T5A8 base pair. The binding of the Z-5 molecule in the minor grove of the d(CACGTG)(2) duplex appears to have a profound effect in conferring stability to a Z-DNA conformation via electrostatic complementarity and hydrogen bonding interactions. The successive base stacking geometry in d(CACGTG)(2) is similar to the corresponding steps in d(CG)(3). These results suggest that specific polyamines such as Z-5 could serve as powerful inducers of Z-type conformation in unmodified DNA sequences with AT base pairs. This structure provides a molecular basis for stabilizing AT base pairs incorporated into an alternating d(CG) sequence.

DNA-Mediated Charge Transport Devices for Protein Detection

Detection of biologically relevant targets, including small molecules, proteins, DNA, and RNA, is vital for fundamental research as well as clinical diagnostics. Sensors with biological elements provide a natural foundation for such devices because of the inherent recognition capabilities of biomolecules. Electrochemical DNA platforms are simple, sensitive, and do not require complex target labeling or expensive instrumentation. Sensitivity and specificity are added to DNA electrochemical platforms when the physical properties of DNA are harnessed. The inherent structure of DNA, with its stacked core of aromatic bases, enables DNA to act as a wire via DNA-mediated charge transport (DNA CT). DNA CT is not only robust over long molecular distances of at least 34 nm, but is also especially sensitive to anything that perturbs proper base stacking, including DNA mismatches, lesions, or DNA-binding proteins that distort the π-stack. Electrochemical sensors based on DNA CT have previously been used for single-nucleotide polymorphism detection, hybridization assays, and DNA-binding protein detection. Here, improvements to (i) the structure of DNA monolayers and (ii) the signal amplification with DNA CT platforms for improved sensitivity and detection are described.

First, improvements to the control over DNA monolayer formation are reported through the incorporation of copper-free click chemistry into DNA monolayer assembly. As opposed to conventional film formation involving the self-assembly of thiolated DNA, copper-free click chemistry enables DNA to be tethered to a pre-formed mixed alkylthiol monolayer. The total amount of DNA in the final film is directly related to the amount of azide in the underlying alkylthiol monolayer. DNA monolayers formed with this technique are significantly more homogeneous and lower density, with a larger amount of individual helices exposed to the analyte solution. With these improved monolayers, significantly more sensitive detection of the transcription factor TATA binding protein (TBP) is achieved.

Using low-density DNA monolayers, two-electrode DNA arrays were designed and fabricated to enable the placement of multiple DNA sequences onto a single underlying electrode. To pattern DNA onto the primary electrode surface of these arrays, a copper precatalyst for click chemistry was electrochemically activated at the secondary electrode. The location of the secondary electrode relative to the primary electrode enabled the patterning of up to four sequences of DNA onto a single electrode surface. As opposed to conventional electrochemical readout from the primary, DNA-modified electrode, a secondary microelectrode, coupled with electrocatalytic signal amplification, enables more sensitive detection with spatial resolution on the DNA array electrode surface. Using this two-electrode platform, arrays have been formed that facilitate differentiation between well-matched and mismatched sequences, detection of transcription factors, and sequence-selective DNA hybridization, all with the incorporation of internal controls.

For effective clinical detection, the two working electrode platform was multiplexed to contain two complementary arrays, each with fifteen electrodes. This platform, coupled with low density DNA monolayers and electrocatalysis with readout from a secondary electrode, enabled even more sensitive detection from especially small volumes (4 μL per well). This multiplexed platform has enabled the simultaneous detection of two transcription factors, TBP and CopG, with surface dissociation constants comparable to their solution dissociation constants.

With the sensitivity and selectivity obtained from the multiplexed, two working electrode array, an electrochemical signal-on assay for activity of the human methyltransferase DNMT1 was incorporated. DNMT1 is the most abundant human methyltransferase, and its aberrant methylation has been linked to the development of cancer. However, current methods to monitor methyltransferase activity are either ineffective with crude samples or are impractical to develop for clinical applications due to a reliance on radioactivity. Electrochemical detection of methyltransferase activity, in contrast, circumvents these issues. The signal-on detection assay translates methylation events into electrochemical signals via a methylation-specific restriction enzyme. Using the two working electrode platform combined with this assay, DNMT1 activity from tumor and healthy adjacent tissue lysate were evaluated. Our electrochemical measurements revealed significant differences in methyltransferase activity between tumor tissue and healthy adjacent tissue.

As differential activity was observed between colorectal tumor tissue and healthy adjacent tissue, ten tumor sets were subsequently analyzed for DNMT1 activity both electrochemically and by tritium incorporation. These results were compared to expression levels of DNMT1, measured by qPCR, and total DNMT1 protein content, measured by Western blot. The only trend detected was that hyperactivity was observed in the tumor samples as compared to the healthy adjacent tissue when measured electrochemically. These advances in DNA CT-based platforms have propelled this class of sensors from the purely academic realm into the realm of clinically relevant detection.

Studies of dinucleoside monophosphates and monomer-polynucleotide interactions by proton magnetic resonance

The nature of the intra- and intermolecular base-stacking interactions involving several dinucleoside monophosphates in aqueous solution have been investigated by proton magnetic resonance spectrosocopy, and this method has been applied to a study of the interaction of polyuridylic acid with purine and adenosine monomers.

The pmr spectra of adenylyl (3' → 5') cytidine (ApC) and cytidylyl (3' → 5') adenosine (CpA) have been studied as a function of concentration and temperature. The results of these studies indicate that the intramolecular base-stacking interactions between the adenine and cytosine bases of these dinucleoside monophosphates are rather strong, and that the stacking tendencies are comparable for the two sequence isomers. The chemical shifts of the cytosine H₅ and adenine H₂ protons, and their variations with temperature, were shown to be consistent with stacked conformations in which both bases of the dinucleoside monophosphates are preferentially oriented in the anti conformation as in similar dApdC, and dCpdA (dA = deoxyadenosine; dC = deoxycytidine) segments in double helical DNA. The intramolecular stacking interaction was found to have a pronounced effect on the conformations of the ribose moieties, and these conformational changes are discussed. The concentration studies indicate extensive self-association of these dinucleoside monophosphates, and analysis of the concentration data facilitated determination of the dimerization constant for the association process as well as the nature of the intermolecular complexes.

The dependence of the ribose conformation upon the extent of intramolecular base-stacking was used to demonstrate that the base-base interaction in cytidylyl (3' → 5') cytidine (CpC) is rather strong, while there appears to be little interaction between the two uracil bases of uridylyl (3' → 5') uridine (UpU).

Studies of the binding of purine to several ribose and deoxyribose dinucleoside monophosphates show that the mode of interaction is base-stacking, and evidence for the formation of a purine-dinucleoside monophosphate intercalated complex is presented. The purine proton resonances are markedly broadened in this complex, and estimates of the purine linewidths in the complex and the equilibrium constant for purine intercalation are obtained.

A study of the interaction of unsubstitued purine with polyuridylic acid at 29°C by pmr indicated that purine binds to the uracil bases of the polymer by base-stacking. The severe broadening of the purine proton resonances observed provides strong evidence for the intercalation of purine between adjacent uracil bases of poly U. This interaction does not result in a more rigid or ordered structure for the polymer.

Investigation of the interaction between adenosine and polyuridylic acid revealed two modes of interaction between the monomer and the polymer, depending on the temperature. At temperatures above 26°C or so, monomeric adenosine binds to poly U by noncooperative A-U base stacking. Below this temperature, a rigid triple-stranded 1A:2U complex is formed, presumably via cooperative hydrogen-bonding as has previously been reported.

These results clearly illustrate the importance of base-stacking in non-specific interactions between bases, nucleosides and nucleotides, and also reveal the important role of the base-stacking interactions in cooperatively for med structures involving specific base-pairing where both types of interaction are possible.

Proton magnetic resonance studies of ribonucleic acid complexes. I. Complexes of biological bases and oligonucleotides with RNA. II. Template recognition and the degeneracy of the genetic code

Part I. Complexes of Biological Bases and Oligonucleotides with RNA

The physical nature of complexes of several biological bases and oligonucleotides with single-stranded ribonucleic acids have been studied by high resolution proton magnetic resonance spectroscopy. The importance of various forces in the stabilization of these complexes is also discussed.

Previous work has shown that purine forms an intercalated complex with single-stranded nucleic acids. This complex formation led to severe and stereospecific broadening of the purine resonances. From the field dependence of the linewidths, T₁ measurements of the purine protons and nuclear Overhauser enhancement experiments, the mechanism for the line broadening was ascertained to be dipole-dipole interactions between the purine protons and the ribose protons of the nucleic acid.

The interactions of ethidium bromide (EB) with several RNA residues have been studied. EB forms vertically stacked aggregates with itself as well as with uridine, 3'-uridine monophosphate and 5'-uridine monophosphate and forms an intercalated complex with uridylyl (3' → 5') uridine and polyuridylic acid (poly U). The geometry of EB in the intercalated complex has also been determined.

The effect of chain length of oligo-A-nucleotides on their mode of interaction with poly U in D₂0 at neutral pD have also been studied. Below room temperatures, ApA and ApApA form a rigid triple-stranded complex involving a stoichiometry of one adenine to two uracil bases, presumably via specific adenine-uracil base pairing and cooperative base stacking of the adenine bases. While no evidence was obtained for the interaction of ApA with poly U above room temperature, ApApA exhibited complex formation of a 1:1 nature with poly U by forming Watson-Crick base pairs. The thermodynamics of these systems are discussed.

Part II. Template Recognition and the Degeneracy of the Genetic Code

The interaction of ApApG and poly U was studied as a model system for the codon-anticodon interaction of tRNA and mRNA in vivo. ApApG was shown to interact with poly U below ~20°C. The interaction was of a 1:1 nature which exhibited the Hoogsteen bonding scheme. The three bases of ApApG are in an anti conformation and the guanosine base appears to be in the lactim tautomeric form in the complex.

Due to the inadequacies of previous models for the degeneracy of the genetic code in explaining the observed interactions of ApApG with poly U, the "tautomeric doublet" model is proposed as a possible explanation of the degenerate interactions of tRNA with mRNA during protein synthesis in vivo.

Comparative studies on parallel and antiparallel duplex and triplex DNA

Parallel strand models for base sequences d(A)(10). d(T)(10), d(AT)(5) . d(TA)(5), d(G(5)C(5)). d(C(5)G(5)), d(GC)(5) . d(CG)(5) and d(CTATAGGGAT). d(GATATCCCTA), where reverse Watson-Crick A-T pairing with two H-bonds and reverse Watson-Crick G-C pairing with one H-bond or with two H-bonds were adopted, and three models of d(T)(14). d(A)(14). d(T)(14) triple helix with different strand orientations were built up by molecular architecture and energy minimization. Comparisons of parallel duplex models with their corresponding B-DNA models and comparisons among the three triple helices showed: (i) conformational energies of parallel AT duplex models were a little lower, while for GC duplex models they were about 8% higher than that of their corresponding B-DNA models; (ii) the energy differences between parallel and B-type duplex models and among the three triple helices arose mainly from base stacking energies, especially for GC base pairing; (iii) the parallel duplexes with one H-bond G-C pairs were less stable than those with two H-bonds G-C pairs. The present paper includes a brief discussion about the effect of base stacking and base sequences on DNA conformations. (C) 1997 Academic Press Limited.

Spectral and electrochemical detection of protonated triplex formation by a small-molecule anticancer agent

Triplex helical formation has been the focus of considerable interest because of possible applications in developing new molecular biology tools as well as therapeutic agents and the possible relevance of H-DNA structures in biology system. We report here that a small-molecule anticancer agent, coralyne, has binding preference to the less stable protonated triplex d(C+-T)(6):d(A-G)(6).d(C-T)(6) over duplex d(A-G)(6).d(C-T)(6) and shows different spectral and electrochemical characteristics when binding to triplex and duplex DNA, indicating that electrochemical technique can detect the less stable protonated triplex formation.

Effect of DNA flanking sequence on charge transport in short DNA duplexes

Several factors can influence charge transport (CT)-mediated DNA, such as sequence, distance, base stacking, base pair mismatch, conformation, tether length, etc. However, the DNA context effect or how flanking sequences influence redox active drugs in the DNA CT reaction and later in DNA enzymatic repair and synthesis is still not well understood. The set of seven DNA molecules in this study have been characterized well for the study of flanking sequence effects. These DNA duplexes are formed from self-complementary strands and contain the common central four-base sequence 5'-A-G-C-T-3', flanked on both sides by either (AT)(n) or (AA)(n) (n = 2, 3, or 4) or AA(AT)(2). UV-vis, fluorescence, UV melting, circular dichroism, and cyclic voltammetry experiments were used to study the flanking sequence effect on CT-mediated DNA by using daunomycin or adriamycin cross-linked with these seven DNA molecules. Our results showed that charge transport was related to the flanking sequence, DNA melting free energy, and ionic strength. For (AA)(n) or (AT)(n) species of the same length, (AA)(n) series were more stable and more efficient CT was observed through the (AA)(n) series. The same trend was observed for (AA)(n) and (AT)(n) series at different ionic strengths, further supporting the idea that flanking sequence can result in different base stacking and modulate charge transport through these seven DNA molecules.

Caractérisation structurale et thermodynamique de la reconnaissance du substrat par le ribozyme VS de Neurospora

Les interactions ARN/ARN de type kissing-loop sont des éléments de structure tertiaire qui jouent souvent des rôles clés chez les ARN, tant au niveau fonctionnel que structural. En effet, ce type d’interaction est crucial pour plusieurs processus dépendant des ARN, notamment pour l’initiation de la traduction, la reconnaissance des ARN antisens et la dimérisation de génome rétroviral. Les interactions kissing-loop sont également importantes pour le repliement des ARN, puisqu’elles permettent d’établir des contacts à longue distance entre différents ARN ou encore entre les domaines éloignés d’un même ARN. Ce type d’interaction stabilise aussi les structures complexes des ARN fonctionnels tels que les ARNt, les riborégulateurs et les ribozymes. Comme d’autres ARN fonctionnels, le ribozyme VS de Neurospora contient une interaction kissing-loop importante. Celle-ci est impliquée dans la reconnaissance du substrat et se forme entre la tige-boucle I (stem-loop I, SLI) du substrat et la tige-boucle V (stem-loop V, SLV) du domaine catalytique. Des études biochimiques ont démontré que l’interaction kissing-loop I/V, dépendante du magnésium, implique trois paires de bases Watson-Crick (W-C). De plus, cette interaction est associée à un réarrangement de la structure du substrat, le faisant passer d’une conformation inactive dite unshifted à une conformation active dite shifted. Les travaux présentés dans cette thèse consistent en une caractérisation structurale et thermodynamique de l’interaction kissing-loop I/V du ribozyme VS, laquelle est formée de fragments d’ARN représentant les tige-boucles I et V dérivées du ribozyme VS (SLI et SLV). Cette caractérisation a été réalisée principalement par spectroscopie de résonance magnétique nucléaire (RMN) et par titrage calorimétrique isotherme (isothermal titration calorimetry, ITC) en utilisant différents complexes SLI/SLV dans lesquels l’ARN SLV est commun à tous les complexes, alors que différentes variations de l’ARN SLI ont été utilisées, soit en conformation shiftable ou preshifted. Les données d’ITC ont permis de démontrer qu’en présence d’une concentration saturante de magnésium, l’affinité d’un substrat SLI preshifted pour SLV est extrêmement élevée, rendant cette interaction plus stable que ce qui est prédit pour un duplexe d’ARN équivalent. De plus, l’étude effectuée par ITC montre que des ARN SLI preshifted présentent une meilleure affinité pour SLV que des ARN SLI shiftable, ce qui a permis de calculer le coût énergétique associé au réarrangement de structure du substrat. En plus de confirmer la formation des trois paires de bases W-C prédites à la jonction I/V, les études de RMN ont permis d’obtenir une preuve structurale directe du réarrangement structural des substrats SLI shiftable en présence de magnésium et de l’ARN SLV. La structure RMN d’un complexe SLI/SLV de grande affinité démontre que les boucles terminales de SLI et SLV forment chacune un motif U-turn, ce qui facilite l’appariement W-C intermoléculaire. Plusieurs autres interactions ont été définies à l’interface I/V, notamment des triplets de bases, ainsi que des empilements de bases. Ces interactions contribuent d’ailleurs à la création d’une structure présentant un empilement continu, c’est-à-dire qui se propage du centre de l’interaction jusqu’aux bouts des tiges de SLI et SLV. Ces études de RMN permettent donc de mieux comprendre la stabilité exceptionnelle de l’interaction kissing-loop I/V au niveau structural et mènent à l’élaboration d’un modèle cinétique de l’activation du substrat par le ribozyme VS. En considérant l’ensemble des données d’ITC et de RMN, l’étonnante stabilité de l’interaction I/V s’explique probablement par une combinaison de facteurs, dont les motifs U-turn, la présence d’un nucléotide exclu de la boucle de SLV (U700), la liaison de cations magnésium et l’empilement de bases continu à la jonction I/V.

Long-lived excited-state dynamics of i-motif structures probed by time-resolved infrared spectroscopy

UV-generated excited states of cytosine (C) nucleobases are precursors to mutagenic photoproduct formation. The i-motif formed from C-rich sequences is known to exhibit high yields of long-lived excited states following UV absorption. Here the excited states of several i-motif structures have been characterized following 267 nm laser excitation using time-resolved infrared spectroscopy (TRIR). All structures possess a long-lived excited state of ~300 ps and notably in some cases decays greater than 1 ns are observed. These unusually long-lived lifetimes are attributed to the interdigitated DNA structure which prevents direct base stacking overlap.

From Rational Design of Drug Crystals to Understanding of Nucleic Acid Structures: Lamivudine Duplex

A DNA-like duplex of nucleosides is probable to exist even without the 5`-phosphate groups needed to assemble the chain backbone. However, double-stranded helical structures of nucleosides are unknown. Here, we report a duplex of nucleoside analogs that is spontaneously assembled due to stacking of the neutral and protonated molecules of lamivudine, a nucleoside reverse transcriptase inhibitor (NTRI) widely used in anti-HIV drug combinatory medication. The left-handed lamivudine duplex has features similar to those of i-motif DNA, as the face-to-face base stacking and the helix rise per base pair. Furthermore, the protonation pattern on alternate bases expected for it DNA-like duplex stabilized by pairing of neutral and protonated cytosine fragments was observed for the first time in the lamivudine double-stranded helix. This structure demonstrates that hydrogen bonds can substitute for covalent phosphodiester linkage in the stabilization of the duplex backbone. This interesting example of spontaneous molecular self-organization indicates that the 5`-phosphate group could not be a requirement for duplex assembly.

Further Quantum Mechanical Evidence that Difluorotoluene does not Hydrogen Bond

Calculations were run on the methylated DNA base pairs adenine:thymine and adenine:difluorotoluene to further investigate the hydrogen-bonding properties of difluorotoluene (F). Geometries were optimized using hybrid density functional theory. Single-point calculations at the MP2(full) level were performed to obtain more rigorous energies. The functional counterpoise method was used to correct for the basis set superposition error (BSSE), and the interaction energies were also corrected for fragment relaxation. These corrections brought the B3LYP and MP2 interaction energies into excellent agreement. In the gas phase, the Gibbs free energies calculated at the B3LYP and MP2 levels of theory predict that A and T will spontaneously form an A:T pair while A:F spontaneously dissociates into A and F. Solvation effects on the pairing of the bases were explored using implicit solvent models for water and chloroform. In aqueous solution, both A:T and A:F are predicted to dissociate into their component monomers. Semiempirical calculations were performed on small sections of B-form DNA containing the two pairs, and the results provide support for the concept that base stacking is more important than hydrogen bonding for the stability of the A:F pair within a DNA helix.

Molecular Dynamics Simulation of a PNA·DNA·PNA Triple Helix in Aqueous Solution

Molecular dynamics simulations have been used to explore the conformational flexibility of a PNA·DNA·PNA triple helix in aqueous solution. Three 1.05 ns trajectories starting from different but reasonable conformations have been generated and analyzed in detail. All three trajectories converge within about 300 ps to produce stable and very similar conformational ensembles, which resemble the crystal structure conformation in many details. However, in contrast to the crystal structure, there is a tendency for the direct hydrogen-bonds observed between the amide hydrogens of the Hoogsteen-binding PNA strand and the phosphate oxygens of the DNA strand to be replaced by water-mediated hydrogen bonds, which also involve pyrimidine O2 atoms. This structural transition does not appear to weaken the triplex structure but alters groove widths and so may relate to the potential for recognition of such structures by other ligands (small molecules or proteins). Energetic analysis leads us to conclude that the reason that the hybrid PNA/DNA triplex has quite different helical characteristics from the all-DNA triplex is not because the additional flexibility imparted by the replacement of sugar−phosphate by PNA backbones allows motions to improve base-stacking but rather that base-stacking interactions are very similar in both types of triplex and the driving force comes from weak but definate conformational preferences of the PNA strands.