Self-replication of Aß42 aggregates occurs on small and isolated fibril sites.

Self-replication of amyloid fibrils via secondary nucleation is an intriguing physicochemical phenomenon in which existing fibrils catalyze the formation of their own copies. The molecular events behind this fibril surface-mediated process remain largely inaccessible to current structural and imaging techniques. Using statistical mechanics, computer modeling, and chemical kinetics, we show that the catalytic structure of the fibril surface can be inferred from the aggregation behavior in the presence and absence of a fibril-binding inhibitor. We apply our approach to the case of Alzheimer's A[Formula: see text] amyloid fibrils formed in the presence of proSP-C Brichos inhibitors. We find that self-replication of A[Formula: see text] fibrils occurs on small catalytic sites on the fibril surface, which are far apart from each other, and each of which can be covered by a single Brichos inhibitor.

Temperature protocols to guide selective self-assembly of competing structures.

Multicomponent self-assembly mixtures offer the possibility of encoding multiple target structures with the same set of interacting components. Selective retrieval of one of the stored structures has been attempted by preparing an initial state that favors the assembly of the required target, through seeding, concentration patterning, or specific choices of interaction strengths. This may not be possible in an experiment where on-the-fly reconfiguration of the building blocks to switch functionality may be required. In this paper, we explore principles of inverse design of a multicomponent, self-assembly mixture capable of encoding two competing structures that can be selected through simple temperature protocols. We design the target structures to realize the generic situation in which one of the targets has the lower nucleation barrier, while the other is globally more stable. We observe that, to avoid the formation of spurious or chimeric aggregates, the number of neighboring component pairs that occur in both structures should be minimal. Our design also requires the inclusion of components that are part of only one of the target structures. We observe, however, that to maximize the selectivity of retrieval, the component library itself should be maximally shared by the two targets, within such a constraint. We demonstrate that temperature protocols can be designed that lead to the formation of either one of the target structures with high selectivity. We discuss the important role played by secondary aggregation products in improving selectivity, which we term "vestigial aggregates."

Liquid network connectivity regulates the stability and composition of biomolecular condensates with many components.

One of the key mechanisms used by cells to control the spatiotemporal organization of their many components is the formation and dissolution of biomolecular condensates through liquid-liquid phase separation (LLPS). Using a minimal coarse-grained model that allows us to simulate thousands of interacting multivalent proteins, we investigate the physical parameters dictating the stability and composition of multicomponent biomolecular condensates. We demonstrate that the molecular connectivity of the condensed-liquid network-i.e., the number of weak attractive protein-protein interactions per unit of volume-determines the stability (e.g., in temperature, pH, salt concentration) of multicomponent condensates, where stability is positively correlated with connectivity. While the connectivity of scaffolds (biomolecules essential for LLPS) dominates the phase landscape, introduction of clients (species recruited via scaffold-client interactions) fine-tunes it by transforming the scaffold-scaffold bond network. Whereas low-valency clients that compete for scaffold-scaffold binding sites decrease connectivity and stability, those that bind to alternate scaffold sites not required for LLPS or that have higher-than-scaffold valencies form additional scaffold-client-scaffold bridges increasing stability. Proteins that establish more connections (via increased valencies, promiscuous binding, and topologies that enable multivalent interactions) support the stability of and are enriched within multicomponent condensates. Importantly, proteins that increase the connectivity of multicomponent condensates have higher critical points as pure systems or, if pure LLPS is unfeasible, as binary scaffold-client mixtures. Hence, critical points of accessible systems (i.e., with just a few components) might serve as a unified thermodynamic parameter to predict the composition of multicomponent condensates.

Computational design of probes to detect bacterial genomes by multivalent binding.

Rapid methods for diagnosis of bacterial infections are urgently needed to reduce inappropriate use of antibiotics, which contributes to antimicrobial resistance. In many rapid diagnostic methods, DNA oligonucleotide probes, attached to a surface, bind to specific nucleotide sequences in the DNA of a target pathogen. Typically, each probe binds to a single target sequence; i.e., target-probe binding is monovalent. Here we show using computer simulations that the detection sensitivity and specificity can be improved by designing probes that bind multivalently to the entire length of the pathogen genomic DNA, such that a given probe binds to multiple sites along the target DNA. Our results suggest that multivalent targeting of long pieces of genomic DNA can allow highly sensitive and selective binding of the target DNA, even if competing DNA in the sample also contains binding sites for the same probe sequences. Our results are robust to mild fragmentation of the bacterial genome. Our conclusions may also be relevant for DNA detection in other fields, such as disease diagnostics more broadly, environmental management, and food safety.

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.

Using sequence data to predict the self-assembly of supramolecular collagen structures.

Collagen fibrils are the major constituents of the extracellular matrix, which provides structural support to vertebrate connective tissues. It is widely assumed that the superstructure of collagen fibrils is encoded in the primary sequences of the molecular building blocks. However, the interplay between large-scale architecture and small-scale molecular interactions makes the ab initio prediction of collagen structure challenging. Here, we propose a model that allows us to predict the periodic structure of collagen fibers and the axial offset between the molecules, purely on the basis of simple predictive rules for the interaction between amino acid residues. With our model, we identify the sequence-dependent collagen fiber geometries with the lowest free energy and validate the predicted geometries against the available experimental data. We propose a procedure for searching for optimal staggering distances. Finally, we build a classification algorithm and use it to scan 11 data sets of vertebrate fibrillar collagens, and predict the periodicity of the resulting assemblies. We analyzed the experimentally observed variance of the optimal stagger distances across species, and find that these distances, and the resulting fibrillar phenotypes, are evolutionary well preserved. Moreover, we observed that the energy minimum at the optimal stagger distance is broad in all cases, suggesting a further evolutionary adaptation designed to improve the assembly kinetics. Our periodicity predictions are not only in good agreement with the experimental data on collagen molecular staggering for all collagen types analyzed, but also for synthetic peptides. We argue that, with our model, it becomes possible to design tailor-made, periodic collagen structures, thereby enabling the design of novel biomimetic materials based on collagen-mimetic trimers.

Numerical Test of the Onsager Relations in a Driven System.

The Onsager reciprocity relations were formulated in the context of irreversible thermodynamics, but they are based on assumptions that have a wider applicability. Here, we present simulations testing the Onsager relations between surface-coupled diffusive and bulk fluxes in a system prepared in a nonequilibrium steady state. The system consists of a mixture of two identical species maintained at different temperatures inside a channel. In order to tune the friction of the two species with the walls independently, while keeping the particle-wall interaction potentials the same, we allow the kinematics of particle-wall collisions to be different: "bounce-back" (B) or "specular" (S). In the BB case, diffusio-capillary transport can only take place if the two species have different temperatures. We find that the Onsager reciprocity relations are obeyed in the linear regime, even in the BB case where all fluxes are the result of perturbing the system from a nonequilibrium steady state in a way that does not satisfy time-reversal symmetry. Our Letter provides a direct, numerical illustration of the validity of the Onsager relations outside their original range of application, and suggests their relevance for transport in driven or active systems.

Thermodynamics and kinetics of phase separation of protein-RNA mixtures by a minimal model.

Intracellular liquid-liquid phase separation enables the formation of biomolecular condensates, such as ribonucleoprotein granules, which play a crucial role in the spatiotemporal organization of biomolecules (e.g., proteins and RNAs). Here, we introduce a patchy-particle polymer model to investigate liquid-liquid phase separation of protein-RNA mixtures. We demonstrate that at low to moderate concentrations, RNA enhances the stability of RNA-binding protein condensates because it increases the molecular connectivity of the condensed-liquid phase. Importantly, we find that RNA can also accelerate the nucleation stage of phase separation. Additionally, we assess how the capacity of RNA to increase the stability of condensates is modulated by the relative protein-protein/protein-RNA binding strengths. We find that phase separation and multiphase organization of multicomponent condensates is favored when the RNA binds with higher affinity to the lower-valency proteins in the mixture than to the cognate higher-valency proteins. Collectively, our results shed light on the roles of RNA in ribonucleoprotein granule formation and the internal structuring of stress granules.

Nano-pump based on exothermic surface reactions.

We present simulations indicating that it should be possible to construct a switchable nano-scale fluid pump, driven by exothermic surface reactions. Such a pump could, for instance, be controlled electro-chemically. In our simulations we explore a simple illustration of such a pump. We argue that the simplicity of the pump design could make it attractive for micro/nano-fluidics applications.

Computation of the chemical potential and solubility of amorphous solids.

Using a recently developed technique to estimate the equilibrium free energy of glassy materials, we explore if equilibrium simulation methods can be used to estimate the solubility of amorphous solids. As an illustration, we compute the chemical potentials of the constituent particles of a two-component Kob-Andersen model glass former. To compute the chemical potential for different components, we combine the calculation of the overall free energy of the glass with a calculation of the chemical potential difference of the two components. We find that the standard method to compute chemical potential differences by thermodynamic integration yields not only a wide scatter in the chemical potential values, but also, more seriously, the average of the thermodynamic integration results is well above the extrapolated value for the supercooled liquid. However, we find that if we compute the difference in the chemical potential of the components with the non-equilibrium free-energy expression proposed by Jarzynski, we obtain a good match with the extrapolated value of the supercooled liquid. The extension of the Jarzynski method that we propose opens a potentially powerful route to compute the free-energy related equilibrium properties of glasses. We find that the solubility estimate of amorphous materials obtained from direct-coexistence simulations is only in fair agreement with the solubility prediction based on the chemical potential calculations of a hypothetical "well-equilibrated glass." In direct-coexistence simulations, we find that, in qualitative agreement with experiments, the amorphous solubility decreases with time and attains a low solubility value.

Estimation of the equilibrium free energy for glasses using the Jarzynski equality.

The free energy of glasses cannot be estimated using thermodynamic integration as glasses are intrinsically not in equilibrium. We present numerical simulations showing that, in contrast, plausible free-energy estimates of a Kob-Andersen glass can be obtained using the Jarzynski relation. Using the Jarzynski relation, we also compute the chemical potential difference of the two components of this system and find that, in the glassy regime, the Jarzynski estimate matches well with the extrapolated value of the supercooled liquid. Our findings are of broader interest as they show that the Jarzynski method can be used under conditions where the thermodynamic integration approach, which is normally more accurate, breaks down completely. Systems where such an approach might be useful are gels and jammed glassy structures formed by compression.

Thermodynamics and kinetics of crystallization in deeply supercooled Stillinger-Weber silicon.

We study the kinetics of crystallization in deeply supercooled liquid silicon employing computer simulations and the Stillinger-Weber three-body potential. The free energy barriers to crystallization are computed using umbrella sampling Monte Carlo simulations and from unconstrained molecular dynamics simulations using a mean first passage time formulation. We focus on state points that have been described in earlier work [S. Sastry and C. A. Angell, Nat. Mater. 2, 739 (2003)] as straddling a liquid-liquid phase transition (LLPT) between two metastable liquid states. It was argued subsequently [Ricci et al., Mol. Phys. 117, 3254 (2019)] that the apparent transition is due to the loss of metastability of the liquid state with respect to the crystalline state. The presence of a barrier to crystallization for these state points is therefore of importance to ascertain, which we investigate, with due attention to ambiguities that may arise from the choice of order parameters. We find a well-defined free energy barrier to crystallization and demonstrate that both umbrella sampling and mean first passage time methods yield results that agree quantitatively. Our results thus provide strong evidence against the possibility that the liquids at state points close to the reported LLPT exhibit slow, spontaneous crystallization, but they do not address the existence of a LLPT (or lack thereof). We also compute the free energy barriers to crystallization at other state points over a broad range of temperatures and pressures and discuss the effect of changes in the microscopic structure of the metastable liquid on the free energy barrier heights.

Free energies of crystals computed using Einstein crystal with fixed center of mass and differing spring constants.

Free energies of crystals computed using a center of mass constraint require a finite-size correction, as shown in previous work by Polson et al. [J. Chem. Phys. 112, 5339-5342 (2000)]. Their reference system is an Einstein crystal with equal spring constants. In this paper, we extend the work of Polson et al. [J. Chem. Phys. 112, 5339-5342 (2000)] to the case of differing spring constants. The generalization is convenient for constraining the center of mass in crystals with atoms of differing masses, and it helps to optimize the free energy calculations. To test the theory, we compare the free energies of LiI and NaCl crystals from calculations with differing spring constants to those computed using equal spring constants. Using these center of mass finite size corrections, we compute the true free energies of these crystals for different system sizes to eliminate the intrinsic finite-size effects. These calculations help demonstrate the size of these finite-size corrections relative to other contributions to the absolute free energy of the crystals.

Computing the Heat Conductivity of Fluids from Density Fluctuations.

Equilibrium molecular dynamics simulations, in combination with the Green-Kubo (GK) method, have been extensively used to compute the thermal conductivity of liquids. However, the GK method relies on an ambiguous definition of the microscopic heat flux, which depends on how one chooses to distribute energies over atoms. This ambiguity makes it problematic to employ the GK method for systems with nonpairwise interactions. In this work, we show that the hydrodynamic description of thermally driven density fluctuations can be used to obtain the thermal conductivity of a bulk fluid unambiguously, thereby bypassing the need to define the heat flux. We verify that, for a model fluid with only pairwise interactions, our method yields estimates of thermal conductivity consistent with the GK approach. We apply our approach to compute the thermal conductivity of a nonpairwise additive water model at supercritical conditions, and of a liquid hydrogen system described by a machine-learning interatomic potential, at 33 GPa and 2000 K.

Reactive Momentum Transfer Contributes to the Self-Propulsion of Janus Particles.

We report simulations of a spherical Janus particle undergoing exothermic surface reactions around one pole only. Our model excludes self-phoretic transport by design. Nevertheless, net motion occurs from direct momentum transfer between solvent and colloid, with speed scaling as the square root of the energy released during the reaction. We find that such propulsion is dominated by the system's short-time response, when neither the time dependence of the flow around the colloid nor the solvent compressibility can be ignored. Our simulations agree reasonably well with previous experiments.

Effect of the interaction strength and anisotropy on the diffusio-phoresis of spherical colloids.

Gradients in temperature, concentration or electrostatic potential cannot exert forces on a bulk fluid; they can, however, exert forces on a fluid in a microscopic boundary layer surrounding a (nano)colloidal solute, resulting in so-called phoretic flow. Here we present a simulation study of phoretic flow around a spherical colloid held fixed in a concentration gradient. We show that the resulting flow velocity depends non-monotonically on the strength of the colloid-fluid interaction. The reason for this non-monotonic dependence is that solute particles are effectively trapped in a shell around the colloid and cannot contribute to diffusio-phoresis. We also observe that the flow depends sensitively on the anisotropy of solute-colloid interaction.

The Lennard-Jones potential: when (not) to use it.

The Lennard-Jones 12-6 potential (LJ) is arguably the most widely used pair potential in molecular simulations. In fact, it is so popular that the question is rarely asked whether it is fit for purpose. In this paper, we argue that, whilst the LJ potential was designed for noble gases such as argon, it is often used for systems where it is not expected to be particularly realistic. Under those circumstances, the disadvantages of the LJ potential become relevant: most important among these is that in simulations the LJ potential is always modified such that it has a finite range. More seriously, there is by now a whole family of different potentials that are all called Lennard-Jones 12-6, and that are all different - and that may have very different macroscopic properties. In this paper, we consider alternatives to the LJ 12-6 potential that could be employed under conditions where the LJ potential is only used as a typical short-ranged potential with attraction. We construct a class of potentials that are, in many respects LJ-like but that are by construction finite ranged, vanishing quadratically at the cut-off distance, and that are designed to be computationally cheap. Below, we present this potential and report numerical data for its thermodynamic and transport properties, for the most important cases (cut-off distance rc = 2σ ("LJ-like") and rc = 1.2σ (a typical "colloidal" potential)).

Optimal multivalent targeting of membranes with many distinct receptors.

Cells can often be recognized by the concentrations of receptors expressed on their surface. For better (targeted drug treatment) or worse (targeted infection by pathogens), it is clearly important to be able to target cells selectively. A good targeting strategy would result in strong binding to cells with the desired receptor profile and barely binding to other cells. Using a simple model, we formulate optimal design rules for multivalent particles that allow them to distinguish target cells based on their receptor profile. We find the following: (i) It is not a good idea to aim for very strong binding between the individual ligands on the guest (delivery vehicle) and the receptors on the host (cell). Rather, one should exploit multivalency: High sensitivity to the receptor density on the host can be achieved by coating the guest with many ligands that bind only weakly to the receptors on the cell surface. (ii) The concentration profile of the ligands on the guest should closely match the composition of the cognate membrane receptors on the target surface. And (iii) irrespective of all details, the effective strength of the ligand-receptor interaction should be of the order of the thermal energy [Formula: see text], where [Formula: see text] is the absolute temperature and [Formula: see text] is Boltzmann's constant. We present simulations that support the theoretical predictions. We speculate that, using the above design rules, it should be possible to achieve targeted drug delivery with a greatly reduced incidence of side effects.

Monte Carlo sampling for stochastic weight functions.

Conventional Monte Carlo simulations are stochastic in the sense that the acceptance of a trial move is decided by comparing a computed acceptance probability with a random number, uniformly distributed between 0 and 1. Here, we consider the case that the weight determining the acceptance probability itself is fluctuating. This situation is common in many numerical studies. We show that it is possible to construct a rigorous Monte Carlo algorithm that visits points in state space with a probability proportional to their average weight. The same approach may have applications for certain classes of high-throughput experiments and the analysis of noisy datasets.

Numerical evidence for thermally induced monopoles.

Electric charges are conserved. The same would be expected to hold for magnetic charges, yet magnetic monopoles have never been observed. It is therefore surprising that the laws of nonequilibrium thermodynamics, combined with Maxwell's equations, suggest that colloidal particles heated or cooled in certain polar or paramagnetic solvents may behave as if they carry an electric/magnetic charge. Here, we present numerical simulations that show that the field distribution around a pair of such heated/cooled colloidal particles agrees quantitatively with the theoretical predictions for a pair of oppositely charged electric or magnetic monopoles. However, in other respects, the nonequilibrium colloidal particles do not behave as monopoles: They cannot be moved by a homogeneous applied field. The numerical evidence for the monopole-like fields around heated/cooled colloidal particles is crucial because the experimental and numerical determination of forces between such colloidal particles would be complicated by the presence of other effects, such as thermophoresis.