A-like guanine-guanine stacking in the aqueous DNA duplex of d(GGGGCCCC)

Assessment of A- to B- DNA Transitions Utilizing the Drude Polarizable Force Field

In addition to the well-characterized B-form of DNA, duplex DNA can adopt various conformations, such as A or Z-DNA. Though less common, these structures can be induced biologically through protein or ligand interactions or experimentally with niche environmental conditions, such as high salt concentrations or in mixed water-ethanol. Reproducing these alternate structures through molecular dynamics simulations in recent years has been quite challenging with the currently available force fields, simulation techniques, and time scales. In this study, the Drude polarizable force field is tested for its ability to facilitate transitions between A-DNA and B-DNA or maintain A-DNA. Though transitions away from B-DNA were observed in high concentrations of ethanol, the resulting structures had hybrid properties taken from both B-DNA and A-DNA structures. This was also true for A-DNA in ethanol, which lost some of the A-DNA properties that it was expected to maintain. When B-DNA was tested in high salt environments, the resulting B-DNA structures showed no distinguishable differences with the increasing salt concentrations tested. These results with the Drude FF and recent results with additive force fields suggest that at present the current additive and polarizable force fields do not facilitate a complete transition between B- to A-DNA conformations under the conditions simulated. At present, the Drude FF favors A-B DNA hybrid structures when simulated in nonphysiological conditions.

Development of Mn2+-Specific Biosensor Using G-Quadruplex-Based DNA

Metal ions are used in various situations in living organisms and as a part of functional materials. Since the excessive intake of metal ions can cause health hazards and environmental pollution, the development of new molecules that can monitor metal ion concentrations with high sensitivity and selectivity is strongly desired. DNA can form various structures, and these structures and their properties have been used in a wide range of fields, including materials, sensors, and drugs. Guanine-rich sequences respond to metal ions and form G-quadruplex structures and G-wires, which are the self-assembling macromolecules of G-quadruplex structures. Therefore, guanine-rich DNA can be applied to a metal ion-detection sensor and functional materials. In this study, the IRDAptamer library originally designed based on G-quadruplex structures was used to screen for Mn²⁺, which is known to induce neurodegenerative diseases. Circular dichroism and fluorescence analysis using Thioflavin T showed that the identified IRDAptamer sequence designated MnG4C1 forms a non-canonical G-quadruplex structure in response to low concentrations of Mn²⁺. A serum resistance and thermostability analysis revealed that MnG4C1 acquired stability in a Mn²⁺-dependent manner. A Förster resonance energy transfer (FRET) system using fluorescent molecules attached to the termini of MnG4C1 showed that FRET was effectively induced based on Mn²⁺-dependent conformational changes, and the limit of detection (LOD) was 0.76 µM for Mn²⁺. These results suggested that MnG4C1 can be used as a novel DNA-based Mn²⁺-detecting molecule.

Thermodynamics of d(GGGGCCCC) Binding to Neomycin-Class Aminoglycosides

DNA adopts a number of conformations that can affect its binding to other macromolecules. The conformations (A, B, Z) can be sequence- and/or solution-dependent. While AT-rich DNA sequences generally adopt a Canonical B-form structure, GC-rich sequences are more promiscuous. Recognition of GC-rich nucleic acids by small molecules has been much more challenging than the recognition of AT-rich duplexes. Spectrophotometric and calorimetric techniques were used to characterize the binding of neomycin-class aminoglycosides to a GC-rich DNA duplex, G₄C₄, in various ionic and pH conditions. Our results reveal that binding enhances the thermal stability of G₄C₄, with thermal enhancement decreasing with increasing pH and/or Na⁺ concentration. Although G₄C₄ bound to aminoglycosides demonstrated a mixed A- and B-form conformation, circular dichroism studies indicate that binding induces a conformational shift toward A-form DNA. Isothermal titration calorimetry studies reveal that aminoglycoside binding to G₄C₄ is linked to the uptake of protons at pH = 7.0 and that this uptake is pH-dependent. Increased pH and/or Na⁺ concentration results in a decrease in G₄C₄ affinity for the aminoglycosides. The binding affinities of the aminoglycosides follow the expected hierarchy: neomycin > paromomycin > ribostamycin. The salt dependence of DNA binding affinities of aminoglycosides is consistent with at least two drug NH₃⁺ groups participating in electrostatic interactions with G₄C₄. These studies further embellish our understanding of the many factors facilitating recognition of GC-rich DNA structures as guided by their optimum charge and shape complementarity for small-molecule amino sugars.

Structural properties and binding mechanism of DNA aptamers sensing saliva melatonin for diagnosis and monitoring of circadian clock and sleep disorders

Circadian desynchrony with the external light-dark cycle influences the rhythmic secretion of melatonin which is among the first signs of circadian rhythm sleep disorders. An accurate dim light melatonin onset (established indicator of circadian rhythm sleep disorders) measurement requires lengthy assays, and antibody affinities alterations, especially in patients with circadian rhythm disorders whose melatonin salivary levels vary significantly, making antibodies detection mostly inadequate. In contrast, aptamers with their numerous advantages (e.g., target selectivity, structural flexibility in tuning binding affinities, small size, etc.) can become preferable biorecognition molecules for salivary melatonin detection with high sensitivity and specificity. This study thoroughly characterizes the structural property and binding mechanism of a single-stranded DNA aptamer full sequence (MLT-C-1) and its truncated versions (MLT-A-2, MLT-A-4) to decipher its optimal characteristics for saliva melatonin detection. We use circular dichroism spectroscopy to determine aptamers' conformational changes under different ionic strengths and showed that aptamers display a hairpin loop structure where few base pairs in the stem play a significant role in melatonin binding and formation of aptamer stabilized structure. Through microscale thermophoresis, aptamers demonstrated a high binding affinity in saliva samples (MLT-C-1F K_d = 12.5 ± 1.7 nM; MLT-A-4F K_d = 11.2 ± 1.6 nM; MLT-A-2F K_d = 2.4 ± 2.8 nM; limit-of-detection achieved in pM, highest sensitivity attained for MLT-A-2F aptamer with the lowest detection limit of 1.35 pM). Our data suggest that aptamers are promising as biorecognition molecules and provide the baseline parameters for the development of an aptamer-based point-of-care diagnostic system for melatonin detection and accurate profiling of its fluctuations in saliva.

Mapping DNA Conformations Using Single-Molecule Conductance Measurements

DNA is an attractive material for a range of applications in nanoscience and nanotechnology, and it has recently been demonstrated that the electronic properties of DNA are uniquely sensitive to its sequence and structure, opening new opportunities for the development of electronic DNA biosensors. In this report, we examine the origin of multiple conductance peaks that can occur during single-molecule break-junction (SMBJ)-based conductance measurements on DNA. We demonstrate that these peaks originate from the presence of multiple DNA conformations within the solutions, in particular, double-stranded B-form DNA (dsDNA) and G-quadruplex structures. Using a combination of circular dichroism (CD) spectroscopy, computational approaches, sequence and environmental controls, and single-molecule conductance measurements, we disentangle the conductance information and demonstrate that specific conductance values come from specific conformations of the DNA and that the occurrence of these peaks can be controlled by controlling the local environment. In addition, we demonstrate that conductance measurements are uniquely sensitive to identifying these conformations in solutions and that multiple configurations can be detected in solutions over an extremely large concentration range, opening new possibilities for examining low-probability DNA conformations in solutions.

Variability of the DNA Backbone Geometry in DNA-Protein Complexes: Experimental Data Analysis

We have analyzed and compared the available experimental data (PDB) on the backbone geometry of the DNA in solution (NMR), in crystals (X-rays), and in complexes with proteins (X-rays and cryo-electron microscopy). The deoxyribose (pseudorotational angle τ₀) and ε/ζ (BI-BII transition in phosphates) flexibilities are practically the same in the four samples. The α/γ mobility is minimal in crystalline DNA: on the histograms, there is one canonical and one noncanonical t/t peak. The α/γ mobility increases in DNA solutions (three more noncanonical peaks) and is maximal in DNA-protein complexes (another additional peak). On a large amount of data, we have confirmed that the three main degrees of freedom of the sugar-phosphate backbone are "orthogonal": changes in any of the angles τ₀, (ζ-ε), and (γ-α) occur, as a rule, at a constant (usually canonical) value of any other. In the DNA-protein complexes, none of the geometrical parameters commonly used to distinguish the A and B forms of DNA, except for Zp and its simpler analog Zp', show an unambiguous correlation with τ₀. Proteins, binding to DNA, in 59% of cases change the local shape of the helix up to the characteristic of the A-form without switching the deoxyribose conformation from south to north. However, we have found simple local characteristics of one nucleotide that correlate with the angles τ₀ and (ζ-ε). These are the angles C3'C1'N* and C4'C3'P(2), respectively. They are orthogonal in DNA-protein complexes exactly as the pair τ₀ and (ζ-ε). Most characteristics of DNA in complexes with proteins are the same in X-ray and in cryo-EM data, except for the histogram for the angle τ₀. We offer a possible explanation for this difference. We also discuss the artifacts on the ε/ζ histogram for DNA in solutions caused by the currently used NMR refinement protocols.

Harnessing intrinsic fluorescence for typing of secondary structures of DNA

High-throughput investigation of structural diversity of nucleic acids is hampered by the lack of suitable label-free methods, combining fast and cheap experimental workflow with high information content. Here, we explore the use of intrinsic fluorescence emitted by nucleic acids for this scope. After a preliminary assessment of suitability of this phenomenon for tracking conformational changes of DNA, we examined steady-state emission spectra of an 89-membered set of oligonucleotides with reported conformation (G-quadruplexes (G4s), i-motifs, single- and double-strands) by means of multivariate analysis. Principal component analysis of emission spectra resulted in successful clustering of oligonucleotides into three corresponding conformational groups, without discrimination between single- and double-stranded structures. Linear discriminant analysis was exploited for the assessment of novel sequences, allowing the evaluation of their G4-forming propensity. Our method does not require any labeling agent or dye, avoiding the related bias, and can be utilized to screen novel sequences of interest in a high-throughput and cost-effective manner. In addition, we observed that left-handed (Z-) G4 structures were systematically more fluorescent than most other G4 structures, almost reaching the quantum yield of 5'-d[(G3T)3G3]-3' (G3T, the most fluorescent G4 structure reported to date).

Discrimination between G/C Binding Sites by Olivomycin A Is Determined by Kinetics of the Drug-DNA Interaction

Olivomycin A (OA) exerts its cytotoxic potency due to binding to the minor groove of the G/C-rich DNA and interfering with replication and transcription. Screening of the complete set of tetranucleotide G/C sites by electrophoretic mobility gel shift assay (EMSA) revealed that the sites containing central GC or GG dinucleotides were able to bind OA, whereas the sites with the central CG dinucleotide were not. However, studies of equilibrium OA binding in solution by fluorescence, circular dichroism and isothermal titration calorimetry failed to confirm the sequence preference of OA, indicating instead a similar type of complex and comparable affinity of OA to all G/C binding sites. This discrepancy was resolved by kinetics analysis of the drug–DNA interaction: the dissociation rate significantly differed between SGCS, SGGS and SCGS sites (S stands for G or C), thereby explaining the disintegration of the complexes during EMSA. The functional relevance of the revealed differential kinetics of OA–DNA interaction was demonstrated in an in vitro transcription assay. These findings emphasize the crucial role of kinetics in the mechanism of OA action and provide an important approach to the screening of new drug candidates.

Concerted dynamics of metallo-base pairs in an A/B-form helical transition

Metal-mediated base pairs expand the repertoire of nucleic acid structures and dynamics. Here we report solution structures and dynamics of duplex DNA containing two all-natural C-Hg^II-T metallo base pairs separated by six canonical base pairs. NMR experiments reveal a 3:1 ratio of well-resolved structures in dynamic equilibrium. The major species contains two (N3)T-Hg^II-(N3)C base pairs in a predominantly B-form helix. The minor species contains (N3)T-Hg^II-(N4)C base pairs and greater A-form characteristics. Ten-fold different ¹J coupling constants (¹⁵N,¹⁹⁹Hg) are observed for (N3)C-Hg^II (114 Hz) versus (N4)C-Hg^II (1052 Hz) connectivities, reflecting differences in cytosine ionization and metal-bonding strengths. Dynamic interconversion between the two types of C-Hg^II-T base pairs are coupled to a global conformational exchange between the helices. These observations inspired the design of a repetitive DNA sequence capable of undergoing a global B-to-A-form helical transition upon adding Hg^II, demonstrating that C-Hg^II-T has unique switching potential in DNA-based materials and devices.

Piperine Modulates Protein Mediated Toxicity in Fragile X-Associated Tremor/Ataxia Syndrome through Interacting Expanded CGG Repeat (r(CGG)exp) RNA

An expansion of CGG tandem repeats in the 5' untranslated region (5'-UTR) of fragile X mental retardation 1 (FMR1) gene causes fragile X-associated tremor/ataxia syndrome (FXTAS). The transcripts of these expanded repeats r(CGG)^exp either form RNA foci or undergo the repeat-associated non-ATG (RAN) translation that produces toxic homopolymeric proteins in neuronal cells. The discovery of small molecule modulators that possess a strong binding affinity and high selectivity to these toxic expanded repeats RNA could be a promising therapeutic approach to cure the expanded repeat-associated neurological diseases. Therefore, here we sought to test the therapeutic potential of a natural alkaloid, piperine, by assessing its ability to bind and neutralize the toxicity of r(CGG)^exp RNA motif. To accomplish this first, we have determined the affinity of piperine to r(CGG)^exp RNA using fluorescence-based binding assay and isothermal titration calorimetry assay. These assays showed that piperine forms a thermodynamically favorable interaction with r(CGG)^exp RNA with high selectivity to the G-rich RNA motif. Interaction of piperine with r(CGG)^exp motif was further validated using several biophysical techniques such as CD, CD melting, NMR spectroscopy, and gel retardation assay. Moreover, piperine was also found to be effective for improving the r(CGG)^exp associated splicing defects and RAN translation in a FXTAS cell model system. Our results effectively provided the evidence that piperine strongly interacts with r(CGG)^exp RNA and could be used as a suitable candidate for therapeutic development against FXTAS.

CD Study of the G-Quadruplex Conformation

Circular Dichroic (CD) spectroscopy is one of the most frequently used methods for guanine quadruplex studies and in general for studies of conformational properties of nucleic acids. The reason is its high sensitivity to even slight changes in mutual orientation of absorbing bases of DNA. CD can reveal formation of particular structural DNA arrangements and can be used to search for the conditions stabilizing the structures, to follow the transitions between various structural states, to explore kinetics of their appearance, to determine thermodynamic parameters, and also to detect formation of higher order structures. CD spectroscopy is an important complementary technique to NMR spectroscopy and X-ray diffraction in quadruplex studies due to its sensitivity, easy manipulation of studied samples, and relative inexpensiveness. In this part, we present the protocol for the use of CD spectroscopy in the study of guanine quadruplexes, together with practical advice and cautions about various, particularly interpretation, difficulties.

A common intronic variant of PARP1 confers melanoma risk and mediates melanocyte growth via regulation of MITF

Previous genome-wide association studies have identified a melanoma-associated locus at 1q42.1 that encompasses a ∼100-kb region spanning the PARP1 gene. Expression quantitative trait locus (eQTL) analysis in multiple cell types of the melanocytic lineage consistently demonstrated that the 1q42.1 melanoma risk allele (rs3219090[G]) is correlated with higher PARP1 levels. In silico fine-mapping and functional validation identified a common intronic indel, rs144361550 (-/GGGCCC; r2 = 0.947 with rs3219090), as displaying allele-specific transcriptional activity. A proteomic screen identified RECQL as binding to rs144361550 in an allele-preferential manner. In human primary melanocytes, PARP1 promoted cell proliferation and rescued BRAFV600E-induced senescence phenotypes in a PARylation-independent manner. PARP1 also transformed TERT-immortalized melanocytes expressing BRAFV600E. PARP1-mediated senescence rescue was accompanied by transcriptional activation of the melanocyte-lineage survival oncogene MITF, highlighting a new role for PARP1 in melanomagenesis.

A bouquet of DNA structures: Emerging diversity

Structural polymorphism of DNA has constantly been evolving from the time of illustration of the double helical model of DNA by Watson and Crick. A variety of non-canonical DNA structures have constantly been documented across the globe. DNA attracted worldwide attention as a carrier of genetic information. In addition to the classical Watson–Crick duplex, DNA can actually adopt diverse structures during its active participation in cellular processes like replication, transcription, recombination and repair. Structures like hairpin, cruciform, triplex, G-triplex, quadruplex, i-motif and other alternative non-canonical DNA structures have been studied at length and have also shown their in vivo occurrence. This review mainly focuses on non-canonical structures adopted by DNA oligonucleotides which have certain prerequisites for their formation in terms of sequence, its length, number and orientation of strands along with varied solution conditions. This conformational polymorphism of DNA might be the basis of different functional properties of a specific set of DNA sequences, further giving some insights for various extremely complicated biological phenomena. Many of these structures have already shown their linkages with diseases like cancer and genetic disorders, hence making them an extremely striking target for structure-specific drug designing and therapeutic applications.

•

DNA can adopt diverse range of structures other than classical Watson–Crick duplex.

•

Discussion of alternate structures like hairpin, cruciform, triplex, quadruplex etc.

•

This review gives some insights for the biological relevance of DNA structures.

Telomeric repeats act as nucleosome-disfavouring sequences in vivo

Telomeric DNAs consist of tandem repeats of G-clusters such as TTAGGG and TG_1-3, which are the human and yeast repeat sequences, respectively. In the yeast Saccharomyces cerevisiae, the telomeric repeats are non-nucleosomal, whereas in humans, they are organized in tightly packaged nucleosomes. However, previous in vitro studies revealed that the binding affinities of human and yeast telomeric repeat sequences to histone octamers in vitro were similar, which is apparently inconsistent with the differences in the human and yeast telomeric chromatin structures. To further investigate the relationship between telomeric sequences and chromatin structure, we examined the effect of telomeric repeats on the formation of positioned nucleosomes in vivo by indirect end-label mapping, primer extension mapping and nucleosome repeat analyses, using a defined minichromosome in yeast cells. We found that the human and yeast telomeric repeat sequences both disfavour nucleosome assembly and alter nucleosome positioning in the yeast minichromosome. We further demonstrated that the G-clusters in the telomeric repeats are required for the nucleosome-disfavouring properties. Thus, our results suggest that this inherent structural feature of the telomeric repeat sequences is involved in the functional dynamics of the telomeric chromatin structure.

A short GC-rich palindrome of human mannose receptor gene coding region displays a conformational switch

Conformational switching in DNA is fundamental to biological processes. The structural status of a palindromic GC-rich dodecamer DNA sequence, integral part of human MRC2 coding region, and a related sequence of opposite polarity from human FDX1 gene were characterized and compared. UV-melting, circular dichroism, and gel electrophoresis experiments demonstrated the formation of intermolecular structures. Although stability and molecularity of both the oligomeric structures were found to be almost identical, their secondary structures differed remarkably as A1 MRC2 sequence showed A-like and B-like DNA conformation, whereas the A2 FDX1 sequence exhibited only the A-like signatures. The study is relevant for understanding structural polymorphism at genomic locations depending on DNA sequence and solution environment.

Circular dichroism spectroscopy of DNA: from duplexes to quadruplexes

Nucleic acids bear the genetic information and participate in its expression and evolution during replication, repair, recombination, transcription, and translation. These phenomena are mostly based on recognition of nucleic acids by proteins. The major factor enabling the specific recognition is structure. Circular dichroism (CD) spectroscopy is very useful to study secondary structures of nucleic acids, in general, and DNA, in particular. CD sensitively reflects isomerizations among distinct conformational states. The isomerizations may operate as molecular switches regulating various physiological or pathological processes. Here, we review CD spectra of nucleic acids, beginning with early studies on natural DNA molecules through analyses of synthetic polynucleotides to study of selected genomic fragments.

Conformational changes in DNA upon ligand binding monitored by circular dichroism

Circular dichroism (CD) spectroscopy is an optical technique that measures the difference in the absorption of left and right circularly polarized light. This technique has been widely employed in the studies of nucleic acids structures and the use of it to monitor conformational polymorphism of DNA has grown tremendously in the past few decades. DNA may undergo conformational changes to B-form, A-form, Z-form, quadruplexes, triplexes and other structures as a result of the binding process to different compounds. Here we review the recent CD spectroscopic studies of the induction of DNA conformational changes by different ligands, which includes metal derivative complex of aureolic family drugs, actinomycin D, neomycin, cisplatin, and polyamine. It is clear that CD spectroscopy is extremely sensitive and relatively inexpensive, as compared with other techniques. These studies show that CD spectroscopy is a powerful technique to monitor DNA conformational changes resulting from drug binding and also shows its potential to be a drug-screening platform in the future.

Thermodynamics of nucleic acid "shape readout" by an aminosugar

Recognition of nucleic acids is important for our understanding of nucleic acid structure as well as for our understanding of nucleic acid-protein interactions. In addition to the direct readout mechanisms of nucleic acids such as H-bonding, shape recognition of nucleic acids is being increasingly recognized as playing an equally important role in DNA recognition. Competition dialysis, UV, flourescent intercalator displacement (FID), computational docking, and calorimetry studies were conducted to study the interaction of neomycin with a variety of nucleic acid conformations (shapes). At pH 5.5, the results suggest the following. (1) Neomycin binds three RNA structures [16S A site rRNA, poly(rA)·poly(rA), and poly(rA)·poly(rU)] with high affinities (K(a) ~ 10(7) M(-1)). (2) The binding of neomycin to A-form GC-rich oligomer d(A(2)G(15)C(15)T(2))(2) has an affinity comparable to those of RNA structures. (3) The binding of neomycin to DNA·RNA hybrids shows a 3-fold variance that can be attributed to their structural differences [for poly(dA)·poly(rU), K(a) = 9.4 × 10(6) M(-1), and for poly(rA)·poly(dT), K(a) = 3.1 × 10(6) M(-1)]. (4) The interaction of neomycin with DNA triplex poly(dA)·2poly(dT) yields a binding affinity (K(a)) of 2.4 × 10(5) M(-1). (5) Poly(dA-dT)(2) shows the lowest association constant for all nucleic acids studied (K(a) < 10(5)). (6) Neomycin binds to G-quadruplexes with K(a) values of ~10(4)-10(5) M(-1). (7) Computational studies show that the decrease in major groove width in the B to A transition correlates with increasing neomycin affinity. Neomycin's affinity for various nucleic acid structures can be ranked as follows: RNAs and GC-rich d(A(2)G(15)C(15)T(2))(2) structures > poly(dA)·poly(rU) > poly(rA)·poly(dT) > T·A-T triplex, G-quadruplex, B-form AT-rich, or GC-rich DNA sequences. The results illustrate the first example of a small molecule-based "shape readout" of different nucleic acid conformations.

Sole and stable RNA duplexes of G-rich sequences located in the 5'-untranslated region of protooncogenes

Guanine- (G-) rich nucleic acid sequences can form four-stranded structures called G-quadruplexes. It is widely held that the formation of a G-quadruplex in RNA is more feasible than in DNA because of the lack of a complementary strand in mRNA. Here, we analyzed sequences of 5'-untranslated regions of protooncogenes and surprisingly found that these regions showed an enrichment of not only guanine (G) but also cytosine (C) nucleotides. Since neighboring cytosine- (C-) rich regions can affect the formation and stability of a G-quadruplex structure, we further investigated the properties of DNA and RNA structures of G-rich and GC-rich regions. We selected typical GC-rich RNA sequences from protooncogenes and corresponding DNA sequences and investigated their structures. It was found that the GC-rich RNA sequences formed stable A-form duplexes as their major structure independent of the surrounding conditions, including the presence of different cations (Na(+), K(+), or Li(+)) or molecular crowding with 40 wt % poly(ethylene glycol) with an average molecular mass of 200 Da although there are a few exceptions in which only a combination of K(+) and molecular crowding induced a G-quadruplex structure of an extremely G-rich RNA sequence. In contrast, structural polymorphisms involving duplexes, G-quadruplexes, and i-motifs were observed for GC-rich DNA sequences depending on the surrounding factors. These results demonstrate the considerable structural and functional differences in GC-rich sequences of the genome (DNA) and transcriptosome (mRNA) with respect to the nucleic acid backbone. Moreover, it was suggested that structural study for a G-rich RNA sequence should be carried out under cell-mimicking condition where K(+) and crowding cosolutes exist.

Circular dichroism and conformational polymorphism of DNA

Here we review studies that provided important information about conformational properties of DNA using circular dichroic (CD) spectroscopy. The conformational properties include the B-family of structures, A-form, Z-form, guanine quadruplexes, cytosine quadruplexes, triplexes and other less characterized structures. CD spectroscopy is extremely sensitive and relatively inexpensive. This fast and simple method can be used at low- as well as high-DNA concentrations and with short- as well as long-DNA molecules. The samples can easily be titrated with various agents to cause conformational isomerizations of DNA. The course of detected CD spectral changes makes possible to distinguish between gradual changes within a single DNA conformation and cooperative isomerizations between discrete structural states. It enables measuring kinetics of the appearance of particular conformers and determination of their thermodynamic parameters. In careful hands, CD spectroscopy is a valuable tool for mapping conformational properties of particular DNA molecules. Due to its numerous advantages, CD spectroscopy significantly participated in all basic conformational findings on DNA.

Hot-spot consensus of fluoroquinolone-mediated DNA cleavage by Gram-negative and Gram-positive type II DNA topoisomerases

Bacterial DNA gyrase and topoisomerase IV are selective targets of fluoroquinolones. Topoisomerase IV versus gyrase and Gram-positive versus Gram-negative behavior was studied based on the different recognition of DNA sequences by topoisomerase-quinolone complexes. A careful statistical analysis of preferred bases was performed on a large number (>400) of cleavage sites. We found discrete preferred sequences that were similar when using different enzymes (i.e. gyrase and topoisomerase IV) from the same bacterial source, but in part diverse when employing enzymes from different origins (i.e. Escherichia coli and Streptococcus pneumoniae). Subsequent analysis on the wild-type and mutated consensus sequences showed that: (i) Gn/Cn-rich sequences at and around the cleavage site are hot spots for quinolone-mediated strand breaks, especially for E. coli topoisomerases: we elucidated positions required for quinolone and enzyme recognition; (ii) for S. pneumoniae enzymes only, A and T at positions -2 and +6 are discriminating cleavage determinants; (iii) symmetry of the target sequence is a key trait to promote cleavage and (iv) the consensus sequence adopts a heteronomous A/B conformation, which may trigger DNA processing by the enzyme-drug complex.

Intramolecular and intermolecular guanine quadruplexes of DNA in aqueous salt and ethanol solutions

DNA guanine quadruplexes are all based on stacks of guanine tetrads, but they can be of many types differing by mutual strand orientation, topology, position and structure of loops, and the number of DNA molecules constituting their structure. Here we have studied a series of nine DNA fragments (G(3)Xn)(3)G(3), where X = A, C or T, and n = 1, 2 or 3, to find how the particular bases and their numbers enable folding of the molecule into quadruplex and what type of quadruplex is formed. We show that any single base between G(3) blocks gives rise to only four-molecular parallel-stranded quadruplexes in water solutions. In contrast to previous models, even two Ts in potential loops lead to tetramolecular parallel quadruplexes and only three consecutive Ts lead to an intramolecular quadruplex, which is antiparallel. Adenines make the DNA less prone to quadruplex formation. (G(3)A(2))(3)G(3) folds into an intramolecular antiparallel quadruplex. The same is true with (G(3)A(3))(3)G(3) but only in KCl. In NaCl or LiCl, (G(3)A(3))(3)G(3) prefers to generate homoduplexes. Cytosine still more interferes with the quadruplex, which only is generated by (G(3)C)(3)G(3), whereas (G(3)C(2))(3)G(3) and (G(3)C(3))(3)G(3) generate hairpins and/or homoduplexes. Ethanol is a more potent DNA guanine quadruplex inducer than are ions in water solutions. It promotes intramolecular folding and parallel orientation of quadruplex strands, which rather corresponds to quadruplex structures observed in crystals.

Vibrational dynamics of DNA. III. Molecular dynamics simulations of DNA in water and theoretical calculations of the two-dimensional vibrational spectra

A theoretical description of the vibrational excitons in DNA is presented by using the vibrational basis mode theory developed in Papers I and II. The parameters obtained from the density functional theory calculations, such as vibrational coupling constants and basis mode frequencies, are used to numerically simulate two-dimensional (2D) IR spectra of dG(n):dC(n) and dA(n):dT(n) double helices with n varying from 1 to 10. From the molecular dynamics simulations of dG(5)C(5) and dA(5)T(5) double helices in D(2)O solution, it is found that the thermally driven internal motions of these systems in an aqueous solution do not induce strong fluctuations of basis mode frequencies nor vibrational couplings. In order to construct the two-exciton Hamiltonian, the vibrational anharmonicities of eight basis modes are obtained by carrying out B3LYP6-31G(*) calculations for the nine basis modes. The simulated 2D IR spectra of dG(n):dC(n) double helix in D(2)O solution are directly compared with closely related experimental results. The 2D IR spectra of dG(n):dC(n) and dA(n):dT(n) are found to be weakly dependent on the number of base pairs. The present work demonstrates that the computational procedure combining quantum chemistry calculation and molecular dynamics simulation methods can be of use to predict 2D IR spectra of nucleic acids in solutions.

Vibrational dynamics of DNA: IV. Vibrational spectroscopic characteristics of A-, B-, and Z-form DNA’s

Linear and nonlinear IR spectroscopic studies of nucleic acids can provide crucial information on solution conformations of DNA double helix and its complex with other molecules. Carrying out density functional theory calculations of A-, B-, and Z-form DNA's, the authors obtained vibrational spectroscopic properties as well as coupling constants between different basis modes. The vibrational couplings that determine the extent of exciton delocalization are strongly dependent on DNA conformation mainly because the interlayer distance between two neighboring base pairs changes with respect to the DNA conformation. The Z-DNA has comparatively small interlayer vibrational coupling constants so that its vibrational spectrum depends little on the number of base pairs, whereas the A-DNA shows a notable dependency on the size. Furthermore, it is shown that a few distinctively different line shape changes in both IR and two-dimensional IR spectra as the DNA conformation changes from B to A or from B to Z can be used as marker bands and characteristic features distinguishing different DNA conformations.

Interpreting DNA Vibrational Circular Dichroism Spectra Using a Coupling Model from Two-Dimensional Infrared Spectroscopy

Two-dimensional infrared spectroscopy was recently used to measure the vibrational couplings between carbonyl bonds located on DNA nucleobases (Krummel, A. T.; Mukherjee, P.; Zanni, M. T. J. Phys. Chem. B 2003, 107, 9165 and Krummel, A. T.; Zanni, M. T. J. Phys. Chem. B 2006, 110, 13991). Here, we extend the coupling model derived from these 2D IR experiments to simulate the vibrational absorption and vibrational circular dichroism (VCD) spectra of three double-stranded DNA oligomers: poly(dG)-poly(dC), poly(dG-dC), and dGGCC. Using this model, we determine that the VCD spectrum of A-form poly(dG)-poly(dC) is dominated by interactions between stacked bases, whereas the coupling between base pairs and stacked bases carries equal importance in the VCD spectrum of B-form poly(dG-dC). We also simulate the absorption and VCD spectra of dGGCC, which is a combination of A- and B-form configurations. These simulations give insight into the structural interpretation of VCD and absorption spectroscopies that have long been used to monitor DNA secondary structure and kinetics.

Photochemical probing of the B--a conformational transition in a linearized pUC19 DNA and its polylinker region

We induced the B-to-A conformational transition by ethanol in a linearized pUC19 DNA. A primer extension method was used in combination with UV light irradiation to follow the transition, based on pausing of DNA synthesis due to the presence of damaged bases in the template. Primer extension data highly correlated with the results of another method monitoring the B-A transition, i.e. inhibition of restriction endonuclease cleavage of UV light-irradiated DNA. Primer extension enabled us to locate damaged nucleotides within the region of interest. Most damaged nucleotides were located in B-form trimers, exclusively containing both pyrimidine bases (TTC, TCT, CTC, and CTT), and in a cytosine tetramer. The amount of damaged bases decreased in the course of B-A transition. Some of the damage even disappeared in the A-form, which mainly concerns the C(4) and C(3) blocks. The cleavage was nearly restored in the A-form within this region (Eco88I). On the contrary the decrease of damage was less significant with thymine dimers, only dropping to 50-60% of the B-form level. Consequently, the cleavage with EcoRI and HindIII remained mostly as before the transition (75% and 60% of uncleaved DNA preserved). We found significant differences in the B- and A-form pattern of UV light-damaged bases within the same region (polylinker) of DNA embedded within long (plasmid) or short (127 bp fragment) DNA molecules. The B-A transition of the fragment was found less cooperative than with linearized plasmid, which was confirmed by both CD spectroscopy and restriction cleavage inhibition.

Structural polymorphism exhibited by a quasipalindrome present in the locus control region (LCR) of the human beta-globin gene cluster

Structural polymorphism of DNA is a widely accepted property. A simple addition to this perception has been our recent finding, where a single nucleotide polymorphism (SNP) site present in a quasipalindromic sequence of beta-globin LCR exhibited a hairpin-duplex equilibrium. Our current studies explore that secondary structures adopted by individual complementary strands compete with formation of a perfect duplex. Using gel-electrophoresis, ultraviolet (UV)-thermal denaturation, circular dichroism (CD) techniques, we have demonstrated the structural transitions within a perfect duplex containing 11 bp quasipalindromic stretch (TGGGG(G/C)CCCCA), to hairpins and bulge duplex forms. The extended version of the 11 bp duplex, flanked by 5 bp on both sides also demonstrated conformational equilibrium between duplex and hairpin species. Gel-electrophoresis confirms that the duplex coexists with hairpin and bulge duplex/cruciform species. Further, in CD spectra of duplexes, presence of two overlapping positive peaks at 265 and 285 nm suggest the features of A- as well as B-type DNA conformation and show oligomer concentration dependence, manifested in A --> B transition. This indicates the possibility of an architectural switching at quasipalindromic region between linear duplex to a cruciform structure. Such DNA structural variations are likely to be found in the mechanics of molecular recognition and manipulation by proteins.

DNA Vibrational Coupling Revealed with Two-Dimensional Infrared Spectroscopy: Insight into Why Vibrational Spectroscopy Is Sensitive to DNA Structure

Two-dimensional infrared (2D IR) spectroscopy was used to study the carbonyl vibrational modes of guanine and cytosine bases in A- and B-form DNA. Located between 1600 and 1700 cm(-1), these modes are often used to monitor DNA secondary structure with traditional infrared spectroscopies such as FTIR, but traditional spectroscopies lack the necessary observables to unravel the coupling mechanisms that make these modes sensitive to secondary structure. By using 2D IR spectroscopy and electronic structure calculations on d(G(5)C(5)) and d(GC)(8) model nucleic acids, we find that hydrogen-bonded guanine/cytosine base pairs are primarily electrostatically coupled and that the coupling between these modes can be modeled with a transition dipole density approach. In comparison, electrostatics is insufficient to model stacked bases because of cooperative charge-sharing effects, but the coupling can be accurately calculated using a finite difference method. We find that the coupling is very strong for both hydrogen-bonded and stacked base geometries, creating vibrational modes that extend both across the base pairs and along the lengths of the helices. Our results provide a physical basis for understanding how strong coupling gives rise to the empirically established relationship between infrared spectroscopy and DNA/RNA secondary structure.

Molecular basis for temperature sensing by an RNA thermometer

Regulatory RNA elements, like riboswitches, respond to intracellular signals by three-dimensional (3D) conformational changes. RNA thermometers employ a similar strategy to sense temperature changes in the cell and regulate the translational machinery. We present here the first 3D NMR structure of the functional domain of a highly conserved bacterial RNA thermometer containing the ribosome binding site that remains occluded at normal temperatures (30 degrees C). We identified a region adjacent to the Shine-Dalgarno sequence that has a network of weak hydrogen bonds within the RNA helix. With the onset of heat shock at 42 degrees C, destabilisation of the RNA structure initiates at this region and favours the release of the ribosome binding site and of the start codon. Deletion of a highly conserved G residue leads to the formation of a stable regular RNA helix that loses thermosensing ability. Our results indicate that RNA thermometers are able to sense temperature changes without the aid of accessory factors.

Interaction of bulged DNA with leucine-containing mimics of NCS-chrom

Synthesis of chiral spirocyclic helical compounds containing leucine that mimic the molecular architecture of the potent DNA bulge binder obtained from the natural product metabolite NCSi-gb has been accomplished. The interaction between the compounds and DNA was studied by circular dichroism (CD) method. The results suggested that the two synthetic diastereoisomers specifically targeted the bulge site of DNA and induced conformational change of bulged DNA greatly.

NMR structure of the apoB mRNA stem-loop and its interaction with the C to U editing APOBEC1 complementary factor

We have solved the NMR structure of the 31-nucleotide (nt) apoB mRNA stem-loop, a substrate of the cytidine deaminase APOBEC1. We found that the edited base located at the 5' end of the octa-loop is stacked between two adenosines in both the unedited (cytidine 6666) and the edited (uridine 6666) forms and that the rest of the loop is unstructured. The 11-nt "mooring" sequence essential for editing is partially flexible although it is mostly in the stem of the RNA. The octa-loop and the internal loop in the middle of the stem confer this flexibility. These findings shed light on why APOBEC1 alone cannot edit efficiently the cytidine 6666 under physiological conditions, the editing base being buried in the loop and not directly accessible. We also show that APOBEC1 does not specifically bind apoB mRNA and requires the auxiliary factor, APOBEC1 complementary factor (ACF), to edit specifically cytidine 6666. The binding of ACF to both the mooring sequence and APOBEC1 explains the specificity of the reaction. Our NMR study lead us to propose a mechanism in which ACF recognizes first the flexible nucleotides of the mooring sequence (the internal loop and the 3' end octa-loop) and subsequently melts the stem-loop, exposing the amino group of the cytidine 6666 to APOBEC1. Thus, the flexibility of the mooring sequence plays a central role in the RNA recognition by ACF.

Mapping the B-A conformational transition along plasmid DNA

A simple method is presented to monitor conformational isomerizations along genomic DNA. We illustrate properties of the method with the B-A conformational transition induced by ethanol in linearized pUC19 plasmid DNA. At various ethanol concentrations, the DNA was irradiated with ultraviolet light, transferred to a restriction endonuclease buffer and the irradiated DNA was cleaved by 17 restriction endonucleases. The irradiation damaged DNA and the damage blocked the restrictase cleavage. The amount of uncleaved, i.e. damaged, DNA depended on the concentration of ethanol in a characteristic S-shape way typical of the cooperative B-A transition. The transition beginning and midpoint were determined for each restriction endonuclease. These data map the B-A transition along the whole polylinker of pUC19 DNA and six evenly distributed recognition sequences within the rest of the plasmid. The transition midpoints fell within the B-A transition region of the plasmid simultaneously determined by CD spectroscopy. The present method complements the previous methods used to study the B-A transition. It can be employed to analyze multikilobase regions of genomic DNA whose restriction endonuclease cleavage fragments can be separated and quantified on agarose gels.

Hairpin-duplex equilibrium reflected in the A-->B transition in an undecamer quasi-palindrome present in the locus control region of the human beta-globin gene cluster

Our recent work on an A-->G single nucleotide polymorphism (SNP) at the quasi-palindromic sequence d(TGGGG[A/G]CCCCA) of HS4 of the human beta-globin locus control region in an Indian population showed a significant association between the G allele and the occurrence of beta-thalassemia. Using UV-thermal denaturation, gel assay, circular dichroism (CD) and nuclease digestion experiments we have demonstrated that the undecamer quasi- palindromic sequence d(TGGGGACCCCA) (HPA11) and its reported polymorphic (SNP) version d(TGG GGGCCCCA) (HPG11) exist in hairpin-duplex equilibria. The biphasic nature of the melting profiles for both the oligonucleotides persisted at low as well as high salt concentrations. The HPG11 hairpin showed a higher T(m) than HPA11. The presence of unimolecular and bimolecular species was also shown by non-denaturating gel electrophoresis experiments. The CD spectra of both oligonucleotides showed features of the A- as well as B-type conformations and, moreover, exhibited a concentration dependence. The disappearance of the 265 nm positive CD signal in an oligomer concentration-dependent manner is indicative of an A-->B transition. The results give unprecedented insight into the in vitro structure of the quasi-palindromic sequence and provide the first report in which a hairpin-duplex equilibrium has been correlated with an A-->B interconversion of DNA. The nuclease-dependent degradation suggests that HPG11 is more resistant to nuclease than HPA11. Multiple sequence alignment of the HS4 region of the beta-globin gene cluster from different organisms revealed that this quasi-palindromic stretch is unique to Homo sapiens. We propose that quasi-palindromic sequences may form stable mini- hairpins or cruciforms in the HS4 region and might play a role in regulating beta-globin gene expression by affecting the binding of transcription factors.

Formation pathways of a guanine-quadruplex DNA revealed by molecular dynamics and thermodynamic analysis of the substates

The formation of a cation-stabilized guanine quadruplex (G-DNA) stem is an exceptionally slow process involving complex kinetics that has not yet been characterized at atomic resolution. Here, we investigate the formation of a parallel stranded G-DNA stem consisting of four strands of d(GGGG) using molecular dynamics simulations with explicit inclusion of counterions and solvent. Due to the limitations imposed by the nanosecond timescale of the simulations, rather than watching for the spontaneous formation of G-DNA, our approach probes the stability of possible supramolecular intermediates (including two-, three-, and four-stranded assemblies with out-of-register base pairing between guanines) on the formation pathway. The simulations suggest that "cross-like" two-stranded assemblies may serve as nucleation centers in the initial formation of parallel stranded G-DNA quadruplexes, proceeding through a series of rearrangements involving trapping of cations, association of additional strands, and progressive slippage of strands toward the full stem. To supplement the analysis, approximate free energies of the models are obtained with explicit consideration of the integral cations. The approach applied here serves as a prototype for qualitatively investigating other G-DNA molecules using molecular dynamics simulation and free-energy analysis.

Aminoglycoside (neomycin) preference is for A-form nucleic acids, not just RNA: results from a competition dialysis study

Aminoglycosides have been at the forefront of antimicrobial therapy for the past 50 years. Their specificity is believed to lie in binding duplex RNAs (rRNA). Competition dialysis studies of various nucleic acid forms with 9-aminoacridine, quinacrine, and a neomycin-acridine conjugate were carried out. Our results suggest a strong preference for aminoglycoside binding to nucleic acids that can adopt an A-type conformation. These results challenge the common belief that aminoglycoside specificity is simply for duplex RNAs.

Circular dichroism spectroscopy of conformers of (guanine + adenine) repeat strands of DNA

(Guanine+adenine) strands of DNA are known to associate into guanine tetraplexes, homodimerize into parallel or antiparallel duplexes, and fold into a cooperatively melting single strand resembling the protein alpha helix. Using CD spectroscopy and other methods, we studied how this conformational polymorphism depended on the primary structure of DNA. The study showed that d(GGGA)(5) and d(GGA)(7) associated into homoduplexes at low salt or in the presence of LiCl but were prone to guanine tetraplex formation, especially in the presence of KCl. In addition, they yielded essentially the same CD spectrum in the presence of ethanol as observed with the ordered single strand of d(GA)(10). Strands of d(GA)(10), d(GGAA)(5), d(GAA)(7), and d(GAAA)(5) associated into homoduplexes in both LiCl and KCl solutions, but not into guanine tetraplexes. d(GAAA)(5) and d(GAA)(7) further failed to form the single-stranded conformer in aqueous ethanol. Adenine protonation, however, stabilized the single-stranded conformer even in these adenine-rich fragments. The ordered single strands, homoduplexes as well as the guanine tetraplexes, all provided strikingly similar CD spectra, indicating that all of the conformers shared similar base stacking geometries. The increasing adenine content only decreased the conformer thermostability.

The structure of DNA within cationic lipid/DNA complexes

The structure of DNA within CLDCs used for gene delivery is controversial. Previous studies using CD have been interpreted to indicate that the DNA is converted from normal B to C form in complexes. This investigation reexamines this interpretation using CD of model complexes, FTIR as well as Raman spectroscopy and molecular dynamics simulations to address this issue. CD spectra of supercoiled plasmid DNA undergo a significant loss of rotational strength in the signal near 275 nm upon interaction with either the cationic lipid dimethyldioctadecylammonium bromide or 1,2-dioleoyltrimethylammonium propane. This loss of rotational strength is shown, however, by both FTIR and Raman spectroscopy to occur within the parameters of the B-type conformation. Contributions of absorption flattening and differential scattering to the CD spectra of complexes are unable to account for the observed spectra. Model studies of the CD of complexes prepared from synthetic oligonucleotides of varying length suggest that significant reductions in rotational strength can occur within short stretches of DNA. Furthermore, some alteration in the hydrogen bonding of bases within CLDCs is indicated in the FTIR and Raman spectroscopy results. In addition, alterations in base stacking interactions as well as hydrogen bonding are suggested by molecular dynamics simulations. A global interpretation of all of the data suggests the DNA component of CLDCs remains in a variant B form in which base/base interactions are perturbed.

Genomic occurrence of microsatellites containing integral and non-integral repeat numbers

We calculated occurrences of all dinucleotide and trinucleotide microsatellites in the human, mouse, and yeast genomes. The microsatellites were considered separately not only according to the repeated dinucleotide or trinucleotide and the microsatellite length but also according to the starting/terminal nucleotide. The analysis showed that dramatically non-equal amounts occurred in the human genome of microsatellites that differed only by the terminal nucleotides. For example, the 23-mer (TTG)(7)TT occurs 635 times in the human genome whereas (GTT)(7)GT is present only three times in the human genome though the two 23-mers share a 22 nucleotide sequence. The dramatically non-equal occurrences of microsatellites differing only by the terminal nucleotides are observed for most dinucleotide and trinucleotide microsatellites and in all analyzed genomes. We suppose that the strikingly non-equal genomic occurrences of these closely related microsatellites originate from conformational properties of DNA.

The S. cerevisiae architectural HMGB protein NHP6A complexed with DNA: DNA and protein conformational changes upon binding

NHP6A is a non-sequence-specific DNA-binding protein from Saccharomyces cerevisiae which belongs to the HMGB protein family. Previously, we have solved the structure of NHP6A in the absence of DNA and modeled its interaction with DNA. Here, we present the refined solution structures of the NHP6A-DNA complex as well as the free 15bp DNA. Both the free and bound forms of the protein adopt the typical L-shaped HMGB domain fold. The DNA in the complex undergoes significant structural rearrangement from its free form while the protein shows smaller but significant conformational changes in the complex. Structural and mutational analysis as well as comparison of the complex with the free DNA provides insight into the factors that contribute to binding site selection and DNA deformations in the complex. Further insight into the amino acid determinants of DNA binding by HMGB domain proteins is given by a correlation study of NHP6A and 32 other HMGB domains belonging to both the DNA-sequence-specific and non-sequence-specific families of HMGB proteins. The resulting correlations can be rationalized by comparison of solved structures of HMGB proteins.

Mediation of the A/B-DNA helix transition by G-tracts in the crystal structure of duplex CATGGGCCCATG

The crystal structure of the DNA dodecamer duplex CATGGGCCCATG lies on a structural continuum along the transition between A- and B-DNA. The dodecamer possesses the normal vector plot and inclination values typical of B-DNA, but has the crystal packing, helical twist, groove width, sugar pucker, slide and x-displacement values typical of A-DNA. The structure shows highly ordered water structures, such as a double spine of water molecules against each side of the major groove, stabilizing the GC base pairs in an A-like conformation. The different hydration of GC and AT base pairs provides a physical basis for solvent-dependent facilitation of the A<-->B helix transition by GC base pairs. Crystal structures of CATGGGCCCATG and other A/B-DNA intermediates support a 'slide first, roll later' mechanism for the B-->A helix transition. In the distribution of helical parameters in protein-DNA crystal structures, GpG base steps show A-like properties, reflecting their innate predisposition for the A conformation.

Circular dichroism spectroscopy reveals invariant conformation of guanine runs in DNA

We demonstrate that the characteristic circular dichroism (CD) features of the parallel-stranded DNA tetraplex of d(G4), especially the strong band at 260 nm, are characteristic for the B and A forms of the antiparallel duplex of d(C4G4). Hence, this band evidently originates from intrastrand guanine-guanine stacking, which is therefore very similar in the duplex and tetraplex DNA. In addition, the same type of the CD spectrum is provided by the ordered single strand of d(GA)10. This observation suggests that the ordered single strand of d(GA)10 is stabilized by a core of guanines stacked like in the parallel tetraplex. This view is used to start the modeling of the molecular structure of the ordered d(GA)10 single strand. Our studies suggest that guanine itself is strong enough to stabilize various secondary structures of DNA, which is a property relevant to thinking about the origin and evolution of molecular replicators.

Simulations of nucleic acids and their complexes

Recent years have seen considerable progress in simulations of nucleic acids. Improvements in force fields, simulation techniques and protocols, and increasing computer power have all contributed to making nanosecond-scale simulations of both DNA and RNA commonplace. The results are already helping to explain how nucleic acids respond to their environment and to their base sequence and to reveal the factors underlying recognition processes by probing biologically important nucleic acid-protein interactions and medically important nucleic acid-drug complexation. This Account summarizes methodological progress and applications of molecular dynamics to nucleic acids over the past few years and tries to identify remaining challenges.

Structural basis of DNA recognition by anticancer antibiotics, chromomycin A(3), and mithramycin: roles of minor groove width and ligand flexibility

Anticancer antibiotics, chromomycin A(3) (CHR) and mithramycin (MTR), inhibit cellular processes like transcription and replication, by binding reversibly to double-stranded DNA via minor groove, in the presence of bivalent metal ions like Mg(2+) with GC base specificity. Here, we have attempted to assess the roles of two parameters-namely DNA groove dimension and flexibility of the ligand-in the structural recognition between the ligands, (drug)(2)Mg(2+) and DNA. For the purpose we have employed three synthetic oligonucleotides with minor groove width lying between B- and A-type structures as model DNA sequences: d(GCGCGCGC)(2) in B-form, d(CCGGCGCCGG)(2) in B-form with unusual wide minor groove, and (GGGGCCCC)(2) in A-form. Association of the (drug)(2)Mg(2+) with the oligomers have been probed using spectroscopic techniques like absorbance, fluorescence, and CD. The binding and thermodynamic parameters for the different association processes have also been characterized. Major conclusions from the above studies are as follows. Groove size of the oligomers influences the conformation of the bound ligand. A saccharide dependent variation in structural rigidity of the ligands, (MTR)(2)Mg(2+) and (CHR)(2)Mg(2+), has been observed that leads to differences in the energetics of recognition of the same DNA sequence by the two ligands. In contrast to (CHR)(2)Mg(2+), higher flexibility in (MTR)(2)Mg(2+) makes its conformation in the DNA bound form less sensitive to the groove dimension of DNA.