A Molecular View of the Dynamics of dsDNA Packing Inside Viral Capsids in the Presence of Ions.

Genome packing in viruses and prokaryotes relies on positively charged ions to reduce electrostatic repulsions, and induce attractions that can facilitate DNA condensation. Here we present molecular dynamics simulations spanning several microseconds of dsDNA packing inside nanometer-sized viral capsids. We use a detailed molecular model of DNA that accounts for molecular structure, basepairing, and explicit counterions. The size and shape of the capsids studied here are based on the 30-nanometer-diameter gene transfer agents of bacterium Rhodobacter capsulatus that transfer random 4.5-kbp (1.5 µm) DNA segments between bacterial cells. Multivalent cations such as spermidine and magnesium induce attraction between packaged DNA sites that can lead to DNA condensation. At high concentrations of spermidine, this condensation significantly increases the shear stresses on the packaged DNA while also reducing the pressure inside the capsid. These effects result in an increase in the packing velocity and the total amount of DNA that can be packaged inside the nanometer-sized capsids. In the simulation results presented here, high concentrations of spermidine3+ did not produce the premature stalling observed in experiments. However, a small increase in the heterogeneity of packing velocities was observed in the systems with magnesium and spermidine ions compared to the system with only salt. The results presented here indicate that the effect of multivalent cations and of spermidine, in particular, on the dynamics of DNA packing, increases with decreasing packing velocities.

Mechanical Response of DNA-Nanoparticle Crystals to Controlled Deformation.

The self-assembly of DNA-conjugated nanoparticles represents a promising avenue toward the design of engineered hierarchical materials. By using DNA to encode nanoscale interactions, macroscale crystals can be formed with mechanical properties that can, at least in principle, be tuned. Here we present in silico evidence that the mechanical response of these assemblies can indeed be controlled, and that subtle modifications of the linking DNA sequences can change the Young's modulus from 97 kPa to 2.1 MPa. We rely on a detailed molecular model to quantify the energetics of DNA-nanoparticle assembly and demonstrate that the mechanical response is governed by entropic, rather than enthalpic, contributions and that the response of the entire network can be estimated from the elastic properties of an individual nanoparticle. The results here provide a first step toward the mechanical characterization of DNA-nanoparticle assemblies, and suggest the possibility of mechanical metamaterials constructed using DNA.

Holliday Junction Thermodynamics and Structure: Coarse-Grained Simulations and Experiments.

Holliday junctions play a central role in genetic recombination, DNA repair and other cellular processes. We combine simulations and experiments to evaluate the ability of the 3SPN.2 model, a coarse-grained representation designed to mimic B-DNA, to predict the properties of DNA Holliday junctions. The model reproduces many experimentally determined aspects of junction structure and stability, including the temperature dependence of melting on salt concentration, the bias between open and stacked conformations, the relative populations of conformers at high salt concentration, and the inter-duplex angle (IDA) between arms. We also obtain a close correspondence between the junction structure evaluated by all-atom and coarse-grained simulations. We predict that, for salt concentrations at physiological and higher levels, the populations of the stacked conformers are independent of salt concentration, and directly observe proposed tetrahedral intermediate sub-states implicated in conformational transitions. Our findings demonstrate that the 3SPN.2 model captures junction properties that are inaccessible to all-atom studies, opening the possibility to simulate complex aspects of junction behavior.

Coarse-Grained Ions for Nucleic Acid Modeling.

We present a general coarse-grained model of sodium, magnesium, spermidine, and chlorine in implicit solvent. The effective potentials between ions are systematically parametrized using a relative entropy coarse-graining approach [Carmichael, S. P. and M. S. Shell, J. Phys. Chem. B, 116, 8383-93 (2012)] that maximizes the information retained in a coarse-grained model. We describe the local distribution of ions in the vicinity of a recently published coarse-grained DNA model and demonstrate a dependence of persistence length on ionic strength that differs from that predicted by Odijk-Skolnick-Fixman theory. Consistent with experimental observations, we show that spermidine induces DNA condensation whereas magnesium and sodium do not. This model can be used alongside any coarse-grained DNA model that has explicit charges and an accurate reproduction of the excluded volume of dsDNA.

A molecular view of DNA-conjugated nanoparticle association energies.

Nanoparticles functionalized with short sequences of DNA represent a promising platform for customizable self assembly. Though much recent research has focused on the phase behavior and assembly of these structures, little has been done to precisely characterize the pairwise interaction between particles. Here we present a detailed calculation of the association between DNA-nanoparticle conjugates using 3SPN.2, a coarse-grained model of DNA that accounts for molecular structure and base-pairing. We compare our results to those obtained experimentally using µm sized particles and analyze the free energy surfaces that characterize interparticle hybridization. Next, we study the importance of three-body effects and their impact on particle association and melting. Lastly, we explore the observation by Park et al. [Nature, 451, 553 (2008)] that DNA-nanoparticle crystallization can be inhibited by the deletion of a single nucleotide. Using our model, we suggest that the role of this nucleotide is to disrupt frustration.

DNA shape dominates sequence affinity in nucleosome formation.

Nucleosomes provide the basic unit of compaction in eukaryotic genomes, and the mechanisms that dictate their position at specific locations along a DNA sequence are of central importance to genetics. In this Letter, we employ molecular models of DNA and proteins to elucidate various aspects of nucleosome positioning. In particular, we show how DNA's histone affinity is encoded in its sequence-dependent shape, including subtle deviations from the ideal straight B-DNA form and local variations of minor groove width. By relying on high-precision simulations of the free energy of nucleosome complexes, we also demonstrate that, depending on DNA's intrinsic curvature, histone binding can be dominated by bending interactions or electrostatic interactions. More generally, the results presented here explain how sequence, manifested as the shape of the DNA molecule, dominates molecular recognition in the problem of nucleosome positioning.

Coarse-grained modeling of DNA curvature.

The interaction of DNA with proteins occurs over a wide range of length scales, and depends critically on its local structure. In particular, recent experimental work suggests that the intrinsic curvature of DNA plays a significant role on its protein-binding properties. In this work, we present a coarse grained model of DNA that is capable of describing base-pairing, hybridization, major and minor groove widths, and local curvature. The model represents an extension of the recently proposed 3SPN.2 description of DNA [D. M. Hinckley, G. S. Freeman, J. K. Whitmer, and J. J. de Pablo, J. Chem. Phys. 139, 144903 (2013)], into which sequence-dependent shape and mechanical properties are incorporated. The proposed model is validated against experimental data including melting temperatures, local flexibilities, dsDNA persistence lengths, and minor groove width profiles.

Coarse-grained modeling of DNA oligomer hybridization: length, sequence, and salt effects.

A recently published coarse-grained DNA model [D. M. Hinckley, G. S. Freeman, J. K. Whitmer, and J. J. de Pablo, J. Chem. Phys. 139, 144903 (2013)] is used to study the hybridization mechanism of DNA oligomers. Forward flux sampling is used to construct ensembles of reactive trajectories from which the effects of sequence, length, and ionic strength are revealed. Heterogeneous sequences are observed to hybridize via the canonical zippering mechanism. In contrast, homogeneous sequences hybridize through a slithering mechanism, while more complex base pair displacement processes are observed for repetitive sequences. In all cases, the formation of non-native base pairs leads to an increase in the observed hybridization rate constants beyond those observed in sequences where only native base pairs are permitted. The scaling of rate constants with length is captured by extending existing hybridization theories to account for the formation of non-native base pairs. Furthermore, that scaling is found to be similar for oligomeric and polymeric systems, suggesting that similar physics is involved.

An experimentally-informed coarse-grained 3-Site-Per-Nucleotide model of DNA: structure, thermodynamics, and dynamics of hybridization.

A new 3-Site-Per-Nucleotide coarse-grained model for DNA is presented. The model includes anisotropic potentials between bases involved in base stacking and base pair interactions that enable the description of relevant structural properties, including the major and minor grooves. In an improvement over available coarse-grained models, the correct persistence length is recovered for both ssDNA and dsDNA, allowing for simulation of non-canonical structures such as hairpins. DNA melting temperatures, measured for duplexes and hairpins by integrating over free energy surfaces generated using metadynamics simulations, are shown to be in quantitative agreement with experiment for a variety of sequences and conditions. Hybridization rate constants, calculated using forward-flux sampling, are also shown to be in good agreement with experiment. The coarse-grained model presented here is suitable for use in biological and engineering applications, including nucleosome positioning and DNA-templated engineering.

A coarse-grain three-site-per-nucleotide model for DNA with explicit ions.

The "three sites per nucleotide" (3SPN) model provides a coarse-grained representation of nucleic acids for simulation of molecular processes. Previously, this model has relied on an implicit representation of the surrounding ionic environment at the level of Debye-Hückel theory. In this work, we eliminate this limitation and present an explicit representation of ions, both monovalent and divalent. The coarse-grain ion-ion and ion-phosphate potential energy functions are inferred from all-atom simulations and parameterized to reproduce key features of the local structure and organization of ions in bulk water and in the presence of DNA. The resulting model, 3SPN.1-I, is capable of reproducing the local structure observed in detailed atomistic simulations, as well as the experimental melting temperature of DNA for a range of DNA oligonucleotide lengths, CG-content, Na(+) concentration, and Mg(2+) concentration.