Indicator of the Stacking Interaction in the DNA Double-Helical Structure: ChiraSac Study

Characterization of Z-DNA Using Circular Dichroism

The B-DNA to Z-DNA transition is a remarkable conformational change in DNA, which was originally observed in poly-GC DNA in the presence of high salt concentration. This eventually prompted the observation of the crystal structure of Z-DNA, a left-handed double-helical DNA, at atomic resolution. Despite advances in Z-DNA research, the application of circular dichroism (CD) spectroscopy as the fundamental technique to characterize this unique DNA conformation has remained constant. In this chapter, we describe a CD spectroscopic method for characterizing the B-DNA to Z-DNA transition of a CG-repeat double-stranded DNA fragment formed from a protein or chemical inducer.

Evaluating Geometric Definitions of Stacking for RNA Dinucleoside Monophosphates Using Molecular Mechanics Calculations

RNA modulation via small molecules is a novel approach in pharmacotherapies, where the determination of the structural properties of RNA motifs is considered a promising way to develop drugs capable of targeting RNA structures to control diseases. However, due to the complexity and dynamic nature of RNA molecules, the determination of RNA structures using experimental approaches is not always feasible, and computational models employing force fields can provide important insight. The quality of the force field will determine how well the predictions are compared to experimental observables. Stacking in nucleic acids is one such structural property, originating mainly from London dispersion forces, which are quantum mechanical and are included in molecular mechanics force fields through nonbonded interactions. Geometric descriptions are utilized to decide if two residues are stacked and hence to calculate the stacking free energies for RNA dinucleoside monophosphates (DNMPs) through statistical mechanics for comparison with experimental thermodynamics data. Here, we benchmark four different stacking definitions using molecular dynamics (MD) trajectories for 16 RNA DNMPs produced by two different force fields (RNA-IL and ff99OL3) and show that our stacking definition better correlates with the experimental thermodynamics data. While predictions within an accuracy of 0.2 kcal/mol at 300 K were observed in RNA CC, CU, UC, AG, GA, and GG, stacked states of purine-pyrimidine and pyrimidine-purine DNMPs, respectively, were typically underpredicted and overpredicted. Additionally, population distributions of RNA UU DNMPs were poorly predicted by both force fields, implying a requirement for further force field revisions. We further discuss the differences predicted by each RNA force field. Finally, we show that discrete path sampling (DPS) calculations can provide valuable information and complement the MD simulations. We propose the use of experimental thermodynamics data for RNA DNMPs as benchmarks for testing RNA force fields.

Sequence preference and structural heterogeneity of BZ junctions

BZ junctions, which connect B-DNA to Z-DNA, are necessary for local transformation of B-DNA to Z-DNA in the genome. However, the limited information on the junction-forming sequences and junction structures has led to a lack of understanding of the structural diversity and sequence preferences of BZ junctions. We determined three crystal structures of BZ junctions with diverse sequences followed by spectroscopic validation of DNA conformation. The structural features of the BZ junctions were well conserved regardless of sequences via the continuous base stacking through B-to-Z DNA with A-T base extrusion. However, the sequence-dependent structural heterogeneity of the junctions was also observed in base step parameters that are correlated with steric constraints imposed during Z-DNA formation. Further, circular dichroism and fluorescence-based analysis of BZ junctions revealed that a base extrusion was only found at the A-T base pair present next to a stable dinucleotide Z-DNA unit. Our findings suggest that Z-DNA formation in the genome is influenced by the sequence preference for BZ junctions.

Probing helical transitions in a DNA duplex

The complex conformational change from B-DNA to Z-DNA requires inversion of helix-handedness. Multiple degrees of freedom are intricately coupled during this transition, and formulating an appropriate reaction coordinate that captures the underlying complexity would be problematic. In this contribution, we adopt an alternative approach, based on the potential energy landscape perspective, to construct a kinetic transition network. Microscopic insight into the B → Z transition is provided in terms of geometrically defined discrete paths consisting of local minima and the transition states that connect them. We find that the inversion of handedness can occur via two competing mechanisms, either involving stretched intermediates, or a B-Z junction, in agreement with previous predictions. The organisation of the free energy landscape further suggests that this process is likely to be slow under physiological conditions. Our results represent a key step towards decoding the more intriguing features of the B → Z transition, such as the role of ionic strength and negative supercoiling in reshaping the landscape.

Accuracy of Td-DFT in the Ultraviolet and Circular Dichroism Spectra of Deoxyguanosine and Uridine

Accuracy of the time-dependent density functional theory (Td-DFT) was examined for the ultraviolet (UV) and circular dichroism (CD) spectra of deoxyguanosine (dG) and uridine, using 11 different DFT functionals and two different basis sets. The Td-DFT results of the UV and CD spectra were strongly dependent on the functionals used. The basis-set dependence was observed only for the CD spectral calculations. For the UV spectra, the B3LYP and PBE0 functionals gave relatively good results. For the CD spectra, the B3LYP and PBE0 with 6-311G(d,p) basis gave relatively permissible result only for dG. The results of other functionals were difficult to be used for the studies of the UV and CD spectra, though the symmetry adapted cluster-configuration interaction (SAC-CI) method reproduced well the experimental spectra of these molecules. To obtain valuable information from the theoretical calculations of the UV and CD spectra, the theoretical tool must be able to reproduce correctly both of the intensities and peak positions of the UV and CD spectra. Then, we can analyze the reasons of the changes of the intensity and/or the peak position to clarify the chemistry involved. It is difficult to recommend Td-DFT as such tools of science, at least from the examinations using dG and uridine.

Similarities and Differences between RNA and DNA Double-Helical Structures in Circular Dichroism Spectroscopy: A SAC-CI Study

The helical structures of DNA and RNA are investigated experimentally using circular dichroism (CD) spectroscopy. The signs and the shapes of the CD spectra are much different between the right- and left-handed structures as well as between DNA and RNA. The main difference lies in the sign at around 295 nm of the CD spectra: it is positive for the right-handed B-DNA and the left-handed Z-RNA but is negative for the left-handed Z-DNA and the right-handed A-RNA. We calculated the SAC-CI CD spectra of DNA and RNA using the tetramer models, which include both hydrogen-bonding and stacking interactions that are important in both DNA and RNA. The SAC-CI results reproduced the features at around 295 nm of the experimental CD spectra of each DNA and RNA, and elucidated that the strong stacking interaction between the two base pairs is the origin of the negative peaks at 295 nm of the CD spectra for both DNA and RNA. On the basis of these facts, we discuss the similarities and differences between RNA and DNA double-helical structures in the CD spectroscopy based on the ChiraSac methodology.

Structural and functional studies of a large winged Z-DNA-binding domain of Danio rerio protein kinase PKZ

The Z-DNA-binding domain of PKZ from zebrafish (Danio rerio; drZαPKZ ) contains the largest β-wing among known Z-DNA-binding domains. To elucidate the functional implication of the β-wing, we solved the crystal structure of apo-drZαPKZ . Structural comparison with its Z-DNA-bound form revealed a large conformational change within the β-wing during Z-DNA binding. Biochemical studies of protein mutants revealed that two basic residues in the β-wing are responsible for Z-DNA recognition as well as fast B-Z transition. Therefore, the extra basic residues in the β-wing of drZαPKZ are necessary for the fast B-Z transition activity.