BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jonathan P. K. Doye

Jonathan P. K. Doye appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.softBiomoleculesBiological Physicsphysics.chem-phcond-mat.stat-mechcond-mat.mtrl-sciPopulations and Evolutioncond-mat.otherMolecular Networks

Catalog footprint

What is connected

44works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

44 published item(s)

preprint2025arXiv

A patchy-particle 3-dimensional octagonal quasicrystal

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The free-energy landscape of a mechanically bistable DNA origami

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Characterizing the free-energy landscapes of DNA origamis

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Characterising DNA T-motifs by Simulation and Experiment

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Computing the elastic mechanical properties of rod-like DNA nanostructures

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

CO Oxidation Catalysed by Pd-based Bimetallic Nanoalloys

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Coarse-grained modelling of supercoiled RNA

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

We introduce an extended version of oxDNA, a coarse-grained model of 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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A nucleotide-level coarse-grained model of RNA

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Characterizing the bending and flexibility induced by bulges in DNA duplexes

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

DNA hairpins primarily promote duplex melting rather than inhibiting hybridization

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Heterogeneous ice nucleation on silver-iodide-like surfaces

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Modelling toehold-mediated RNA strand displacement

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Ni-based nanoalloys: Towards thermally stable highly magnetic materials

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Plectoneme tip bubbles: Coupled denaturation and writhing in supercoiled DNA

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Coarse-grained simulations of DNA overstretching

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

DNA hybridization kinetics: zippering, internal displacement and sequence dependence

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Homogeneous TIP4P/2005 ice nucleation at low supercooling

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Sequence-dependent thermodynamics of a coarse-grained DNA model

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Evolutionary Dynamics in a Simple Model of Self-Assembly

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Reentrant phase behaviour for systems with competition between phase separation and self-assembly

In patchy particle systems where there is competition between the self-assembly of finite clusters and liquid-vapour phase separation, reentrant phase behaviour is observed, with the system passing from a monomeric vapour phase to a region of liquid-vapour phase coexistence and then to a vapour phase of clusters as the temperature is decreased at constant density. Here, we present a classical statistical mechanical approach to the determination of the complete phase diagram of such a system. We model the system as a van der Waals fluid, but one where the monomers can assemble into monodisperse clusters that have no attractive interactions with any of the other species. The resulting phase diagrams show a clear region of reentrance. However, for the most physically reasonable parameter values of the model, this behaviour is restricted to a certain range of density, with phase separation still persisting at high densities.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Structural, mechanical and thermodynamic properties of a coarse-grained DNA model

We explore in detail the structural, mechanical and thermodynamic properties of a coarse-grained model of DNA similar to that introduced in Thomas E. Ouldridge, Ard A. Louis, Jonathan P.K. Doye, Phys. Rev. Lett. 104 178101 (2010). Effective interactions are used to represent chain connectivity, excluded volume, base stacking and hydrogen bonding, naturally reproducing a range of DNA behaviour. We quantify the relation to experiment of the thermodynamics of single-stranded stacking, duplex hybridization and hairpin formation, as well as structural properties such as the persistence length of single strands and duplexes, and the torsional and stretching stiffness of double helices. We also explore the model's representation of more complex motifs involving dangling ends, bulged bases and internal loops, and the effect of stacking and fraying on the thermodynamics of the duplex formation transition.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Templated self-assembly of patchy particles

We explore the use of templated self-assembly to facilitate the formation of complex target structures made from patchy particles. First, we consider the templating of high-symmetry shell structures around a spherical core particle. We find that nucleation around the core particle can inhibit aggregate formation, a process which often hinders self-assembly. In particular, this new assembly pathway allows dodecahedral shells to form readily, whereas these structures never form in the absence of the template. Secondly, we consider the self-assembly of multi-shell structures, where the central icosahedral core is known to form readily on its own, and which could then template the growth of further layers. We are able to find conditions under which two- and three-shell structures successfully assemble, illustrating the power of the templating approach.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The effect of scale-free topology on the robustness and evolvability of genetic regulatory networks

We investigate how scale-free (SF) and Erdos-Renyi (ER) topologies affect the interplay between evolvability and robustness of model gene regulatory networks with Boolean threshold dynamics. In agreement with Oikonomou and Cluzel (2006) we find that networks with SFin topologies, that is SF topology for incoming nodes and ER topology for outgoing nodes, are significantly more evolvable towards specific oscillatory targets than networks with ER topology for both incoming and outgoing nodes. Similar results are found for networks with SFboth and SFout topologies. The functionality of the SFout topology, which most closely resembles the structure of biological gene networks (Babu et al., 2004), is compared to the ER topology in further detail through an extension to multiple target outputs, with either an oscillatory or a non-oscillatory nature. For multiple oscillatory targets of the same length, the differences between SFout and ER networks are enhanced, but for non-oscillatory targets both types of networks show fairly similar evolvability. We find that SF networks generate oscillations much more easily than ER networks do, and this may explain why SF networks are more evolvable than ER networks are for oscillatory phenotypes. In spite of their greater evolvability, we find that networks with SFout topologies are also more robust to mutations than ER networks. Furthermore, the SFout topologies are more robust to changes in initial conditions (environmental robustness). For both topologies, we find that once a population of networks has reached the target state, further neutral evolution can lead to an increase in both the mutational robustness and the environmental robustness to changes in initial conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The stability of a crystal with diamond structure for patchy particles with tetrahedral symmetry

The phase diagram of model anisotropic particles with four attractive patches in a tetrahedral arrangement has been computed at two different values for the range of the potential, with the aim of investigating the conditions under which a diamond crystal can be formed. We find that the diamond phase is never stable for our longer-ranged potential. At low temperatures and pressures, the fluid freezes into a body-centred-cubic solid that can be viewed as two interpenetrating diamond lattices with a weak interaction between the two sublattices. Upon compression, an orientationally ordered face-centred-cubic crystal becomes more stable than the body-centred-cubic crystal, and at higher temperatures a plastic face-centered-cubic phase is stabilized by the increased entropy due to orientational disorder. A similar phase diagram is found for the shorter-ranged potential, but at low temperatures and pressures, we also find a region over which the diamond phase is thermodynamically favored over the body-centred-cubic phase. The higher vibrational entropy of the diamond structure with respect to the body-centred-cubic solid explains why it is stable even though the enthalpy of the latter phase is lower. Some preliminary studies on the growth of the diamond structure starting from a crystal seed were performed. Even though the diamond phase is never thermodynamically stable for the longer-ranged model, direct coexistence simulations of the interface between the fluid and the body-centred-cubic crystal and between the fluid and the diamond crystal show that, at sufficiently low pressures, it is quite probable that in both cases the solid grows into a diamond crystal, albeit involving some defects. These results highlight the importance of kinetic effects in the formation of diamond crystals in systems of patchy particles.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

DNA nanotweezers studied with a coarse-grained model of DNA

We introduce a coarse-grained rigid nucleotide model of DNA that reproduces the basic thermodynamics of short strands: duplex hybridization, single-stranded stacking and hairpin formation, and also captures the essential structural properties of DNA: the helical pitch, persistence length and torsional stiffness of double-stranded molecules, as well as the comparative flexibility of unstacked single strands. We apply the model to calculate the detailed free-energy landscape of one full cycle of DNA 'tweezers', a simple machine driven by hybridization and strand displacement.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Extracting bulk properties of self-assembling systems from small simulations

For systems that self assemble into finite-sized objects, it is sometimes convenient to compute the thermodynamics for a small system where a single assembly can form. However, we show that in the canonical ensemble the use of small systems can lead to significant finite-size effects due to the suppression of concentration fluctuations. We introduce methods to estimate the bulk-yields from simulations of small systems and to follow the convergence of yields with system size, under the assumptions that the various species behave ideally. We also propose an extension to the umbrella sampling technique that allows the formation of multiple finite-sized objects.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Modelling the Self-Assembly of Virus Capsids

We use computer simulations to study a model, first proposed by Wales [1], for the reversible and monodisperse self-assembly of simple icosahedral virus capsid structures. The success and efficiency of assembly as a function of thermodynamic and geometric factors can be qualitatively related to the potential energy landscape structure of the assembling system. Even though the model is strongly coarse-grained, it exhibits a number of features also observed in experiments, such as sigmoidal assembly dynamics, hysteresis in capsid formation and numerous kinetic traps. We also investigate the effect of macromolecular crowding on the assembly dynamics. Crowding agents generally reduce capsid yields at optimal conditions for non-crowded assembly, but may increase yields for parameter regimes away from the optimum. Finally, we generalize the model to a larger triangulation number T = 3, and observe more complex assembly dynamics than that seen for the original T = 1 model.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Monodisperse self-assembly in a model with protein-like interactions

We study the self-assembly behaviour of patchy particles with `protein-like' interactions that can be considered as a minimal model for the assembly of viral capsids and other shell-like protein complexes. We thoroughly explore the thermodynamics and dynamics of self assembly as a function of the parameters of the model and find robust assembly of all target structures considered. Optimal assembly occurs in the region of parameter space where a free energy barrier regulates the rate of nucleation, thus preventing the premature exhaustion of the supply of monomers that can lead to the formation of incomplete shells. The interactions also need to be specific enough to prevent the assembly of malformed shells, but whilst maintaining kinetic accessibility. Free-energy landscapes computed for our model have a funnel-like topography guiding the system to form the target structure, and show that the torsional component of the interparticle interactions prevents the formation of disordered aggregates that would otherwise act as kinetic traps.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Self-assembly of monodisperse clusters: Dependence on target geometry

We apply a simple model system of patchy particles to study monodisperse self-assembly, using the Platonic solids as target structures. We find marked differences between the assembly behaviours of the different systems. Tetrahedra, octahedra and icosahedra assemble easily, while cubes are more challenging and dodecahedra do not assemble. We relate these differences to the kinetics and thermodynamics of assembly, with the formation of large disordered aggregates a particular important competitor to correct assembly. In particular, the free energy landscapes of those targets that are easy to assemble are funnel-like, whereas for the dodecahedral system the landscape is relatively flat with little driving force to facilitate escape from disordered aggregates.

0 citations0 reviewsNo code linkedOpen access