Multiscale equilibration of highly entangled isotropic model polymer melts.

We present a computationally efficient multiscale method for preparing equilibrated, isotropic long-chain model polymer melts. As an application, we generate Kremer-Grest melts of 1000 chains with 200 entanglements and 25 000-2000 beads/chain, which cover the experimentally relevant bending rigidities up to and beyond the limit of the isotropic-nematic transition. In the first step, we employ Monte Carlo simulations of a lattice model to equilibrate the large-scale chain structure above the tube scale while ensuring a spatially homogeneous density distribution. We then use theoretical insight from a constrained mode tube model to introduce the bead degrees of freedom together with random walk conformational statistics all the way down to the Kuhn scale of the chains. This is followed by a sequence of simulations with carefully parameterized force-capped bead-spring models, which slowly introduce the local bead packing while reproducing the larger-scale chain statistics of the target Kremer-Grest system at all levels of force-capping. Finally, we can switch to the full Kremer-Grest model without perturbing the structure. The resulting chain statistics is in excellent agreement with literature results on all length scales accessible in brute-force simulations of shorter chains.

Single-molecule stretching experiments of flexible (wormlike) chain molecules in different ensembles: Theory and a potential application of finite chain length effects to nick-counting in DNA.

We propose a formalism for deriving force-elongation and elongation-force relations for flexible chain molecules from analytical expressions for their radial distribution function, which provides insight into the factors controlling the asymptotic behavior and finite chain length corrections. In particular, we apply this formalism to our previously developed interpolation formula for the wormlike chain end-to-end distance distribution. The resulting expression for the asymptotic limit of infinite chain length is of similar quality to the numerical evaluation of Marko and Siggia's variational theory and considerably more precise than their interpolation formula. A comparison to numerical data suggests that our analytical finite chain length corrections achieve a comparable accuracy. As an application of our results, we discuss the possibility of inferring the time-dependent number of nicks in single-molecule stretching experiments on double-stranded DNA from the accompanying changes in the effective chain length.

Local loop opening in untangled ring polymer melts: a detailed "Feynman test" of models for the large scale structure.

The conformational statistics of ring polymers in melts or dense solutions is strongly affected by their quenched microscopic topological state. The effect is particularly strong for untangled (i.e. non-concatenated and unknotted) rings, which are known to crumple and segregate. Here we study these systems using a computationally efficient multi-scale approach, where we combine massive simulations on the fiber level with the explicit construction of untangled ring melt configurations based on theoretical ideas for their large scale structure. We find (i) that topological constraints may be neglected on scales below the standard entanglement length, Le, (ii) that rings with a size 1 ≤ Lr/Le ≤ 30 exhibit nearly ideal lattice tree behavior characterized by primitive paths which are randomly branched on the entanglement scale, and (iii) that larger rings are compact with gyration radii Rg2(Lr) ∝ Lr2/3. The detailed comparison between equilibrated and constructed ensembles allows us to perform a "Feynman test" of our understanding of untangled rings: can we convert ideas for the large scale ring structure into algorithms for constructing (nearly) equilibrated ring melt samples? We show that most structural observables are quantitatively reproduced by two different construction schemes: hierarchical crumpling and ring melts derived from the analogy to interacting branched polymers. However, the latter fail the "Feynman test" with respect to the magnetic radius, Rm, which we have defined based on an analogy to magnetostatics. While Rm is expected to vanish for double-folded structures, the observed values of Rm2(Lr) ∝ Rg2(Lr) provide a simple and computationally convenient measure of the presence of a non-negligible amount of local loop opening in crumpled rings.

Conformational statistics of randomly branching double-folded ring polymers.

The conformations of topologically constrained double-folded ring polymers can be described as wrappings of randomly branched primitive trees. We extend previous work on the tree statistics under different (solvent) conditions to explore the conformational statistics of double-folded rings in the limit of tight wrapping. In particular, we relate the exponents characterizing the ring statistics to those describing the primitive trees and discuss the distribution functions [Formula: see text] and [Formula: see text] for the spatial distance, [Formula: see text], and tree contour distance, L, between monomers as a function of their ring contour distance, [Formula: see text].

Flory theory of randomly branched polymers.

Randomly branched polymer chains (or trees) are a classical subject of polymer physics with connections to the theory of magnetic systems, percolation and critical phenomena. More recently, the model has been reconsidered for RNA, supercoiled DNA and the crumpling of topologically-constrained polymers. While solvable in the ideal case, little is known exactly about randomly branched polymers with volume interactions. Flory theory provides a simple, unifying description for a wide range of branched systems, including isolated trees in good and θ-solvent, and tree melts. In particular, the approach provides a common framework for the description of randomly branched polymers with quenched connectivity and for randomly branching polymers with annealed connectivity. Here we review the Flory theory for interacting trees in the asymptotic limit of very large polymerization degree for good solvent, θ-solutions and melts, and report its predictions for annealed connectivity in θ-solvents. We compare the predictions of Flory theory for randomly branched polymers to a wide range of available analytical and numerical results and conclude that they are qualitatively excellent and quantitatively good in most cases.

Computer simulations of melts of randomly branching polymers.

Randomly branching polymers with annealed connectivity are model systems for ring polymers and chromosomes. In this context, the branched structure represents transient folding induced by topological constraints. Here we present computer simulations of melts of annealed randomly branching polymers of 3 ≤ N ≤ 1800 segments in d = 2 and d = 3 dimensions. In all cases, we perform a detailed analysis of the observed tree connectivities and spatial conformations. Our results are in excellent agreement with an asymptotic scaling of the average tree size of R ∼ N1/d, suggesting that the trees behave as compact, territorial fractals. The observed swelling relative to the size of ideal trees, R ∼ N1/4, demonstrates that excluded volume interactions are only partially screened in melts of annealed trees. Overall, our results are in good qualitative agreement with the predictions of Flory theory. In particular, we find that the trees swell by the combination of modified branching and path stretching. However, the former effect is subdominant and difficult to detect in d = 3 dimensions.

Ring polymers in the melt state: the physics of crumpling.

The conformational statistics of ring polymers in melts or dense solutions is strongly affected by their quenched microscopic topological state. The effect is particularly strong for nonconcatenated unknotted rings, which are known to crumple and segregate and which have been implicated as models for the generic behavior of interphase chromosomes. Here we use a computationally efficient multiscale approach to show that melts of rings of total contour length Lr can be quantitatively mapped onto melts of interacting lattice trees with gyration radii ⟨R(g)(2)(Lr)⟩ ∝ L(r)(2ν) and ν = 0.32 ± 0.01.

Temperature dependence of the DNA double helix at the nanoscale: structure, elasticity, and fluctuations.

Biological organisms exist over a broad temperature range of -15°C to +120°C, where many molecular processes involving DNA depend on the nanoscale properties of the double helix. Here, we present results of extensive molecular dynamics simulations of DNA oligomers at different temperatures. We show that internal basepair conformations are strongly temperature-dependent, particularly in the stretch and opening degrees of freedom whose harmonic fluctuations can be considered the initial steps of the DNA melting pathway. The basepair step elasticity contains a weaker, but detectable, entropic contribution in the roll, tilt, and rise degrees of freedom. To extend the validity of our results to the temperature interval beyond the standard melting transition relevant to extremophiles, we estimate the effects of superhelical stress on the stability of the basepair steps, as computed from the Benham model. We predict that although the average twist decreases with temperature in vitro, the stabilizing external torque in vivo results in an increase of ∼1°/bp (or a superhelical density of Δσ ≃ +0.03) in the interval 0-100°C. In the final step, we show that the experimentally observed apparent bending persistence length of torsionally unconstrained DNA can be calculated from a hybrid model that accounts for the softening of the double helix and the presence of transient denaturation bubbles. Although the latter dominate the behavior close to the melting transition, the inclusion of helix softening is important around standard physiological temperatures.

Hamiltonian adaptive resolution simulation for molecular liquids.

Adaptive resolution schemes allow the simulation of a molecular fluid treating simultaneously different subregions of the system at different levels of resolution. In this work we present a new scheme formulated in terms of a global Hamiltonian. Within this approach equilibrium states corresponding to well-defined statistical ensembles can be generated making use of all standard molecular dynamics or Monte Carlo methods. Models at different resolutions can thus be coupled, and thermodynamic equilibrium can be modulated keeping each region at desired pressure or density without disrupting the Hamiltonian framework.

Monte carlo adaptive resolution simulation of multicomponent molecular liquids.

Complex soft matter systems can be efficiently studied with the help of adaptive resolution simulation methods, concurrently employing two levels of resolution in different regions of the simulation domain. The nonmatching properties of high- and low-resolution models, however, lead to thermodynamic imbalances between the system's subdomains. Such inhomogeneities can be healed by appropriate compensation forces, whose calculation requires nontrivial iterative procedures. In this work we employ the recently developed Hamiltonian adaptive resolution simulation method to perform Monte Carlo simulations of a binary mixture, and propose an efficient scheme, based on Kirkwood thermodynamic integration, to regulate the thermodynamic balance of multicomponent systems.

From crystal and NMR structures, footprints and cryo-electron-micrographs to large and soft structures: nanoscale modeling of the nucleosomal stem.

The interaction of histone H1 with linker DNA results in the formation of the nucleosomal stem structure, with considerable influence on chromatin organization. In a recent paper [Syed,S.H., Goutte-Gattat,D., Becker,N., Meyer,S., Shukla,M.S., Hayes,J.J., Everaers,R., Angelov,D., Bednar,J. and Dimitrov,S. (2010) Single-base resolution mapping of H1-nucleosome interactions and 3D organization of the nucleosome. Proc. Natl Acad. Sci. USA, 107, 9620-9625], we published results of biochemical footprinting and cryo-electron-micrographs of reconstituted mono-, di- and tri-nucleosomes, for H1 variants with different lengths of the cationic C-terminus. Here, we present a detailed account of the analysis of the experimental data and we include thermal fluctuations into our nano-scale model of the stem structure. By combining (i) crystal and NMR structures of the nucleosome core particle and H1, (ii) the known nano-scale structure and elasticity of DNA, (iii) footprinting information on the location of protected sites on the DNA backbone and (iv) cryo-electron micrographs of reconstituted tri-nucleosomes, we arrive at a description of a polymorphic, hierarchically organized stem with a typical length of 20 ± 2 base pairs. A comparison to linker conformations inferred for poly-601 fibers with different linker lengths suggests, that intra-stem interactions stabilize and facilitate the formation of dense chromatin fibers.

Single-base resolution mapping of H1-nucleosome interactions and 3D organization of the nucleosome.

Despite the key role of the linker histone H1 in chromatin structure and dynamics, its location and interactions with nucleosomal DNA have not been elucidated. In this work we have used a combination of electron cryomicroscopy, hydroxyl radical footprinting, and nanoscale modeling to analyze the structure of precisely positioned mono-, di-, and trinucleosomes containing physiologically assembled full-length histone H1 or truncated mutants of this protein. Single-base resolution *OH footprinting shows that the globular domain of histone H1 (GH1) interacts with the DNA minor groove located at the center of the nucleosome and contacts a 10-bp region of DNA localized symmetrically with respect to the nucleosomal dyad. In addition, GH1 interacts with and organizes about one helical turn of DNA in each linker region of the nucleosome. We also find that a seven amino acid residue region (121-127) in the COOH terminus of histone H1 was required for the formation of the stem structure of the linker DNA. A molecular model on the basis of these data and coarse-grain DNA mechanics provides novel insights on how the different domains of H1 interact with the nucleosome and predicts a specific H1-mediated stem structure within linker DNA.

DNA supercoiling in bacteria: state of play and challenges from a viewpoint of physics based modeling.

DNA supercoiling is central to many fundamental processes of living organisms. Its average level along the chromosome and over time reflects the dynamic equilibrium of opposite activities of topoisomerases, which are required to relax mechanical stresses that are inevitably produced during DNA replication and gene transcription. Supercoiling affects all scales of the spatio-temporal organization of bacterial DNA, from the base pair to the large scale chromosome conformation. Highlighted in vitro and in vivo in the 1960s and 1970s, respectively, the first physical models were proposed concomitantly in order to predict the deformation properties of the double helix. About fifteen years later, polymer physics models demonstrated on larger scales the plectonemic nature and the tree-like organization of supercoiled DNA. Since then, many works have tried to establish a better understanding of the multiple structuring and physiological properties of bacterial DNA in thermodynamic equilibrium and far from equilibrium. The purpose of this essay is to address upcoming challenges by thoroughly exploring the relevance, predictive capacity, and limitations of current physical models, with a specific focus on structural properties beyond the scale of the double helix. We discuss more particularly the problem of DNA conformations, the interplay between DNA supercoiling with gene transcription and DNA replication, its role on nucleoid formation and, finally, the problem of scaling up models. Our primary objective is to foster increased collaboration between physicists and biologists. To achieve this, we have reduced the respective jargon to a minimum and we provide some explanatory background material for the two communities.

Dynamics of fluid bilayer vesicles: Soft meshes and robust curvature energy discretization.

Continuum models like the Helfrich Hamiltonian are widely used to describe fluid bilayer vesicles. Here we study the molecular dynamics compatible dynamics of the vertices of two-dimensional meshes representing the bilayer, whose in-plane motion is only weakly constrained. We show (i) that Jülicher's discretization of the curvature energy offers vastly superior robustness for soft meshes compared to the commonly employed expression by Gommper and Kroll and (ii) that for sufficiently soft meshes, the typical behavior of fluid bilayer vesicles can emerge even if the mesh connectivity remains fixed throughout the simulations. In particular, soft meshes can accommodate large shape transformations, and the model can generate the typical ℓ^{-4} signal for the amplitude of surface undulation modes of nearly spherical vesicles all the way up to the longest wavelength modes. Furthermore, we compare results for Newtonian, Langevin, and Brownian dynamics simulations of the mesh vertices to demonstrate that the internal friction of the membrane model is negligible, making it suitable for studying the internal dynamics of vesicles via coupling to hydrodynamic solvers or particle-based solvent models.

Multi-scale modeling of diffusion-controlled reactions in polymers: renormalisation of reactivity parameters.

The quantitative description of polymeric systems requires hierarchical modeling schemes, which bridge the gap between the atomic scale, relevant to chemical or biomolecular reactions, and the macromolecular scale, where the longest relaxation modes occur. Here, we use the formalism for diffusion-controlled reactions in polymers developed by Wilemski, Fixman, and Doi to discuss the renormalisation of the reactivity parameters in polymer models with varying spatial resolution. In particular, we show that the adjustments are independent of chain length. As a consequence, it is possible to match reactions times between descriptions with different resolution for relatively short reference chains and to use the coarse-grained model to make quantitative predictions for longer chains. We illustrate our results by a detailed discussion of the classical problem of chain cyclization in the Rouse model, which offers the simplest example of a multi-scale descriptions, if we consider differently discretized Rouse models for the same physical system. Moreover, we are able to explore different combinations of compact and non-compact diffusion in the local and large-scale dynamics by varying the embedding dimension.

DNA nanomechanics in the nucleosome.

The sequence-dependent mechanics of DNA-histone binding are essential for nucleosome positioning and mobility. Here we reanalyze nucleosome crystal structures in terms of the well-characterized base-pair scale DNA elasticity, extracting the forces and torques acting on all bound DNA base pairs. We find that the strongest forces follow a characteristic repeating pattern that recovers the 12 known DNA backbone-histone contact sites. DNA twist defects and histone point mutations modify this pattern in interpretable ways. Additional, irregular forces between contact sites reveal histone tail-DNA interactions, whereas requiring the absence of external forces leads to a structural refinement of linker DNA. Based on these observations, we propose a simple, structure-based mechanical model of the nucleosome that is able to explain the placement of DNA twist defects in 146 base-pair nucleosomes and allows an estimate of the elastic energy spectrum of nucleosome twist defect states.

Monte Carlo simulation of a lattice model for the dynamics of randomly branching double-folded ring polymers.

Supercoiled DNA, crumpled interphase chromosomes, and topologically constrained ring polymers often adopt treelike, double-folded, randomly branching configurations. Here we study an elastic lattice model for tightly double-folded ring polymers, which allows for the spontaneous creation and deletion of side branches coupled to a diffusive mass transport, which is local both in space and on the connectivity graph of the tree. We use Monte Carlo simulations to study systems falling into three different universality classes: ideal double-folded rings without excluded volume interactions, self-avoiding double-folded rings, and double-folded rings in the melt state. The observed static properties are in good agreement with exact results, simulations, and predictions of Flory theory for randomly branching polymers. For example, in the melt state rings adopt compact configurations and exhibit territorial behavior. In particular, we show that the emergent dynamics is in excellent agreement with a recent scaling theory and illustrate the qualitative differences with the familiar reptation dynamics of linear chains.

Looping probabilities in model interphase chromosomes.

Fluorescence in-situ hybridization (FISH) and chromosome conformation capture (3C) are two powerful techniques for investigating the three-dimensional organization of the genome in interphase nuclei. The use of these techniques provides complementary information on average spatial distances (FISH) and contact probabilities (3C) for specific genomic sites. To infer the structure of the chromatin fiber or to distinguish functional interactions from random colocalization, it is useful to compare experimental data to predictions from statistical fiber models. The current estimates of the fiber stiffness derived from FISH and 3C differ by a factor of 5. They are based on the wormlike chain model and a heuristic modification of the Shimada-Yamakawa theory of looping for unkinkable, unconstrained, zero-diameter filaments. Here, we provide an extended theoretical and computational framework to explain the currently available experimental data for various species on the basis of a unique, minimal model of decondensing chromosomes: a kinkable, topologically constraint, semiflexible polymer with the (FISH) Kuhn length of l(K) = 300 nm, 10 kinks per Mbp, and a contact distance of 45 nm. In particular: 1), we reconsider looping of finite-diameter filaments on the basis of an analytical approximation (novel, to our knowledge) of the wormlike chain radial density and show that unphysically large contact radii would be required to explain the 3C data based on the FISH estimate of the fiber stiffness; 2), we demonstrate that the observed interaction frequencies at short genomic lengths can be explained by the presence of a low concentration of curvature defects (kinks); and 3), we show that the most recent experimental 3C data for human chromosomes are in quantitative agreement with interaction frequencies extracted from our simulations of topologically confined model chromosomes.

Stress relaxation in entangled polymer melts.

We present an extensive set of simulation results for the stress relaxation in equilibrium and step-strained bead-spring polymer melts. The data allow us to explore the chain dynamics and the shear relaxation modulus, G(t), into the plateau regime for chains with Z=40 entanglements and into the terminal relaxation regime for Z=10. Using the known (Rouse) mobility of unentangled chains and the melt entanglement length determined via the primitive path analysis of the microscopic topological state of our systems, we have performed parameter-free tests of several different tube models. We find excellent agreement for the Likhtman-McLeish theory using the double reptation approximation for constraint release, if we remove the contribution of high-frequency modes to contour length fluctuations of the primitive chain.

Prediction of RNA multiloop and pseudoknot conformations from a lattice-based, coarse-grain tertiary structure model.

We present a semiquantitative lattice model of RNA folding, which is able to reproduce complex folded structures such as multiloops and pseudoknots without relying on the frequently employed ad hoc generalization of the Jacobson-Stockmayer loop entropy. We derive the model parameters from the Turner description of simple secondary structural elements and pay particular attention to the unification of mismatch and coaxial stacking parameters as well as of border and nonlocal loop parameters, resulting in a reduced, unified parameter set for simple loops of arbitrary type and size. For elementary structures, the predictive power of the model is comparable to the standard secondary structure approaches, from which its parameters are derived. For complex structures, our approach offers a systematic treatment of generic effects of chain connectivity as well as of excluded volume or attractive interactions between and within all elements of the secondary structure. We reproduce the native structures of tRNA multiloops and of viral frameshift signal pseudoknots.