Stretching semiflexible polymer chains: evidence for the importance of excluded volume effects from Monte Carlo simulation

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.

Generating Chromosome Geometries in a Minimal Cell From Cryo-Electron Tomograms and Chromosome Conformation Capture Maps

JCVI-syn3A is a genetically minimal bacterial cell, consisting of 493 genes and only a single 543 kbp circular chromosome. Syn3A’s genome and physical size are approximately one-tenth those of the model bacterial organism Escherichia coli’s, and the corresponding reduction in complexity and scale provides a unique opportunity for whole-cell modeling. Previous work established genome-scale gene essentiality and proteomics data along with its essential metabolic network and a kinetic model of genetic information processing. In addition to that information, whole-cell, spatially-resolved kinetic models require cellular architecture, including spatial distributions of ribosomes and the circular chromosome’s configuration. We reconstruct cellular architectures of Syn3A cells at the single-cell level directly from cryo-electron tomograms, including the ribosome distributions. We present a method of generating self-avoiding circular chromosome configurations in a lattice model with a resolution of 11.8 bp per monomer on a 4 nm cubic lattice. Realizations of the chromosome configurations are constrained by the ribosomes and geometry reconstructed from the tomograms and include DNA loops suggested by experimental chromosome conformation capture (3C) maps. Using ensembles of simulated chromosome configurations we predict chromosome contact maps for Syn3A cells at resolutions of 250 bp and greater and compare them to the experimental maps. Additionally, the spatial distributions of ribosomes and the DNA-crowding resulting from the individual chromosome configurations can be used to identify macromolecular structures formed from ribosomes and DNA, such as polysomes and expressomes.

Elasticity of a DNA chain dotted with bubbles under force

The flexibility and the extension along the direction of the force are shown to be related to the bubble number fluctuation and the average number of bubbles, respectively, when the strands of the DNA are subjected to a force along the same direction, here called a stretching force. The force-temperature phase diagram shows the existence of a tricritical point, where the first-order force-induced zipping transition becomes continuous. On the other hand, when the forces are being applied in opposite directions, here called an unzipping force, the transition remains first order, with the possibility of vanishing of the low-temperature reentrant phase for a semiflexible DNA. Moreover, we found that the bulk elasticity changes only if an external force penetrates the bound phase and affects the bubble states.

Conformational behavior and self-assembly of disjoint semi-flexible ring polymers adsorbed on solid substrates

The conformational behavior and spatial organization of self-avoiding semi-flexible ring polymers, that are fully adsorbed on solid substrates, are investigated via systematic coarse-grained molecular dynamics simulations. Our results show that both conformations and spatial organization of the polymers depend strongly on their bending stiffness, κ, and on their areal number density, ρ. For ρ < ρ*, where ρ* is the overlap density, and for low values of κ, thermal fluctuations lead to weakly anisotropic instantaneous conformations of the polymers. The interplay between thermal fluctuations and polymer stiffness leads to a non-monotonic dependence of the polymers elongation on κ with a maximum elongation at some intermediate κ. Regardless of κ, the polymers elongation is almost independent of ρ for ρ ⪅ ρ*, then increases with ρ. At ρ ≈ ρ* and high κ, the almost circularly-shaped polymers self-assemble into a triangular lattice with quasi-long range order. For ρ above ρ* and high κ, crowding of the polymers leads to their self-assembly into liquid-crystalline phases. In particular, for ρ moderately above ρ* and high κ, the polymer conformations are obround and self-assemble into domains with smectic-A-like order. At higher densities, the polymer have a biconcave geometry and self-assemble into domains with smectic-C-like order.

The Persistence Length of Semiflexible Polymers in Lattice Monte Carlo Simulations

While applying computer simulations to study semiflexible polymers, it is a primary task to determine the persistence length that characterizes the chain stiffness. One frequently asked question concerns the relationship between persistence length and the bending constant of applied bending potential. In this paper, theoretical persistence lengths of polymers with two different bending potentials were analyzed and examined by using lattice Monte Carlo simulations. We found that the persistence length was consistent with theoretical predictions only in bond fluctuation model with cosine squared angle potential. The reason for this is that the theoretical persistence length is calculated according to a continuous bond angle, which is discrete in lattice simulations. In lattice simulations, the theoretical persistence length is larger than that in continuous simulations.

Rigidity-induced scale invariance in polymer ejection from capsid

While the dynamics of a fully flexible polymer ejecting a capsid through a nanopore has been extensively studied, the ejection dynamics of semiflexible polymers has not been properly characterized. Here we report results from simulations of ejection dynamics of semiflexible polymers ejecting from spherical capsids. Ejections start from strongly confined polymer conformations of constant initial monomer density. We find that, unlike for fully flexible polymers, for semiflexible polymers the force measured at the pore does not show a direct relation to the instantaneous ejection velocity. The cumulative waiting time t(s), that is, the time at which a monomer s exits the capsid the last time, shows a clear change when increasing the polymer rigidity κ. The major part of an ejecting polymer is driven out of the capsid by internal pressure. At the final stage the polymer escapes the capsid by diffusion. For the driven part there is a crossover from essentially exponential growth of t with s of the fully flexible polymers to a scale-invariant form. In addition, a clear dependence of t on polymer length N_{0} was found. These findings combined give the dependence t(s)∝N_{0}^{0.55}s^{1.33} for the strongly rigid polymers. This crossover in dynamics where κ acts as a control parameter is reminiscent of a phase transition. This analogy is further enhanced by our finding a perfect data collapse of t for polymers of different N_{0} and any constant κ.

Optimal Reactivity and Improved Self-Healing Capability of Structurally Dynamic Polymers Grafted on Janus Nanoparticles Governed by Chain Stiffness and Spatial Organization

Structurally dynamic polymers are recognized as a key potential to revolutionize technologies ranging from design of self-healing materials to numerous biomedical applications. Despite intense research in this area, optimizing reactivity and thereby improving self-healing ability at the most fundamental level pose urgent issue for wider applications of such emerging materials. Here, the authors report the first mechanistic investigation of the fundamental principle for the dependence of reactivity and self-healing capabilities on the properties inherent to dynamic polymers by combining large-scale computer simulation, theoretical analysis, and experimental discussion. The results allow to reveal how chain stiffness and spatial organization regulate reactivity of dynamic polymers grafted on Janus nanoparticles and mechanically mediated reaction in their reverse chemistry, and, particularly, identify that semiflexible dynamic polymers possess the optimal reactivity and self-healing ability. The authors also develop an analytical model of blob theory of polymer chains to complement the simulation results and reveal essential scaling laws for optimal reactivity. The findings offer new insights into the physical mechanism in various systems involving reverse/dynamic chemistry. These studies highlight molecular engineering of polymer architecture and intrinsic property as a versatile strategy in control over the structural responses and functionalities of emerging materials with optimized self-healing capabilities.

Semiflexible polymers confined in a slit pore with attractive walls: two-dimensional liquid crystalline order versus capillary nematization

Semiflexible polymers under good solvent conditions interacting with attractive planar surfaces are investigated by Molecular Dynamics (MD) simulations and classical Density Functional Theory (DFT). A bead-spring type potential complemented by a bending potential is used, allowing variation of chain stiffness from completely flexible coils to rod-like polymers whose persistence length by far exceeds their contour length. Solvent is only implicitly included, monomer-monomer interactions being purely repulsive, while two types of attractive wall-monomer interactions are considered: (i) a strongly attractive Mie-type potential, appropriate for a strictly structureless wall, and (ii) a corrugated wall formed by Lennard-Jones particles arranged on a square lattice. It is found that in dilute solutions the former case leads to the formation of a strongly adsorbed surface layer, and the profile of density and orientational order in the z-direction perpendicular to the wall is predicted by DFT in nice agreement with MD. While for very low bulk densities a Kosterlitz-Thouless type transition from the isotropic phase to a phase with power-law decay of nematic correlations is suggested to occur in the strongly adsorbed layer, for larger densities a smectic-C phase in the surface layer is detected. No "capillary nematization" effect at higher bulk densities is found in this system, unlike systems with repulsive walls. This finding is attributed to the reduction of the bulk density (in the center of the slit pore) due to polymer adsorption on the attractive wall, for a system studied in the canonical ensemble. Consequently in a system with two attractive walls nematic order in the slit pore can occur only at a higher density than for a bulk system.

Stiffer double-stranded DNA in two-dimensional confinement due to bending anisotropy

Using analytical approach and Monte Carlo (MC) simulations, we study the elastic behavior of the intrinsically twisted elastic ribbons with bending anisotropy, such as double-stranded DNA (dsDNA), in two-dimensional (2D) confinement. We show that, due to the bending anisotropy, the persistence length of dsDNA in 2D conformations is always greater than three-dimensional (3D) conformations. This result is in consistence with the measured values for DNA persistence length in 2D and 3D in equal biological conditions. We also show that in two dimensions, an anisotropic, intrinsically twisted polymer exhibits an implicit twist-bend coupling, which leads to the transient curvature increasing with a half helical turn periodicity along the bent polymer.

Semiflexible macromolecules in quasi-one-dimensional confinement: Discrete versus continuous bond angles

The conformations of semiflexible polymers in two dimensions confined in a strip of width D are studied by computer simulations, investigating two different models for the mechanism by which chain stiffness is realized. One model (studied by molecular dynamics) is a bead-spring model in the continuum, where stiffness is controlled by a bond angle potential allowing for arbitrary bond angles. The other model (studied by Monte Carlo) is a self-avoiding walk chain on the square lattice, where only discrete bond angles (0° and ±90°) are possible, and the bond angle potential then controls the density of kinks along the chain contour. The first model is a crude description of DNA-like biopolymers, while the second model (roughly) describes synthetic polymers like alkane chains. It is first demonstrated that in the bulk the crossover from rods to self-avoiding walks for both models is very similar, when one studies average chain linear dimensions, transverse fluctuations, etc., despite their differences in local conformations. However, in quasi-one-dimensional confinement two significant differences between both models occur: (i) The persistence length (extracted from the average cosine of the bond angle) gets renormalized for the lattice model when D gets less than the bulk persistence length, while in the continuum model it stays unchanged. (ii) The monomer density near the repulsive walls for semiflexible polymers is compatible with a power law predicted for the Kratky-Porod model in the case of the bead-spring model, while for the lattice case it tends to a nonzero constant across the strip. However, for the density of chain ends, such a constant behavior seems to occur for both models, unlike the power law observed for flexible polymers. In the regime where the bulk persistence length ℓp is comparable to D, hairpin conformations are detected, and the chain linear dimensions are discussed in terms of a crossover from the Daoud/De Gennes "string of blobs"-picture to the flexible rod picture when D decreases and/or the chain stiffness increases. Introducing a suitable further coarse-graining of the chain contours of the continuum model, direct estimates for the deflection length and its distribution could be obtained.

Scattering and Gaussian Fluctuation Theory for Semiflexible Polymers

The worm-like chain is one of the best theoretical models of the semiflexible polymer. The structure factor, which can be obtained by scattering experiment, characterizes the density correlation in different length scales. In the present review, the numerical method to compute the static structure factor of the worm-like chain model and its general properties are demonstrated. Especially, the chain length and persistence length involved multi-scale nature of the worm-like chain model are well discussed. Using the numerical structure factor, Gaussian fluctuation theory of the worm-like chain model can be developed, which is a powerful tool to analyze the structure stability and to predict the spinodal line of the system. The microphase separation of the worm-like diblock copolymer is considered as an example to demonstrate the usage of Gaussian fluctuation theory.

Conformational Properties of Active Semiflexible Polymers

The conformational properties of flexible and semiflexible polymers exposed to active noise are studied theoretically. The noise may originate from the interaction of the polymer with surrounding active (Brownian) particles or from the inherent motion of the polymer itself, which may be composed of active Brownian particles. In the latter case, the respective monomers are independently propelled in directions changing diffusively. For the description of the polymer, we adopt the continuous Gaussian semiflexible polymer model. Specifically, the finite polymer extensibility is taken into account, which turns out to be essential for the polymer conformations. Our analytical calculations predict a strong dependence of the relaxation times on the activity. In particular, semiflexible polymers exhibit a crossover from a bending elasticity-dominated dynamics to the flexible polymer dynamics with increasing activity. This leads to a significant activity-induced polymer shrinkage over a large range of self-propulsion velocities. For large activities, the polymers swell and their extension becomes comparable to the contour length. The scaling properties of the mean square end-to-end distance with respect to the polymer length and monomer activity are discussed.

Effect of Uniformly Applied Force and Molecular Characteristics of a Polymer Chain on Its Adhesion to Graphene Substrates

The force-induced desorption of a polymer chain from a graphene substrate is studied with molecular dynamics (MD). A critical force needs to be exceeded before detachment of the polymer from the substrate. It is found that for a chain to exhibit good adhesive properties the chain configuration should consist of fibrils-elongated, aligned sections of polymers and cavities which dissipate the applied energy. A fibrillation index is defined to quantify the quality of fibrils. We focus on the molecular properties of the polymer chain, which can lead to large amounts of fibrillation, and find that both strong attraction between the polymer and substrate and good solvency conditions are important conditions for this. We also vary the stiffness of the chain and find that for less stiff chains a plateau in the stress-strain curve gives rise to good adhesion however for very stiff chains there is limited elongation of the chain but the chain can still exhibit good fibrillation by a lamella-like rearrangement. Finally, it is found that the detachment time, t, of a polymer from the adsorbed substrate is inversely proportional to force, F (i.e., t ∝ F(-γ)), where exponent γ depends on the solvent quality, polymer-substrate attraction, and chain stiffness.

Effect of excluded volume on the force-extension of wormlike chains in slit confinement

We use pruned-enriched Rosenbluth method simulations to develop a quantitative phase diagram for the stretching of a real wormlike chain confined in a slit. Our simulations confirm the existence of a "confined Pincus" regime in slit confinement, analogous to the Pincus regime in free solution, where excluded volume effects are sensible. The lower bound for the confined Pincus regime in the force-molecular weight plane, as well as the scaling of the extension with force and slit size, agree with an existing scaling theory for this regime. The upper bound of the confined Pincus regime depends on the strength of the confinement. For strong confinement, the confined Pincus regime ends when the contour length in the Pincus blob is too short to have intrablob excluded volume. As a result, the chain statistics become ideal and the confined Pincus regime at low forces is connected directly to ideal chain stretching at large forces. In contrast, for weak confinement, the confined Pincus regime ends when the Pincus blobs no longer fit inside the slit, even though there is sufficient contour length to have excluded volume inside the Pincus blob. As a result, weak confinement leads to a free-solution Pincus regime intervening between the confined Pincus regime for weak forces and ideal chain stretching at strong forces. Our results highlight shortcomings in existing models for the stretching of wormlike chains in slits.

Chain stiffness regulates entropy-templated perfect mixing at single-nanoparticle level

The mixing on a single-particle level of chemically incompatible nanoparticles is an outstanding challenge for many applications. Burgeoning research activity suggests that entropic templating is a potential strategy to address this issue. Herein, using systematic computer simulations of model nanoparticle systems, we show that the entropy-templated interfacial organization of nanoparticles significantly depends on the stiffness of tethered chains. Unexpectedly, the optimal chain stiffness can be identified wherein a system exhibits the most perfect mixing for a certain compression ratio. Our simulations demonstrate that entropic templating regulated by chain stiffness precisely reflects various entropic repulsion states that arise from typical conformation regimes of semiflexible chains. The physical mechanism of the chain stiffness effect is revealed by analyzing the entropic repulsion states of tethered chains and quantitatively estimating the resulting entropy penalties, which provides direct evidence that supports the key role of entropic transition in the entropic templating strategy, as suggested in experiments. Moreover, the model nanoparticle systems are found to evolve into binary nanoparticle superlattices by remixing at extremely high stiffness. The findings facilitate the wide application of the entropic templating strategy in creating interfacially reactive nanomaterials with ordered structures on the single-nanoparticle level as well as mechanomutable responses.

Flexible polyelectrolyte chain in a strong electrolyte solution: Insight into equilibrium properties and force-extension behavior from mesoscale simulation

Macromolecules with ionizable groups are ubiquitous in biological and synthetic systems. Due to the complex interaction between chain and electrostatic decorrelation lengths, both equilibrium properties and micro-mechanical response of dilute solutions of polyelectrolytes (PEs) are more complex than their neutral counterparts. In this work, the bead-rod micromechanical description of a chain is used to perform hi-fidelity Brownian dynamics simulation of dilute PE solutions to ascertain the self-similar equilibrium behavior of PE chains with various linear charge densities, scaling of the Kuhn step length (lE) with salt concentration cs and the force-extension behavior of the PE chain. In accord with earlier theoretical predictions, our results indicate that for a chain with n Kuhn segments, lE ∼ cs (-0.5) as linear charge density approaches 1/n. Moreover, the constant force ensemble simulation results accurately predict the initial non-linear force-extension region of PE chain recently measured via single chain experiments. Finally, inspired by Cohen's extraction of Warner's force law from the inverse Langevin force law, a novel numerical scheme is developed to extract a new elastic force law for real chains from our discrete set of force-extension data similar to Padè expansion, which accurately depicts the initial non-linear region where the total Kuhn length is less than the thermal screening length.

Modeling the stretching of wormlike chains in the presence of excluded volume

We propose an interpolation formula (the EV-WLC relation) for the force-extension behavior of wormlike chains in the presence of hard-core excluded volume interactions, analogous to the classic interpolation formula from Marko and Siggia for ideal wormlike chains. Using pruned-enriched Rosenbluth method (PERM) simulations of asymptotically long, discrete wormlike chains in an external force, we show that the error in the EV-WLC interpolation formula to describe discrete wormlike chains is systematically smaller than the error in the Marko-Siggia interpolation formula, except for the saturation region in which both formulas have the same limiting behavior. We anticipate that the EV-WLC interpolation formula will prove useful in the coarse-graining of wormlike chain models for dynamic simulations. Related results for the excess free energy due to excluded volume provide strong support for the physical basis of the Pincus regime.

Monte Carlo simulations of lattice models for single polymer systems

Single linear polymer chains in dilute solutions under good solvent conditions are studied by Monte Carlo simulations with the pruned-enriched Rosenbluth method up to the chain length N~O(10(4)). Based on the standard simple cubic lattice model (SCLM) with fixed bond length and the bond fluctuation model (BFM) with bond lengths in a range between 2 and √10, we investigate the conformations of polymer chains described by self-avoiding walks on the simple cubic lattice, and by random walks and non-reversible random walks in the absence of excluded volume interactions. In addition to flexible chains, we also extend our study to semiflexible chains for different stiffness controlled by a bending potential. The persistence lengths of chains extracted from the orientational correlations are estimated for all cases. We show that chains based on the BFM are more flexible than those based on the SCLM for a fixed bending energy. The microscopic differences between these two lattice models are discussed and the theoretical predictions of scaling laws given in the literature are checked and verified. Our simulations clarify that a different mapping ratio between the coarse-grained models and the atomistically realistic description of polymers is required in a coarse-graining approach due to the different crossovers to the asymptotic behavior.

Force-extension curves for broken-rod macromolecules: Dramatic effects of different probing methods for two and three rods

We study force-extension curves of a single semiflexible chain consisting of several rigid rods connected by flexible spacers. The atomic force microscopy and laser optical or magnetic tweezers apparatus stretching these rod-coil macromolecules are discussed. In addition, the stretching by external isotropic force is analyzed. The main attention is focused on computer simulation and analytical results. We demonstrate that the force-extension curves for rod-coil chains composed of two or three rods of equal length differ not only quantitatively but also qualitatively in different probe methods. These curves have an anomalous shape for a chain of two rods. End-to-end distributions of rod-coil chains are calculated by Monte Carlo method and compared with analytical equations. The influence of the spacer's length on the force-extension curves in different probe methods is analyzed. The results can be useful for interpreting experiments on the stretching of rod-coil block-copolymers.

Semiflexible polymer brushes and the brush-mushroom crossover

Semiflexible polymers end-grafted to a repulsive planar substrate under good solvent conditions are studied by scaling arguments, computer simulations, and self-consistent field theory. Varying the chain length N, persistence length lp, and grafting density σg, the chain linear dimensions and distribution functions of all monomers and of the free chain ends are studied. Particular attention is paid to the limit of very small σg, where the grafted chains behave as "mushrooms" no longer interacting with each other. Unlike a flexible mushroom, which has a self-similar structure from the size (a) of an effective monomer up to the mushroom height (h/a ∝ N(v), ν ≈ 3/5), a semiflexible mushroom (like a free semiflexible chain) exhibits three different scaling regimes, h/a ∝ N for contour length L = Na < lp, a Gaussian regime, h/a ∝ (Llp)(1/2)/a for lp ≪ L ≪ R* ∝ (lp(2)/a), and a regime controlled by excluded volume, h/a ∝ (lp/a)(1/5)N(ν). The semiflexible brush is predicted to scale as h/a ∝ (lpaσg)(1/3)N in the excluded volume regime, and h/a ∝ (lpa(3)σ(2))(1/4)N in the Gaussian regime. Since in the volume taken by a semiflexible mushroom excluded-volume interactions are much weaker in comparison to a flexible mushroom, there occurs an additional regime where semiflexible mushrooms overlap without significant chain stretching. Moreover, since the size of a semiflexible mushroom is much larger than the size of a flexible mushroom with the same N, the crossover from mushroom to brush behavior is predicted to take place at much smaller densities than for fully flexible chains. The numerical results, however, confirm the scaling predictions only qualitatively; for chain lengths that are relevant for experiments, often intermediate effective exponents are observed due to extended crossovers.

Unconventional ordering behavior of semi-flexible polymers in dense brushes under compression

Using a coarse-grained bead-spring model for semi-flexible macromolecules which form a polymer brush, the structure and dynamics of the polymers were investigated, varying the chain stiffness and the grafting density. The anchoring conditions for the grafted chains were chosen such that their first bonds were oriented along the normal to the substrate plane. The compression of such a semi-flexible brush by a planar piston was observed to be a two-stage process: for a small compression the chains were shown to contract by "buckling" deformation whereas for a larger compression the chains exhibited a collective (almost uniform) bending deformation. Thus, the stiff polymer brush underwent a 2nd order phase transition of collective bond reorientation. The pressure, required to keep the stiff brush at a given degree of compression, was thereby significantly smaller than for an otherwise identical brush made of entirely flexible polymer chains! While both the brush height and the chain linear dimensions in the z-direction perpendicular to the substrate increased monotonically with an increase in the chain stiffness, the lateral (xy) chain linear dimensions exhibited a maximum at an intermediate chain stiffness. Increasing the grafting density led to a strong decrease of these lateral dimensions which is compatible with an exponential decay. Also the recovery kinetics after removal of the compressing piston were studied, and were found to follow a power-law/exponential decay with time. A simple mean-field theoretical consideration, accounting for the buckling/bending behavior of semi-flexible polymer brushes under compression was suggested.

Driven translocation of a semi-flexible chain through a nanopore: a Brownian dynamics simulation study in two dimensions

We study translocation dynamics of a semi-flexible polymer chain through a nanoscopic pore in two dimensions using Langevin dynamics simulation in presence of an external bias F inside the pore. For chain length N and stiffness parameter κb considered in this paper, we observe that the mean first passage time <τ> increases as <τ(κb)>~<τ(κb=0)>lp(aN) , where κb and lp are the stiffness parameter and persistence length, respectively, and aN is a constant that has a weak N dependence. We monitor the time dependence of the last monomer xN(t) at the cis compartment and calculate the tension propagation time (TP) ttp directly from simulation data for <xN(t)> ~ t as alluded in recent nonequlibrium TP theory [T. Sakaue, Phys. Rev. E 76, 021803 (2007)] and its modifications to Brownian dynamics tension propagation theory [T. Ikonen, A. Bhattacharya, T. Ala-Nissila, and W. Sung, Phys. Rev. E 85, 051803 (2012); and J. Chem. Phys. 137, 085101 (2012)] originally developed to study translocation of a fully flexible chain. We also measure ttp from peak position of the waiting time distribution W(s) of the translocation coordinate s (i.e., the monomer inside the pore), and explicitly demonstrate the underlying TP picture along the chain backbone of a translocating chain to be valid for semi-flexible chains as well. From the simulation data, we determine the dependence of ttp on chain persistence length lp and show that the ratio ttp∕<τ> is independent of the bias F.

Is DNA a Good Model Polymer?

The details surrounding the cross-over from wormlike-specific to universal polymeric behavior has been the subject of debate and confusion even for the simple case of a dilute, unconfined wormlike chain. We have directly computed the polymer size, form factor, free energy and Kirkwood diffusivity for unconfined wormlike chains as a function of molecular weight, focusing on persistence lengths and effective widths that represent single-stranded and double-stranded DNA in a high ionic strength buffer. To do so, we use a chain-growth Monte Carlo algorithm, the Pruned-Enriched Rosenbluth Method (PERM), which allows us to estimate equilibrium and near-equilibrium dynamic properties of wormlike chains over an extremely large range of contour lengths. From our calculations, we find that very large DNA chains (≈ 1,000,000 base pairs depending on the choice of size metric) are required to reach flexible, swollen non-draining coils. Furthermore, our results indicate that the commonly used model polymer λ-DNA (48,500 base pairs) does not exhibit "ideal" scaling, but exists in the middle of the transition to long-chain behavior. We subsequently conclude that typical DNA used in experiments are too short to serve as an accurate model of long-chain, universal polymer behavior.

Stretching tethered polymer chains: density functional approach

We propose application of density functional theory to calculate the force acting on a selected segment of a tethered polymer chain that leads to stretching the chain. The density functional allows one to determine the effects due to the presence of other chains and solvent molecules. For high and moderate solvent densities the plot of the force versus the distance of the segment from the surface exhibits oscillatory behavior that has not been predicted by other approaches.

Scattering function of semiflexible polymer chains under good solvent conditions

Using the pruned-enriched Rosenbluth Monte Carlo algorithm, the scattering functions of semiflexible macromolecules in dilute solution under good solvent conditions are estimated both in d = 2 and d = 3 dimensions, considering also the effect of stretching forces. Using self-avoiding walks of up to N = 25,600 steps on the square and simple cubic lattices, variable chain stiffness is modeled by introducing an energy penalty ε(b) for chain bending; varying q(b) = exp (-ε(b)∕k(B)T) from q(b) = 1 (completely flexible chains) to q(b) = 0.005, the persistence length can be varied over two orders of magnitude. For unstretched semiflexible chains, we test the applicability of the Kratky-Porod worm-like chain model to describe the scattering function and discuss methods for extracting persistence length estimates from scattering. While in d = 2 the direct crossover from rod-like chains to self-avoiding walks invalidates the Kratky-Porod description, it holds in d = 3 for stiff chains if the number of Kuhn segments n(K) does not exceed a limiting value n(K)(*) (which depends on the persistence length). For stretched chains, the Pincus blob size enters as a further characteristic length scale. The anisotropy of the scattering is well described by the modified Debye function, if the actual observed chain extension <X> (end-to-end distance in the direction of the force) as well as the corresponding longitudinal and transverse linear dimensions <X(2)> - <X>(2), <R(g,⊥)(2)> are used.

Force spectroscopy of complex biopolymers with heterogeneous elasticity

Cellular biopolymers can exhibit significant compositional heterogeneities as a result of the non-uniform binding of associated proteins, the formation of microstructural defects during filament assembly, or the imperfect bundling of filaments into composite structures of variable diameter. These can lead to significant variations in the local mechanical properties of biopolymers along their length. Existing spectral analysis methods assume filament homogeneity and therefore report only a single average stiffness for the entire filament. However, understanding how local effects modulate biopolymer mechanics in a spatially resolved manner is essential to understanding how binding and bundling proteins regulate biopolymer stiffness and function in cellular contexts. Here, we present a new method to determine the spatially varying material properties of individual complex biopolymers from the observation of passive thermal fluctuations of the filament conformation. We develop new statistical mechanics-based approaches for heterogeneous filaments that estimate local bending elasticities as a function of the filament arc-length. We validate this methodology using simulated polymers with known stiffness distributions, and find excellent agreement between derived and expected values. We then determine the bending elasticity of microtubule filaments of variable composition generated by repeated rounds of tubulin polymerization using either GTP or GMPCPP, a nonhydrolyzable GTP analog. Again, we find excellent agreement between mechanical and compositional heterogeneities.

From flexible to stiff: systematic analysis of structural phases for single semiflexible polymers

Inspired by recent studies revealing unexpected pliability of semiflexible biomolecules like RNA and DNA, we systematically investigate the range of structural phases by means of a simple generic polymer model. Using a two-dimensional variant of Wang-Landau sampling to explore the conformational space in energy and stiffness within a single simulation, we identify the entire diversity of structures existing from the well-studied limit of flexible polymers to that of wormlike chains. We also discuss, in detail, the influence of finite-size effects in the formation of crystalline structures that are virtually inaccessible via conventional computational approaches.