Base Pairs In Duplex DNA Research Articles

Sequence Dependence in Nucleosome Dynamics.

The basic packaging unit of eukaryotic chromatin is the nucleosome that contains 145-147 base pair duplex DNA wrapped around an octameric histone protein. While the DNA sequence plays a crucial role in controlling the positioning of the nucleosome, the molecular details behind the interplay between DNA sequence and nucleosome dynamics remain relatively unexplored. This study analyzes this interplay in detail by performing all-atom molecular dynamics simulations of nucleosomes, comparing the human α-satellite palindromic (ASP) and the strong positioning "Widom-601" DNA sequence at time scales of 12 μs. The simulations are performed at salt concentrations 10-20 times higher than physiological salt concentrations to screen the electrostatic interactions and promote unwrapping. These microsecond-long simulations give insight into the molecular-level sequence-dependent events that dictate the pathway of DNA unwrapping. We find that the "ASP" sequence forms a loop around SHL ± 5 for three sets of simulations. Coincident with loop formation is a cooperative increase in contacts with the neighboring N-terminal H2B tail and C-terminal H2A tail and the release of neighboring counterions. We find that the Widom-601 sequence exhibits a strong breathing motion of the nucleic acid ends. Coincident with the breathing motion is the collapse of the full N-terminal H3 tail and formation of an α-helix that interacts with the H3 histone core. We postulate that the dynamics of these histone tails and their modification with post-translational modifications (PTMs) may play a key role in governing this dynamics.

Importance of isothermal titration calorimetry for the detection of the direct binding of metal ions to mismatched base pairs in duplex DNA.

Metal ion-nucleic acid interactions contribute substantially to the structure and biological activity of nucleic acids and have a wide range of potential applications in nanotechnology. In this study, we examined the interactions between metal ions and mismatched base pairs in duplex DNA to reveal the underlying molecular mechanism. UV melting analyses showed that the melting temperature (Tm) of a 21-base pair duplex DNA with each of the C-A, C-C and C-T mismatched base pairs increased upon the addition of Ag+. However, isothermal titration calorimetry (ITC) demonstrated that Ag+ only bound to the C-C mismatched base pair of the duplex DNA to form C-Ag-C bonds, without binding to the C-A and C-T mismatches. These results showed that Tm increased even when metal ions did not bind to the mismatched base pairs of the duplex DNA. Although the increase in Tm upon the addition of the metal ions is often used to detect metal ion binding to mismatched base pairs of duplex DNA, these results indicated that UV melting analyses are unable to detect the direct binding of metal ions to the mismatched base pairs. Because ITC analyses directly detect the heat derived from metal ion binding to mismatched base pairs of duplex DNA, we concluded that this may be an effective detection approach.

Hoogsteen Base Pairings

In contrast to Watson-Crick (WC) base pairing, Hoogsteen (HG) base pairing involves flipping a purine base 180° between its anti and syn conformation. Recent studies have shown that HG pairs coexist in dynamical equilibrium, and several biological functions depend on them. This significance has stirred computational research on this base-pairing transition. However, a methodical reproduction of sequence variations has continued to be out of reach. It is becoming increasingly clear that Hoogsteen base pairs play a crucial role in DNA replication, recognition, damage repair, and incorrect sequence repair. The Protein Data Bank contains a variety of Hoogsteen base pairing modes that include the preference for A–T versus G–C bps, TA versus GG steps, and a preference for 5'-purines at terminal ends. RNA A-form duplexes are strongly disfavored by Hoogsteen base pairs, in stark contrast to B-form DNA. Therefore, N1-methyl adenosine and N1-methyl guanosine, which transpire in DNA as alkylation impairment and in RNA as posttranscriptional adjustments, have great differences in effects. They create G–C+ and A–U Hoogsteen base pairs in duplex DNA that preserve the structural integrity of the double helix, but obstruct base pairing altogether and induce local duplex melting in RNA, providing a mechanism for potently disrupting RNA structure through posttranscriptional modifications. In duplex DNA, they maintain the structural integrity of the double helix by creating G–C+ and A–U Hoogsteen base pairs, but block base pairing altogether and cause local duplex melting in RNA, thus providing a potent means for disrupting RNA structure post transcriptionally. As a result of the markedly different inclinations for B-DNA and A-RNA to form Hoogsteen base pairs, they may be able to balance the opposing demands of maintaining genome stability and dynamically modulating the epitranscriptome. This review examines the occurrence of Hoogsteen base pairs in DNA and RNA duplexes.

Recognition of 5-methyl-CG and CG base pairs in duplex DNA with high stability using antiparallel-type triplex-forming oligonucleotides with 2-guanidinoethyl-2'-deoxynebularine.

The formation of triplex DNA is a site-specific recognition method that directly targets duplex DNA. However, triplex DNA formation is generally formed for the GC and AT base pairs of duplex DNA, and there are no natural nucleotides that recognize the CG and TA base pairs, or even the 5-methyl-CG (5mCG) base pair. Moreover, duplex DNA, including 5mCG base pairs, epigenetically regulates gene expression in vivo, and thus targeting strategies are of biological importance. Therefore, the development of triplex-forming oligonucleotides (TFOs) with artificial nucleosides that selectively recognize these base pairs with high affinity is needed. We recently reported that 2'-deoxy-2-aminonebularine derivatives exhibited the ability to recognize 5mCG and CG base pairs in triplex formation; however, this ability was dependent on sequences. Therefore, we designed and synthesized new nucleoside derivatives based on the 2'-deoxy-nebularine (dN) skeleton to shorten the linker length connecting to the hydrogen-bonding unit in formation of the antiparallel motif triplex. We successfully demonstrated that TFOs with 2-guanidinoethyl-2'-deoxynebularine (guanidino-dN) recognized 5mCG and CG base pairs with very high affinity in all four DNA sequences with different adjacent nucleobases of guanidino-dN as well as in the promoter sequences of human genes containing 5mCG base pairs with a high DNA methylation frequency.

Structural basis of the XPB helicase-Bax1 nuclease complex interacting with the repair bubble DNA.

Nucleotide excision repair (NER) removes various DNA lesions caused by UV light and chemical carcinogens. The DNA helicase XPB plays a key role in DNA opening and coordinating damage incision by nucleases during NER, but the underlying mechanisms remain unclear. Here, we report crystal structures of XPB from Sulfurisphaera tokodaii (St) bound to the nuclease Bax1 and their complex with a bubble DNA having one arm unwound in the crystal. StXPB and Bax1 together spirally encircle 10 base pairs of duplex DNA at the double-/single-stranded (ds–ss) junction. Furthermore, StXPB has its ThM motif intruding between the two DNA strands and gripping the 3′-overhang while Bax1 interacts with the 5′-overhang. This ternary complex likely reflects the state of repair bubble extension by the XPB and nuclease machine. ATP binding and hydrolysis by StXPB could lead to a spiral translocation along dsDNA and DNA strand separation by the ThM motif, revealing an unconventional DNA unwinding mechanism. Interestingly, the DNA is kept away from the nuclease domain of Bax1, potentially preventing DNA incision by Bax1 during repair bubble extension.

Single-molecule level structural dynamics of DNA unwinding by human mitochondrial Twinkle helicase

Knowledge of the molecular events in mitochondrial DNA (mtDNA) replication is crucial to understanding the origins of human disorders arising from mitochondrial dysfunction. Twinkle helicase is an essential component of mtDNA replication. Here, we employed atomic force microscopy imaging in air and liquids to visualize ring assembly, DNA binding, and unwinding activity of individual Twinkle hexamers at the single-molecule level. We observed that the Twinkle subunits self-assemble into hexamers and higher-order complexes that can switch between open and closed-ring configurations in the absence of DNA. Our analyses helped visualize Twinkle loading onto and unloading from DNA in an open-ringed configuration. They also revealed that closed-ring conformers bind and unwind several hundred base pairs of duplex DNA at an average rate of ∼240 bp/min. We found that the addition of mitochondrial single-stranded (ss) DNA-binding protein both influences the ways Twinkle loads onto defined DNA substrates and stabilizes the unwound ssDNA product, resulting in a ∼5-fold stimulation of the apparent DNA-unwinding rate. Mitochondrial ssDNA-binding protein also increased the estimated translocation processivity from 1750 to >9000 bp before helicase disassociation, suggesting that more than half of the mitochondrial genome could be unwound by Twinkle during a single DNA-binding event. The strategies used in this work provide a new platform to examine Twinkle disease variants and the core mtDNA replication machinery. They also offer an enhanced framework to investigate molecular mechanisms underlying deletion and depletion of the mitochondrial genome as observed in mitochondrial diseases.

Investigating Nucleosome Stability via FRET

DNA's susceptibility to oxidative damage leads to the formation of the mutagenic Spiroiminodihydantoin (Sp) lesion. The Sp lesion contains a chiral center that disrupts the DNA double helix. When the Sp lesion is not repaired, it may lead to permanent mutations, apoptosis, cellular aging, or cancer. Previous work in our lab suggests that the Sp lesion can cause changes in nucleosome stability and DNA positioning. A nucleosome is the fundamental level of DNA packaging into chromosomes and forms when 146 base pair DNA duplex wraps around a histone octamer. This structure regulates access to DNA by repair proteins, transcription factors, and polymerases. With the use of Florescence Resonance Energy Transfer (FRET), we can examine how the Sp lesion affects nucleosome stability. We synthesized 146-mer duplexes with Cy3 and Cy5 fluorophores with Sp lesions at three potential locations. Then, we assembled nucleosome core particles with fluorescent donor-acceptor pairs on the DNA. FRET is a function of the proximity of the donor and acceptor fluorophores. Thus, when nucleosomes are gradually titrated with NaCl, DNA ends begin peeling and leads to a loss in acceptor fluorescence. With FRET, we are quantitatively measuring the correlation between nucleosome stability and location of the Sp lesion. We hope to advance our understanding of how the Sp lesion influences DNA compaction into nucleosomes.

Direct evidence for (G)O6···H2-N4(C)+ hydrogen bonding in transient G(syn)-C+ and G(syn)-m5C+ Hoogsteen base pairs in duplex DNA from cytosine amino nitrogen off-resonance R1ρ relaxation dispersion measurements.

NMR relaxation dispersion studies have shown that Watson-Crick G-C and A-T base pairs in duplex DNA exist in dynamic equilibrium with their Hoogsteen counterparts. Hoogsteen base pairs form through concurrent rotation of the purine base about the glycosidic bond from an anti to a syn conformation and constriction of the C1'-C1' distance across the base pair by ∼2 Å to allow Hoogsteen type hydrogen bonding. Owing to their unique structure, Hoogsteen base pairs can play important roles in DNA recognition, the accommodation, recognition, and repair of DNA damage, and in DNA replication. NMR relaxation dispersion experiments targeting imino nitrogen and protonated base and sugar carbons have provided insights into many structural features of transient Hoogsteen base pairs, including one of two predicted hydrogen bonds involving (G)N7···H-N3(C)+ and (A)N7···H-N3(T). Here, through measurement of cytosine amino (N4) R1ρ relaxation dispersion, we provide direct evidence for the second (G)O6···H2-N4(C)+ hydrogen bond in G(syn)-C+ transient Hoogsteen base pairs. The utility of cytosine N4 R1ρ relaxation dispersion as a new sensitive probe of transient Hoogsteen base pairs, and cytosine dynamics in general, is further demonstrated by measuring G(syn)-C+ Hoogsteen exchange near neutral pH and in the context of the naturally occurring DNA modification 5-methyl cytosine (m5C), in DNA samples prepared using chemical synthesis and a 15N labeled m5C phosphoramidite.

Half-Intercalation Stabilizes Slipped Mispairing and Explains Genome Vulnerability to Frameshift Mutagenesis by Endogenous "Molecular Bookmarks".

Some 60 years ago chemicals that intercalate between base pairs of duplex DNA were found to amplify frameshift mutagenesis. Surprisingly, the robust induction of frameshifts by intercalators still lacks a mechanistic model, leaving this classic phenomenon annoyingly intractable. A promising idea of asymmetric half-intercalation-stabilizing frameshift intermediates during DNA synthesis has never been developed into a model. Instead, researchers of frameshift mutagenesis embraced the powerful slipped-mispairing concept that unexpectedly struggled with the role of intercalators in frameshifting. It is proposed that the slipped mispairing and the half-intercalation ideas are two sides of the same coin. Further, existing findings are reviewed to test predictions of the combined "half-intercalation into the slipped-mispairing intermediate" model against accumulated knowledge. The existence of potential endogenous intercalators and the phenomenon of "DNA bookmarks" reveal ample possibilities for natural frameshift mutagenisis in the cell. From this alarming perspective, it is discussed how the cell could prevent genome deterioration from frameshift mutagenesis.

Tracking Photoinduced Charge Separation in DNA: from Start to Finish.

The initial studies of the dynamics of photoinduced charge separation conducted in our laboratories 20 years ago found strongly distance-dependent rate constants over short distances but failed to detect intermediates in the transport of positive charge (holes). These observations were consistent with the single-step superexchange or tunneling mechanism that had been observed for numerous donor-bridge-acceptor systems at that time. Subsequent studies found weak distance dependence for hole transport over longer distances in DNA, characteristic of incoherent hopping of either localized or delocalized holes. The introduction of synthetic DNA capped hairpin constructs possessing hole donor and acceptor chromophores (or purine bases) at opposite ends of a base-pair domain made it possible to determine the time required for transit of charge from one chromophore to the other and, in some cases, to distinguish between the transit time and the much faster initial charge injection time. These studies eliminated conventional tunneling as a viable mechanism for charge transport in DNA except at very short donor-acceptor separations; however, they did not establish the presence or nature of intermediates in the charge separation process. Recent studies in our laboratories have succeeded in identifying key intermediates as well as untangling the dynamics and efficiency of the charge separation process from start to finish. The dynamics of the initial charge injection process is dependent upon both its free energy and the stacking of the hole donor chromophore and adjacent purine base. The transport of positive charge (holes) over multiple base pairs in duplex DNA occurs most efficiently via repeating adenine bases, known as A-tracts. The transit time across an A-tract is strongly dependent upon the free energy for hole injection, whereas the efficiency of charge separation depends on the competition between charge delocalization and charge recombination in the contact radical ion pair. The guanine cation radical has been detected both by femtosecond transient absorption and by stimulated Raman spectroscopies when the guanine is located near the chromophore employed for hole injection into an A-tract. Replacement of guanine by its derivative 8-phenylethynylguanine (EG), permits tracking of hole transport across longer poly(purine) sequences as a consequence of the stronger transient absorption and stimulated Raman scattering for EG+• vs G+•. We have recently obtained evidence based on femtosecond transient absorption spectroscopy for the formation of delocalized A-polarons in A-tracts possessing four or more A-T base pairs. Similar methods have been used to track hole transport across less-common DNA structures including diblock and triblock poly(purines), locked nucleic acids, three-way junctions, and G-quadruplexes. Similar methods are have been applied to the study of photoinduced electron transport in DNA.

Atomic structures of excited state A-T Hoogsteen base pairs in duplex DNA by combining NMR relaxation dispersion, mutagenesis, and chemical shift calculations.

NMR relaxation dispersion studies indicate that in canonical duplex DNA, Watson-Crick base pairs (bps) exist in dynamic equilibrium with short-lived low abundance excited state Hoogsteen bps. N1-methylated adenine (m1A) and guanine (m1G) are naturally occurring forms of damage that stabilize Hoogsteen bps in duplex DNA. NMR dynamic ensembles of DNA duplexes with m1A-T Hoogsteen bps reveal significant changes in sugar pucker and backbone angles in and around the Hoogsteen bp, as well as kinking of the duplex towards the major groove. Whether these structural changes also occur upon forming excited state Hoogsteen bps in unmodified duplexes remains to be established because prior relaxation dispersion probes provided limited information regarding the sugar-backbone conformation. Here, we demonstrate measurements of C3' and C4' spin relaxation in the rotating frame (R1ρ) in uniformly 13C/15N labeled DNA as sensitive probes of the sugar-backbone conformation in DNA excited states. The chemical shifts, combined with structure-based predictions using an automated fragmentation quantum mechanics/molecular mechanics method, show that the dynamic ensemble of DNA duplexes containing m1A-T Hoogsteen bps accurately model the excited state Hoogsteen conformation in two different sequence contexts. Formation of excited state A-T Hoogsteen bps is accompanied by changes in sugar-backbone conformation that allow the flipped syn adenine to form hydrogen-bonds with its partner thymine and this in turn results in overall kinking of the DNA toward the major groove. Results support the assignment of Hoogsteen bps as the excited state observed in canonical duplex DNA, provide an atomic view of DNA dynamics linked to formation of Hoogsteen bps, and lay the groundwork for a potentially general strategy for solving structures of nucleic acid excited states.

Monitoring the Intercalation of Acridine Orange (AO) Molecules into 6-Methyl isoxanthopterin (6-MI)-Labeled DNA by Circular Dichroism and Single-Molecule FRET

Intercalation into duplex DNA occurs when a planar and aromatic small molecule ligand binds between adjacent base pairs of duplex DNA, forming a sandwich-like structure. This insertion significantly perturbs the double-stranded DNA structure. As a result, such intercalation may also perturb DNA-structure-dependent biological processes, such as DNA transcription, replication and repair. Based on such functional perturbations, extensive studies have been performed that focus on the development of intercalative drugs with possible clinical potential. Among the intercalators that have been examined, acridine derivatives have been much investigated. However, more work is required to develop effective drugs with the desired selectivities and reduced side effects. In an attempt to provide more incisive molecular approaches to understanding intercalation mechanisms, we have shown that 6-methyl isoxanthopterin (6-MI) and acridine orange (AO) can act as a donor-acceptor Forster Resonance Energy Transfer (FRET) pair. 6-MI is a fluorescent guanine analogue that can serve as a local reporter of DNA conformational changes, while acridine orange is a well-studied DNA intercalator. We show how these two optical probes can be paired to monitor local AO-DNA intercalation dynamics with multi-millisecond time resolution. Applying single-molecule FRET and linear dichroism methods developed in our group, we are using these approaches to study the roles of intermediate conformational states during the intercalation process.

M(1)A and m(1)G disrupt A-RNA structure through the intrinsic instability of Hoogsteen base pairs.

The B-DNA double helix can dynamically accommodate G–C and A–T base pairs in either Watson-Crick or Hoogsteen configurations. Here, we show that G–C+ and A–U Hoogsteen base pairs are strongly disfavored in A-RNA. As a result, N1-methyl adenosine and N1-methyl guanosine, which occur in DNA as a form of alkylation damage, and in RNA as a posttranscriptional modification, have dramatically different consequences. They create G–C+ and A–U Hoogsteen base pairs in duplex DNA that maintain the structural integrity of the double helix, but block base pairing all together and induce local duplex melting in RNA, providing a mechanism for potently disrupting RNA structure through posttranscriptional modifications. The markedly different propensities to form Hoogsteen base pairs in B-DNA and A-RNA may help meet the opposing requirements of maintaining genome stability on one hand, and dynamically modulating the structure of the epitranscriptome on the other.

Extracting water and ion distributions from solution x-ray scattering experiments.

Small-angle X-ray scattering measurements can provide valuable information about the solvent environment around biomolecules, but it can be difficult to extract solvent-specific information from observed intensity profiles. Intensities are proportional to the square of scattering amplitudes, which are complex quantities. Amplitudes in the forward direction are real, and the contribution from a solute of known structure (and from the waters it excludes) can be estimated from theory; hence, the amplitude arising from the solvent environment can be computed by difference. We have found that this "square root subtraction scheme" can be extended to non-zero q values, out to 0.1 Å(-1) for the systems considered here, since the phases arising from the solute and from the water environment are nearly identical in this angle range. This allows us to extract aspects of the water and ion distributions (beyond their total numbers), by combining experimental data for the complete system with calculations for the solutes. We use this approach to test molecular dynamics and integral-equation (3D-RISM (three-dimensional reference interaction site model)) models for solvent structure around myoglobin, lysozyme, and a 25 base-pair duplex DNA. Comparisons can be made both in Fourier space and in terms of the distribution of interatomic distances in real space. Generally, computed solvent distributions arising from the MD simulations fit experimental data better than those from 3D-RISM, even though the total small-angle X-ray scattering patterns are very similar; this illustrates the potential power of this sort of analysis to guide the development of computational models.

Hidden Hoogsteens in the Data

Although extremely rare in the crystallographic database, Hoogsteen (HG) base pairs have been modeled, especially at duplex termini, protein-DNA interfaces, damaged DNA and DNA mismatches. Formation of HG base pairs in duplex DNA requires a 180 degree rotation along the glycosidic bond of the purine relative to the canonical Watson-Crick (WC) base-pair. In a recent survey of the PDB we identified 75 HS base pairs in 2,910 x-ray derived models. Recent evidence suggests that the prevalence of HG base pairs in nature to be greater than previously thought. Due to electron density density ambiguity at mid to low resolution, it is possible for bases to be modeled incorrectly. We have developed a high-throughput method which examines the difference density around DNA to identify mis-modeled bases. Here we show where bases were incorrectly modeled and our subsequent fixes which show pairs modeled as WC that really are HG, vice versa, and evidence of WC-HG alternates. This research shows that base pairing conformations in DNA are more heterogeneous than previously thought.

A historical account of hoogsteen base‐pairs in duplex DNA

In 1957, a unique pattern of hydrogen bonding between N3 and O4 on uracil and N7 and N6 on adenine was proposed to explain how poly(rU) strands can associate with poly(rA)-poly(rU) duplexes to form triplexes. Two years later, Karst Hoogsteen visualized such a noncanonical A-T base-pair through X-ray analysis of co-crystals containing 9-methyladenine and 1-methylthymine. Subsequent X-ray analyses of guanine and cytosine derivatives yielded the expected Watson-Crick base-pairing, but those of adenine and thymine (or uridine) did not yield Watson-Crick base-pairs, instead favoring "Hoogsteen" base-pairing. More than two decades ensued without experimental "proof" for A-T Watson-Crick base-pairs, while Hoogsteen base-pairs continued to surface in AT-rich sequences, closing base-pairs of apical loops, in structures of DNA bound to antibiotics and proteins, damaged and chemically modified DNA, and in polymerases that replicate DNA via Hoogsteen pairing. Recently, NMR studies have shown that base-pairs in duplex DNA exist as a dynamic equilibrium between Watson-Crick and Hoogsteen forms. There is now little doubt that Hoogsteen base-pairs exist in significant abundance in genomic DNA, where they can expand the structural and functional versatility of duplex DNA beyond that which can be achieved based only on Watson-Crick base-pairing. Here, we provide a historical account of the discovery and characterization of Hoogsteen base-pairs, hoping that this will inform future studies exploring the occurrence and functional importance of these alternative base-pairs.

Sensing of a nucleic acid binding protein via a label-free perylene probe fluorescence recovery assay

A novel label-free fluorescence recovery assay for the sensing of a DNA binding protein has been developed. A transcription factor c-Jun protein, and a 21 base pair duplex DNA containing the c-Jun protein binding site (J-DNA) were selected. J-DNA was mixed with a cationic fluorescent perylene probe (compound 1), and induced aggregation of the probe. Quenching of the probe's fluorescence was observed. However, when c-Jun protein was mixed with the J-DNA, c-Jun bound to the duplex DNA, which reduced the degree of the induced perylene probe aggregation, and a turn on fluorescence signal was observed. The recovered fluorescence intensity was directly related to the amount of c-Jun added. The method is highly selective, six non-DNA binding proteins and one randomly selected 21 base pair duplex DNA (con-1) were tested. No noticeable compound 1 fluorescence recovery was observed. Mutations were also introduced to the c-Jun recognition sequence and much reduced fluorescence recovery was observed. Our assay is label-free, convenient, inexpensive, and fast. It can be used in biomedical research such as high throughput screening of drugs targeted at DNA-binding proteins.

Simulated single-cycle kinetics improves the design of surface plasmon resonance assays

Instruments based on the surface plasmon resonance (SPR) principle are widely used to monitor in real time molecular interactions between a partner, immobilized on a sensor chip surface and another one injected in a continuous flow of buffer. In a classical SPR experiment, several cycles of binding and regeneration of the surface are performed in order to determine the rate and the equilibrium constants of the reaction. In 2006, Karlsson and co-workers introduced a new method named single-cycle kinetics (SCK) to perform SPR assays. The method consists in injecting sequentially increasing concentrations of the partner in solution, with only one regeneration step performed at the end of the complete binding cycle. A 10 base-pair DNA duplex was characterized kinetically to show how simulated sensorgrams generated by the BiaEvaluation software provided by Biacore™ could really improve the design of SPR assays performed with the SCK method. The DNA duplex was investigated at three temperatures, 10, 20 and 30°C, to analyze fast and slow rate constants. The results show that after a short obligatory preliminary experiment, simulations provide users with the best experimental conditions to be used, in particular, the maximum concentration used to reach saturation, the dilution factor for the serial dilutions of the sample injected and the duration of the dissociation and association phases. The use of simulated single-cycle kinetics saves time and reduces sample consumption. Simulations can also be used to design SPR experiments with ternary complexes.

Structural, Electronic, and Optical Properties of Metallo Base Pairs in Duplex DNA: A Theoretical Insight

Using density functional theory calculations, we investigated the structural, energetic, electronic, and optical properties of recently synthesized duplex DNA containing metal-mediated base pairs. The studied duplex DNA consists of three imidazole (Im) units linked through metal (Im-M-Im, M = metal) and four flanking A:T base pairs (two on each side). We examined the role of artificial base pairing in the presence of two distinctive metal ions, diamagnetic Ag(+) and magnetic Cu(2+) ions, on the stability of duplex DNA. We found that metal-mediated base pairs form stable duplex DNA by direct metal ion coordination to the Im bases. Our results suggest a higher binding stability of base pairing mediated by Cu(2+) ions than by Ag(+) ions, which is attributed to a larger extent of orbital hybridization. We furthermore found that DNA modified with Im-Ag(+)-Im shows the low-energy optical absorption characteristic of π-π*orbital transition of WC A:T base pairs. On the other hand, we found that the low-energy optical absorption peaks for DNA modified with Im-Cu(2+)-Im originate from spin-spin interactions. Additionally, this complex exhibits weak ferromagnetic coupling between Cu(2+) ions and strong spin polarization, which could be used for memory devices. Moreover, analyzing the role of counter ions (Na(+)) and the presence of explicit water molecules on the structural stability and electronic properties of the DNA duplex modified with Im-Ag(+)-Im, we found that the impact of these two factors is negligible. Our results are fruitful for understanding the experimental data and suggest a potential route for constructing effective metal-mediated base pairs in duplex DNA for optoelectronic applications.

Thermodynamic and structural properties of the specific binding between Ag+ ion and C:C mismatched base pair in duplex DNA to form C–Ag–C metal-mediated base pair

Metal ion-nucleic acid interactions have attracted considerable interest for their involvement in structure formation and catalytic activity of nucleic acids. Although interactions between metal ion and mismatched base pair duplex are important to understand mechanism of gene mutations related to heavy metal ions, they have not been well-characterized. We recently found that the Ag(+) ion stabilized a C:C mismatched base pair duplex DNA. A C-Ag-C metal-mediated base pair was supposed to be formed by the binding between the Ag(+) ion and the C:C mismatched base pair to stabilize the duplex. Here, we examined specificity, thermodynamics and structure of possible C-Ag-C metal-mediated base pair. UV melting indicated that only the duplex with the C:C mismatched base pair, and not of the duplexes with the perfectly matched and other mismatched base pairs, was specifically stabilized on adding the Ag(+) ion. Isothermal titration calorimetry demonstrated that the Ag(+) ion specifically bound with the C:C base pair at 1:1 molar ratio with a binding constant of 10(6) M(-1), which was significantly larger than those for nonspecific metal ion-DNA interactions. Electrospray ionization mass spectrometry also supported the specific 1:1 binding between the Ag(+) ion and the C:C base pair. Circular dichroism spectroscopy and NMR revealed that the Ag(+) ion may bind with the N3 positions of the C:C base pair without distorting the higher-order structure of the duplex. We conclude that the specific formation of C-Ag-C base pair with large binding affinity would provide a binding mode of metal ion-DNA interactions, similar to that of the previously reported T-Hg-T base pair. The C-Ag-C base pair may be useful not only for understanding of molecular mechanism of gene mutations related to heavy metal ions but also for wide variety of potential applications of metal-mediated base pairs in various fields, such as material, life and environmental sciences.