Showing 1–50 of 97 results for author: Doye, J P K

A patchy-particle 3-dimensional octagonal quasicrystal

Authors: Akie Kowaguchi, Savan Mehta, Jonathan P. K. Doye, Eva G. Noya

Abstract: We devise an ideal 3-dimensional octagonal quasicrystal that is based upon the 2-dimensional Ammann-Beenker tiling and that is potentially suitable for realization with patchy particles. Based on an analysis of its local environments we design a binary system of 8- and 5-patch particles that in simulations assembles into a 3-dimensional octagonal quasicrystal. The local structure is subtly differe… ▽ More We devise an ideal 3-dimensional octagonal quasicrystal that is based upon the 2-dimensional Ammann-Beenker tiling and that is potentially suitable for realization with patchy particles. Based on an analysis of its local environments we design a binary system of 8- and 5-patch particles that in simulations assembles into a 3-dimensional octagonal quasicrystal. The local structure is subtly different from the original ideal quasicrystal possessing a narrower coordination-number distribution; in fact, the 8-patch particles are not needed and a one-component system of the 5-patch particles assembles into an essentially identical octagonal quasicrystal. We also consider a one-component system of the 8-patch particles; this assembles into a cluster with a number of crystalline domains, but which, because of the coherent boundaries between the crystallites, has approximate eight-fold order. We envisage that these systems could be realized using DNA origami or protein design. △ Less

Submitted 9 August, 2024; originally announced August 2024.

Comments: main: 6 pages, 4 figures; Supplementary: 12 pages, 6 figures

A one-component patchy-particle icosahedral quasicrystal

Authors: Eva G. Noya, Jonathan P. K. Doye

Abstract: Designing particles that are able to form icosahedral quasicrystals (IQCs) and that are as simple as possible is not only of fundamental interest but is also important to the potential realization of IQCs in materials other than metallic alloys. Here we introduce one-component patchy-particle systems that in simulations are able to form face-centred IQCs that are made up of interconnected icosahed… ▽ More Designing particles that are able to form icosahedral quasicrystals (IQCs) and that are as simple as possible is not only of fundamental interest but is also important to the potential realization of IQCs in materials other than metallic alloys. Here we introduce one-component patchy-particle systems that in simulations are able to form face-centred IQCs that are made up of interconnected icosahedra. The directional bonding of the particles facilitates the formation of a network of bonds with icosahedral orientational order and hence quasiperiodic positional order. The assembled quasicrystals have similar energies to periodic approximants but are entropically stabilized by phason disorder. Their long-range quasiperiodic order is confirmed by a higher-dimensional analysis. Promising routes to realize these IQCs experimentally include via protein design and DNA origami particles. △ Less

Submitted 24 July, 2024; originally announced July 2024.

Comments: main: 8 pages, 3 figures; Supplementary: 24 pages, 15 figures

Coarse-grained modelling of DNA-RNA hybrids

Authors: Eryk J. Ratajczyk, Petr Šulc, Andrew J. Turberfield, Jonathan P. K. Doye, Ard A. Louis

Abstract: We introduce oxNA, a new model for the simulation of DNA-RNA hybrids which is based on two previously developed coarse-grained models$\unicode{x2014}$oxDNA and oxRNA. The model naturally reproduces the physical properties of hybrid duplexes including their structure, persistence length and force-extension characteristics. By parameterising the DNA-RNA hydrogen bonding interaction we fit the model'… ▽ More We introduce oxNA, a new model for the simulation of DNA-RNA hybrids which is based on two previously developed coarse-grained models$\unicode{x2014}$oxDNA and oxRNA. The model naturally reproduces the physical properties of hybrid duplexes including their structure, persistence length and force-extension characteristics. By parameterising 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$\unicode{x2014}$calculating the free energy profiles of hybrid strand displacement reactions, studying the resolution of a short R-loop and simulating RNA-scaffolded wireframe origami. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Comments: 15 pages, 10 figures

Journal ref: J. Chem. Phys. 160, 115101 (2024)

Pipelined information flow in molecular mechanical circuits leads to increased error and irreversibility

Authors: Ian Seet, Thomas E. Ouldridge, Jonathan P. K. Doye

Abstract: Pipelining is a design technique for logical circuits that allows for higher throughput than circuits in which multiple computations are fed through the system one after the other. It allows for much faster computation than architectures in which inputs must pass through every layer of the circuit before the next computation can begin (phased chaining). We explore the hypothesis that these advanta… ▽ More Pipelining is a design technique for logical circuits that allows for higher throughput than circuits in which multiple computations are fed through the system one after the other. It allows for much faster computation than architectures in which inputs must pass through every layer of the circuit before the next computation can begin (phased chaining). We explore the hypothesis that these advantages may be offset by a higher error rate, logical irreversibility and greater thermodynamic costs by simulating pipelined molecular mechanical circuits using an explicit physical model. We observe the emergent logical irreversibility, and see that the simultaneous action of multiple components indeed leads to a higher error rate than in phase-chained circuits. The thermodynamic costs of operating the gates are much larger than in equivalent phase-chained circuits, and these costs do not appear to tend to zero in the limit of slow gate operation. Redesigning the gates to eliminate errors and artificially enforcing logical reversibility reduces the thermodynamic costs and recovers thermodynamically reversible behaviour in the limit of slow gate operation. The breakdown of logical reversibility and accuracy are both associated with a breakdown of the digital behaviour of the device, likely contributing to thermodynamic costs that are large relative to the scale of the information being processed. △ Less

Submitted 20 June, 2023; originally announced June 2023.

Simulation of reversible molecular mechanical logic gates and circuits

Authors: Ian Seet, Thomas E. Ouldridge, Jonathan P. K. Doye

Abstract: Landauer's principle places a fundamental lower limit on the work required to perform a logically irreversible operation. Logically reversible gates provide a way to avoid these work costs, and also simplify the task of making the computation as a whole thermodynamically reversible. The inherent reversibility of mechanical logic gates would make them good candidates for the design of practical log… ▽ More Landauer's principle places a fundamental lower limit on the work required to perform a logically irreversible operation. Logically reversible gates provide a way to avoid these work costs, and also simplify the task of making the computation as a whole thermodynamically reversible. The inherent reversibility of mechanical logic gates would make them good candidates for the design of practical logically reversible computing systems if not for the relatively large size and mass of such systems. In this paper, we outline the design and simulation of reversible molecular mechanical logic gates that come close to the limits of thermodynamic reversibility even under the effects of thermal noise, and outline associated circuit components from which arbitrary combinatorial reversible circuits can be constructed and simulated. We demonstrate that isolated components can be operated in a thermodynamically reversible manner, and explore the complexities of combining components to implement more complex computations. Finally, we demonstrate a method to construct arbitrarily large reversible combinatorial circuits using multiple external controls and signal boosters with a working half-adder circuit. △ Less

Submitted 16 November, 2022; v1 submitted 2 May, 2022; originally announced May 2022.

Journal ref: Phys. Rev. E 107, 024134 (2023)

The free-energy landscape of a mechanically bistable DNA origami

Authors: Chak Kui Wong, Jonathan P. K. Doye

Abstract: Molecular simulations using coarse-grained models allow the structure, dynamics and mechanics of DNA origamis to be comprehensively characterized. Here, we focus on the free-energy landscape of a jointed DNA origami that has been designed to exhibit two mechanically stable states and for which a bistable landscape has been inferred from ensembles of structures visualized by electron microscopy. Su… ▽ More Molecular simulations using coarse-grained models allow the structure, dynamics and mechanics of DNA origamis to be comprehensively characterized. Here, we focus on the free-energy landscape of a jointed DNA origami that has been designed to exhibit two mechanically stable states and for which a bistable landscape has been inferred from ensembles of structures visualized by electron microscopy. Surprisingly, simulations using the oxDNA model predict that the defect-free origami has a single free-energy minimum. The expected second state is not stable because the hinge joints do not simply allow free angular motion but instead lead to increasing free-energetic penalties as the joint angles relevant to the second state are approached. This raises interesting questions about the cause of this difference between simulations and experiment, such as how assembly defects might affect the ensemble of structures observed experimentally. △ Less

Submitted 21 January, 2022; originally announced January 2022.

Comments: main: 7 pages, 6 figures; supplementary material: 6 pages, 8 figures

Journal ref: Appl. Sci. 12, 5875 (2022)

Free-energy landscapes of DNA and its assemblies: Perspectives from coarse-grained modelling

Authors: Jonathan P. K. Doye, Ard A. Louis, John S. Schreck, Flavio Romano, Ryan M. Harrison, Majid Mosayebi, Megan C. Engel, Thomas E. Ouldridge

Abstract: This chapter will provide an overview of how characterizing free-energy landscapes can provide insights into the biophysical properties of DNA, as well as into the behaviour of the DNA assemblies used in the field of DNA nanotechnology. The landscapes for these complex systems are accessible through the use of accurate coarse-grained descriptions of DNA. Particular foci will be the landscapes asso… ▽ More This chapter will provide an overview of how characterizing free-energy landscapes can provide insights into the biophysical properties of DNA, as well as into the behaviour of the DNA assemblies used in the field of DNA nanotechnology. The landscapes for these complex systems are accessible through the use of accurate coarse-grained descriptions of DNA. Particular foci will be the landscapes associated with DNA self-assembly and mechanical deformation, where the latter can arise from either externally imposed forces or internal stresses. △ Less

Submitted 19 November, 2021; originally announced November 2021.

Comments: 20 pages, 5 figures

Journal ref: Frontiers of Nanoscience, Vol. 21, Ch. 9, pp 195-210 (2022)

Characterizing the free-energy landscapes of DNA origamis

Authors: Chak Kui Wong, Chuyan Tang, John S. Schreck, Jonathan P. K. Doye

Abstract: We show how coarse-grained modelling combined with umbrella sampling using distance-based order parameters can be applied to compute the free-energy landscapes associated with mechanical deformations of large DNA nanostructures. We illustrate this approach for the strong bending of DNA nanotubes and the potentially bistable landscape of twisted DNA origami sheets. The homogeneous bending of the DN… ▽ More We show how coarse-grained modelling combined with umbrella sampling using distance-based order parameters can be applied to compute the free-energy landscapes associated with mechanical deformations of large DNA nanostructures. We illustrate this approach for the strong bending of DNA nanotubes and the potentially bistable landscape of twisted DNA origami sheets. The homogeneous bending of the DNA nanotubes is well described by the worm-like chain model; for more extreme bending the nanotubes reversibly buckle with the bending deformations localized at one or two "kinks". For a twisted one-layer DNA origami, the twist is coupled to the bending of the sheet giving rise to a free-energy landscape that has two nearly-degenerate minima that have opposite curvatures. By contrast, for a two-layer origami, the increased stiffness with respect to bending leads to a landscape with a single free-energy minimum that has a saddle-like geometry. The ability to compute such landscapes is likely to be particularly useful for DNA mechanotechnology and for understanding stress accumulation during the self-assembly of origamis into higher-order structures. △ Less

Submitted 14 August, 2021; originally announced August 2021.

Comments: main (17 pages, 7 figures); Supporting Information (14 pages, 11 figures)

Journal ref: Nanoscale 14, 2638-2648 (2022)

Programming patchy particles to form three-dimensional dodecagonal quasicrystals

Authors: Daniel F. Tracey, Eva G. Noya, Jonathan P. K. Doye

Abstract: 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 quasicrystal… ▽ More 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 7-patch particles but with wide patches that allow them to adopt both 7- and 8-coordinated environments. The second is a ternary system that contains a mixture of 7- and 8-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. △ Less

Submitted 9 March, 2021; originally announced March 2021.

Comments: 10 pages, 6 figures

Journal ref: J. Chem. Phys. 154, 194505 (2021)

Measuring internal forces in single-stranded DNA: Application to a DNA force clamp

Authors: Megan C. Engel, Flavio Romano, Ard A. Louis, Jonathan P. K. Doye

Abstract: We present a new method for calculating internal forces in DNA structures using coarse-grained models and demonstrate its utility with the oxDNA model. The instantaneous forces on individual nucleotides are explored and related to model potentials, and using our framework, internal forces are calculated for two simple DNA systems and for a recently-published nanoscopic force clamp. Our results hig… ▽ More We present a new method for calculating internal forces in DNA structures using coarse-grained models and demonstrate its utility with the oxDNA model. The instantaneous forces on individual nucleotides are explored and related to model potentials, and using our framework, internal forces are calculated for two simple DNA systems and for a recently-published nanoscopic force clamp. Our results highlight some pitfalls associated with conventional methods for estimating internal forces, which are based on elastic polymer models, and emphasise the importance of carefully considering secondary structure and ionic conditions when modelling the elastic behaviour of single-stranded DNA. Beyond its relevance to the DNA nanotechnological community, we expect our approach to be broadly applicable to calculations of internal force in a variety of structures -- from DNA to protein -- and across other coarse-grained simulation models. △ Less

Submitted 27 July, 2020; originally announced July 2020.

Comments: 39 pages, 9 figures

Journal ref: J. Chem. Theor. Comput. 16, 7764-7775 (2020)

Computing the elastic mechanical properties of rod-like DNA nanostructures

Authors: Hemani Chhabra, Garima Mishra, Yijing Cao, Domen Prešern, Enrico Skoruppa, Maxime M. C. Tortora, Jonathan P. K. Doye

Abstract: To study the elastic properties of rod-like DNA nanostructures, we perform long simulations of these structure using the oxDNA coarse-grained model. By analysing the fluctuations in these trajectories we obtain estimates of the bend and twist persistence lengths, and the underlying bend and twist elastic moduli and couplings between them. Only on length scales beyond those associated with the spac… ▽ More To study the elastic properties of rod-like DNA nanostructures, we perform long simulations of these structure using the oxDNA coarse-grained model. By analysing the fluctuations in these trajectories we obtain estimates of the bend and twist persistence lengths, and the underlying bend and twist elastic moduli and couplings between them. Only on length scales beyond those associated with the spacings between the interhelix crossovers do the bending fluctuations behave like those of a worm-like chain. The obtained bending persistence lengths are much larger than that for double-stranded DNA and increase non-linearly with the number of helices, whereas the twist moduli increase approximately linearly. To within the numerical error in our data, the twist-bend coupling constants are of order zero. That the bending persistence lengths we obtain are generally somewhat higher than in experiment probably reflects both that the simulated origami have no assembly defects and that the oxDNA extensional modulus for double-stranded DNA is too large. △ Less

Submitted 26 June, 2020; originally announced June 2020.

Journal ref: J. Chem. Theory Comput. 16, 7748-7763 (2020)

Characterising DNA T-motifs by Simulation and Experiment

Authors: Behnam Najafi, Katherine G. Young, Jonathan Bath, Ard A. Louis, Jonathan P. K. Doye, Andrew J. Turberfield

Abstract: The success of DNA nanotechnology has been driven by the discovery of novel structural motifs with a wide range of shapes and uses. We present a comprehensive study of the T-motif, a 3-armed, planar, right-angled junction that has been used in the self-assembly of DNA polyhedra and periodic structures. The motif is formed through the interaction of a bulge loop in one duplex and a sticky end of an… ▽ More The success of DNA nanotechnology has been driven by the discovery of novel structural motifs with a wide range of shapes and uses. We present a comprehensive study of the T-motif, a 3-armed, planar, right-angled junction that has been used in the self-assembly of DNA polyhedra and periodic structures. The motif is formed through the interaction of a bulge loop in one duplex and a sticky end of another. The polarity of the sticky end has significant consequences for the thermodynamic and geometrical properties of the T-motif: different polarities create junctions spanning different grooves of the duplex. We compare experimental binding strengths with predictions of oxDNA, a coarse-grained model of DNA, for various loop sizes. We find that, although both sticky-end polarities can create stable junctions, junctions resulting from 5$'$ sticky ends are stable over a wider range of bulge loop sizes. We highlight the importance of possible coaxial stacking interactions within the motif and investigate how each coaxial stacking interaction stabilises the structure and favours a particular geometry. △ Less

Submitted 23 May, 2020; originally announced May 2020.

Comments: 8 pages; 8 figures

The oxDNA coarse-grained model as a tool to simulate DNA origami

Authors: Jonathan P. K. Doye, Hannah Fowler, Domen Prešern, Joakim Bohlin, Lorenzo Rovigatti, Flavio Romano, Petr Šulc, Chak Kui Wong, Ard A. Louis, John S. Schreck, Megan C. Engel, Michael Matthies, Erik Benson, Erik Poppleton, Benedict E. K. Snodin

Abstract: This chapter introduces how to run molecular dynamics simulations for DNA origami using the oxDNA coarse-grained model. This chapter introduces how to run molecular dynamics simulations for DNA origami using the oxDNA coarse-grained model. △ Less

Submitted 10 April, 2020; originally announced April 2020.

Comments: 17 pages, 5 figures

Programming patchy particles to form complex periodic structures

Authors: Daniel F. Tracey, Eva G. Noya, Jonathan P. K. Doye

Abstract: 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,… ▽ More 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. △ Less

Submitted 23 September, 2019; originally announced September 2019.

Comments: 12 pages, 6 figures

Journal ref: J. Chem. Phys. 151, 224506 (2019)

Chiral shape fluctuations and the origin of chirality in cholesteric phases of DNA origamis

Authors: Maxime M. C. Tortora, Garima Mishra, Domen Presern, Jonathan P. K. Doye

Abstract: Lyotropic cholesteric liquid crystal phases are ubiquitously observed in biological and synthetic polymer solutions, characterized by a complex interplay between thermal fluctuations, entropic and enthalpic forces. The elucidation of the link between microscopic features and macroscopic chiral structure, and of the relative roles of these competing contributions on phase organization, remains a to… ▽ More Lyotropic cholesteric liquid crystal phases are ubiquitously observed in biological and synthetic polymer solutions, characterized by a complex interplay between thermal fluctuations, entropic and enthalpic forces. The elucidation of the link between microscopic features and macroscopic chiral structure, and of the relative roles of these competing contributions on phase organization, remains a topical issue. Here we provide theoretical evidence of a novel mechanism of chirality amplification in lyotropic liquid crystals, whereby phase chirality is governed by fluctuation-stabilized helical deformations in the conformations of their constituent molecules. Our results compare favorably to recent experimental studies of DNA origami assemblies and demonstrate the influence of intra-molecular mechanics on chiral supra-molecular order, with potential implications for a broad class of experimentally-relevant colloidal systems. △ Less

Submitted 1 August, 2020; v1 submitted 29 November, 2018; originally announced November 2018.

Journal ref: Sci. Adv. 6 (2020) eaaw8331

Coarse-grained modelling of the structural properties of DNA origami

Authors: Benedict E. K. Snodin, John S. Schreck, Flavio Romano, Ard A. Louis, Jonathan P. K. Doye

Abstract: We use the oxDNA coarse-grained model to provide a detailed characterization of the fundamental structural properties of DNA origamis, focussing on archetypal 2D and 3D origamis. 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… ▽ More We use the oxDNA coarse-grained model to provide a detailed characterization of the fundamental structural properties of DNA origamis, focussing on archetypal 2D and 3D origamis. 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 Angstrom, 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 origamis, and has the potential to be used to routinely pre-screen putative origami designs. △ Less

Submitted 22 September, 2018; originally announced September 2018.

Comments: 14 pages, 10 figures

Journal ref: Nucleic Acids Res. 47, 1585-1597 (2019)

Incorporating particle flexibility in a density functional description of nematics and cholesterics

Authors: Maxime M. C. Tortora, Jonathan P. K. Doye

Abstract: We describe a general implementation of the Fynewever-Yethiraj density functional theory (DFT) for the investigation of nematic and cholesteric self-assembly in arbitrary solutions of semi-flexible polymers. The basic assumptions of the theory are discussed in the context of other generalised Onsager descriptions for flexible polyatomic systems. The location of the isotropic-to-nematic phase trans… ▽ More We describe a general implementation of the Fynewever-Yethiraj density functional theory (DFT) for the investigation of nematic and cholesteric self-assembly in arbitrary solutions of semi-flexible polymers. The basic assumptions of the theory are discussed in the context of other generalised Onsager descriptions for flexible polyatomic systems. The location of the isotropic-to-nematic phase transition is found to be in good agreement with molecular simulations for elongated chains up to relatively high polymer flexibilities, although the predictions of the theory in the nematic regime lead to gradual underestimations of order parameters with decreasing particle stiffness. This shortcoming is attributed to increasing overestimations of the molecular conformational entropy in higher-density phases, which may not be easily addressed in the formalism of DFT for realistic particle models. Practical consequences of these limitations are illustrated through the application of DFT to systems of near-persistence-length DNA duplexes, whose cholesteric behaviour is found to be strongly contingent on their detailed accessible conformational space in concentrated solutions. △ Less

Submitted 30 January, 2019; v1 submitted 31 January, 2018; originally announced January 2018.

Journal ref: Mol. Phys. 116, 2773-2791 (2018)

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

Authors: Pedro Fonseca, Flavio Romano, John S. Schreck, Thomas E. Ouldridge, Jonathan P. K. Doye, Ard A. Louis

Abstract: 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 perform… ▽ More 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. △ Less

Submitted 28 March, 2018; v1 submitted 6 December, 2017; originally announced December 2017.

Journal ref: J. Chem. Phys 148, 134910 (2018)

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

Authors: Maxime M. C. Tortora, Jonathan P. K. Doye

Abstract: 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 parti… ▽ More 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 behaviour 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^\circ$. 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. △ Less

Submitted 13 December, 2017; v1 submitted 30 August, 2017; originally announced August 2017.

Journal ref: J. Chem. Phys. 147 (2017) 224504

Perturbative density functional methods for cholesteric liquid crystals

Authors: Maxime M. C. Tortora, Jonathan P. K. Doye

Abstract: We introduce a comprehensive numerical framework to generically infer the emergent macroscopic properties of uniaxial nematic and cholesteric phases from that of their microscopic constituent mesogens. This approach, based on the full numerical resolution of the Poniewierski-Stecki equations in the weak chirality limit, may expediently handle a wide range of particle models through the use of Mont… ▽ More We introduce a comprehensive numerical framework to generically infer the emergent macroscopic properties of uniaxial nematic and cholesteric phases from that of their microscopic constituent mesogens. This approach, based on the full numerical resolution of the Poniewierski-Stecki equations in the weak chirality limit, may expediently handle a wide range of particle models through the use of Monte-Carlo sampling for all virial-type integrals. Its predictions in terms of equilibrium cholesteric structures are found to be in excellent agreement with previous full-functional descriptions, thereby demonstrating the quantitative validity of the perturbative treatment of chirality for pitch lengths as short as a few dozen particle diameters. Furthermore, the use of the full angle-dependent virial coefficients in the Onsager-Parsons-Lee formalism increases its numerical efficiency by several orders of magnitude over that of these previous methods. The comparison of our results with numerical simulations however reveals some shortcomings of the Parsons-Lee approximation for systems of strongly non-convex particles, notwithstanding the accurate inclusion of their full effective molecular volume. Further potential limitations of our theory in terms of phase symmetry assumptions are also examined, and prospective directions for future improvements discussed. △ Less

Submitted 20 May, 2018; v1 submitted 21 February, 2017; originally announced February 2017.

Journal ref: J. Chem. Phys. 146 (2017) 184504

DNA-stabilized Ag-Au bimetallic clusters: The effects of alloying and embedding on optical properties

Authors: Dennis Palagin, Jonathan P. K. Doye

Abstract: Global geometry optimization and time-dependent density functional theory calculations have been used to study the structural evolution and optical properties of (AgAu)n (n=2-6) nanoalloys both as individual clusters and as clusters stabilized with the fragments of DNA of different size. We show that alloying can be used to control and tune the level of interaction between the metal atoms of the c… ▽ More Global geometry optimization and time-dependent density functional theory calculations have been used to study the structural evolution and optical properties of (AgAu)n (n=2-6) nanoalloys both as individual clusters and as clusters stabilized with the fragments of DNA of different size. We show that alloying can be used to control and tune the level of interaction between the metal atoms of the cluster and the organic fragments of the DNA ligands. For instance, gold and silver atoms are shown to exhibit synergistic effects in the process of charge transfer from the nucleobase to the cluster, with the silver atoms directly connected to the nitrogen atoms of cytosine increasing their positive partial charge, while their more electronegative neighbouring gold atoms host the excess negative charge. This allows the geometrical structures and optical absorption spectra of small bimetallic clusters to retain many of their main features upon aggregation with relatively large DNA fragments, such as a cytosine-based 9-nucleotide hairpin loop, which suggests a potential synthetic route to such hybrid metal-organic compounds, and opens up the possibility of bringing the unique tunable properties of bimetallic nanoalloys to biological applications. △ Less

Submitted 25 July, 2016; originally announced July 2016.

Journal ref: Phys. Chem. Chem. Phys. 18, 22311-22322 (2016)

Self-assembly of two-dimensional binary quasicrystals: A possible route to a DNA quasicrystal

Authors: Aleks Reinhardt, John S. Schreck, Flavio Romano, Jonathan P. K. Doye

Abstract: We use Monte Carlo simulations and free-energy techniques to show that binary solutions of penta- and hexavalent two-dimensional patchy particles can form thermodynamically stable quasicrystals even at very narrow patch widths, provided their patch interactions are chosen in an appropriate way. Such patchy particles can be thought of as a coarse-grained representation of DNA multi-arm `star' motif… ▽ More We use Monte Carlo simulations and free-energy techniques to show that binary solutions of penta- and hexavalent two-dimensional patchy particles can form thermodynamically stable quasicrystals even at very narrow patch widths, provided their patch interactions are chosen in an appropriate way. Such patchy particles can be thought of as a coarse-grained representation of DNA multi-arm `star' motifs, which can be chosen to bond with one another very specifically by tuning the DNA sequences of the protruding arms. We explore several possible design strategies and conclude that DNA star tiles that are designed to interact with one another in a specific but not overly constrained way could potentially be used to construct soft quasicrystals in experiment. We verify that such star tiles can form stable dodecagonal motifs using oxDNA, a realistic coarse-grained model of DNA. △ Less

Submitted 22 July, 2016; originally announced July 2016.

Journal ref: J. Phys.: Condens. Matter 29, 014006 (2017)

Coarse-grained modelling of strong DNA bending II: Cyclization

Authors: Ryan M. Harrison, Flavio Romano, Thomas E. Ouldridge, Ard A. Louis, Jonathan P. K. Doye

Abstract: DNA cyclization is a powerful technique to gain insight into the nature of DNA bending. The worm-like chain model provides a good description of small to moderate bending fluctuations, but some experiments on strongly-bent shorter molecules suggest enhanced flexibility over and above that expected from the worm-like chain. Here, we use a coarse-grained model of DNA to investigate the thermodynamic… ▽ More DNA cyclization is a powerful technique to gain insight into the nature of DNA bending. The worm-like chain model provides a good description of small to moderate bending fluctuations, but some experiments on strongly-bent shorter molecules suggest enhanced flexibility over and above that expected from the worm-like chain. Here, we use a coarse-grained model of DNA to investigate the thermodynamics of DNA cyclization for molecules with less than 210 base pairs. As the molecules get shorter we find increasing deviations between our computed equilibrium j-factor and the worm-like chain predictions of Shimada and Yamakawa. These deviations are due to sharp kinking, first at nicks, and only subsequently in the body of the duplex. At the shortest lengths, substantial fraying at the ends of duplex domains is the dominant method of relaxation. We also estimate the dynamic j-factor measured in recent FRET experiments. We find that the dynamic j-factor is systematically larger than its equilibrium counterpart, with the deviation larger for shorter molecules, because not all the stress present in the fully cyclized state is present in the transition state. These observations are important for the interpretation of recent experiments, as only kinking within the body of the duplex is genuinely indicative of non-worm-like chain behaviour. △ Less

Submitted 30 June, 2015; originally announced June 2015.

Coarse-grained modelling of strong DNA bending I: Thermodynamics and comparison to an experimental "molecular vice"

Authors: Ryan M. Harrison, Flavio Romano, Thomas E. Ouldridge, Ard A. Louis, Jonathan P. K. Doye

Abstract: DNA bending is biologically important for genome regulation and is relevant to a range of nanotechnological systems. Recent results suggest that sharp bending is much easier than implied by the widely-used worm-like chain model; many of these studies, however, remain controversial. We use a coarse-grained model, previously fitted to DNA's basic thermodynamic and mechanical properties, to explore s… ▽ More DNA bending is biologically important for genome regulation and is relevant to a range of nanotechnological systems. Recent results suggest that sharp bending is much easier than implied by the widely-used worm-like chain model; many of these studies, however, remain controversial. We use a coarse-grained model, previously fitted to DNA's basic thermodynamic and mechanical properties, to explore strongly bent systems. We find that as the end-to-end distance is decreased sufficiently short duplexes undergo a transition to a state in which the bending strain is localized at a flexible kink that involves disruption of base-pairing and stacking. This kinked state, which is not well-described by the worm-like chain model, allows the duplex to more easily be sharply bent. It is not completely flexible, however, due to constraints arising from the connectivity of both DNA backbones. We also perform a detailed comparison to recent experiments on a "molecular vice" that probes highly bent DNA. Close agreement between simulations and experiments strengthens the hypothesis that localised bending via kinking occurs in the molecular vice and causes enhanced flexibility of duplex DNA. Our calculations therefore suggests that the cost of kinking implied by this experiment is consistent with the known thermodynamic and mechanical properties of DNA. △ Less

Submitted 30 June, 2015; originally announced June 2015.

Coarse-grained modelling of supercoiled RNA

Authors: Christian Matek, Petr Šulc, Ferdinando Randisi, Jonathan P. K. Doye, Ard A. Louis

Abstract: We study the behaviour of double-stranded RNA under twist and tension using oxRNA, a recently developed coarse-grained model of RNA. Introducing explicit salt-dependence into the model allows us to directly compare our results to data from recent single-molecule experiments. The model reproduces extension curves as a function of twist and stretching force, including the buckling transition and the… ▽ More We study the behaviour of double-stranded RNA under twist and tension using oxRNA, a recently developed coarse-grained model of RNA. Introducing explicit salt-dependence into the model allows us to directly compare our results to data from recent single-molecule experiments. The model reproduces extension curves as a function of twist and stretching force, including the buckling transition and the behaviour of plectoneme structures. For negative supercoiling, we predict denaturation bubble formation in plectoneme end-loops, suggesting preferential plectoneme localisation in weak base sequences. OxRNA exhibits a positive twist-stretch coupling constant, in agreement with recent experimental observations. △ Less

Submitted 8 June, 2015; originally announced June 2015.

Comments: 8 pages + 5 pages Supplementary Material

Journal ref: J. Chem. Phys. 143, 243122 (2015)

Introducing Improved Structural Properties and Salt Dependence into a Coarse-Grained Model of DNA

Authors: Benedict E. K. Snodin, Ferdinando Randisi, Majid Mosayebi, Petr Sulc, John S. Schreck, Flavio Romano, Thomas E. Ouldridge, Roman Tsukanov, Eyal Nir, Ard A. Louis, Jonathan P. K. Doye

Abstract: We introduce an extended version of oxDNA, a coarse-grained model of 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… ▽ More We introduce an extended version of oxDNA, a coarse-grained model of 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. △ Less

Submitted 19 May, 2015; v1 submitted 3 April, 2015; originally announced April 2015.

Journal ref: J. Chem. Phys. 142, 234901 (2015)

Force-induced rupture of a DNA duplex

Authors: Majid Mosayebi, Ard A. Louis, Jonathan P. K. Doye, Thomas E. Ouldridge

Abstract: The rupture of double-stranded DNA under stress is a key process in biophysics and nanotechnology. In this article we consider the shear-induced rupture of short DNA duplexes, a system that has been given new importance by recently designed force sensors and nanotechnological devices. We argue that rupture must be understood as an activated process, where the duplex state is metastable and the str… ▽ More The rupture of double-stranded DNA under stress is a key process in biophysics and nanotechnology. In this article we consider the shear-induced rupture of short DNA duplexes, a system that has been given new importance by recently designed force sensors and nanotechnological devices. We argue that rupture must be understood as an activated process, where the duplex state is metastable and the strands will separate in a finite time that depends on the duplex length and the force applied. Thus, the critical shearing force required to rupture a duplex within a given experiment depends strongly on the time scale of observation. We use simple models of DNA to demonstrate that this approach naturally captures the experimentally observed dependence of the critical force on duplex length for a given observation time. In particular, the critical force is zero for the shortest duplexes, before rising sharply and then plateauing in the long length limit. The prevailing approach, based on identifying when the presence of each additional base pair within the duplex is thermodynamically unfavorable rather than allowing for metastability, does not predict a time-scale-dependent critical force and does not naturally incorporate a critical force of zero for the shortest duplexes. Additionally, motivated by a recently proposed force sensor, we investigate application of stress to a duplex in a mixed mode that interpolates between shearing and unzipping. As with pure shearing, the critical force depends on the time scale of observation; at a fixed time scale and duplex length, the critical force exhibits a sigmoidal dependence on the fraction of the duplex that is subject to shearing. △ Less

Submitted 12 February, 2015; originally announced February 2015.

Comments: 10 pages, 6 figures

Journal ref: ACS Nano 9, 11993-12003 (2015)

CO Oxidation Catalysed by Pd-based Bimetallic Nanoalloys

Authors: Dennis Palagin, Jonathan P. K. Doye

Abstract: Density functional theory based global geometry optimization has been used to demonstrate the crucial influence of the geometry of the catalytic cluster on the energy barriers for the CO oxidation reaction over Pd-based bimetallic nanoalloys. We show that dramatic geometry change between the reaction intermediates can lead to very high energy barriers and thus be prohibitive for the whole process.… ▽ More Density functional theory based global geometry optimization has been used to demonstrate the crucial influence of the geometry of the catalytic cluster on the energy barriers for the CO oxidation reaction over Pd-based bimetallic nanoalloys. We show that dramatic geometry change between the reaction intermediates can lead to very high energy barriers and thus be prohibitive for the whole process. This introduces challenges for both the design of new catalysts, and theoretical methods employed. On the theory side, a careful choice of geometric configurations of all reaction intermediates is crucial for an adequate description of a possible reaction path. From the point of view of the catalyst design, the cluster geometry can be controlled by adjusting the level of interaction between the cluster and the dopant metal, as well as between the adsorbate molecules and the catalyst cluster by mixing different metals in a single nanoalloy particle. We show that substitution of a Pd atom in the Pd$_{5}$ cluster with a single Ag atom to form Pd$_{4}$Ag$_{1}$ leads to a potential improvement of the catalytic properties of the cluster for the CO oxidation reaction. On the other hand, a single Au atom does not enhance the properties of the catalyst, which is attributed to a weaker hybridization between the cluster's constituent metals and the adsorbate molecules. Such flexibility of properties of bimetallic nanoalloy clusters illustrates the possibility of fine-tuning, which might be used for design of novel efficient catalytic materials. △ Less

Submitted 12 February, 2015; originally announced February 2015.

Comments: 12 pages, 8 figures

Journal ref: Phys. Chem. Chem. Phys. 17, 28010-28021 (2015)

Characterizing the bending and flexibility induced by bulges in DNA duplexes

Authors: John S. Schreck, Thomas E. Ouldridge, Flavio Romano, Ard A. Louis, Jonathan P. K. Doye

Abstract: Advances in DNA nanotechnology have stimulated the search for simple motifs that can be used to control the properties of DNA nanostructures. One such motif, which has been used extensively in structures such as polyhedral cages, two-dimensional arrays, and ribbons, is a bulged duplex, that is two helical segments that connect at a bulge loop. We use a coarse-grained model of DNA to characterize s… ▽ More Advances in DNA nanotechnology have stimulated the search for simple motifs that can be used to control the properties of DNA nanostructures. One such motif, which has been used extensively in structures such as polyhedral cages, two-dimensional arrays, and ribbons, is a bulged duplex, that is two helical segments that connect at a bulge loop. We use a coarse-grained model of DNA to characterize such bulged duplexes. We find that this motif can adopt structures belonging to two main classes: one where the stacking of the helices at the center of the system is preserved, the geometry is roughly straight and the bulge is on one side of the duplex, and the other where the stacking at the center is broken, thus allowing this junction to act as a hinge and increasing flexibility. Small loops favor states where stacking at the center of the duplex is preserved, with loop bases either flipped out or incorporated into the duplex. Duplexes with longer loops show more of a tendency to unstack at the bulge and adopt an open structure. The unstacking probability, however, is highest for loops of intermediate lengths, when the rigidity of single-stranded DNA is significant and the loop resists compression. The properties of this basic structural motif clearly correlate with the structural behavior of certain nano-scale objects, where the enhanced flexibility associated with larger bulges has been used to tune the self-assembly product as well as the detailed geometry of the resulting nanostructures. △ Less

Submitted 19 December, 2014; originally announced December 2014.

Comments: 12 pages + 4 pages of supplemental materials

Journal ref: J. Chem. Phys. 142, 165101 (2015)

Modelling toehold-mediated RNA strand displacement

Authors: Petr Šulc, Thomas E. Ouldridge, Flavio Romano, Jonathan P. K. Doye, Ard A. Louis

Abstract: We study the thermodynamics and kinetics of an RNA toehold-mediated strand displacement reaction with a recently developed coarse-grained model of RNA. Strand displacement, during which a single strand displaces a different strand previously bound to a complementary substrate strand, is an essential mechanism in active nucleic acid nanotechnology and has also been hypothesized to occur in vivo. We… ▽ More We study the thermodynamics and kinetics of an RNA toehold-mediated strand displacement reaction with a recently developed coarse-grained model of RNA. Strand displacement, during which a single strand displaces a different strand previously bound to a complementary substrate strand, is an essential mechanism in active nucleic acid nanotechnology and has also been hypothesized to occur in vivo. We study the rate of displacement reactions as a function of the length of the toehold and temperature and make two experimentally testable predictions: that the displacement is faster if the toehold is placed at the 5' end of the substrate and that the displacement slows down with increasing temperature for longer toeholds. △ Less

Submitted 12 November, 2014; originally announced November 2014.

Journal ref: Biophysical Journal , Volume 108 , Issue 5 , 1238 - 1247, 2015

The role of loop stacking in the dynamics of DNA hairpin formation

Authors: Majid Mosayebi, Flavio Romano, Thomas E. Ouldridge, Ard A. Louis, Jonathan P. K. Doye

Abstract: We study the dynamics of DNA hairpin formation using oxDNA, a nucleotide-level coarse-grained model of DNA. In particular, we explore the effects of the loop stacking interactions and non-native base pairing on the hairpin closing times. We find a non-monotonic variation of the hairpin closing time with temperature, in agreement with the experimental work of Wallace et al. [Proc. Nat. Acad. Sci. U… ▽ More We study the dynamics of DNA hairpin formation using oxDNA, a nucleotide-level coarse-grained model of DNA. In particular, we explore the effects of the loop stacking interactions and non-native base pairing on the hairpin closing times. We find a non-monotonic variation of the hairpin closing time with temperature, in agreement with the experimental work of Wallace et al. [Proc. Nat. Acad. Sci. USA 2001, 98, 5584-5589]. The hairpin closing process involves the formation of an initial nucleus of one or two bonds between the stems followed by a rapid zippering of the stem. At high temperatures, typically above the hairpin melting temperature, an effective negative activation enthalpy is observed because the nucleus has a lower enthalpy than the open state. By contrast, at low temperatures, the activation enthalpy becomes positive mainly due to the increasing energetic cost of bending a loop that becomes increasingly highly stacked as the temperature decreases. We show that stacking must be very strong to induce this experimentally observed behavior, and that the existence of just a few weak stacking points along the loop can substantially suppress it. Non-native base pairs are observed to have only a small effect, slightly accelerating hairpin formation. △ Less

Submitted 5 October, 2014; originally announced October 2014.

Comments: 24 pages, 9 figures

Journal ref: J. Phys. Chem. B, 2014, 118 (49), 14326-14335

Ni-based nanoalloys: Towards thermally stable highly magnetic materials

Authors: Dennis Palagin, Jonathan P. K. Doye

Abstract: Molecular dynamics simulations and density functional theory calculations have been used to demonstrate the possibility of preserving high spin states of the magnetic cores within Ni-based core-shell bimetallic nanoalloys over a wide range of temperatures. We show that, unlike the case of Ni-Al clusters, Ni-Ag clusters preserve high spin states (up to 8 $μ_{\mathrm{B}}$ in case of Ni$_{13}$Ag… ▽ More Molecular dynamics simulations and density functional theory calculations have been used to demonstrate the possibility of preserving high spin states of the magnetic cores within Ni-based core-shell bimetallic nanoalloys over a wide range of temperatures. We show that, unlike the case of Ni-Al clusters, Ni-Ag clusters preserve high spin states (up to 8 $μ_{\mathrm{B}}$ in case of Ni$_{13}$Ag$_{32}$ cluster) due to small hybridization between the electronic levels of two species. Intriguingly, such clusters are also able to maintain geometrical and electronic integrity of their cores at temperatures up to 1000 K (e.g. for Ni$_{7}$Ag$_{27}$ cluster). Furthermore, we also show the possibility of creating ordered arrays of such magnetic clusters on a suitable support by soft-landing pre-formed clusters on the surface, without introducing much disturbance in geometrical and electronic structure of the cluster. We illustrate this approach with the example of Ni$_{13}$Ag$_{38}$ clusters adsorbed on the Si(111)-(7$\times$7) surface, which, having two distinctive halves to the unit cell, acts as a selective template for cluster deposition. △ Less

Submitted 30 September, 2014; originally announced September 2014.

Comments: 9 pages, 10 figures

Journal ref: J. Chem. Phys. 141, 214302 (2014)

DNA hairpins primarily promote duplex melting rather than inhibiting hybridization

Authors: John S. Schreck, Thomas E. Ouldridge, Flavio Romano, Petr Sulc, Liam Shaw, Ard A. Louis, Jonathan P. K. Doye

Abstract: The effect of secondary structure on DNA duplex formation is poorly understood. We use a coarse-grained model of DNA to show that specific 3- and 4-base pair hairpins reduce hybridization rates by factors of 2 and 10 respectively, in good agreement with experiment. By contrast, melting rates are accelerated by factors of ~100 and ~2000. This surprisingly large speed-up occurs because hairpins form… ▽ More The effect of secondary structure on DNA duplex formation is poorly understood. We use a coarse-grained model of DNA to show that specific 3- and 4-base pair hairpins reduce hybridization rates by factors of 2 and 10 respectively, in good agreement with experiment. 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, stabilizing partially melted states, and facilitating dissociation. These results may help guide the design of DNA devices that use hairpins to modulate hybridization and dissociation pathways and rates. △ Less

Submitted 19 August, 2014; originally announced August 2014.

Comments: 5 pages + 14 pages of appendices

Journal ref: Nucl. Acids Res. 43, 6181-6190 (2015)

Heterogeneous ice nucleation on silver-iodide-like surfaces

Authors: Guillaume Fraux, Jonathan P. K. Doye

Abstract: We attempt to simulate the heterogeneous nucleation of ice at model silver-iodide surfaces and find relatively facile ice nucleation and growth at the Ag+ termi nated basal face, but never see nucleation at the I- terminated basal face or the prism and normal faces. Water molecules strongly adsorb onto the Ag+ terminate d face to give a well-ordered hexagonal ice-like bilayer that then acts as a t… ▽ More We attempt to simulate the heterogeneous nucleation of ice at model silver-iodide surfaces and find relatively facile ice nucleation and growth at the Ag+ termi nated basal face, but never see nucleation at the I- terminated basal face or the prism and normal faces. Water molecules strongly adsorb onto the Ag+ terminate d face to give a well-ordered hexagonal ice-like bilayer that then acts as a template for further ice growth. △ Less

Submitted 13 July, 2014; originally announced July 2014.

Comments: 2 pages, 1 figure

Journal ref: J. Chem. Phys. 141, 216101 (2014)

Effects of surface interactions on heterogeneous ice nucleation for a monatomic water model

Authors: Aleks Reinhardt, Jonathan P. K. Doye

Abstract: Despite its importance in atmospheric science, much remains unknown about the microscopic mechanism of heterogeneous ice nucleation. In this work, we perform hybrid Monte Carlo simulations of the heterogeneous nucleation of ice on a range of generic surfaces, both flat and structured, in order to probe the underlying factors affecting the nucleation process. The structured surfaces we study compri… ▽ More Despite its importance in atmospheric science, much remains unknown about the microscopic mechanism of heterogeneous ice nucleation. In this work, we perform hybrid Monte Carlo simulations of the heterogeneous nucleation of ice on a range of generic surfaces, both flat and structured, in order to probe the underlying factors affecting the nucleation process. The structured surfaces we study comprise one basal plane bilayer of ice with varying lattice parameters and interaction strengths. We show that what determines the propensity for nucleation is not just the surface attraction, but also the orientational ordering imposed on liquid water near a surface. In particular, varying the ratio of the surface's attraction and orientational ordering can change the mechanism by which nucleation occurs: ice can nucleate on the structured surface even when the orientational ordering imposed by the surface is weak, as the water molecules that interact strongly with the surface are themselves a good template for further growth. We also show that lattice matching is important for heterogeneous nucleation on the structured surface we study. We rationalise these brute-force simulation results by explicitly calculating the interfacial free energies of ice and liquid water in contact with the nucleating surface and their variation with surface interaction parameters. △ Less

Submitted 22 April, 2014; originally announced April 2014.

Journal ref: J. Chem. Phys. 141, 084501 (2014)

Plectoneme tip bubbles: Coupled denaturation and writhing in supercoiled DNA

Authors: Christian Matek, Thomas E. Ouldridge, Jonathan P. K. Doye, Ard A. Louis

Abstract: Biological information is not only stored in the digital chemical sequence of double helical DNA, but is also encoded in the mechanical properties of the DNA strands, which can influence biochemical processes involving its readout. For example, loop formation in the Lac operon can regulate the expression of key genes, and DNA supercoiling is closely correlated to rhythmic circardian gene expressio… ▽ More Biological information is not only stored in the digital chemical sequence of double helical DNA, but is also encoded in the mechanical properties of the DNA strands, which can influence biochemical processes involving its readout. For example, loop formation in the Lac operon can regulate the expression of key genes, and DNA supercoiling is closely correlated to rhythmic circardian gene expression in cyanobacteria. Supercoiling is also important for large scale organisation of the genome in both eukaryotic and prokaryotic cells. DNA can respond to torsional stress by writhing to form looped structures called plectonemes, thus transferring energy stored as twist into energy stored in bending. Denaturation bubbles can also relax torsional stress, with the enthalpic cost of breaking bonds being compensated by their ability to absorb undertwist. Here we predict a novel regime where bubbles form at the tips of plectonemes, and study its properties using coarse-grained simulations. These tip bubbles can occur for both positive and negative supercoiling and greatly reduce plectoneme diffusion by a pinning mechanism. They can cause plectonemes to preferentially localise to AT rich regions, because bubbles more easily form there. The tip-bubble regime occurs for supercoiling densities and forces that are typically encountered for DNA in vivo, and may be exploited for biological control of genomic processes. △ Less

Submitted 10 April, 2014; originally announced April 2014.

Comments: 5 pages + 26 pages Supplementary Material

Journal ref: Scientific Reports, 5, 7655 (2015)

A nucleotide-level coarse-grained model of RNA

Authors: Petr Šulc, Flavio Romano, Thomas E. Ouldridge, Jonathan P. K. Doye, Ard A. Louis

Abstract: We present a new, nucleotide-level model for RNA, oxRNA, based on the coarse-graining methodology recently developed for the oxDNA model of DNA. The model is designed to reproduce structural, mechanical and thermodynamic properties of RNA, and the coarse-graining level aims to retain the relevant physics for RNA hybridization and the structure of single- and double-stranded RNA. In order to explor… ▽ More We present a new, nucleotide-level model for RNA, oxRNA, based on the coarse-graining methodology recently developed for the oxDNA model of DNA. The model is designed to reproduce structural, mechanical and thermodynamic properties of RNA, and the coarse-graining level aims to retain the relevant physics for RNA hybridization and the structure of single- and double-stranded RNA. In order to explore its strengths and weaknesses, we test the model in a range of nanotechnological and biological settings. Applications explored include the folding thermodynamics of a pseudoknot, the formation of a kissing loop complex, the structure of a hexagonal RNA nanoring, and the unzipping of a hairpin motif. We argue that the model can be used for efficient simulations of the structure of systems with thousands of base pairs, and for the assembly of systems of up to hundreds of base pairs. The source code implementing the model is released for public use. △ Less

Submitted 17 March, 2014; originally announced March 2014.

Comments: 25 pages

Journal ref: J. Chem. Phys. 140, 235102 (2014)

Coarse-graining DNA for simulations of DNA nanotechnology

Authors: Jonathan P. K. Doye, Thomas E. Ouldridge, Ard A. Louis, Flavio Romano, Petr Sulc, Christian Matek, Benedict E. K. Snodin, Lorenzo Rovigatti, John S. Schreck, Ryan M. Harrison, William P. J. Smith

Abstract: 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 sui… ▽ More 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. △ Less

Submitted 18 August, 2013; originally announced August 2013.

Comments: 20 pages, 9 figures

Journal ref: Phys. Chem. Chem. Phys. 15, 20395-20414 (2013)

Homogeneous TIP4P/2005 ice nucleation at low supercooling

Authors: Aleks Reinhardt, Jonathan P. K. Doye

Abstract: We present a partial free energy profile for the homogeneous nucleation of ice using an all-atom model of water at low supercooling, at which ice growth dynamics are reasonably accessible to simulation. We demonstrate that the free energy profile is well described by classical nucleation theory, and that the nucleation barrier is entropic in origin. We also estimate to first order the temperature… ▽ More We present a partial free energy profile for the homogeneous nucleation of ice using an all-atom model of water at low supercooling, at which ice growth dynamics are reasonably accessible to simulation. We demonstrate that the free energy profile is well described by classical nucleation theory, and that the nucleation barrier is entropic in origin. We also estimate to first order the temperature dependence of the interfacial free energy. △ Less

Submitted 21 June, 2013; originally announced June 2013.

Journal ref: J. Chem. Phys. 139, 096102 (2013)

DNA hybridization kinetics: zippering, internal displacement and sequence dependence

Authors: Thomas E. Ouldridge, Petr Šulc, Flavio Romano, Jonathan P. K. Doye, Ard A. Louis

Abstract: While the thermodynamics of DNA hybridization is well understood, much less is known about the kinetics of this classic system. Filling this gap in our understanding has new urgency because DNA nanotechnology often depends critically on binding rates. Here we use a coarse-grained model to explore the hybridization kinetics of DNA oligomers, finding that strand association proceeds through a comple… ▽ More While the thermodynamics of DNA hybridization is well understood, much less is known about the kinetics of this classic system. Filling this gap in our understanding has new urgency because DNA nanotechnology often depends critically on binding rates. Here we use a coarse-grained model to explore the hybridization kinetics of DNA oligomers, finding that strand association proceeds through a complex set of intermediate states. Successful binding events start with the formation of a few metastable base-pairing interactions, followed by zippering of the remaining bonds. However, despite reasonably strong interstrand interactions, initial contacts frequently fail to lead to zippering because the typical configurations in which they form differ from typical states of similar enthalpy in the double-stranded equilibrium ensemble. Therefore, if the association process is analyzed on the base-pair (secondary structure) level, it shows non-Markovian behavior. Initial contacts must be stabilized by two or three base pairs before full zippering is likely, resulting in negative effective activation enthalpies. Non-Arrhenius behavior is observed as the number of base pairs in the effective transition state increases with temperature. In addition, we find that alternative pathways involving misbonds can increase association rates. For repetitive sequences, misaligned duplexes frequently rearrange to form fully paired duplexes by two distinct processes which we label `pseudoknot' and `inchworm' internal displacement. We show how the above processes can explain why experimentally observed association rates of GC-rich oligomers are higher than rates of AT-rich equivalents. More generally, we argue that association rates can be modulated by sequence choice. △ Less

Submitted 14 March, 2013; originally announced March 2013.

Journal ref: Nucleic Acids Research, 41, 8886-8895 (2013)

Computing phase diagrams for a quasicrystal-forming patchy-particle system

Authors: Aleks Reinhardt, Flavio Romano, Jonathan P. K. Doye

Abstract: We introduce an approach to computing the free energy of quasicrystals, which we use to calculate phase diagrams for systems of two-dimensional patchy particles with five regularly arranged patches that have previously been shown to form dodecagonal quasicrystals. We find that the quasicrystal is a thermodynamically stable phase for a wide range of conditions and remains a robust feature of the sy… ▽ More We introduce an approach to computing the free energy of quasicrystals, which we use to calculate phase diagrams for systems of two-dimensional patchy particles with five regularly arranged patches that have previously been shown to form dodecagonal quasicrystals. We find that the quasicrystal is a thermodynamically stable phase for a wide range of conditions and remains a robust feature of the system as the potential's parameters are varied. We also demonstrate that the quasicrystal is entropically stabilised over its crystalline approximants. △ Less

Submitted 8 April, 2013; v1 submitted 11 February, 2013; originally announced February 2013.

Journal ref: Phys. Rev. Lett. 110, 255503 (2013)

Simulating a burnt-bridges DNA motor with a coarse-grained DNA model

Authors: Petr Šulc, Thomas E. Ouldridge, Flavio Romano, Jonathan P. K. Doye, Ard A. Louis

Abstract: We apply a recently-developed coarse-grained model of DNA, designed to capture the basic physics of nanotechnological DNA systems, to the study of a `burnt-bridges' DNA motor consisting of a single-stranded cargo that steps processively along a track of single-stranded stators. We demonstrate that the model is able to simulate such a system, and investigate the sensitivity of the stepping process… ▽ More We apply a recently-developed coarse-grained model of DNA, designed to capture the basic physics of nanotechnological DNA systems, to the study of a `burnt-bridges' DNA motor consisting of a single-stranded cargo that steps processively along a track of single-stranded stators. We demonstrate that the model is able to simulate such a system, and investigate the sensitivity of the stepping process to the spatial separation of stators, finding that an increased distance can suppress successful steps due to the build up of unfavourable tension. The mechanism of suppression suggests that varying the distance between stators could be used as a method for improving signal-to-noise ratios for motors that are required to make a decision at a junction of stators. △ Less

Submitted 18 December, 2012; originally announced December 2012.

Comments: 12 pages

Journal ref: Natural Computing, 13, 535-547 (2014)

Coarse-grained simulations of DNA overstretching

Authors: Flavio Romano, Debayan Chakraborty, Jonathan P. K. Doye, Thomas E. Ouldridge, Ard. A. Louis

Abstract: We use a recently developed coarse-grained model to simulate the overstretching of duplex DNA. Overstretching at 23C occurs at 74 pN in the model, about 6-7 pN higher than the experimental value at equivalent salt conditions. Furthermore, the model reproduces the temperature dependence of the overstretching force well. The mechanism of overstretching is always force-induced melting by unpeeling fr… ▽ More We use a recently developed coarse-grained model to simulate the overstretching of duplex DNA. Overstretching at 23C occurs at 74 pN in the model, about 6-7 pN higher than the experimental value at equivalent salt conditions. Furthermore, the model reproduces the temperature dependence of the overstretching force well. The mechanism of overstretching is always force-induced melting by unpeeling from the free ends. That we never see S-DNA (overstretched duplex DNA), even though there is clear experimental evidence for this mode of overstretching under certain conditions, suggests that S-DNA is not simply an unstacked but hydrogen-bonded duplex, but instead probably has a more exotic structure. △ Less

Submitted 29 January, 2013; v1 submitted 26 September, 2012; originally announced September 2012.

Comments: 11 pages, 11 figures

Journal ref: J. Chem. Phys. 138, 085101 (2013)

Local order parameters for use in driving homogeneous ice nucleation with all-atom models of water

Authors: Aleks Reinhardt, Jonathan P. K. Doye, Eva G. Noya, Carlos Vega

Abstract: We present a local order parameter based on the standard Steinhardt-Ten Wolde approach that is capable both of tracking and of driving homogeneous ice nucleation in simulations of all-atom models of water. We demonstrate that it is capable of forcing the growth of ice nuclei in supercooled liquid water simulated using the TIP4P/2005 model using overbiassed umbrella sampling Monte Carlo simulations… ▽ More We present a local order parameter based on the standard Steinhardt-Ten Wolde approach that is capable both of tracking and of driving homogeneous ice nucleation in simulations of all-atom models of water. We demonstrate that it is capable of forcing the growth of ice nuclei in supercooled liquid water simulated using the TIP4P/2005 model using overbiassed umbrella sampling Monte Carlo simulations. However, even with such an order parameter, the dynamics of ice growth in deeply supercooled liquid water in all-atom models of water are shown to be very slow, and so the computation of free energy landscapes and nucleation rates remains extremely challenging. △ Less

Submitted 8 November, 2012; v1 submitted 29 August, 2012; originally announced August 2012.

Comments: This version incorporates the minor changes made to the paper following peer review

Journal ref: J. Chem. Phys. 137, 194504 (2012)

Sequence-dependent thermodynamics of a coarse-grained DNA model

Authors: Petr Šulc, Flavio Romano, Thomas E. Ouldridge, Lorenzo Rovigatti, Jonathan P. K. Doye, Ard A. Louis

Abstract: We introduce a sequence-dependent parametrization for a coarse-grained DNA model [T. E. Ouldridge, A. A. Louis, and J. P. K. Doye, J. Chem. Phys. 134, 085101 (2011)] originally designed to reproduce the properties of DNA molecules with average sequences. The new parametrization introduces sequence-dependent stacking and base-pairing interaction strengths chosen to reproduce the melting temperature… ▽ More We introduce a sequence-dependent parametrization for a coarse-grained DNA model [T. E. Ouldridge, A. A. Louis, and J. P. K. Doye, J. Chem. Phys. 134, 085101 (2011)] originally designed to reproduce the properties of DNA molecules with average sequences. The new parametrization introduces sequence-dependent stacking and base-pairing interaction strengths chosen to reproduce the melting temperatures of short duplexes. By developing a histogram reweighting technique, we are able to fit our parameters to the melting temperatures of thousands of sequences. To demonstrate the flexibility of the model, we study the effects of sequence on: (a) the heterogeneous stacking transition of single strands, (b) the tendency of a duplex to fray at its melting point, (c) the effects of stacking strength in the loop on the melting temperature of hairpins, (d) the force-extension properties of single strands and (e) the structure of a kissing-loop complex. Where possible we compare our results with experimental data and find a good agreement. A simulation code called oxDNA, implementing our model, is available as free software. △ Less

Submitted 14 July, 2012; originally announced July 2012.

Comments: 15 pages

Journal ref: J. Chem. Phys. 137, 135101 (2012)

DNA cruciform arms nucleate through a correlated but non-synchronous cooperative mechanism

Authors: Christian Matek, Thomas E. Ouldridge, Adam Levy, Jonathan P. K. Doye, Ard A. Louis

Abstract: Inverted repeat (IR) sequences in DNA can form non-canonical cruciform structures to relieve torsional stress. We use Monte Carlo simulations of a recently developed coarse-grained model of DNA to demonstrate that the nucleation of a cruciform can proceed through a cooperative mechanism. Firstly, a twist-induced denaturation bubble must diffuse so that its midpoint is near the centre of symmetry o… ▽ More Inverted repeat (IR) sequences in DNA can form non-canonical cruciform structures to relieve torsional stress. We use Monte Carlo simulations of a recently developed coarse-grained model of DNA to demonstrate that the nucleation of a cruciform can proceed through a cooperative mechanism. Firstly, a twist-induced denaturation bubble must diffuse so that its midpoint is near the centre of symmetry of the IR sequence. Secondly, bubble fluctuations must be large enough to allow one of the arms to form a small number of hairpin bonds. Once the first arm is partially formed, the second arm can rapidly grow to a similar size. Because bubbles can twist back on themselves, they need considerably fewer bases to resolve torsional stress than the final cruciform state does. The initially stabilised cruciform therefore continues to grow, which typically proceeds synchronously, reminiscent of the S-type mechanism of cruciform formation. By using umbrella sampling techniques we calculate, for different temperatures and superhelical densities, the free energy as a function of the number of bonds in each cruciform along the correlated but non-synchronous nucleation pathways we observed in direct simulations. △ Less

Submitted 13 September, 2012; v1 submitted 13 June, 2012; originally announced June 2012.

Comments: 12 pages main paper + 11 pages supplementary data

Journal ref: J. Phys. Chem. B 116, 11616-11625 (2012)

The effect of topology on the structure and free energy landscape of DNA kissing complexes

Authors: Flavio Romano, Alex Hudson, Jonathan P. K. Doye, Thomas E. Ouldridge, Ard A. Louis

Abstract: We use a recently developed coarse-grained model for DNA to study kissing complexes formed by hybridization of complementary hairpin loops. The binding of the loops is topologically constrained because their linking number must remain constant. By studying systems with linking numbers -1, 0 or 1 we show that the average number of interstrand base pairs is larger when the topology is more favourabl… ▽ More We use a recently developed coarse-grained model for DNA to study kissing complexes formed by hybridization of complementary hairpin loops. The binding of the loops is topologically constrained because their linking number must remain constant. By studying systems with linking numbers -1, 0 or 1 we show that the average number of interstrand base pairs is larger when the topology is more favourable for the right-handed wrapping of strands around each other. The thermodynamic stability of the kissing complex also decreases when the linking number changes from -1 to 0 to 1. The structures of the kissing complexes typically involve two intermolecular helices that coaxially stack with the hairpin stems at a parallel four-way junction. △ Less

Submitted 15 March, 2012; originally announced March 2012.

Journal ref: J. Chem. Phys. 136, 215102 (2012)

Formation of dodecagonal quasicrystals in two-dimensional systems of patchy particles

Authors: Marjolein N. van der Linden, Jonathan P. K. Doye, Ard A. Louis

Abstract: The behaviour of two-dimensional patchy particles with 5 and 7 regularly-arranged patches is investigated by computer simulation. For higher pressures and wider patch widths, hexagonal crystals have the lowest enthalpy, whereas at lower pressures and for narrower patches, lower-density crystals with five nearest neighbours and that are based on the (3^2,4,3,4) tiling of squares and triangles becom… ▽ More The behaviour of two-dimensional patchy particles with 5 and 7 regularly-arranged patches is investigated by computer simulation. For higher pressures and wider patch widths, hexagonal crystals have the lowest enthalpy, whereas at lower pressures and for narrower patches, lower-density crystals with five nearest neighbours and that are based on the (3^2,4,3,4) tiling of squares and triangles become lower in enthalpy. Interestingly, in regions of parameter space near to that where the hexagonal crystals become stable, quasicrystalline structures with dodecagonal symmetry form on cooling from high temperature. These quasicrystals can be considered as tilings of squares and triangles, and are probably stabilized by the large configurational entropy associated with all the different possible such tilings. The potential for experimentally realizing such structures using DNA multi-arm motifs are discussed. △ Less

Submitted 24 November, 2011; originally announced November 2011.

Comments: 12 pages, 12 figures

Journal ref: J. Chem. Phys. 136, 054904 (2012)

Free energy landscapes for homogeneous nucleation of ice for a monatomic water model

Authors: Aleks Reinhardt, Jonathan P. K. Doye

Abstract: We simulate the homogeneous nucleation of ice from supercooled liquid water at 220 K in the isobaric-isothermal ensemble using the MW monatomic water potential. Monte Carlo simulations using umbrella sampling are performed in order to determine the nucleation free energy barrier. We find the Gibbs energy profile to be relatively consistent with that predicted by classical nucleation theory; the fr… ▽ More We simulate the homogeneous nucleation of ice from supercooled liquid water at 220 K in the isobaric-isothermal ensemble using the MW monatomic water potential. Monte Carlo simulations using umbrella sampling are performed in order to determine the nucleation free energy barrier. We find the Gibbs energy profile to be relatively consistent with that predicted by classical nucleation theory; the free energy barrier to nucleation was determined to be ~18 kT and the critical nucleus comprised ~85 ice particles. Growth from the supercooled liquid gives clusters that are predominantly cubic, whilst starting with a pre-formed subcritical nucleus of cubic or hexagonal ice results in the growth of predominantly that phase of ice only. △ Less

Submitted 21 December, 2011; v1 submitted 30 September, 2011; originally announced September 2011.

Comments: 11 pages, 6 figures; updated with nucleation rates and additional comparisons with some newly published papers

Journal ref: J. Chem. Phys. 136, 054501 (2012)

Evolutionary Dynamics in a Simple Model of Self-Assembly

Authors: Iain G. Johnston, Sebastian A. Ahnert, Jonathan P. K. Doye, Ard A. Louis

Abstract: We investigate the evolutionary dynamics of an idealised model for the robust self-assembly of two-dimensional structures called polyominoes. The model includes rules that encode interactions between sets of square tiles that drive the self-assembly process. The relationship between the model's rule set and its resulting self-assembled structure can be viewed as a genotype-phenotype map and incorp… ▽ More We investigate the evolutionary dynamics of an idealised model for the robust self-assembly of two-dimensional structures called polyominoes. The model includes rules that encode interactions between sets of square tiles that drive the self-assembly process. The relationship between the model's rule set and its resulting self-assembled structure can be viewed as a genotype-phenotype map and incorporated into a genetic algorithm. The rule sets evolve under selection for specified target structures. The corresponding, complex fitness landscape generates rich evolutionary dynamics as a function of parameters such as the population size, search space size, mutation rate, and method of recombination. Furthermore, these systems are simple enough that in some cases the associated model genome space can be completely characterised, shedding light on how the evolutionary dynamics depends on the detailed structure of the fitness landscape. Finally, we apply the model to study the emergence of the preference for dihedral over cyclic symmetry observed for homomeric protein tetramers. △ Less

Submitted 28 February, 2011; originally announced February 2011.

Journal ref: Phys. Rev. E 83, 066105 (2011)