DNA Cyclization: Suppression or Enhancement by Electrostatic Repulsions?

Stability of DNA and RNA hairpins: a comparative study based on ox-DNA

Advances in single-molecule experiments on macromolecular crowding urgently need an efficient simulation method to resolve their discrepancies quantitatively. Ox-DNA model has been since reworked to treat the thermodynamics and mechanical properties of DNA/RNA hairpin at a stretching force. In hopping experiments, the critical forces of RNA hairpins at different temperatures are greater than those of DNA hairpins, in addition, the Gibbs free energy at a fixed temperature required to convert an RNA hairpin into a single-stranded molecule at zero force is obviously greater than that of DNA hairpin and gradually decreases by increasing the temperature. As far as force-ramping experiments are concerned, the first-rupture forces of RNA/DNA hairpins corresponding to the maximum probability density linearly pertain to the force-loading rate, with those of RNA hairpins being greater. The extended ox-DNA model could potentially identify the interaction between biologically inert polymer and RNA/DNA hairpins in crowded environments.

Charged elastic rings: deformation and dynamics

We report the counter-intuitive instability of charged elastic rings, and the persistence of sinusoidal deformations in the lowest-energy configurations by the combination of high-precision numerical simulations and analytical perturbation calculation. We also study the dynamical evolution of the charged ring under random disturbance, and reveal the modulation of the dominant frequencies by the electrostatic force. The purely mechanical analysis of the classical ring system presented in this work yields insights into the subtlety of long-range forces in the organization and dynamics of matter.

Mutual effects between single-stranded DNA conformation and Na+-Mg2+ ion competition in mixed salt solutions

The ion-dependence of single-stranded DNA (ssDNA) conformational changes has attracted growing attention because of its biological and technological importance. Although single-species ion effects have been extensively explored, it is challenging to study the ssDNA conformational properties under mixed monovalent/divalent ion conditions due to the complications of ssDNA flexibility and ion-ion competition. In this study, we apply Langevin dynamics simulations to investigate mixed Na⁺/Mg²⁺ ion-dependent ssDNA conformations. The ssDNA structure is described using a coarse-grained model, in which the phosphate, base, and sugar of each nucleotide are represented by three different beads. A novel improvement in our simulation model is that mixed-salt-related electrostatic interactions are computed via combining Manning counterion condensation (MCC) theory with the Monte Carlo tightly bound ion (MCTBI) model. Based on this MCC-MCTBI combination, we report new empirical functions to describe the ion-concentration-dependent and ssDNA conformation/structure-dependent electrostatic effects. The calculation results relating to the ion binding properties and the simulation results relating to the ssDNA conformational properties are validated against experimental results. In addition, our simulation results suggest a quantitative relationship between the ssDNA conformation and Na⁺-Mg²⁺ competition; this in turn reveals their mutual impact in the ion atmosphere.

Pseudomonas-tailed lytic phages: genome mechanical analysis and putative correlation with virion morphogenesis yield

Aim: To unveil a putative correlation between phage genome flexibility and virion morphogenesis yield. Materials & methods: A deeper analysis of the mechanical properties of three Pseudomonas aeruginosa lytic phage genomes was undertaken, together with full genome cyclizability calculations. Results & conclusion: A putative correlation was established among phage genome flexibility, eclipse timeframe and virion particle morphogenesis yield, with a more flexible phage genome leading to a higher burst size and a more rigid phage genome leading to lower burst sizes. The results obtained are highly relevant to understand the influence of the phage genome plasticity on the virion morphogenesis yield inside the infected bacterial host cells and assumes particular relevance in the actual context of bacterial resistance to antibiotics.

Mechanical Flexibility of DNA: A Quintessential Tool for DNA Nanotechnology

The mechanical properties of DNA have enabled it to be a structural and sensory element in many nanotechnology applications. While specific base-pairing interactions and secondary structure formation have been the most widely utilized mechanism in designing DNA nanodevices and biosensors, the intrinsic mechanical rigidity and flexibility are often overlooked. In this article, we will discuss the biochemical and biophysical origin of double-stranded DNA rigidity and how environmental and intrinsic factors such as salt, temperature, sequence, and small molecules influence it. We will then take a critical look at three areas of applications of DNA bending rigidity. First, we will discuss how DNA's bending rigidity has been utilized to create molecular springs that regulate the activities of biomolecules and cellular processes. Second, we will discuss how the nanomechanical response induced by DNA rigidity has been used to create conformational changes as sensors for molecular force, pH, metal ions, small molecules, and protein interactions. Lastly, we will discuss how DNA's rigidity enabled its application in creating DNA-based nanostructures from DNA origami to nanomachines.

Predicting Monovalent Ion Correlation Effects in Nucleic Acids

Ion correlation and fluctuation can play a potentially significant role in metal ion-nucleic acid interactions. Previous studies have focused on the effects for multivalent cations. However, the correlation and fluctuation effects can be important also for monovalent cations around the nucleic acid surface. Here, we report a model, gMCTBI, that can explicitly treat discrete distributions of both monovalent and multivalent cations and can account for the correlation and fluctuation effects for the cations in the solution. The gMCTBI model enables investigation of the global ion binding properties as well as the detailed discrete distributions of the bound ions. Accounting for the ion correlation effect for monovalent ions can lead to more accurate predictions, especially in a mixed monovalent and multivalent salt solution, for the number and location of the bound ions. Furthermore, although the monovalent ion-mediated correlation does not show a significant effect on the number of bound ions, the correlation may enhance the accumulation of monovalent ions near the nucleic acid surface and hence affect the ion distribution. The study further reveals novel ion correlation-induced effects in the competition between the different cations around nucleic acids.

The role of solvent quality and chain stiffness on the end-to-end contact kinetics of semiflexible polymers

Motivated by loop closure during protein folding and DNA packing, we systemically studied the effects of the solvent quality and chain stiffness on the thermodynamics and kinetics of the end-to-end contact formation for semiflexible polymer chains with reactive ends by Langevin dynamics simulations. In thermodynamics, a rich variety of products of the end-to-end contact have been discovered, such as loop, hairpin, toroid, and rodlike bundle, the populations of which are dependent on the solvent quality and chain stiffness. In kinetics, the overall pathways to form the end-to-end contact have been identified. The change of solvent quality and chain stiffness can tune the roughness of energy landscape and modulate the kinetic partitioning of the end-to-end contact formation pathways, leading to differing kinetic behaviors. In good or poor solvents, the first end-to-end contact rate k _c decreases with increasing the strength of bending stiffness k _θ monotonically. In very poor solvents, however, the dependence of the logarithm of the first end-to-end contact rate ln k _c on k _θ exhibits erratic behavior, which stems from more rugged energy landscape due to the polymer chain getting trapped into the intermediate state composed of the rodlike bundles with two ends in separation. For semiflexible chains, with increasing chain length N, the rate k _c increases initially and then decreases: in good solvents, the rate k _c exhibits a power-law relationship to chain length N with an exponent of ∼-1.50 in the region of long chains, which is in good agreement with the value derived from the experiment in the asymptotic limit of large N; and in poor solvents, the rate k _c exhibits a significantly stronger chain length dependence than those observed in good solvents in the region of long chains due to frustration to form the end-to-end contact along a specific path, especially the scaling exponent between the rate k _c and chain length N is ∼-3.62 for the case of polymer chains with k _θ = 4 at the solvent quality ε _ij = 1, in accord with the value obtained from the experiments.

Predicting RNA-Metal Ion Binding with Ion Dehydration Effects

Metal ions play essential roles in nucleic acids folding and stability. The interaction between metal ions and nucleic acids can be highly complicated because of the interplay between various effects such as ion correlation, fluctuation, and dehydration. These effects may be particularly important for multivalent ions such as Mg²⁺ ions. Previous efforts to model ion correlation and fluctuation effects led to the development of the Monte Carlo tightly bound ion model. Here, by incorporating ion hydration/dehydration effects into the Monte Carlo tightly bound ion model, we develop a, to our knowledge, new approach to predict ion binding. The new model enables predictions for not only the number of bound ions but also the three-dimensional spatial distribution of the bound ions. Furthermore, the new model reveals several intriguing features for the bound ions such as the mutual enhancement/inhibition in ion binding between the fully hydrated (diffuse) ions, the outer-shell dehydrated ions, and the inner-shell dehydrated ions and novel features for the monovalent-divalent ion interplay due to the hydration effect.

Salt-Induced Assembly Transformation of DNA-AuNP Conjugates Based on RCA Origami: From Linear Arrays to Nanorings

We developed a simple method to adjust the structural transformation of DNA-gold nanoparticles assemblies from linear arrays to nanorings by increasing salt concentrations. A DNA nanoladder constructing from RCA origami acted as templates to assemble periodic AuNPs arrays by a terminal thiol located on staple oligonucleotides. The linear AuNPs arrays could be transformed into nanorings only by changing the concentration of NaCl aqueous solution during the assembly process. It was proven that the electrostatic repulsion, being asymmetrically diminished by the high concentration of NaCl, caused the formation of nanoring architectures.

Extreme bendability of DNA double helix due to bending asymmetry

Experimental data of the DNA cyclization (J-factor) at short length scales exceed the theoretical expectation based on the wormlike chain (WLC) model by several orders of magnitude. Here, we propose that asymmetric bending rigidity of the double helix in the groove direction can be responsible for extreme bendability of DNA at short length scales and it also facilitates DNA loop formation at these lengths. To account for the bending asymmetry, we consider the asymmetric elastic rod (AER) model which has been introduced and parametrized in an earlier study [B. Eslami-Mossallam and M. R. Ejtehadi, Phys. Rev. E 80, 011919 (2009)]. Exploiting a coarse grained representation of the DNA molecule at base pair (bp) level and using the Monte Carlo simulation method in combination with the umbrella sampling technique, we calculate the loop formation probability of DNA in the AER model. We show that the DNA molecule has a larger J-factor compared to the WLC model which is in excellent agreement with recent experimental data.

Flexibility of short DNA helices under mechanical stretching

The flexibility of short DNA fragments is studied by a Hamiltonian model which treats the inter-strand and intra-strand forces at the level of the base pair.

Cylindrical confinement of semiflexible polymers

Equilibrium states of a closed semiflexible polymer binding to a cylinder are described. This may be either by confinement or by constriction. Closed completely bound states are labeled by two integers: the number of oscillations, n, and the number of times it winds the cylinder, p, the latter being a topological invariant. We examine the behavior of these states as the length of the loop is increased by evaluating the energy, the conserved axial torque, and the contact force. The ground state for a given p is the state with n=1; a short loop with p=1 is an elliptic deformation of a parallel circle; as its length increases it elongates along the cylinder axis with two hairpin ends. Excited states with n≥2 and p=1 possess n-fold axial symmetry. Short (long) loops possess energies ≈pE(0)(nE(0)), with E(0) the energy of a circular loop with same radius as the cylinder; in long loops the axial torque vanishes. Confined bound excited states are initially unstable; however, above a critical length each n-fold state becomes stable: The folded hairpin cannot be unfolded. The ground state for each p is also initially unstable with respect to deformations rotating the loop off the surface into the interior. A closed planar elastic curve aligned along the cylinder axis making contact with the cylinder on its two sides is identified as the ground state of a confined loop. Exterior bound states behave very differently, if free to unbind, as signaled by the reversal in the sign of the contact force. If p=1, all such states are unstable. If p≥2, however, a topological obstruction to complete unbinding exists. If the loop is short, the bound state with p=2 and n=1 provides a stable constriction of the cylinder, partially unbinding as the length is increased. This motif could be relevant to an understanding of the process of membrane fission mediated by dynamin rings.

Effects of the internal friction and the solvent quality on the dynamics of a polymer chain closure

Using 3D Langevin dynamics simulations, we investigate the effects of the internal friction and the solvent quality on the dynamics of a polymer chain closure. We show that the chain closure in good solvents is a purely diffusive process. By extrapolation to zero solvent viscosity, we find that the internal friction of a chain plays a non-ignorable role in the dynamics of the chain closure. When the solvent quality changes from good to poor, the mean closure time τc decreases by about 1 order of magnitude for the chain length 20 ≤ N ≤ 100. Furthermore, τc has a minimum as a function of the solvent quality. With increasing the chain length N, the minimum of τc occurs at a better solvent. Finally, the single exponential distributions of the closure time in poor solvents suggest that the negative excluded volume of segments does not alter the nearly Poisson statistical characteristics of the process of the chain closure.

Kinetics of polymer looping with macromolecular crowding: effects of volume fraction and crowder size

The looping of polymers such as DNA is a fundamental process in the molecular biology of living cells, whose interior is characterised by a high degree of molecular crowding. We here investigate in detail the looping dynamics of flexible polymer chains in the presence of different degrees of crowding. From the analysis of the looping-unlooping rates and the looping probabilities of the chain ends we show that the presence of small crowders typically slows down the chain dynamics but larger crowders may in fact facilitate the looping. We rationalise these non-trivial and often counterintuitive effects of the crowder size on the looping kinetics in terms of an effective solution viscosity and standard excluded volume. It is shown that for small crowders the effect of an increased viscosity dominates, while for big crowders we argue that confinement effects (caging) prevail. The tradeoff between both trends can thus result in the impediment or facilitation of polymer looping, depending on the crowder size. We also examine how the crowding volume fraction, chain length, and the attraction strength of the contact groups of the polymer chain affect the looping kinetics and hairpin formation dynamics. Our results are relevant for DNA looping in the absence and presence of protein mediation, DNA hairpin formation, RNA folding, and the folding of polypeptide chains under biologically relevant high-crowding conditions.

Modeling DNA condensation on freestanding cationic lipid membranes

Motivated by recent experimental observations of a rapid spontaneous DNA coil-globule transition on freestanding cationic lipid bilayers, we propose simple theoretical models for DNA condensation on cationic lipid membranes. First, for a single DNA rod, we examine the conditions of full wrapping of a cylindrical DNA-like semi-flexible polyelectrolyte by an oppositely charged membrane. Then, for two parallel DNA rods, we self-consistently analyze the shape and the extent of the membrane enveloping them, focusing on membrane elastic deformations and the membrane-DNA embracing angle, which enables us to compute the membrane-mediated DNA-DNA interactions. We examine the effects of the membrane composition and its charge density, which are the experimentally tunable parameters. We show that membrane-driven rod-rod attraction is more pronounced for higher charge densities and for smaller surface tensions of the membrane. Thus, we demonstrate that for a long DNA chain adhered to a cationic lipid membrane, such membrane-induced DNA-DNA attraction can trigger compaction of DNA.

Electrostatic effect of H1-histone protein binding on nucleosome repeat length

Within a simple biophysical model we describe the effect of electrostatic binding of H1 histone proteins on the nucleosome repeat length in chromatin. The length of wrapped DNA optimizes its binding energy to the histone core and the elastic energy penalty of DNA wrapping. The magnitude of the effect predicted from our model is in agreement with the systematic experimental data on the linear variation of nucleosome repeat lengths with H1/nucleosome ratio (Woodcock C L et al 2006 Chromos. Res. 14 17-25). We compare our model to the data for different cell types and organisms, with a widely varying ratio of bound H1 histones per nucleosome. We underline the importance of this non-specific histone-DNA charge-balance mechanism in regulating the positioning of nucleosomes and the degree of compaction of chromatin fibers in eukaryotic cells.

Understanding nucleic acid-ion interactions

Ions surround nucleic acids in what is referred to as an ion atmosphere. As a result, the folding and dynamics of RNA and DNA and their complexes with proteins and with each other cannot be understood without a reasonably sophisticated appreciation of these ions' electrostatic interactions. However, the underlying behavior of the ion atmosphere follows physical rules that are distinct from the rules of site binding that biochemists are most familiar and comfortable with. The main goal of this review is to familiarize nucleic acid experimentalists with the physical concepts that underlie nucleic acid-ion interactions. Throughout, we provide practical strategies for interpreting and analyzing nucleic acid experiments that avoid pitfalls from oversimplified or incorrect models. We briefly review the status of theories that predict or simulate nucleic acid-ion interactions and experiments that test these theories. Finally, we describe opportunities for going beyond phenomenological fits to a next-generation, truly predictive understanding of nucleic acid-ion interactions.

Mechanical properties of DNA-like polymers

The molecular structure of the DNA double helix has been known for 60 years, but we remain surprisingly ignorant of the balance of forces that determine its mechanical properties. The DNA double helix is among the stiffest of all biopolymers, but neither theory nor experiment has provided a coherent understanding of the relative roles of attractive base stacking forces and repulsive electrostatic forces creating this stiffness. To gain insight, we have created a family of double-helical DNA-like polymers where one of the four normal bases is replaced with various cationic, anionic or neutral analogs. We apply DNA ligase-catalyzed cyclization kinetics experiments to measure the bending and twisting flexibilities of these polymers under low salt conditions. Interestingly, we show that these modifications alter DNA bending stiffness by only 20%, but have much stronger (5-fold) effects on twist flexibility. We suggest that rather than modifying DNA stiffness through a mechanism easily interpretable as electrostatic, the more dominant effect of neutral and charged base modifications is their ability to drive transitions to helical conformations different from canonical B-form DNA.