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

We use Monte Carlo simulations and freeenergy techniques to show that binary
solutions of penta and hexavalent twodimensional 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 coarsegrained representation of DNA multiarm
`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 coarsegrained model of DNA.

We introduce an extended version of oxDNA, a coarsegrained model of DNA
designed to capture the thermodynamic, structural and mechanical properties of
single and doublestranded DNA. By including explicit major and minor grooves,
and by slightly modifying the coaxial stacking and backbonebackbone
interactions, we improve the ability of the model to treat large
(kilobasepair) 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 singlestranded 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 saltdependent properties of
DNA.

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,
twodimensional arrays, and ribbons, is a bulged duplex, that is two helical
segments that connect at a bulge loop. We use a coarsegrained 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 singlestranded
DNA is significant and the loop resists compression. The properties of this
basic structural motif clearly correlate with the structural behavior of
certain nanoscale objects, where the enhanced flexibility associated with
larger bulges has been used to tune the selfassembly product as well as the
detailed geometry of the resulting nanostructures.

The effect of secondary structure on DNA duplex formation is poorly
understood. We use a coarsegrained model of DNA to show that specific 3 and
4base 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 speedup
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.

Protein aggregation is an important field of investigation because it is
closely related to the problem of neurodegenerative diseases, to the
development of biomaterials, and to the growth of cellular structures such as
cytoskeleton. Selfaggregation of protein amyloids, for example, is a
complicated process involving many species and levels of structures. This
complexity, however, can be dealt with using statistical mechanical tools, such
as free energies, partition functions, and transfer matrices. In this article,
we review general strategies for studying protein aggregation using statistical
mechanical approaches and show that canonical and grand canonical ensembles can
be used in such approaches. The grand canonical approach is particularly
convenient since competing pathways of assembly and disassembly can be
considered simultaneously. Another advantage of using statistical mechanics is
that numerically exact solutions can be obtained for all of the thermodynamic
properties of fibrils, such as the amount of fibrils formed, as a function of
initial protein concentration. Furthermore, statistical mechanics models can be
used to fit experimental data when they are available for comparison.

We propose a kinetic model for the selfaggregation by amyloid proteins. By
extending several wellknown models for protein aggregation, the time evolution
of aggregate concentrations containing $r$ proteins, denoted $c_r(t)$, can be
written in terms of generalized Smoluchowski kinetics. With this approach we
take into account all possible aggregation and fragmentation reactions
involving clusters of any size. Correspondingly, an aggregate of size x+y could
be formed by or breakup into two smaller constituent aggregates of sizes x and
y. The rates of each aggregation or fragmentation reaction, called kernels, are
specified in terms of the aggregate size, and we solve $c_r(t)$ for large
cluster sizes using numerical techniques. We show that by using Smoluchowski
kinetics many pathways to fibrillation are possible and quantities, such as the
aggregate length distribution at an arbitrary time, can be calculated. We show
that the predicted results of the model are in agreement with the experimental
observations.

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

We develop a theory of aggregation using statistical mechanical methods. An
example of a complicated aggregation system with several levels of structures
is peptide/protein selfassembly. The problem of protein aggregation is
important for the understanding and treatment of neurodegenerative diseases and
also for the development of biomacromolecules as new materials. We write the
effective Hamiltonian in terms of interaction energies between protein
monomers, protein and solvent, as well as between protein filaments. The grand
partition function can be expressed in terms of a ZimmBragglike transfer
matrix, which is calculated exactly and all thermodynamic properties can be
obtained. We start with twostate and threestate descriptions of protein
monomers using Potts models that can be generalized to include qstates, for
which the exactly solvable feature of the model remains. We focus on n X N
lattice systems, corresponding to the ordered structures observed in some real
fibrils. We have obtained results on nucleation processes and phase diagrams,
in which a protein property such as the sheet content of aggregates is
expressed as a function of the number of proteins on the lattice and
interprotein or interfacial interaction energies. We have applied our methods
to A{\beta}(140) and Curli fibrils and obtained results in good agreement with
experiments.

In view of the important role helixsheet transitions play in protein
aggregation, we introduce a simple model to study secondary structural
transitions of helixcoilsheet systems using a Potts model starting with an
effective Hamiltonian. This energy function depends on four parameters that
approximately describe entropic and enthalpic contributions to the stability of
a polypeptide in helical and sheet conformations. The sheet structures involve
longrange interactions between residues which are far in sequence, but are in
contact in real space. Such contacts are included in the Hamiltonian. Using
standard statistical mechanical techniques, the partition function is solved
exactly using transfer matrices. Based on this model, we study thermodynamic
properties of polypeptides, including phase transitions between helix, sheet,
and coil structures.