Salt contribution to the flexibility of single-stranded nucleic acid offinite length

Conformational Transition of Semiflexible Ring Polyelectrolyte in Tetravalent Salt Solutions: A Simple Numerical Modeling without the Effect of Twisting

In this work, the conformational behaviors of ring polyelectrolyte in tetravalent salt solutions are discussed in detail through molecular dynamics simulation. For simplification, here we have neglected the effect of the twisting interaction, although it has been well known that both bending and twisting interactions play a deterministic in the steric conformation of a semiflexible ring polymer. The salt concentration CS and the bending energy b take a decisive role in the conformation of the ring polyelectrolyte (PE). Throughout our calculations, the b varies from b = 0 (freely joint chain) to b = 120. The salt concentration CS changes in the range of 3.56 × 10-4 M ≤ CS ≤ 2.49 × 10-1 M. Upon the addition of salt, ring PE contracts at first, subsequently re-expands. More abundant conformations are observed for a semiflexible ring PE. For b = 10, the conformation of semiflexible ring PE shifts from the loop to two-racquet-head spindle, then it condenses into toroid, finally arranges into coil with the increase of CS. As b increases further, four phase transitions are observed. The latter two phase transitions are different. The semiflexible ring PE experiences transformation from toroid to two racquet head spindle, finally to loop in the latter two phase transitions. Its conformation is determined by the competition among the bending energy, cation-bridge, and entropy. Combined, our findings indicate that the conformations of semiflexible ring PE can be controlled by changing the salt concentration and chain stiffness.

RNA 3D Structure Prediction: Progress and Perspective

Ribonucleic acid (RNA) molecules play vital roles in numerous important biological functions such as catalysis and gene regulation. The functions of RNAs are strongly coupled to their structures or proper structure changes, and RNA structure prediction has been paid much attention in the last two decades. Some computational models have been developed to predict RNA three-dimensional (3D) structures in silico, and these models are generally composed of predicting RNA 3D structure ensemble, evaluating near-native RNAs from the structure ensemble, and refining the identified RNAs. In this review, we will make a comprehensive overview of the recent advances in RNA 3D structure modeling, including structure ensemble prediction, evaluation, and refinement. Finally, we will emphasize some insights and perspectives in modeling RNA 3D structures.

Machine Learning Color Feature Analysis of a High Throughput Nanoparticle Conjugate Sensing Assay

Plasmonic nanoparticles are finding applications within the single molecule sensing field in a "dimer" format, where interaction of the target with hairpin DNA causes a decrease in the interparticle distance, leading to a localized surface plasmon resonance shift. While this shift may be detected using spectroscopy, achieving statistical relevance requires the measurement of thousands of nanoparticle dimers and the timescales required for spectroscopic analysis are incompatible with point-of-care devices. However, using dark-field imaging of the dimer structures, simultaneous digital analysis of the plasmonic resonance shift after target interaction of thousands of dimer structures may be achieved in minutes. The main challenge of this digital analysis on the single-molecule scale was the occurrence of false signals caused by non-specifically bound clusters of nanoparticles. This effect may be reduced by digitally separating dimers from other nanoconjugate types. Variation in image intensity was observed to have a discernible impact on the color analysis of the nanoconjugate constructs and thus the accuracy of the digital separation. Color spaces wherein intensity may be uncoupled from the color information (hue, saturation, and value (HSV) and luminance, a* vector, and b* vector (LAB) were contrasted to a color space which cannot uncouple intensity (RGB) to train a classifier algorithm. Each classifier algorithm was validated to determine which color space produced the most accurate digital separation of the nanoconjugate types. The LAB-based learning classifier demonstrated the highest accuracy for digitally separating nanoparticles. Using this classifier, nanoparticle conjugates were monitored for their plasmonic color shift after interaction with a synthetic RNA target, resulting in a platform with a highly accurate yes/no response with a true positive rate of 88% and a true negative rate of 100%. The sensor response of tested single stranded RNA (ssRNA) samples was well above control responses for target concentrations in the range of 10 aM-1 pM.

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.

Modulation of DNA conformation in electrolytic nanodroplets

The behavior of deoxyribonucleic acid (DNA) molecules in confinement is of profound importance in various bioengineering and medical applications. In the present study, all-atom molecular dynamics simulation is utilized to investigate the transition of the double-strand DNA (dsDNA) conformation in the electrolytic nanodroplet. Three typical conformations, i.e., C-shaped, folded S-shaped, and double C-shaped, are observed for different droplet sizes and ionic concentrations. To reveal the physics underlying this phenomenon, the characteristics of the dsDNA molecules, such as the overcharging intensity, the end-to-end distance, the radius of gyration, etc. are analyzed in detail based on the numerical results. It is found that the transition can be ascribed to the buckling of the polymer molecules under the compression due to the confinement of the nanodroplet, and it can be modulated by the ionic concentration in the electrolyte. Generally, nanoscale confinement dominates dsDNA behavior over the electrostatic effects in smaller nanodroplets, while the latter becomes more important for larger nanodroplets. This competition results in the persistence length increasing with the nanodroplet radii. Based on these discussions, a non-dimensional elasto-capillary number μ is proposed to classify the dsDNA conformations into three regions.

Effect of mono- and multi-valent ionic environments on the in-lattice nanoparticle-grafted single-stranded DNA

Polyelectrolyte (PE) chains respond in a complex manner to multivalent salt environments, and this behavior depends on pH, temperature, and the presence of specific counter ions. Although much work has been done to understand the behaviour of free PE chains, it is important to reveal their behaviour on a nanoparticle's surface, where surface constraints, particle geometry, and multi-chain environment can affect their behaviour and contribute to particles' assembly states. Our work investigates, using in situ small-angle X-ray scattering (SAXS), the morphology of PE (single-stranded DNA) chains grafted onto the surface of spherical gold nanoparticles assembled in a lattice in the presence of monovalent, divalent and trivalent salts. For divalent salts, the DNA brush length was found to decrease at a faster rate with salt concentration than in the monovalent salt environment, while trivalent salts led to chain collapse. Using a power law analysis and the modified Daoud-Cotton model, we have obtained insight into the mechanism of a nanoparticle-grafted chain's response to ionic environments. Our analysis suggests that the decrease in brush length is due to the conventional electrostatic screening for monovalent systems, whereas for divalent systems both electrostatic screening and divalent ion bridging must be considered.

Resistive amplitude fingerprints during translocation of linear molecules through charged solid-state nanopores

We report the first analytical theory on the amplitude of resistive signals during molecular translocation through charged solid-state nanopores with variable cross-sectional area and piecewise-constant surface charge densities. By providing closed-form explicit algebraic expressions for the concentration profiles inside charged nanopores, this theory allows the prediction of baseline and translocation resistive signals without the need for numerical simulation of the electrokinetic phenomena. A transversely homogenized theory and an asymptotic expansion for weakly charged pores capture DC or quasi-static rectification due to field-induced intrapore concentration polarization (as a result of pore charge inhomogeneity or a translocating molecule). This theory, validated by simulations and experiments, is then used to explain why the amplitude of a single stranded DNA molecule can be twice as high as the amplitude of its double stranded counterpart. It also suggests designs for intrapore concentration polarization and volume exclusion effects that can produce biphasic and other amplitude fingerprints for high-throughput and yet discriminating molecular identification.

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.

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.

Predicting 3D structure and stability of RNA pseudoknots in monovalent and divalent ion solutions

RNA pseudoknots are a kind of minimal RNA tertiary structural motifs, and their three-dimensional (3D) structures and stability play essential roles in a variety of biological functions. Therefore, to predict 3D structures and stability of RNA pseudoknots is essential for understanding their functions. In the work, we employed our previously developed coarse-grained model with implicit salt to make extensive predictions and comprehensive analyses on the 3D structures and stability for RNA pseudoknots in monovalent/divalent ion solutions. The comparisons with available experimental data show that our model can successfully predict the 3D structures of RNA pseudoknots from their sequences, and can also make reliable predictions for the stability of RNA pseudoknots with different lengths and sequences over a wide range of monovalent/divalent ion concentrations. Furthermore, we made comprehensive analyses on the unfolding pathway for various RNA pseudoknots in ion solutions. Our analyses for extensive pseudokonts and the wide range of monovalent/divalent ion concentrations verify that the unfolding pathway of RNA pseudoknots is mainly dependent on the relative stability of unfolded intermediate states, and show that the unfolding pathway of RNA pseudoknots can be significantly modulated by their sequences and solution ion conditions.

Potential of mean force between oppositely charged nanoparticles: A comprehensive comparison between Poisson-Boltzmann theory and Monte Carlo simulations

Ion-mediated interactions between like-charged polyelectrolytes have been paid much attention, and the Poisson-Boltzmann (PB) theory has been shown to fail in qualitatively predicting multivalent ion-mediated like-charge attraction. However, inadequate attention has been paid to the ion-mediated interactions between oppositely charged polyelectrolytes. In this work, the potentials of mean force (PMF) between oppositely charged nanoparticles in 1:1 and 2:2 salt solutions were investigated by Monte Carlo simulations and the PB theory. Our calculations show that the PMFs between oppositely charged nanoparticles are generally attractive in 1:1 and 2:2 salt solutions and that such attractive PMFs become weaker at higher 1:1 or 2:2 salt concentrations. The comprehensive comparisons show that the PB theory can quantitatively predict the PMFs between oppositely charged nanoparticles in 1:1 salt solutions, except for the slight deviation at very high 1:1 salt concentration. However, for 2:2 salt solutions, the PB theory generally overestimates the attractive PMF between oppositely charged nanoparticles, and this overestimation becomes more pronounced for nanoparticles with higher charge density and for higher 2:2 salt concentration. Our microscopic analyses suggest that the overestimation of the PB theory on the attractive PMFs for 2:2 salt solutions is attributed to the underestimation of divalent ions bound to nanoparticles.

DNA-Assisted Dispersion of Carbon Nanotubes and Comparison with Other Dispersing Agents

Separation and sorting of pristine carbon nanotubes (CNTs) from bundle geometry is a very challenging task due to the insoluble and nondispersive nature of CNTs in aqueous medium. Recently, many studies have been performed to address this problem using various organic and inorganic solutions, surfactant molecules, and biomolecules as dispersing agents. Recent experimental studies have reported the DNA to be highly efficient in dispersing CNTs from bundle geometry. However, there is no microscopic study and also quantitative estimation of the dispersion efficiency of the DNA. Using all-atom molecular dynamics simulation, we study the structure and stability of single-stranded DNA (ssDNA)-single-walled carbon nanotube (SWNT) (6,5) complex. To quantify the dispersion efficiency of various DNA sequences, we perform potential of mean forces (PMF) calculation between two bare SWNTs as well ssDNA-wrapped CNTs for different base sequences. From the PMF calculation, we find the PMF between two bare (6,5) SWNTs to be approximately -29 kcal/mol. For the ssDNA-wrapped SWNTs, the PMF reduces significantly and becomes repulsive. In the presence of ssDNA of different polynucleotide bases (A, T, G, and C), we present a microscopic picture of the ssDNA-SWNT (6,5) complex and also a quantitative estimate of the interaction strength between nanotubes from PMF calculation. From PMF, we show the sequence of dispersion efficiency for four different nucleic bases to be T > A > C > G. We have also presented a comparison of the dispersion efficiencies of ssDNA, flavin mononucleotide surfactant, and poly(amidoamine) (PAMAM) dendrimer by comparing their respective PMF values.

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.

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.

Scaling analysis of random walks with persistence lengths: Application to self-avoiding walks

We develop an approach for performing scaling analysis of N-step random walks (RWs). The mean square end-to-end distance, 〈R[over ⃗]_{N}^{2}〉, is written in terms of inner persistence lengths (IPLs), which we define by the ensemble averages of dot products between the walker's position and displacement vectors, at the jth step. For RW models statistically invariant under orthogonal transformations, we analytically introduce a relation between 〈R[over ⃗]_{N}^{2}〉 and the persistence length, λ_{N}, which is defined as the mean end-to-end vector projection in the first step direction. For self-avoiding walks (SAWs) on 2D and 3D lattices we introduce a series expansion for λ_{N}, and by Monte Carlo simulations we find that λ_{∞} is equal to a constant; the scaling corrections for λ_{N} can be second- and higher-order corrections to scaling for 〈R[over ⃗]_{N}^{2}〉. Building SAWs with typically 100 steps, we estimate the exponents ν_{0} and Δ_{1} from the IPL behavior as function of j. The obtained results are in excellent agreement with those in the literature. This shows that only an ensemble of paths with the same length is sufficient for determining the scaling behavior of 〈R[over ⃗]_{N}^{2}〉, being that the whole information needed is contained in the inner part of the paths.

Exploring ion induced folding of a single-stranded DNA oligomer from molecular simulation studies

One crucial issue in DNA hydration is the effect of salts on its conformational features. This has relevance in biology as cations present in the cellular environment shield the negative charges on the DNA backbone, thereby reducing the repulsive force between them. By screening the negative charges along the backbone, cations stabilize the folded structure of DNA. To study the effect of the added salt on single-stranded DNA (ss-DNA) conformations, we have performed room temperature molecular dynamics simulations of an aqueous solution containing the ss-DNA dodecamer with the 5'-CGCGAATTCGCG-3' sequence in the presence of 0.2, 0.5, and 0.8 M NaCl. Our calculations reveal that in the presence of the salt, the DNA molecule forms more collapsed coil-like conformations due to the screening of negative charges along the backbone. Additionally, we demonstrated that the formation of an octahedral inner-sphere complex by the strongly bound ion plays an important role in the stabilization of such folded conformation of DNA. Importantly, it is found that ion-DNA interactions can also explain the formation of non-sequential base stackings with longer lifetimes. Such non-sequential base stackings further stabilize the collapsed coil-like folded form of the DNA oligomer.

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

Flexibility of short DNA helices is important for the biological functions such as nucleosome formation and DNA-protein recognition. Recent experiments suggest that short DNAs of tens of base pairs (bps) may have apparently higher flexibility than those of kilo bps, while there is still the debate on such high flexibility. In the present work, we have studied the flexibility of short DNAs with finite-length of 5-50 bps by the all-atomistic molecular dynamics simulations and Monte Carlo simulations with the worm-like chain model. Our microscopic analyses reveal that short DNAs have apparently high flexibility which is attributed to the significantly strong bending and stretching flexibilities of ∼6 bps at each helix end. Correspondingly, the apparent persistence length lp of short DNAs increases gradually from ∼29 nm to ∼45 nm as DNA length increases from 10 to 50 bps, in accordance with the available experimental data. Our further analyses show that the short DNAs with excluding ∼6 bps at each helix end have the similar flexibility with those of kilo bps and can be described by the worm-like chain model with lp ∼ 50 nm.

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.

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.

Excluded volume and ion-ion correlation effects on the ionic atmosphere around B-DNA: theory, simulations, and experiments

The ionic atmosphere around a nucleic acid regulates its stability in aqueous salt solutions. One major source of complexity in biological activities involving nucleic acids arises from the strong influence of the surrounding ions and water molecules on their structural and thermodynamic properties. Here, we implement a classical density functional theory for cylindrical polyelectrolytes embedded in aqueous electrolytes containing explicit (neutral hard sphere) water molecules at experimental solvent concentrations. Our approach allows us to include ion correlations as well as solvent and ion excluded volume effects for studying the structural and thermodynamic properties of highly charged cylindrical polyelectrolytes. Several models of size and charge asymmetric mixtures of aqueous electrolytes at physiological concentrations are studied. Our results are in good agreement with Monte Carlo simulations. Our numerical calculations display significant differences in the ion density profiles for the different aqueous electrolyte models studied. However, similar results regarding the excess number of ions adsorbed to the B-DNA molecule are predicted by our theoretical approach for different aqueous electrolyte models. These findings suggest that ion counting experimental data should not be used alone to validate the performance of aqueous DNA-electrolyte models.

Molecular dynamics simulation of a single-stranded DNA with heterogeneous distribution of nucleobases in aqueous medium

The DNA metabolic processes often involve single-stranded DNA (ss-DNA) molecules as important intermediates. In the absence of base complementarity, ss-DNAs are more flexible and interact strongly with water in aqueous media. Ss-DNA-water interactions are expected to control the conformational flexibility of the DNA strand, which in turn should influence the properties of the surrounding water molecules. We have performed room temperature molecular dynamics simulation of an aqueous solution containing the ss-DNA dodecamer, 5'-CGCGAATTCGCG-3'. The conformational flexibility of the DNA strand and the microscopic structure and ordering of water molecules around it have been explored. The simulation reveals transformation of the initial base-stacked form of the ss-DNA to a fluctuating collapsed coil-like conformation with the formation of a few non-sequentially stacked base pairs. A preliminary analysis shows further collapse of the DNA conformation in presence of additional salt (NaCl) due to screening of negative charges along the backbone by excess cations. Additionally, higher packing of water molecules within a short distance from the DNA strand is found to be associated with realignment of water molecules by breaking their regular tetrahedral ordering.