Molecular structure of nicked DNA: a substrate for DNA repair enzymes

Cracked actin filaments as mechanosensitive receptors

Actin filament networks are exposed to mechanical stimuli, but the effect of strain on actin filament structure has not been well established in molecular detail. This is a critical gap in understanding because the activity of a variety of actin-binding proteins has recently been determined to be altered by actin filament strain. We therefore used all-atom molecular dynamics simulations to apply tensile strains to actin filaments and find that changes in actin subunit organization are minimal in mechanically strained, but intact, actin filaments. However, a conformational change disrupts the critical D-loop to W-loop connection between longitudinal neighboring subunits, which leads to a metastable cracked conformation of the actin filament whereby one protofilament is broken prior to filament severing. We propose that the metastable crack presents a force-activated binding site for actin regulatory factors that specifically associate with strained actin filaments. Through protein-protein docking simulations, we find that 43 evolutionarily diverse members of the dual zinc-finger-containing LIM-domain family, which localize to mechanically strained actin filaments, recognize two binding sites exposed at the cracked interface. Furthermore, through its interactions with the crack, LIM domains increase the length of time damaged filaments remain stable. Our findings propose a new molecular model for mechanosensitive binding to actin filaments.

In Addition to Damaging the Plasma Membrane, Phenolic Monoterpenoid Carvacrol Can Bind to the Minor Groove of DNA of Phytopathogenic Fungi to Potentially Control Tea Leaf Spot Caused by Lasiodiplodia theobromae

Carvacrol expresses a wide range of biological activities, but the studies of its mechanisms focused on bacteria, mainly involving the destruction of the plasma membrane. In this study, carvacrol exhibited strong activities against several phytopathogenic fungi and demonstrated a novel antifungal mechanism against Lasiodiplodia theobromae. RNA sequencing indicated that many genes of L. theobromae hyphae were predominately induced by carvacrol, particularly those involved in replication and transcription. Hyperchromic, hypsochromic, and bathochromic effects in the UV-visible absorption spectrum were observed following titration of calf thymus DNA (ctDNA) and carvacrol, which indicated the formation of a DNA-carvacrol complex. Circular dichroism (CD) spectroscopy indicated that the response of DNA to carvacrol was similar to that of 4',6-diamidino-2-phenylindole (DAPI) but different from that of ethidium bromide (EB), implying the ionic bonds between carvacrol and ctDNA. Fluorescence spectrum (FS) analysis indicated that carvacrol quenched the fluorescence of double-stranded DNA (dsDNA) more than single-stranded DNA, indicating that carvacrol mainly bound to dsDNA. A displacement assay showed that carvacrol reduced the fluorescence intensity of the DNA-DAPI complex through competition with DAPI, but this did not occur for DNA-EB. The FS assay revealed that carvacrol bound to the AAA sequence on the minor groove of ds-oligonucleotides. The hydroxyl of carvacrol was verified to bind to ctDNA through a comparative test in which structural analogs of carvacrol, including thymol and 4-ethyl-1,2-dimethyl, were analyzed. The current study indicated carvacrol can destruct plasma membranes and bind to the minor groove of DNA, inhibiting fungal proliferation by disturbing the stability of dsDNA.

Scattering of Attosecond Laser Pulses on a DNA Molecule during Its Nicking and Bending

It is well known that X-ray crystallography is based on X-ray diffraction (XRD) for atoms and molecules. The diffraction pattern arises as a result of scattering of incident radiation, which makes it possible to determine the structure of the scattering substance. With the advent of ultrashort radiation sources, the theory and interpretation of X-ray diffraction analysis have remained the same. This work shows that when an attosecond laser pulse is scattered on a DNA molecule, including during its nicking and bending, the pulse duration is an important characteristic of the scattering. In this case, the diffraction pattern changes significantly compared to the previously known scattering theory. The results obtained must be used in XRD theory to study DNA structures, their mutations and damage, since the previously known theory can produce large errors and, therefore, the DNA structure can be "decoding" incorrectly.

Cracked actin filaments as mechanosensitive receptors

Actin filament networks are exposed to mechanical stimuli, but the effect of strain on actin filament structure has not been well-established in molecular detail. This is a critical gap in understanding because the activity of a variety of actin-binding proteins have recently been determined to be altered by actin filament strain. We therefore used all-atom molecular dynamics simulations to apply tensile strains to actin filaments and find that changes in actin subunit organization are minimal in mechanically strained, but intact, actin filaments. However, a conformational change disrupts the critical D-loop to W-loop connection between longitudinal neighboring subunits, which leads to a metastable cracked conformation of the actin filament, whereby one protofilament is broken prior to filament severing. We propose that the metastable crack presents a force-activated binding site for actin regulatory factors that specifically associate with strained actin filaments. Through protein-protein docking simulations, we find that 43 evolutionarily-diverse members of the dual zinc finger containing LIM domain family, which localize to mechanically strained actin filaments, recognize two binding sites exposed at the cracked interface. Furthermore, through its interactions with the crack, LIM domains increase the length of time damaged filaments remain stable. Our findings propose a new molecular model for mechanosensitive binding to actin filaments.

Toward Force Fields with Improved Base Stacking Descriptions

Recent DNA force fields indicate good performance in describing flexibility and structural stability of double-stranded B-DNA. However, it is not clear how accurately base stacking interactions are represented that are critical for simulating structure formation processes and conformational changes. Based on the equilibrium nucleoside association and base pair nicking, we find that the recent Tumuc1 force field improves the description of base stacking compared to previous state-of-the-art force fields. Nevertheless, base pair stacking is still overstabilized compared to experiment. We propose a rapid method to reweight calculated free energies of stacking upon force field modifications in order to generate improved parameters. A decrease of the Lennard-Jones attraction between nucleo-bases alone appears insufficient; however, adjustments in the partial charge distribution on base atoms could help to further improve the force field description of base stacking.

Crystal structure of an RNA/DNA strand exchange junction

Short segments of RNA displace one strand of a DNA duplex during diverse processes including transcription and CRISPR-mediated immunity and genome editing. These strand exchange events involve the intersection of two geometrically distinct helix types-an RNA:DNA hybrid (A-form) and a DNA:DNA homoduplex (B-form). Although previous evidence suggests that these two helices can stack on each other, it is unknown what local geometric adjustments could enable A-on-B stacking. Here we report the X-ray crystal structure of an RNA-5'/DNA-3' strand exchange junction at an anisotropic resolution of 1.6 to 2.2 Å. The structure reveals that the A-to-B helical transition involves a combination of helical axis misalignment, helical axis tilting and compression of the DNA strand within the RNA:DNA helix, where nucleotides exhibit a mixture of A- and B-form geometry. These structural principles explain previous observations of conformational stability in RNA/DNA exchange junctions, enabling a nucleic acid architecture that is repeatedly populated during biological strand exchange events.

Supercoiling and looping promote DNA base accessibility and coordination among distant sites

DNA in cells is supercoiled and constrained into loops and this supercoiling and looping influence every aspect of DNA activity. We show here that negative supercoiling transmits mechanical stress along the DNA backbone to disrupt base pairing at specific distant sites. Cooperativity among distant sites localizes certain sequences to superhelical apices. Base pair disruption allows sharp bending at superhelical apices, which facilitates DNA writhing to relieve torsional strain. The coupling of these processes may help prevent extensive denaturation associated with genomic instability. Our results provide a model for how DNA can form short loops, which are required for many essential processes, and how cells may use DNA loops to position nicks to facilitate repair. Furthermore, our results reveal a complex interplay between site-specific disruptions to base pairing and the 3-D conformation of DNA, which influences how genomes are stored, replicated, transcribed, repaired, and many other aspects of DNA activity.

DNA punch cards for storing data on native DNA sequences via enzymatic nicking

Synthetic DNA-based data storage systems have received significant attention due to the promise of ultrahigh storage density and long-term stability. However, all known platforms suffer from high cost, read-write latency and error-rates that render them noncompetitive with modern storage devices. One means to avoid the above problems is using readily available native DNA. As the sequence content of native DNA is fixed, one can modify the topology instead to encode information. Here, we introduce DNA punch cards, a macromolecular storage mechanism in which data is written in the form of nicks at predetermined positions on the backbone of native double-stranded DNA. The platform accommodates parallel nicking on orthogonal DNA fragments and enzymatic toehold creation that enables single-bit random-access and in-memory computations. We use Pyrococcus furiosus Argonaute to punch files into the PCR products of Escherichia coli genomic DNA and accurately reconstruct the encoded data through high-throughput sequencing and read alignment.

Current synthetic DNA-based data storage systems have high recording costs, read-write latency and error-rates that make them uncompetitive compared to traditional digital storage. The authors use nicks in native DNA to encode data in parallel and create access sites for in-memory computations.

A Strategy Employing a TF-Splinting Duplex Nanoswitch to Achieve Single-Step, Enzyme-Free, Signal-On Detection of l-Tryptophan

Transcription factor (TF)-based metabolite detection mainly depends on TF-regulated gene expression in cells. From TF activation to gene transcription and translation, the signal travels a relatively long way before it is received. Here, we propose a TF-splinting duplex DNA nanoswitch to detect metabolites. We show its feasibility using tryptophan repressor (TrpR) to detect l-tryptophan as a model. The assay has been optimized and characterized after obtaining a proof of concept, and the detection of l-tryptophan in complex biological samples is feasible. Unlike an equivalent gene expression approach, the whole process is a single-step, enzyme-free, and signal-on method. It can be completed within 20 min. This proposed TF-splinting duplex has the potential to be applied to the quick and convenient detection of other metabolites or even TFs.

DNA clamp function of the monoubiquitinated Fanconi anaemia ID complex

The FANCI-FANCD2 (ID) complex, mutated in the Fanconi Anemia (FA) cancer predisposition syndrome, is required for the repair of interstrand crosslinks (ICL) and related lesions¹. The FA pathway is activated when a replication fork stalls at an ICL², triggering the mono-ubiquitination of the ID complex. ID mono-ubiquitination is essential for ICL repair by excision, translesion synthesis and homologous recombination, but its function was hitherto unknown^1,3. Here, the 3.5 Å cryo-EM structure of mono-ubiquitinated ID (ID^Ub) bound to DNA reveals that it forms a closed ring that encircles the DNA. Compared to the cryo-EM structure of the non-ubiquitinated ID complex bound to ICL DNA, described here as well, mono-ubiquitination triggers a complete re-arrangement of the open, trough-like ID structure through the ubiquitin of one protomer binding to the other protomer in a reciprocal fashion. The structures, in conjunction with biochemical data, indicate the mono-ubiquitinated ID complex looses its preference for ICL and related branched DNA structures, becoming a sliding DNA clamp that can coordinate the subsequent repair reactions. Our findings also reveal how mono-ubiquitination in general can induce an alternate structure with a new function.

Critical Sites of DNA Backbone Integrity for Damaged Base Removal by Formamidopyrimidine-DNA Glycosylase

DNA glycosylases, the enzymes that initiate base excision DNA repair, recognize damaged bases through a series of precisely orchestrated movements. Most glycosylases sharply kink the DNA axis at the lesion site and extrude the target base from the DNA double helix into the enzyme's active site. Little attention has been paid so far to the role of the physical continuity of the DNA backbone in allowing the required conformational distortion. Here, we analyze base excision by formamidopyrimidine-DNA glycosylase (Fpg) from substrates keeping all phosphates but containing a nick within three nucleotides of the lesion in either DNA strand. Four phosphoester linkages at the damaged nucleotide and two nucleotides 3' to it were essential for Fpg activity, while the breakage of the others, even at the same critical phosphates, had no effect or even stimulated the reaction. Reduction of the likelihood of hydrogen bonding at the nicks by using dideoxynucleotides as their 3'-terminal groups was more detrimental for the activity. All phosphoester bonds in the complementary strand were dispensable for base excision, but nicks close to the orphaned nucleotide caused early termination of damaged strand cleavage. Elastic network analysis of Fpg-DNA structures showed that the vibrational motions of the critical phosphates are strongly correlated, in part due to the presence of the protein. Overall, our results suggest that mechanical forces propagating along the DNA backbone play a critical role in the correct conformational distortion of DNA by Fpg and possibly by other target base-everting DNA glycosylases.

A DNA structural alphabet provides new insight into DNA flexibility

DNA is a structurally plastic molecule, and its biological function is enabled by adaptation to its binding partners. To identify the DNA structural polymorphisms that are possible in such adaptations, the dinucleotide structures of 60 000 DNA steps from sequentially nonredundant crystal structures were classified and an automated protocol assigning 44 distinct structural (conformational) classes called NtC (for Nucleotide Conformers) was developed. To further facilitate understanding of the DNA structure, the NtC were assembled into the DNA structural alphabet CANA (Conformational Alphabet of Nucleic Acids) and the projection of CANA onto the graphical representation of the molecular structure was proposed. The NtC classification was used to define a validation score called confal, which quantifies the conformity between an analyzed structure and the geometries of NtC. NtC and CANA assignment were applied to analyze the structural properties of typical DNA structures such as Dickerson-Drew dodecamers, guanine quadruplexes and structural models based on fibre diffraction. NtC, CANA and confal assignment, which is accessible at the website https://dnatco.org, allows the quantitative assessment and validation of DNA structures and their subsequent analysis by means of pseudo-sequence alignment. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:Acta_Cryst_D:2.

Orientation-dependent FRET system reveals differences in structures and flexibilities of nicked and gapped DNA duplexes

Differences in structures and flexibilities of DNA duplexes play important roles on recognition by DNA-binding proteins. We herein describe a novel method for structural analyses of DNA duplexes by using orientation dependence of Förster resonance energy transfer (FRET). We first analyzed canonical B-form duplex and correct structural parameters were obtained. The experimental FRET efficiencies were in excellent agreement with values theoretically calculated by using determined parameters. We then investigated DNA duplexes with nick and gaps, which are key intermediates in DNA repair systems. Effects of gap size on structures and flexibilities were successfully revealed. Since our method is facile and sensitive, it could be widely used to analyze DNA structures containing damages and non-natural molecules.

Biomimetic Autonomous Enzymatic Nanowalker of High Fuel Efficiency

Replicating efficient chemical energy utilization of biological nanomotors is one ultimate goal of nanotechnology and energy technology. Here, we report a rationally designed autonomous bipedal nanowalker made of DNA that achieves a fuel efficiency of less than two fuel molecules decomposed per productive forward step, hence breaking a general threshold for chemically powered machines invented to date. As a genuine enzymatic nanomotor without changing itself nor the track, the walker demonstrates a sustained motion on an extended double-stranded track at a speed comparable to previous burn-bridge motors. Like its biological counterparts, this artificial nanowalker realizes multiple chemomechanical gatings, especially a bias-generating product control unique to chemically powered nanomotors. This study yields rich insights into how pure physical effects facilitate harvest of chemical energy at the single-molecule level and provides a rarely available motor system for future development toward replicating the efficient, repeatable, automatic, and mechanistically sophisticated transportation seen in biomotor-based intracellular transport but beyond the capacity of the current burn-bridge motors.

Key determinants of target DNA recognition by retroviral intasomes

Background

Retroviral integration favors weakly conserved palindrome sequences at the sites of viral DNA joining and generates a short (4–6 bp) duplication of host DNA flanking the provirus. We previously determined two key parameters that underlie the target DNA preference for prototype foamy virus (PFV) and human immunodeficiency virus type 1 (HIV-1) integration: flexible pyrimidine (Y)/purine (R) dinucleotide steps at the centers of the integration sites, and base contacts with specific integrase residues, such as Ala188 in PFV integrase and Ser119 in HIV-1 integrase. Here we examined the dinucleotide preference profiles of a range of retroviruses and correlated these findings with respect to length of target site duplication (TSD).

Results

Integration datasets covering six viral genera and the three lengths of TSD were accessed from the literature or generated in this work. All viruses exhibited significant enrichments of flexible YR and/or selection against rigid RY dinucleotide steps at the centers of integration sites, and the magnitude of this enrichment inversely correlated with TSD length. The DNA sequence environments of in vivo-generated HIV-1 and PFV sites were consistent with integration into nucleosomes, however, the local sequence preferences were largely independent of target DNA chromatinization. Integration sites derived from cells infected with the gammaretrovirus reticuloendotheliosis virus strain A (Rev-A), which yields a 5 bp TSD, revealed the targeting of global chromatin features most similar to those of Moloney murine leukemia virus, which yields a 4 bp duplication. In vitro assays revealed that Rev-A integrase interacts with and is catalytically stimulated by cellular bromodomain containing 4 protein.

Conclusions

Retroviral integrases have likely evolved to bend target DNA to fit scissile phosphodiester bonds into two active sites for integration, and viruses that cut target DNA with a 6 bp stagger may not need to bend DNA as sharply as viruses that cleave with 4 bp or 5 bp staggers. For PFV and HIV-1, the selection of signature bases and central flexibility at sites of integration is largely independent of chromatin structure. Furthermore, global Rev-A integration is likely directed to chromatin features by bromodomain and extraterminal domain proteins.

Electronic supplementary material

The online version of this article (doi:10.1186/s12977-015-0167-3) contains supplementary material, which is available to authorized users.

Syntheses and photophysical properties of 5'-6-locked fluorescent nucleosides

Nine fluorescent 5'-6-locked nucleosides were synthesized by condensation of various 1,2-diketones with 5-amino-2'-deoxycytidine. The nucleosides have different substituents on the pyrazine core structure, ranging from two methyl groups to polyaromatic rings. The photophysical properties of each nucleoside were determined, with the nucleosides displaying diverse absorption and emission maxima, extinction coefficients and quantum yields. The nucleoside with the highest fluorescence brightness was phosphitylated and incorporated into an oligonucleotide by means of automated oligonucleotide synthesis. The labelled oligonucleotide in aqueous buffer exhibited a substantially lowered extinction coefficient and quantum yield compared to the nucleoside in THF. The photophysical properties of the nucleoside were also compared in different DNA structural contexts, a single strand, a 14-mer duplex, a 14-mer duplex with an 11-mer overhang, and a 25-mer nicked duplex labelled at the nick site. Circular dichroism and melting temperature studies verified that the nucleoside did not perturb or destabilize the DNA helixes. In fact, when incorporated at the nick site, the nucleoside was found to stabilize the nicked duplex notably compared to its unmodified counterpart. The brightness of the fluorescent nucleoside in DNA increased as the polarity of its surroundings decreased, being highest in the 25-mer nicked duplex where exposure to the polar solvent is minimized by stacking to the adjacent bases on both the 3'- and 5'-side. The nucleosides brightness in the nicked duplex was also found to increase with lowered temperature, in accordance with expected temperature-dependent changes in the stacked-unstacked equilibrium at the nick site.

Conformations of the sugar-phosphate backbone in helical DNA crystal structures

The DNA backbone geometry was analyzed for 96 crystal structures of oligodeoxynucleotides. The ranges and mean values of the torsion angles for the observed subclasses of the A-, B-, and Z-DNA conformational types were determined by analyzing distributions of the torsion angles and scattergrams relating pairs of angles.

First principles effective electronic couplings for hole transfer in natural and size-expanded DNA

Hole transfer processes between base pairs in natural DNA and size-expanded DNA (xDNA) are studied and compared, by means of an accurate first principles evaluation of the effective electronic couplings (also known as transfer integrals), in order to assess the effect of the base augmentation on the efficiency of charge transport through double-stranded DNA. According to our results, the size expansion increases the average electronic coupling, and thus the CT rate, with potential implications in molecular biology and in the implementation of molecular nanoelectronics. Our analysis shows that the effect of the nucleobase expansion on the charge-transfer (CT) rate is sensitive to the sequence of base pairs. Furthermore, we find that conformational variability is an important factor for the modulation of the CT rate. From a theoretical point of view, this work offers a contribution to the CT chemistry in pi-stacked arrays. Indeed, we compare our methodology against other standard computational frameworks that have been adopted to tackle the problem of CT in DNA, and unravel basic principles that should be accounted for in selecting an appropriate theoretical level.

Synthesis and properties of nicked dumbbell and dumbbell DNA conjugates having stilbenedicarboxamide linkers

The synthesis and properties of nicked dumbbell and dumbbell DNA conjugates having A-tract base pair domains connected by rod-like stilbenedicarboxamide linkers are reported. The nicked dumbbells have one to eight dA-dT base pairs and are missing a sugar-phosphate bond either between the linker and a thymine nucleoside residue or between two thymine residues. Chemical ligation of all of the nicked dumbbells with cyanogen bromide affords the dumbbell conjugates in good yield, providing the smallest mini-dumbbells prepared to date. The dumbbells have exceptionally high thermal stability, whereas the nicked dumbbells are only marginally more stable than the hairpin structures on either side of the nick. The structures of the nicked dumbbells and dumbbells have been investigated using a combination of circular dichroism spectroscopy and molecular modeling. The base pair domains are found to adopt normal B'-DNA geometry and thus provide a helical ruler for studies of the distance and angular dependence of electronic interactions between the chromophore linkers.

Base-stacking and base-pairing contributions into thermal stability of the DNA double helix

Two factors are mainly responsible for the stability of the DNA double helix: base pairing between complementary strands and stacking between adjacent bases. By studying DNA molecules with solitary nicks and gaps we measure temperature and salt dependence of the stacking free energy of the DNA double helix. For the first time, DNA stacking parameters are obtained directly (without extrapolation) for temperatures from below room temperature to close to melting temperature. We also obtain DNA stacking parameters for different salt concentrations ranging from 15 to 100 mM Na⁺. From stacking parameters of individual contacts, we calculate base-stacking contribution to the stability of A•T- and G•C-containing DNA polymers. We find that temperature and salt dependences of the stacking term fully determine the temperature and the salt dependence of DNA stability parameters. For all temperatures and salt concentrations employed in present study, base-stacking is the main stabilizing factor in the DNA double helix. A•T pairing is always destabilizing and G•C pairing contributes almost no stabilization. Base-stacking interaction dominates not only in the duplex overall stability but also significantly contributes into the dependence of the duplex stability on its sequence.

Stacked-unstacked equilibrium at the nick site of DNA

Stability of duplex DNA with respect to separation of complementary strands is crucial for DNA executing its major functions in the cell and it also plays a central role in major biotechnology applications of DNA: DNA sequencing, polymerase chain reaction, and DNA microarrays. Two types of interaction are well known to contribute to DNA stability: stacking between adjacent base-pairs and pairing between complementary bases. However, their contribution into the duplex stability is yet to be determined. Now we fill this fundamental gap in our knowledge of the DNA double helix. We have prepared a series of 32, 300 bp-long DNA fragments with solitary nicks in the same position differing only in base-pairs flanking the nick. Electrophoretic mobility of these fragments in the gel has been studied. Assuming the equilibrium between stacked and unstacked conformations at the nick site, all 32 stacking free energy parameters have been obtained. Only ten of them are essential and they govern the stacking interactions between adjacent base-pairs in intact DNA double helix. A full set of DNA stacking parameters has been determined for the first time. From these data and from a well-known dependence of DNA melting temperature on G.C content, the contribution of base-pairing into duplex stability has been estimated. The obtained energy parameters of the DNA double helix are of paramount importance for understanding sequence-dependent DNA flexibility and for numerous biotechnology applications.

DNA-Mediated Exciton Coupling and Electron Transfer between Donor and Acceptor Stilbenes Separated by a Variable Number of Base Pairs

The synthesis, steady-state spectroscopy, and transient absorption spectroscopy of DNA conjugates possessing both stilbene electron donor and electron acceptor chromophores are described. These conjugates are proposed to form nicked DNA dumbbell structures in which a stilbenedicarboxamide acceptor and stilbenediether donor are separated by variable numbers of A-T or G-C base pairs. The nick is located either adjacent to one of the chromophores or between two of the bases. Thermal dissociation profiles indicate that stable structures are formed possessing as few as two A-T base pairs. Circular dichroism (CD) spectra in the base pair region are characteristic of B-DNA duplex structures, whereas CD spectra at longer wavelengths display two bands attributed to exciton coupling between the two stilbenes. The sign and intensity of these bands are dependent upon both the distance between the chromophores and the dihedral angle between their transition dipoles [Deltaepsilon approximately Rda(-2) sin(2theta)]. Pulsed laser excitation of the stilbenediamide results in creation of the acceptor-donor radical ion pair, which decays via charge recombination. The dynamics of charge separation and charge recombination display an exponential distance dependence, similar to that observed previously for systems in which guanine serves as the electron donor. Unlike exciton coupling between the stilbenes, there is no apparent dependence of the charge-transfer rates upon the dihedral angle between donor and acceptor stilbenes. The introduction of a single G-C base pair between the donor and acceptor results in a change in the mechanism for charge separation from single step superexchange to hole hopping.

High-throughput approach for detection of DNA bending and flexibility based on cyclization

We have developed a high-throughput approach to the labor-intensive problems of DNA cyclization, which we use to characterize DNA curvature and mechanical properties. The method includes a combinatorial approach to make the DNA constructs needed and automated real-time measurement of the kinetics using fluorescence. We validated the approach and investigated the flexibility of two kinds of nicked DNA and AT dinucleotide repeats. We found that, although the nicks hardly alter the bending flexibility, they significantly increase the torsional flexibility, and that the AT repeat has 28% (+/-12%) lower bending rigidity than a generic DNA sequence.

Monitoring of single nicks in duplex DNA by gel electrophoretic mobility-shift assay

We demonstrate that the gel electrophoretic mobility-shift assay (EMSA) can be used for site-selective and quantitative monitoring of nicks in linear double-stranded DNA (dsDNA) thus allowing to expediently follow the nicking activity of enzymes or other agents targeted to a designated dsDNA site. At elevated temperature and/or in the presence of urea, DNA fragments carrying a single nick produced by the nicking enzyme N.BstNBI exhibit a well-detectable gel retardation effect. On the basis of permutation analysis, the decreased electrophoretic mobility of nicked dsDNA fragments is attributed to a bend (or hinge) in the DNA double helix sequence-specifically generated by a nick. Since nick-induced DNA bending depends on interaction between base pairs adjacent to a nick, the change in mobility is different for nicked DNA sites with different sequences. Therefore, EMSA monitoring of differential mobility change caused by nicks within various DNA sequences could be useful for studying the differential base stacking and nearest-neighbor energetics.

Inhibitory module of Ets-1 allosterically regulates DNA binding through a dipole-facilitated phosphate contact

DNA binding of the transcription factor Ets-1 is negatively regulated by three inhibitory helices that lie near the ETS domain. The current model suggests that this negative regulation, termed autoinhibition, is caused by the energetic expense of a DNA-induced structural transition that includes the unfolding of one inhibitory helix. This report investigates the role of helix H1 of the ETS domain in the autoinhibition mechanism. Previous structural studies modeled the inhibitory helices packing together and connecting with helix H1, suggesting a role of this helix in the configuration of an inhibitory module. Recently, high-resolution structures of the ETS domain-DNA interface indicate that the N terminus of helix H1 directly contacts DNA. The contact, which is augmented by the macrodipole of helix H1, consists of a hydrogen bond between the amide NH of leucine 337 in helix H1 and the oxygen of a corresponding phosphate. We propose that this hydrogen bond positions helix H1 to be a link between autoinhibition and DNA binding. Four independent approaches tested this hypothesis. First, the hydrogen bond was disrupted by removal of the phosphate in a missing phosphate analysis. Second, base pairs that surround the helix H1-contacting phosphate and appear to dictate DNA backbone conformation were mutated. Next, a hydrophobic residue in helix H1 that is expected to position the N terminus of the helix was altered. Finally, a residue on the surface of helix H1 that may contact the inhibitory elements was changed. In each case DNA binding and autoinhibition was affected. Taken together, the results demonstrate the role of the dipole-facilitated phosphate contact in DNA binding. Furthermore, the findings support a model in which helix H1 links the inhibitory elements to the ETS domain. We speculate that this helix, which is conserved in all Ets proteins, provides a common route to regulation.

Thermodynamic analysis of stacking hybridization of oligonucleotides with DNA template

Contiguous stacking hybridization of oligodeoxyribonucleotides with DNA as template was investigated using three types of complexes: oligonucleotide contiguously stacked with the stem of the preformed minihairpin (complexes I), oligonucleotide tandems containing two (complexes II) or three (complexes III) short oligomers with a common DNA template. Enthalpy Delta H degrees and entropy Delta S degrees of the coaxial stacking of adjacent duplexes were determined for GC/G*pC, GT/A*pC, AC/G*pT, AT/A*pT, CT/A*pG, AG/C*pT, AA/T*pT and TT/A*pA nicked (*) dinucleotide base pairs. The maximal efficiency of co-operative interaction was found for the GC/G*pC interface (Delta G degrees(NN/N*pN)=-2.7 kcal/mol) and the minimal one for the AA/T*pT interface (Delta G degrees(NN/N*pN)=-1.2 kcal/mol) at 37 degrees C. As a whole, the efficiency of the base pairs interaction Delta G degrees(NN/N*pN) in the nick is not lower than that within the intact DNA helix (Delta G degrees(NN/NN)). These observed Delta G degrees(NN/N*pN) values are proposed may include the effect of the partial removal of fraying at the adjacent helix ends additionally to the effect of the direct stacking of the terminal base pairs in the duplex junction (Delta G degrees(NN/NN). The thermodynamic parameters have been found to describe adequately the formation of all tandem complexes of the II and III types with oligonucleotides of various length and hybridization properties. The performed thermodynamic analysis reveals features of stacking oligonucleotide hybridization which allow one to predict the temperature dependence of association of oligonucleotides and the DNA template within tandem complexes as well as to determine optimal concentration for formation of these complexes characterized by high co-operativity level.

Intrinsic bending and deformability at the T-A step of CCTTTAAAGG: a comparative analysis of T-A and A-T steps within A-tracts

Introduction of a T-A or pyrimidine-purine step into a straight and rigid A-tract can cause a positive roll deformation that kinks the DNA helix at that step. In CCTTTAAAGG, the central T-A step has an 8.6 degrees bend toward the major groove. We report the structural analysis of CCTTTAAAGG and a comparison with 25 other representative crystal structures from the NDB containing at least four consecutive A or T bases. On average, more local bending occurs at the disruptive T-A step (8.21 degrees ) than at an A-T step (5.71 degrees ). In addition, A-tracts containing an A-T step are more bent than are pure A-tracts, and hence A-A and A-T steps are not equivalent. All T-A steps examined exhibit positive roll, bending towards the major groove, while A-T steps display negative roll and bend slightly towards the minor groove. This illustrates how inherent negative and positive roll are, respectively, at A-T and T-A steps within A-tracts. T-A steps are more deformable, showing larger and more variable deformations of minor groove width, rise, cup, twist, and buckle. Standard deviations of twist, rise, and cup for T-A steps are 6.66 degrees, 0.55 A, and 15.90 degrees, versus 2.28 degrees, 0.21 A, and 2.99 degrees for A-T steps. Packing constraints determine which local values of these helical parameters an individual T-A step will adopt. For instance, with CCTTTAAAGG and three isomorphous structures, CGATTAATCG, CGATATATCG, and CGATCGATCG, crystal packing forces lead to a series of correlated changes: widened minor groove, large slide, low twist, and large rise. The difference in helical parameters between A-T steps lying within A-tracts, versus A-T steps within alternating AT sequences, demonstrates the importance of neighboring steps on the conformation of a given dinucleotide step.

A molecular dynamics simulation study of coaxial stacking in RNA

We report on unrestrained molecular dynamics simulations of an RNA tetramer binding to a tetra-nucleotide overhang at the 5'-end of an RNA hairpin (nicked structure) and of the corresponding continuous hairpin with Na+ as counterions. The simulations lead to stable structures and in this way a structural model for the coaxially stacked RNA hairpin is generated. The stacking interface in the coaxially stacked nicked hairpin structure is characterized by a reduced twist and shift and a slightly increased propeller twist as compared to the continuous system. This leads to an increased overlap between C22 and G23 in the stacking interface of the nicked structure. In the simulations the continuous RNA hairpin has an almost straight helical axis. On the other hand, the corresponding axis for the nicked structure exhibits a marked kink of 39 degrees. The stacking interface exhibits no increased flexibility as compared to the corresponding base pair step in the continuous structure.

Crystal studies of B-DNA: the answers and the questions

The history of DNA crystallography is reviewed and is followed by discussion of the methods used for structure determination. The features of B-DNA molecular and crystal structures are described.

Refined high-field NMR solution structure of a binary-addressed pyrene/perfluoro-azide complementary DNA oligonucleotide system shows extensive distortion in the central nick region

The structural analysis of the photoactivated binary system of complementary-addressing nucleic acid sequences (1:2:3) by high-resolution NMR spectroscopy and restrained molecular dynamics is reported. The binary system comprised a 12 base-pair target DNA sequence, pdGTATCAGTTTCT (1), and two hexanucleotides, (dAGAAACp-L-Az (2) and Pyr-pdTGATAC (3)), complementary to neighbouring sites in the target DNA. Oligonucleotide (2) is conjugated with a p-azidotetrafludrobenzyl group (Az) via a linker group (L), and the other oligonucleotide (3) is equipped with the photosensitizing pyrenyl-1-methylamino group (Pyr). We now extend the structural analysis of 1:2:3, which was previously based on qualitative 2D 1H-NMR data and thermodynamic analysis of complex formation from UV-visible thermal denaturation experiments. In the current work structural refinement was performed by separate molecular dynamics runs for six different starting structures based on 318 proton-proton distance-range constraints, evaluated from the 1H-NOESY spectrum (tau(mix) = 200 ms, 600 MHz) using complete relaxation matrix analysis (NMR/TRIAD/MARDIGRAS). Additional Watson-Crick hydrogen bond restraints were included in the calculations based on the detected signals from the exchangeable protons, using REFOPT(NY) methods. The final averaged structure obtained from the six refined co-ordinate sets showed a considerable degree of axis bend (62.5 degrees) with the bending point in the middle of the duplex in the region of the backbone nick between the two short oligonucleotides. The complex behaves dynamically as the equivalent of two short B-DNA-like duplexes displaying a hinge-like flexing at their junction. In all final structures the Pyr function location was very restricted, the aromatic group lying in the duplex minor groove near residues 4T, 5C and 2T. In contrast, the location of the perfluoroazido group was different in all the final structures, indicating the high flexibility of this group in the duplex. The only feature common to all six final azido group orientations was the outside location on the side of the major groove. The distance between the Pyr and Az groups varied from 11 A to 24 A for the six final structures (17 A, final average structure). The dynamics of duplex denaturation for 1:2:3 was probed by monitoring the temperature-induced NMR line broadening of the imino protons in a 1D variable temperature NMR experiment. The melting of 1:2:3 starts both from the ends and from the middle part of the duplex at the backbone break between the two short oligonucleotides reflecting the destabilisation of the pyrene-arylazido nick region in the duplex.

Internal correlations in minor groove profiles of experimental and computed B-DNA conformations

It has been noticed that converged conformations of B-DNA oligomers obtained in MD calculations often have very small atom position rmsd values from the canonical B-DNA and all helical parameters close to the standard values, but their minor grooves tend to be somewhat narrower. This apparent bias disappears, however, when C5' rather than phosphorus atoms are used for measuring the groove width. At the origin of this effect is the specific orientation of phosphate groups in the canonical B-DNA model which maximizes their separation across the minor groove. When measured by C5' traces, minor groove profiles of experimental structures available in the Nucleic Acids Database show much less tendency to narrow below the canonical width. Correlation analysis reveals a high degree of correspondence in shapes of minor grooves of calculated and experimental single-crystal structures of B-DNA oligomers.

SELEX and missing phosphate contact analyses reveal flexibility within the AP-2[alpha] protein: DNA binding complex

The AP-2 family of transcription factors are defined by the presence of a novel DNA binding domain, termed a 'basic helix-span-helix' motif. The AP-2 genes regulate important aspects of vertebrate embryogenesis and have also been linked to the control of cell proliferation and tumorigenesis, but the cellular targets that the AP-2 proteins control are largely undefined. In particular, since only a limited number of sequences have previously been utilized to define the nature of the AP-2 binding site, the range of DNA sequences recognized by the AP-2 proteins remains unknown. We have therefore utilized a SELEX analysis to identify multiple new AP-2[alpha] binding sites. Moreover, we have devised a novel missing phosphate and nucleotide competition analysis to characterize the residues in the binding site required for AP-2[alpha] protein:DNA contact. These studies suggest that the AP-2[alpha] protein:DNA complex is flexible and indicate that AP-2[alpha] can bind three related sequence motifs: GCC N3 GGC, GCC N4 GGC and GCC N3/4 GGG. The availability of these more refined consensus sequences should assist in the identification of target genes for this critical transcription factor.

DNA structure and flexibility in the sequence-specific binding of papillomavirus E2 proteins

The papillomavirus E2 proteins are transcriptional regulators that bind to a consensus DNA sequence ACCG NNNN CGGT. Multiple copies of this binding site are found in the viral genomes. The affinities of the naturally occurring binding sites for the E2 proteins are predominantly dependent upon the sequence of the NNNN spacer. The hierarchies of binding site affinities among the sites present in the viral genomes result in differential occupancy during the viral life-cycle. In turn, this differential binding regulates transcription from viral promoters, including those for the oncogenes E6 and E7. Structural and biochemical studies have shown that E2 proteins bend the DNA to which they specifically bind. Atomic resolution structures of complexes of the bovine papillomavirus strain 1 (BPV-1) E2 protein and DNA show that the protein does not contact the spacer DNA. A direct comparison of the binding of the DNA-binding domains of the E2 proteins from BPV-1 and human papillomavirus strain 16 (HPV-16) to a series of binding sites as a function of the sequence of their central spacer and/or the presence of a nick or gap in one strand of the spacer DNA is presented in this paper. The BPV-1 E2 DNA-binding domain is only moderately sensitive to the nature of the central spacer; less than several fold differences in affinity were observed when the DNA sequence of the spacer was varied and/or a nick or gap was introduced. In contrast, the HPV-16 E2 DNA-binding domain binds to sites containing A:T-rich central spacers with significantly increased affinity. The introduction of a nick or gap into the spacer of these high affinity sequences is very detrimental to HPV-16 E2 binding while comparable nicks or gaps have only small effects in the low affinity sequences. These results suggest that the HPV-16 E2 protein recognizes the structure of the DNA spacer and that the mechanism of DNA-sequence specific binding of the homologous HPV-16 E2 and BPV-1 E2 proteins is significantly different.

Structure of a DNA analog of the primer for HIV-1 RT second strand synthesis

The non-self-complementary DNA decamer C-A-A-A-G-A-A-A-A-G/C-T-T-T-T-C-T-T-T-G is a DNA/DNA analogue of a portion of the polypurine tract or PPT, which is a RNA/DNA hybrid that serves as a primer for synthesis of the (+) DNA strand by HIV reverse transcriptase (RT), and which is not digested by the RNase H domain of reverse transcriptase following (-) strand synthesis. The same unusual conformation that eludes RNase H, thought to be a change in width of minor groove, may also be responsible for the inhibition of HIV RT by minor groove binding drugs such as distamycin and their bis-linked derivatives. The present X-ray crystal structure of this DNA decamer exhibits the usual properties of A-tract B-DNA under biologically relevant conditions: large propeller twist of base-pairs, narrowed minor groove, and a straight helix axis. Groove narrowing is fully developed in the A-A-A-A region, but not in the A-A-A region, which previous investigators have proposed as being too short to exhibit typical A-tract properties. The RNA/DNA hybrid produced by HIV reverse transcriptase during (-) strand synthesis presumably forms a "heteromerous" or H-helix with narrower minor groove than an A-helical RNA/RNA duplex. If the narrowing of minor groove in A-tract H-helices is comparable to that seen in A-tract B-helices, then the narrowed minor groove of the polypurine tract could make the second primer site both (1) impervious to RNase H digestion, and (2) susceptible to inhibition by minor groove binding drugs.

A DNA decamer with a sticky end: the crystal structure of d-CGACGATCGT

The crystal structure of d-CGACGATCGT has been determined to a resolution of 2.6 A. The molecule was synthesized by standard phosphoramidite procedures, and purified by anion-exchange HPLC. Crystals are monolclinic, space group P2(1), with unit cell dimensions, a = 26.45 A, b = 34.66 A, c = 32.17 A, beta = 113.45 degrees and Z = 4, containing a B-DNA double helix in each crystallographic asymmetric unit. The structure was solved using molecular replacement, aided by an isomorphous derivative, in which a bromine atom was attached to the 5 position of cytosine 8. Problems of fit between the search model and the structure ultimately obtained necessitated the use of Patterson correlation procedures between the determination of the orientation and the translation of the molecule. In all, 69 solvent molecules have been identified, and the structure has been refined to an R-factor of 0.214, using the 1421 reflections with F > 2sigma(F), collected at -120 degrees C. The sequence produces a molecule containing eight Watson-Crick base-pairs and a two-nucleotide 5'-sticky end at each end of the duplex. The sticky ends cohere with one another, so the molecules form continuous 10-fold double helices throughout the crystal, with each strand being interrupted by inherent staggered nicks. The relative angular relationships between helices in the structure differ from each other; most of the arrangements differ from Holliday junctions, whose rotational orientations are phased by a crossover and which are modeled to contain double helices that are exactly parallel or antiparallel. However, one helical juxtaposition in this crystal is similar to the alignment of double helices in parallel Holliday junctions. A survey of DNA decamers that also form infinite helices in crystals reveals relationships that approximate both parallel and antiparallel Holliday junction alignments.

Analysis of local helix bending in crystal structures of DNA oligonucleotides and DNA-protein complexes

Sequence-dependent bending of the helical axes in 112 oligonucleotide duplex crystal structures resident in the Nucleic Acid Database have been analyzed and compared with the use of bending dials, a computer graphics tool. Our analysis includes structures of both A and B forms of DNA and considers both uncomplexed forms of the double helix as well as those bound to drugs and proteins. The patterns in bending preferences in the crystal structures are analyzed by base pair steps, and emerging trends are noted. Analysis of the 66 B-form structures in the Nucleic Acid Database indicates that uniform trends within all pyrimidine-purine and purine-pyrimidine steps are not necessarily observed but are found particularly at CG and GC steps of dodecamers. The results support the idea that AA steps are relatively straight and that larger roll bends occur at or near the junctions of these A-tracts with their flanking sequences. The data on 16 available crystal structures of protein-DNA complexes indicate that the majority of the DNA bends induced via protein binding are sharp localized kinks. The analysis of the 30 available A-form DNA structures indicates that these structures are also bent and show a definitive preference for bending into the deep major groove over the shallow minor groove.

Stable loop in the crystal structure of the intercalated four-stranded cytosine-rich metazoan telomere

In most metazoans, the telomeric cytosine-rich strand repeating sequence is d(TAACCC). The crystal structure of this sequence was solved to 1.9-A resolution. Four strands associate via the cytosine-containing parts to form a four-stranded intercalated structure held together by C.C+ hydrogen bonds. The base-paired strands are parallel to each other, and the two duplexes are intercalated into each other in opposite orientations. One TAA end forms a highly stabilized loop with the 5' thymine Hoogsteen-base-paired to the third adenine. The 5' end of this loop is in close proximity to the 3' end of one of the other intercalated cytosine strands. Instead of being entirely in a DNA duplex, this structure suggests the possibility of an alternative conformation for the cytosine-rich telomere strands.

Crystal and solution structures of d(CGC[e6G]AATTCGCG)-drug complexes reveal conformational polymorphism of O6-ethyl-guanine:cytosine base pair

O6-ethyl-guanine (e6G) is a relatively persistent alkylation lesion caused by the exposure of DNA to carcinogen N-ethyl-N-nitrosourea. We have studied the structural consequences of the e6G incorporation in DNA by X-ray crystallography and NMR. We have obtained crystals of the modified DNA dodecamer d(CGC[e6G]AATTCGCG) complexed to several minor groove binding drugs including Hoechst 33258, Hoechst 33342, netropsin, and SN6999. The space group of the crystals from those complexes is P2(1)2(1)2(1). However the crystal structure of the SN6999 complex is not isomorphous to that from the other three complexes. In all four refined crystal structures the drugs bind in the narrow minor groove at or close to the central AATT region of the dodecamer B-DNA duplex. The DNA conformation is influenced by the binding of drugs. The eight independent e6G:C base pairs have a conformation ranging from one with three-centered hydrogen bonds between the bases to a wobble conformation with two hydrogen bonds. The ethyl group of the eight e6G bases is mostly in the proximal orientation to N7. Our 1D and 2D-NMR studies of the same (free) dodecamer reveal that the e6G:C base pairs in the duplex are likely to adopt a wobble conformation in solution. Those results suggest that the e6G:C base pair has a dynamic equilibrium among various conformations, which may present an ambiguous signal to cells. In contrast, the e6G:T base pair adopts a Watson-Crick-like conformation. This may be a plausible explanation of why thymine is found preferentially incorporated across the e6G during replication.

Nonenzymatic plasmid ligation mediated by minor groove-binding molecules

DNA ligation is the weak link in the chain of gene cloning. We have developed a straightforward nonenzymatic alternative to this reaction that employs easily available commercial reagents. The method uses the affinity of distamycin for the minor groove to join DNA ends together. Phosphodiester bonds are formed after cyanoimidazole-promoted phosphate activation in the presence of manganese(II) cations. When transfected into eukaryotic cells, the chemically ligated plasmid is transcribed even more efficiently than after enzymatic ligation. At present this technique compares favorably with T4 ligation for AT-rich cohesive termini, but in principle it could be extended to any restriction site.

Electrophoretic evidence that single-stranded regions of one or more nucleotides dramatically increase the flexibility of DNA

The influence of single-stranded nicks and gaps on the flexibility of DNA has been investigated by subjecting to gel electrophoresis sets of molecules containing single-stranded regions of defined position and length. The DNA molecules were produced by ligating together synthetic oligomers that contained either nicks or single-stranded gaps of 1-4 nucleotides; the oligomer repeat lengths were 20, 21, 22, 23, or 26 bp, in order to produce nicks or gaps that were either in- or out-of-phase with the helix repeat of DNA. Nick-containing DNA molecules displayed nearly normal electrophoretic behavior, with maximum reductions in gel mobility (41 degrees C; 12% polyacrylamide gels) of approximately 10% for 230-bp molecules containing 10 nicks. In contrast, molecules containing gaps of 2-4 nucleotides demonstrated dramatic reductions in mobility, approaching one-half of the values of their full-duplex counterparts; molecules containing 1-nucleotide gaps displayed intermediate behavior. The observed (relative) mobilities of molecules containing gaps of more than 1 nucleotide were remarkably insensitive to temperature and to the presence of magnesium ions in the electrophoresis buffer. The central conclusion of the current study is that single-stranded gaps represent points of swivel-like character, whereas nicks retain much of the rigid character of double-helical DNA.

Minor groove binding of SN6999 to an alkylated DNA: molecular structure of d(CGC[e6G]AATTCGCG)-SN6999 complex

The interaction between a potent synthetic antitumor and antiviral minor groove binding drug 1-methyl-4-[4-[4-(4-(1-methylquinolinium)amino)benzamido]anilino] pyridinium dichloride (SN6999) and an alkylated DNA d(CGC[e6G]AATTCGCG) dodecamer has been studied by X-ray crystallography. The complex forms a new crystal lattice in the space group P2(1)2(1)2(1) with unit cell dimensions of a = 28.48 A, b = 36.11 A, and c = 69.60 A. The structure has been solved by the molecular replacement method and refined to an R-factor of 17.0% at approximately 2.5 A resolution using 1618 reflections. In the complex, the SN6999 covers almost six base pairs in the narrow minor groove with the 1-methylquinolinium (Q) ring near T8-A17 and the 1-methylpyridinium (P) ring near the C3-G22 base pair. The central benzamido (BQ) and anilino (BP) rings are essentially coplanar, with the Q and P rings having large dihedral angles of 38 degrees and 39 degrees, respectively, to the plane of BQ/BP. There is only one direct hydrogen bond between the amide NH of SN6999 to T20O2 of DNA. The drug-DNA interaction is stabilized by stacking interaction of sugar oxygens from T20O4' to BQ and C21O4' to BP. There is charge-induced dipole interaction between the positively charged nitrogen atom of 1-methylquinolinium with C9O4' and that of 1-methylpyridinium with G22O4'. The crystal structure of the complex can be used to explain the NMR results. SN6999 lacks the crescent shape observed in other minor groove binding drugs and distorts the DNA duplex upon binding. The complex packs in the lattice using the G-N2:G-N3 interlocking base pairs at both ends of the helix. As in earlier cases, the two independent e6G:C base pairs adopt different base pairing schemes. The e6G16:C9 base pair adopts a previously observed bifurcated configuration involving three-centered hydrogen bonds and is similar to a Watson-Crick pairing. In contrast, the e6G4:C21 base pair adopts a novel "reverse wobble" configuration with C21 being pushed toward the major groove side. The ethyl group is in the proximal orientation (to N7) in both base pairs. Taken together with the observations found in the same DNA complexed to Hoechst 33258, Hoechst 33342, and retropsin from different crystal lattices, the results suggest that the e6G:C base pairing is weak and polymorphic when compared to a normal G:C base pair and the DNA duplex containing this lesion is readily distorted.