How to design an icosahedral quasicrystal through directional bonding.

Icosahedral quasicrystals (IQCs) are materials that exhibit long-range order but lack periodicity in any direction. Although IQCs were the first reported quasicrystals1, they have been experimentally observed only in metallic alloys2, not in other materials. By contrast, quasicrystals with other symmetries (particularly dodecagonal) have now been found in several soft-matter systems3-5. Here we introduce a class of IQCs built from model patchy colloids that could be realized experimentally using DNA origami particles. Our rational design strategy leads to systems that robustly assemble in simulations into a target IQC through directional bonding. This is illustrated for both body-centred and primitive IQCs, with the simplest systems involving just two particle types. The key design feature is the geometry of the interparticle interactions favouring the propagation of an icosahedral network of bonds, despite this leading to many particles not being fully bonded. As well as furnishing model systems in which to explore the fundamental physics of IQCs, our approach provides a potential route towards functional quasicrystalline materials.

Symmetry and simplicity spontaneously emerge from the algorithmic nature of evolution.

SignificanceWhy does evolution favor symmetric structures when they only represent a minute subset of all possible forms? Just as monkeys randomly typing into a computer language will preferentially produce outputs that can be generated by shorter algorithms, so the coding theorem from algorithmic information theory predicts that random mutations, when decoded by the process of development, preferentially produce phenotypes with shorter algorithmic descriptions. Since symmetric structures need less information to encode, they are much more likely to appear as potential variation. Combined with an arrival-of-the-frequent mechanism, this algorithmic bias predicts a much higher prevalence of low-complexity (high-symmetry) phenotypes than follows from natural selection alone and also explains patterns observed in protein complexes, RNA secondary structures, and a gene regulatory network.

Coarse-grained modeling of DNA-RNA hybrids.

We introduce oxNA, a new model for the simulation of DNA-RNA hybrids that is based on two previously developed coarse-grained models-oxDNA and oxRNA. The model naturally reproduces the physical properties of hybrid duplexes, including their structure, persistence length, and force-extension characteristics. By parameterizing the DNA-RNA hydrogen bonding interaction, we fit the model's thermodynamic properties to experimental data using both average-sequence and sequence-dependent parameters. To demonstrate the model's applicability, we provide three examples of its use-calculating the free energy profiles of hybrid strand displacement reactions, studying the resolution of a short R-loop, and simulating RNA-scaffolded wireframe origami.

The interplay of supercoiling and thymine dimers in DNA.

Thymine dimers are a major mutagenic photoproduct induced by UV radiation. While they have been the subject of extensive theoretical and experimental investigations, questions of how DNA supercoiling affects local defect properties, or, conversely, how the presence of such defects changes global supercoiled structure, are largely unexplored. Here, we introduce a model of thymine dimers in the oxDNA forcefield, parametrized by comparison to melting experiments and structural measurements of the thymine dimer induced bend angle. We performed extensive molecular dynamics simulations of double-stranded DNA as a function of external twist and force. Compared to undamaged DNA, the presence of a thymine dimer lowers the supercoiling densities at which plectonemes and bubbles occur. For biologically relevant supercoiling densities and forces, thymine dimers can preferentially segregate to the tips of the plectonemes, where they enhance the probability of a localized tip-bubble. This mechanism increases the probability of highly bent and denatured states at the thymine dimer site, which may facilitate repair enzyme binding. Thymine dimer-induced tip-bubbles also pin plectonemes, which may help repair enzymes to locate damage. We hypothesize that the interplay of supercoiling and local defects plays an important role for a wider set of DNA damage repair systems.

Designer protein assemblies with tunable phase diagrams in living cells.

Protein self-organization is a hallmark of biological systems. Although the physicochemical principles governing protein-protein interactions have long been known, the principles by which such nanoscale interactions generate diverse phenotypes of mesoscale assemblies, including phase-separated compartments, remain challenging to characterize. To illuminate such principles, we create a system of two proteins designed to interact and form mesh-like assemblies. We devise a new strategy to map high-resolution phase diagrams in living cells, which provide self-assembly signatures of this system. The structural modularity of the two protein components allows straightforward modification of their molecular properties, enabling us to characterize how interaction affinity impacts the phase diagram and material state of the assemblies in vivo. The phase diagrams and their dependence on interaction affinity were captured by theory and simulations, including out-of-equilibrium effects seen in growing cells. Finally, we find that cotranslational protein binding suffices to recruit a messenger RNA to the designed micron-scale structures.

Massively Parallelized Molecular Force Manipulation with On-Demand Thermal and Optical Control.

In single-molecule force spectroscopy (SMFS), a tethered molecule is stretched using a specialized instrument to study how macromolecules extend under force. One problem in SMFS is the serial and slow nature of the measurements, performed one molecule at a time. To address this long-standing challenge, we report on the origami polymer force clamp (OPFC) which enables parallelized manipulation of the mechanical forces experienced by molecules without the need for dedicated SMFS instruments or surface tethering. The OPFC positions target molecules between a rigid nanoscale DNA origami beam and a responsive polymer particle that shrinks on demand. As a proof-of-concept, we record the steady state and time-resolved mechanical unfolding dynamics of DNA hairpins using the fluorescence signal from ensembles of molecules and confirm our conclusion using modeling.

Programming patchy particles to form three-dimensional dodecagonal quasicrystals.

Model patchy particles have been shown to be able to form a wide variety of structures, including symmetric clusters, complex crystals, and even two-dimensional quasicrystals. Here, we investigate whether we can design patchy particles that form three-dimensional quasicrystals, in particular targeting a quasicrystal with dodecagonal symmetry that is made up of stacks of two-dimensional quasicrystalline layers. We obtain two designs that are able to form such a dodecagonal quasicrystal in annealing simulations. The first is a one-component system of seven-patch particles but with wide patches that allow them to adopt both seven- and eight-coordinated environments. The second is a ternary system that contains a mixture of seven- and eight-patch particles and is likely to be more realizable in experiments, for example, using DNA origami. One interesting feature of the first system is that the resulting quasicrystals very often contain a screw dislocation.

Coarse-grained modelling of the structural properties of DNA origami.

We use the oxDNA coarse-grained model to provide a detailed characterization of the fundamental structural properties of DNA origami, focussing on archetypal 2D and 3D origami. The model reproduces well the characteristic pattern of helix bending in a 2D origami, showing that it stems from the intrinsic tendency of anti-parallel four-way junctions to splay apart, a tendency that is enhanced both by less screened electrostatic interactions and by increased thermal motion. We also compare to the structure of a 3D origami whose structure has been determined by cryo-electron microscopy. The oxDNA average structure has a root-mean-square deviation from the experimental structure of 8.4 Å, which is of the order of the experimental resolution. These results illustrate that the oxDNA model is capable of providing detailed and accurate insights into the structure of DNA origami, and has the potential to be used to routinely pre-screen putative origami designs and to investigate the molecular mechanisms that regulate the properties of DNA origami.

Design and synthesis of pleated DNA origami nanotubes with adjustable diameters.

DNA origami allows for the synthesis of nanoscale structures and machines with nanometre precision and high yields. Tubular DNA origami nanostructures are particularly useful because their geometry facilitates a variety of applications including nanoparticle encapsulation, the construction of artificial membrane pores and as structural scaffolds that can uniquely spatially arrange nanoparticles in circular, linear and helical arrays. Here we report a system of parametrization for the design of radially symmetric DNA origami nanotubes with adjustable diameter, length, crossover density, pleat angle and chirality. The system is implemented into a computational algorithm that provides a practical means to navigate the complex geometry of DNA origami nanotube design. We apply this in the design, synthesis and characterization of novel DNA origami nanotubes. These include structures with pleated walls where the same number of duplexes can form nanotubes with different diameters, and to vary the diameter within the same structure. We also construct nanotubes that can be reconfigured into different chiral shapes. Finally, we explore the effect of strain on the local and global geometry of DNA origami nanotubes and demonstrate how pleated walls can provide a strategy to rigidify nanotubes and to construct closely packed parallel duplexes.

Substrate conformational dynamics facilitate structure-specific recognition of gapped DNA by DNA polymerase.

DNA-binding proteins utilise different recognition mechanisms to locate their DNA targets; some proteins recognise specific DNA sequences, while others interact with specific DNA structures. While sequence-specific DNA binding has been studied extensively, structure-specific recognition mechanisms remain unclear. Here, we study structure-specific DNA recognition by examining the structure and dynamics of DNA polymerase I Klenow Fragment (Pol) substrates both alone and in DNA-Pol complexes. Using a docking approach based on a network of 73 distances collected using single-molecule FRET, we determined a novel solution structure of the single-nucleotide-gapped DNA-Pol binary complex. The structure resembled existing crystal structures with regards to the downstream primer-template DNA substrate, and revealed a previously unobserved sharp bend (∼120°) in the DNA substrate; this pronounced bend was present in living cells. MD simulations and single-molecule assays also revealed that 4-5 nt of downstream gap-proximal DNA are unwound in the binary complex. Further, experiments and coarse-grained modelling showed the substrate alone frequently adopts bent conformations with 1-2 nt fraying around the gap, suggesting a mechanism wherein Pol recognises a pre-bent, partially-melted conformation of gapped DNA. We propose a general mechanism for substrate recognition by structure-specific enzymes driven by protein sensing of the conformational dynamics of their DNA substrates.

Rolling circle amplification shows a sinusoidal template length-dependent amplification bias.

Biophysical properties of DNA such as its longitudinal and torsional persistence length govern many processes and phenomena in biology, DNA nanotechnology and biotechnology. It has, for example, long been known that the circularization efficiency of short DNA fragments shows a periodic pattern where fragments with integer helical turns circularize much more efficiently than those with odd helical half turns due to stronger stacking of duplex ends. Small DNA circles can serve as templates for rolling circle amplification (RCA), which is a common and extremely robust amplification mechanism for nucleic acids. We discovered a strong template length-dependent amplification efficiency bias of RCA with the same periodicity as B-DNA. However, stacking cannot explain the mechanism behind this bias as the presence of the polymerase in the bifurcation fork inhibits base stacking of ends. Instead, coarse-grained molecular dynamics simulations imply that different amplification efficiencies come from a varying fraying probability of the last two downstream base pairs. We conclude that an increased strain-promoted fraying probability can increase the polymerization rate compared to a relaxed template.

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

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.

Reconfigurable T-junction DNA Origami.

DNA self-assembly allows the construction of nanometre-scale structures and devices. Structures with thousands of unique components are routinely assembled in good yield. Experimental progress has been rapid, based largely on empirical design rules. Herein, we demonstrate a DNA origami technique designed as a model system with which to explore the mechanism of assembly. The origami fold is controlled through single-stranded loops embedded in a double-stranded DNA template and is programmed by a set of double-stranded linkers that specify pairwise interactions between loop sequences. Assembly is via T-junctions formed by hybridization of single-stranded overhangs on the linkers with the loops. The sequence of loops on the template and the set of interaction rules embodied in the linkers can be reconfigured with ease. We show that a set of just two interaction rules can be used to assemble simple T-junction origami motifs and that assembly can be performed at room temperature.

TacoxDNA: A user-friendly web server for simulations of complex DNA structures, from single strands to origami.

Simulations of nucleic acids at different levels of structural details are increasingly used to complement and interpret experiments in different fields, from biophysics to medicine and materials science. However, the various structural models currently available for DNA and RNA and their accompanying suites of computational tools can be very rarely used in a synergistic fashion. The tacoxDNA webserver and standalone software package presented here are a step toward a long-sought interoperability of nucleic acids models. The webserver offers a simple interface for converting various common input formats of DNA structures and setting up molecular dynamics (MD) simulations. Users can, for instance, design DNA rings with different topologies, such as knots, with and without supercoiling, by simply providing an XYZ coordinate file of the DNA centre-line. More complex DNA geometries, as designable in the cadnano, CanDo and Tiamat tools, can also be converted to all-atom or oxDNA representations, which can then be used to run MD simulations. Though the latter are currently geared toward the native and LAMMPS oxDNA representations, the open-source package is designed to be further expandable. TacoxDNA is available at http://tacoxdna.sissa.it. © 2019 Wiley Periodicals, Inc.

Programming patchy particles to form complex periodic structures.

We introduce a scheme to design patchy particles so that a given target crystal is the global free-energy minimum at sufficiently low temperature. A key feature is a torsional component to the potential that only allows binding when particles have the correct relative orientations. In all examples studied, the target crystal structures readily assembled on annealing from a low-density fluid phase, albeit with the simpler target structures assembling more rapidly. The most complex example was a clathrate with 46 particles in its primitive unit cell. We also explored whether the structural information encoded in the particle interactions could be further reduced. For example, removing the torsional restrictions led to the assembly of an alternative crystal structure for the BC8-forming design, but the more complex clathrate design was still able to assemble because of the greater remaining specificity.

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.

Hierarchical bounding structures for efficient virial computations: Towards a realistic molecular description of cholesterics.

We detail the application of bounding volume hierarchies to accelerate second-virial evaluations for arbitrary complex particles interacting through hard and soft finite-range potentials. This procedure, based on the construction of neighbour lists through the combined use of recursive atom-decomposition techniques and binary overlap search schemes, is shown to scale sub-logarithmically with particle resolution in the case of molecular systems with high aspect ratios. Its implementation within an efficient numerical and theoretical framework based on classical density functional theory enables us to investigate the cholesteric self-assembly of a wide range of experimentally relevant particle models. We illustrate the method through the determination of the cholesteric behavior of hard, structurally resolved twisted cuboids, and report quantitative evidence of the long-predicted phase handedness inversion with increasing particle thread angles near the phenomenological threshold value of 45°. Our results further highlight the complex relationship between microscopic structure and helical twisting power in such model systems, which may be attributed to subtle geometric variations of their chiral excluded-volume manifold.

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.