DNA-DNA interactions in bacteriophage capsids are responsible for the observed DNA knotting

The structure and physical properties of a packaged bacteriophage particle

A string of nucleotides confined within a protein capsid contains all the instructions necessary to make a functional virus particle, a virion. Although the structure of the protein capsid is known for many virus species1,2, the three-dimensional organization of viral genomes has mostly eluded experimental probes3,4. Here we report all-atom structural models of an HK97 virion5, including its entire 39,732 base pair genome, obtained through multiresolution simulations. Mimicking the action of a packaging motor6, the genome was gradually loaded into the capsid. The structure of the packaged capsid was then refined through simulations of increasing resolution, which produced a 26 million atom model of the complete virion, including water and ions confined within the capsid. DNA packaging occurs through a loop extrusion mechanism7 that produces globally different configurations of the packaged genome and gives each viral particle individual traits. Multiple microsecond-long all-atom simulations characterized the effect of the packaged genome on capsid structure, internal pressure, electrostatics and diffusion of water, ions and DNA, and revealed the structural imprints of the capsid onto the genome. Our approach can be generalized to obtain complete all-atom structural models of other virus species, thereby potentially revealing new drug targets at the genome-capsid interface.

Quantum-inspired encoding enhances stochastic sampling of soft matter systems

Quantum advantage in solving physical problems is still hard to assess due to hardware limitations. However, algorithms designed for quantum computers may engender transformative frameworks for modeling and simulating paradigmatically hard systems. Here, we show that the quadratic unconstrained binary optimization encoding enables tackling classical many-body systems that are challenging for conventional Monte Carlo. Specifically, in self-assembled melts of rigid lattice ring polymers, the combination of high density, chain stiffness, and topological constraints results in divergent autocorrelation times for real-space Monte Carlo. Our quantum-inspired encoding overcomes this problem and enables sampling melts of lattice rings with fixed curvature and compactness, unveiling counterintuitive topological effects. Tackling the same problems with the D-Wave quantum annealer leads to substantial performance improvements and advantageous scaling of sampling computational cost with the size of the self-assembled ring melts.

Why Na+ has higher propensity than K+ to condense DNA in a crowded environment

Experimentally, in the presence of the crowding agent polyethylene glycol (PEG), sodium ions compact double-stranded DNA more readily than potassium ions. Here, we have used molecular dynamics simulations and the "ion binding shells model" of DNA condensation to provide an explanation for the observed variations in condensation of short DNA duplexes in solutions containing different monovalent cations and PEG; several predictions are made. According to the model we use, externally bound ions contribute the most to the ion-induced aggregation of DNA duplexes. The simulations reveal that for two adjacent DNA duplexes, the number of externally bound Na+ ions is larger than the number of K+ ions over a wide range of chloride concentrations in the presence of PEG, providing a qualitative explanation for the higher propensity of sodium ions to compact DNA under crowded conditions. The qualitative picture is confirmed by an estimate of the corresponding free energy of DNA aggregation that is at least 0.2kBT per base pair more favorable in solution with NaCl than with KCl at the same ion concentration. The estimated attraction free energy of DNA duplexes in the presence of Na+ depends noticeably on the DNA sequence; we predict that AT-rich DNA duplexes are more readily condensed than GC-rich ones in the presence of Na+. Counter-intuitively, the addition of a small amount of a crowding agent with high affinity for the specific condensing ion may lead to the weakening of the ion-mediated DNA-DNA attraction, shifting the equilibrium away from the DNA condensed phase.

Role of DNA-DNA sliding friction and nonequilibrium dynamics in viral genome ejection and packaging

Many viruses eject their DNA via a nanochannel in the viral shell, driven by internal forces arising from the high-density genome packing. The speed of DNA exit is controlled by friction forces that limit the molecular mobility, but the nature of this friction is unknown. We introduce a method to probe the mobility of the tightly confined DNA by measuring DNA exit from phage phi29 capsids with optical tweezers. We measure extremely low initial exit velocity, a regime of exponentially increasing velocity, stochastic pausing that dominates the kinetics and large dynamic heterogeneity. Measurements with variable applied force provide evidence that the initial velocity is controlled by DNA-DNA sliding friction, consistent with a Frenkel-Kontorova model for nanoscale friction. We confirm several aspects of the ejection dynamics predicted by theoretical models. Features of the pausing suggest that it is connected to the phenomenon of 'clogging' in soft matter systems. Our results provide evidence that DNA-DNA friction and clogging control the DNA exit dynamics, but that this friction does not significantly affect DNA packaging.

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, 3₁, 5₁, and 5₂, demonstrates that topological mirror symmetry is broken─knots formed by helices with a given sense prefer one handedness over the other. For the 3₁ and 5₁ knots, positive helical sense favors positive handedness. Intriguingly, an opposite trend is observed for 5₂ 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.

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.

DNA compaction and dynamic observation in a nanopore gated sub-attoliter silicon nanocavity

Confinement of biopolymers inside volumes with micro- or nanoscale lateral dimensions is ubiquitous in nature. Investigating the behavior of biopolymers in a confined environment is essential to improve our basic understanding in life sciences. In this work, we present a nanopore gated sub-attoliter silicon nanocavity device, which allows DNA compaction similar to that in virus capsids. Single DNA molecules can be electrically driven into the nanocavity, and then get compacted inside the nanocavity under certain conditions. The dynamic fluctuations of the compacted DNA can be monitored via ionic current measurements. The mechanism for the DNA compaction is elucidated by varying the DNA length or concentration, voltage polarity, nanocavity dimensions and ionic strength. Furthermore, Brownian dynamics simulations reveal the dynamic fluctuations of the compacted DNA, which are reflected in the measured ionic current. Our nanocavity device is anticipated to provide a controlled environment in extremely small volumes for investigating the physics of confined biopolymers.

Statistics and topology of fluctuating ribbons

Ribbons are a class of slender structures whose length, width, and thickness are widely separated from each other. This scale separation gives a ribbon unusual mechanical properties in athermal macroscopic settings, for example, it can bend without twisting, but cannot twist without bending. Given the ubiquity of ribbon-like biopolymers in biology and chemistry, here we study the statistical mechanics of microscopic inextensible, fluctuating ribbons loaded by forces and torques. We show that these ribbons exhibit a range of topologically and geometrically complex morphologies exemplified by three phases-a twist-dominated helical phase (HT), a writhe-dominated helical phase (HW), and an entangled phase-that arise as the applied torque and force are varied. Furthermore, the transition from HW to HT phases is characterized by the spontaneous breaking of parity symmetry and the disappearance of perversions (that correspond to chirality-reversing localized defects). This leads to a universal response curve of a topological quantity, the link, as a function of the applied torque that is similar to magnetization curves in second-order phase transitions.

The Local Topological Free Energy of the SARS-CoV-2 Spike Protein

The novel coronavirus SARS-CoV-2 infects human cells using a mechanism that involves binding and structural rearrangement of its Spike protein. Understanding protein rearrangement and identifying specific amino acids where mutations affect protein rearrangement has attracted much attention for drug development. In this manuscript, we use a mathematical method to characterize the local topology/geometry of the SARS-CoV-2 Spike protein backbone. Our results show that local conformational changes in the FP, HR1, and CH domains are associated with global conformational changes in the RBD domain. The SARS-CoV-2 variants analyzed in this manuscript (alpha, beta, gamma, delta Mink, G614, N501) show differences in the local conformations of the FP, HR1, and CH domains as well. Finally, most mutations of concern are either in or in the vicinity of high local topological free energy conformations, suggesting that high local topological free energy conformations could be targets for mutations with significant impact of protein function. Namely, the residues 484, 570, 614, 796, and 969, which are present in variants of concern and are targeted as important in protein function, are predicted as such from our model.

Geometric Predictors of Knotted and Linked Arcs

Inspired by how certain proteins “sense” knots and entanglements in DNA molecules, here, we ask if local geometric features that may be used as a readout of the underlying topology of generic polymers exist. We perform molecular simulations of knotted and linked semiflexible polymers and study four geometric measures to predict topological entanglements: local curvature, local density, local 1D writhe, and nonlocal 3D writhe. We discover that local curvature is a poor predictor of entanglements. In contrast, segments with maximum local density or writhe correlate as much as 90% of the time with the shortest knotted and linked arcs. We find that this accuracy is preserved across different knot types and also under significant spherical confinement, which is known to delocalize essential crossings in knotted polymers. We further discover that nonlocal 3D writhe is the best geometric readout of the knot location. Finally, we discuss how these geometric features may be used to computationally analyze entanglements in generic polymer melts and gels.

Interplay between nematic and cholesteric interactions in self-consistent field theory

Chirality is a design feature of a number of biomolecules (e.g., collagen). In these molecules, cholesteric (chiral-nematic) behavior emerges from a combination of the tendency for the biopolymers to align (nematic interactions) and for the alignment direction to change with position, rotating around an axis normal to the alignment direction. This paper presents self-consistent field theory (SCFT) of chiral-nematic polymers, which takes into account polymer flexibility and the orientational degrees of freedom of polymer segments. Using the resulting SCFT, we construct a phase diagram showing regions of stability for isotropic, nematic, and cholesteric phases. Furthermore, we find that nematic interactions can stabilize the cholesteric phase, pushing the isotropic-cholesteric phase transition to lower cholesteric interaction strength, until the isotropic-nematic-cholesteric triple point is reached.

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.

Quantifying the effects of slit confinement on polymer knots using the tube model

Knots can spontaneously form in DNA, proteins, and other polymers and affect their properties. These knots often experience spatial confinement in biological systems and experiments. While confinement dramatically affects the knot behavior, the physical mechanisms underlying the confinement effects are not fully understood. In this work, we provide a simple physical picture of the polymer knots in slit confinement using the tube model. In the tube model, the polymer segments in the knot core are assumed to be confined in a virtual tube due to the topological restriction. We first perform Monte Carlo simulation of a flexible knotted chain confined in a slit. We find that with the decrease of the slit height from H=+∞ (the 3D case) to H=2a (the 2D case), the most probable knot size L_{knot}^{*} dramatically shrinks from (L_{knot}^{*})_{3D}≈140a to (L_{knot}^{*})_{2D}≈26a, where a is the monomer diameter of the flexible chain. Then we quantitatively explain the confinement-induced knot shrinking and knot deformation using the tube model. Our results for H=2a can be applied to a polymer knot on a surface, which resembles DNA knots measured by atomic force microscopy under the conditions that DNA molecules are weakly absorbed on the surface and reach equilibrium 2D conformations. This work demonstrates the effectiveness of the tube model in understanding polymer knots.

Topological and physical links in soft matter systems

Linking, or multicomponent topological entanglement, is ubiquitous in soft matter systems, from mixtures of polymers and DNA filaments packedin vivoto interlocked line defects in liquid crystals and intertwined synthetic molecules. Yet, it is only relatively recently that theoretical and experimental advancements have made it possible to probe such entanglements and elucidate their impact on the physical properties of the systems. Here, we review the state-of-the-art of this rapidly expanding subject and organize it as follows. First, we present the main concepts and notions, from topological linking to physical linking and then consider the salient manifestations of molecular linking, from synthetic to biological ones. We next cover the main physical models addressing mutual entanglements in mixtures of polymers, both linear and circular. Finally, we consider liquid crystals, fluids and other non-filamentous systems where topological or physical entanglements are observed in defect or flux lines. We conclude with a perspective on open challenges.

Vassiliev measures of complexity of open and closed curves in 3-space

In this article, we define Vassiliev measures of complexity for open curves in 3-space. These are related to the coefficients of the enhanced Jones polynomial of open curves in 3-space. These Vassiliev measures are continuous functions of the curve coordinates; as the ends of the curve tend to coincide, they converge to the corresponding Vassiliev invariants of the resulting knot. We focus on the second Vassiliev measure from the enhanced Jones polynomial for closed and open curves in 3-space. For closed curves, this second Vassiliev measure can be computed by a Gauss code diagram and it has an integral formulation, the double alternating self-linking integral. The double alternating self-linking integral is a topological invariant of closed curves and a continuous function of the curve coordinates for open curves in 3-space. For polygonal curves, the double alternating self-linking integral obtains a simpler expression in terms of geometric probabilities.

The local topological free energy of proteins

Protein folding, the process by which proteins attain a 3-dimensional conformation necessary for their function, remains an important unsolved problem in biology. A major gap in our understanding is how local properties of proteins relate to their global properties. In this manuscript, we use the Writhe and Torsion to introduce a new local topological/geometrical free energy that can be associated to 4 consecutive amino acids along the protein backbone. By analyzing a culled protein dataset from the PDB, our results show that high local topological free energy conformations are independent of sequence and may be involved in the rate limiting step in protein folding. By analyzing a set of 2-state single domain proteins, we find that the total local topological free energy of these proteins correlates with the experimentally observed folding rates reported in Plaxco et al. (2000).

Ion-dependent DNA configuration in bacteriophage capsids

Bacteriophages densely pack their long double-stranded DNA genome inside a protein capsid. The conformation of the viral genome inside the capsid is consistent with a hexagonal liquid crystalline structure. Experiments have confirmed that the details of the hexagonal packing depend on the electrochemistry of the capsid and its environment. In this work, we propose a biophysical model that quantifies the relationship between DNA configurations inside bacteriophage capsids and the types and concentrations of ions present in a biological system. We introduce an expression for the free energy that combines the electrostatic energy with contributions from bending of individual segments of DNA and Lennard-Jones-type interactions between these segments. The equilibrium points of this energy solve a partial differential equation that defines the distributions of DNA and the ions inside the capsid. We develop a computational approach that allows us to simulate much larger systems than what is possible using the existing molecular-level methods. In particular, we are able to estimate bending and repulsion between the DNA segments as well as the full electrochemistry of the solution, both inside and outside of the capsid. The numerical results show good agreement with existing experiments and with molecular dynamics simulations for small capsids.

Chromonic liquid crystals and packing configurations of bacteriophage viruses

We study equilibrium configurations of hexagonal columnar liquid crystals in the context of characterizing packing structures of bacteriophage viruses in a protein capsid. These are viruses that infect bacteria and are currently the focus of intense research efforts, with the goal of finding new therapies for bacteria-resistant antibiotics. The energy that we propose consists of the Oseen-Frank free energy of nematic liquid crystals that penalizes bending of the columnar directions, in addition to the cross-sectional elastic energy accounting for distortions of the transverse hexagonal structure; we also consider the isotropic contribution of the core and the energy of the unknown interface between the outer ordered region of the capsid and the inner disordered core. The problem becomes of free boundary type, with constraints. We show that the concentric, azimuthal, spool-like configuration is the absolute minimizer. Moreover, we present examples of toroidal structures formed by DNA in free solution and compare them with the analogous ones occurring in experiments with other types of lyotropic liquid crystals, such as food dyes and additives. This article is part of the theme issue 'Topics in mathematical design of complex materials'.

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.

Buckling transitions and soft-phase invasion of two-component icosahedral shells

What is the optimal distribution of two types of crystalline phases on the surface of icosahedral shells, such as of many viral capsids? We here investigate the distribution of a thin layer of soft material on a crystalline convex icosahedral shell. We demonstrate how the shapes of spherical viruses can be understood from the perspective of elasticity theory of thin two-component shells. We develop a theory of shape transformations of an icosahedral shell upon addition of a softer, but still crystalline, material onto its surface. We show how the soft component "invades" the regions with the highest elastic energy and stress imposed by the 12 topological defects on the surface. We explore the phase diagram as a function of the surface fraction of the soft material, the shell size, and the incommensurability of the elastic moduli of the rigid and soft phases. We find that, as expected, progressive filling of the rigid shell by the soft phase starts from the most deformed regions of the icosahedron. With a progressively increasing soft-phase coverage, the spherical segments of domes are filled first (12 vertices of the shell), then the cylindrical segments connecting the domes (30 edges) are invaded, and, ultimately, the 20 flat faces of the icosahedral shell tend to be occupied by the soft material. We present a detailed theoretical investigation of the first two stages of this invasion process and develop a model of morphological changes of the cone structure that permits noncircular cross sections. In conclusion, we discuss the biological relevance of some structures predicted from our calculations, in particular for the shape of viral capsids.

Topological Disentanglement of Linear Polymers under Tension

We develop a theoretical description of the topological disentanglement occurring when torus knots reach the ends of a semiflexible polymer under tension. These include decays into simpler knots and total unknotting. The minimal number of crossings and the minimal knot contour length are the topological invariants playing a key role in the model. The crossings behave as particles diffusing along the chain and the application of appropriate boundary conditions at the ends of the chain accounts for the knot disentanglement. Starting from the number of particles and their positions, suitable rules allow reconstructing the type and location of the knot moving on the chain Our theory is extensively benchmarked with corresponding molecular dynamics simulations and the results show a remarkable agreement between the simulations and the theoretical predictions of the model.

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.

Quantitative Study of the Chiral Organization of the Phage Genome Induced by the Packaging Motor

Molecular motors that translocate DNA are ubiquitous in nature. During morphogenesis of double-stranded DNA bacteriophages, a molecular motor drives the viral genome inside a protein capsid. Several models have been proposed for the three-dimensional geometry of the packaged genome, but very little is known of the signature of the molecular packaging motor. For instance, biophysical experiments show that in some systems, DNA rotates during the packaging reaction, but most current biophysical models fail to incorporate this property. Furthermore, studies including rotation mechanisms have reached contradictory conclusions. In this study, we compare the geometrical signatures imposed by different possible mechanisms for the packaging motors: rotation, revolution, and rotation with revolution. We used a previously proposed kinetic Monte Carlo model of the motor, combined with Brownian dynamics simulations of DNA to simulate deterministic and stochastic motor models. We find that rotation is necessary for the accumulation of DNA writhe and for the chiral organization of the genome. We observe that although in the initial steps of the packaging reaction, the torsional strain of the genome is released by rotation of the molecule, in the later stages, it is released by the accumulation of writhe. We suggest that the molecular motor plays a key role in determining the final structure of the encapsidated genome in bacteriophages.

Identifying knot types of polymer conformations by machine learning

We investigate the use of artificial neural networks (NNs) as an alternative tool to current analytical methods for recognizing knots in a given polymer conformation. The motivation is twofold. First, it is of interest to examine whether NNs are effective at learning the global and sequential properties that uniquely define a knot. Second, knot classification is an important and unsolved problem in mathematical and physical sciences, and NNs may provide insights into this problem. Motivated by these points, we generate millions of polymer conformations for five knot types: 0, 3_{1}, 4_{1}, 5_{1}, and 5_{2}, and we design various NN models for classification. Our best model achieves a five-class classification accuracy of above 99% on a polymer of 100 monomers. We find that the sequential modeling ability of recurrent NNs is crucial for this result, as it outperforms feed-forward NNs and successfully generalizes to differently sized conformations as well. We present our methods and suggest that deep learning may be used in specific applications of knot detection where some error is permissible. Hopefully, with further development, NNs can offer an alternative computational method for knot identification and facilitate knot research in mathematical and physical sciences.

Fine structure of viral dsDNA encapsidation

Unraveling the mechanisms of packing of DNA inside viral capsids is of fundamental importance to understanding the spread of viruses. It could also help develop new applications to targeted drug delivery devices for a large range of therapies. In this article, we present a robust, predictive mathematical model and its numerical implementation to aid the study and design of bacteriophage viruses for application purposes. Exploiting the analogies between the columnar hexagonal chromonic phases of encapsidated viral DNA and chromonic aggregates formed by plank-shaped molecular compounds, we develop a first-principles effective mechanical model of DNA packing in a viral capsid. The proposed expression of the packing energy, which combines relevant aspects of the liquid crystal theory, is developed from the model of hexagonal columnar phases, together with that describing configurations of polymeric liquid crystals. The method also outlines a parameter selection strategy that uses available data for a collection of viruses, aimed at applications to viral design. The outcome of the work is a mathematical model and its numerical algorithm, based on the method of finite elements, and computer simulations to identify and label the ordered and disordered regions of the capsid and calculate the inner pressure. It also presents the tools for the local reconstruction of the DNA "scaffolding" and the center curve of the filament within the capsid.

Field-driven tracer diffusion through curved bottlenecks: fine structure of first passage events

Using scaling arguments and extensive numerical simulations, we study the dynamics of a tracer particle in a corrugated channel represented by a periodic sequence of broad chambers and narrow funnel-like bottlenecks enclosed by a hard-wall boundary.

Computational methods in the study of self-entangled proteins: a critical appraisal

The existence of self-entangled proteins, the native structure of which features a complex topology, unveils puzzling, and thus fascinating, aspects of protein biology and evolution. The discovery that a polypeptide chain can encode the capability to self-entangle in an efficient and reproducible way during folding, has raised many questions, regarding the possible function of these knots, their conservation along evolution, and their role in the folding paradigm. Understanding the function and origin of these entanglements would lead to deep implications in protein science, and this has stimulated the scientific community to investigate self-entangled proteins for decades by now. In this endeavour, advanced experimental techniques are more and more supported by computational approaches, that can provide theoretical guidelines for the interpretation of experimental results, and for the effective design of new experiments. In this review we provide an introduction to the computational study of self-entangled proteins, focusing in particular on the methodological developments related to this research field. A comprehensive collection of techniques is gathered, ranging from knot theory algorithms, that allow detection and classification of protein topology, to Monte Carlo or molecular dynamics strategies, that constitute crucial instruments for investigating thermodynamics and kinetics of this class of proteins.

Spontaneous Domain Formation in Spherically Confined Elastic Filaments

Although the free energy of a genome packing into a virus is dominated by DNA-DNA interactions, ordering of the DNA inside the capsid is elasticity driven, suggesting general solutions with DNA organized into spool-like domains. Using analytical calculations and computer simulations of a long elastic filament confined to a spherical container, we show that the ground state is not a single spool as assumed hitherto, but an ordering mosaic of multiple homogeneously ordered domains. At low densities, we observe concentric spools, while at higher densities, other morphologies emerge, which resemble topological links. We discuss our results in the context of metallic wires, viral DNA, and flexible polymers.

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.

Effects of spherical confinement and backbone stiffness on flexible polymer jamming

We use molecular simulations to study jamming of a crumpled bead-spring model polymer in a finite container and compare to jamming of repulsive spheres. After proper constraint counting, the onset of rigidity is seen to occur isostatically as in the case of repulsive spheres. Despite this commonality, the presence of the curved container wall and polymer backbone bonds introduce new mechanical properties. Notably, these include additional bands in the vibrational density of states that reflect the material structure as well as oscillations in local contact number and density near the wall but with lower amplitude for polymers. Polymers have fewer boundary contacts, and this low-density surface layer strongly reduces the global bulk modulus. We further show that bulk-modulus dependence on backbone stiffness can be described by a model of stiffnesses in series and discuss potential experimental and biological applications.

Topological Disentanglement Dynamics of Torus Knots on Open Linear Polymers

We simulate and study the topological disentanglement occurring when torus knots reach the ends of a semiflexible open polymer (decay into simpler knots or unknotting). Through a rescaling procedure and the application of appropriate boundary conditions, we show that the full unknotting process can be understood in terms of point-like particles representing essential crossings, diffusing on the support [0, 1]. We address the bending and configurational free energy drives on the diffusion process, together with the scaling properties of the effective diffusion and friction coefficients. Agreement with simulations suggests universal features for these two model parameters.

Synergy of topoisomerase and structural-maintenance-of-chromosomes proteins creates a universal pathway to simplify genome topology

Topological entanglements severely interfere with important biological processes. For this reason, genomes must be kept unknotted and unlinked during most of a cell cycle. Type II topoisomerase (TopoII) enzymes play an important role in this process but the precise mechanisms yielding systematic disentanglement of DNA in vivo are not clear. Here we report computational evidence that structural-maintenance-of-chromosomes (SMC) proteins-such as cohesins and condensins-can cooperate with TopoII to establish a synergistic mechanism to resolve topological entanglements. SMC-driven loop extrusion (or diffusion) induces the spatial localization of essential crossings, in turn catalyzing the simplification of knots and links by TopoII enzymes even in crowded and confined conditions. The mechanism we uncover is universal in that it does not qualitatively depend on the specific substrate, whether DNA or chromatin, or on SMC processivity; we thus argue that this synergy may be at work across organisms and throughout the cell cycle.

Topological entanglement of interlocked knotted-unknotted polymer rings

Topological entanglements in biopolymers could drive them to certain internal statics and dynamics with important implications for biological functions. In this study, by means of molecular dynamics simulations, we demonstrate that the minimal crossing pattern of a braid plays a major role in its structural and dynamical properties; the braid consists of a knotted ring and an interlocked entwined unknotted polymer ring. In particular, we show that depending on the bending rigidity of the chains, the conformational energy of the braid can be either lower or higher than the unlocked polymer rings. Additionally, we find that a non-identical crossing pattern in the braid could distinctly enforce concerted internal conformational fluctuations between the interlocked rings.

Polymer Translocation Across a Corrugated Channel: Fick⁻Jacobs Approximation Extended Beyond the Mean First-Passage Time

Polymer translocation across a corrugated channel is a paradigmatic stochastic process encountered in diverse systems. The instance of time when a polymer first arrives to some prescribed location defines an important characteristic time-scale for various phenomena, which are triggered or controlled by such an event. Here we discuss the translocation dynamics of a Gaussian polymer in a periodically-corrugated channel using an appropriately generalized Fick⁻Jacobs approach. Our main aim is to probe an effective broadness of the first-passage time distribution (FPTD), by determining the so-called coefficient of variation γ of the FPTD, defined as the ratio of the standard deviation versus the mean first-passage time (MFPT). We present a systematic analysis of γ as a function of a variety of system's parameters. We show that γ never significantly drops below 1 and, in fact, can attain very large values, implying that the MFPT alone cannot characterize the first-passage statistics of the translocation process exhaustively well.

Topological Sieving of Rings According to Their Rigidity

We present a novel mechanism for resolving the mechanical rigidity of nanoscopic circular polymers that flow in a complex environment. The emergence of a regime of negative differential mobility induced by topological interactions between the rings and the substrate is the key mechanism for selective sieving of circular polymers with distinct flexibilities. A simple model accurately describes the sieving process observed in molecular dynamics simulations and yields experimentally verifiable analytical predictions, which can be used as a reference guide for improving filtration procedures of circular filaments. The topological sieving mechanism we propose ought to be relevant also in probing the microscopic details of complex substrates.

Dynamics of supercoiled DNA with complex knots: large-scale rearrangements and persistent multi-strand interlocking

Knots and supercoiling are both introduced in bacterial plasmids by catalytic processes involving DNA strand passages. While the effects on plasmid organization has been extensively studied for knotting and supercoiling taken separately, much less is known about their concurrent action. Here, we use molecular dynamics simulations and oxDNA, an accurate mesoscopic DNA model, to study the kinetic and metric changes introduced by complex (five-crossing) knots and supercoiling in 2 kbp-long DNA rings. We find several unexpected results. First, the conformational ensemble is dominated by two distinct states, differing in branchedness and knot size. Secondly, fluctuations between these states are as fast as the metric relaxation of unknotted rings. In spite of this, certain boundaries of knotted and plectonemically-wound regions can persist over much longer timescales. These pinned regions involve multiple strands that are interlocked by the cooperative action of topological and supercoiling constraints. Their long-lived character may be relevant for the simplifying action of topoisomerases.

Characterising knotting properties of polymers in nanochannels

Using a lattice model of polymers in a tube, we define one way to characterise different configurations of a given knot as either "local" or "non-local", based on a standard approach for measuring the "size" of a knot within a knotted polymer chain. The method involves associating knot-types to subarcs of the chain, and then identifying a knotted subarc with minimal arclength; this arclength is then the knot-size. If the resulting knot-size is small relative to the whole length of the chain, then the knot is considered to be localised or "local"; otherwise, it is "non-local". Using this definition, we establish that all but exponentially few sufficiently long self-avoiding polygons (closed chains) in a tubular sublattice of the simple cubic lattice are "non-locally" knotted. This is shown to also hold for the case when the same polygons are subject to an external tensile force, as well as in the extreme case when they are as compact as possible (no empty lattice sites). We also provide numerical evidence for small tube sizes that at equilibrium non-local knotting is more likely than local knotting, regardless of the strength of the stretching or compressing force. The relevance of these results to other models and recent experiments involving DNA knots is also discussed.

KymoKnot: A web server and software package to identify and locate knots in trajectories of linear or circular polymers

The KymoKnot software package and web server identifies and locates physical knots or proper knots in a series of polymer conformations. It is mainly intended as an analysis tool for trajectories of linear or circular polymers, but it can be used on single instances too, e.g. protein structures in PDB format. A key element of the software package is the so-called minimally interfering chain closure algorithm that is used to detect physical knots in open chains and to locate the knotted region in both open and closed chains. The web server offers a user-friendly graphical interface that identifies the knot type and highlights the knotted region on each frame of the trajectory, which the user can visualize interactively from various viewpoints. The dynamical evolution of the knotted region along the chain contour is presented as a kymograph. All data can be downloaded in text format. The KymoKnot package is licensed under the BSD 3-Clause licence. The server is publicly available at http://kymoknot.sissa.it/kymoknot/interactive.php .