Detailed scaling analysis of low-force polyelectrolyte elasticity

Noncooperative thermodynamics and kinetic models of ligand binding to polymers: Connecting McGhee-von Hippel model with the Tonks gas model

Ligand binding to polymers modifies the physical and chemical properties of the polymers, leading to physical, chemical, and biological implications. McGhee and von Hippel obtained the equilibrium coverage as a function of the ligand affinity, through the computation of the possible binding sites for the ligand. Here, we complete this theory deriving the kinetic model for the ligand-binding dynamics and the associated equilibrium chemical potential, which turns out to be of the Tonks gas model type. At low coverage, the Tonks chemical potential becomes the Fermi chemical potential and even the ideal gas chemical potential. We also discuss kinetic models associated with these chemical potentials. These results clarify the kinetic models of ligand binding, their relations with the chemical potentials, and their range of validity. Our results highlight the inaccuracy of ideal and simplified kinetic approaches for medium and high coverages.

Polyelectrolyte Conformation Controlled by a Trivalent-Rich Ion Jacket

The configuration of charged polymers is heavily dependent on interactions with surrounding salt ions, typically manifesting as a sensitivity to the bulk ionic strength. Here, we use single-molecule mechanical measurements to show that a charged polysaccharide, hyaluronic acid, shows a surprising regime of insensitivity to ionic strength in the presence of trivalent ions. Using simulations and theory, we propose that this is caused by the formation of a "jacket" of ions, tightly associated with the polymer, whose charge (and thus effect on configuration) is robust against changes in solution composition.

DNA Mechanics and Topology

We review the current understanding of the mechanics of DNA and DNA-protein complexes, from scales of base pairs up to whole chromosomes. Mechanics of the double helix as revealed by single-molecule experiments will be described, with an emphasis on the role of polymer statistical mechanics. We will then discuss how topological constraints- entanglement and supercoiling-impact physical and mechanical responses. Models for protein-DNA interactions, including effects on polymer properties of DNA of DNA-bending proteins will be described, relevant to behavior of protein-DNA complexes in vivo. We also discuss control of DNA entanglement topology by DNA-lengthwise-compaction machinery acting in concert with topoisomerases. Finally, the chapter will conclude with a discussion of relevance of several aspects of physical properties of DNA and chromatin to oncology.

Force-Induced Unravelling of DNA Origami

The mechanical properties of DNA nanostructures are of widespread interest as applications that exploit their stability under constant or intermittent external forces become increasingly common. We explore the force response of DNA origami in comprehensive detail by combining AFM single molecule force spectroscopy experiments with simulations using oxDNA, a coarse-grained model of DNA at the nucleotide level, to study the unravelling of an iconic origami system: the Rothemund tile. We contrast the force-induced melting of the tile with simulations of an origami 10-helix bundle. Finally, we simulate a recently proposed origami biosensor, whose function takes advantage of origami behavior under tension. We observe characteristic stick-slip unfolding dynamics in our force-extension curves for both the Rothemund tile and the helix bundle and reasonable agreement with experimentally observed rupture forces for these systems. Our results highlight the effect of design on force response: we observe regular, modular unfolding for the Rothemund tile that contrasts with strain-softening of the 10-helix bundle which leads to catastropic failure under monotonically increasing force. Further, unravelling occurs straightforwardly from the scaffold ends inward for the Rothemund tile, while the helix bundle unfolds more nonlinearly. The detailed visualization of the yielding events provided by simulation allows preferred pathways through the complex unfolding free-energy landscape to be mapped, as a key factor in determining relative barrier heights is the extensional release per base pair broken. We shed light on two important questions: how stable DNA nanostructures are under external forces and what design principles can be applied to enhance stability.

Sequence-Dependent Three Interaction Site Model for Single- and Double-Stranded DNA

We develop a robust coarse-grained model for single- and double-stranded DNA by representing each nucleotide by three interaction sites (TIS) located at the centers of mass of sugar, phosphate, and base. The resulting TIS model includes base-stacking, hydrogen bond, and electrostatic interactions as well as bond-stretching and bond angle potentials that account for the polymeric nature of DNA. The choices of force constants for stretching and the bending potentials were guided by a Boltzmann inversion procedure using a large representative set of DNA structures extracted from the Protein Data Bank. Some of the parameters in the stacking interactions were calculated using a learning procedure, which ensured that the experimentally measured melting temperatures of dimers are faithfully reproduced. Without any further adjustments, the calculations based on the TIS model reproduce the experimentally measured salt and sequence-dependence of the size of single-stranded DNA (ssDNA), as well as the persistence lengths of poly(dA) and poly(dT) chains. Interestingly, upon application of mechanical force, the extension of poly(dA) exhibits a plateau, which we trace to the formation of stacked helical domains. In contrast, the force-extension curve (FEC) of poly(dT) is entropic in origin and could be described by a standard polymer model. We also show that the persistence length of double-stranded DNA, formed from two complementary ssDNAs, is consistent with the prediction based on the worm-like chain. The persistence length, which decreases with increasing salt concentration, is in accord with the Odijk-Skolnick-Fixman theory intended for stiff polyelectrolyte chains near the rod limit. Our model predicts the melting temperatures of DNA hairpins with excellent accuracy, and we are able to recover the experimentally known sequence-specific trends. The range of applications, which did not require adjusting any parameter after the initial construction based solely on PDB structures and melting profiles of dimers, attests to the transferability and robustness of the TIS model for ssDNA and dsDNA.

Influence of ligand-receptor interactions on force-extension behavior within the freely jointed chain model

We study the influence of receptor-ligand interactions on the force response of single polymer chains theoretically. The extension of the chain is modeled in terms of freely jointed chain or elastic freely jointed chain (EFJC) models. The situation involving noninteracting bonds is solved exactly, while effects of interactions are treated within a mean-field approximation. The form with shorter bonds governs the low force situation, while the form with longer bonds is relevant in the high force regime. We further discuss the accuracy of approximate relations, which were used to describe the response of the EFJC model.

Mechanics, thermodynamics, and kinetics of ligand binding to biopolymers

Ligands binding to polymers regulate polymer functions by changing their physical and chemical properties. This ligand regulation plays a key role in many biological processes. We propose here a model to explain the mechanical, thermodynamic, and kinetic properties of the process of binding of small ligands to long biopolymers. These properties can now be measured at the single molecule level using force spectroscopy techniques. Our model performs an effective decomposition of the ligand-polymer system on its covered and uncovered regions, showing that the elastic properties of the ligand-polymer depend explicitly on the ligand coverage of the polymer (i.e., the fraction of the polymer covered by the ligand). The equilibrium coverage that minimizes the free energy of the ligand-polymer system is computed as a function of the applied force. We show how ligands tune the mechanical properties of a polymer, in particular its length and stiffness, in a force dependent manner. In addition, it is shown how ligand binding can be regulated applying mechanical tension on the polymer. Moreover, the binding kinetics study shows that, in the case where the ligand binds and organizes the polymer in different modes, the binding process can present transient shortening or lengthening of the polymer, caused by changes in the relative coverage by the different ligand modes. Our model will be useful to understand ligand-binding regulation of biological processes, such as the metabolism of nucleic acid. In particular, this model allows estimating the coverage fraction and the ligand mode characteristics from the force extension curves of a ligand-polymer system.

Spider Silk Peptide Is a Compact, Linear Nanospring Ideal for Intracellular Tension Sensing

Recent development and applications of calibrated, fluorescence resonance energy transfer (FRET)-based tension sensors have led to a new understanding of single molecule mechanotransduction in a number of biological systems. To expand the range of accessible forces, we systematically measured FRET versus force trajectories for 25, 40, and 50 amino acid peptide repeats derived from spider silk. Single molecule fluorescence-force spectroscopy showed that the peptides behaved as linear springs instead of the nonlinear behavior expected for a disordered polymer. Our data are consistent with a compact, rodlike structure that measures 0.26 nm per 5 amino acid repeat that can stretch by 500% while maintaining linearity, suggesting that the remarkable elasticity of spider silk proteins may in part derive from the properties of individual chains. We found the shortest peptide to have the widest range of force sensitivity: between 2 pN and 11 pN. Live cell imaging of the three tension sensor constructs inserted into vinculin showed similar force values around 2.4 pN. We also provide a lookup table for force versus intracellular FRET for all three constructs.

Sequence-dependent elasticity and electrostatics of single-stranded DNA: signatures of base-stacking

Base-stacking is a key factor in the energetics that determines nucleic acid structure. We measure the tensile response of single-stranded DNA as a function of sequence and monovalent salt concentration to examine the effects of base-stacking on the mechanical and thermodynamic properties of single-stranded DNA. By comparing the elastic response of highly stacked poly(dA) and that of a polypyrimidine sequence with minimal stacking, we find that base-stacking in poly(dA) significantly enhances the polymer's rigidity. The unstacking transition of poly(dA) at high force reveals that the intrinsic electrostatic tension on the molecule varies significantly more weakly on salt concentration than mean-field predictions. Further, we provide a model-independent estimate of the free energy difference between stacked poly(dA) and unstacked polypyrimidine, finding it to be ∼-0.25 kBT/base and nearly constant over three orders of magnitude in salt concentration.

The snakelike chain character of unstructured RNA

In the absence of base-pairing and tertiary structure, ribonucleic acid (RNA) assumes a random-walk conformation, modulated by the electrostatic self-repulsion of the charged, flexible backbone. This behavior is often modeled as a Kratky-Porod "wormlike chain" (WLC) with a Barrat-Joanny scale-dependent persistence length. In this study we report measurements of the end-to-end extension of poly(U) RNA under 0.1 to 10 pN applied force and observe two distinct elastic-response regimes: a low-force, power-law regime characteristic of a chain of swollen blobs on long length scales and a high-force, salt-valence-dependent regime consistent with ion-stabilized crumpling on short length scales. This short-scale structure is additionally supported by force- and salt-dependent quantification of the RNA ion atmosphere composition, which shows that ions are liberated under stretching; the number of ions liberated increases with increasing bulk salt concentration. Both this result and the observation of two elastic-response regimes directly contradict the WLC model, which predicts a single elastic regime across all forces and, when accounting for scale-dependent persistence length, the opposite trend in ion release with salt concentration. We conclude that RNA is better described as a "snakelike chain," characterized by smooth bending on long length scales and ion-stabilized crumpling on short length scales. In monovalent salt, these two regimes are separated by a characteristic length that scales with the Debye screening length, highlighting the determining importance of electrostatics in RNA conformation.

Elastic properties and secondary structure formation of single-stranded DNA at monovalent and divalent salt conditions

Single-stranded DNA (ssDNA) plays a major role in several biological processes. It is therefore of fundamental interest to understand how the elastic response and the formation of secondary structures are modulated by the interplay between base pairing and electrostatic interactions. Here we measure force-extension curves (FECs) of ssDNA molecules in optical tweezers set up over two orders of magnitude of monovalent and divalent salt conditions, and obtain its elastic parameters by fitting the FECs to semiflexible models of polymers. For both monovalent and divalent salts, we find that the electrostatic contribution to the persistence length is proportional to the Debye screening length, varying as the inverse of the square root of cation concentration. The intrinsic persistence length is equal to 0.7 nm for both types of salts, and the effectivity of divalent cations in screening electrostatic interactions appears to be 100-fold as compared with monovalent salt, in line with what has been recently reported for single-stranded RNA. Finally, we propose an analysis of the FECs using a model that accounts for the effective thickness of the filament at low salt condition and a simple phenomenological description that quantifies the formation of non-specific secondary structure at low forces.

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.

Fluctuations in the coil-stretch transition of flexible polymers in good solvents: a peak due to nonlinear force relation

Long flexible polymers undergo a coil to stretch transition (CST) in an elongational flow. Near the CST, a peak can be observed in the fluctuations of the size of a molecule (|R|). Solvent effects on the fluctuations are studied using Brownian dynamics simulations of a nonlinear spring force relation that can represent real molecules. Ignoring the influence of hydrodynamic interactions, a linear region in the spring force relation is known to cause the peak in |R| fluctuations. In contrast, we find that a peak in the fluctuations can be obtained even for the nonlinear spring force relation. We analyze the influence of hydrodynamic interactions on the fluctuations using a dumbbell model with a conformation-dependent drag coefficient.

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.

Quantifying screening ion excesses in single-molecule force-extension experiments

We derive a thermodynamic identity that allows one to infer the change in the number of screening ions that are associated with a charged macromolecule as the macromolecule is continuously stretched. Applying this identity to force-extension data on both single-stranded and double-stranded DNA, we find that the number of polymer-associated ions depends nontrivially on both the bulk salt concentration and the bare rigidity of the polymer, with single-stranded DNA exhibiting a relatively large decrease in ion excess upon stretching. We rationalize these observations using simple models for polyelectrolyte extension.

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

Semiflexible macromolecules in dilute solution under very good solvent conditions are modeled by self-avoiding walks on the simple cubic lattice (d = 3 dimensions) and square lattice (d = 2 dimensions), varying chain stiffness by an energy penalty ε(b) for chain bending. In the absence of excluded volume interactions, the persistence length l(p) of the polymers would then simply be l(p) = l(b)(2d - 2)(-1)q(b) (-1) with q(b) = exp(-ε(b)/k(B)T), the bond length l(b) being the lattice spacing, and k(B)T is the thermal energy. Using Monte Carlo simulations applying the pruned-enriched Rosenbluth method (PERM), both q(b) and the chain length N are varied over a wide range (0.005 ≤ q(b) ≤ 1, N ≤ 50,000), and also a stretching force f is applied to one chain end (fixing the other end at the origin). In the absence of this force, in d = 2 a single crossover from rod-like behavior (for contour lengths less than l(p)) to swollen coils occurs, invalidating the Kratky-Porod model, while in d = 3 a double crossover occurs, from rods to Gaussian coils (as implied by the Kratky-Porod model) and then to coils that are swollen due to the excluded volume interaction. If the stretching force is applied, excluded volume interactions matter for the force versus extension relation irrespective of chain stiffness in d = 2, while theories based on the Kratky-Porod model are found to work in d = 3 for stiff chains in an intermediate regime of chain extensions. While for q(b) ≪ 1 in this model a persistence length can be estimated from the initial decay of bond-orientational correlations, it is argued that this is not possible for more complex wormlike chains (e.g., bottle-brush polymers). Consequences for the proper interpretation of experiments are briefly discussed.

Single-molecule elasticity measurements of the onset of excluded volume in poly(ethylene glycol)

A polymer must reach a certain size to exhibit significant excluded-volume interactions and adopt a swollen random-walk configuration. We show that single-molecule measurements can sense the onset of swelling by modulating the effective chain size with force: as the force is reduced from a large value, the polymer is first highly aligned, then a Gaussian coil, then finally a swollen chain, with each regime exhibiting a distinct elasticity. We use this approach to quantify the structural parameters of poly(ethylene glycol) and show that they vary in the expected manner with changes in solvent.

Quantitative High-Resolution Sensing of DNA Hybridization Using Magnetic Tweezers with Evanescent Illumination

We applied the combined approach of evanescent nanometry and force spectroscopy using magnetic tweezers to quantify the degree of hybridization of a single synthetic single-stranded DNA oligomer to a resolution approaching a single-base. In this setup, the 200 nucleotide long DNA was covalently attached to the surface of an optically transparent solid support at one end and to the surface of a superparamagnetic fluorescent microsphere (force probe) at the other end. The force was applied to the probes using an electromagnet. The end-to-end molecular distance (i.e. out-of-image-plane position of the force probe) was determined from the intensity of the probe fluorescent image observed with total-internal reflectance microscopy. An equation of state for single stranded DNA molecules under tension (extensible freely jointed chain) was used to derive the penetration depth of the evanescent field and to calibrate the magnetic properties of the force probes. The parameters of the magnetic response of the force probes obtained from the equation of state remained constant when changing the penetration depth, indicating a robust calibration procedure. The results of such a calibration were also confirmed using independently measured probe-surface distances for probes mounted onto cantilevers of an atomic force microscope. Upon hybridization of the complementary 50 nucleotide-long oligomer to the surface-bound 200-mer, the changes in the force-distance curves were consistent with the quantitative conversion of 25% of the original single-stranded DNA to its double-stranded form, which was modeled as an elastic rod. The method presented here for quantifying the hybridization state of the single DNA molecules has potential for determining the degree of hybridization of individual molecules in a single molecule array with high accuracy.

Direct observation of single flexible polymers using single stranded DNA()

Over the last 15 years, double stranded DNA (dsDNA) has been used as a model polymeric system for nearly all single polymer dynamics studies. However, dsDNA is a semiflexible polymer with markedly different molecular properties compared to flexible chains, including synthetic organic polymers. In this work, we report a new system for single polymer studies of flexible chains based on single stranded DNA (ssDNA). We developed a method to synthesize ssDNA for fluorescence microscopy based on rolling circle replication, which generates long strands (>65 kb) of ssDNA containing "designer" sequences, thereby preventing intramolecular base pair interactions. Polymers are synthesized to contain amine-modified bases randomly distributed along the backbone, which enables uniform labelling of polymer chains with a fluorescent dye to facilitate fluorescence microscopy and imaging. Using this approach, we synthesized ssDNA chains with long contour lengths (>30 μm) and relatively low dye loading ratios (~1 dye per 100 bases). In addition, we used epifluorescence microscopy to image single ssDNA polymer molecules stretching in flow in a microfluidic device. Overall, we anticipate that ssDNA will serve as a useful model system to probe the dynamics of polymeric materials at the molecular level.

Model systems for single molecule polymer dynamics

Double stranded DNA (dsDNA) has long served as a model system for single molecule polymer dynamics. However, dsDNA is a semiflexible polymer, and the structural rigidity of the DNA double helix gives rise to local molecular properties and chain dynamics that differ from flexible chains, including synthetic organic polymers. Recently, we developed single stranded DNA (ssDNA) as a new model system for single molecule studies of flexible polymer chains. In this work, we discuss model polymer systems in the context of "ideal" and "real" chain behavior considering thermal blobs, tension blobs, hydrodynamic drag and force-extension relations. In addition, we present monomer aspect ratio as a key parameter describing chain conformation and dynamics, and we derive dynamical scaling relations in terms of this molecular-level parameter. We show that asymmetric Kuhn segments can suppress monomer-monomer interactions, thereby altering global chain dynamics. Finally, we discuss ssDNA in the context of a new model system for single molecule polymer dynamics. Overall, we anticipate that future single polymer studies of flexible chains will reveal new insight into the dynamic behavior of "real" polymers, which will highlight the importance of molecular individualism and the prevalence of non-linear phenomena.