Flexibility of short DNA helices with finite-length effect: From base pairs to tens of base pairs

Mechanical response of double-stranded DNA: Bend, twist, and overwind

We employed all-atom molecular dynamics simulations to explore the mechanical response of bending, twisting, and overwinding for double-stranded DNA (dsDNA). We analyzed the bending and twisting deformations, as well as their stiffnesses, using the tilt, roll, and twist modes under stretching force. Findings indicate that the roll and twist angles vary linearly with the stretching force but show opposite trends. The tilt, roll, and twist elastic moduli are considered constants, while the coupling between roll and twist modes slightly decreases under stretching force. The effect of the stretching force on the roll and twist modes, including both their deformations and elasticities, exhibits sequence-dependence, with symmetry around the base pair step. Furthermore, we examined the overwinding path and mechanism of dsDNA from the perspective of the stiffness matrix, based on the tilt, roll, and twist modes. The correlations among tilt, roll, and twist angles imply an alternative overwinding pathway via twist-roll coupling when dsDNA is stretched, wherein entropic contribution prevails.

Effect of ethanol on the elasticities of double-stranded RNA and DNA revealed by magnetic tweezers and simulations

The elasticities of double-stranded (ds) DNA and RNA, which are critical to their biological functions and applications in materials science, can be significantly modulated by solution conditions such as ions and temperature. However, there is still a lack of a comprehensive understanding of the role of solvents in the elasticities of dsRNA and dsDNA in a comparative way. In this work, we explored the effect of ethanol solvent on the elasticities of dsRNA and dsDNA by magnetic tweezers and all-atom molecular dynamics simulations. We found that the bending persistence lengths and contour lengths of dsRNA and dsDNA decrease monotonically with the increase in ethanol concentration. Furthermore, the addition of ethanol weakens the positive twist-stretch coupling of dsRNA, while promotes the negative twist-stretch coupling of dsDNA. Counter-intuitively, the lower dielectric environment of ethanol causes a significant re-distribution of counterions and enhanced ion neutralization, which overwhelms the enhanced repulsion along dsRNA/dsDNA, ultimately leading to the softening in bending for dsRNA and dsDNA. Moreover, for dsRNA, ethanol causes slight ion-clamping across the major groove, which weakens the major groove-mediated twist-stretch coupling, while for dsDNA, ethanol promotes the stretch-radius correlation due to enhanced ion binding and consequently enhances the helical radius-mediated twist-stretch coupling.

Effect of temperature on anisotropic bending elasticity of dsRNA: an all-atom molecular dynamics simulation

Employing all-atom molecular dynamics simulations, we examined the temperature-dependent behavior of bending elasticity in double-stranded RNA (dsRNA). Specifically, we focused on the bending persistence length and its constituent components, namely, the tilt and roll stiffness. Our results revealed a near-linear decrease in these stiffness components as a function of temperature, thereby highlighting the increased flexibility of dsRNA at elevated temperatures. Furthermore, our data revealed a significant anisotropy in dsRNA bending elasticity, which diminished with increasing temperature, attributable to marked disparities in tilt and roll stiffness components. We delineated the underlying biophysical mechanisms and corroborated our findings with extant literature. These observations offer salient implications for advancing our understanding of nucleic acid elasticity, and are pertinent to potential medical applications.

Influence of temperature on bend, twist and twist-bend coupling of dsDNA

The temperature-dependent bend and twist elasticities of dsDNA, as well as their couplings, were explored through all-atom molecular dynamics simulations. Three rotational parameters, tilt, roll, and twist, were employed to assess the bend and twist elasticities through their stiffness matrix. Our analysis indicates that the bend and twist stiffnesses decrease as the temperature rises, primarily owing to entropic influences stemming from thermodynamic fluctuations. Furthermore, the couplings between these rotational parameters also exhibit a decline with increasing temperature, although the roll-twist coupling displays greater strength than the tilt-roll and tilt-twist couplings, attributed to its more robust correction component. We elucidated the influence of temperature on bend and twist elasticities based on the comparisons between various models and existing data.

Visualizing the 4D Impact of Gold Nanoparticles on DNA

The genotoxicity of AuNPs has sparked a scientific debate, with one perspective attributing it to direct DNA damage and another to oxidative damage through reactive oxygen species (ROS) activation. This controversy poses challenges for the widespread use of AuNPs in biomedical applications. To address this debate, we employed four-dimensional atomic force microscopy (4DAFM) to examine the ability of AuNPs to damage DNA in vitro in the absence of ROS. To further examine whether the size and chemical coupling of these AuNPs are properties that control their toxicity, we exposed individual DNA molecules to three different types of AuNPs: small (average diameter = 10 nm), large (average diameter = 22 nm), and large conjugated (average diameter = 39 nm) AuNPs. We found that all types of AuNPs caused rapid (within minutes) and direct damage to the DNA molecules without the involvement of ROS. This research holds significant promise for advancing nanomedicines in diverse areas like viral therapy (including COVID-19), cancer treatment, and biosensor development for detecting DNA damage or mutations by resolving the ongoing debate regarding the genotoxicity mechanism. Moreover, it actively contributes to the continuous endeavors aimed at fully harnessing the capabilities of AuNPs across diverse biomedical fields, promising transformative healthcare solutions.

Exploration of the Character Representation of DNA Chiral Conformations and Deformations via a Curved Surface Discrete Frenet Frame

While undergoing structural deformation, DNA experiences changes in the interactions between its internal base pairs, presenting challenges to conventional elastic methods. To address this, we propose the Discrete Critical State (DCS) model in this paper. This model combines surface discrete frame theory with gauge theory and Landau phase transition theory to investigate DNA's structural deformation, phase transitions, and chirality. Notably, the DCS model considers both the internal interactions within DNA and formulates an overall equation using unified physical and geometric parameters. By employing the discrete frame, we derive the evolution of physical quantities along the helical axis of DNA, including geodesic curvature, geodesic torsion, and others. Our findings indicate that B-DNA has a significantly lower free energy density compared to Z-DNA, which is in agreement with experimental observations. This research reveals that the direction of base pairs is primarily governed by the geodesic curve within the helical plane, aligning closely with the orientation of the base pairs. Moreover, the geodesic curve has a profound influence on the arrangement of base pairs at the microscopic level and effectively regulates the configuration and geometry of DNA through macroscopic-level free energy considerations.

Twist-stretch relations in nucleic acids

Nucleic acids are highly deformable helical molecules constantly stretched, twisted and bent in their biological functioning. Single molecule experiments have shown that double stranded (ds)-RNA and standard ds-DNA have opposite twist-stretch patterns and stretching properties when overwound under a constant applied load. The key structural features of the A-form RNA and B-form DNA helices are here incorporated in a three-dimensional mesoscopic Hamiltonian model which accounts for the radial, bending and twisting fluctuations of the base pairs. Using path integral techniques which sum over the ensemble of the base pair fluctuations, I compute the average helical repeat of the molecules as a function of the load. The obtained twist-stretch relations and stretching properties, for short A- and B-helical fragments, are consistent with the opposite behaviors observed in kilo-base long molecules.

Asymmetric Periodic Boundary Conditions for All-Atom Molecular Dynamics and Coarse-Grained Simulations of Nucleic Acids

Periodic boundary conditions are commonly applied in molecular dynamics simulations in the microcanonical (NVE), canonical (NVT), and isothermal-isobaric (NpT) ensembles. In their simplest application, a biological system of interest is placed in the middle of a solvation box, which is chosen 'sufficiently large' to minimize any numerical artifacts associated with the periodic boundary conditions. This practical approach brings limitations to the size of biological systems that can be simulated. Here, we study simulations of effectively infinitely long nucleic acids, which are solvated in the directions perpendicular to the polymer chain, while periodic boundary conditions are also applied along the polymer chain. We study the effects of these asymmetric periodic boundary conditions (APBC) on the simulated results, including the mechanical properties of biopolymers and the properties of the surrounding solvent. To get some further insights into the advantages of using the APBC, a coarse-grained worm-like chain model is first studied, illustrating how the persistence length can be extracted from the local properties of the polymer chain, which are less affected by the APBC than some global averages. This is followed by all-atom molecular dynamics simulations of DNA in ionic solutions, where we use the APBC to investigate sequence-dependent properties of DNA molecules and properties of the surrounding solvent.

Mechanical codes of chemical-scale specificity in DNA motifs

In gene transcription, certain sequences of double-stranded (ds)DNA play a vital role in nucleosome positioning and expression initiation. That dsDNA is deformed to various extents in these processes leads us to ask: Could the genomic DNA also have sequence specificity in its chemical-scale mechanical properties? We approach this question using statistical machine learning to determine the rigidity between DNA chemical moieties. What emerges for the polyA, polyG, TpA, and CpG sequences studied here is a unique trigram that contains the quantitative mechanical strengths between bases and along the backbone. In a way, such a sequence-dependent trigram could be viewed as a DNA mechanical code. Interestingly, we discover a compensatory competition between the axial base-stacking interaction and the transverse base-pairing interaction, and such a reciprocal relationship constitutes the most discriminating feature of the mechanical code. Our results also provide chemical-scale understanding for experimental observables. For example, the long polyA persistence length is shown to have strong base stacking while its complement (polyAc) exhibits high backbone rigidity. The mechanical code concept enables a direct reading of the physical interactions encoded in the sequence which, with further development, is expected to shed new light on DNA allostery and DNA-binding drugs.

Computational Modeling of DNA 3D Structures: From Dynamics and Mechanics to Folding

DNA carries the genetic information required for the synthesis of RNA and proteins and plays an important role in many processes of biological development. Understanding the three-dimensional (3D) structures and dynamics of DNA is crucial for understanding their biological functions and guiding the development of novel materials. In this review, we discuss the recent advancements in computer methods for studying DNA 3D structures. This includes molecular dynamics simulations to analyze DNA dynamics, flexibility, and ion binding. We also explore various coarse-grained models used for DNA structure prediction or folding, along with fragment assembly methods for constructing DNA 3D structures. Furthermore, we also discuss the advantages and disadvantages of these methods and highlight their differences.

Temperature dependence of DNA elasticity: An all-atom molecular dynamics simulation study

We used all-atom molecular dynamics simulation to investigate the elastic properties of double-stranded DNA (dsDNA). We focused on the influences of temperature on the stretch, bend, and twist elasticities, as well as the twist-stretch coupling, of the dsDNA over a wide range of temperature. The results showed that the bending and twist persistence lengths, together with the stretch and twist moduli, decrease linearly with temperature. However, the twist-stretch coupling behaves in a positive correction and enhances as the temperature increases. The potential mechanisms of how temperature affects dsDNA elasticity and coupling were investigated by using the trajectories from atomistic simulation, in which thermal fluctuations in structural parameters were analyzed in detail. We analyzed the simulation results by comparing them with previous simulation and experimental data, which are in good agreement. The prediction about the temperature dependence of dsDNA elastic properties provides a deeper understanding of DNA elasticities in biological environments and potentially helps in the further development of DNA nanotechnology.

Mechanical Constraint Effect on DNA Persistence Length

Persistence length is a significant criterion to characterize the semi-flexibility of DNA molecules. The mechanical constraints applied on DNA chains in new single-molecule experiments play a complex role in measuring DNA persistence length; however, there is a difficulty in quantitatively characterizing the mechanical constraint effects due to their complex interactions with electrostatic repulsions and thermal fluctuations. In this work, the classical buckling theory of Euler beam and Manning's statistical theories of electrostatic force and thermal fluctuation force are combined for an isolated DNA fragment to formulate a quantitative model, which interprets the relationship between DNA persistence length and critical buckling length. Moreover, this relationship is further applied to identify the mechanical constraints in different DNA experiments by fitting the effective length factors of buckled fragments. Then, the mechanical constraint effects on DNA persistence lengths are explored. A good agreement among the results by theoretical models, previous experiments, and present molecular dynamics simulations demonstrates that the new superposition relationship including three constraint-dependent terms can effectively characterize changes in DNA persistence lengths with environmental conditions, and the strong constraint-environment coupling term dominates the significant changes of persistence lengths; via fitting effective length factors, the weakest mechanical constraints on DNAs in bulk experiments and stronger constraints on DNAs in single-molecule experiments are identified, respectively. Moreover, the consideration of DNA buckling provides a new perspective to examine the bendability of short-length DNA.

Structural dynamics of DNA strand break sensing by PARP-1 at a single-molecule level

Single-stranded breaks (SSBs) are the most frequent DNA lesions threatening genomic integrity. A highly kinked DNA structure in complex with human PARP-1 domains led to the proposal that SSB sensing in Eukaryotes relies on dynamics of both the broken DNA double helix and PARP-1's multi-domain organization. Here, we directly probe this process at the single-molecule level. Quantitative smFRET and structural ensemble calculations reveal how PARP-1's N-terminal zinc fingers convert DNA SSBs from a largely unperturbed conformation, via an intermediate state into the highly kinked DNA conformation. Our data suggest an induced fit mechanism via a multi-domain assembly cascade that drives SSB sensing and stimulates an interplay with the scaffold protein XRCC1 orchestrating subsequent DNA repair events. Interestingly, a clinically used PARP-1 inhibitor Niraparib shifts the equilibrium towards the unkinked DNA conformation, whereas the inhibitor EB47 stabilizes the kinked state.

Ab initio predictions for 3D structure and stability of single- and double-stranded DNAs in ion solutions

The three-dimensional (3D) structure and stability of DNA are essential to understand/control their biological functions and aid the development of novel materials. In this work, we present a coarse-grained (CG) model for DNA based on the RNA CG model proposed by us, to predict 3D structures and stability for both dsDNA and ssDNA from the sequence. Combined with a Monte Carlo simulated annealing algorithm and CG force fields involving the sequence-dependent base-pairing/stacking interactions and an implicit electrostatic potential, the present model successfully folds 20 dsDNAs (≤52nt) and 20 ssDNAs (≤74nt) into the corresponding native-like structures just from their sequences, with an overall mean RMSD of 3.4Å from the experimental structures. For DNAs with various lengths and sequences, the present model can make reliable predictions on stability, e.g., for 27 dsDNAs with/without bulge/internal loops and 24 ssDNAs including pseudoknot, the mean deviation of predicted melting temperatures from the corresponding experimental data is only ~2.0°C. Furthermore, the model also quantificationally predicts the effects of monovalent or divalent ions on the structure stability of ssDNAs/dsDNAs.

To determine 3D structures and quantify stability of single- (ss) and double-stranded (ds) DNAs is essential to unveil the mechanisms of their functions and to further guide the production and development of novel materials. Although many DNA models have been proposed to reproduce the basic structural, mechanical, or thermodynamic properties of dsDNAs based on the secondary structure information or preset constraints, there are very few models can be used to investigate the ssDNA folding or dsDNA assembly from the sequence. Furthermore, due to the polyanionic nature of DNAs, metal ions (e.g., Na⁺ and Mg²⁺) in solutions can play an essential role in DNA folding and dynamics. Nevertheless, ab initio predictions for DNA folding in ion solutions are still an unresolved problem. In this work, we developed a novel coarse-grained model to predict 3D structures and thermodynamic stabilities for both ssDNAs and dsDNAs in monovalent/divalent ion solutions from their sequences. As compared with the extensive experimental data and available existing models, we showed that the present model can successfully fold simple DNAs into their native-like structures, and can also accurately reproduce the effects of sequence and monovalent/divalent ions on structure stability for ssDNAs including pseudoknot and dsDNAs with/without bulge/internal loops.

Multivalent Cations Reverse the Twist-Stretch Coupling of RNA

When stretched, both DNA and RNA duplexes change their twist angles through twist-stretch coupling. The coupling is negative for DNA but positive for RNA, which is not yet completely understood. Here, our magnetic tweezers experiments show that the coupling of RNA reverses from positive to negative by multivalent cations. Combining with the previously reported tension-induced negative-to-positive coupling reversal of DNA, we propose a unified mechanism of the couplings of both RNA and DNA based on molecular dynamics simulations. Two deformation pathways are competing when stretched: shrinking the radius causes positive couplings but widening the major groove causes negative couplings. For RNA whose major groove is clamped by multivalent cations and canonical DNA, their radii shrink when stretched, thus exhibiting positive couplings. For elongated DNA whose radius already shrinks to the minimum and canonical RNA, their major grooves are widened when stretched, thus exhibiting negative couplings.

Nucleic Acid Identity, Structure, and Flexibility Affect the Electrochemical Signal of Tethered Redox Molecules upon Biopolymer Collapse

We demonstrate that cation condensation can induce the collapse of surface-bound nucleic acids and that the electrochemical signal from a tethered redox molecule (methylene blue) upon collapse reports on nucleic acid identity, structure, and flexibility. Furthermore, the correlation of the electrochemical signal and structure is consistent with theoretical considerations of nucleic acid collapse. Changes in solution dielectric permittivity or the concentration of trivalent cations cause the structure of nucleic acids to become more compact due to an increase in attractive electrostatic interactions between the charged biopolymer backbone and multivalent ions in the solution. Consequently, the compaction of nucleic acids results in a change in the dynamics and location of the terminally appended redox marker, which is reflected in the faradaic current measured using cyclic voltammetry. In comparison to ssDNA, nucleic acid duplexes (dsDNA, DNA/peptide nucleic acid, and dsRNA) require nucleic-acid-composition-specific solution conditions for the collapse to occur. Moreover, the magnitude of current increase observed after the collapse is different for each nucleic structure, and we find here that these changes are dictated by physical parameters of the nucleic acids including the axial charge spacing and the periodicity of the helix. The work here aims to provide quantitative and predicative measures of the effects of the nucleic acid structure on the electrochemical signal produced from distal-end appended redox markers. This architecture is commonly employed in functional nucleic acid sensors and a better understanding of structure-to-signal correlations will enable the rational design of sensitive sensing architectures.

Using capillary electrophoresis to characterize the hydrodynamic and electrostatic properties of DNA in solutions containing various monovalent cations

This review describes the results obtained by using free-solution capillary electrophoresis to probe the electrostatic and hydrodynamic properties of DNA in solutions containing various monovalent cations. In brief, we found that the mobilities of double-stranded DNAs (dsDNAs) increase with increasing molecular weight before leveling off and becoming constant at molecular weights ≥400 bp. The mobilities of single-stranded DNAs (ssDNAs) go through a maximum at ∼10-20 nucleotides before decreasing and becoming constant for oligomers larger than ∼30-50 bases. The mobilities of both ss- and dsDNAs increase linearly with the logarithm of increasing charge per unit length and decrease linearly with the logarithm of increasing ionic strength. Surprisingly, ss- and dsDNA mobilities level off and become nearly constant at ionic strengths ≥0.6 M. The thermal stabilities of dsDNAs decrease linearly with increasing solution viscosity. The diffusion coefficients of dsDNA are modulated by the diffusion coefficients of their counterions because of electrostatic DNA-cation coupling interactions. Finally, the anomalously slow mobilities observed for A-tract-containing DNAs can be attributed both to differences in shape and to the preferential localization of small cations in the A-tract minor groove. Since many of these results are mirrored in other polyion-counterion systems, free-solution electrophoresis can be viewed as a reporter of the electrostatics and hydrodynamics of highly charged polyions. New results describing the mobilities of dsDNA analogues of a microRNA-messenger RNA complex are also presented.

Sequence-dependent mechanics of collagen reflect its structural and functional organization

Extracellular matrix mechanics influence diverse cellular functions, yet surprisingly little is known about the mechanical properties of their constituent collagen proteins. In particular, network-forming collagen IV, an integral component of basement membranes, has been far less studied than fibril-forming collagens. A key feature of collagen IV is the presence of interruptions in the triple-helix-defining (Gly-X-Y) sequence along its collagenous domain. Here, we used atomic force microscopy to determine the impact of sequence heterogeneity on the local flexibility of collagen IV and of the fibril-forming collagen III. Our extracted flexibility profile of collagen IV reveals that it possesses highly heterogeneous mechanics, ranging from semiflexible regions as found for fibril-forming collagens to a lengthy region of high flexibility toward its N-terminus. A simple model in which flexibility is dictated only by the presence of interruptions fit the extracted profile reasonably well, providing insight into the alignment of chains and demonstrating that interruptions, particularly when coinciding in multiple chains, significantly enhance local flexibility. To a lesser extent, sequence variations within the triple helix lead to variable flexibility, as seen along the continuously triple-helical collagen III. We found this fibril-forming collagen to possess a high-flexibility region around its matrix-metalloprotease binding site, suggesting a unique mechanical fingerprint of this region that is key for matrix remodeling. Surprisingly, proline content did not correlate with local flexibility in either collagen type. We also found that physiologically relevant changes in pH and chloride concentration did not alter the flexibility of collagen IV, indicating such environmental changes are unlikely to control its compaction during secretion. Although extracellular chloride ions play a role in triggering collagen IV network formation, they do not appear to modulate the structure of its collagenous domain.

Ion-mediated interactions between like-charged polyelectrolytes with bending flexibility

Ion-mediated interactions between polyelectrolytes (PEs) are crucial to the properties of flexible biopolymers such as nucleic acids and proteins but the effect of PE flexibility on such interactions has not been explicitly addressed until now. In this work, the potentials of mean force (PMFs) between like-charged PEs with different bending flexibility have been investigated by Monte Carlo simulations and a cylindrical confinement around each PE was involved to model two PEs in an array. We found that in the absence of trivalent salt, the PMFs between like-charged PEs in an array are apparently repulsive while the bending flexibility can visibly decrease the repulsive PMFs. With the addition of high trivalent salt, the PMFs become significantly attractive whereas the attractive PMFs can be apparently weakened by the bending flexibility. Our analyses reveal that the effect of bending flexibility is attributed to the increased PE conformational space, which allows the PEs to fluctuate away to decrease the monovalent ion-mediated repulsion or to weaken the trivalent ion-mediated attraction through disrupting trivalent ion-bridging configuration. Additionally, our further calculations show that the effect of bending flexibility on the ion-mediated interactions is less apparent for PEs without cylindrical confinement.

Determinants of cyclization-decyclization kinetics of short DNA with sticky ends

Cyclization of DNA with sticky ends is commonly used to measure DNA bendability as a function of length and sequence, but how its kinetics depend on the rotational positioning of the sticky ends around the helical axis is less clear. Here, we measured cyclization (looping) and decyclization (unlooping) rates (k_loop and k_unloop) of DNA with sticky ends over three helical periods (100-130 bp) using single-molecule fluorescence resonance energy transfer (FRET). k_loop showed a nontrivial undulation as a function of DNA length whereas k_unloop showed a clear oscillation with a period close to the helical turn of DNA (∼10.5 bp). The oscillation of k_unloop was almost completely suppressed in the presence of gaps around the sticky ends. We explain these findings by modeling double-helical DNA as a twisted wormlike chain with a finite width, intrinsic curvature, and stacking interaction between the end base pairs. We also discuss technical issues for converting the FRET-based cyclization/decyclization rates to an equilibrium quantity known as the J factor that is widely used to characterize DNA bending mechanics.

Terahertz Wave Accelerates DNA Unwinding: A Molecular Dynamics Simulation Study

Unwinding the double helix of the DNA molecule is the basis of gene duplication and gene editing, and the acceleration of this unwinding process is crucial to the rapid detection of genetic information. Based on the unwinding of six-base-pair DNA duplexes, we demonstrate that a terahertz stimulus at a characteristic frequency (44.0 THz) can serve as an efficient, nonthermal, and long-range method to accelerate the unwinding process of DNA duplexes. The average speed of the unwinding process increased by 20 times at least, and its temperature was significantly reduced. The mechanism was revealed to be the resonance between the terahertz stimulus and the vibration of purine connected by the weak hydrogen bond and the consequent break in hydrogen bond connections between these base pairs. Our findings potentially provide a promising application of terahertz technology for the rapid detection of nucleic acids, biomedicine, and therapy.

Electrophoretic Mobility of DNA in Solutions of High Ionic Strength

The free-solution mobilities of small single-stranded DNA (ssDNA) and double-stranded DNA (dsDNA) have been measured by capillary electrophoresis in solutions containing 0.01-1.0 M sodium acetate. The mobility of dsDNA is greater than that of ssDNA at all ionic strengths because of the greater charge density of dsDNA. The mobilities of both ssDNA and dsDNA decrease with increasing ionic strength until approaching plateau values at ionic strengths greater than ∼0.6 M. Hence, ssDNA and dsDNA appear to interact in a similar manner with the ions in the background electrolyte. For dsDNA, the mobilities predicted by the Manning electrophoresis equation are reasonably close to the observed mobilities, using no adjustable parameters, if the average distance between phosphate residues (the b parameter) is taken to be 1.7 Å. For ssDNA, the predicted mobilities are close to the observed mobilities at ionic strengths ≤0.01 M if the b-value is taken to be 4.1 Å. The predicted and observed mobilities diverge strongly at higher ionic strengths unless the b-value is reduced significantly. The results suggest that ssDNA strands exist as an ensemble of relatively compact conformations at high ionic strengths, with b-values corresponding to the relatively short phosphate-phosphate distances through space.

Strongly Bent Double-Stranded DNA: Reconciling Theory and Experiment

The strong bending of polymers is poorly understood. We propose a general quantitative framework of polymer bending that includes both the weak and strong bending regimes on the same footing, based on a single general physical principle. As the bending deformation increases beyond a certain (polymer-specific) point, the change in the convexity properties of the effective bending energy of the polymer makes the harmonic deformation energetically unfavorable: in this strong bending regime the energy of the polymer varies linearly with the average bending angle as the system follows the convex hull of the deformation energy function. For double-stranded DNA, the effective bending deformation energy becomes non-convex for bends greater than ~ 2° per base-pair, equivalent to the curvature of a closed circular loop of ~ 160 base pairs. A simple equation is derived for the polymer loop energy that covers both the weak and strong bending regimes. The theory shows quantitative agreement with recent DNA cyclization experiments on short DNA fragments, while maintaining the expected agreement with experiment in the weak bending regime. Counter-intuitively, cyclization probability (j-factor) of very short DNA loops is predicted to increase with decreasing loop length; the j-factor reaches its minimum for loops of ≃ 45 base pairs. Atomistic simulations reveal that the attractive component of the short-range Lennard-Jones interaction between the backbone atoms can explain the underlying non-convexity of the DNA effective bending energy, leading to the linear bending regime. Applicability of the theory to protein-DNA complexes, including the nucleosome, is discussed.

The Contribution of Backbone Electrostatic Repulsion to DNA Mechanical Properties is Length-Scale-Dependent

The mechanics of DNA bending is crucially related to many vital biological processes. Recent experiments reported anomalous flexibility for DNA on short length scales, calling into doubt the validity of the harmonic worm-like chain (WLC) model in this region. In the present work, we systematically probed the bending dynamics of DNA at different length scales. In contrast to the remarkable deviation from the WLC description for DNA duplexes of less than three helical turns, our atomistic studies indicate that the neutral "null isomer" behaves in accord with the ideal elastic WLC and exhibits a uniform decay for the directional correlation of local bending. The backbone neutralization weakens the anisotropy in the effective bending preference and the helical periodicity of bend correlation that have previously been observed for normal DNA. The contribution of electrostatic repulsion to stretching cooperativity and the mechanical properties of DNA strands is length-scale-dependent: the phosphate neutralization increases the stiffness of DNA below two helical turns, but it is decreased for longer strands. We find that DNA rigidity is largely determined by base pair stacking, with electrostatic interactions contributing only around 10% of the total persistence length.

Structural Flexibility of DNA-RNA Hybrid Duplex: Stretching and Twist-Stretch Coupling

DNA-RNA hybrid (DRH) duplexes play essential roles during the replication of DNA and the reverse transcription of RNA viruses, and their flexibility is important for their biological functions. Recent experiments indicated that A-form RNA and B-form DNA have a strikingly different flexibility in stretching and twist-stretch coupling, and the structural flexibility of DRH duplex is of great interest, especially in stretching and twist-stretch coupling. In this work, we performed microsecond all-atom molecular dynamics simulations with new AMBER force fields to characterize the structural flexibility of DRH duplex in stretching and twist-stretch coupling. We have calculated all the helical parameters, stretch modulus, and twist-stretch coupling parameters for the DRH duplex. First, our analyses on structure suggest that the DRH duplex exhibits an intermediate conformation between A- and B-forms and closer to A-form, which can be attributed to the stronger rigidity of the RNA strand than the DNA strand. Second, our calculations show that the DRH duplex has the stretch modulus of 834 ± 34 pN and a very weak twist-stretch coupling. Our quantitative analyses indicate that, compared with DNA and RNA duplexes, the different flexibility of the DRH duplex in stretching and twist-stretch coupling is mainly attributed to its apparently different basepair inclination in the helical structure.

Competitive Binding of Mg2+ and Na+ Ions to Nucleic Acids: From Helices to Tertiary Structures

Nucleic acids generally reside in cellular aqueous solutions with mixed divalent/monovalent ions, and the competitive binding of divalent and monovalent ions is critical to the structures of nucleic acids because of their polyanionic nature. In this work, we first proposed a general and effective method for simulating a nucleic acid in mixed divalent/monovalent ion solutions with desired bulk ion concentrations via molecular dynamics (MD) simulations and investigated the competitive binding of Mg²⁺/Na⁺ ions to various nucleic acids by all-atom MD simulations. The extensive MD-based examinations show that single MD simulations conducted using the proposed method can yield desired bulk divalent/monovalent ion concentrations for various nucleic acids, including RNA tertiary structures. Our comprehensive analyses show that the global binding of Mg²⁺/Na⁺ to a nucleic acid is mainly dependent on its structure compactness, as well as Mg²⁺/Na⁺ concentrations, rather than the specific structure of the nucleic acid. Specifically, the relative global binding of Mg²⁺ over Na⁺ is stronger for a nucleic acid with higher effective surface charge density and higher relative Mg²⁺/Na⁺ concentrations. Furthermore, the local binding of Mg²⁺/Na⁺ to a phosphate of a nucleic acid mainly depends on the local phosphate density in addition to Mg²⁺/Na⁺ concentrations.

DNA Confined in a Nanodroplet: A Molecular Dynamics Study

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

Modeling Structure, Stability, and Flexibility of Double-Stranded RNAs in Salt Solutions

Double-stranded (ds) RNAs play essential roles in many processes of cell metabolism. The knowledge of three-dimensional (3D) structure, stability, and flexibility of dsRNAs in salt solutions is important for understanding their biological functions. In this work, we further developed our previously proposed coarse-grained model to predict 3D structure, stability, and flexibility for dsRNAs in monovalent and divalent ion solutions through involving an implicit structure-based electrostatic potential. The model can make reliable predictions for 3D structures of extensive dsRNAs with/without bulge/internal loops from their sequences, and the involvement of the structure-based electrostatic potential and corresponding ion condition can improve the predictions for 3D structures of dsRNAs in ion solutions. Furthermore, the model can make good predictions for thermal stability for extensive dsRNAs over the wide range of monovalent/divalent ion concentrations, and our analyses show that the thermally unfolding pathway of dsRNA is generally dependent on its length as well as its sequence. In addition, the model was employed to examine the salt-dependent flexibility of a dsRNA helix, and the calculated salt-dependent persistence lengths are in good accordance with experiments.

Ionic liquids make DNA rigid

Persistence length of double-stranded DNA (dsDNA) is known to decrease with an increase in ionic concentration of the solution. In contrast to this, here we show that the persistence length of dsDNA increases dramatically as a function of ionic liquid (IL) concentration. Using all atom explicit solvent molecular dynamics simulations and theoretical models, we present, for the first time, a systematic study to determine the mechanical properties of dsDNA in various hydrated ILs at different concentrations. We find that dsDNA in 50 wt % ILs have lower persistence length and stretch modulus in comparison to 80 wt % ILs. We further observe that both the persistence length and stretch modulus of dsDNA increase as we increase the concentration of ILs. The present trend of the stretch modulus and persistence length of dsDNA with IL concentration supports the predictions of the macroscopic elastic theory, in contrast to the behavior exhibited by dsDNA in monovalent salt. Our study further suggests the preferable ILs that can be used for maintaining DNA stability during long-term storage.

Radial distribution function of semiflexible oligomers with stretching flexibility

The radial distribution of the end-to-end distance R_ee is crucial for quantifying the global size and flexibility of a linear polymer. For semiflexible polymers, several analytical formulas have been derived for the radial distribution of R_ee ignoring the stretching flexibility. However, for semiflexible oligomers, such as DNA or RNA, the stretching flexibility can be rather pronounced and can significantly affect the radial distribution of R_ee. In this study, we obtained an extended formula that includes the stretch modulus to describe the distribution of R_ee for semiflexible oligomers on the basis of previous formulas for semiflexible polymers without stretching flexibility. The extended formula was validated by extensive Monte Carlo simulations over wide ranges of the stretch modulus and persistence length, as well as all-atom molecular dynamics simulations of short DNAs and RNAs. Additionally, our analyses showed that the effect of stretching flexibility on the distribution of R_ee becomes negligible for DNAs longer than ∼130 base pairs and RNAs longer than ∼240 base pairs.

Molecular rigidity and enthalpy-entropy compensation in DNA melting

Enthalpy-entropy compensation (EEC) is observed in diverse molecular binding processes of importance to living systems and manufacturing applications, but this widely occurring phenomenon is not sufficiently understood from a molecular physics standpoint. To gain insight into this fundamental problem, we focus on the melting of double-stranded DNA (dsDNA) since measurements exhibiting EEC are extensive for nucleic acid complexes and existing coarse-grained models of DNA allow us to explore the influence of changes in molecular parameters on the energetic parameters by using molecular dynamics simulations. Previous experimental and computational studies have indicated a correlation between EEC and changes in molecular rigidity in certain binding-unbinding processes, and, correspondingly, we estimate measures of DNA molecular rigidity under a wide range of conditions, along with resultant changes in the enthalpy and entropy of binding. In particular, we consider variations in dsDNA rigidity that arise from changes of intrinsic molecular rigidity such as varying the associative interaction strength between the DNA bases, the length of the DNA chains, and the bending stiffness of the individual DNA chains. We also consider extrinsic changes of molecular rigidity arising from the addition of polymer additives and geometrical confinement of DNA between parallel plates. All our computations confirm EEC and indicate that this phenomenon is indeed highly correlated with changes in molecular rigidity. However, two distinct patterns relating to how DNA rigidity influences the entropy of association emerge from our analysis. Increasing the intrinsic DNA rigidity increases the entropy of binding, but increases in molecular rigidity from external constraints decreases the entropy of binding. EEC arises in numerous synthetic and biological binding processes and we suggest that changes in molecular rigidity might provide a common origin of this ubiquitous phenomenon in the mutual binding and unbinding of complex molecules.

Kinetic theory for DNA melting with vibrational entropy

By treating DNA as a vibrating nonlinear lattice, an activated kinetic theory for DNA melting is developed to capture the breakage of the hydrogen bonds and subsequent softening of torsional and bending vibration modes. With a coarse-grained lattice model, we identify a key bending mode with GHz frequency that replaces the hydrogen vibration modes as the dominant out-of-phase phonon vibration at the transition state. By associating its bending modulus to a universal in-phase bending vibration modulus at equilibrium, we can hence estimate the entropic change in the out-of-phase vibration from near-equilibrium all-atom simulations. This and estimates of torsional and bending entropy changes lead to the first predictive and sequence-dependent theory with good quantitative agreement with experimental data for the activation energy of melting of short DNA molecules without intermediate hairpin structures.

Divalent Ion-Mediated DNA-DNA Interactions: A Comparative Study of Triplex and Duplex

Ion-mediated interaction between DNAs is essential for DNA condensation, and it is generally believed that monovalent and nonspecifically binding divalent cations cannot induce the aggregation of double-stranded (ds) DNAs. Interestingly, recent experiments found that alkaline earth metal ions such as Mg²⁺ can induce the aggregation of triple-stranded (ts) DNAs, although there is still a lack of deep understanding of the surprising findings at the microscopic level. In this work, we employed all-atom dynamic simulations to directly calculate the potentials of mean force (PMFs) between tsDNAs, between dsDNAs, and between tsDNA and dsDNA in Mg²⁺ solutions. Our calculations show that the PMF between tsDNAs is apparently attractive and becomes more strongly attractive at higher [Mg²⁺], although the PMF between dsDNAs cannot become apparently attractive even at high [Mg²⁺]. Our analyses show that Mg²⁺ internally binds into grooves and externally binds to phosphate groups for both tsDNA and dsDNA, whereas the external binding of Mg²⁺ is much stronger for tsDNA. Such stronger external binding of Mg²⁺ for tsDNA favors more apparent ion-bridging between helices than for dsDNA. Furthermore, our analyses illustrate that bridging ions, as a special part of external binding ions, are tightly and positively coupled to ion-mediated attraction between DNAs.

Understanding the Relative Flexibility of RNA and DNA Duplexes: Stretching and Twist-Stretch Coupling

The flexibility of double-stranded (ds) RNA and dsDNA is crucial for their biological functions. Recent experiments have shown that the flexibility of dsRNA and dsDNA can be distinctively different in the aspects of stretching and twist-stretch coupling. Although various studies have been performed to understand the flexibility of dsRNA and dsDNA, there is still a lack of deep understanding of the distinctive differences in the flexibility of dsRNA and dsDNA helices as pertains to their stretching and twist-stretch coupling. In this work, we have explored the relative flexibility in stretching and twist-stretch coupling between dsRNA and dsDNA by all-atom molecular dynamics simulations. The calculated stretch modulus and twist-stretch coupling are in good accordance with the existing experiments. Our analyses show that the differences in stretching and twist-stretch coupling between dsRNA and dsDNA helices are mainly attributed to their different (A- and B-form) helical structures. Stronger basepair inclination and slide in dsRNA is responsible for the apparently weaker stretching rigidity versus that of dsDNA, and the opposite twist-stretch coupling for dsRNA and dsDNA is also attributed to the stronger basepair inclination in dsRNA than in dsDNA. Our calculated macroscopic elastic parameters and microscopic analyses are tested and validated by different force fields for both dsRNA and dsDNA.

Twist-stretch profiles of DNA chains

Helical molecules change their twist number under the effect of a mechanical load. We study the twist-stretch relation for a set of short DNA molecules modeled by a mesoscopic Hamiltonian. Finite temperature path integral techniques are applied to generate a large ensemble of possible configurations for the base pairs of the sequence. The model also accounts for the bending and twisting fluctuations between adjacent base pairs along the molecules stack. Simulating a broad range of twisting conformation, we compute the helix structural parameters by averaging over the ensemble of base pairs configurations. The method selects, for any applied force, the average twist angle which minimizes the molecule's free energy. It is found that the chains generally over-twist under an applied stretching and the over-twisting is physically associated to the contraction of the average helix diameter, i.e. to the damping of the base pair fluctuations. Instead, assuming that the maximum amplitude of the bending fluctuations may decrease against the external load, the DNA molecule first over-twists for weak applied forces and then untwists above a characteristic force value. Our results are discussed in relation to available experimental information albeit for kilo-base long molecules.

Predicting 3D Structure, Flexibility, and Stability of RNA Hairpins in Monovalent and Divalent Ion Solutions

A full understanding of RNA-mediated biology would require the knowledge of three-dimensional (3D) structures, structural flexibility, and stability of RNAs. To predict RNA 3D structures and stability, we have previously proposed a three-bead coarse-grained predictive model with implicit salt/solvent potentials. In this study, we further develop the model by improving the implicit-salt electrostatic potential and including a sequence-dependent coaxial stacking potential to enable the model to simulate RNA 3D structure folding in divalent/monovalent ion solutions. The model presented here can predict 3D structures of RNA hairpins with bulges/internal loops (<77 nucleotides) from their sequences at the corresponding experimental ion conditions with an overall improved accuracy compared to the experimental data; the model also makes reliable predictions for the flexibility of RNA hairpins with bulge loops of different lengths at several divalent/monovalent ion conditions. In addition, the model successfully predicts the stability of RNA hairpins with various loops/stems in divalent/monovalent ion solutions.

Nucleic acid polymeric properties and electrostatics: Directly comparing theory and simulation with experiment

Nucleic acids are biopolymers that carry genetic information and are also involved in various gene regulation functions such as gene silencing and protein translation. Because of their negatively charged backbones, nucleic acids are polyelectrolytes. To adequately understand nucleic acid folding and function, we need to properly describe its i) polymer/polyelectrolyte properties and ii) associating ion atmosphere. While various theories and simulation models have been developed to describe nucleic acids and the ions around them, many of these theories/simulations have not been well evaluated due to complexities in comparison with experiment. In this review, I discuss some recent experiments that have been strategically designed for straightforward comparison with theories and simulation models. Such data serve as excellent benchmarks to identify limitations in prevailing theories and simulation parameters.

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.

Probing DNA Stiffness through Optical Fluctuation Analysis of Plasmon Rulers

The distance-dependent plasmon coupling between biopolymer tethered gold or silver nanoparticles forms the foundation for the so-called plasmon rulers. While conventional plasmon ruler applications focus on the detection of singular events in the far-field spectrum, we perform in this Letter a ratiometric analysis of the continuous spectral fluctuations arising from thermal interparticle separation variations in plasmon rulers confined to fluid lipid membranes. We characterized plasmon rulers with different DNA tethers and demonstrate the ability to detect and quantify differences in the plasmon ruler potential and tether stiffness. The influence of the nature of the tether (single-stranded versus double-stranded DNA) and the length of the tether is analyzed. The characterization of the continuous variation of the interparticle separation in individual plasmon rulers through optical fluctuation analysis provides additional information about the conformational flexibility of the tether molecule(s) located in the confinement of the deeply subdiffraction limit interparticle gap and enhances the versatility of plasmon rulers as a tool in Biophysics and Nanotechnology.

Multivalent ion-mediated nucleic acid helix-helix interactions: RNA versus DNA

Ion-mediated interaction is critical to the structure and stability of nucleic acids. Recent experiments suggest that the multivalent ion-induced aggregation of double-stranded (ds) RNAs and DNAs may strongly depend on the topological nature of helices, while there is still lack of an understanding on the relevant ion-mediated interactions at atomistic level. In this work, we have directly calculated the potentials of mean force (PMF) between two dsRNAs and between two dsDNAs in Co(NH3)6 (3+) (Co-Hex) solutions by the atomistic molecular dynamics simulations. Our calculations show that at low [Co-Hex], the PMFs between B-DNAs and between A-RNAs are both (strongly) repulsive. However, at high [Co-Hex], the PMF between B-DNAs is strongly attractive, while those between A-RNAs and between A-DNAs are still (weakly) repulsive. The microscopic analyses show that for A-form helices, Co-Hex would become 'internal binding' into the deep major groove and consequently cannot form the evident ion-bridge between adjacent helices, while for B-form helices without deep grooves, Co-Hex would exhibit 'external binding' to strongly bridge adjacent helices. In addition, our further calculations show that, the PMF between A-RNAs could become strongly attractive either at very high [Co-Hex] or when the bottom of deep major groove is fixed with a layer of water.