Measuring inter-DNA potentials in solution

Effects of monovalent and divalent cations on the rheology of entangled DNA

In this paper we investigate the effects of varying cation valency and concentration on the rheology of entangled λDNA solutions. We show that monovalent cations moderately increase the viscoelasticty of the solutions mainly by stabilising linear concatenation of λDNA "monomers" via hybridisation of their sticky ends. On the contrary, divalent cations have a far more complex and dramatic effect on the rheology of the solution and we observe evidence of inter-molecular DNA-DNA bridging by Mg2+. We argue that these results may be interesting in the context of dense solutions of single and double stranded DNA, e.g. in vivo or in biotechnology applications such as DNA origami and DNA hydrogels.

Programming the morphology of DNA origami crystals by magnesium ion strength

Harnessing the programmable nature of DNA origami for controlling structural features in crystalline materials affords opportunities to bring crystal engineering to a remarkable level. However, the challenge of crystallizing a single type of DNA origami unit into varied structural outcomes remains, given the requirement for specific DNA designs for each targeted structure. Here, we show that crystals with distinct equilibrium phases and shapes can be realized using a single DNA origami morphology with an allosteric factor to modulate the binding coordination. As a result, origami crystals undergo phase transitions from a simple cubic lattice to a simple hexagonal (SH) lattice and eventually to a face-centered cubic (FCC) lattice. After selectively removing internal nanoparticles from DNA origami building blocks, the body-centered tetragonal and chalcopyrite lattice are derived from the SH and FCC lattices, respectively, revealing another phase transition involving crystal system conversions. The rich phase space was realized through the de novo synthesis of crystals under varying solution environments, followed by the individual characterizations of the resulting products. Such phase transitions can lead to associated transitions in the shape of the resulting products. Hexagonal prism crystals, crystals characterized by triangular facets, and twinned crystals are observed to form from SH and FCC systems, which have not previously been experimentally realized by DNA origami crystallization. These findings open a promising pathway toward accessing a rich phase space with a single type of building block and wielding other instructions as tools to develop crystalline materials with tunable properties.

In situ small-angle X-ray scattering reveals strong condensation of DNA origami during silicification

Silicification of DNA origami structures increases their stability and provides chemical protection. Yet, it is unclear whether the whole DNA framework is embedded or if silica just forms an outer shell and how silicification affects the origami's internal structure. Employing in situ small-angle X-ray scattering (SAXS), we show that addition of silica precursors induces substantial condensation of the DNA origami at early reaction times by almost 10 %. Subsequently, the overall size of the silicified DNA origami increases again due to increasing silica deposition. We further identify the SAXS Porod invariant as a reliable, model-free parameter for the evaluation of the amount of silica formation at a given time. Contrast matching of the DNA double helix Lorentzian peak reveals silica growth also inside the origami. The less polar silica forming within the origami structure, replacing more than 40 % of the internal hydration water, causes a hydrophobic effect: condensation. DNA origami objects with flat surfaces show a strong tendency towards aggregation during silicification, presumably driven by the same entropic forces causing condensation. Maximally condensed origami displayed thermal stability up to 60 °C. Our studies provide insights into the silicification reaction allowing for the formulation of optimized reaction protocols.

Opto-Electrostatic Determination of Nucleic Acid Double-Helix Dimensions and the Structure of the Molecule–Solvent Interface

A DNA molecule is highly electrically charged in solution. The electrical potential at the molecular surface is known to vary strongly with the local geometry of the double helix and plays a pivotal role in DNA–protein interactions. Further out from the molecular surface, the electrical field propagating into the surrounding electrolyte bears fingerprints of the three-dimensional arrangement of the charged atoms in the molecule. However, precise extraction of the structural information encoded in the electrostatic “far field” has remained experimentally challenging. Here, we report an optical microscopy-based approach that detects the field distribution surrounding a charged molecule in solution, revealing geometric features such as the radius and the average rise per basepair of the double helix with up to sub-Angstrom precision, comparable with traditional molecular structure determination techniques like X-ray crystallography and nuclear magnetic resonance. Moreover, measurement of the helical radius furnishes an unprecedented view of both hydration and the arrangement of cations at the molecule–solvent interface. We demonstrate that a probe in the electrostatic far field delivers structural and chemical information on macromolecules, opening up a new dimension in the study of charged molecules and interfaces in solution.

The Effects of Flexibility on dsDNA–dsDNA Interactions

A detailed understanding of the physical mechanism of ion-mediated dsDNA interactions is important in biological functions such as DNA packaging and homologous pairing. We report the potential of mean force (PMF) or the effective solvent mediated interactions between two parallel identical dsDNAs as a function of interhelical separation in 0.15 M NaCl solution. Here, we study the influence of flexibility of dsDNAs on the effective interactions by comparing PMFs between rigid models and flexible ones. The role of flexibility of dsDNA pairs in their association is elucidated by studying the energetic properties of Na⁺ ions as well as the fluctuations of ions around dsDNAs. The introduction of flexibility of dsDNAs softens the vdW contact wall and induces more counterion fluctuations around dsDNAs. In addition, flexibility facilitates the Na⁺ ions dynamics affecting their distribution. The results quantify the extent of attraction influenced by dsDNA flexibility and further emphasize the importance of non-continuum solvation approaches.

Electrostatic free energies carry structural information on nucleic acid molecules in solution

Over the last several decades, a range of experimental techniques from x-ray crystallography and atomic force microscopy to nuclear magnetic resonance and small angle x-ray scattering have probed nucleic acid structure and conformation with high resolution both in the condensed state and in solution. We present a computational study that examines the prospect of using electrostatic free energy measurements to detect 3D conformational properties of nucleic acid molecules in solution. As an example, we consider the conformational difference between A- and B-form double helices whose structures differ in the values of two key parameters—the helical radius and rise per basepair. Mapping the double helix onto a smooth charged cylinder reveals that electrostatic free energies for molecular helices can, indeed, be described by two parameters: the axial charge spacing and the radius of a corresponding equivalent cylinder. We show that electrostatic free energies are also sensitive to the local structure of the molecular interface with the surrounding electrolyte. A free energy measurement accuracy of 1%, achievable using the escape time electrometry (ET e) technique, could be expected to offer a measurement precision on the radius of the double helix of approximately 1 Å. Electrostatic free energy measurements may, therefore, not only provide information on the structure and conformation of biomolecules but could also shed light on the interfacial hydration layer and the size and arrangement of counterions at the molecular interface in solution.

Divalent metal ions and intermolecular interactions facilitate DNA network formation

The interactions between divalent metal ions and DNA are crucial for basic life processes. These interactions are also important in advanced technological products such as DNA-based ion sensors. Current polyelectrolyte theories cannot describe these interactions well and do not consider the corresponding dynamics. In this study, we report the single-molecule dynamics of the binding of divalent metal ions to a single DNA molecule and the morphology characterization of the complex. We found that most of the divalent metal ions (Mn²⁺, Zn²⁺, Co²⁺, Ni²⁺, and Cd²⁺), except Mg²⁺ and Ca²⁺, could cause monomolecular DNA condensation. For transition metal ions, different ionic strengths were required to induce the compaction, and different shortening speeds were displayed in the dynamics, indicating ionic specificity. Atomic force microscopy revealed that the morphologies of the metal ion-DNA complexes were affected by the ionic strength of the metal ion, DNA chain length, and DNA concentration. At low metal ion concentration, DNA tended to adopt a random coil conformation. Increasing the ionic strength led to network-like condensed structures, suggesting that divalent metal ions can induce attraction between DNA molecules. Furthermore, higher DNA concentration and longer chain length enhanced intermolecular interactions and consequently resulted in network structures with a higher degree of interconnectivity.

Visualizing Disordered Single-Stranded RNA: Connecting Sequence, Structure, and Electrostatics

Disordered homopolymeric regions of single-stranded RNA, such as U or A tracts, are found within functional RNAs where they play distinct roles in defining molecular structure and facilitating recognition by partners. Despite this prominence, details of conformational and biophysical properties of these regions have not yet been resolved. We apply a number of experimental techniques to investigate the conformations of these biologically important motifs and provide quantitative measurements of their ion atmospheres. Single strands of RNA display pronounced sequence-dependent conformations that relate to the unique ion atmospheres each attracts. Chains of rU bases are relatively unstructured under all conditions, while chains of rA bases display distinct ordering through stacking or clustering motifs, depending on the composition of the surrounding solution. These dramatic structural differences are consistent with the measured disparity in ion composition and atmospheres around each homopolymer, revealing a complex interplay of base, ion, and single-strand ordering. The unique structural and ionic signatures of homopolymer ssRNAs explains their role(s) in folding structured RNAs and may explain their distinct recognition by protein partners.

Counterion-Dependent Mechanisms of DNA Origami Nanostructure Stabilization Revealed by Atomistic Molecular Simulation

The DNA origami technique has proven to have tremendous potential for therapeutic and diagnostic applications like drug delivery, but the relatively low concentrations of cations in physiological fluids cause destabilization and degradation of DNA origami constructs preventing in vivo applications. To reveal the mechanisms behind DNA origami stabilization by cations, we performed atomistic molecular dynamics simulations of a DNA origami rectangle in aqueous solvent with varying concentrations of magnesium and sodium as well as polyamines like oligolysine and spermine. We explored the binding of these ions to DNA origami in detail and found that the mechanism of stabilization differs between ion types considerably. While sodium binds weakly and quickly exchanges with the solvent, magnesium and spermine bind close to the origami with spermine also located in between helices, stabilizing the crossovers characteristic for DNA origami and reducing repulsion of parallel helices. In contrast, oligolysine of length ten prevents helix repulsion by binding to adjacent helices with its flexible side chains, spanning the gap between the helices. Shorter oligolysine molecules with four subunits are weak stabilizers as they lack both the ability to connect helices and to prevent helix repulsion. This work thus shows how the binding modes of ions influence the stabilization of DNA origami nanostructures on a molecular level.

Resolving solution conformations of the model semi-flexible polyelectrolyte homogalacturonan using molecular dynamics simulations and small-angle x-ray scattering

The conformation of polyelectrolytes in the solution state has long been of interest in polymer science. Herein we utilize all atom molecular dynamics simulations (MD) and small-angle x-ray scattering experiments (SAXS) to elucidate the molecular structure of the model polyelectrolyte homogalacturonan. Several degrees of polymerization were studied and in addition partial methylesterification of the otherwise charge-carrying carboxyl groups was used in order to generate samples with varying intra-chain charge distributions. It is shown that at length scales above around 1nm the conformation of isolated chains has surprisingly little dependence on the charge distribution or the concentration of attendant monovalent salts, reflective of the intrinsic stiffness of the saccharide rings and the dynamical constraints of the glycosidic linkage. Indeed the conformation of isolated chains over all accessible length scales is well described by the atomic coordinates available from fibre diffraction studies. Furthermore, in more concentrated systems it is shown that, after careful analysis of the SAXS data, the form of the inter-particle effects heralded by the emergence of a so-called polyelectrolyte peak, can be extracted, and that this phenomena can be reproduced by multiple chain MD simulations.

"End-to-end" stacking of small dsRNA

PELDOR (pulsed electron-electron double resonance) is an established method to study intramolecular distances and can give evidence for conformational changes and flexibilities. However, it can also be used to study intermolecular interactions as for example oligerimization. Here, we used PELDOR to study the "end-to-end" stacking of small double-stranded (ds) RNAs. For this study, the dsRNA molecules were only singly labeled with the spin label TPA to avoid multispin effects and to measure only the intermolecular stacking interactions. It can be shown that small dsRNAs tend to assemble to rod-like structures due to π-π interactions between the base pairs at the end of the strands. On the one hand, these interactions can influence or complicate measurements aimed at the determining of the structure and dynamics of the dsRNA molecule itself. On the other hand, it can be interesting to study such intermolecular stacking interactions in more detail, as for example their dependence on ion concentration. We quantitatively determined the stacking probability as a function of the monovalent NaCl salt and the dsRNA concentration. From these data, the dissociation constant K_d was deduced and found to depend on the ratio between the NaCl salt and dsRNA concentrations. Additionally, the distances and distance distributions obtained predict a model for the stacking geometry of dsRNAs. Introducing a nucleotide overhangs at one end of the dsRNA molecule restricts the stacking to the other end, leading only to dimer formations. Introducing such an overhang at both ends of the dsRNA molecule fully suppresses stacking, as we demonstrate by PELDOR experiments quantitatively.

Direct Demonstration of DNA Compaction Mediated by Divalent Counterions

We unambiguously demonstrated DNA attraction and its regulation mediated by divalent cations Mg²⁺ and Ca²⁺ by tethering a DNA single chain at various pH solutions. It is found that DNA is compacted when the pH of the solution containing these divalent counterions is decreased below 5. When the pH of the medium is ∼4, DNA is in an unstable transition state, being able to switch between compact and extensible states. We can also regulate the DNA attraction through a cyclic process of DNA compaction and unraveling by alternating the pH of the solution between 3 and 8. The corresponding change of morphology of DNA modulated by pH is also confirmed by atomic force microscopy (AFM). In the theoretical aspect, the present experimental finding is consistent with the coarse-grained simulation of Langevin dynamics on the effect of pH on DNA in a solution of divalent counterions.

A General Approach to Access Morphologies of Polyoxometalates in Solution by Using SAXS: An Ab Initio Modeling Protocol

Herein, we reported a general protocol for an ab initio modeling approach to deduce structure information of polyoxometalates (POMs) in solutions from scattering data collected by the small-angle X-ray scattering (SAXS) technique. To validate the protocol, the morphologies of a serious of known POMs in either aqueous or organic solvents were analyzed. The obtained particle morphologies were compared and confirmed with previous reported crystal structures. To extend the feasibility of the protocol to an unknown system of aqueous solutions of Na₂ MoO₄ with the pH ranging from -1 to 8.35, the formation of {Mo₃₆ } clusters was probed, identified, and confirmed by SAXS. The approach was further optimized with a multi-processing capability to achieve fast analysis of experimental data, thereby, facilitating in situ studies of formations of POMs in solutions. The advantage of this approach is to generate intuitive 3D models of POMs in solutions without confining information such as symmetries and possible sizes.

Potential of mean force between oppositely charged nanoparticles: A comprehensive comparison between Poisson-Boltzmann theory and Monte Carlo simulations

Ion-mediated interactions between like-charged polyelectrolytes have been paid much attention, and the Poisson-Boltzmann (PB) theory has been shown to fail in qualitatively predicting multivalent ion-mediated like-charge attraction. However, inadequate attention has been paid to the ion-mediated interactions between oppositely charged polyelectrolytes. In this work, the potentials of mean force (PMF) between oppositely charged nanoparticles in 1:1 and 2:2 salt solutions were investigated by Monte Carlo simulations and the PB theory. Our calculations show that the PMFs between oppositely charged nanoparticles are generally attractive in 1:1 and 2:2 salt solutions and that such attractive PMFs become weaker at higher 1:1 or 2:2 salt concentrations. The comprehensive comparisons show that the PB theory can quantitatively predict the PMFs between oppositely charged nanoparticles in 1:1 salt solutions, except for the slight deviation at very high 1:1 salt concentration. However, for 2:2 salt solutions, the PB theory generally overestimates the attractive PMF between oppositely charged nanoparticles, and this overestimation becomes more pronounced for nanoparticles with higher charge density and for higher 2:2 salt concentration. Our microscopic analyses suggest that the overestimation of the PB theory on the attractive PMFs for 2:2 salt solutions is attributed to the underestimation of divalent ions bound to nanoparticles.

Like-charge attraction and opposite-charge decomplexation between polymers and DNA molecules

We scrutinize the effect of polyvalent ions on polymer-DNA interactions. We extend a recently developed test-charge theory [S. Buyukdagli et al., Phys. Rev. E 94, 042502 (2016)1539-375510.1103/PhysRevE.94.042502] to the case of a stiff polymer interacting with a DNA molecule in an electrolyte mixture. The theory accounts for one-loop level electrostatic correlation effects such as the ionic cloud deformation around the strongly charged DNA molecule as well as image-charge forces induced by the low DNA permittivity. Our model can reproduce and explain various characteristics of the experimental phase diagrams for polymer solutions. First, the addition of polyvalent cations to the electrolyte solution results in the attraction of the negatively charged polymer by the DNA molecule. The glue of the like-charge attraction is the enhanced shielding of the polymer charges by the dense counterion layer at the DNA surface. Second, through the shielding of the DNA-induced electrostatic potential, mono- and polyvalent cations of large concentration both suppress the like-charge attraction. Within the same formalism, we also predict a new opposite-charge repulsion effect between the DNA molecule and a positively charged polymer. In the presence of polyvalent anions such as sulfate or phosphate, their repulsion by the DNA charges leads to the charge screening deficiency of the region around the DNA molecule. This translates into a repulsive force that results in the decomplexation of the polymer from DNA. This opposite-charge repulsion phenomenon can be verified by current experiments and the underlying mechanism can be beneficial to gene therapeutic applications where the control over polymer-DNA interactions is the key factor.

Electrolyte-Mediated Assembly of Charged Nanoparticles

Solutions at high salt concentrations are used to crystallize or segregate charged colloids, including proteins and polyelectrolytes via a complex mechanism referred to as "salting-out". Here, we combine small-angle X-ray scattering (SAXS), molecular dynamics (MD) simulations, and liquid-state theory to show that salting-out is a long-range interaction, which is controlled by electrolyte concentration and colloid charge density. As a model system, we analyze Au nanoparticles coated with noncomplementary DNA designed to prevent interparticle assembly via Watson-Crick hybridization. SAXS shows that these highly charged nanoparticles undergo "gas" to face-centered cubic (FCC) to "glass-like" transitions with increasing NaCl or CaCl2 concentration. MD simulations reveal that the crystallization is concomitant with interparticle interactions changing from purely repulsive to a "long-range potential well" condition. Liquid-state theory explains this attraction as a sum of cohesive and depletion forces that originate from the interelectrolyte ion and electrolyte-ion-nanoparticle positional correlations. Our work provides fundamental insights into the effect of ionic correlations in the salting-out mechanism and suggests new routes for the crystallization of colloids and proteins using concentrated salts.

Salt contribution to the flexibility of single-stranded nucleic acid offinite length

Nucleic acids are negatively charged macromolecules and their structure properties are strongly coupled to metal ions in solutions. In this article, the salt effects on the flexibility of single-stranded (ss) nucleic acid chain ranging from 12 to 120 nucleotides are investigated systematically by the coarse-grained Monte Carlo simulations where the salt ions are considered explicitly and the ss chain is modeled with the virtual-bond structural model. Our calculations show that, the increase of ion concentration causes the structural collapse of ss chain and multivalent ions are much more efficient in causing such collapse, and both trivalent/small divalent ions can induce more compact state than a random relaxation state. We found that monovalent, divalent, and trivalent ions can all overcharge ss chain, and the dominating source for such overcharging changes from ion-exclusion-volume effect to ion Coulomb correlations. In addition, the predicted Na(+) and Mg(2+)-dependent persistence length l(p)'s of ss nucleic acid are in accordance with the available experimental data, and through systematic calculations, we obtained the empirical formulas for l(p) as a function of [Na(+)], [Mg(2+)] and chain length.

Polyelectrolyte properties of single stranded DNA measured using SAXS and single-molecule FRET: Beyond the wormlike chain model

Nucleic acids are highly charged polyelectrolytes that interact strongly with salt ions. Rigid, base-paired regions are successfully described with wormlike chain models, but nonbase-paired single stranded regions have fundamentally different polymer properties because of their greater flexibility. Recently, attention has turned to single stranded nucleic acids due to the growing recognition of their biological importance, as well as the availability of sophisticated experimental techniques sensitive to the conformation of individual molecules. We investigate polyelectrolyte properties of poly(dT), an important and widely studied model system for flexible single stranded nucleic acids, in physiologically important mixed mono- and divalent salt. We report measurements of the form factor and interparticle interactions using SAXS, end-to-end distances using smFRET, and number of excess ions using ASAXS. We present a coarse-grained model that accounts for flexibility, excluded volume, and electrostatic interactions in these systems. Predictions of the model are validated against experiment. We also discuss the state of all-atom, explicit solvent molecular dynamics simulations of poly(dT), the next step in understanding the complexities of ion interactions with these highly charged and flexible polymers.

Determining the Locations of Ions and Water around DNA from X-Ray Scattering Measurements

Nucleic acids carry a negative charge, attracting salt ions and water. Interactions with these components of the solvent drive DNA to condense, RNA to fold, and proteins to bind. To understand these biological processes, knowledge of solvent structure around the nucleic acids is critical. Yet, because they are often disordered, ions and water evade detection by x-ray crystallography and other high-resolution methods. Small-angle x-ray scattering (SAXS) is uniquely sensitive to the spatial correlations between solutes and the surrounding solvent. Thus, SAXS provides an experimental constraint to guide or test emerging solvation theories. However, the interpretation of SAXS profiles is nontrivial because of the difficulty in separating the scattering signals of each component: the macromolecule, ions, and hydration water. Here, we demonstrate methods for robustly deconvoluting these signals, facilitating a more straightforward comparison with theory. Using SAXS data collected on an absolute intensity scale for short DNA duplexes in solution with Na(+), K(+), Rb(+), or Cs(+) counterions, we mathematically decompose the scattering profiles into components (DNA, water, and ions) and validate the decomposition using anomalous scattering measurements. In addition, we generate a library of physically motivated ion atmosphere models and rank them by agreement with the scattering data. The best-fit models have relatively compact ion atmospheres when compared to predictions from the mean-field Poisson-Boltzmann theory of electrostatics. Thus, the x-ray scattering methods presented here provide a valuable measurement of the global structure of the ion atmosphere that can be used to test electrostatics theories that go beyond the mean-field approximation.

Potential of mean force between like-charged nanoparticles: Many-body effect

Ion-mediated interaction is important for the properties of polyelectrolytes such as colloids and nucleic acids. The effective pair interactions between two polyelectrolytes have been investigated extensively, but the many-body effect for multiple polyelectrolytes still remains elusive. In this work, the many-body effect in potential of mean force (PMF) between like-charged nanoparticles in various salt solutions has been comprehensively examined by Monte Carlo simulation and the nonlinear Poisson-Boltzmann theory. Our calculations show that, at high 1:1 salt, the PMF is weakly repulsive and appears additive, while at low 1:1 salt, the additive assumption overestimates the repulsive many-body PMF. At low 2:2 salt, the pair PMF appears weakly repulsive while the many-body PMF can become attractive. In contrast, at high 2:2 salt, the pair PMF is apparently attractive while the many-body effect can cause a weaker attractive PMF than that from the additive assumption. Our microscopic analyses suggest that the elusive many-body effect is attributed to ion-binding which is sensitive to ion concentration, ion valence, number of nanoparticles and charges on nanoparticles.

DNA under Force: Mechanics, Electrostatics, and Hydration

Quantifying the basic intra- and inter-molecular forces of DNA has helped us to better understand and further predict the behavior of DNA. Single molecule technique elucidates the mechanics of DNA under applied external forces, sometimes under extreme forces. On the other hand, ensemble studies of DNA molecular force allow us to extend our understanding of DNA molecules under other forces such as electrostatic and hydration forces. Using a variety of techniques, we can have a comprehensive understanding of DNA molecular forces, which is crucial in unraveling the complex DNA functions in living cells as well as in designing a system that utilizes the unique properties of DNA in nanotechnology.

Ion competition in condensed DNA arrays in the attractive regime

Physical origin of DNA condensation by multivalent cations remains unsettled. Here, we report quantitative studies of how one DNA-condensing ion (Cobalt(3+) Hexammine, or Co(3+)Hex) and one nonDNA-condensing ion (Mg(2+)) compete within the interstitial space in spontaneously condensed DNA arrays. As the ion concentrations in the bath solution are systematically varied, the ion contents and DNA-DNA spacings of the DNA arrays are determined by atomic emission spectroscopy and x-ray diffraction, respectively. To gain quantitative insights, we first compare the experimentally determined ion contents with predictions from exact numerical calculations based on nonlinear Poisson-Boltzmann equations. Such calculations are shown to significantly underestimate the number of Co(3+)Hex ions, consistent with the deficiencies of nonlinear Poisson-Boltzmann approaches in describing multivalent cations. Upon increasing the concentration of Mg(2+), the Co(3+)Hex-condensed DNA array expands and eventually redissolves as a result of ion competition weakening DNA-DNA attraction. Although the DNA-DNA spacing depends on both Mg(2+) and Co(3+)Hex concentrations in the bath solution, it is observed that the spacing is largely determined by a single parameter of the DNA array, the fraction of DNA charges neutralized by Co(3+)Hex. It is also observed that only ∼20% DNA charge neutralization by Co(3+)Hex is necessary for spontaneous DNA condensation. We then show that the bath ion conditions can be reduced to one variable with a simplistic ion binding model, which is able to describe the variations of both ion contents and DNA-DNA spacings reasonably well. Finally, we discuss the implications on the nature of interstitial ions and cation-mediated DNA-DNA interactions.

Elucidating internucleosome interactions and the roles of histone tails

The nucleosome is the first level of genome organization and regulation in eukaryotes where negatively charged DNA is wrapped around largely positively charged histone proteins. Interaction between nucleosomes is dominated by electrostatics at long range and guided by specific contacts at short range, particularly involving their flexible histone tails. We have thus quantified how internucleosome interactions are modulated by salts (KCl, MgCl2) and histone tail deletions (H3, H4 N-terminal), using small-angle x-ray scattering and theoretical modeling. We found that measured effective charges at low salts are ∼1/5th of the theoretically predicted renormalized charges and that H4 tail deletion suppresses the attraction at high salts to a larger extent than H3 tail deletion.

Helical structure determines different susceptibilities of dsDNA, dsRNA, and tsDNA to counterion-induced condensation

Recent studies of counterion-induced condensation of nucleic acid helices into aggregates produced several puzzling observations. For instance, trivalent cobalt hexamine ions condensed double-stranded (ds) DNA oligomers but not their more highly charged dsRNA counterparts. Divalent alkaline earth metal ions condensed triple-stranded (ts) DNA oligomers but not dsDNA. Here we show that these counterintuitive experimental results can be rationalized within the electrostatic zipper model of interactions between molecules with helical charge motifs. We report statistical mechanical calculations that reveal dramatic and nontrivial interplay between the effects of helical structure and thermal fluctuations on electrostatic interaction between oligomeric nucleic acids. Combining predictions for oligomeric and much longer helices, we also interpret recent experimental studies of the role of counterion charge, structure, and chemistry. We argue that an electrostatic zipper attraction might be a major or even dominant force in nucleic acid condensation.

DNA methylation regulated nucleosome dynamics

A strong correlation between nucleosome positioning and DNA methylation patterns has been reported in literature. However, the mechanistic model accounting for the correlation remains elusive. In this study, we evaluated the effects of specific DNA methylation patterns on modulating nucleosome conformation and stability using FRET and SAXS. CpG dinucleotide repeats at 10 bp intervals were found to play different roles in nucleosome stability dependent on their methylation states and their relative nucleosomal locations. An additional (CpG)₅ stretch located in the nucleosomal central dyad does not alter the nucleosome conformation, but significant conformational differences were observed between the unmethylated and methylated nucleosomes. These findings suggest that the correlation between nucleosome positioning and DNA methylation patterns can arise from the variations in nucleosome stability dependent on their sequence and epigenetic content. This knowledge will help to reveal the detailed role of DNA methylation in regulating chromatin packaging and gene transcription.

Fast-SAXS-pro: a unified approach to computing SAXS profiles of DNA, RNA, protein, and their complexes

A generalized method, termed Fast-SAXS-pro, for computing small angle x-ray scattering (SAXS) profiles of proteins, nucleic acids, and their complexes is presented. First, effective coarse-grained structure factors of DNA nucleotides are derived using a simplified two-particle-per-nucleotide representation. Second, SAXS data of a 18-bp double-stranded DNA are measured and used for the calibration of the scattering contribution from excess electron density in the DNA solvation layer. Additional test on a 25-bp DNA duplex validates this SAXS computational method and suggests that DNA has a different contribution from its hydration surface to the total scattering compared to RNA and protein. To account for such a difference, a sigmoidal function is implemented for the treatment of non-uniform electron density across the surface of a protein/nucleic-acid complex. This treatment allows differential scattering from the solvation layer surrounding protein/nucleic-acid complexes. Finally, the applications of this Fast-SAXS-pro method are demonstrated for protein/DNA and protein/RNA complexes.

DNA self-assembly: from chirality to evolution

Transient or long-term DNA self-assembly participates in essential genetic functions. The present review focuses on tight DNA-DNA interactions that have recently been found to play important roles in both controlling DNA higher-order structures and their topology. Due to their chirality, double helices are tightly packed into stable right-handed crossovers. Simple packing rules that are imposed by DNA geometry and sequence dictate the overall architecture of higher order DNA structures. Close DNA-DNA interactions also provide the missing link between local interactions and DNA topology, thus explaining how type II DNA topoisomerases may sense locally the global topology. Finally this paper proposes that through its influence on DNA self-assembled structures, DNA chirality played a critical role during the early steps of evolution.

Salt dependence of the radius of gyration and flexibility of single-stranded DNA in solution probed by small-angle x-ray scattering

Short single-stranded nucleic acids are ubiquitous in biological processes; understanding their physical properties provides insights to nucleic acid folding and dynamics. We used small-angle x-ray scattering to study 8-100 residue homopolymeric single-stranded DNAs in solution, without external forces or labeling probes. Poly-T's structural ensemble changes with increasing ionic strength in a manner consistent with a polyelectrolyte persistence length theory that accounts for molecular flexibility. For any number of residues, poly-A is consistently more elongated than poly-T, likely due to the tendency of A residues to form stronger base-stacking interactions than T residues.

Distribution of counterions and interaction between two similarly charged dielectric slabs: roles of charge discreteness and dielectric inhomogeneity

The distribution of counterions and the electrostatic interaction between two similarly charged dielectric slabs is studied in the strong coupling limit. Dielectric inhomogeneities and discreteness of charge on the slabs have been taken into account. It is found that the amount of dielectric constant difference between the slabs and the environment, and the discreteness of charge on the slabs have opposing effects on the equilibrium distribution of the counterions. At small interslab separations, increasing the amount of dielectric constant difference increases the tendency of the counterions toward the middle of the intersurface space between the slabs and the discreteness of charge pushes them to the surfaces of the slabs. In the limit of point charges, independent of the strength of dielectric inhomogeneity, counterions distribute near the surfaces of the slabs. The interaction between the slabs is attractive at low temperatures and its strength increases with the dielectric constant difference. At room temperature, the slabs may completely attract each other, reach to an equilibrium separation, or have two equilibrium separations with a barrier in between, depending on the system parameters.

Understanding thermodynamic competitivity between biopolymer folding and misfolding under large-scale intermolecular interactions

Cooperativity is a hallmark of spontaneous biopolymer folding. The presence of intermolecular interactions could create off-pathway misfolding structures and suppress folding cooperativity. This raises the hypothesis that thermodynamic competitivity between off-pathway misfolding and on-pathway folding may intervene with cooperativity and govern biopolymer folding dynamics under conditions permitting large-scale intermolecular interactions. Here we report direct imaging and theoretical modeling of thermodynamic competitivity between biopolymer folding and misfolding under such conditions, using a two-dimensional array of proton-fueled DNA molecular motors packed at the maximal density as a model system. Time-resolved liquid-phase atomic force microscopy with enhanced phase contrast revealed that the misfolding and folding intermediates transiently self-organize into spatiotemporal patterns on the nanoscale in thermodynamic states far away from equilibrium as a result of thermodynamic competitivity. Computer simulations using a novel cellular-automaton network model provide quantitative insights into how large-scale intermolecular interactions correlate the structural dynamics of individual biomolecules together at the systems level.

Counterion condensation theory of attraction between like charges in the absence of multivalent counterions

There is abundant experimental evidence suggesting the existence of attractive interactions among identically charged polyelectrolytes in ordinary salt solutions. The presence of multivalent counterions is not required. We review the relevant literature in detail and conclude that it merits more attention than it has received. We discuss also some recent observations of a low ionic strength attraction of negatively charged DNA to the region of a negatively charged glass nanoslit where the floor of the nanoslit meets the walls, again in the absence of multivalent ions. On the theoretical side, it has become clear that purely electrostatic interactions require the presence of multivalent counterions if they are to generate like-charge attraction. Any theory of like-charge attraction in the absence of multivalent counterions must therefore contain a non-electrostatic component. We point out that counterion condensation theory, which has predicted like-charge polyelectrolyte attraction in an intermediate range of distances in ordinary 1:1 salt conditions, contains both electrostatic and non-electrostatic elements. The non-electrostatic component of the theory is the modeling constraint that the counterions fall into two explicit populations, condensed and uncondensed. As reviewed in the paper, this physically motivated constraint is supported by strong experimental evidence. We proceed to offer an explanation of the nanoslit observations by showing in an idealized model that the line of intersection of two intersecting planes is a virtual polyelectrolyte. Since we have previously developed a counterion condensation theory of attraction of two like-charged polyelectrolytes, our suggestion is that the DNA is attracted to the virtual polyelectrolytes that may be located in the nanoslit where floor meets walls. We present the detailed calculations needed to document this suggestion: an extension of previous theory to the case of polyelectrolytes with like but not identical charges; the demonstration of counterion condensation on a plane with bare charge density greater than an explicitly exhibited critical value; a calculation of the free energy of the plane; a calculation of the interaction of a line charge polyelectrolyte with a like-charged plane; and the detailed demonstration that the line of intersection of two planes is a virtual polyelectrolyte.

Effects of a protecting osmolyte on the ion atmosphere surrounding DNA duplexes

Osmolytes are small, chemically diverse, organic solutes that function as an essential component of cellular stress response. Protecting osmolytes enhance protein stability via preferential exclusion, and nonprotecting osmolytes, such as urea, destabilize protein structures. Although much is known about osmolyte effects on proteins, less is understood about osmolyte effects on nucleic acids and their counterion atmospheres. Nonprotecting osmolytes destabilize nucleic acid structures, but effects of protecting osmolytes depend on numerous factors including the type of nucleic acid and the complexity of the functional fold. To begin quantifying protecting osmolyte effects on nucleic acid interactions, we used small-angle X-ray scattering (SAXS) techniques to monitor DNA duplexes in the presence of sucrose. This protecting osmolyte is a commonly used contrast matching agent in SAXS studies of protein-nucleic acid complexes; thus, it is important to characterize interaction changes induced by sucrose. Measurements of interactions between duplexes showed no dependence on the presence of up to 30% sucrose, except under high Mg(2+) conditions where stacking interactions were disfavored. The number of excess ions associated with DNA duplexes, reported by anomalous small-angle X-ray scattering (ASAXS) experiments, was sucrose independent. Although protecting osmolytes can destabilize secondary structures, our results suggest that ion atmospheres of individual duplexes remain unperturbed by sucrose.

SAXS studies of ion-nucleic acid interactions

Positively charged ions, atoms, or molecules compensate the high negative charge of the nucleic acid backbone. Their presence is critical to the biological function of DNA and RNA. This review focuses on experimental studies probing (a) interactions between small ions and nucleic acids and (b) ion-mediated interactions between nucleic acid duplexes. Experimental results on these simple model systems can be compared with specific theoretical models to validate their predictions. Small angle X-ray scattering (SAXS) provides unique insight into these interactions. Anomalous SAXS reports the spatial correlations of condensed (e.g., locally concentrated) counterions to individual DNA or RNA duplexes. SAXS very effectively reports interactions between nucleic acid helices, which range from strongly repulsive to strongly attractive depending on the ionic species present. The sign and strength of interparticle interactions are easily deduced from dramatic changes in the scattering profiles of interacting duplexes.

Theoretical assessment of the oligolysine model for ionic interactions in protein-DNA complexes

The observed salt dependence of charged ligand binding to polyelectrolytes, such as proteins to DNA or antithrombin to heparin, is often interpreted by means of the "oligolysine model." We review this model as derived entirely within the framework of the counterion condensation theory of polyelectrolytes with no introduction of outside assumptions. We update its comparison with experimental data on the structurally simple systems for which it was originally intended. We then compute the salt dependence of the binding free energy for a variety of protein-DNA complexes with nonlinear Poisson-Boltzmann (NLPB) simulation methods. The results of the NLPB calculations confirm the central prediction of the oligolysine model that the net charge density of DNA remains invariant to protein binding. Specifically, when a cationic protein residue penetrates the layer of counterions condensed on DNA, a counterion is released to bulk solution, and when an anionic protein residue penetrates the condensed counterion layer, an additional counterion is condensed from bulk solution. We also conclude, however, that the cumulative effect of charged protein residues distant from the binding interface makes a significant contribution to the salt dependence of binding, an observation not accommodated by the oligolysine model.

Double-stranded RNA resists condensation

Much attention has been focused on DNA condensation because of its fundamental biological importance. The recent discovery of new roles for RNA duplexes demands efficient packaging of double-stranded RNA for therapeutics. Here we report measurements of short DNA and RNA duplexes in the presence of trivalent ions. Under conditions where UV spectroscopy indicates condensation of DNA duplexes into (insoluble) precipitates, RNA duplexes remain soluble. Small angle x-ray scattering results suggest that the differing surface topologies of RNA and DNA may be crucial in generating the attractive forces that result in precipitation.

DNA-DNA interactions in tight supercoils are described by a small effective charge density

DNA-DNA interactions are important for genome compaction and transcription regulation. In studies of such complex processes, DNA is often modeled as a homogeneously charged cylinder and its electrostatic interactions are calculated within the framework of the Poisson-Boltzmann equation. Commonly, a charge adaptation factor is used to address limitations of this theoretical approach. Despite considerable theoretical and experimental efforts, a rigorous quantitative assessment of this parameter is lacking. Here, we comprehensively characterized DNA-DNA interactions in the presence of monovalent ions by analyzing the supercoiling behavior of single DNA molecules held under constant tension. Both a theoretical model and coarse-grained simulations of this process revealed a surprisingly small effective DNA charge of 40% of the nominal charge density, which was additionally supported by all-atom molecular dynamics simulations.

Electrostatics of strongly charged biological polymers: ion-mediated interactions and self-organization in nucleic acids and proteins

Charges on biological polymers in physiologically relevant solution conditions are strongly screened by water and salt solutions containing counter-ions. However, the entropy of these counterions can result in surprisingly strong interactions between charged objects in water despite short screening lengths, via coupling between osmotic and electrostatic interactions. Widespread work in theory, experiment, and computation has been carried out to gain a fundamental understanding of the rich, yet sometimes counterintuitive, behavior of these polyelectrolyte systems. Examples of polyelectrolyte association in biology include DNA packaging and RNA folding, as well as aggregation and self-organization phenomena in different disease states.

Differential stability of DNA crossovers in solution mediated by divalent cations

The assembly of DNA duplexes into higher-order structures plays a major role in many vital cellular functions such as recombination, chromatin packaging and gene regulation. However, little is currently known about the molecular structure and stability of direct DNA–DNA interactions that are required for such functions. In nature, DNA helices minimize electrostatic repulsion between double helices in several ways. Within crystals, B-DNA forms either right-handed crossovers by groove–backbone interaction or left-handed crossovers by groove–groove juxtaposition. We evaluated the stability of such crossovers at various ionic concentrations using large-scale atomistic molecular dynamics simulations. Our results show that right-handed DNA crossovers are thermodynamically stable in solution in the presence of divalent cations. Attractive forces at short-range stabilize such crossover structures with inter-axial separation of helices less than 20 Å. Right-handed crossovers, however, dissociate swiftly in the presence of monovalent ions only. Surprisingly, left-handed crossovers, assembled by sequence-independent juxtaposition of the helices, appear unstable even at the highest concentration of Mg²⁺studied here. Our study provides new molecular insights into chiral association of DNA duplexes and highlights the unique role divalent cations play in differential stabilization of crossover structures. These results may serve as a rational basis to understand the role DNA crossovers play in biological processes.

Probing the conformational distributions of subpersistence length DNA

We have measured the bending elasticity of short double-stranded DNA (dsDNA) chains through small-angle x-ray scattering from solutions of dsDNA-linked dimers of gold nanoparticles. This method, which does not require exertion of external forces or binding to a substrate, reports on the equilibrium distribution of bending fluctuations, not just an average value (as in ensemble fluorescence resonance energy transfer) or an extreme value (as in cyclization), and in principle provides a more robust data set for assessing the suitability of theoretical models. Our experimental results for dsDNA comprising 42-94 basepairs are consistent with a simple wormlike chain model of dsDNA elasticity, whose behavior we have determined from Monte Carlo simulations that explicitly represent nanoparticles and their alkane tethers. A persistence length of 50 nm (150 basepairs) gave a favorable comparison, consistent with the results of single-molecule force-extension experiments on much longer dsDNA chains, but in contrast to recent suggestions of enhanced flexibility at these length scales.

Both helix topology and counterion distribution contribute to the more effective charge screening in dsRNA compared with dsDNA

The recent discovery of the RNA interference mechanism emphasizes the biological importance of short, isolated, double-stranded (ds) RNA helices and calls for a complete understanding of the biophysical properties of dsRNA. However, most previous studies of the electrostatics of nucleic acid duplexes have focused on DNA. Here, we present a comparative investigation of electrostatic effects in RNA and DNA. Using resonant (anomalous) and non-resonant small-angle X-ray scattering, we characterized the charge screening efficiency and counterion distribution around short (25 bp) dsDNA and RNA molecules of comparable sequence. Consistent with theoretical predictions, we find counterion mediated screening to be more efficient for dsRNA than dsDNA. Furthermore, the topology of the RNA A-form helix alters the spatial distribution of counterions relative to B-form DNA. The experimental results reported here agree well with ion-size-corrected non-linear Poisson-Boltzmann calculations. We propose that differences in electrostatic properties aid in selective recognition of different types of short nucleic acid helices by target binding partners.

A mesoscale model of DNA and its renaturation

A mesoscale model of DNA is presented (3SPN.1), extending the scheme previously developed by our group. Each nucleotide is mapped onto three interaction sites. Solvent is accounted for implicitly through a medium-effective dielectric constant and electrostatic interactions are treated at the level of Debye-Hückel theory. The force field includes a weak, solvent-induced attraction, which helps mediate the renaturation of DNA. Model parameterization is accomplished through replica exchange molecular dynamics simulations of short oligonucleotide sequences over a range of composition and chain length. The model describes the melting temperature of DNA as a function of composition as well as ionic strength, and is consistent with heat capacity profiles from experiments. The dependence of persistence length on ionic strength is also captured by the force field. The proposed model is used to examine the renaturation of DNA. It is found that a typical renaturation event occurs through a nucleation step, whereby an interplay between repulsive electrostatic interactions and colloidal-like attractions allows the system to undergo a series of rearrangements before complete molecular reassociation occurs.

Sequence effects in the melting and renaturation of short DNA oligonucleotides: structure and mechanistic pathways

The renaturation/denaturation of DNA oligonucleotides is characterized in the context of expanded ensemble (EXE) and transition path sampling (TPS) simulations. Free energy profiles have been determined from EXE for DNA sequences of varying composition, chain length, and ionic strength. TPS simulations within a Langevin dynamics formalism have been carried out to obtain further information of the transition state for renaturation. Simulation results reveal that free energy profiles are strikingly similar for the various DNA sequences considered in this work. Taking intact double-stranded DNA to have an extent of reaction ξ = 1.0, the maximum of the free energy profile appears at ξ≈0.15, corresponding to ∼2 base pairs. In terms of chain length, the free energy barrier of longer oligonucleotides (30 versus 15 base pairs) is higher and slightly narrower, due to increased sharpness associated with the transition. Low ionic strength tends to decrease free energy barriers, whereby increasing strand rigidity facilitates reassociation. Two mechanisms for DNA reassociation emerge from our analysis of the transition state ensemble. Repetitive sequences tend to reassociate through a non-specific pathway involving molecular slithering. In contrast, random sequences associate through a more restrictive pathway involving the formation of specific contacts, which then leads to overall molecular zippering. In both random and repetitive sequences, the distribution of contacts suggests that nucleation is favored for sites located within the middle region of the chain. The prevalent extent of reaction for the transition state is ξ≈0.25, and the critical size of the nucleus as obtained from our analysis involves ∼4 base pairs.

Nonmonotonic density dependence of the diffusion of DNA fragments in low-salt suspensions

The high linear charge density of 20-base-pair oligomers of DNA is shown to lead to a striking nonmonotonic dependence of the long-time self-diffusion on the concentration of DNA in low-salt conditions. This generic nonmonotonic behavior results from the strong coupling between the electrostatic and solvent-mediated hydrodynamic interactions, from the renormalization of these electrostatic interactions at large separations, and specifically from the dominance of the far-field hydrodynamic interactions caused by the strong repulsion between the DNA fragments.

Abrupt transition from a free, repulsive to a condensed, attractive DNA phase, induced by multivalent polyamine cations

We have investigated the energetics of DNA condensation by multivalent polyamine cations. Solution small angle x-ray scattering was used to monitor interactions between short 25 base pair dsDNA strands in the free supernatant DNA phase that coexists with the condensed DNA phase. Interestingly, when tetravalent spermine is used, significant inter-DNA repulsion is observed in the free phase, in contrast with the presumed inter-DNA attraction in the coexisting condensed phase. DNA condensation thus appears to be a discrete, first-order-like, transition from a repulsive gaseous to an attractive condensed solid phase, in accord with the reported all-or-none condensation of giant DNA. We further quantify the electrostatic repulsive potentials in the free DNA phase and estimate the number of additional spermine cations that bind to DNA upon condensation.

How DNA coiling enhances target localization by proteins

Many genetic processes depend on proteins interacting with specific sequences on DNA. Despite the large excess of nonspecific DNA in the cell, proteins can locate their targets rapidly. After initial nonspecific binding, they are believed to find the target site by 1D diffusion ("sliding") interspersed by 3D dissociation/reassociation, a process usually referred to as facilitated diffusion. The 3D events combine short intrasegmental "hops" along the DNA contour, intersegmental "jumps" between nearby DNA segments, and longer volume "excursions." The impact of DNA conformation on the search pathway is, however, still unknown. Here, we show direct evidence that DNA coiling influences the specific association rate of EcoRV restriction enzymes. Using optical tweezers together with a fast buffer exchange system, we obtained association times of EcoRV on single DNA molecules as a function of DNA extension, separating intersegmental jumping from other search pathways. Depending on salt concentration, targeting rates almost double when the DNA conformation is changed from fully extended to a coiled configuration. Quantitative analysis by an extended facilitated diffusion model reveals that only a fraction of enzymes are ready to bind to DNA. Generalizing our results to the crowded environment of the cell we predict a major impact of intersegmental jumps on target localization speed on DNA.

Reentrant condensation of proteins in solution induced by multivalent counterions

Negatively charged globular proteins in solution undergo a condensation upon adding trivalent counterions between two critical concentrations C and C, C <C. This reentrant condensation behavior above C is caused by short-ranged electrostatic interactions between multivalent cations and acidic residues, mechanistically different from the case of DNA. Small-angle x-ray scattering indicates a short-ranged attraction between counterion-bound proteins near C and C. Monte Carlo simulations (under these strong electrostatic coupling conditions) support an effective inversion of charge on surface side chains through binding of the multivalent counterions.

RNA folding: conformational statistics, folding kinetics, and ion electrostatics

RNA folding is a remarkably complex problem that involves ion-mediated electrostatic interaction, conformational entropy, base pairing and stacking, and noncanonical interactions. During the past decade, results from a variety of experimental and theoretical studies pointed to (a) the potential ion correlation effect in Mg2+-RNA interactions, (b) the rugged energy landscapes and multistate RNA folding kinetics even for small RNA systems such as hairpins and pseudoknots, (c) the intraloop interactions and sequence-dependent loop free energy, and (d) the strong nonadditivity of chain entropy in RNA pseudoknot and other tertiary folds. Several related issues, which have not been thoroughly resolved, require combined approaches with thermodynamic and kinetic experiments, statistical mechanical modeling, and all-atom computer simulations.