Twist-induced defects of the P-SSP7 genome revealed by modeling the cryo-EM density

Supercoiling and looping promote DNA base accessibility and coordination among distant sites

DNA in cells is supercoiled and constrained into loops and this supercoiling and looping influence every aspect of DNA activity. We show here that negative supercoiling transmits mechanical stress along the DNA backbone to disrupt base pairing at specific distant sites. Cooperativity among distant sites localizes certain sequences to superhelical apices. Base pair disruption allows sharp bending at superhelical apices, which facilitates DNA writhing to relieve torsional strain. The coupling of these processes may help prevent extensive denaturation associated with genomic instability. Our results provide a model for how DNA can form short loops, which are required for many essential processes, and how cells may use DNA loops to position nicks to facilitate repair. Furthermore, our results reveal a complex interplay between site-specific disruptions to base pairing and the 3-D conformation of DNA, which influences how genomes are stored, replicated, transcribed, repaired, and many other aspects of DNA activity.

A Primer on the oxDNA Model of DNA: When to Use it, How to Simulate it and How to Interpret the Results

The oxDNA model of Deoxyribonucleic acid has been applied widely to systems in biology, biophysics and nanotechnology. It is currently available via two independent open source packages. Here we present a set of clearly documented exemplar simulations that simultaneously provide both an introduction to simulating the model, and a review of the model's fundamental properties. We outline how simulation results can be interpreted in terms of-and feed into our understanding of-less detailed models that operate at larger length scales, and provide guidance on whether simulating a system with oxDNA is worthwhile.

Effects of Model Shape, Volume, and Softness of the Capsid for DNA Packaging of phi29

Double-stranded DNA is under extreme confinement when packed in phage phi29 with osmotic pressures approaching 60 atm and densities near liquid crystalline. The shape of the capsid determined from experiment is elongated. We consider the effects of the capsid shape and volume on the DNA distribution. We propose simple models for the capsid of phage phi29 to capture volume, shape, and wall flexibility, leading to an accurate DNA density profile. The effect of the packaging motor twisting the DNA on the resulting density distribution has been explored. We find packing motor induced twisting leads to a greater numbers of defects formed. The emergence of defects such as bubbles or large roll angles along the DNA shows a sequence dependence, and the resulting flexibility leads to an inhomogeneous distribution of defects occurring more often at TpA steps and AT-rich regions. In conjunction with capsid elongation, this has effects on the global DNA packing structures.

Structure and the role of filling rate on model dsDNA packed in a phage capsid

The conformation of DNA inside bacteriophages is of paramount importance for understanding packaging and ejection mechanisms. Models describing the structure of the confined macromolecule have depicted highly ordered conformations, such as spooled or toroidal arrangements that focus on reproducing experimental results obtained by averaging over thousands of configurations. However, it has been seen that more disordered states, including DNA kinking and the presence of domains with different DNA orientation can also accurately reproduce many of the structural experiments. In this work we have compared the results obtained through different simulated filling rates. We find a rate dependence for the resulting constrained states showing different anisotropic configurations. We present a quantitative analysis of the density distribution and the DNA orientation across the capsid showing excellent agreement with structural experiments. Second, we have analyzed the correlations within the capsid, finding evidence of the presence of domains characterized by aligned segments of DNA characterized by the structure factor. Finally, we have measured the number and distribution of DNA defects such as the emergence of bubbles and kinks as function of the filling rate. We find the slower the rate the fewer kink defects that appear and they would be unlikely at experimental filling rates with our model parameters. DNA domains of various orientation get larger with slower rates.

Influence of DNA sequence on the structure of minicircles under torsional stress

The sequence dependence of the conformational distribution of DNA under various levels of torsional stress is an important unsolved problem. Combining theory and coarse-grained simulations shows that the DNA sequence and a structural correlation due to topology constraints of a circle are the main factors that dictate the 3D structure of a 336 bp DNA minicircle under torsional stress. We found that DNA minicircle topoisomers can have multiple bend locations under high torsional stress and that the positions of these sharp bends are determined by the sequence, and by a positive mechanical correlation along the sequence. We showed that simulations and theory are able to provide sequence-specific information about individual DNA minicircles observed by cryo-electron tomography (cryo-ET). We provided a sequence-specific cryo-ET tomogram fitting of DNA minicircles, registering the sequence within the geometric features. Our results indicate that the conformational distribution of minicircles under torsional stress can be designed, which has important implications for using minicircle DNA for gene therapy.

Phage-like packing structures with mean field sequence dependence

Packing of double-stranded DNA in phages must overcome both electrostatic repulsions and the problem of persistence length. We consider coarse-grained models with the ability to kink and with randomly generated disorder. We show that the introduction of kinking into configurations of the DNA polymer packaged within spherical confinement results in significant reductions of the overall energies and pressures. We use a kink model which has the ability to deform every 24 bp, close to the average length predicted from phage sequence. The introduction of such persistence length defects even with highly random packing models increases the local nematic ordering of the packed DNA polymer segments. Such local ordering allowed by kinking not only reduces the total bending energy of confined DNA due to nonlinear elasticity but also reduces the electrostatic component of the energy and pressure. We show that a broad ensemble of polymer configurations is consistent with the structural data. © 2016 Wiley Periodicals, Inc.

Sequence Affects the Cyclization of DNA Minicircles

Understanding how the sequence of a DNA molecule affects its dynamic properties is a central problem affecting biochemistry and biotechnology. The process of cyclizing short DNA, as a critical step in molecular cloning, lacks a comprehensive picture of the kinetic process containing sequence information. We have elucidated this process by using coarse-grained simulations, enhanced sampling methods, and recent theoretical advances. We are able to identify the types and positions of structural defects during the looping process at a base-pair level. Correlations along a DNA molecule dictate critical sequence positions that can affect the looping rate. Structural defects change the bending elasticity of the DNA molecule from a harmonic to subharmonic potential with respect to bending angles. We explore the subelastic chain as a possible model in loop formation kinetics. A sequence-dependent model is developed to qualitatively predict the relative loop formation time as a function of DNA sequence.