Aromatic and arginine content drives multiphasic condensation of protein-RNA mixtures.

Multiphasic architectures are found ubiquitously in biomolecular condensates and are thought to have important implications for the organization of multiple chemical reactions within the same compartment. Many of these multiphasic condensates contain RNA in addition to proteins. Here, we investigate the importance of different interactions in multiphasic condensates comprising two different proteins and RNA using computer simulations with a residue-resolution coarse-grained model of proteins and RNA. We find that in multilayered condensates containing RNA in both phases, protein-RNA interactions dominate, with aromatic residues and arginine forming the key stabilizing interactions. The total aromatic and arginine content of the two proteins must be appreciably different for distinct phases to form, and we show that this difference increases as the system is driven toward greater multiphasicity. Using the trends observed in the different interaction energies of this system, we demonstrate that we can also construct multilayered condensates with RNA preferentially concentrated in one phase. The "rules" identified can thus enable the design of synthetic multiphasic condensates to facilitate further study of their organization and function.

Phase diagrams-Why they matter and how to predict them.

Understanding the thermodynamic stability and metastability of materials can help us to, for example, gauge whether crystalline polymorphs in pharmaceutical formulations are likely to be durable. It can also help us to design experimental routes to novel phases with potentially interesting properties. In this Perspective, we provide an overview of how thermodynamic phase behavior can be quantified both in computer simulations and machine-learning approaches to determine phase diagrams, as well as combinations of the two. We review the basic workflow of free-energy computations for condensed phases, including some practical implementation advice, ranging from the Frenkel-Ladd approach to thermodynamic integration and to direct-coexistence simulations. We illustrate the applications of such methods on a range of systems from materials chemistry to biological phase separation. Finally, we outline some challenges, questions, and practical applications of phase-diagram determination which we believe are likely to be possible to address in the near future using such state-of-the-art free-energy calculations, which may provide fundamental insight into separation processes using multicomponent solvents.

A streamlined molecular-dynamics workflow for computing solubilities of molecular and ionic crystals.

Computing the solubility of crystals in a solvent using atomistic simulations is notoriously challenging due to the complexities and convergence issues associated with free-energy methods, as well as the slow equilibration in direct-coexistence simulations. This paper introduces a molecular-dynamics workflow that simplifies and robustly computes the solubility of molecular or ionic crystals. This method is considerably more straightforward than the state-of-the-art, as we have streamlined and optimised each step of the process. Specifically, we calculate the chemical potential of the crystal using the gas-phase molecule as a reference state, and employ the S0 method to determine the concentration dependence of the chemical potential of the solute. We use this workflow to predict the solubilities of sodium chloride in water, urea polymorphs in water, and paracetamol polymorphs in both water and ethanol. Our findings indicate that the predicted solubility is sensitive to the chosen potential energy surface. Furthermore, we note that the harmonic approximation often fails for both molecular crystals and gas molecules at or above room temperature, and that the assumption of an ideal solution becomes less valid for highly soluble substances.

Simulations of DNA-Origami Self-Assembly Reveal Design-Dependent Nucleation Barriers.

Nucleation is the rate-determining step in the kinetics of many self-assembly processes. However, the importance of nucleation in the kinetics of DNA-origami self-assembly, which involves both the binding of staple strands and the folding of the scaffold strand, is unclear. Here, using Monte Carlo simulations of a lattice model of DNA origami, we find that some, but not all, designs can have a nucleation barrier and that this barrier disappears at lower temperatures, rationalizing the success of isothermal assembly. We show that the height of the nucleation barrier depends primarily on the coaxial stacking of staples that are adjacent on the same helix, a parameter that can be modified with staple design. Creating a nucleation barrier to DNA-origami assembly could be useful in optimizing assembly times and yields, while eliminating the barrier may allow for fast molecular sensors that can assemble/disassemble without hysteresis in response to changes in the environment.

Targeted modulation of protein liquid-liquid phase separation by evolution of amino-acid sequence.

Rationally and efficiently modifying the amino-acid sequence of proteins to control their ability to undergo liquid-liquid phase separation (LLPS) on demand is not only highly desirable, but can also help to elucidate which protein features are important for LLPS. Here, we propose a computational method that couples a genetic algorithm to a sequence-dependent coarse-grained protein model to evolve the amino-acid sequences of phase-separating intrinsically disordered protein regions (IDRs), and purposely enhance or inhibit their capacity to phase-separate. We validate the predicted critical solution temperatures of the mutated sequences with ABSINTH, a more accurate all-atom model. We apply the algorithm to the phase-separating IDRs of three naturally occurring proteins, namely FUS, hnRNPA1 and LAF1, as prototypes of regions that exist in cells and undergo homotypic LLPS driven by different types of intermolecular interaction, and we find that the evolution of amino-acid sequences towards enhanced LLPS is driven in these three cases, among other factors, by an increase in the average size of the amino acids. However, the direction of change in the molecular driving forces that enhance LLPS (such as hydrophobicity, aromaticity and charge) depends on the initial amino-acid sequence. Finally, we show that the evolution of amino-acid sequences to modulate LLPS is strongly coupled to the make-up of the medium (e.g. the presence or absence of RNA), which may have significant implications for our understanding of phase separation within the many-component mixtures of biological systems.

Direct observation and rational design of nucleation behavior in addressable self-assembly.

To optimize a self-assembly reaction, it is essential to understand the factors that govern its pathway. Here, we examine the influence of nucleation pathways in a model system for addressable, multicomponent self-assembly based on a prototypical "DNA-brick" structure. By combining temperature-dependent dynamic light scattering and atomic force microscopy with coarse-grained simulations, we show how subtle changes in the nucleation pathway profoundly affect the yield of the correctly formed structures. In particular, we can increase the range of conditions over which self-assembly occurs by using stable multisubunit clusters that lower the nucleation barrier for assembling subunits in the interior of the structure. Consequently, modifying only a small portion of a structure is sufficient to optimize its assembly. Due to the generality of our coarse-grained model and the excellent agreement that we find with our experimental results, the design principles reported here are likely to apply generically to addressable, multicomponent self-assembly.

Predicting the phase diagram of titanium dioxide with random search and pattern recognition.

Predicting phase stabilities of crystal polymorphs is central to computational materials science and chemistry. Such predictions are challenging because they first require searching for potential energy minima and then performing arduous free-energy calculations to account for entropic effects at finite temperatures. Here, we develop a framework that facilitates such predictions by exploiting all the information obtained from random searches of crystal structures. This framework combines automated clustering, classification and visualisation of crystal structures with machine-learning estimation of their enthalpy and entropy. We demonstrate the framework on the technologically important system of TiO2, which has many polymorphs, without relying on prior knowledge of known phases. We find a number of new phases and predict the phase diagram and metastabilities of crystal polymorphs at 1600 K, benchmarking the results against full free-energy calculations.

Numerical calculation of free-energy barriers for entangled polymer nucleation.

The crystallization of entangled polymers from their melt is investigated using computer simulation with a coarse-grained model. Using hybrid Monte Carlo simulations enables us to probe the behavior of long polymer chains. We identify solid-like beads with a centrosymmetry local order parameter and compute the nucleation free-energy barrier at relatively high supercooling with adaptive-bias windowed umbrella sampling. Our results demonstrate that the critical nucleus sizes and the heights of free-energy barriers do not significantly depend on the molecular weight of the polymer; however, the nucleation rate decreases with the increase in molecular weight. Moreover, an analysis of the composition of the critical nucleus suggests that intra-molecular growth of the nucleated cluster does not contribute significantly to crystallization for this system.

Programming the Nucleation of DNA Brick Self-Assembly with a Seeding Strand.

Recently, the DNA brick strategy has provided a highly modular and scalable approach for the construction of complex structures, which can be used as nanoscale pegboards for the precise organization of molecules and nanoparticles for many applications. Despite the dramatic increase of structural complexity provided by the DNA brick method, the assembly pathways are still poorly understood. Herein, we introduce a "seed" strand to control the crucial nucleation and assembly pathway in DNA brick assembly. Through experimental studies and computer simulations, we successfully demonstrate that the regulation of the assembly pathways through seeded growth can accelerate the assembly kinetics and increase the optimal temperature by circa 4-7 °C for isothermal assembly. By improving our understanding of the assembly pathways, we provide new guidelines for the design of programmable pathways to improve the self-assembly of DNA nanostructures.

Microscopic analysis of thermo-orientation in systems of off-centre Lennard-Jones particles.

When fluids of anisotropic molecules are placed in temperature gradients, the molecules may align themselves along the gradient: this is called thermo-orientation. We discuss the theory of this effect in a fluid of particles that interact by a spherically symmetric potential, where the particles' centres of mass do not coincide with their interaction centres. Starting from the equations of motion of the molecules, we show how a simple assumption of local equipartition of energy can be used to predict the thermo-orientation effect, recovering the result of Wirnsberger et al. [Phys. Rev. Lett. 120, 226001 (2018)]. Within this approach, we show that for particles with a single interaction centre, the thermal centre of the molecule must coincide with the interaction centre. The theory also explains the coupling between orientation and kinetic energy that is associated with this non-Boltzmann distribution. We discuss deviations from this local equipartition assumption, showing that these can occur in linear response to a temperature gradient. We also present numerical simulations showing significant deviations from the local equipartition predictions, which increase as the centre of mass of the molecule is displaced further from its interaction centre.

Lattice models and Monte Carlo methods for simulating DNA origami self-assembly.

The optimal design of DNA origami systems that assemble rapidly and robustly is hampered by the lack of a model for self-assembly that is sufficiently detailed yet computationally tractable. Here, we propose a model for DNA origami that strikes a balance between these two criteria by representing these systems on a lattice at the level of binding domains. The free energy of hybridization between individual binding domains is estimated with a nearest-neighbour model. Double helical segments are treated as being rigid, but we allow flexibility at points where the backbone of one of the strands is interrupted, which provides a reasonably realistic representation of partially and fully assembled states. Particular attention is paid to the constraints imposed by the double helical twist, as they determine where strand crossovers between adjacent helices can occur. To improve the efficiency of sampling configuration space, we develop Monte Carlo methods for sampling scaffold conformations in near-assembled states, and we carry out simulations in the grand canonical ensemble, enabling us to avoid considering states with unbound staples. We demonstrate that our model can quickly sample assembled configurations of a small origami design previously studied with the oxDNA model, as well as a design with staples that span longer segments of the scaffold. The sampling ability of our method should allow for good statistics to be obtained when studying the assembly pathways and is suited to investigating, in particular, the effects of design and assembly conditions on these pathways and their resulting final assembled structures.

Rational design of self-assembly pathways for complex multicomponent structures.

The field of complex self-assembly is moving toward the design of multiparticle structures consisting of thousands of distinct building blocks. To exploit the potential benefits of structures with such "addressable complexity," we need to understand the factors that optimize the yield and the kinetics of self-assembly. Here we use a simple theoretical method to explain the key features responsible for the unexpected success of DNA-brick experiments, which are currently the only demonstration of reliable self-assembly with such a large number of components. Simulations confirm that our theory accurately predicts the narrow temperature window in which error-free assembly can occur. Even more strikingly, our theory predicts that correct assembly of the complete structure may require a time-dependent experimental protocol. Furthermore, we predict that low coordination numbers result in nonclassical nucleation behavior, which we find to be essential for achieving optimal nucleation kinetics under mild growth conditions. We also show that, rather surprisingly, the use of heterogeneous bond energies improves the nucleation kinetics and in fact appears to be necessary for assembling certain intricate 3D structures. This observation makes it possible to sculpt nucleation pathways by tuning the distribution of interaction strengths. These insights not only suggest how to improve the design of structures based on DNA bricks, but also point the way toward the creation of a much wider class of chemical or colloidal structures with addressable complexity.

Investigating the role of boundary bricks in DNA brick self-assembly.

In the standard DNA brick set-up, distinct 32-nucleotide strands of single-stranded DNA are each designed to bind specifically to four other such molecules. Experimentally, it has been demonstrated that the overall yield is increased if certain bricks which occur on the outer faces of target structures are merged with adjacent bricks. However, it is not well understood by what mechanism such 'boundary bricks' increase the yield, as they likely influence both the nucleation process and the final stability of the target structure. Here, we use Monte Carlo simulations with a patchy particle model of DNA bricks to investigate the role of boundary bricks in the self-assembly of complex multicomponent target structures. We demonstrate that boundary bricks lower the free-energy barrier to nucleation and that boundary bricks on edges stabilize the final structure. However, boundary bricks are also more prone to aggregation, as they can stabilize partially assembled intermediates. We explore some design strategies that permit us to benefit from the stabilizing role of boundary bricks whilst minimizing their ability to hinder assembly; in particular, we show that maximizing the total number of boundary bricks is not an optimal strategy.

Effects of co-ordination number on the nucleation behaviour in many-component self-assembly.

We report canonical and grand-canonical lattice Monte Carlo simulations of the self-assembly of addressable structures comprising hundreds of distinct component types. The nucleation behaviour, in the form of free-energy barriers to nucleation, changes significantly as the co-ordination number of the building blocks is changed from 4 to 8 to 12. Unlike tetrahedral structures - which roughly correspond to DNA bricks that have been studied in experiments - the shapes of the free-energy barriers of higher co-ordination structures depend strongly on the supersaturation, and such structures require a very significant driving force for structure growth before nucleation becomes thermally accessible. Although growth at high supersaturation results in more defects during self-assembly, we show that high co-ordination number structures can still be assembled successfully in computer simulations and that they exhibit self-assembly behaviour analogous to DNA bricks. In particular, the self-assembly remains modular, enabling in principle a wide variety of nanostructures to be assembled, with a greater spatial resolution than is possible in low co-ordination structures.

DNA brick self-assembly with an off-lattice potential.

We report Monte Carlo simulations of a simple off-lattice patchy-particle model for DNA 'bricks'. We relate the parameters that characterise this model with the binding free energy of pairs of single-stranded DNA molecules. We verify that an off-lattice potential parameterised in this way reproduces much of the behaviour seen with a simpler lattice model we introduced previously, although the relaxation of the geometric constraints leads to a more error-prone self-assembly pathway. We investigate the self-assembly process as a function of the strength of the non-specific interactions. We show that our off-lattice model for DNA bricks results in robust self-assembly into a variety of target structures.

Communication: theoretical prediction of free-energy landscapes for complex self-assembly.

We present a technique for calculating free-energy profiles for the nucleation of multicomponent structures that contain as many species as building blocks. We find that a key factor is the topology of the graph describing the connectivity of the target assembly. By considering the designed interactions separately from weaker, incidental interactions, our approach yields predictions for the equilibrium yield and nucleation barriers. These predictions are in good agreement with corresponding Monte Carlo simulations. We show that a few fundamental properties of the connectivity graph determine the most prominent features of the assembly thermodynamics. Surprisingly, we find that polydispersity in the strengths of the designed interactions stabilizes intermediate structures and can be used to sculpt the free-energy landscape for self-assembly. Finally, we demonstrate that weak incidental interactions can preclude assembly at equilibrium due to the combinatorial possibilities for incorrect association.

Numerical evidence for nucleated self-assembly of DNA brick structures.

The observation by Ke et al. [Science 338, 1177 (2012)] that large numbers of short, predesigned DNA strands can assemble into three-dimensional target structures came as a great surprise, as no colloidal self-assembling system has ever achieved the same degree of complexity. That failure seemed easy to rationalize: the larger the number of distinct building blocks, the higher the expected error rate for self-assembly. The experiments of Ke et al. have disproved this argument. Here, we report Monte Carlo simulations of the self-assembly of a DNA brick cube, comprising approximately 1000 types of DNA strand, using a simple model. We model the DNA strands as lattice tetrahedra with attractive patches, the interaction strengths of which are computed using a standard thermodynamic model. We find that, within a narrow temperature window, the target structure assembles with high probability. Our simulations suggest that misassembly is disfavored because of a slow nucleation step. As our model incorporates no aspect of DNA other than its binding properties, these simulations suggest that, with proper design of the building blocks, other systems, such as colloids, may also assemble into truly complex structures.

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.

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 stabilized over its crystalline approximants.

Note: 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.