Discrimination among individual Watson-Crick base pairs at the termini of single DNA hairpin molecules

Nanoscale α-hemolysin pores can be used to analyze individual DNA or RNA molecules. Serial examination of hundreds to thousands of molecules per minute is possible using ionic current impedance as the measured property. In a recent report, we showed that a nanopore device coupled with machine learning algorithms could automatically discriminate among the four combinations of Watson–Crick base pairs and their orientations at the ends of individual DNA hairpin molecules. Here we use kinetic analysis to demonstrate that ionic current signatures caused by these hairpin molecules depend on the number of hydrogen bonds within the terminal base pair, stacking between the terminal base pair and its nearest neighbor, and 5′ versus 3′ orientation of the terminal bases independent of their nearest neighbors. This report constitutes evidence that single Watson–Crick base pairs can be identified within individual unmodified DNA hairpin molecules based on their dynamic behavior in a nanoscale pore.     

Fast Inactivation in Shaker K+ Channels : Properties of Ionic and Gating Currents

Fast inactivating Shaker H4 potassium channels and nonconducting pore mutant Shaker H4 W434F channels have been used to correlate the installation and recovery of the fast inactivation of ionic current with changes in the kinetics of gating current known as “charge immobilization” (Armstrong, C.M., and F. Bezanilla. 1977. J. Gen. Physiol. 70:567–590.). Shaker H4 W434F gating currents are very similar to those of the conducting clone recorded in potassium-free solutions. This mutant channel allows the recording of the total gating charge return, even when returning from potentials that would largely inactivate conducting channels. As the depolarizing potential increased, the OFF gating currents decay phase at −90 mV return potential changed from a single fast component to at least two components, the slower requiring ∼200 ms for a full charge return. The charge immobilization onset and the ionic current decay have an identical time course. The recoveries of gating current (Shaker H4 W434F) and ionic current (Shaker H4) in 2 mM external potassium have at least two components. Both recoveries are similar at −120 and −90 mV. In contrast, at higher potentials (−70 and −50 mV), the gating charge recovers significantly more slowly than the ionic current. A model with a single inactivated state cannot account for all our data, which strongly support the existence of “parallel” inactivated states. In this model, a fraction of the charge can be recovered upon repolarization while the channel pore is occupied by the NH₂-terminus region.     

Temperature and Electrolyte Optimization of the α-Hemolysin Latch Sensing Zone for Detection of Base Modification in Double-Stranded DNA

The latch region of the wild-type protein pore α-hemolysin (α-HL) constitutes a sensing zone for individual abasic sites (and furan analogs) in double-stranded DNA (dsDNA). The presence of an abasic site or furan within a DNA duplex, electrophoretically captured in the α-HL vestibule and positioned at the latch region, can be detected based on the current blockage prior to duplex unzipping. We investigated variations in blockage current as a function of temperature (12–35°C) and KCl concentration (0.15–1.0 M) to understand the origin of the current signature and to optimize conditions for identifying the base modification. In 1 M KCl solution, substitution of a furan for a cytosine base in the latch region results in an ∼8 kJ mol⁻¹ decrease in the activation energy for ion transport through the protein pore. This corresponds to a readily measured ∼2 pA increase in current at room temperature. Optimal resolution for detecting the presence of a furan in the latch region is achieved at lower KCl concentrations, where the noise in the measured blockage current is significantly lower. The noise associated with the blockage current also depends on the stability of the duplex (as measured from the melting temperature), where a greater noise in the measured blockage current is observed for less stable duplexes.     

Gating pore currents in DIIS4 mutations of NaV1.4 associated with periodic paralysis: saturation of ion flux and implications for disease pathogenesis

S4 voltage–sensor mutations in CaV1.1 and NaV1.4 channels cause the human muscle disorder hypokalemic periodic paralysis (HypoPP). The mechanism whereby these mutations predispose affected sarcolemma to attacks of sustained depolarization and loss of excitability is poorly understood. Recently, three HypoPP mutations in the domain II S4 segment of NaV1.4 were shown to create accessory ionic permeation pathways, presumably extending through the aqueous gating pore in which the S4 segment resides. However, there are several disparities between reported gating pore currents from different investigators, including differences in ionic selectivity and estimates of current amplitude, which in turn have important implications for the pathological relevance of these aberrant currents. To clarify the features of gating pore currents arising from different DIIS4 mutants, we recorded gating pore currents created by HypoPP missense mutations at position R666 in the rat isoform of Nav1.4 (the second arginine from the outside, at R672 in human NaV1.4). Extensive measurements were made for the index mutation, R666G, which created a gating pore that was permeable to K⁺ and Na⁺. This current had a markedly shallow slope conductance at hyperpolarized voltages and robust inward rectification, even when the ionic gradient strongly favored outward ionic flow. These characteristics were accounted for by a barrier model incorporating a voltage-gated permeation pathway with a single cation binding site oriented near the external surface of the electrical field. The amplitude of the R666G gating pore current was similar to the amplitude of a previously described proton-selective current flowing through the gating pore in rNaV1.4-R663H mutant channels. Currents with similar amplitude and cation selectivity were also observed in R666S and R666C mutant channels, while a proton-selective current was observed in R666H mutant channels. These results add support to the notion that HypoPP mutations share a common biophysical profile comprised of a low-amplitude inward current at the resting potential that may contribute to the pathological depolarization during attacks of weakness.     

Rapid discrimination among individual DNA hairpin molecules at single-nucleotide resolution using an ion channel

RNA and DNA strands produce ionic current signatures when driven through an alpha-hemolysin channel by an applied voltage. Here we combine this nanopore detector with a support vector machine (SVM) to analyze DNA hairpin molecules on the millisecond time scale. Measurable properties include duplex stem length, base pair mismatches, and loop length. This nanopore instrument can discriminate between individual DNA hairpins that differ by one base pair or by one nucleotide.     

Important Role of Asparagines in Coupling the Pore and Votage-Sensor Domain in Voltage-Gated Sodium Channels

Voltage-gated sodium (Na_V) channels contain an α-subunit incorporating the channel’s pore and gating machinery composed of four homologous domains (DI–DIV), with a pore domain formed by the S5 and S6 segments and a voltage-sensor domain formed by the S1–S4 segments. During a membrane depolarization movement, the S4s in the voltage-sensor domains exert downstream effects on the S6 segments to control ionic conductance through the pore domain. We used lidocaine, a local anesthetic and antiarrhythmic drug, to probe the role of conserved Asn residues in the S6s of DIII and DIV in Na_V1.5 and Na_V1.4. Previous studies have shown that lidocaine binding to the pore domain causes a decrease in the maximum gating (Q_max) charge of ∼38%, and three-fourths of this decrease results from the complete stabilization of DIII-S4 (contributing a 30% reduction in Q_max) and one-fourth is due to partial stabilization of DIV-S4 (a reduction of 8–10%). Even though substitutions for the Asn in DIV-S6 in Na_V1.5, N1764A and N1764C, produce little ionic current in transfected mammalian cells, they both express robust gating currents. Anthopleurin-A toxin, which inhibits movement of DIV-S4, still reduced Q_max by nearly 30%, a value similar to that observed in wild-type channels, in both N1764A and N1764C. By applying lidocaine and measuring the gating currents, we demonstrated that Asn residues in the S6s of DIII and DIV are important for coupling their pore domains to their voltage-sensor domains, and that Ala and Cys substitutions for Asn in both S6s result in uncoupling of the pore domains from their voltage-sensor domains. Similar observations were made for Na_V1.4, although substitutions for Asn in DIII-S6 showed somewhat less uncoupling.     

Hybridization of single-stranded DNA targets to immobilized complementary DNA probes: comparison of hairpin versus linear capture probes

A microtiter-based assay system is described in which DNA hairpin probes with dangling ends and single-stranded, linear DNA probes were immobilized and compared based on their ability to capture single-strand target DNA. Hairpin probes consisted of a 16 bp duplex stem, linked by a T₂-biotin·dT-T₂ loop. The third base was a biotinylated uracil (U_B) necessary for coupling to avidin coated microtiter wells. The capture region of the hairpin was a 3′ dangling end composed of either 16 or 32 bases. Fundamental parameters of the system, such as probe density and avidin adsorption capacity of the plates were characterized. The target DNA consisted of 65 bases whose 3′ end was complementary to the dangling end of the hairpin or to the linear probe sequence. The assay system was employed to measure the time dependence and thermodynamic stability of target hybridization with hairpin and linear probes. Target molecules were labeled with either a 5′-FITC, or radiolabeled with [γ-³³P]ATP and captured by either linear or hairpin probes affixed to the solid support. Over the range of target concentrations from 10 to 640 pmol hybridization rates increased with increasing target concentration, but varied for the different probes examined. Hairpin probes displayed higher rates of hybridization and larger equilibrium amounts of captured targets than linear probes. At 25 and 45°C, rates of hybridization were better than twice as great for the hairpin compared with the linear capture probes. Hairpin–target complexes were also more thermodynamically stable. Binding free energies were evaluated from the observed equilibrium constants for complex formation. Results showed the order of stability of the probes to be: hairpins with 32 base dangling ends > hairpin probes with l6 base dangling ends > 16 base linear probes > 32 base linear probes. The physical characteristics of hairpins could offer substantial advantages as nucleic acid capture moieties in solid support based hybridization systems.     

Solution structure of the DNA decamer duplex containing a 3'-T x T basepair of the cis-syn cyclobutane pyrimidine dimer: implication for the mutagenic property of the cis-syn dimer

The cis–syn dimer is the major DNA photoproduct produced by UV irradiation. In order to determine the origin of the mutagenic property of the cis–syn dimer, we used NMR restraints and molecular dynamics to determine the solution structure of a DNA decamer duplex containing a wobble pair between the 3′-T of the cis–syn dimer and the opposite T residue (CS/TA duplex). The solution structure of the CS/TA duplex revealed that the 3′-T·T base pair of the cis–syn dimer had base pair geometry that was significantly different from the canonical Watson–Crick base pair and caused destabilization and conformational distortion of its 3′-region. However, a 3′-T·A base pair at the cis–syn dimer within this related DNA decamer maintains the normal Watson–Crick base pair geometry and causes little distortion in the conformation of its 3′-side. Our results show that in spite of its stable hydrogen bonding, the insertion of a T residue opposite the 3′-T of the cis–syn dimer is inhibited by structural distortion caused by the 3′-T·T base pair. This may explain why the frequency of the 3′-T→A transversion, which is the major mutation produced by the cis–syn dimer, is only 4%.     

Ionic current blockades from DNA and RNA molecules in the alpha-hemolysin nanopore

We characterize the substate structure of current blockades produced when single-stranded polynucleotide molecules were electrophoretically driven into the alpha-hemolysin protein pore. We frequently observe substates where the ionic current is reduced by approximately 50%. Most of these substates can be associated with a molecular configuration where a polymer occupies only the vestibule region of the pore, though a few appear related to a polymer occupying only the transmembrane beta-barrel region of the pore. The duration of the vestibule configuration depends on polymer composition and on which end of the polymer, 3' or 5', subsequently threads into the narrowest constriction and initiates translocation. Below approximately 140 mV a polymer is more likely to escape from the vestibule against the applied voltage gradient, while at higher voltages a polymer is more likely to follow the voltage gradient by threading through the narrowest constriction and translocating through the pore. Increasing the applied voltage also increases the duration of the vestibule configuration. A semiquantitative model of these trends suggests that escape has stronger voltage dependence than threading, and that threading is sensitive to polymer orientation while escape is not. These results emphasize the utility of alpha-hemolysin as a model system to study biologically relevant physical and chemical processes at the single-molecule level.     

Immunostimulatory properties of phosphorothioate CpG DNA containing both 3'-5'- and 2'-5'-internucleotide linkages

Synthetic oligodeoxyribonucleotides containing CpG-dinucleotides (CpG DNA) in specific sequence contexts activate the vertebrate immune system. We have examined the effect of 3′-deoxy-2′–5′-ribonucleoside (3′-deoxynucleoside) incorporation into CpG DNA on the immunostimulatory activity. Incorporation of 3′-deoxynucleosides results in the formation of 2′–5′-internucleotide linkages in an otherwise 3′–5′-linked CpG DNA. In studies, both in vitro and in vivo, CpG DNA containing unnatural 3′-deoxynucleoside either within the CpG-dinucleotide or adjacent to the CpG-dinucleotide failed to induce immunostimulatory activity, suggesting that the modification was not recognized by the receptors. Incorporation of the same modification distal to the CpG-dinucleotide in the 5′-flanking sequence potentiated the immunostimulatory activity of the CpG DNA. The same modification when incorporated in the 3′-flanking sequence had an insignificant effect on immunostimulatory activity of CpG DNA. Interestingly, substitution of a 3′-deoxynucleoside in the 5′-flanking sequence distal to the CpG-dinucleotide resulted in increased IL-6 and IL-10 secretion with similar levels of IL-12 compared with parent CpG DNA. The incorporation of the same modification in the 3′-flanking sequence resulted in lower IL-6 and IL-10 secretion with similar levels of IL-12 compared with parent CpG DNA. These results suggest that site-specific incorporation of 3′-deoxynucleotides in CpG DNA modulates immunostimulatory properties.     

A novel single-stranded DNA-specific 3′–5′ exonuclease,Thermus thermophilus exonuclease I,is involved in several DNA repair pathways

Single-stranded DNA (ssDNA)-specific exonucleases (ssExos) are expected to be involved in a variety of DNA repair pathways corresponding to their cleavage polarities; however, the relationship between the cleavage polarity and the respective DNA repair pathways is only partially understood. To understand the cellular function of ssExos in DNA repair better, genes encoding ssExos were disrupted in Thermus thermophilus HB8 that seems to have only a single set of 5′–3′ and 3′–5′ ssExos unlike other model organisms. Disruption of the tthb178 gene, which was expected to encode a 3′–5′ ssExo, resulted in significant increase in the sensitivity to H₂O₂ and frequency of the spontaneous mutation rate, but scarcely affected the sensitivity to ultraviolet (UV) irradiation. In contrast, disruption of the recJ gene, which encodes a 5′–3′ ssExo, showed little effect on the sensitivity to H₂O₂, but caused increased sensitivity to UV irradiation. In vitro characterization revealed that TTHB178 possessed 3′–5′ ssExo activity that degraded ssDNAs containing deaminated and methylated bases, but not those containing oxidized bases or abasic sites. Consequently, we concluded that TTHB178 is a novel 3′–5′ ssExo that functions in various DNA repair systems in cooperation with or independently of RecJ. We named TTHB178 as T. thermophilus exonuclease I.     

Ion permeation and block of the gating pore in the voltage sensor of NaV1.4 channels with hypokalemic periodic paralysis mutations

Hypokalemic periodic paralysis and normokalemic periodic paralysis are caused by mutations of the gating charge–carrying arginine residues in skeletal muscle Na_V1.4 channels, which induce gating pore current through the mutant voltage sensor domains. Inward sodium currents through the gating pore of mutant R666G are only ∼1% of central pore current, but substitution of guanidine for sodium in the extracellular solution increases their size by 13- ± 2-fold. Ethylguanidine is permeant through the R666G gating pore at physiological membrane potentials but blocks the gating pore at hyperpolarized potentials. Guanidine is also highly permeant through the proton-selective gating pore formed by the mutant R666H. Gating pore current conducted by the R666G mutant is blocked by divalent cations such as Ba²⁺ and Zn²⁺ in a voltage-dependent manner. The affinity for voltage-dependent block of gating pore current by Ba²⁺ and Zn²⁺ is increased at more negative holding potentials. The apparent dissociation constant (K_d) values for Zn²⁺ block for test pulses to −160 mV are 650 ± 150 µM, 360 ± 70 µM, and 95.6 ± 11 µM at holding potentials of 0 mV, −80 mV, and −120 mV, respectively. Gating pore current is blocked by trivalent cations, but in a nearly voltage-independent manner, with an apparent K_d for Gd³⁺ of 238 ± 14 µM at −80 mV. To test whether these periodic paralyses might be treated by blocking gating pore current, we screened several aromatic and aliphatic guanidine derivatives and found that 1-(2,4-xylyl)guanidinium can block gating pore current in the millimolar concentration range without affecting normal Na_V1.4 channel function. Together, our results demonstrate unique permeability of guanidine through Na_V1.4 gating pores, define voltage-dependent and voltage-independent block by divalent and trivalent cations, respectively, and provide initial support for the concept that guanidine-based gating pore blockers could be therapeutically useful.     

Biochemical analysis of human Dna2

Yeast Dna2 helicase/nuclease is essential for DNA replication and assists FEN1 nuclease in processing a subset of Okazaki fragments that have long single-stranded 5′ flaps. It is also involved in the maintenance of telomeres. DNA2 is a gene conserved in eukaryotes, and a putative human ortholog of yeast DNA2 (ScDNA2) has been identified. Little is known about the role of human DNA2 (hDNA2), although complementation experiments have shown that it can function in yeast to replace ScDNA2. We have now characterized the biochemical properties of hDna2. Recombinant hDna2 has single-stranded DNA-dependent ATPase and DNA helicase activity. It also has 5′–3′ nuclease activity with preference for single-stranded 5′ flaps adjacent to a duplex DNA region. The nuclease activity is stimulated by RPA and suppressed by steric hindrance at the 5′ end. Moreover, hDna2 shows strong 3′–5′ nuclease activity. This activity cleaves single-stranded DNA in a fork structure and, like the 5′–3′ activity, is suppressed by steric hindrance at the 3′-end, suggesting that the 3′–5′ nuclease requires a 3′ single-stranded end for activation. These biochemical specificities are very similar to those of the ScDna2 protein, but suggest that the 3′–5′ nuclease activity may be more important than previously thought.     

Kinetics of DNA Unwinding by the RecD2 Helicase from Deinococcus radiodurans

RecD2 from Deinococcus radiodurans is a superfamily 1 DNA helicase that is homologous to the Escherichia coli RecD protein but functions outside the context of RecBCD enzyme. We report here on the kinetics of DNA unwinding by RecD2 under single and multiple turnover conditions. There is little unwinding of 20-bp substrates by preformed RecD2-dsDNA complexes when excess ssDNA is present to trap enzyme molecules not bound to the substrate. A shorter 12-bp substrate is unwound rapidly under single turnover conditions. The 12-bp unwinding reaction could be simulated with a mechanism in which the DNA is unwound in two kinetic steps with rate constant of k_unw = 5.5 s⁻¹ and a dissociation step from partially unwound DNA of k_off = 1.9 s⁻¹. These results indicate a kinetic step size of about 3–4 bp, unwinding rate of about 15–20 bp/s, and low processivity (p = 0.74). The reaction time courses with 20-bp substrates, determined under multiple turnover conditions, could be simulated with a four-step mechanism and rate constant values very similar to those for the 12-bp substrate. The results indicate that the faster unwinding of a DNA substrate with a forked end versus only a 5′-terminal single-stranded extension can be accounted for by a difference in the rate of enzyme binding to the DNA substrates. Analysis of reactions done with different RecD2 concentrations indicates that the enzyme forms an inactive dimer or other oligomer at high enzyme concentrations. RecD2 oligomers can be detected by glutaraldehyde cross-linking but not by size exclusion chromatography.     

Conductance and block of hair-cell mechanotransducer channels in transmembrane channel–like protein mutants

Transmembrane channel–like (TMC) proteins TMC1 and TMC2 are crucial to the function of the mechanotransducer (MT) channel of inner ear hair cells, but their precise function has been controversial. To provide more insight, we characterized single MT channels in cochlear hair cells from wild-type mice and mice with mutations in Tmc1, Tmc2, or both. Channels were recorded in whole-cell mode after tip link destruction with BAPTA or after attenuating the MT current with GsMTx-4, a peptide toxin we found to block the channels with high affinity. In both cases, the MT channels in outer hair cells (OHCs) of wild-type mice displayed a tonotopic gradient in conductance, with channels from the cochlear base having a conductance (110 pS) nearly twice that of those at the apex (62 pS). This gradient was absent, with channels at both cochlear locations having similar small conductances, with two different Tmc1 mutations. The conductance of MT channels in inner hair cells was invariant with cochlear location but, as in OHCs, was reduced in either Tmc1 mutant. The gradient of OHC conductance also disappeared in Tmc1/Tmc2 double mutants, in which a mechanically sensitive current could be activated by anomalous negative displacements of the hair bundle. This “reversed stimulus–polarity” current was seen with two different Tmc1/Tmc2 double mutants, and with Tmc1/Tmc2/Tmc3 triple mutants, and had a pharmacological sensitivity comparable to that of native MT currents for most antagonists, except dihydrostreptomycin, for which the affinity was less, and for curare, which exhibited incomplete block. The existence in the Tmc1/Tmc2 double mutants of MT channels with most properties resembling those of wild-type channels indicates that proteins other than TMCs must be part of the channel pore. We suggest that an external vestibule of the MT channel may partly account for the channel’s large unitary conductance, high Ca²⁺ permeability, and pharmacological profile, and that this vestibule is disrupted in Tmc mutants.     

Nonindependent K+ Movement through the Pore in IRK1 Potassium Channels

We measured unidirectional K⁺ in- and efflux through an inward rectifier K channel (IRK1) expressed in Xenopus oocytes. The ratio of these unidirectional fluxes differed significantly from expectations based on independent ion movement. In an extracellular solution with a K⁺ concentration of 25 mM, the data were described by a Ussing flux-ratio exponent, n′, of ∼2.2 and was constant over a voltage range from −50 to −25 mV. This result indicates that the pore of IRK1 channels may be simultaneously occupied by at least three ions. The IRK1 n′ value of 2.2 is significantly smaller than the value of 3.5 obtained for Shaker K channels under identical conditions. To determine if other permeation properties that reflect multi-ion behavior differed between these two channel types, we measured the conductance (at 0 mV) of single IRK1 channels as a function of symmetrical K⁺ concentration. The conductance could be fit by a saturating hyperbola with a half-saturation K⁺ activity of 40 mM, substantially less than the reported value of 300 mM for Shaker K channels. We investigated the ability of simple permeation models based on absolute reaction rate theory to simulate IRK1 current–voltage, conductance, and flux-ratio data. Certain classes of four-barrier, three-site permeation models are inconsistent with the data, but models with high lateral barriers and a deep central well were able to account for the flux-ratio and single channel data. We conclude that while the pore in IRK1 and Shaker channels share important similarities, including K⁺ selectivity and multi-ion occupancy, they differ in other properties, including the sensitivity of pore conductance to K⁺ concentration, and may differ in the number of K⁺ ions that can simultaneously occupy the pore: IRK1 channels may contain three ions, but the pore in Shaker channels can accommodate four or more ions.     

New insights into Hoogsteen base pairs in DNA duplexes from a structure-based survey

Hoogsteen (HG) base pairs (bps) provide an alternative pairing geometry to Watson–Crick (WC) bps and can play unique functional roles in duplex DNA. Here, we use structural features unique to HG bps (syn purine base, HG hydrogen bonds and constricted C1′–C1′ distance across the bp) to search for HG bps in X-ray structures of DNA duplexes in the Protein Data Bank. The survey identifies 106 A•T and 34 G•C HG bps in DNA duplexes, many of which are undocumented in the literature. It also uncovers HG-like bps with syn purines lacking HG hydrogen bonds or constricted C1′–C1′ distances that are analogous to conformations that have been proposed to populate the WC-to-HG transition pathway. The survey reveals HG preferences similar to those observed for transient HG bps in solution by nuclear magnetic resonance, including stronger preferences for A•T versus G•C bps, TA versus GG steps, and also suggests enrichment at terminal ends with a preference for 5′-purine. HG bps induce small local perturbations in neighboring bps and, surprisingly, a small but significant degree of DNA bending (∼14°) directed toward the major groove. The survey provides insights into the preferences and structural consequences of HG bps in duplex DNA.     

Electron attachment-induced DNA single-strand breaks at the pyrimidine sites

To elucidate the contribution of pyrimidine in DNA strand breaks caused by low-energy electrons (LEEs), theoretical investigations of the LEE attachment-induced C_3′–O_3′, and C_5′–O_5′ σ bond as well as N-glycosidic bond breaking of 2′-deoxycytidine-3′,5′-diphosphate and 2′-deoxythymidine-3′,5′-diphosphate were performed using the B3LYP/DZP++ approach. The base-centered radical anions are electronically stable enough to assure that either the C–O or glycosidic bond breaking processes might compete with the electron detachment and yield corresponding radical fragments and anions. In the gas phase, the computed glycosidic bond breaking activation energy (24.1 kcal/mol) excludes the base release pathway. The low-energy barrier for the C_3′–O_3′ σ bond cleavage process (∼6.0 kcal/mol for both cytidine and thymidine) suggests that this reaction pathway is the most favorable one as compared to other possible pathways. On the other hand, the relatively low activation energy barrier (∼14 kcal/mol) for the C_5′–O_5′ σ bond cleavage process indicates that this bond breaking pathway could be possible, especially when the incident electrons have relatively high energy (a few electronvolts). The presence of the polarizable medium greatly increases the activation energies of either C–O σ bond cleavage processes or the N-glycosidic bond breaking process. The only possible pathway that dominates the LEE-induced DNA single strands in the presence of the polarizable surroundings (such as in an aqueous solution) is the C_3′–O_3′ σ bond cleavage (the relatively low activation energy barrier, ∼13.4 kcal/mol, has been predicted through a polarizable continuum model investigation). The qualitative agreement between the ratio for the bond breaks of C_5′–O_5′, C_3′–O_3′ and N-glycosidic bonds observed in the experiment of oligonucleotide tetramer CGAT and the theoretical sequence of the bond breaking reaction pathways have been found. This consistency between the theoretical predictions and the experimental observations provides strong supportive evidences for the base-centered radical anion mechanism of the LEE-induced single-strand bond breaking around the pyrimidine sites of the DNA single strands.     

Structural basis of replication origin recognition by the DnaA protein

Escherichia coli DnaA binds to 9 bp sequences (DnaA boxes) in the replication origin, oriC, to form a complex initiating chromosomal DNA replication. In the present study, we determined the crystal structure of its DNA-binding domain (domain IV) complexed with a DnaA box at 2.1 Å resolution. DnaA domain IV contains a helix–turn–helix motif for DNA binding. One helix and a loop of the helix– turn–helix motif are inserted into the major groove and 5 bp (3′ two-thirds of the DnaA box sequence) are recognized through base-specific hydrogen bonds and van der Waals contacts with the C5-methyl groups of thymines. In the minor groove, Arg399, located in the loop adjacent to the motif, recognizes three more base pairs (5′ one-third of the DnaA box sequence) by base-specific hydrogen bonds. DNA bending by ~28° was also observed in the complex. These base-specific interactions explain how DnaA exhibits higher affinity for the strong DnaA boxes (R1, R2 and R4) than the weak DnaA boxes (R3 and M) in the replication origin.