Atomic-resolution crystal structures of B-DNA reveal specific influences of divalent metal ions on conformation and packing

In silico studies of the interaction of the minor groove binder Hoechst 33258 with B-DNA

Interaction of the minor groove binder, Hoechst 33258, with the Dickerson-Drew DNA dodecamer sequence has been investigated using docking, MM/QM, MM/GBSA and molecular dynamics computations to study the modes of binding and the interactions responsible for the binding. Besides the original Hoechst 33258 ligand (HT), a total of 12 ionization and stereochemical states for the ligand are obtained at the physiological pH and have been docked into B-DNA. These states have one or the other or both benzimidazole rings in protonated states, apart from the piperazine nitrogen, which has a quaternary nitrogen in all the states. Most of these states are found to exhibit good docking scores and free energy of binding with B-DNA. The best docked state has been taken further for molecular dynamics simulations and compared with the original HT. This state is protonated at both benzimidazole rings besides the piperazine ring and hence has very highly negative coulombic interaction energy. In both cases, there are strong coulombic interactions, but these are offset by the almost equally unfavorable solvation energies. Thus, the nonpolar forces, particularly van der Waals contacts, dominate the interaction, and the polar interactions highlight subtle changes in the binding energies, leading to more highly protonated states having more negative binding energies.Communicated by Ramaswamy H. Sarma.

On Water Arrangements in Right- and Left-Handed DNA Structures

DNA requires hydration to maintain its structural integrity. Crystallographic analyses have enabled patterns of water arrangements to be visualized. We survey these water motifs in this review, focusing on left- and right-handed duplex and quadruplex DNAs, together with the i-motif. Common patterns of linear spines of water organization in grooves have been identified and are widely prevalent in right-handed duplexes and quadruplexes. By contrast, a left-handed quadruplex has a distinctive wheel of hydration populating the almost completely circular single groove in this structure.

Transport of double-stranded ribonucleic acids (dsRNA) and deoxyribonucleic acids (DNA) in sand and iron oxide-coated sand columns under varying solution chemistries

Assessing ecological risks associated with the use of genetically modified RNA interference crops demands an understanding of the fate of crop-released insecticidal double-stranded RNA (dsRNA) molecules in soils. We studied the adsorption of one dsRNA and two double-stranded DNA as model nucleic acids (NAs) during transport through sand- and iron oxide-coated sand (IOCS)-filled columns over a range of solution pH and ionic compositions. Consistent with NA-sand electrostatic repulsion, we observed only slight retention of NAs in sand columns. Conversely, pronounced NA retention in IOCS columns is consistent with strong and irreversible NA adsorption involving electrostatic attraction to and inner-sphere complex formation of NAs with iron oxide coatings. Adsorption of NAs to iron oxides revealed a fast and a slow kinetic adsorption regime, possibly caused by the excluded-area effect. Adsorption of NAs to sand and IOCS increased in the presence of dissolved Mg2+ and with increasing ionic strength, reflecting cation-bridging and screening of repulsive electrostatics, respectively. The co-solute phosphate and a pre-adsorbed dissolved organic matter isolate competitively suppressed dsRNA adsorption to IOCS. Similar adsorption characteristics of dsRNA and similarly sized DNA suggest that existing information on DNA adsorption to soil particles helps in predicting adsorption and fate of dsRNA molecules in soils.

Conformational Morphing by a DNA Analogue Featuring 7-Deazapurines and 5-Halogenpyrimidines and the Origins of Adenine-Tract Geometry

Several efforts are currently directed at the creation and cellular implementation of alternative genetic systems composed of pairing components that are orthogonal to the natural dA/dT and dG/dC base pairs. In an alternative approach, Watson-Crick-type pairing is conserved, but one or all of the four letters of the A, C, G, and T alphabet are substituted by modified components. Thus, all four nucleobases were altered to create halogenated deazanucleic acid (DZA): dA was replaced by 7-deaza-2'-deoxyadenosine (dzA), dG by 7-deaza-2'-deoxyguanosine (dzG), dC by 5-fluoro-2'-deoxycytidine (FdC), and dT by 5-chloro-2'-deoxyuridine (CldU). This base-pairing system was previously shown to retain function in Escherichia coli. Here, we analyze the stability, hydration, structure, and dynamics of a DZA Dickerson-Drew Dodecamer (DDD) of sequence 5'-FdC-dzG-FdC-dzG-dzA-dzA-CldU-CldU-FdC-dzG-FdC-dzG-3'. Contrary to similar stabilities of DDD and DZA-DDD, osmotic stressing revealed a dramatic loss of hydration for the DZA-DDD relative to that for the DDD. The parent DDD 5'-d(CGCGAATTCGCG)-3' features an A-tract, a run of adenosines uninterrupted by a TpA step, and exhibits a hallmark narrow minor groove. Crystal structures─in the presence of RNase H─and MD simulations show increased conformational plasticity ("morphing") of DZA-DDD relative to that of the DDD. The narrow dzA-tract minor groove in one structure widens to resemble that in canonical B-DNA in a second structure. These changes reflect an indirect consequence of altered DZA major groove electrostatics (less negatively polarized compared to that in DNA) and hydration (reduced compared to that in DNA). Therefore, chemical modifications outside the minor groove that lead to collapse of major groove electrostatics and hydration can modulate A-tract geometry.

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.

Molecular dynamics study of collective water vibrations in a DNA hydration shell

The structure of DNA double helix is stabilized by water molecules and metal counterions that form the ion-hydration shell around the macromolecule. Understanding the role of the ion-hydration shell in the physical mechanisms of the biological functioning of DNA requires detailed studies of its structure and dynamics at the atomistic level. In the present work, the study of collective vibrations of water molecules around the DNA double helix was performed within the framework of classical all-atom molecular dynamics methods. Calculating the vibrational density of states, the vibrations of water molecules in the low-frequency spectra ranged from [Formula: see text]30 to [Formula: see text]300 [Formula: see text] were analyzed for the case of different regions of the DNA double helix (minor groove, major groove, and phosphate groups). The analysis revealed significant differences in the collective vibrations behavior of water molecules in the DNA hydration shell, compared to the vibrations of bulk water. All low-frequency modes of the DNA ion-hydration shell are shifted by about 15-20 [Formula: see text] towards higher frequencies, which is more significant for water molecules in the minor groove region of the double helix. The interactions of water molecules with the atoms of the macromolecule induce intensity decrease of the mode of hydrogen-bond symmetrical stretching near 150 [Formula: see text], leading to the disappearance of this mode in the DNA spectra. The obtained results can provide an interpretation of the experimental data for DNA low-frequency spectra and may be important for the understanding of the processes of indirect protein-nucleic recognition.

Conformation and Pairing Properties of an O6-Methyl-2'-deoxyguanosine-Directed Benzimidazole Nucleoside Analog in Duplex DNA

O6-Methyl-2'-deoxyguanosine (O6-MeG) is one of the most common DNA lesions and arises as a consequence of both xenobiotic carcinogens and endogenous methylation by S-adenosylmethionine. O6-MeG frequently causes G-to-A mutations during DNA replication due to the misincorporation of dTTP and continued DNA synthesis. Efforts to detect DNA adducts such as O6-MeG, and to understand their impacts on DNA structure and function, have motivated the creation of nucleoside analogs with altered base moieties to afford a more favorable interaction with the adduct as compared to the unmodified nucleotide. Such analogs directed at O6-MeG include benzimidazolinone and benzimidazole nucleotides, as well as their extended π surface analogs naphthimidazolinone and napthimidazole derivatives. These analogs form a more stable pair with O6-MeG than with G, most likely due to a combination of H-bonding and stacking. While extending the π surface of the analogs enhances their performance as adduct-directed probes, the precise origins of the increased affinity between the synthetic analogs and O6-MeG remain unclear. To better understand relevant conformational and pairing properties, we used X-ray crystallography and analyzed the structures of the DNA duplexes with naphthimidazolinone inserted opposite G or O6-MeG. The structures reveal a complex interaction of the analog found either in an anti orientation and stacked inside the duplex, either above or below G or O6-MeG, or in a syn orientation and paired opposite G with formation of a single H-bond. The experimental structural data are consistent with the stabilizing effect of the synthetic analog observed in UV melting experiments and calculations and moreover reveal that the origin of these observations appears to be superior stacking between O6-MeG and the extended π system of the synthetic probe.

The Multiple States of Environmental DNA and What Is Known about Their Persistence in Aquatic Environments

Increased use of environmental DNA (eDNA) analysis for indirect species detection has spurred the need to understand eDNA persistence in the environment. Understanding the persistence of eDNA is complex because it exists in a mixture of different states (e.g., dissolved, particle adsorbed, intracellular, and intraorganellar), and each state is expected to have a specific decay rate that depends on environmental parameters. Thus, improving knowledge about eDNA conversion rates between states and the reactions that degrade eDNA in different states is needed. Here, we focus on eukaryotic extraorganismal eDNA, outline how water chemistry and suspended mineral particles likely affect conversion among each eDNA state, and indicate how environmental parameters affect persistence of states in the water column. On the basis of deducing these controlling parameters, we synthesized the eDNA literature to assess whether we could already derive a general understanding of eDNA states persisting in the environment. However, we found that these parameters are often not being measured or reported when measured, and in many cases very few experimental data exist from which to draw conclusions. Therefore, further study of how environmental parameters affect eDNA state conversion and eDNA decay in aquatic environments is needed. We recommend analytic controls that can be used during the processing of water to assess potential losses of different eDNA states if all were present in a water sample, and we outline future experimental work that would help determine the dominant eDNA states in water.

Peptide-Decorated DNA Nanostructures Promote Site-Specific Hydroxyapatite Growth

The guiding principle for mineralized tissue formation is that mineral growth occurs through the interaction of Ca²⁺ and phosphate ions with extracellular matrix (ECM) proteins. Recently, nanoengineered DNA structures have been proposed as mimics to ECM scaffolds. However, these principles have not been applied to mineralized tissues. Here, we describe DNA nanostructures, namely, a DNA nanotube and a DNA origami rectangle that are site specifically functionalized with a mineral-promoting "SSEE" peptide derived from ECM proteins present in mineralized tissues. In the presence of Ca²⁺ and phosphate ions (mineralizing conditions), site-specific calcium phosphate formation occurred on the DNA nanostructures. Amorphous calcium phosphate or hydroxyapatite was formed depending on the incubation time, shape of the DNA nanostructure, and amount of Ca²⁺ and phosphate ions present. The ability to design and control the growth of hydroxyapatite through nanoengineered scaffolds provides insights into the mechanisms that may occur during crystal nucleation and growth of mineralized tissues and can inspire mineralized tissue regeneration strategies.

Extended magnesium and calcium force field parameters for accurate ion-nucleic acid interactions in biomolecular simulations

Magnesium and calcium play an essential role in the folding and function of nucleic acids. To correctly describe their interactions with DNA and RNA in biomolecular simulations, an accurate parameterization is crucial. In most cases, the ion parameters are optimized based on a set of experimental solution properties such as solvation free energies, radial distribution functions, water exchange rates, and activity coefficient derivatives. However, the transferability of such bulk-optimized ion parameters to quantitatively describe biomolecular systems is limited. Here, we extend the applicability of our previous bulk-optimized parameters by including experimental binding affinities toward the phosphate oxygen on nucleic acids. In particular, we systematically adjust the combination rules that are an integral part of the pairwise interaction potentials of classical force fields. This allows us to quantitatively describe specific ion binding to nucleic acids without changing the solution properties in the most simple and efficient way. We show the advancement of the optimized Lorentz combination rule for two representative nucleic acid systems. For double-stranded DNA, the optimized combination rule for Ca²⁺ significantly improves the agreement with experiments, while the standard combination rule leads to unrealistically distorted DNA structures. For the add A-riboswitch, the optimized combination rule for Mg²⁺ improves the structure of two specifically bound Mg²⁺ ions as judged by the experimental distance to the binding site. Including experimental binding affinities toward specific ion binding sites on biomolecules, therefore, provides a promising perspective to develop a more accurate description of metal cations for biomolecular simulations.

Adsorption of double-stranded ribonucleic acids (dsRNA) to iron (oxyhydr-)oxide surfaces: comparative analysis of model dsRNA molecules and deoxyribonucleic acids (DNA)

Double-stranded ribonucleic acid (dsRNA) molecules are novel plant-incorporated protectants expressed in genetically modified RNA interference (RNAi) crops. Ecological risk assessment (ERA) of RNAi crops requires a heretofore-missing detailed understanding of dsRNA adsorption in soils, a key fate process. Herein, we systematically study the adsorption of a model dsRNA molecule and of two double-stranded deoxyribonucleic acid (DNA) molecules of varying lengths to three soil iron (oxyhydr-)oxides - goethite, lepidocrocite, and hematite - over a range of solution pH (4.5-10), ionic strength (I = 10-100 mM NaCl) and composition (0.5, 1, and 3 mM MgCl2) and in the absence and presence of phosphate (0.05-5 mM) as co-adsorbate. We hypothesized comparable adsorption characteristics of dsRNA and DNA based on their structural similarities. Consistently, the three nucleic acids (NAs) showed high adsorption affinities to the iron (oxyhydr-)oxides with decreasing adsorption in the order goethite, lepidocrocite, and hematite, likely reflecting a decrease in the hydroxyl group density and positive charges of the oxide surfaces in the same order. NA adsorption also decreased with increasing solution pH, consistent with weakening of NA electrostatic attraction to and inner-sphere complex formation with the iron (oxyhydr-)oxides surfaces as pH increased. Adsorbed NA concentrations increased with increasing I and in the presence of Mg2+, consistent with adsorbed NA molecules adopting more compact conformations. Strong NA-phosphate adsorption competition demonstrates that co-adsorbates need consideration in assessing dsRNA fate in soils. Comparable adsorption characteristics of dsRNA and DNA molecules to iron (oxyhydr-)oxides imply that information on DNA adsorption to soil particle surfaces can inform dsRNA ERA.

Targeting the ALS/FTD-associated A-DNA kink with anthracene-based metal complex causes DNA backbone straightening and groove contraction

The use of a small molecule compound to reduce toxic repeat RNA transcripts or their translated aberrant proteins to target repeat-expanded RNA/DNA with a G4C2 motif is a promising strategy to treat C9orf72-linked disorders. In this study, the crystal structures of DNA and RNA–DNA hybrid duplexes with the -GGGCCG- region as a G4C2 repeat motif were solved. Unusual groove widening and sharper bending of the G4C2 DNA duplex A-DNA conformation with B-form characteristics inside was observed. The G4C2 RNA–DNA hybrid duplex adopts a more typical rigid A form structure. Detailed structural analysis revealed that the G4C2 repeat motif of the DNA duplex exhibits a hydration shell and greater flexibility and serves as a ‘hot-spot’ for binding of the anthracene-based nickel complex, Ni^II(Chro)₂ (Chro = Chromomycin A3). In addition to the original GGCC recognition site, Ni^II(Chro)₂ has extended specificity and binds the flanked G:C base pairs of the GGCC core, resulting in minor groove contraction and straightening of the DNA backbone. We have also shown that Chro-metal complexes inhibit neuronal toxicity and suppresses locomotor deficits in a Drosophila model of C9orf72-associated ALS. The approach represents a new direction for drug discovery against ALS and FTD diseases by targeting G4C2 repeat motif DNA.

Water spines and networks in G-quadruplex structures

Quadruplex DNAs can fold into a variety of distinct topologies, depending in part on loop types and orientations of individual strands, as shown by high-resolution crystal and NMR structures. Crystal structures also show associated water molecules. We report here on an analysis of the hydration arrangements around selected folded quadruplex DNAs, which has revealed several prominent features that re-occur in related structures. Many of the primary-sphere water molecules are found in the grooves and loop regions of these structures. At least one groove in anti-parallel and hybrid quadruplex structures is long and narrow and contains an extensive spine of linked primary-sphere water molecules. This spine is analogous to but fundamentally distinct from the well-characterized spine observed in the minor groove of A/T-rich duplex DNA, in that every water molecule in the continuous quadruplex spines makes a direct hydrogen bond contact with groove atoms, principally phosphate oxygen atoms lining groove walls and guanine base nitrogen atoms on the groove floor. By contrast, parallel quadruplexes do not have extended grooves, but primary-sphere water molecules still cluster in them and are especially associated with the loops, helping to stabilize loop conformations.

DNA nanostructures as templates for biomineralization

Nature uses extracellular matrix scaffolds to organize biominerals into hierarchical structures over various length scales. This has inspired the design of biomimetic mineralization scaffolds, with DNA nanostructures being among the most promising. DNA nanotechnology makes use of molecular recognition to controllably give 1D, 2D and 3D nanostructures. The control we have over these structures makes them attractive templates for the synthesis of mineralized tissues, such as bones and teeth. In this Review, we first summarize recent work on the crystallization processes and structural features of biominerals on the nanoscale. We then describe self-assembled DNA nanostructures and come to the intersection of these two themes: recent applications of DNA templates in nanoscale biomineralization, a crucial process to regenerate mineralized tissues.

The influence of Ni2+ and other ions on the trigonal structure of DNA

We present a new structure of a DNA dodecamer obtained in the presence of Ni²⁺ ions. The DNA forms Ni-guanine cross-links between neighboring molecules. Our results show that an adequate dosage of Ni²⁺ may help to form well-defined DNA nanostructures. We also compare our structure with other dodecamers which present unique features and also crystallize in trigonal unit cells, strongly influenced by the counterions associated with DNA. In all cases, the DNA duplexes form parallel pseudo-helical columns in the crystal, similar to DNA-protamine and native DNA fibers.

Functionalized DNA nanostructures as scaffolds for guided mineralization

The field of DNA nanotechnology uses synthetic DNA strands as building blocks for designing complex shapes in one-, two- and three-dimensions. Here, we investigate whether DNA nanostructures are feasible platforms for the precise organization of polyaspartic acid (pAsp), a known mineral carrier, with a goal towards biomimetic mineralization for enamel regeneration. We describe the preparation of DNA-pAsp conjugates and their subsequent assembly into ordered nanostructures. Covalent attachment of pAsp to DNA was noted to hinder DNA nanostructure formation past a certain threshold (50% pAsp) when tested on a previously published DNA system. However, a simplified double stranded DNA system (3sDH system) was more robust and efficient in its pAsp incorporation. In addition, the 3sDH system was successful in organizing mineral inducing groups in one dimension at repeating intervals of 28.7 ± 4.0 nm, as determined by atomic force microscopy. Our results demonstrate that DNA nanostructures can be functionalized with pAsp and act as a platform to investigate guided mineralization.

Molecular Dynamics Simulation of Homo-DNA: The Role of Crystal Packing in Duplex Conformation

The (4′→6′)-linked DNA homolog 2′,3′-dideoxy-β-D-glucopyranosyl nucleic acid (dideoxy-glucose nucleic acid or homo-DNA) exhibits stable self-pairing of the Watson–Crick and reverse-Hoogsteen types, but does not cross-pair with DNA. Molecular modeling and NMR solution studies of homo-DNA duplexes pointed to a conformation that was nearly devoid of a twist and a stacking distance in excess of 4.5 Å. By contrast, the crystal structure of the homo-DNA octamer dd(CGAATTCG) revealed a right-handed duplex with average values for helical twist and rise of ca. 15° and 3.8 Å, respectively. Other key features of the structure were strongly inclined base-pair and backbone axes in the duplex with concomitant base-pair slide and cross-strand stacking, and the formation of a dimer across a crystallographic dyad with inter-duplex base swapping. To investigate the conformational flexibility of the homo-DNA duplex and a potential influence of lattice interactions on its geometry, we used molecular dynamics (MD) simulations of the crystallographically observed dimer of duplexes and an isolated duplex in the solution state. The dimer of duplexes showed limited conformational flexibility, and key parameters such as helical rise, twist, and base-pair slide exhibited only minor fluctuations. The single duplex was clearly more flexible by comparison and underwent partial unwinding, albeit without significant lengthening. Thus, base stacking was preserved in the isolated duplex and two adenosines extruded from the stack in the dimer of duplexes were reinserted into the duplex and pair with Ts in a Hoogsteen mode. Our results confirmed that efficient stacking in homo-DNA seen in the crystal structure of a dimer of duplexes was maintained in the separate duplex. Therefore, lattice interactions did not account for the different geometries of the homo-DNA duplex in the crystal and earlier models that resembled inclined ladders with large base-pair separations that precluded efficient stacking.

Effect of the chemical environment of the DNA guanine quadruplex on the free energy of binding of Na and K ions

Guanine quadruplex (G-quadruplex) structures play a vital role in stabilizing the DNA genome and in protecting healthy cells from transforming into cancer cells. The structural stability of G-quadruplexes is greatly enhanced by the binding of monovalent cations such as Na⁺ or K⁺ into the interior axial channel. We computationally study the free energy of binding of Na⁺ and K⁺ ions to two intramolecular G-quadruplexes that differ considerably in their degree of rigidity and the presence or absence of terminal nucleotides. The goal of our study is two-fold. On the one hand, we study the free energy of binding every ion, which complements the experimental findings that report the average free energy for replacing Na⁺ with K⁺ ions. On the other hand, we examine the role of the G-quadruplex structure in the binding free energy. In the study, we employ all-atom molecular dynamics simulations and the alchemical transformation method for the computation of the free energies. To compare the cation-dependent contribution to the structural stability of G-quadruplexes, we use a two-step approach to calculate the individual free energy difference ΔG of binding two Na⁺ and two K⁺ to two G-quadruplexes: the unimolecular DNA d[T₂GT₂(G₃T)₃] (Protein Data Bank ID 2M4P) and the human telomeric DNA d[AGGG(TTAGGG)₃] (PDB ID 1KF1). In contrast to the experimental studies that estimate the average free energy of binding, we find a varying difference of approximately 2-9 kcal/mol between the free energy contribution of binding the first and second cation, Na⁺ or K⁺. Furthermore, we found that the free energy of binding K⁺ is not affected by the chemical nature of the two quadruplexes. By contrast, Na⁺ showed dependency on the G-quadruplex structure; the relatively small size allows Na⁺ to explore larger configurational space than K⁺. Numerical results presented here may offer reference values for future design of cationic drug-like ligands that replace the metal ions in G-quadruplexes.

Comparative analysis of inosine-substituted duplex DNA by circular dichroism and X-ray crystallography

Leveraging structural biology tools, we report the results of experiments seeking to determine if the different mechanical properties of DNA polymers with base analog substitutions can be attributed, at least in part, to induced changes from classical B-form DNA. The underlying hypothesis is that different inherent bending and twisting flexibilities may characterize non-canonical B-DNA, so that it is inappropriate to interpret mechanical changes caused by base analog substitution as resulting simply from 'electrostatic' or 'base stacking' influences without considering the larger context of altered helical geometry. Circular dichroism spectra of inosine-substituted oligonucleotides and longer base-substituted DNAs in solution indicated non-canonical helical conformations, with the degree of deviation from a standard B-form geometry depending on the number of I⋅C pairs. X-ray diffraction of a highly inosine-substituted DNA decamer crystal (eight I⋅C and two A⋅T pairs) revealed an A-tract-like conformation with a uniformly narrow minor groove, reduced helical rise, and the majority of sugars adopting a C1'-exo (southeastern) conformation. This contrasts with the standard B-DNA geometry with C2'-endo sugar puckers (south conformation). In contrast, the crystal structure of a decamer with only four I⋅C pairs has a geometry similar to that of the reference duplex with eight G⋅C and two A⋅T pairs. The unique crystal geometry of the inosine-rich duplex is noteworthy given its unusual CD signature in solution and the altered mechanical properties of some inosine-containing DNAs.

Photochemical analysis of structural transitions in DNA liquid crystals reveals differences in spatial structure of DNA molecules organized in liquid crystalline form

The anisotropic shape of DNA molecules allows them to form lyotropic liquid crystals (LCs) at high concentrations. This liquid crystalline arrangement is also found in vivo (e.g., in bacteriophage capsids, bacteria or human sperm nuclei). However, the role of DNA liquid crystalline organization in living organisms still remains an open question. Here we show that in vitro, the DNA spatial structure is significantly changed in mesophases compared to non-organized DNA molecules. DNA LCs were prepared from pBluescript SK (pBSK) plasmid DNA and investigated by photochemical analysis of structural transitions (PhAST). We reveal significant differences in the probability of UV-induced pyrimidine dimer photoproduct formation at multiple loci on the DNA indicative of changes in major groove architecture.

Metal-DNA Interactions Improve signal in High-Resolution Melting of DNA for Species Differentiation of Plasmodium Parasite

The success of high-resolution melting (HRM) analysis for distinguishing similar DNAs with minor base mismatch differences is limited. Here, metal-mediated structural change in DNA has been exploited to amplify HRM signals leading to differentiation of target DNAs in an orthologous gene corresponding to four Plasmodium species. Conserved 26-mer ssDNAs from ldh gene of the four Plasmodium species were employed as targets. A capture probe (CP) that is fully complementary to the Plasmodium falciparum target (FT) and has two base mismatches each, with the targets of Plasmodium vivax (VT), Plasmodium malariae, (MT), and Plasmodium ovale (OT), was considered. The DNA duplexes were treated with metal ions for structural perturbation and analyzed by HRM. Distinct resolution of melting fluorescence signal in otherwise identical HRM profiles for each of the DNA duplexes was achieved by using Ca⁺² or Mg⁺² ions, where, Ca⁺² conferred higher resolution. The increase in resolution for CP-FT versus CP-OT, CP-FT versus CP-VT, CP-FT versus CP-MT, CP-VT versus CP-OT, and CP-MT versus CP-OT with Ca-DNA as compared to control was 67.3-, 20.4-, 22.0-, 10.9-, and 8.3-fold, respectively. The signal resolution was the highest at pH 8. The method could detect 0.25 pmol/µl of the target DNA. Structural analysis showed that Ca⁺² and Mg⁺² ions perturbed the structure of DNA. This perturbation helped to improve HRM signal resolution among DNA targets corresponding to the orthologous gene of four Plasmodium species. This novel approach has potential application not only for Plasmodium species-specific diagnosis but also for differentiation of DNAs with minor sequence variation.

Photochromic switching of the DNA helicity induced by azobenzene derivatives

The photochromic properties of azobenzene, involving conformational changes occurring upon interaction with light, provide an excellent tool to establish new ways of selective regulation applied to biosystems. We report here on the binding of two water-soluble 4-(phenylazo)benzoic acid derivatives (Azo-2N and Azo-3N) with double stranded DNA and demonstrate that the photoisomerization of Azo-3N leads to changes in DNA structure. In particular, we show that stabilization and destabilization of the B-DNA secondary structure can be photochemically induced in situ by light. This photo-triggered process is fully reversible and could be an alternative pathway to control a broad range of biological processes. Moreover, we found that the bicationic Azo-3N exhibited a higher DNA-binding constant than the monocationic Azo-2N pointing out that the number of positive charges along the photosensitive polyamines chain plays a pivotal role in stabilizing the photochrome-DNA complex.

Stacking interactions in RNA and DNA: Roll-slide energy hyperspace for ten unique dinucleotide steps

Understanding dinucleotide sequence directed structures of nuleic acids and their variability from experimental observation remained ineffective due to unavailability of statistically meaningful data. We have attempted to understand this from energy scan along twist, roll, and slide degrees of freedom which are mostly dependent on dinucleotide sequence using ab initio density functional theory. We have carried out stacking energy analysis in these dinucleotide parameter phase space for all ten unique dinucleotide steps in DNA and RNA using DFT-D by ωB97X-D/6-31G(2d,2p), which appears to satisfactorily explain conformational preferences for AU/AU step in our recent study. We show that values of roll, slide, and twist of most of the dinucleotide sequences in crystal structures fall in the low energy region. The minimum energy regions with large twist values are associated with the roll and slide values of B-DNA, whereas, smaller twist values correspond to higher stability to RNA and A-DNA like conformations. Incorporation of solvent effect by CPCM method could explain the preference shown by some sequences to occur in B-DNA or A-DNA conformations. Conformational preference of BII sub-state in B-DNA is preferentially displayed mainly by pyrimidine-purine steps and partly by purine-purine steps. The purine-pyrimidine steps show largest effect of 5-methyl group of thymine in stacking energy and the introduction of solvent reduces this effect significantly. These predicted structures and variabilities can explain the effect of sequence on DNA and RNA functionality.

Generating Crystallographic Models of DNA Dodecamers from Structures of RNase H:DNA Complexes

The DNA dodecamer 5'-d(CGCGAATTCGCG)-3' is arguably the best studied oligonucleotide and crystal structures of duplexes with this sequence account for a considerable portion of the total number of oligo-2'-deoxynucleotide structures determined over the last 30 years. The dodecamer has commonly served as a template to analyze the effects of sequence on DNA conformation, the conformational properties of chemically modified nucleotides, DNA-ligand interactions as well as water structure and DNA-cation binding. Although molecular replacement is the phasing method of choice given the large number of available models of the dodecamer, this strategy often fails as a result of conformational changes caused by chemical modification, mismatch pairs, or differing packing modes. Here, we describe an alternative approach to determine crystal structures of the dodecamer in cases where molecular replacement does not produce a solution or when crystals of the DNA alone cannot be grown. It is based on the discovery that many dodecamers of the above sequence can be readily co-crystallized with Bacillus halodurans RNase H, whereby the enzyme is unable to cleave the DNA. Determination of the structure of the complex using the protein portion as the search model yields a structural model of the DNA. Provided crystals of the DNA alone are also available, the DNA model from the complex then enables phasing their structures by molecular replacement.

In situ NMR measurement of macromolecule-bound metal ion concentrations

Many nucleic acids and proteins require divalent metal ions such as Mg(2+) and Ca(2+) for folding and function. The lipophilic alignment media frequently used as membrane mimetics also bind these divalent metals. Here we demonstrate the use of (31)P NMR spectrum of a metal ion chelator (deoxycytidine diphosphate) to measure the bound [Mg(2+)] and [Ca(2+)] in situ for several biological model systems at relatively high divalent ion concentrations (1-10 mM). This method represents a general approach to measuring divalent metal ion binding in NMR samples where the amount and type of metal ion added to the system is known.

The intrinsic mechanics of B-DNA in solution characterized by NMR

Experimental characterization of the structural couplings in free B-DNA in solution has been elusive, because of subtle effects that are challenging to tackle. Here, the exploitation of the NMR measurements collected on four dodecamers containing a substantial set of dinucleotide sequences provides new, consistent correlations revealing the DNA intrinsic mechanics. The difference between two successive residual dipolar couplings (ΔRDCs) involving C6/8-H6/8, C3′-H3′ and C4′-H4′ vectors are correlated to the ³¹P chemical shifts (δP), which reflect the populations of the BI and BII backbone states. The δPs are also correlated to the internucleotide distances (D_inter) involving H6/8, H2′ and H2″ protons. Calculations of NMR quantities on high resolution X-ray structures and controlled models of DNA enable to interpret these couplings: the studied ΔRDCs depend mostly on roll, while D_inter are mainly sensitive to twist or slide. Overall, these relations demonstrate how δP measurements inform on key inter base parameters, in addition to probe the BI↔BII backbone equilibrium, and shed new light into coordinated motions of phosphate groups and bases in free B-DNA in solution. Inspection of the 5′ and 3′ ends of the dodecamers also supplies new information on the fraying events, otherwise neglected.

Can A Denaturant Stabilize DNA? Pyridine Reverses DNA Denaturation in Acidic pH

The stability of DNA is highly dependent on the properties of the surrounding solvent, such as ionic strength, pH, and the presence of denaturants and osmolytes. Addition of pyridine is known to unfold DNA by replacing π-π stacking interactions between bases, stabilizing conformations in which the nucleotides are solvent exposed. We show here experimental and theoretical evidences that pyridine can change its role and in fact stabilize the DNA under acidic conditions. NMR spectroscopy and MD simulations demonstrate that the reversal in the denaturing role of pyridine is specific, and is related to its character as pseudo groove binder. The present study sheds light on the nature of DNA stability and on the relationship between DNA and solvent, with clear biotechnological implications.

Antibody-drug conjugates as novel anti-cancer chemotherapeutics

Over the past couple of decades, antibody-drug conjugates (ADCs) have revolutionized the field of cancer chemotherapy. Unlike conventional treatments that damage healthy tissues upon dose escalation, ADCs utilize monoclonal antibodies (mAbs) to specifically bind tumour-associated target antigens and deliver a highly potent cytotoxic agent. The synergistic combination of mAbs conjugated to small-molecule chemotherapeutics, via a stable linker, has given rise to an extremely efficacious class of anti-cancer drugs with an already large and rapidly growing clinical pipeline. The primary objective of this paper is to review current knowledge and latest developments in the field of ADCs. Upon intravenous administration, ADCs bind to their target antigens and are internalized through receptor-mediated endocytosis. This facilitates the subsequent release of the cytotoxin, which eventually leads to apoptotic cell death of the cancer cell. The three components of ADCs (mAb, linker and cytotoxin) affect the efficacy and toxicity of the conjugate. Optimizing each one, while enhancing the functionality of the ADC as a whole, has been one of the major considerations of ADC design and development. In addition to these, the choice of clinically relevant targets and the position and number of linkages have also been the key determinants of ADC efficacy. The only marketed ADCs, brentuximab vedotin and trastuzumab emtansine (T-DM1), have demonstrated their use against both haematological and solid malignancies respectively. The success of future ADCs relies on improving target selection, increasing cytotoxin potency, developing innovative linkers and overcoming drug resistance. As more research is conducted to tackle these issues, ADCs are likely to become part of the future of targeted cancer therapeutics.

Molecular modeling of nucleic Acid structure: electrostatics and solvation

This unit presents an overview of computer simulation techniques as applied to nucleic acid systems, ranging from simple in vacuo molecular modeling techniques to more complete all-atom molecular dynamics treatments that include an explicit representation of the environment. The third in a series of four units, this unit focuses on critical issues in solvation and the treatment of electrostatics. UNITS 7.5 & 7.8 introduced the modeling of nucleic acid structure at the molecular level. This included a discussion of how to generate an initial model, how to evaluate the utility or reliability of a given model, and ultimately how to manipulate this model to better understand its structure, dynamics, and interactions. Subject to an appropriate representation of the energy, such as a specifically parameterized empirical force field, the techniques of minimization and Monte Carlo simulation, as well as molecular dynamics (MD) methods, were introduced as a way of sampling conformational space for a better understanding of the relevance of a given model. This discussion highlighted the major limitations with modeling in general. When sampling conformational space effectively, difficult issues are encountered, such as multiple minima or conformational sampling problems, and accurately representing the underlying energy of interaction. In order to provide a realistic model of the underlying energetics for nucleic acids in their native environments, it is crucial to include some representation of solvation (by water) and also to properly treat the electrostatic interactions. These subjects are discussed in detail in this unit.

Detection of an en masse and reversible B- to A-DNA conformational transition in prokaryotes in response to desiccation

The role that DNA conformation plays in the biochemistry of cells has been the subject of intensive research since DNA polymorphism was discovered. B-DNA has long been considered the native form of DNA in cells although alternative conformations of DNA are thought to occur transiently and along short tracts. Here, we report the first direct observation of a fully reversible en masse conformational transition between B- and A-DNA within live bacterial cells using Fourier transform infrared (FTIR) spectroscopy. This biospectroscopic technique allows for non-invasive and reagent-free examination of the holistic biochemistry of samples. For this reason, we have been able to observe the previously unknown conformational transition in all four species of bacteria investigated. Detection of this transition is evidence of a previously unexplored biological significance for A-DNA and highlights the need for new research into the role that A-DNA plays as a cellular defence mechanism and in stabilizing the DNA conformation. Such studies are pivotal in understanding the role of A-DNA in the evolutionary pathway of nucleic acids. Furthermore, this discovery demonstrates the exquisite capabilities of FTIR spectroscopy and opens the door for further investigations of cell biochemistry with this under-used technique.

DNA conformational flexibility study using phosphate backbone neutralization model

Due to the critical role of DNA in the processes of the cell cycle, the structural and physicochemical properties of DNA have long been of concern. In the present work, the effect of interplay between the DNA duplex and metal ions in solution on the DNA structure and conformational flexibility is studied by comparing the structure and dynamic conformational behavior of a duplex in a normal form and its “null isomer” using molecular dynamics methods. It was found that the phosphate neutralization changes the cation atmosphere around the DNA duplex greatly, increases the major groove width, decreases the minor groove width, and reduces the global bending direction preference. We also noted that the probability of BI phosphate linkages increases significantly because of the charge reduction in the backbone phosphate groups. More importantly, we found that the electrostatic effect on the DNA conformational flexibility is dependent on the sequence; that is, the phosphate backbone neutralization induces the global dynamic bending to be less flexible for GC-rich sequences but more flexible for AT-rich sequences.

Small-molecule binding to the DNA minor groove is mediated by a conserved water cluster

High-resolution crystal structures of the DNA duplex sequence d(CGCGAATTCGCG)(2) complexed with three minor-groove ligands are reported. A highly conserved cluster of 11 linked water molecules has been found in the native and all 3 ligand-bound structures, positioned at the boundary of the A/T and G/C regions where the minor groove widens. This cluster appears to play a key structural role in stabilizing noncovalently binding small molecules in the AT region of the B-DNA minor groove. The cluster extends from the backbone phosphate groups along the mouth of the groove and links to DNA and ligands by a network of hydrogen bonds that help to maintain the ligands in position. This arrangement of water molecules is distinct from, but linked by, hydrogen bonding to the well-established spine of hydration, which is displaced by bound ligands. Features of the water cluster and observed differences in binding modes help to explain the measured binding affinities and thermodynamic characteristics of these ligands on binding to AT sites in DNA.

Interactions of Mn2+ with a non-self-complementary Z-type DNA duplex

Crystal structures of the hexanucleotide d(CACGCG)·d(CGCGTG) were determined in two crystal lattices when different concentrations of the counterion Mn2+ were used in crystallization. The availability of Mn2+ during the crystallization process appears to play an important role in inducing different crystal packings that lead to crystals belonging to the two space groups P2(1) and P6(5). Analysis of the molecular interactions of Mn2+ with the Z-form duplexes shows direct coordination to the purine residues G and A.

Effect of initial ion positions on the interactions of monovalent and divalent ions with a DNA duplex as revealed with atomistic molecular dynamics simulations

Monovalent (Na(+)) and divalent (Mg(2+)) ion distributions around the Dickerson-Drew dodecamer were studied by atomistic molecular dynamics (MD) simulations with AMBER molecular modeling software. Different initial placements of ions were tried and the resulting effects on the ion distributions around DNA were investigated. For monovalent ions, results were found to be nearly independent of initial cation coordinates. However, Mg(2+) ions demonstrated a strong initial coordinate dependent behavior. While some divalent ions initially placed near the DNA formed essentially permanent direct coordination complexes with electronegative DNA atoms, Mg(2+) ions initially placed further away from the duplex formed a full, nonexchanging, octahedral first solvation shell. These fully solvated cations were still capable of binding with DNA with events lasting up to 20 ns, and in comparison were bound much longer than Na(+) ions. Force field parameters were also investigated with modest and little differences arising from ion (ions94 and ions08) and nucleic acid description (ff99, ff99bsc0, and ff10), respectively. Based on known Mg(2+) ion solvation structure, we conclude that in most cases Mg(2+) ions retain their first solvation shell, making only solvent-mediated contacts with DNA duplex. The proper way to simulate Mg(2+) ions around DNA duplex, therefore, should begin with ions placed in the bulk water.

Structure, Stiffness and Substates of the Dickerson-Drew Dodecamer

The Dickerson-Drew dodecamer (DD) d-[CGCGAATTCGCG]₂ is a prototypic B-DNA molecule whose sequence-specific structure and dynamics have been investigated by many experimental and computational studies. Here, we present an analysis of DD properties based on extensive atomistic molecular dynamics (MD) simulations using different ionic conditions and water models. The 0.6-2.4-µs-long MD trajectories are compared to modern crystallographic and NMR data. In the simulations, the duplex ends can adopt an alternative base-pairing, which influences the oligomer structure. A clear relationship between the BI/BII backbone substates and the basepair step conformation has been identified, extending previous findings and exposing an interesting structural polymorphism in the helix. For a given end pairing, distributions of the basepair step coordinates can be decomposed into Gaussian-like components associated with the BI/BII backbone states. The nonlocal stiffness matrices for a rigid-base mechanical model of DD are reported for the first time, suggesting salient stiffness features of the central A-tract. The Riemann distance and Kullback-Leibler divergence are used for stiffness matrix comparison. The basic structural parameters converge very well within 300 ns, convergence of the BI/BII populations and stiffness matrices is less sharp. Our work presents new findings about the DD structural dynamics, mechanical properties, and the coupling between basepair and backbone configurations, including their statistical reliability. The results may also be useful for optimizing future force fields for DNA.

Mg2+ in the major groove modulates B-DNA structure and dynamics

This study investigates the effect of Mg(2+) bound to the DNA major groove on DNA structure and dynamics. The analysis of a comprehensive dataset of B-DNA crystallographic structures shows that divalent cations are preferentially located in the DNA major groove where they interact with successive bases of (A/G)pG and the phosphate group of 5'-CpA or TpG. Based on this knowledge, molecular dynamics simulations were carried out on a DNA oligomer without or with Mg(2+) close to an ApG step. These simulations showed that the hydrated Mg(2+) forms a stable intra-strand cross-link between the two purines in solution. ApG generates an electrostatic potential in the major groove that is particularly attractive for cations; its intrinsic conformation is well-adapted to the formation of water-mediated hydrogen bonds with Mg(2+). The binding of Mg(2+) modulates the behavior of the 5'-neighboring step by increasing the BII (ε-ζ>0°) population of its phosphate group. Additional electrostatic interactions between the 5'-phosphate group and Mg(2+) strengthen both the DNA-cation binding and the BII character of the 5'-step. Cation binding in the major groove may therefore locally influence the DNA conformational landscape, suggesting a possible avenue for better understanding how strong DNA distortions can be stabilized in protein-DNA complexes.