Reconstructing diploid 3D chromatin structures from single cell Hi-C data with a polymer-based approach.

Detailed understanding of the 3D structure of chromatin is a key ingredient to investigate a variety of processes inside the cell. Since direct methods to experimentally ascertain these structures lack the desired spatial fidelity, computational inference methods based on single cell Hi-C data have gained significant interest. Here, we develop a progressive simulation protocol to iteratively improve the resolution of predicted interphase structures by maximum-likelihood association of ambiguous Hi-C contacts using lower-resolution predictions. Compared to state-of-the-art methods, our procedure is not limited to haploid cell data and allows us to reach a resolution of up to 5,000 base pairs per bead. High resolution chromatin models grant access to a multitude of structural phenomena. Exemplarily, we verify the formation of chromosome territories and holes near aggregated chromocenters as well as the inversion of the CpG content for rod photoreceptor cells.

Growth and aging in a few phase-separating active matter systems.

Via computer simulations we study evolution dynamics in systems of continuously moving active Brownian particles. The obtained results are discussed against those from the passive 2D Ising case. Following sudden quenches of random configurations to state points lying within the miscibility gaps and to the critical points, we investigate the far-from-steady-state dynamics by calculating quantities associated with structure and characteristic length scales. We also study aging for quenches into the miscibility gap and provide a quantitative picture for the scaling behavior of the two-time order-parameter correlation function. The overall structure and dynamics are consistent with expectations from the Ising model. This remains true for certain active lattice models as well, for which we present results for quenches to the critical points.

300-Times-Increased Diffusive Skyrmion Dynamics and Effective Pinning Reduction by Periodic Field Excitation.

Thermally induced skyrmion dynamics, as well as skyrmion pinning effects, in thin films have attracted significant interest. While pinning poses challenges in deterministic skyrmion devices and slows down skyrmion diffusion, for applications in non-conventional computing, both pinning of an appropriate strength and skyrmion diffusion speed are key. Here, periodic field excitations are employed to realize an increase of the skyrmion diffusion by more than two orders of magnitude. Amplifying the excitation, a drastic reduction of the effective skyrmion pinning, is reported, and a transition from pinning-dominated diffusive hopping to dynamics approaching free diffusion is observed. By tailoring the field oscillation frequency and amplitude, a continuous tuning of the effective pinning and skyrmion dynamics is demonstrated, which is a key asset and enabler for non-conventional computing applications. It is found that the periodic excitations additionally allow stabilization of skyrmions at different sizes for field values that are inaccessible in static systems, opening up new approaches to ultrafast skyrmion motion by transiently exciting moving skyrmions.

Can Polymer Helicity Affect Topological Chirality of Polymer Knots?

We investigate the effect of helicity in isolated polymers on the topological chirality of their knots with computer simulations. Polymers are described by generic worm-like chains (WLC), where helical conformations are promoted by chiral coupling between segments that are neighbors along the chain contour. The sign and magnitude of the coupling coefficient u determine the sense and strength of helicity. The corrugation of the helix is adjusted via the radius R of a spherical, hard excluded volume around each WLC segment. Open and compact helices are, respectively, obtained for R that is either zero or smaller than the length of the WLC bond, and R that is a few times larger than the bond length. We use a Monte Carlo algorithm to sample polymer conformations for different values of u, spanning the range from achiral polymers to chains with well-developed helices. Monitoring the average helix torsion and fluctuations of chiral order as a function of u, for two very different chain lengths, demonstrates that the coil-helix transition in this model is not a phase transition but a crossover. Statistical analysis of conformations forming the simplest chiral knots, 31, 51, and 52, demonstrates that topological mirror symmetry is broken─knots formed by helices with a given sense prefer one handedness over the other. For the 31 and 51 knots, positive helical sense favors positive handedness. Intriguingly, an opposite trend is observed for 52 knots, where positive helical sense promotes negative handedness. We argue that this special coupling between helicity and topological chirality stems from a generic mechanism: conformations where some of the knot crossings are found in "braids" formed by two tightly interwoven sections of the polymer.

AlphaFold predicts the most complex protein knot and composite protein knots.

The computer artificial intelligence system AlphaFold has recently predicted previously unknown three-dimensional structures of thousands of proteins. Focusing on the subset with high-confidence scores, we algorithmically analyze these predictions for cases where the protein backbone exhibits rare topological complexity, that is, knotting. Amongst others, we discovered a 71 -knot, the most topologically complex knot ever found in a protein, as well several six-crossing composite knots comprised of two methyltransferase or carbonic anhydrase domains, each containing a simple trefoil knot. These deeply embedded composite knots occur evidently by gene duplication and interconnection of knotted dimers. Finally, we report two new five-crossing knots including the first 51 -knot. Our list of analyzed structures forms the basis for future experimental studies to confirm these novel-knotted topologies and to explore their complex folding mechanisms.

Skyrmion pinning energetics in thin film systems.

A key issue for skyrmion dynamics and devices are pinning effects present in real systems. While posing a challenge for the realization of conventional skyrmionics devices, exploiting pinning effects can enable non-conventional computing approaches if the details of the pinning in real samples are quantified and understood. We demonstrate that using thermal skyrmion dynamics, we can characterize the pinning of a sample and we ascertain the spatially resolved energy landscape. To understand the mechanism of the pinning, we probe the strong skyrmion size and shape dependence of the pinning. Magnetic microscopy imaging demonstrates that in contrast to findings in previous investigations, for large skyrmions the pinning originates at the skyrmion boundary and not at its core. The boundary pinning is strongly influenced by the very complex pinning energy landscape that goes beyond the conventional effective rigid quasi-particle description. This gives rise to complex skyrmion shape distortions and allows for dynamic switching of pinning sites and flexible tuning of the pinning.

Knot formation of dsDNA pushed inside a nanochannel.

Recent experiments demonstrated that knots in single molecule dsDNA can be formed by compression in a nanochannel. In this manuscript, we further elucidate the underlying molecular mechanisms by carrying out a compression experiment in silico, where an equilibrated coarse-grained double-stranded DNA confined in a square channel is pushed by a piston. The probability of forming knots is a non-monotonic function of the persistence length and can be enhanced significantly by increasing the piston speed. Under compression knots are abundant and delocalized due to a backfolding mechanism from which chain-spanning loops emerge, while knots are less frequent and only weakly localized in equilibrium. Our in silico study thus provides insights into the formation, origin and control of DNA knots in nanopores.

Influence of ionic conditions on knotting in a coarse-grained model for DNA.

We investigate knotting probabilities of long double-stranded DNA strands in a coarse-grained Kratky-Porod model for DNA with Monte Carlo simulations. Various ionic conditions are implemented by adjusting the effective diameter of monomers. We find that the occurrence of knots in DNA can be reinforced considerably by high salt conditions and confinement between plates. Likewise, knots can almost be dissolved completely in a low salt scenario. Comparisons with recent experiments confirm that the coarse-grained model is able to capture and quantitatively predict topological features of DNA and can be used for guiding future experiments on DNA knots.

Shear-Thinning in Oligomer Melts-Molecular Origins and Applications.

We investigate the molecular origin of shear-thinning in melts of flexible, semiflexible and rigid oligomers with coarse-grained simulations of a sheared melt. Entanglements, alignment, stretching and tumbling modes or suppression of the latter all contribute to understanding how macroscopic flow properties emerge from the molecular level. In particular, we identify the rise and decline of entanglements with increasing chain stiffness as the major cause for the non-monotonic behaviour of the viscosity in equilibrium and at low shear rates, even for rather small oligomeric systems. At higher shear rates, chains align and disentangle, contributing to shear-thinning. By performing simulations of single chains in shear flow, we identify which of these phenomena are of collective nature and arise through interchain interactions and which are already present in dilute systems. Building upon these microscopic simulations, we identify by means of the Irving-Kirkwood formula the corresponding macroscopic stress tensor for a non-Newtonian polymer fluid. Shear-thinning effects in oligomer melts are also demonstrated by macroscopic simulations of channel flows. The latter have been obtained by the discontinuous Galerkin method approximating macroscopic polymer flows. Our study confirms the influence of microscopic details in the molecular structure of short polymers such as chain flexibility on macroscopic polymer flows.

Critical behavior in active lattice models of motility-induced phase separation.

Lattice models allow for a computationally efficient investigation of motility-induced phase separation (MIPS) compared to off-lattice systems. Simulations are less demanding, and thus, bigger systems can be accessed with higher accuracy and better statistics. In equilibrium, lattice and off-lattice models with comparable interactions belong to the same universality class. Whether concepts of universality also hold for active particles is still a controversial and open question. Here, we examine two recently proposed active lattice systems that undergo MIPS and investigate numerically their critical behavior. In particular, we examine the claim that these systems and MIPS in general belong to the Ising universality class. We also take a more detailed look on the influence and role of rotational diffusion and active velocity in these systems.

Comparing equilibration schemes of high-molecular-weight polymer melts with topological indicators.

Recent theoretical studies have demonstrated that the behaviour of molecular knots is a sensitive indicator of polymer structure. Here, we use knots to verify the ability of two state-of-the-art algorithms-configuration assembly and hierarchical backmapping-to equilibrate high-molecular-weight (MW) polymer melts. Specifically, we consider melts with MWs equivalent to several tens of entanglement lengths and various chain flexibilities, generated with both strategies. We compare their unknotting probability, unknotting length, knot spectra, and knot length distributions. The excellent agreement between the two independent methods with respect to knotting properties provides an additional strong validation of their ability to equilibrate dense high-MW polymeric liquids. By demonstrating this consistency of knotting behaviour, our study opens the way for studying topological properties of polymer melts beyond time and length scales accessible to brute-force molecular dynamics simulations.

Can Soft Models Describe Polymer Knots?

Similar to macroscopic ropes and cables, long polymers create knots. We address the fundamental question whether and under which conditions it is possible to describe these intriguing objects with crude models that capture only mesoscale polymer properties. We focus on melts of long polymers which we describe by a model typical for mesoscopic simulations. A worm-like chain model defines the polymer architecture. To describe nonbonded interactions, we deliberately choose a generic "soft" repulsive potential that leads to strongly overlapping monomers and coarse local liquid structure. The soft model is parametrized to accurately reproduce mesoscopic structure and conformations of reference polymer melts described by a microscopic model. The microscopically resolved samples retain all generic features affecting polymer topology and provide, therefore, reliable reference data on knots. We compare characteristic knotting properties in mesoscopic and microscopically resolved melts for different cases of chain stiffness. We conclude that mesoscopic models can reliably describe knots in those melts, where the length scale characterizing polymer stiffness is substantially larger than the size of monomer-monomer excluded volume. In this case, simplified local liquid structure influences knotting properties only marginally. In contrast, mesoscopic models perform poorly in melts with flexible chains. We qualitatively explain our findings through a free energy model of simple knots available in the literature.

Thermal skyrmion diffusion used in a reshuffler device.

Magnetic skyrmions in thin films can be efficiently displaced with high speed by using spin-transfer torques1,2 and spin-orbit torques3-5 at low current densities. Although this favourable combination of properties has raised expectations for using skyrmions in devices6,7, only a few publications have studied the thermal effects on the skyrmion dynamics8-10. However, thermally induced skyrmion dynamics can be used for applications11 such as unconventional computing approaches12, as they have been predicted to be useful for probabilistic computing devices13. In our work, we uncover thermal diffusive skyrmion dynamics by a combined experimental and numerical study. We probed the dynamics of magnetic skyrmions in a specially tailored low-pinning multilayer material. The observed thermally excited skyrmion motion dominates the dynamics. Analysing the diffusion as a function of temperature, we found an exponential dependence, which we confirmed by means of numerical simulations. The diffusion of skyrmions was further used in a signal reshuffling device as part of a skyrmion-based probabilistic computing architecture. Owing to its inherent two-dimensional texture, the observation of a diffusive motion of skyrmions in thin-film systems may also yield insights in soft-matter-like characteristics (for example, studies of fluctuation theorems, thermally induced roughening and so on), which thus makes it highly desirable to realize and study thermal effects in experimentally accessible skyrmion systems.

Proteins' Knotty Problems.

Knots in proteins are increasingly being recognized as an important structural concept, and the folding of these peculiar structures still poses considerable challenges. From a functional point of view, most protein knots discovered so far are either enzymes or DNA-binding proteins. Our comprehensive topological analysis of the Protein Data Bank reveals several novel structures including knotted mitochondrial proteins and the most deeply embedded protein knot discovered so far. For the latter, we propose a novel folding pathway based on the idea that a loose knot forms at a terminus and slides to its native position. For the mitochondrial proteins, we discuss the folding problem from the perspective of transport and suggest that they fold inside the mitochondria. We also discuss the evolutionary origin of a novel class of knotted membrane proteins and argue that a novel knotted DNA-binding protein constitutes a new fold. Finally, we have also discovered a knot in an artificially designed protein structure.

KnotGenome: a server to analyze entanglements of chromosomes.

The KnotGenome server enables the topological analysis of chromosome model data using three-dimensional coordinate files of chromosomes as input. In particular, it detects prime and composite knots in single chromosomes, and links between chromosomes. The knotting complexity of the chromosome is presented in the form of a matrix diagram that reveals the knot type of the entire polynucleotide chain and of each of its subchains. Links are determined by means of the Gaussian linking integral and the HOMFLY-PT polynomial. Entangled chromosomes are presented graphically in an intuitive way. It is also possible to relax structure with short molecular dynamics runs before the analysis. KnotGenome is freely available at http://knotgenom.cent.uw.edu.pl/.

Do the contact angle and line tension of surface-attached droplets depend on the radius of curvature?

Results from Monte Carlo simulations of wall-attached droplets in the three-dimensional Ising lattice gas model and in a symmetric binary Lennard-Jones fluid, confined by antisymmetric walls, are analyzed, with the aim to estimate the dependence of the contact angle [Formula: see text] on the droplet radius [Formula: see text] of curvature. Sphere-cap shape of the wall-attached droplets is assumed throughout. An approach, based purely on 'thermodynamic' observables, e.g. chemical potential, excess density due to the droplet, etc, is used, to avoid ambiguities in the decision which particles belong (or do not belong, respectively) to the droplet. It is found that the results are compatible with a variation [Formula: see text], [Formula: see text] being the contact angle in the thermodynamic limit ([Formula: see text]). The possibility to use such results to estimate the excess free energy related to the contact line of the droplet, namely the line tension, at the wall, is discussed. Various problems that hamper this approach and were not fully recognized in previous attempts to extract the line tension are identified. It is also found that the dependence of wall tensions on the difference of chemical potential of the droplet from that at the bulk coexistence provides effectively a change of the contact angle of similar magnitude. The simulation approach yields precise estimates for the excess density due to wall-attached droplets and the corresponding free energy excess, relative to a system without a droplet at the same chemical potential. It is shown that this information suffices to estimate nucleation barriers, not affected by ambiguities on droplet shape, contact angle and line tension.

Equilibrium between a Droplet and Surrounding Vapor: A Discussion of Finite Size Effects.

In a theoretical description of homogeneous nucleation one frequently assumes an "equilibrium" coexistence of a liquid droplet with surrounding vapor of a density exceeding that of a saturated vapor at bulk vapor-liquid two-phase coexistence. Thereby one ignores the caveat that in the thermodynamic limit, for which the vapor would be called supersaturated, such states will at best be metastable with finite lifetime, and thus not be well-defined within equilibrium statistical mechanics. In contrast, in a system of finite volume stable equilibrium coexistence of droplet and supersaturated vapor at constant total density is perfectly possible, and numerical analysis of equilibrium free energies of finite systems allows to obtain physically relevant results. In particular, such an analysis can be used to derive the dependence of the droplet surface tension γ( R) on the droplet radius R by computer simulations. Unfortunately, however, the precision of the results produced by this approach turns out to be seriously affected by a hitherto unexplained spurious dependence of γ( R) on the total volume V of the simulation box. These finite size effects are studied here for the standard Ising/lattice gas model in d = 2 dimensions and an Ising model on the face-centered cubic lattice with 3-spin interaction, lacking symmetry between vapor and liquid phases. There also the analogous case of bubbles surrounded by undersaturated liquid is treated. It is argued that (at least a large part of) the finite size effects result from the translation entropy of the droplet or bubble in the system. This effect has been shown earlier to occur also for planar interfaces for simulations in the slab geometry. Consequences for the estimation of the Tolman length are briefly discussed. In particular, we find clear evidence that in d = 2 the leading correction of the curvature-dependent interface tension is a logarithmic term, compatible with theoretical expectations, and we show that then the standard Tolman-style analysis is inapplicable.

Mapping onto Ideal Chains Overestimates Self-Entanglements in Polymer Melts.

We analyze the occurrence of knots, their spectrum, and sizes in polymer melts. Surprisingly, the number of knots in melt conformations is much lower than expected from a mapping to a random walk with the same Kuhn segment length. The effective random walk severely overrates the occurrence of knots and their complexity, particularly when compared to melts of flexible chains, indicating that nontrivial effects due to remnants of self-avoidance play a role for the chain lengths considered in this numerical study. For melt chains with higher persistence length, the effect is less pronounced. In addition, we find that chains in a melt have a knot structure very similar to dilute single chains close to the collapse transition. We finally show that typical equilibration procedures are well-suited to relax the topology in melts.

Estimation of the critical behavior in an active colloidal system with Vicsek-like interactions.

We study numerically the critical behavior of a modified, active Asakura-Oosawa model for colloid-polymer mixtures. The colloids are modeled as self-propelled particles with Vicsek-like interactions. This system undergoes phase separation between a colloid-rich and a polymer-rich phase, whereby the phase diagram depends on the strength of the Vicsek-like interactions. Employing a subsystem-block-density distribution analysis, we determine the critical point and make an attempt to estimate the critical exponents. In contrast to the passive model, we find that the critical point is not located on the rectilinear diameter. A first estimate of the critical exponents ß and ν is consistent with the underlying 3d-Ising universality class observed for the passive model.

Phase behavior of active Brownian disks, spheres, and dimers.

In this paper we provide high precision estimates of the phase diagram of active Brownian particles. We extract coexisting densities from simulations of phase separated states in an elongated box (slab geometry) which minimizes finite-size effects and allows for precise determination of points on the binodal lines. Using this method, we study the influence of both shape and dimensionality on the two-phase region. Active spheres and dimers of active particles are compared to the known phase diagram of active Brownian disks. In the case of dimers, both correlated and uncorrelated propulsion of the two beads are studied. The influence of correlations is discussed through a simple mapping.