John S. Schreck

Coarse-graining DNA for simulations of DNA nanotechnology

To simulate long time and length scale processes involving DNA it is necessary to use a coarse-grained description. Here we provide an overview of different approaches to such coarse-graining, focussing on those at the nucleotide level that allow the self-assembly processes associated with DNA nanotechnology to be studied. OxDNA, our recently-developed coarse-grained DNA model, is particularly suited to this task, and has opened up this field to systematic study by simulations. We illustrate some of the range of DNA nanotechnology systems to which the model is being applied, as well as the insights it can provide into fundamental biophysical properties of DNA.

Introducing improved structural properties and salt dependence into a coarse-grained model of DNA

We introduce an extended version of oxDNA, a coarse-grained model of deoxyribonucleic acid (DNA) designed to capture the thermodynamic, structural, and mechanical properties of single- and double-stranded DNA. By including explicit major and minor grooves and by slightly modifying the coaxial stacking and backbone-backbone interactions, we improve the ability of the model to treat large (kilobase-pair) structures, such as DNA origami, which are sensitive to these geometric features. Further, we extend the model, which was previously parameterised to just one salt concentration ([Na(+)] = 0.5M), so that it can be used for a range of salt concentrations including those corresponding to physiological conditions. Finally, we use new experimental data to parameterise the oxDNA potential so that consecutive adenine bases stack with a different strength to consecutive thymine bases, a feature which allows a more accurate treatment of systems where the flexibility of single-stranded regions is important. We illustrate the new possibilities opened up by the updated model, oxDNA2, by presenting results from simulations of the structure of large DNA objects and by using the model to investigate some salt-dependent properties of DNA.

DNA hairpins destabilize duplexes primarily by promoting melting rather than by inhibiting hybridization

The effect of secondary structure on DNA duplex formation is poorly understood. Using oxDNA, a nucleotide level coarse-grained model of DNA, we study how hairpins influence the rate and reaction pathways of DNA hybridzation. We compare to experimental systems studied by Gao et al. (1) and find that 3-base pair hairpins reduce the hybridization rate by a factor of 2, and 4-base pair hairpins by a factor of 10, compared to DNA with limited secondary structure, which is in good agreement with experiments. By contrast, melting rates are accelerated by factors of ∼100 and ∼2000. This surprisingly large speed-up occurs because hairpins form during the melting process, and significantly lower the free energy barrier for dissociation. These results should assist experimentalists in designing sequences to be used in DNA nanotechnology, by putting limits on the suppression of hybridization reaction rates through the use of hairpins and offering the possibility of deliberately increasing dissociation rates by incorporating hairpins into single strands.

Characterizing DNA Star-Tile-Based Nanostructures Using a Coarse-Grained Model.

We use oxDNA, a coarse-grained model of DNA at the nucleotide level, to simulate large nanoprisms that are composed of multi-arm star tiles, in which the size of bulge loops that have been incorporated into the tile design is used to control the flexibility of the tiles. The oxDNA model predicts equilibrium structures for several different nanoprism designs that are in excellent agreement with the experimental structures as measured by cryoTEM. In particular we reproduce the chiral twisting of the top and bottom faces of the nanoprisms, as the bulge sizes in these structures are varied due to the greater flexibility of larger bulges. We are also able to follow how the properties of the star tiles evolve as the prisms are assembled. Individual star tiles are very flexible, but their structures become increasingly well-defined and rigid as they are incorporated into larger assemblies. oxDNA also finds that the experimentally observed prisms are more stable than their inverted counterparts, but interestingly this preference for the arms of the tiles to bend in a given direction only emerges after they are part of larger assemblies. These results show the potential for oxDNA to provide detailed structural insight as well as to predict the properties of DNA nanostructures and hence to aid rational design in DNA nanotechnology.

A statistical mechanical approach to protein aggregation

We develop a theory of aggregation using statistical mechanical methods. An example of a complicated aggregation system with several levels of structures is peptide/protein self-assembly. The problem of protein aggregation is important for the understanding and treatment of neurodegenerative diseases and also for the development of bio-macromolecules as new materials. We write the effective Hamiltonian in terms of interaction energies between protein monomers, protein and solvent, as well as between protein filaments. The grand partition function can be expressed in terms of a Zimm-Bragg-like transfer matrix, which is calculated exactly and all thermodynamic properties can be obtained. We start with two-state and three-state descriptions of protein monomers using Potts models that can be generalized to include q-states, for which the exactly solvable feature of the model remains. We focus on n × N lattice systems, corresponding to the ordered structures observed in some real fibrils. We have obtained results on nucleation processes and phase diagrams, in which a protein property such as the sheet content of aggregates is expressed as a function of the number of proteins on the lattice and inter-protein or interfacial interaction energies. We have applied our methods to Aβ(1-40) and Curli fibrils and obtained results in good agreement with experiments.

Characterizing the Motion of Jointed DNA Nanostructures Using a Coarse-Grained Model

As detailed structural characterizations of large complex DNA nanostructures are hard to obtain experimentally, particularly if they have substantial flexibility, coarse-grained modeling can potentially provide an important complementary role. Such modeling can provide a detailed view of both the average structure and the structural fluctuations, as well as providing insight into how the nanostructures design determines its structural properties. Here, we present a case study of jointed DNA nanostructures using the oxDNA model. In particular, we consider archetypal hinge and sliding joints, as well as more complex structures involving a number of such coupled joints. Our results highlight how the nature of the motion in these structures can sensitively depend on the precise details of the joints. Furthermore, the generally good agreement with experiments illustrates the power of this approach and suggests the use of such modeling to prescreen the properties of putative designs.

DNA bipedal motor walking dynamics: an experimental and theoretical study of the dependency on step size

Abstract We present a detailed coarse-grained computer simulation and single molecule fluorescence study of the walking dynamics and mechanism of a DNA bipedal motor striding on a DNA origami. In particular, we study the dependency of the walking efficiency and stepping kinetics on step size. The simulations accurately capture and explain three different experimental observations. These include a description of the maximum possible step size, a decrease in the walking efficiency over short distances and a dependency of the efficiency on the walking direction with respect to the origami track. The former two observations were not expected and are non-trivial. Based on this study, we suggest three design modifications to improve future DNA walkers. Our study demonstrates the ability of the oxDNA model to resolve the dynamics of complex DNA machines, and its usefulness as an engineering tool for the design of DNA machines that operate in the three spatial dimensions.

Multi-scale coarse-graining for the study of assembly pathways in DNA-brick self-assembly

Inspired by recent successes using single-stranded DNA tiles to produce complex structures, we develop a two-step coarse-graining approach that uses detailed thermodynamic calculations with oxDNA, a nucleotide-based model of DNA, to parametrize a coarser kinetic model that can reach the time and length scales needed to study the assembly mechanisms of these structures. We test the model by performing a detailed study of the assembly pathways for a two-dimensional target structure made up of 334 unique strands each of which are 42 nucleotides long. Without adjustable parameters, the model reproduces a critical temperature for the formation of the assembly that is close to the temperature at which assembly first occurs in experiments. Furthermore, the model allows us to investigate in detail the nucleation barriers and the distribution of critical nucleus shapes for the assembly of a single target structure. The assembly intermediates are compact and highly connected (although not maximally so), and classical nucleation theory provides a good fit to the height and shape of the nucleation barrier at temperatures close to where assembly first occurs.

Research data supporting “Self-assembly of two-dimensional binary quasicrystals: A possible route to a DNA quasicrystal”

Configurations, patch positions and diffraction patterns for the various types of patchy particles used, for the quasicryalline and the approximant phases as well as representative kinetic products, studied in the paper. Initial and final configurations, control files and energy outputs for the oxDNA simulations reported in the paper.