Knotting probability of DNA molecules confined in restricted volumes: DNA knotting in phage capsids

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.

Following the entangled state of filaments

California blackworms serve as a template for the topological design of active matter.

Energetic preference and topological constraint effects on the formation of DNA twisted toroidal bundles

DNA toroids are compact torus-shaped bundles formed by one or multiple DNA molecules being condensed from the solution due to various condensing agents. It has been shown that the DNA toroidal bundles are twisted. However, the global conformations of DNA inside these bundles are still not well understood. In this study, we investigate this issue by solving different models for the toroidal bundles and performing replica-exchange molecular dynamics (REMD) simulations for self-attractive stiff polymers of various chain lengths. We find that a moderate degree of twisting is energetically favorable for toroidal bundles, yielding optimal configurations of lower energies than for other bundles corresponding to spool-like and constant radius of curvature arrangements. The REMD simulations show that the ground states of the stiff polymers are twisted toroidal bundles with the average twist degrees close to those predicted by the theoretical model. Constant-temperature simulations show that twisted toroidal bundles can be formed through successive processes of nucleation, growth, quick tightening, and slow tightening of the toroid, with the two last processes facilitating the polymer threading through the toroid's hole. A relatively long chain of 512 beads has an increased dynamical difficulty to access the twisted bundle states due to the polymer's topological constraint. Interestingly, we also observed significantly twisted toroidal bundles with a sharp U-shaped region in the polymer conformation. It is suggested that this U-shaped region makes the formation of twisted bundles easier by effectively reducing the polymer length. This effect can be equivalent to having multiple chains in the toroid.

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.

Expansion of Single Chains Released from a Spherical Cavity

A two-stage model is developed to explain the phenomena of chain expansion, released from a confining cavity. In the first stage, the chain is assumed to expand as a sphere, while in the second stage it expands like a coil. The kinetic equations for the variation of chain size are derived in the two stages by balancing the rate of the free energy change with the rate of the energy dissipation. Langevin dynamics simulations are then performed to examine the theory. We find that the expansion process is dominated by the second stage and the evolution of chain size follows, mainly, the predicted curve for coil expansion, which depends on the chain length and is not sensitive to the confining volume fraction. It permits to define the expansion time for the process. Further study reveals that the chain does undergo a spherical expansion in the first stage with the characteristic time much shorter than the one for the second stage. As a consequence, the first-stage variation of chain size can be regarded as an add-on to the principal curve of expansion designated by the second stage. The scaling behaviors and the associated scaling exponents are analyzed in details. The simulation results well support the theory.

Biomolecular Topology: Modelling and Analysis

With the great advancement of experimental tools, a tremendous amount of biomolecular data has been generated and accumulated in various databases. The high dimensionality, structural complexity, the nonlinearity, and entanglements of biomolecular data, ranging from DNA knots, RNA secondary structures, protein folding configurations, chromosomes, DNA origami, molecular assembly, to others at the macromolecular level, pose a severe challenge in their analysis and characterization. In the past few decades, mathematical concepts, models, algorithms, and tools from algebraic topology, combinatorial topology, computational topology, and topological data analysis, have demonstrated great power and begun to play an essential role in tackling the biomolecular data challenge. In this work, we introduce biomolecular topology, which concerns the topological problems and models originated from the biomolecular systems. More specifically, the biomolecular topology encompasses topological structures, properties and relations that are emerged from biomolecular structures, dynamics, interactions, and functions. We discuss the various types of biomolecular topology from structures (of proteins, DNAs, and RNAs), protein folding, and protein assembly. A brief discussion of databanks (and databases), theoretical models, and computational algorithms, is presented. Further, we systematically review related topological models, including graphs, simplicial complexes, persistent homology, persistent Laplacians, de Rham-Hodge theory, Yau-Hausdorff distance, and the topology-based machine learning models.

On the Classification of Polyhedral Links

Knots and links are ubiquitous in chemical systems. Their structure can be responsible for a variety of physical and chemical properties, making them very important in materials development. In this article, we analyze the topological structures of interlocking molecules composed of metal-peptide rings using the concept of polyhedral links. To that end, we discuss the topological classification of alternating polyhedral links.

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.

A statistical approach to knot confinement via persistent homology

In this paper, we study how randomly generated knots occupy a volume of space using topological methods. To this end, we consider the evolution of the first homology of an immersed metric neighbourhood of a knot's embedding for growing radii. Specifically, we extract features from the persistent homology (PH) of the Vietoris-Rips complexes built from point clouds associated with knots. Statistical analysis of our data shows the existence of increasing correlations between geometric quantities associated with the embedding and PH-based features, as a function of the knots' lengths. We further study the variation of these correlations for different knot types. Finally, this framework also allows us to define a simple notion of deviation from ideal configurations of knots.

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.

Keeping intracellular DNA untangled: A new role for condensin?

The DNA-passage activity of topoisomerase II accidentally produces DNA knots and interlinks within and between chromatin fibers. Fortunately, these unwanted DNA entanglements are actively removed by some mechanism. Here we present an outline on DNA knot formation and discuss recent studies that have investigated how intracellular DNA knots are removed. First, although topoisomerase II is able to minimize DNA entanglements in vitro to below equilibrium values, it is unclear whether such capacity performs equally in vivo in chromatinized DNA. Second, DNA supercoiling could bias topoisomerase II to untangle the DNA. However, experimental evidence indicates that transcriptional supercoiling of intracellular DNA boosts knot formation. Last, cohesin and condensin could tighten DNA entanglements via DNA loop extrusion (LE) and force their dissolution by topoisomerase II. Recent observations indicate that condensin activity promotes the removal of DNA knots during interphase and mitosis. This activity might facilitate the spatial organization and dynamics of chromatin.

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).

DNA Knot Malleability in Single-Digit Nanopores

Knots in long DNA molecules are prevalent in biological systems and serve as a model system for investigating static and dynamic properties of biopolymers. We explore the dynamics of knots in double-stranded DNA in a new regime of nanometer-scale confinement, large forces, and short time scales, using solid-state nanopores. We show that DNA knots undergo isomorphic translocation through a nanopore, retaining their equilibrium morphology by swiftly compressing in a lateral direction to fit the constriction. We observe no evidence of knot tightening or jamming, even for single-digit nanopores. We explain the observations as the malleability of DNA, characterized by sharp buckling of the DNA in nanopores, driven by the transient disruption of base pairing. Our molecular dynamics simulations support the model. These results are relevant not only for the understanding of DNA packing and manipulation in living cells but also for the polymer physics of DNA and the development of nanopore-based sequencing technologies.

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.

Time-dependent knotting of agitated chains

Agitated strings serve as macroscale models of spontaneous knotting, providing valuable insight into knotting dynamics at the microscale while allowing explicit analysis of the resulting knot topologies. We present an experimental setup for confined macroscale knot formation via tumbling along with a software interface to process complex knot data. Our setup allows characterization of knotting probability, knot complexity, and knot formation dynamics for knots with as many as 50 crossings. We find that the probability of knotting saturates below 80% within 100 s of the initiation of tumbling and that this saturation probability does not increase for chains above a critical length, an indication of nonequilibrium knot-formation conditions in our experiment. Despite the saturation in knot formation, we show that longer chains, while being more confined, will always tend to form knots of higher complexity since the free end can access a greater number of loops during tumbling.

Investigation of Ring and Star Polymers in Confined Geometries: Theory and Simulations

The calculations of the dimensionless layer monomer density profiles for a dilute solution of phantom ideal ring polymer chains and star polymers with f=4 arms in a Θ-solvent confined in a slit geometry of two parallel walls with repulsive surfaces and for the mixed case of one repulsive and the other inert surface were performed. Furthermore, taking into account the Derjaguin approximation, the dimensionless layer monomer density profiles for phantom ideal ring polymer chains and star polymers immersed in a solution of big colloidal particles with different adsorbing or repelling properties with respect to polymers were calculated. The density-force relation for the above-mentioned cases was analyzed, and the universal amplitude ratio B was obtained. Taking into account the small sphere expansion allowed obtaining the monomer density profiles for a dilute solution of phantom ideal ring polymers immersed in a solution of small spherical particles, or nano-particles of finite size, which are much smaller than the polymer size and the other characteristic mesoscopic length of the system. We performed molecular dynamics simulations of a dilute solution of linear, ring, and star-shaped polymers with N=300, 300 (360), and 1201 (4 × 300 + 1-star polymer with four arms) beads accordingly. The obtained analytical and numerical results for phantom ring and star polymers are compared with the results for linear polymer chains in confined geometries.

Condensin minimizes topoisomerase II-mediated entanglements of DNA in vivo

The juxtaposition of intracellular DNA segments, together with the DNA-passage activity of topoisomerase II, leads to the formation of DNA knots and interlinks, which jeopardize chromatin structure and gene expression. Recent studies in budding yeast have shown that some mechanism minimizes the knotting probability of intracellular DNA. Here, we tested whether this is achieved via the intrinsic capacity of topoisomerase II for simplifying the equilibrium topology of DNA; or whether it is mediated by SMC (structural maintenance of chromosomes) protein complexes like condensin or cohesin, whose capacity to extrude DNA loops could enforce dissolution of DNA knots by topoisomerase II. We show that the low knotting probability of DNA does not depend on the simplification capacity of topoisomerase II nor on the activities of cohesin or Smc5/6 complexes. However, inactivation of condensin increases the occurrence of DNA knots throughout the cell cycle. These results suggest an in vivo role for the DNA loop extrusion activity of condensin and may explain why condensin disruption produces a variety of alterations in interphase chromatin, in addition to persistent sister chromatid interlinks in mitotic chromatin.

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.

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.

Transcriptional supercoiling boosts topoisomerase II-mediated knotting of intracellular DNA

Recent studies have revealed that the DNA cross-inversion mechanism of topoisomerase II (topo II) not only removes DNA supercoils and DNA replication intertwines, but also produces small amounts of DNA knots within the clusters of nucleosomes that conform to eukaryotic chromatin. Here, we examine how transcriptional supercoiling of intracellular DNA affects the occurrence of these knots. We show that although (−) supercoiling does not change the basal DNA knotting probability, (+) supercoiling of DNA generated in front of the transcribing complexes increases DNA knot formation over 25-fold. The increase of topo II-mediated DNA knotting occurs both upon accumulation of (+) supercoiling in topoisomerase-deficient cells and during normal transcriptional supercoiling of DNA in TOP1 TOP2 cells. We also show that the high knotting probability (P^kn ≥ 0.5) of (+) supercoiled DNA reflects a 5-fold volume compaction of the nucleosomal fibers in vivo. Our findings indicate that topo II-mediated DNA knotting could be inherent to transcriptional supercoiling of DNA and other chromatin condensation processes and establish, therefore, a new crucial role of topoisomerase II in resetting the knotting–unknotting homeostasis of DNA during chromatin dynamics.

Quantitative disclosure of DNA knot chirality by high-resolution 2D-gel electrophoresis

The characterization of knots formed in duplex DNA has proved useful to infer biophysical properties and the spatial trajectory of DNA, both in free solution and across its macromolecular interactions. Since knotting, like supercoiling, makes DNA molecules more compact, DNA knot probability and knot complexity can be assessed by the electrophoretic velocity of nicked DNA circles. However, the chirality of the DNA knots has to be determined by visualizing the sign of their DNA crossings by means of electron microscopy. This procedure, which requires purifying the knotted DNA molecules and coating them with protein, is semi-quantitative and it is impracticable in biological samples that contain little amount of knotted DNA forms. Here, we took advantage of an earlier observation that the two chiral forms of a trefoil knot acquire slightly different electrophoretic velocity when the DNA is supercoiled. We introduced a second gel dimension to reveal these chiral forms in DNA mixtures that are largely unknotted. The result is a high-resolution 2D-gel electrophoresis procedure that quantitatively discerns the fractions of positive- and negative-noded trefoil knots formed in vitro and in vivo systems. This development in DNA knot analysis may uncover valuable information toward disclosing the architecture of DNA ensembles.

Complex DNA knots detected with a nanopore sensor

Equilibrium knots are common in biological polymers-their prevalence, size distribution, structure, and dynamics have been extensively studied, with implications to fundamental biological processes and DNA sequencing technologies. Nanopore microscopy is a high-throughput single-molecule technique capable of detecting the shape of biopolymers, including DNA knots. Here we demonstrate nanopore sensors that map the equilibrium structure of DNA knots, without spurious knot tightening and sliding. We show the occurrence of both tight and loose knots, reconciling previous contradictory results from different experimental techniques. We evidence the occurrence of two quantitatively different modes of knot translocation through the nanopores, involving very different tension forces. With large statistics, we explore the complex knots and, for the first time, reveal the existence of rare composite knots. We use parametrized complexity, in concert with simulations, to test the theoretical assumptions of the models, further asserting the relevance of nanopores in future investigation of knots.

Conformational State Hopping of Knots in Tensioned Polymer Chains

We use Brownian dynamics simulations to study the conformational states of knots on tensioned chains. Focusing specifically on the 8₁ knot, we observe knot conformational state hopping and show that the process can be described by a two-state kinetic model in the presence of an external force. The distribution of knot conformational states depends on the applied chain tension, which leads to a force-dependent distribution of knot untying pathways. We generalize our findings by considering the untying pathways of other knots and find that the way knots untie is generally governed by the force applied to the chain. From a broader perspective, being able to influence how a knot unties via external force can potentially be useful for applications of single-molecule techniques in which knots are unwanted.

Knot spectrum of turbulence

Streamlines, vortex lines and magnetic flux tubes in turbulent fluids and plasmas display a great amount of coiling, twisting and linking, raising the question as to whether their topological complexity (continually created and destroyed by reconnections) can be quantified. In superfluid helium, the discrete (quantized) nature of vorticity can be exploited to associate to each vortex loop a knot invariant called the Alexander polynomial whose degree characterizes the topology of that vortex loop. By numerically simulating the dynamics of a tangle of quantum vortex lines, we find that this quantum turbulence always contains vortex knots of very large degree which keep forming, vanishing and reforming, creating a distribution of topologies which we quantify in terms of a knot spectrum and its scaling law. We also find results analogous to those in the wider literature, demonstrating that the knotting probability of the vortex tangle grows with the vortex length, as for macromolecules, and saturates above a characteristic length, as found for tumbled strings.

DNA knots occur in intracellular chromatin

In vivo DNA molecules are narrowly folded within chromatin fibers and self-interacting chromatin domains. Therefore, intra-molecular DNA entanglements (knots) might occur via DNA strand passage activity of topoisomerase II. Here, we assessed the presence of such DNA knots in a variety of yeast circular minichromosomes. We found that small steady state fractions of DNA knots are common in intracellular chromatin. These knots occur irrespective of DNA replication and cell proliferation, though their abundance is reduced during DNA transcription. We found also that in vivo DNA knotting probability does not scale proportionately with chromatin length: it reaches a value of ∼0.025 in domains of ∼20 nucleosomes but tends to level off in longer chromatin fibers. These figures suggest that, while high flexibility of nucleosomal fibers and clustering of nearby nucleosomes facilitate DNA knotting locally, some mechanism minimizes the scaling of DNA knot formation throughout intracellular chromatin. We postulate that regulation of topoisomerase II activity and the fractal architecture of chromatin might be crucial to prevent a potentially massive and harmful self-entanglement of DNA molecules in vivo.

The Rabl configuration limits topological entanglement of chromosomes in budding yeast

The three dimensional organization of genomes remains mostly unknown due to their high degree of condensation. Biophysical studies predict that condensation promotes the topological entanglement of chromatin fibers and the inhibition of function. How organisms balance between functionally active genomes and a high degree of condensation remains to be determined. Here we hypothesize that the Rabl configuration, characterized by the attachment of centromeres and telomeres to the nuclear envelope, helps to reduce the topological entanglement of chromosomes. To test this hypothesis we developed a novel method to quantify chromosome entanglement complexity in 3D reconstructions obtained from Chromosome Conformation Capture (CCC) data. Applying this method to published data of the yeast genome, we show that computational models implementing the attachment of telomeres or centromeres alone are not sufficient to obtain the reduced entanglement complexity observed in 3D reconstructions. It is only when the centromeres and telomeres are attached to the nuclear envelope (i.e. the Rabl configuration) that the complexity of entanglement of the genome is comparable to that of the 3D reconstructions. We therefore suggest that the Rabl configuration is an essential player in the simplification of the entanglement of chromatin fibers.

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.

Calculation of π and Classification of Self-avoiding Lattices via DNA Configuration

Numerical simulation (e.g. Monte Carlo simulation) is an efficient computational algorithm establishing an integral part in science to understand complex physical and biological phenomena related with stochastic problems. Aside from the typical numerical simulation applications, studies calculating numerical constants in mathematics, and estimation of growth behavior via a non-conventional self-assembly in connection with DNA nanotechnology, open a novel perspective to DNA related to computational physics. Here, a method to calculate the numerical value of π, and way to evaluate possible paths of self-avoiding walk with the aid of Monte Carlo simulation, are addressed. Additionally, experimentally obtained variation of the π as functions of DNA concentration and the total number of trials, and the behaviour of self-avoiding random DNA lattice growth evaluated through number of growth steps, are discussed. From observing experimental calculations of π (π_exp) obtained by double crossover DNA lattices and DNA rings, fluctuation of π_exp tends to decrease as either DNA concentration or the number of trials increases. Based upon experimental data of self-avoiding random lattices grown by the three-point star DNA motifs, various lattice configurations are examined and analyzed. This new kind of study inculcates a novel perspective for DNA nanostructures related to computational physics and provides clues to solve analytically intractable problems.

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.

DNA Confined in a Nanodroplet: A Molecular Dynamics Study

As a major genetic material, the configuration and the mechanical properties of a double-stranded DNA (dsDNA) molecule in confinement are crucial for the application of nanotechnology and biological engineering. In the present paper, molecular dynamics simulation is utilized to study the configuration of dsDNA in a nanodroplet on a graphene substrate. The results show that the semiflexible dsDNA molecule changes its configuration with radius of gyration ( R_g) of a few nanometers because of the confined space, that is, the R_g of the dsDNA molecule decreases with the reduction of the nanodroplet size. In comparison, the dsDNA in the bulk usually has a persistent length of tens of nanometers. Especially, if the nanodroplet is small enough, the dsDNA molecule might form a loop structure inside. The dsDNA molecule affects the wetting properties of the graphene substrate. It is found that the graphene becomes more hydrophilic in smaller systems containing the dsDNA molecule, whereas for larger droplets, the changes of the contact angles are not significant with the presence of dsDNA. Moreover, the results indicate that for larger droplets, the line tension of the droplet containing DNA is positive and greater than that without DNA; for smaller droplets, the line tension becomes negative because the dsDNA is compressed and bent in the confinement, and has the potential to expand outwards. The worm-like chain model is used to study the bending energy of a dsDNA molecule in a droplet. The results address that the bending energy of the non-loop-structured dsDNA decreases as the droplet becomes larger, and it is larger than that of loop-structured dsDNA, as the loop structure efficiently prevents the DNA from bending in the vertical direction.

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.

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.

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 .