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

Ard A. Louis

Ard A. Louis 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-phPopulations and Evolutioncond-mat.stat-mechMolecular NetworksComputational ComplexityComputational Geometrycond-mat.mtrl-sciInformation Theorymath.IT

Catalog footprint

What is connected

43works
12topics
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

43 published item(s)

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
preprint2021arXiv

From genotypes to organisms: State-of-the-art and perspectives of a cornerstone in evolutionary dynamics

Understanding how genotypes map onto phenotypes, fitness, and eventually organisms is arguably the next major missing piece in a fully predictive theory of evolution. We refer to this generally as the problem of the genotype-phenotype map. Though we are still far from achieving a complete picture of these relationships, our current understanding of simpler questions, such as the structure induced in the space of genotypes by sequences mapped to molecular structures, has revealed important facts that deeply affect the dynamical description of evolutionary processes. Empirical evidence supporting the fundamental relevance of features such as phenotypic bias is mounting as well, while the synthesis of conceptual and experimental progress leads to questioning current assumptions on the nature of evolutionary dynamics-cancer progression models or synthetic biology approaches being notable examples. This work delves into a critical and constructive attitude in our current knowledge of how genotypes map onto molecular phenotypes and organismal functions, and discusses theoretical and empirical avenues to broaden and improve this comprehension. As a final goal, this community should aim at deriving an updated picture of evolutionary processes soundly relying on the structural properties of genotype spaces, as revealed by modern techniques of molecular and functional analysis.

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

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

Genetic correlations greatly increase mutational robustness and can both reduce and enhance evolvability

Mutational neighbourhoods in genotype-phenotype (GP) maps are widely believed to be more likely to share characteristics than expected from random chance. Such genetic correlations should, as John Maynard Smith famously pointed out, strongly influence evolutionary dynamics. We explore and quantify these intuitions by comparing three GP maps - RNA SS, HP for tertiary, Polyominoes for protein quaternary structure - to a simple random null model that maintains the number of genotypes mapping to each phenotype, but assigns genotypes randomly. The mutational neighbourhood of a genotype in these GP maps is much more likely to contain (mutationally neutral) genotypes mapping to the same phenotype than in the random null model. These neutral correlations can increase the robustness to mutations by orders of magnitude over that of the null model, raising robustness above the critical threshold for the formation of large neutral networks that enhance the capacity for neutral exploration. We also study {\em non-neutral correlations}: Compared to the null model, i) If a particular (non-neutral) phenotype is found once in the 1-mutation neighbourhood of a genotype, then the chance of finding that phenotype multiple times in this neighbourhood is larger than expected; ii) If two genotypes are connected by a single neutral mutation, then their respective non-neutral 1-mutation neighbourhoods are more likely to be similar; iii) If a genotype maps to a folding or self-assembling phenotype, then its non-neutral neighbours are less likely to be a potentially deleterious non-folding or non-assembling phenotype. Non-neutral correlations of type i) and ii) reduce the rate at which new phenotypes can be found by neutral exploration, and so may diminish evolvability, while non-neutral correlations of type iii) may instead facilitate evolutionary exploration and so increase evolvability.

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
preprint2015arXiv

The structure of the genotype-phenotype map strongly constrains the evolution of non-coding RNA

The prevalence of neutral mutations implies that biological systems typically have many more genotypes than phenotypes. But can the way that genotypes are distributed over phenotypes determine evolutionary outcomes? Answering such questions is difficult because the number of genotypes can be hyper-astronomically large. By solving the genotype-phenotype (GP) map for RNA secondary structure for systems up to length $L=126$ nucleotides (where the set of all possible RNA strands would weigh more than the mass of the visible universe) we show that the GP map strongly constrains the evolution of non-coding RNA (ncRNA). Simple random sampling over genotypes predicts the distribution of properties such as the mutational robustness or the number of stems per secondary structure found in naturally occurring ncRNA with surprising accuracy. Since we ignore natural selection, this strikingly close correspondence with the mapping suggests that structures allowing for functionality are easily discovered, despite the enormous size of the genetic spaces. The mapping is extremely biased: the majority of genotypes map to an exponentially small portion of the morphospace of all biophysically possible structures. Such strong constraints provide a non-adaptive explanation for the convergent evolution of structures such as the hammerhead ribozyme. These results presents a particularly clear example of bias in the arrival of variation strongly shaping evolutionary outcomes and may be relevant to Mayr's distinction between proximate and ultimate causes in evolutionary biology.

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

Correlation of Automorphism Group Size and Topological Properties with Program-size Complexity Evaluations of Graphs and Complex Networks

We show that numerical approximations of Kolmogorov complexity (K) applied to graph adjacency matrices capture some group-theoretic and topological properties of graphs and empirical networks ranging from metabolic to social networks. That K and the size of the group of automorphisms of a graph are correlated opens up interesting connections to problems in computational geometry, and thus connects several measures and concepts from complexity science. We show that approximations of K characterise synthetic and natural networks by their generating mechanisms, assigning lower algorithmic randomness to complex network models (Watts-Strogatz and Barabasi-Albert networks) and high Kolmogorov complexity to (random) Erdos-Renyi graphs. We derive these results via two different Kolmogorov complexity approximation methods applied to the adjacency matrices of the graphs and networks. The methods used are the traditional lossless compression approach to Kolmogorov complexity, and a normalised version of a Block Decomposition Method (BDM) measure, based on algorithmic probability theory.

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

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

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

A tractable genotype-phenotype map for the self-assembly of protein quaternary structure

The mapping between biological genotypes and phenotypes is central to the study of biological evolution. Here we introduce a rich, intuitive, and biologically realistic genotype-phenotype (GP) map, that serves as a model of self-assembling biological structures, such as protein complexes, and remains computationally and analytically tractable. Our GP map arises naturally from the self-assembly of polyomino structures on a 2D lattice and exhibits a number of properties: $\textit{redundancy}$ (genotypes vastly outnumber phenotypes), $\textit{phenotype bias}$ (genotypic redundancy varies greatly between phenotypes), $\textit{genotype component disconnectivity}$ (phenotypes consist of disconnected mutational networks) and $\textit{shape space covering}$ (most phenotypes can be reached in a small number of mutations). We also show that the mutational robustness of phenotypes scales very roughly logarithmically with phenotype redundancy and is positively correlated with phenotypic evolvability. Although our GP map describes the assembly of disconnected objects, it shares many properties with other popular GP maps for connected units, such as models for RNA secondary structure or the HP lattice model for protein tertiary structure. The remarkable fact that these important properties similarly emerge from such different models suggests the possibility that universal features underlie a much wider class of biologically realistic GP maps.

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

Viscous fingering at ultralow interfacial tension

We experimentally study the viscous fingering instability in a fluid-fluid phase separated colloid-polymer mixture by means of laser scanning confocal microscopy and microfluidics. We focus on three aspects of the instability. (i) The interface between the two demixed phases has an ultralow surface tension, such that we can address the role of thermal interface fluctuations. (ii) We image the interface in three dimensions allowing us to study the interplay between interface curvature and flow. (iii) The displacing fluid wets all walls completely, in contrast to traditional viscous fingering experiments, in which the displaced fluid wets the walls. We also perform lattice Boltzmann simulations, which help to interpret the experimental observations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Zigzag transitions and nonequilibrium pattern formation in colloidal chains

Paramagnetic colloidal particles that are optically trapped in a linear array can form a zigzag pattern when an external magnetic field induces repulsive interparticle interactions. When the traps are abruptly turned off, the particles form a nonequilibrium expanding pattern with a zigzag symmetry, even when the strength of the magnetic interaction is weaker than that required to break the linear symmetry of the equilibrium state. We show that the transition to the equilibrium zigzag state is always potentially possible for purely harmonic traps. For anharmonic traps that have a finite height, the equilibrium zigzag state becomes unstable above a critical anharmonicity. A normal mode analysis of the equilibrium line configuration demonstrates that increasing the magnetic field leads to a hardening and softening of the spring constants in the longitudinal and transverse directions, respectively. The mode that first becomes unstable is the mode with the zigzag symmetry, which explains the symmetry of nonequilibrium patterns. Our analytically tractable models help to give further insight into the way that the interplay of such factors as the length of the chain, hydrodynamic interactions, thermal fluctuations affect the formation and evolution of the experimentally observed nonequilibrium patterns.

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

Effect of bending rigidity on the knotting of a polymer under tension

A coarse-grained computational model is used to investigate how the bending rigidity of a polymer under tension affects the formation of a trefoil knot. Thermodynamic integration techniques are applied to demonstrate that the free-energy cost of forming a knot has a minimum at non-zero bending rigidity. The position of the minimum exhibits a power-law dependence on the applied tension. For knotted polymers with non-uniform bending rigidity, the knots preferentially localize in the region with a bending rigidity that minimizes the free-energy.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Extracting short-ranged interactions from structure factors

Inverting scattering experiments to obtain effective interparticle interactions is generally a poorly conditioned problem. L. Reatto (Phil. Mag. A 58, 37 (1986)) showed that for atomic liquids close to the triple point, inversions are hard because the structure closely resembles that of an equivalent hard-sphere fluid. Here I demonstrate that at low concentrations and for particles with short-ranged attractive potentials, S(k) also exhibits a very weak dependence on potential shape. Instead, different potentials can generate an S(k) that closely resembles that of the Baxter model with a similar second-virial coefficient. By contrast, in this energetic fluid regime, the inversion of an attractive interaction from real-space correlations such as the radial distribution function g(r) is well conditioned. Nevertheless, one may extract further information from S(k) by measuring isosbestic points, values of k where the scattering intensity I(k) or the structure factor S(k) is invariant to changes in interaction-potential well-depth. These points suggest a new extended corresponding states principle for particles in solution based on the packing fraction, the second osmotic virial coefficient, and a new measure of effective potential range.

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

Epistasis can lead to fragmented neutral spaces and contingency in evolution

In evolution, the effects of a single deleterious mutation can sometimes be compensated for by a second mutation which recovers the original phenotype. Such epistatic interactions have implications for the structure of genome space - namely, that networks of genomes encoding the same phenotype may not be connected by single mutational moves. We use the folding of RNA sequences into secondary structures as a model genotype-phenotype map and explore the neutral spaces corresponding to networks of genotypes with the same phenotype. In most of these networks, we find that it is not possible to connect all genotypes to one another by single point mutations. Instead, a network for a phenotypic structure with $n$ bonds typically fragments into at least $2^n$ neutral components, often of similar size. While components of the same network generate the same phenotype, they show important variations in their properties, most strikingly in their evolvability and mutational robustness. This heterogeneity implies contingency in the evolutionary process.

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

Pattern formation in colloidal explosions

We study the non-equilibrium pattern formation that emerges when magnetically repelling colloids, trapped by optical tweezers, are abruptly released, forming colloidal explosions. For multiple colloids in a single trap we observe a pattern of expanding concentric rings. For colloids individually trapped in a line, we observe explosions with a zigzag pattern that persists even when magnetic interactions are much weaker than those that break the linear symmetry in equilibrium. Theory and computer simulations quantitatively describe these phenomena both in and out of equilibrium. An analysis of the mode spectrum allows us to accurately quantify the non-harmonic nature of the optical traps. Colloidal explosions provide a new way to generate well-characterized non-equilibrium behaviour in colloidal systems.

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
preprint2008arXiv

Hydrodynamics of confined colloidal fluids in two dimensions

We apply a hybrid Molecular Dynamics and mesoscopic simulation technique to study the dynamics of two dimensional colloidal discs in confined geometries. We calculate the velocity autocorrelation functions, and observe the predicted $t^{-1}$ long time hydrodynamic tail that characterizes unconfined fluids, as well as more complex oscillating behavior and negative tails for strongly confined geometries. Because the $t^{-1}$ tail of the velocity autocorrelation function is cut off for longer times in finite systems, the related diffusion coefficient does not diverge, but instead depends logarithmically on the overall size of the system.

0 citations0 reviewsNo code linkedOpen access