Inverse design of a pyrochlore lattice of DNA origami through model-driven experiments.

Sophisticated statistical mechanics approaches and human intuition have demonstrated the possibility of self-assembling complex lattices or finite-size constructs. However, attempts so far have mostly only been successful in silico and often fail in experiment because of unpredicted traps associated with kinetic slowing down (gelation, glass transition) and competing ordered structures. Theoretical predictions also face the difficulty of encoding the desired interparticle interaction potential with the experimentally available nano- and micrometer-sized particles. To overcome these issues, we combine SAT assembly (a patchy-particle interaction design algorithm based on constrained optimization) with coarse-grained simulations of DNA nanotechnology to experimentally realize trap-free self-assembly pathways. We use this approach to assemble a pyrochlore three-dimensional lattice, coveted for its promise in the construction of optical metamaterials, and characterize it with small-angle x-ray scattering and scanning electron microscopy visualization.

Diffusion, viscosity, and linear rheology of valence-limited disordered fluids.

We numerically investigate the dynamics and linear rheology of disordered systems made of patchy particles, focusing on the role of valence, temperature, and bonding mechanism. We demonstrate that the dynamics is enslaved to bonding, giving rise to an activated behavior at low temperatures. By independently computing the diffusion constant and the viscosity from the simulations, we also confirm the validity of the Stokes-Einstein relation in valence-limited systems, with two caveats: (i) the diffusion constant requires a finite-size correction, at least at the intermediate density we investigate, and (ii) there is the onset of a breakdown that appears at the lowest temperatures considered. Finally, our results show that the storage and loss moduli of mixtures of divalent and M-valent particles exhibit an apparent power-law dependence on frequency, hinting at the possibility of using the composition to finely tune the rheological response of these materials. Our results compare well with literature experimental data on valence-limited DNA nanostars. In addition, the wealth of data we present and analyze here will help develop and test theoretical frameworks aimed at describing the dynamics of flexible limited-valence particles that self-assemble into disordered networks.

The oxDNA Coarse-Grained Model as a Tool to Simulate DNA Origami.

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

Structure and elasticity of model disordered, polydisperse, and defect-free polymer networks.

The elasticity of disordered and polydisperse polymer networks is a fundamental problem of soft matter physics that is still open. Here, we self-assemble polymer networks via simulations of a mixture of bivalent and tri- or tetravalent patchy particles, which result in an exponential strand length distribution analogous to that of experimental randomly cross-linked systems. After assembly, the network connectivity and topology are frozen and the resulting system is characterized. We find that the fractal structure of the network depends on the number density at which the assembly has been carried out, but that systems with the same mean valence and same assembly density have the same structural properties. Moreover, we compute the long-time limit of the mean-squared displacement, also known as the (squared) localization length, of the cross-links and of the middle monomers of the strands, showing that the dynamics of long strands is well described by the tube model. Finally, we find a relation connecting these two localization lengths at high density and connect the cross-link localization length to the shear modulus of the system.

A simple solution to the problem of self-assembling cubic diamond crystals.

The self-assembly of colloidal diamond (CD) crystals is considered as one of the most coveted goals of nanotechnology, both from the technological and fundamental points of view. For applications, colloidal diamond is a photonic crystal which can open new possibilities of manipulating light for information processing. From a fundamental point of view, its unique symmetry exacerbates a series of problems that are commonly faced during the self-assembly of target structures, such as the presence of kinetic traps and the formation of crystalline defects and alternative structures (polymorphs). Here we demonstrate that all these problems can be systematically addressed via SAT-assembly, a design framework that converts self-assembly into a Boolean satisfiability problem (SAT). Contrary to previous solutions (requiring four or more components), we prove that the assembly of the CD crystal only requires a binary mixture. Moreover, we use molecular dynamics simulations of a system composed by nearly a million nucleotides to test a DNA nanotechnology design that constitutes a promising candidate for experimental realization.

Designing Enhanced Entropy Binding in Single-Chain Nanoparticles.

Single-chain nanoparticles (SCNPs) are a new class of bio- and soft-matter polymeric objects in which a fraction of the monomers are able to form equivalently intra- or interpolymer bonds. Here we numerically show that a fully entropic gas-liquid phase separation can take place in SCNP systems. Control over the discontinuous (first-order) change-from a phase of independent diluted (fully-bonded) polymers to a phase in which polymers entropically bind to each other to form a (fully-bonded) polymer network-can be achieved by a judicious design of the patterns of reactive monomers along the polymer chain. Such a sensitivity arises from a delicate balance between the distinct entropic contributions controlling the binding.

SAT-assembly: a new approach for designing self-assembling systems.

We propose a general framework for solving inverse self-assembly problems, i.e. designing interactions between elementary units such that they assemble spontaneously into a predetermined structure. Our approach uses patchy particles as building blocks, where the different units bind at specific interaction sites (the patches), and we exploit the possibility of having mixtures with several components. The interaction rules between the patches is determined by transforming the combinatorial problem into a Boolean satisfiability problem (SAT) which searches for solutions where all bonds are formed in the target structure. Additional conditions, such as the non-satisfiability of competing structures (e.g. metastable states) can be imposed, allowing to effectively design the assembly path in order to avoid kinetic traps. We demonstrate this approach by designing and numerically simulating a cubic diamond structure from four particle species that assembles without competition from other polymorphs, including the hexagonal structure.

Hydrodynamic instability and flow reduction in polymer brush coated channels.

A polymer brush is a passive medium. At equilibrium the knowledge of its chemical composition and thickness is enough for a full system characterization. However, when the brush is exposed to fluid flow it reveals a much more intriguing nature, in which filamentous protrusions and the way they interact among themselves and with the surrounding fluid are of outmost importance. Here we investigate such a rich behavior via numerical simulations. We focus on the brush hydrodynamic response at low Reynolds numbers, observing a significant fluid flow reduction inside a polymer-brush coated channel. We find that the reduction of the flow inside the channel is significantly larger than what would happen if the brush effect consisted only in reducing the effective channel width. This amplified reduction is understood as being due to the morphological instability of the brush-liquid interface which is shown to have an elastic origin: the mechanical stress acting on the brush due to the imposed flow is partially released by the interface modulation. In turn, this modulation dissipates more energy than a flat interface in the surrounding fluid, causing a reduction of flow velocity. Our results and interpretations provide an explanation for recent experimental measurements.

Sharing the Load: Stress Redistribution Governs Fracture of Polymer Double Networks.

The stress response of polymer double networks depends not only on the properties of the constituent networks but also on the interactions arising between them. Here, we demonstrate, via coarse-grained simulations, that both their global stress response and their microscopic fracture mechanics are governed by load sharing through these internetwork interactions. By comparing our results with affine predictions, where stress redistribution is by definition homogeneous, we show that stress redistribution is highly inhomogeneous. In particular, the affine prediction overestimates the fraction of broken chains by almost an order of magnitude. Furthermore, homogeneous stress distribution predicts a single fracture process, while in our simulations, fracture of sacrificial chains takes place in two steps governed by load sharing within a network and between networks, respectively. Our results thus provide a detailed microscopic picture of how inhomogeneous stress redistribution after rupture of chains governs the fracture of polymer double networks.

Two-step deswelling in the Volume Phase Transition of thermoresponsive microgels.

Thermoresponsive microgels are one of the most investigated types of soft colloids, thanks to their ability to undergo a Volume Phase Transition (VPT) close to ambient temperature. However, this fundamental phenomenon still lacks a detailed microscopic understanding, particularly regarding the presence and the role of charges in the deswelling process. This is particularly important for the widely used poly(N-isopropylacrylamide)-based microgels, where the constituent monomers are neutral but charged groups arise due to the initiator molecules used in the synthesis. Here, we address this point combining experiments with state-of-the-art simulations to show that the microgel collapse does not happen in a homogeneous fashion, but through a two-step mechanism, entirely attributable to electrostatic effects. The signature of this phenomenon is the emergence of a minimum in the ratio between gyration and hydrodynamic radii at the VPT. Thanks to simulations of microgels with different cross-linker concentrations, charge contents, and charge distributions, we provide evidence that peripheral charges arising from the synthesis are responsible for this behavior and we further build a universal master curve able to predict the two-step deswelling. Our results have direct relevance on fundamental soft condensed matter science and on applications where microgels are involved, ranging from materials to biomedical technologies.

Effect of Internal Architecture on the Assembly of Soft Particles at Fluid Interfaces.

Monolayers of soft colloidal particles confined at fluid interfaces are at the core of a broad range of technological processes, from the stabilization of responsive foams and emulsions to advanced lithographic techniques. However, establishing a fundamental relation between their internal architecture, which is controlled during synthesis, and their structural and mechanical properties upon interfacial confinement remains an elusive task. To address this open issue, which defines the monolayer's properties, we synthesize core-shell microgels, whose soft core can be chemically degraded in a controlled fashion. This strategy allows us to obtain a series of particles ranging from analogues of standard batch-synthesized microgels to completely hollow ones after total core removal. Combined experimental and numerical results show that our hollow particles have a thin and deformable shell, leading to a temperature-responsive collapse of the internal cavity and a complete flattening after adsorption at a fluid interface. Mechanical characterization shows that a critical degree of core removal is required to obtain soft disk-like particles at an oil-water interface, which present a distinct response to compression. At low packing fractions, the mechanical response of the monolayer is dominated by the outer polymer chains forming a corona surrounding the particles within the interfacial plane, regardless of the presence of a core. By contrast, at high compression, the absence of a core enables the particles to deform in the direction orthogonal to the interface and to be continuously compressed without altering the monolayer structure. These findings show how fine, single-particle architectural control during synthesis can be engineered to determine the interfacial behavior of microgels, enabling one to link particle conformation with the resulting material properties.

Effect of Chain Polydispersity on the Elasticity of Disordered Polymer Networks.

Due to their unique structural and mechanical properties, randomly cross-linked polymer networks play an important role in many different fields, ranging from cellular biology to industrial processes. In order to elucidate how these properties are controlled by the physical details of the network (e.g., chain-length and end-to-end distributions), we generate disordered phantom networks with different cross-linker concentrations C and initial densities ρinit and evaluate their elastic properties. We find that the shear modulus computed at the same strand concentration for networks with the same C, which determines the number of chains and the chain-length distribution, depends strongly on the preparation protocol of the network, here controlled by ρinit. We rationalize this dependence by employing a generic stress-strain relation for polymer networks that does not rely on the specific form of the polymer end-to-end distance distribution. We find that the shear modulus of the networks is a nonmonotonic function of the density of elastically active strands, and that this behavior has a purely entropic origin. Our results show that if short chains are abundant, as it is always the case for randomly cross-linked polymer networks, the knowledge of the exact chain conformation distribution is essential for correctly predicting the elastic properties. Finally, we apply our theoretical approach to literature experimental data, qualitatively confirming our interpretations.

OxDNA.org: a public webserver for coarse-grained simulations of DNA and RNA nanostructures.

OxDNA and oxRNA are popular coarse-grained models used by the DNA/RNA nanotechnology community to prototype, analyze and rationalize designed DNA and RNA nanostructures. Here, we present oxDNA.org, a graphical web interface for running, visualizing and analyzing oxDNA and oxRNA molecular dynamics simulations on a GPU-enabled high performance computing server. OxDNA.org automatically generates simulation files, including a multi-step relaxation protocol for structures exported in non-physical states from DNA/RNA design tools. Once the simulation is complete, oxDNA.org provides an interactive visualization and analysis interface using the browser-based visualizer oxView to facilitate the understanding of simulation results for a user's specific structure. This online tool significantly lowers the entry barrier of integrating simulations in the nanostructure design pipeline for users who are not experts in the technical aspects of molecular simulation. The webserver is freely available at oxdna.org.

Charge affinity and solvent effects in numerical simulations of ionic microgels.

Ionic microgel particles are intriguing systems in which the properties of thermo-responsive polymeric colloids are enriched by the presence of charged groups. In order to rationalize their properties and predict the behaviour of microgel suspensions, it is necessary to develop a coarse-graining strategy that starts from the accurate modelling of single particles. Here, we provide a numerical advancement of a recently-introduced model for charged co-polymerized microgels by improving the treatment of ionic groups in the polymer network. We investigate the thermoresponsive properties of the particles, in particular their swelling behaviour and structure, finding that, when charged groups are considered to be hydrophilic at all temperatures, highly charged microgels do not achieve a fully collapsed state, in favorable comparison to experiments. In addition, we explicitly include the solvent in the description and put forward a mapping between the solvophobic potential in the absence of the solvent and the monomer-solvent interactions in its presence, which is found to work very accurately for any charge fraction of the microgel. Our work paves the way for comparing single-particle properties and swelling behaviour of ionic microgels to experiments and to tackle the study of these charged soft particles at a liquid-liquid interface.

Shape control of soft patchy nanoparticles under confinement.

Molecular building blocks undergoing a hierarchical assembly process form nano-scale objects which can further assemble into supramolecular structures. When the intermediate units have a limited valence in bonding, complex structures with tailored properties can be created. Here, we consider a composite, star-shaped particle made of f diblock copolymer chains uniformly grafted on a spherical colloid and investigate its first self-assembly stage both in the bulk and under lateral confinement. By means of numerical simulations, we show that, in the bulk, this system develops aggregates whose number and size depend on the temperature as well as on the relative ratio of solvophobic monomers. The emerging aggregates are referred to as patches and impart directionality in bonding to the complex particle. We further characterize how we can control, by changing the lateral confinement, the shape of the brush and the patch properties as a function of the distance between the confining walls. We find that the number of the patches can be determined by tuning the degree of confinement imposed on the particle. Finally, we employ a continuum mechanics model, known as the Liquid Drop Model, to gain insight into the elastic properties of the system. This theoretical approach allows to connect the patch properties to the elastic response of the composite particle.

Designer protein assemblies with tunable phase diagrams in living cells.

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

Modeling Microgels with a Controlled Structure across the Volume Phase Transition.

Thermoresponsive microgels are soft colloids that find widespread use as model systems for soft matter physics. Their complex internal architecture, made of a disordered and heterogeneous polymer network, has been so far a major challenge for computer simulations. In this work, we put forward a coarse-grained model of microgels whose structural properties are in quantitative agreement with results obtained with small-angle X-ray scattering experiments across a wide range of temperatures, encompassing the volume phase transition. These results bridge the gap between experiments and simulations of individual microgel particles, paving the way to theoretically address open questions about their bulk properties with unprecedented nano- and microscale resolution.

Numerical insights on ionic microgels: structure and swelling behaviour.

Recent progress has been made in the numerical modelling of neutral microgel particles with a realistic, disordered structure. In this work we extend this approach to the case of co-polymerised microgels where a thermoresponsive polymer is mixed with acidic groups. We compare the cases where counterions directly interact with microgel charges or are modelled implicitly through a Debye-Hückel description. We do so by performing extensive numerical simulations of single microgels across the volume phase transition (VPT) varying the temperature and the fraction of charged monomers. We find that the presence of charges considerably alters the microgel structure, quantified by the monomer density profiles and by the form factors of the microgels, particularly close to the VPT. We observe significant deviations between the implicit and explicit models, with the latter comparing more favourably to available experiments. In particular, we observe a shift of the VPT temperature to larger values as the amount of charged monomers increases. We also find that below the VPT the microgel-counterion complex is almost neutral, while it develops a net charge above the VPT. Interestingly, under these conditions the collapsed microgel still retains a large amount of counterions inside its structure. Since these interesting features cannot be captured by the implicit model, our results show that it is crucial to explicitly include the counterions in order to realistically model ionic thermoresponsive microgels.

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

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

q-Independent Slow Dynamics in Atomic and Molecular Systems.

Investigating million-atom systems for very long simulation times, we demonstrate that the collective density-density correlation time (τ_{α}) in simulated supercooled water and silica becomes wave-vector independent (q^{0}) when the probing wavelength is several times larger than the interparticle distance. The q independence of the collective density-density correlation functions, a feature clearly observed in light-scattering studies of some soft-matter systems, is thus a genuine feature of many (but not all) slow-dynamics systems, either atomic, molecular, or colloidal. Indeed, we show that when the dynamics of the density fluctuations includes particle-type diffusion, as in the case of the Lennard-Jones binary-mixture model, the q^{0} regime does not set in and the relaxation time continues to scale as τ_{α}∼q^{-2} even at small q.