
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39232063
70641
10.1038/s41598-024-70641-x
Article
Sensing the structural and conformational properties of single-stranded nucleic acids using electrometry and molecular simulations
Walker-Gibbons Rowan 1
Zhu Xin 1
Behjatian Ali 1
Bennett Timothy J. D. 1
Krishnan Madhavi madhavi.krishnan@chem.ox.ac.uk

12
1 https://ror.org/052gg0110 grid.4991.5 0000 0004 1936 8948 Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford, OX1 3QZ UK
2 grid.4991.5 0000 0004 1936 8948 The Kavli Institute for Nanoscience Discovery, Sherrington Road, Oxford, OX1 3QU UK
4 9 2024
4 9 2024
2024
14 2058226 3 2024
20 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Inferring the 3D structure and conformation of disordered biomolecules, e.g., single stranded nucleic acids (ssNAs), remains challenging due to their conformational heterogeneity in solution. Here, we use escape-time electrometry (ETe) to measure with sub elementary-charge precision the effective electrical charge in solution of short to medium chain length ssNAs in the range of 5–60 bases. We compare measurements of molecular effective charge with theoretically calculated values for simulated molecular conformations obtained from Molecular Dynamics simulations using a variety of forcefield descriptions. We demonstrate that the measured effective charge captures subtle differences in molecular structure in various nucleic acid homopolymers of identical length, and also that the experimental measurements can find agreement with computed values derived from coarse-grained molecular structure descriptions such as oxDNA, as well next generation ssNA force fields. We further show that comparing the measured effective charge with calculations for a rigid, charged rod—the simplest model of a nucleic acid—yields estimates of molecular structural dimensions such as linear charge spacings that capture molecular structural trends observed using high resolution structural analysis methods such as X-ray scattering. By sensitively probing the effective charge of a molecule, electrometry provides a powerful dimension supporting inferences of molecular structural and conformational properties, as well as the validation of biomolecular structural models. The overall approach holds promise for a high throughput, microscopy-based biomolecular analytical approach offering rapid screening and inference of molecular 3D conformation, and operating at the single molecule level in solution.

Keywords

Single stranded nucleic acids
Disordered molecules
Molecular electrostatics
Molecular dynamics
Subject terms

Physical chemistry
Biophysical chemistry
http://dx.doi.org/10.13039/100019180 HORIZON EUROPE European Research Council 724180 http://dx.doi.org/10.13039/501100000268 Biotechnology and Biological Sciences Research Council BB/W017415/1 Zhu Xin issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

The structural and conformational properties of nucleic acids have been thoroughly examined in a range of contexts using a variety of experimental techniques and we now have a detailed picture of their atomic structure, surrounding ion atmospheres in solution, the dynamics of hydration water and the complexes they form with proteins1–6. Whilst double-stranded (ds) DNA has been extensively studied due to its obvious importance as the “genetic blueprint”, far less is known about the structural properties of single stranded nucleic acids (ssNAs) in comparison. This is in large part due to the fact that ssNAs are flexible polymers, with an inherently disordered structure, and exist as heterogeneous conformational ensembles that may not be well represented by single ensemble-averaged 3D structures. Unlike folded biomolecules, e.g., globular proteins, ribozymes and riboswitches, disordered biomolecules often do not lend themselves well to high-resolution structural techniques such as X-ray crystallography and cryo-electron microscopy7,8. Yet knowledge of the molecular level conformation and local atomistic arrangement of ssNAs is important to characterise their interactions, e.g., protein-RNA interactions are important for CRISPR-Cas9 technologies and in driving phase separation and modulating the properties of biomolecular condensates9–12. Furthermore, the conformational properties of biologically significant homopolymeric motifs such as U/A tracts in ssRNA serve important functions in protein recognition, the regulation of gene expression, and the assembly of virus like particles13–18. ssDNA also finds extensive use as a therapeutic agent for a host of diseases, a platform for the development of vaccines and as a polymer in nucleic acid origami nanotechnology19–22. The conformational richness displayed by ssNAs holds the key to their diverse roles and interactions. The challenges and opportunities that lie in the problem of structure-conformation inferences on disordered biomolecules heightens the importance of developing new and orthogonal measurement modalities that in conjunction with other techniques may contribute to multi-parametric structural and conformational analyses15,23.

Experimental investigations into the structure of ssNAs thus far have largely relied upon techniques such as small angle X-ray scattering (SAXS), nuclear magnetic resonance (NMR) and single molecule Förster resonance energy transfer (smFRET)15,24–27. SAXS is a label-free method that has proven very powerful in examining ssNA conformation, revealing both local and global sequence dependent chain conformations that relate to base stacking affinities and the local ion atmospheres around nucleic acids15,23,24. A range of experimental studies have uncovered clear salt- and sequence-dependent trends in conformational properties, e.g., the persistence length, lp, and radius of gyration Rg, of ssDNA were observed to decrease with an increasing salt concentration, and stiffer polymer conformations were inferred for poly-dA structures compared to -dT which was attributed to stronger inter-base stacking between adenine bases compared to thymine26,28–31.

We recently pioneered escape time electrometry (ETe)—a field free, optical microscopy-based technique that measures electrostatic interaction energies between a charged molecule and a like-charged surface with high precision. Because the magnitude of electrostatic interactions carries a strong signature of molecular shape or 3D conformation, such measurements can offer a fundamentally different approach to inferring molecular structural information in solution32–34. The ETe method relies on the electrostatic fluidic trap to spatially confine charged molecules diffusing in solution in a geometrically tailored slit created by two like-charged parallel surfaces (Fig. 1A,C). Measurement of the time spent by individual molecules in the trapped state furnishes a high precision inference of the interaction energy between the molecule and the flanking slit surfaces, which in turn can be parametrized in terms of an effective charge, qeff, of the molecular species, as previously described32–38. The effective charge of a molecule in solution is often smaller than the molecule’s structural charge, qstr, (given by the sum of charges carried by ionized structural groups and any structurally bound ions from the electrolyte), and its departure from the latter is captured by a renormalization factor given by η=qeffqstr. The greater the extent of charge renormalization in the molecule the smaller the value of η. For a highly ionized (acidic) object such as a nucleic acid molecule, η is solely a function of the geometry of the molecular charge distribution or the 3D conformation of the molecule33,35,39,40.Fig. 1 Measurement of the effective charge of ssNAs using ETe. (A) Schematic representation (not to scale) of ssNA molecules confined in an array of electrostatic fluidic traps and imaged using wide-field fluorescence microscopy. (B) A time course of optical images in a single trap displays the duration of a single recorded residence event of duration, Δt. (C) Graphical representation of the cross- section of a fluidic slit of height 2h≈70 nm, with a pocket depth d≈210 nm (top). Below: Calculated electrostatic potential distribution, ϕ, presented on logarithmic scale for a salt concentration c0= 1.2 mM and h= 35 nm (left half-space). Line plot of midplane potential value, ϕmid, in the vicinity of a single pocket (right half-space), displays a value Δϕmid=ϕm≈ −0.06 kBT/e in the slit region, corresponding to an effective surface potential of ϕs≈ −2.2 kBT/e. (D) A representative electrostatic energy landscape (ΔFel) for a single trap in this study showing the spatial states of the molecule in the slit (out: “1”) and pocket (in: “2”). (E) Probability density distributions PΔt of escape times obtained from 3–5 min of imaging for short chain ssDNA fragments (bottom). All measurements were performed in devices with slits of depth 2h= 70 nm and surface potential value ϕs= -1.91kBT/e, under in an electrolyte containing NaCl at a concentration c0= 0.6 mM (κh ≈ 5.5) and pH = 7.5. Measured tesc values for oligomers of length nb= 5–12 reveal clear differences (top panel). The average molecular escape times, tesc, measured effective charge values, qeff, and number of escape events, N, are as follows: nb=5: tesc= 29.16±0.38 ms (qeff=-6.3±0.082e, N=2.2×104), nb=8: tesc= 50.82±0.43 ms (qeff=-9.2±0.077e,N=3.7×104), nb=9: tesc= 65.68±0.89 ms (qeff=-9.8±0.13e,N=2.4×104), nb=11: tesc= 79.24±1.30 ms (qeff=-10.7±0.175e,N=8.0×104), nb=12: tesc= 98.08±2.3 ms (qeff=-11.9±0.28e,N=3.0×104). Widths of presented Gaussian distributions of tesc arise predominantly from the statistical error on the fitted tesc values which are given approximately by tesc/N (see “Methods”).

Our definition of molecular effective charge is founded in electrostatic interaction free energies as previously described, and is in fact in good agreement with a range of other estimates of qeff based on different but related considerations35,40–42. For example qeff as defined in our work agrees reasonably well with analytical estimates for charged spheres and cylinders given in Ref.40 which we invoke here in order to provide physical intuition and understanding of the background of our approach. For charged cylinders of radius R and length l≫κ-1 carrying a uniformly distributed charge of qstr, which describes polyelectrolytes such as DNA, Ref.40 suggests1 qeffe=llB1-ln[qstre(lBl)]2ln(κR)

in the regime κR<1 which corresponds to a typical polyelectrolyte immersed in an electrolyte of ionic strength c0≈ 1 mM. Here e denotes the elementary charge, and κ-1≈0.304/c0 nm is the Debye length which characterizes the spatial rate of decay of the electrostatic potential from a charged interface in an electrolyte containing only monovalent ions. We also have lB≈0.704 nm, the Bjerrum length in water, which denotes the spatial separation at which two identical charges experience kBT worth of interaction energy, where kB is Boltzmann’s constant and T is the absolute temperature. Equation (1) shows that (1) the renormalized linear charge density tends to the Manning limit of 1/lB as κR→0 which e.g., describes a line charge, and (2) since l=nbb, where nb is the number of charged monomers and b the spacing between them, a precise measurement of qeff can be sensitive to the charge spacing of a uniform polyelectrolyte. We recently demonstrated measurements of the rise per basepair and radius of two forms of the DNA double helix using this principle, and will demonstrate the ability to make structural inferences on ssNAs in this study33,34,38. In general, as indicated in Eq. (1), the more compact the 3D charge distribution—or the higher the linear charge density of the molecular backbone in the case of nucleic acids—the greater the influence of charge renormalization, the lower the magnitude of qeff and hence the lower the value of η35.

It is also worth connecting the interpretation of our definition of qeff with readouts of experimental techniques like atomic emission spectroscopy (AES) and anomalous SAXS (ASAXS) which infer an excess concentration of counterions of type i, Γi, in the ion atmosphere, surrounding a charged molecule relative to the bulk solution43. Both qeff and Γi provide descriptions of charge renormalization and derive from the same underlying physical quantity: the electrical potential distribution around a molecule in solution. In order to make a direct link between these two closely related quantities we have conducted an analysis for a negatively charged rod under the conditions of our present experiments and suggest the approximate relationship Γi∝a1+a2η, where a2<0 (see Supplementary Information Sect. S7). In general the parameters a1 and a2 in this equation depend on molecular structural charge, geometry and salt concentration. Importantly the given relationship captures the qualitative expectation that a lower magnitude of qeff (smaller η) entails a larger number of excess counterions in the molecule’s ion atmosphere, and will be often invoked in comparative discussions of results in this study with those from ion counting methods like AES.

Previous studies using electrometry have shown that measured values of qeff agree remarkably well with effective charge values calculated using Poisson-Boltzmann (PB) theory across a diverse range of molecules displaying distinct molecular conformations32–35. PB theory is widely known to be successful in describing the electrostatics of DNA, for instance in capturing the fraction of bound counterions and the nature of the surrounding ion distribution, predicting salt effects on the binding affinities of ligands and proteins to nucleic acids, and in modelling DNA adsorption behaviour to biological membranes and DNA-DNA interactions in supercoils44–52. A recent ETe investigation of DNA nanostructures of identical molecular weight but different 3D conformation illustrated the intimate relationship between 3D molecular conformation and the far-field electrostatic signature of a molecule, permitting various molecular 3D conformations to be readily differentiated according to their measured qeff values33. A related study focussing on the nucleic acid double helix demonstrated that high-precision ETe could further provide a route to determining higher resolution structural information such as the radius and the average rise per basepair of A- and B-form double helices characteristic of dsDNA and dsRNA. These measurements also furnished an experimentally inferred view of the interfacial hydration layer and the arrangement of ions at the molecular interface with the surrounding electrolyte34.

Following in the vein of high-resolution structural modelling techniques that compare experimental readouts with expectations based on computational modelling, in this study we undertake a comparison of measured qeff values with the theoretically expected effective charge, qcalc, for various ssNA species. In particular we employ a range of all-atom force fields including the standard Amber and CHARMM models and some next generation specialised ssNA forcefields such as DES-Amber, as well as a coarse-grained model—oxDNA. In order to relate our measurements of qeff to ssNA structural properties, we select a subset of representative model structures from each simulated conformational ensemble for nucleic acid molecules, and calculate their theoretically expected effective charge, qcalc, using our previously described PB free energy framework35,38. We then compare our measurements with modelling results and discuss the performance of each model and forcefield in capturing both our measurements as well as trends reported using other experimental techniques such as SAXS and AES.

The first part of the study focuses on the experimental approach, demonstrating the ability to measure molecular effective charge on small ssNA fragments. We start by probing an experimental regime where there is not much influence of molecular 3D conformation on qeff and show that this is indeed reflected in computationally modelled values based on MD conformations of short fragments. Next we perform ETe on longer ssNA fragments (30 and 60 base length) and show that differences in effective charge between various homopolymeric stretches of NAs clearly emerge in measurement. In comparing measured effective charge values with those computed from molecular models, in general we found that qeff carried contributions both from global conformation features and local structural information such as charge density along the molecular backbone. Whilst this is expected on qualitative grounds, a systematic quantitative delineation of the relative contributions is likely to be highly challenging. We also found that absolute comparisons between qeff and qcalc are limited by our current level of experimental inaccuracy, which we estimate at about 10% as described later. We therefore do not attempt at this stage to perform ensemble optimization and refinement on the modelled molecular conformations in order to achieve agreement with experiments, as is more common practice in studies involving established structural techniques such as SAXS15,23. Nonetheless, relative trends in qeff values for different species measured under the same conditions are robust to measurement inaccuracies and could be used for comparisons with corresponding trends in computational models. In general our results will show that we did not obtain perfect agreement with any particular molecular model, although we did note that specific models such as DES-Amber and oxDNA performed well in capturing certain measurement trends. Finally, we also demonstrate that line charge models may provide a simple route to parametrising measured differences between molecular species.

We begin by briefly describing the ETe method to measure electrostatic interaction energies and infer corresponding effective charge values for ssNAs in solution. We measure the average residence time or escape time, tesc, of a molecular species in a landscape of electrostatic traps created by the geometric modulation of the gap between two charged plates in solution (Fig. 1). The escape time depends exponentially on the effective charge of the molecule, qeff, which can be then be inferred from measured tesc values as described previously32–34,36,37,53. We work in the electrostatic regime corresponding to κh≈ 5.5–8. Here 2h=70-75 nm is the separation between the charged surfaces in the slit region of an ETe measurement device (Fig. S2). The salt concentration c0 lies in the range 0.6-1.2 mM in this work.

We consider a series of oligomers of length between 5 and 12 bases, and two longer fragments of length 30 and 60 bases. We explore how the molecular effective charge of ssNAs in solution varies as a function of both number of bases, nb, as well as sequence. In particular, for nb = 30 and 60, we examine ssNA species carrying the same amount of total structural charge but differing in nucleotide sequence, namely poly-dT, -dA, -rU and a mixed sequence case. In tandem we use a PB interaction free energy framework described previously to calculate effective charge values, qcalc, for model structures of ssNAs generated using Molecular Dynamics (MD) simulations. We then compare measurements of qeff with the computed results. The conformational properties of the homopolymeric motifs in particular have been extensively studied using different experimental techniques and therefore serve as well characterised test cases with which to compare our results. Comparing our measurements with computationally generated qcalc values, we demonstrate how ETe measurements of the effective molecular charge in solution can serve to infer aspects of molecular structure and conformation, and may also provide an additional physical dimension using which to validate structural models of charged biomolecules in solution.

Results and discussion

Molecular effective charge values for short chain length ssNA fragments

We first measured the escape dynamics of a series of short chain length ssDNA fragments (nb= 5–12) diffusing in a lattice of electrostatic fluidic traps (see Fig. 1). The oligomeric species displayed distinct average escape times, tesc, with the smallest nb = 5 fragment exhibiting an escape time of tesc≈ 29 ms and the longest fragment with nb = 12 a value of tesc≈ 98 ms (Fig. 1E). The measured values of tesc were then converted to values of molecular effective charge, qeff (see “Methods” and Refs.32–34,36,37,53). A calibration measurement using 60 basepair dsDNA—a species that has been previously thoroughly characterized in electrometry—was used to determine the value of the surface potential of the slit, ϕs, which is a priori unknown (see “Methods”)32,34.

We then compared our measurements with effective charge values for simulated molecular structures, qcalc, calculated as described previously33–35,38,54. Note that each molecular species (apart from the nb = 5 case) is characterised by a value of qstr=-nb+3e, which results from accounting for the additional charges of the attached fluorescent dyes (see “Methods” and Figs. S1, S3). We found close agreement (within 5% on average) between the effective charge values computed for model ‘half-helix’ structures, qcalc,HH, (orange data points) and the experimentally measured data (blue data points) (see Fig. 2A). Performing 0.5 μs long MD simulations (employing the DES-Amber forcefield) on these short chain molecules permitted their configurations to depart from their initial helical structure, yielding a conformational ensemble of structures characterised by a spectrum of Rg values and therefore a range of calculated qcalc,MD values. Calculations of the effective charge for conformations corresponding approximately to the minimum and maximum Rg values obtained from MD simulations (see “Methods” and Fig. S2 for details) revealed a small amount of variation in their qcalc,MD values, reflecting the conformational heterogeneity (associated error bars of the MD structure data points in Fig. 2A). In order to enable a quantitative comparison of calculated values with experimental readouts, the effective charge contributed by the covalently attached fluorescent dye molecules was added a posteriori to the effective charge calculated for dye-free molecular conformations (see “Methods” and Fig. S3). Importantly, for the series of short nucleic acid molecules, we found that the measurements revealed |qeff|≈|qstr|. In the regime where a typical length scale describing the extent of the molecule, such as the polymer contour length lc, is smaller than the Debye length, i.e., lc≈5 nm <κ-1≈10 nm, the ‘far-field’ electrostatic potential distribution resembles that of a sphere and charge renormalization vanishes, as recognised in previous theoretical work39,40. Under our present experimental conditions (κlc< 1), the conformational properties of very short fragments of ssDNA are not expected to significantly impact the measured charge renormalization values, therefore implying η≈ 1. Indeed, for 5-base ssDNA, for example, we measure η = 0.94. Thus for short nucleic acids our present measurements are rather insensitive to the precise 3D arrangement of the charge in the molecule, since under the conditions of the measurement, the spatial extent of the molecular charge distribution under consideration is smaller than the Debye length. Figure 2B displays the contours of electrostatic potential around molecular structures corresponding approximately to the minimum and maximum Rg for an nb=12 polymer from MD, which can be seen to resemble one another a short distance away (∼0.5κ-1) from the molecule, with the potential distribution appearing similar to that of a small sphere for smaller values of nb. The two conformations are characterised by a difference in qcalc of 0.7 e or less than 6% of the mean value of 12.6 e. For this reason, qcalc values calculated for model structures consisting of only a half-helix, without any higher-order structural information, find good agreement both with MD-simulated conformations as well as with the experimental data. Thus, although little molecular structural or conformational insight can be gleaned in this electrostatic regime (lc<κ-1), the results demonstrate sensitivity of the measurement to the addition of just a single nucleobase. Note that the measurement of 100–1000 molecules per species in this study yields high-precision measurements of qeff (measurement imprecision ≈ 0.1–1%), which should enable small structural differences between species to be discerned, as described later. Although highly sensitive to relative differences between species, we expect that absolute values of the measurements reported here may reflect an inaccuracy of up to 10%, which arises from a number of different experimental sources such as sample purity, and the influence thereof on the fitting procedure to extract tesc values, as well as accuracy limits on the value of qeff of the calibration molecule (see “Methods”). This can however be improved upon in future such that the accuracy of the method better reflects the precision. The inaccuracy of ≈ 10% on qeff in this work nonetheless does not preclude comparisons with computational estimates of the same quantity under conditions where qeff and qcalc are more sensitive to the 3D charge distribution in the molecule, i.e., lc>κ-1, as described in the following section.Fig. 2 Molecular modelling of short chain ssNAs (nb= 5–12) and comparison of experimental measurements of molecular effective charge (|qeff|) with calculations (|qcalc|). (A) Plot of the molecular effective charge |qeff| as measured by ETe (blue symbols), calculated for model half-helix structures (|qcalc,HH|, orange symbols) and for model conformations taken from MD simulations (|qcalc,MD|, red symbols) vs. molecular structural charge |qstr|. The simulated structures are free of dye molecules but the effective charge of two ATTO dye molecules is added post-calculation to enable comparison with the experimental measurements (see “Methods” and Figs. S1, S3). Error bars on |qeff| arise predominantly from the statistical error on the fitted tesc values. The error bars for the MD data points represent extrema values of the effective charge of structures across the Rg range of the molecule observed in simulations. Dashed grey line represents η=1, presented for reference. Lower panel: Plot of the percentage difference, r, between calculations of effective charge for model half-helix structures, |qcalc,HH|, and experimental measurements |qeff|. (B) Representative electrostatic potential distributions surrounding an nb= 12 oligomer (|qstr|= 15 e) obtained by solving the PB equation for 3D atomic charge distributions of conformations corresponding to extrema values of Rg= 1.65 nm (top) and Rg= 1.1 nm (bottom) taken from an MD simulation using the DES-Amber forcefield. Electrostatic potential contours are shown at ϕ= − 5, − 2.5, − 1.25 and − 0.625 kBT/e. The difference in |qeff| between the two structures is ≈ 0.7 e or < 6% of the mean value of 12.6 e.

Sequence dependence of effective charge values for longer chain ssNAs

Inspired by a range of X-ray scattering studies, smFRET, and force extension measurements that have reported structural differences between homopolymeric sequence ssNAs composed of an identical number of bases, we turned our attention to examining differences in homopolymeric ssNAs. The precision offered by ETe and the sensitivity of the underlying physics to the 3D charge distribution in the molecule may indeed offer the prospect of using effective charge to detect sequence dependent differences in molecular structure and/or conformation. Note that when measurements on various molecular species are performed under identical conditions, as defined by the same measurement device and solution conditions, the measurement precision (≤ 1% in this study) effectively becomes the relevant quantity in determining the ability of the approach to discriminate between different molecular states, conformations or species.

We performed experiments on ssNA species of two different lengths (nb= 30 and nb = 60) and explored the effect of nucleic acid sequence alone on the measured effective molecular charge. The measurement regime corresponds to lc>κ-1 where κlc≈ 1.7 for 30 base and κlc≈ 3.4 for 60 base ssDNA, assuming an approximate contour length per base value bc = 0.5 nm. In this regime we expect substantial charge renormalization (η<1), in contrast to the measurements performed on short ssNA fragments. Indeed we found significant charge renormalization corresponding to η≈ 0.6–0.8 in general (e.g., η≈ 0.8 and η≈ 0.7 for 30 and 60 base poly-dT respectively), and differences in measured qeff of approximately 5% between 30 base poly-dT and poly-dA, and 10% between 60 base poly-dT and poly-dA. For instance, the measured difference in qeff between poly-dT and poly-dA homopolymers is ≈ 2 e and ≈ 4.5 e for the nb= 30 and nb = 60 length chains respectively. These values reflect inter-species disparities in effective charge that are much larger than measurement precisions of approximately 0.5% and 1% for nb = 30 and nb = 60 respectively. Importantly, the measured order of qeff follows poly-dT > mixed sequence > -rU ≥ -dA and implies different average molecular conformations of the homopolymers (Fig. 3A,B). We attribute the sensitivity of the readout to molecular structural detail to the fact that in the regime lc>κ-1, characterising 60 base fragments, higher order multipoles of the molecular charge distribution make a greater contribution to the electrical potential in the far-field, i.e., at large distances from the molecule, h>3κ-155,56. This accentuates the sensitivity of the electrical free energy measurement to molecular 3D conformation and highlights the ability of the ETe technique to sense spatial features of molecular charge density distributions.Fig. 3 ETe detects composition differences between ssNA oligos. (A) Probability density distributions PΔt of residence times for measurements on single stranded homopolymers poly-dT (blue), -dA (red), a mixed DNA sequence (purple) and RNA poly-rU (orange), shown here for the fragment length nb=60, and for the calibrator molecule, 30dsDNA (grey) (see “Methods”). Atomic representations of the nucleotide bases are shown as insets. (B) Measured tesc values for the data presented in (A) (upper panel). Two datasets of |qeff| values inferred from the measured data (solid and dashed vertical lines), for the nb=60 ssNA polymers (lower panel). Widths of presented Gaussian distributions of tesc and error bars on |qeff| arise predominantly from the statistical error on the fitted tesc values in (A). All measurements were performed in devices with slits of depth 2h = 75 nm and surface potential value ϕs = −2.24kBT/e, under similar electrolyte conditions NaCl = 1.2 mM, pH = 7.5. The average molecular escape times, tesc, measured effective charge values, qeff, and number of escape events N for the given measurement dataset are as follows: poly-dT: tesc= 372.01±2.13 ms (qeff=-44.3±0.0923e,N=1.1×105), poly-dA: tesc= 265.79±2.62 ms (qeff=-40.7±0.165e,N=6×104), mixed sequence: tesc= 323.73±6.864 ms (qeff=-42.3±0.342e,N=5.2×103), poly-rU: tesc= 293.11±2.94 ms (qeff=-41.5±0.165e,N=5.8×104). (C) Gel electrophoresis of ssNA samples for nb = 30 and 60 (see Fig. S6 for details). The electrophoretic mobility follows poly-dA > mixed sequence >  -dT >  -rU.

We further noted good qualitative agreement in comparing trends in molecular effective charge and conformation deduced using electrometry with those inferred from other experimental techniques. The measured qeff values for poly-dA and poly-rU are significantly smaller than that of poly-dT indicating higher density charge distributions that are subject to greater charge renormalization (smaller η), compared to poly-dT (e.g. qeff= 44.1 e, 38.7 e and 40.7 e for poly-dT, -dA and -rU, corresponding to η values 0.70, 0.61, 0.65 respectively in a given measurement dataset)24. Techniques such as ASAXS and atomic emission spectroscopy (AES) directly quantify the number of excess ions in the ion atmosphere surrounding nucleic acids and therefore also report on a molecule’s effective charge in solution43,57–59. Plumridge et al. have conducted AES measurements of NA homopolymers and reported values for the excess number of Na+ counterions per phosphate to be ≈ 0.68, 0.71 and 0.83 for poly-dT30, -dA30 and -rU30 respectively, implying a smaller magnitude of qeff for poly-dA and -rU compared to poly-dT (see Eq. (1) and Supplementary Information Sect. S7)15,60. In accordance with these results, SAXS measurements found that poly-dT polymers were associated with the fewest excess ions whereas other oligonucleotides considered, which had T bases substituted with A, G or C, were all associated with a larger number of excess counterions61. Furthermore, an AES study comparing 50 base poly-rU and -dT ssNA homopolymers implied a smaller magnitude of qeff for RNA compared to DNA, and attributed this disparity to an increased polymer charge density brought about by the smaller inter-charge spacing for RNA62. We also note here that the counterion excess Γ+ as reported by AES measurements can be less sensitive to the charge density of a rigid rod model of ssDNA than qeff, particularly at low ionic strengths and in relation to the measurement precision of each method (see Supplementary Information Sect. S7 for further details). For example, the magnitude of the difference in Γ+ between poly-dA and poly-dT of around 4% in AES measurements can be seen to correspond to > 10% difference in η (see Fig. S7)60. Furthermore, our measured differences in qeff between molecular species are much larger than the typical measurement imprecision and immediately suggest differences in the 3D distribution of charge in the various homopolymers that are qualitatively consistent with expectations from ion counting studies. In summary, ETe measurements of molecular effective charge in solution indicate that |qeff| values follow the ordering poly-dT > mixed sequence > -rU ≥ -dA and are in agreement with AES studies. The ssNAs under consideration also display clear but different qualitative trends in their electrophoretic mobilities in polyacrylamide gels which may be more challenging to interpret and are discussed further in Supplementary Information Sect. S5 (Fig. 3C). Furthermore, inferences on molecular conformation from ETe in combination with molecular modelling analyses are in some cases supported both by experimental smFRET and SAXS measurements, as described in the following section24.

Comparing electrometry measurements with the results of computational molecular modelling

General considerations

Because the ETe readout is highly amenable to computational modelling, we embarked on an analysis that sought to quantitatively compare and relate the measurements with predictions from molecular modelling approaches. We therefore generated a host of representative ssNA structures from MD simulations employing various forcefield descriptions. The conformational free energy landscape of the molecule obtained from MD simulations was characterised by the radius of gyration Rg and end-to-end distance parameter R (Fig. 4B). We then determined qcalc values for a pool of representative structures using our PB electrostatics framework (see “Methods”). We typically worked with a subset of structures for each case that included the most likely conformation as well as structures corresponding to the maximum and minimum R values on the 25th percentile contour of structures in (Rg,R) parameter space (black crosses in Fig. 4B), furnishing a range of possible effective charge values for a given molecular species. Although the measured sequence-dependent differences in qeff are more significant for nb= 60 compared to nb= 30, we restricted all-atom simulation studies to 30-base oligomers owing to the significantly lower computational cost.Fig. 4 Molecular modelling and calculation of theoretical effective charge values for long chain ssNAs (fragment lengths nb= 30 and 60) (A) Molecular model of the initial ‘half-helix’ configuration used for MD simulations employing conventional all-atom forcefields Amberbsc1 and CHARMM36, specialised ssNA DES-Amber and HBCUFIX forcefields and the coarse grained oxDNA and oxRNA models. (B) Simulated conformational landscapes of ssNA species (poly-dT—blue, poly-dA—red, poly-rU—orange), parametrised in terms of the polymer radius of gyration, Rg, and end-to-end distance, R. Each contour encloses a 12.5th percentile of structures. R vs. Rg landscapes are shown for the DES-Amber (left) and oxDNA models (right) for nb= 30 and nb= 60 respectively (inset: representative structures with the most likely R, Rg combination). The R vs. Rg landscapes for ssNAs modelled with the other forcefield descriptions are shown in Fig. S5. (C) Comparison of the effective charge values, |qcalc|, for long chain ssNA structures (nb= 30 and 60) obtained via molecular modelling approaches, with experimentally measured |qeff| data (open circles). Measurements (left column, circular symbols and associated measurement precision error bars) represent the average values of two measurement datasets. Horizontal shaded bands depict the range of qeff values across two measurement datasets. qcalc values of representative model ssNA structures (right columns) carry upper and lower bounds that correspond to structures with the maximum and minimum R values contained within the most likely 25th percentile of structures (corresponding to black crosses in (B)). (D) Contour length per base, bc, values calculated from molecular simulations employing different forcefield descriptions, and where bc is calculated as the average inter-phosphate distance for all atom models or as the average inter-backbone site distance for oxDNA models. Black symbols denote values for both poly-dT and poly-dA sequences that share the same bc value.

In keeping with qualitative expectations, we found in general that model structures with larger values of (Rg,R) gave rise to a larger magnitude of qcalc, and vice versa (error bars in Fig. 4C display the range of calculated values for each case). Despite the fact that we are only able to report all-atom simulation results for nb= 30 where the sensitivity of effective charge to molecular properties is poorer under the present experimental conditions, and measured disparities between the molecular species are therefore smaller, it must be noted that percentage changes in qcalc are small compared to the corresponding percentage changes in Rg and R (e.g. ≈ 10% and 20% change in Rg and R respectively resulted in a 1.5% difference in qcalc for poly-dT modelled with DES-Amber). Overall, we found that model structures for different molecular species, e.g., poly-dT, -dA and -rU, that are characterised by very similar values of (Rg,R), also had very a similar qcalc. However, we also found that this is not always the case, and there were instances where we noted more substantial differences (≈ 2–3%) in qcalc for model structures of different species characterised by the same (Rg,R). Therefore, the ETe readout appears to carry information reflecting both the global properties of the charge distribution, as captured for instance by (Rg,R), as well as local properties such as the charge density along the backbone of the molecule, and it is challenging to clearly quantify the relative magnitude of these contributions to the overall molecular effective charge.

Molecular dynamics simulations: the Amber and CHARMM forcefields

Of the available classical all-atom models for simulating nucleic acid dynamics, the parmbcs1 model of the Amber forcefields and the CHARMM36 forcefield represent the state of the art, and have undergone significant rounds of improvement since their inception63. However, these forcefields were originally parametrised for dsDNA and are therefore known to have shortcomings in modelling ssDNA21,64. Here we discuss theoretical molecular effective charge values obtained for ssNA structures taken from all-atom MD simulations performed with the Amber and CHARMM forcefields for nb = 30.

As noted previously, simulations of poly-dT and -dA30 performed with the Amber-bsc1 forcefield, generated structures did not depart from the initial half-helix configuration (Fig. S4)21. This outcome is in contradiction with the SAXS and smFRET experimental data that poly-dT is a somewhat disordered polymer chain, even under low salt conditions21,64. Overall, the simulations indicated very similar conformational properties for poly-dT and poly-dA as well as comparatively narrow Rg vs R distributions for both (Figs. S4, S5). This is reflected in the similarity of their corresponding calculated molecular effective charge values, qcalc≈ 22.6± 0.1 e (blue and red square symbols in Fig. 4C). This result is in obvious qualitative disagreement with our experimentally measured qeff values that displays clear differences of approximately 8% between poly-dT and poly-dA (open circles in Fig. 4C).

The CHARMM36 forcefield, on the other hand, appeared to perform much better at capturing the conformational properties of ssNAs in general. For instance, the results compare more favourably with the experimentally reported Rg and R values for poly-dT30, as well as with inferences on the ‘local’ polymer chain properties from SAXS measurements at 20 mM NaCl (see Supplementary Information Sect. S4 for further discussion)15,23,43. Interestingly, simulations with the CHARMM forcefield revealed a marked difference in conformational ensembles for poly-dT and poly-dA sequences, with the poly-dA structures characterised by collapsed conformations that resided in metastable hairpin structures, as reflected in a small average end-to-end distance (Rg≈2.3 nm and R≈3 nm), well below the reported SAXS values of Rg=2.72 and R=7 nm for poly-dA30 (Fig. S5)23. However, such globular forms of poly-dA have not been observed in SAXS experiments up to salt concentrations as high as 1 M and possibly points to a shortcoming of the CHARMM36 forcefield in modelling poly-dA28. Indeed we found that the simulated conformations for poly-dA, and to a lesser extent for poly-dT, gave calculated values of molecular effective charge that fell significantly short of those measured by ETe by about 20% (blue and red triangle symbols in Fig. 4C). This is not surprising, given that conformational compactness generally entails higher charge renormalization (η< 1) and therefore a lower magnitude of effective charge. Nonetheless, despite the fact that the absolute qcalc values were lower by about 20% than the measured qeff values (e.g., poly-dT30 |qcalc|≈21.5e, |qeff|≈27.1e), the qualitative trend indicated by simulations that poly-dT has a less compact conformation than poly-dA (≈ 7% disparity in |qcalc|) is indeed captured in experimental data that points to an ≈ 8% disparity in |qeff|. However, for poly-rU RNA, the CHARMM36 forcefield resulted in more extended structures and thus higher magnitudes of qcalc values than the DNA homopolymers (see Figs. 4, S5)—a trend clearly not observed in the qeff values measured by ETe. Such disparities between measurement and structural modelling trends was a recurring theme for the force-fields examined in this study, except for DES-Amber discussed next.

The specialised DES-Amber and HB-CUFIX forcefields

Next we analysed results from DES-Amber, a next generation ssNA forcefield. We found that DES-Amber (star symbols in Fig. 4C) was able to better capture experimental sequence dependent trends in qeff for ssNAs, i.e., qeff of poly-dT >  -rU ≥ -dA. The DES-Amber forcefield overcomes some of the deficiencies seen with the Amber and CHARMM36 models, with the ssDNA exhibiting some degree of flexibility whilst not engaging in extensive internal base-pairing that results in the DNA collapsing into globular conformations. Furthermore, DES-Amber indicates minor differences in the conformational ensembles of poly-dT and poly-dA, with poly-dT exhibiting a higher degree of polymer flexibility (able to transiently access more extended and collapsed conformations) and poly-dA adopting more rigid and compact structures with a higher degree of inter-base stacking (Fig. 4B). These trends may be seen to be reflected in the |qcalc| values, which are ≈ 22.8 ± 0.15 e and ≈ 22.6 ± 0.08 e for poly-dT and poly-dA respectively (Fig. 4C). Importantly however, this difference in effective charge between the two species of ≈ 0.2 e (≈ 1%) remains substantially lower than the experimentally measured disparity in qeff of ≈ 2 e (≈ 8%) measured by ETe. Such quantitative discrepancies between experimental readouts and DES-Amber modelling results are also evident in other studies. For example, whilst values of Rg for both poly-dA and -dT modelled with DES-Amber are in good agreement with the reported experimental SAXS value of approximately 3 nm, the DES Amber mean end-to-end length (≈ 9 nm) for both cases appears to be larger than the experimentally reported values of R≈7 nm at the lowest salt concentrations probed23,60.

Finally we turn our attention to modelling results for poly-rU. The DES-Amber forcefield was introduced as being able to model RNA successfully, having been validated to capture the end-to-end length for RNA as measured by smFRET65. Importantly, we found that the DES-Amber forcefield generated structures for DNA and RNA whose calculated |qcalc| values were indeed significantly different (see Fig. 4C). The MD conformational ensemble for DES-RNA reflects a higher degree of polymer flexibility and a much shorter average end-to-end length (Rg≈ 3 nm and R≈ 8 nm) than its DNA counterpart (see Fig. 4B). Furthermore, the average contour length of DES-RNA (calculated as the sum of inter-phosphate distances along the polymer chain) was calculated to be 17.2 nm for nb=30, considerably shorter than poly-dT, which measured 19.5 nm. This ratio in contour lengths of poly-rU to -dT is 0.88, which is in close agreement with the value of ≈ 0.87 reported using a combination of SAXS and smFRET in Ref.24, albeit for measurements performed under different conditions. We expect the suggested reduced spatial extent and inter-backbone charge spacing of the RNA polymer chain to result in a reduced molecular effective charge for RNA compared to DNA. Indeed, we found that representative configurations for 30 base poly-rU and -dT from DES-Amber have average |qcalc| values of 22.1 and 22.8 e respectively, giving a ratio of 0.97 between -rU and -dT (Fig. 4C). Despite the fact that the absolute values are similar but not identical to the measured values (|qeff|≈ 26.2 ± 0.12 e and 27.4 ± 0.1 e), the qualitative trend obtained for the simulated molecular structures, as reflected in the ratio of the two values, is in fact captured in the electrometry measurements that indicate a ratio of 0.96. Thus, experimentally measured disparities between the two molecular species are indeed successfully captured in DES-Amber simulations. It is worth noting again in this context, that whilst the ETe measurements are of high precision (≤ 1% in this work), which enables differences and relative trends in effective charge between molecular species to be measured reliably, the accuracy of the absolute values of qeff can display around 10% uncertainty in the present measurements which precludes rigorous testing against absolute qeff values for computed structures. Furthermore, minor additional uncertainties on modelling parameters enter the picture on the computational side, carrying implications for quantitative comparisons of absolute values between experiment and simulation. One source of uncertainty in the electrostatic modelling of molecular structures concerns the width of the ion accessible region, w, surrounding the molecular structure, which introduces a source of variation in the magnitude of the qcalc values that may well account for about 5% of the presently noted disparity in absolute values. For instance, increasing the value of w from 0.2 to 0.4 nm raises the magnitude of the qeff by ≈ 2% per Å for 30 and 60 base poly-dT. Currently therefore, these considerations in both experiments and computation taken together place constraints on the ability of our present methodology to clearly highlight shortcomings and gaps in molecular simulation models when disparities in absolute values of qeff between experiments and structural models lie in the range of ≈ 10% or less. Despite the noted disparities between absolute values of effective charge from measurement and computation, inferences based on measured differences or relative trends in effective charge between species measured under identical conditions are less subject to experimental inaccuracy and model-parameter uncertainty in computation. Comparative trends in effective charge therefore do nonetheless provide a suitable basis on which to compare experiments with computational indications, and may aid in the fine-tuning of molecular models.

Finally, we employed an alternative ssRNA model, HBCUFIX, which has been demonstrated to better capture the persistence length and Rg of short poly-rU chains compared to DES-Amber66. We found that the conformational landscape sampled by this molecular model was overall similar to that of DES-Amber, exhibiting a substantial exploration of the conformational parameter space, reflecting that of a flexible polymer chain (see Fig. S5). The effective charge value of a representative poly-rU structure using the HB-CUFIX model of ssRNA gave |qcalc|≈ 21.3 e, similar to the DES-Amber value of ≈ 22.3 e. Thus we found that under the conditions of this study, the specialised ssRNA forcefield descriptions generated effective charge values that were in good mutual agreement.

The coarse grained oxDNA model

Having examined and compared our experimental data with a range of atomistic MD models we then performed a comparison of the results of oxDNA2, a coarse-grained DNA model, with our experiments for both nucleic acids lengths nb= 30 and 6067. We proceeded along the same lines, relating computed effective charge values, |qcalc|, with experimental measurements as described previously. Since oxDNA was optimised to predict the self-assembly, structure and mechanical and thermodynamic properties of DNA nanostructures, e.g., DNA origami, it is not entirely clear if oxDNA predictions may generally be relied upon as providing an accurate picture of ssNA molecular 3D conformation67. However, the oxDNA2 model incorporates a simple combination of both sequence-specific inter-base stacking strengths and a salt dependent non-bonded interaction term between backbone sites that could be sufficient to describe the molecular configuration of DNA. In fact we found that computed trends in |qcalc| values for oxDNA structures converted to all atom structures (see “Methods” for details) compared favourably with experimental |qeff| trends (i.e., qeff poly-dT > mixed sequence > -dA) determined with ETe for DNA, for both nb = 30 and nb = 60 (crosses in Fig. 4C). For example, the ratios in |qcalc| values of poly-dT to poly-dA, modelled with oxDNA, were calculated to be ≈ 1.03 and 1.06 for nb = 30 and nb = 60 respectively, which finds qualitative agreement with the ETe measured ratios of 1.08 and 1.11. At the same time, the mean end-to-end length, R≈10 nm, implied by the oxDNA results for poly-dT30 at the salt concentrations of our measurements (ca. 1 mM) is similar to DES-Amber (R≈ 9 nm) but is slightly larger than SAXS measurements, implying some quantitative disparities in comparisons of oxDNA results with other experimental data23,43.

We note however, that although successful for DNA, the oxRNA model, which incorporates the weakest inter-base stacking between adjacent U bases compared to A and T, generated structures with the largest values of Rg and R, and the largest |qcalc| values of all sequences considered, similar to the results obtained with the CHARMM36 forcefield (see Fig. 4C). This trend, reflected in the ratio of |qcalc| values for poly-dT to poly-rU of approximately 0.97, does not agree well with our measured ratio in |qeff| of 1.07 for nb= 60. Thus, of all the models tested for RNA it appears that the DES-Amber forcefield better captures the trend in effective charge measurements. This likely points to the importance of the incorporation of QM-level accurate base stacking and torsional energetics in determining RNA 3D conformation65.

In summary, ETe reveals a lower |qeff| value for poly-dA DNA compared to poly-dT, which we may attribute to a more compact polymer structure by comparing our measured trends with those that are reflected in MD model structures. At the microscopic level, this compaction is believed to arise from an increased inter-base stacking strength of neighbouring adenine bases compared to thymine (see Fig. 4B), further experimental evidence for which can be found from AFM pulling measurements68–70. Many of the other inter-base stacking strengths (including G, C and other pairwise permutations) have been calculated to lie in between these two extrema which can in fact be seen to be reflected in the intermediate value of both qeff and electrophoretic mobility of the mixed sequence case for nb = 60 (see Figs. 4C, 3C)71. However, it is interesting that the ETe measurements place poly-rU as having an effective charge similar to that of poly-dA, whereas arguments on base stacking strength alone would place poly-rU as having the highest |qeff| of all homopolymers examined in this study, as also suggested by coarse grained oxRNA simulations. However, poly-rU structures obtained with the DES-Amber forcefield clearly indicate a shorter average end-to-end distance R and contour length lc for RNA compared to DNA, and capture the ETe readout better.

The rigid-rod model of ssDNA

Although molecular simulations are growing increasingly sophisticated and powerful, field theoretical descriptions remain important in quantitative modelling of experimental data. This is particularly important in descriptions of many-body phenomena where computational expense grows rapidly with system size. For example, in a solution phase electrostatics simulation a typical atomistic computation could involve at least 1010 ions and up to 1000 times as many water molecules.

Thus in order to glean further physical insight into the experimental trends observed for qeff as a function of ssNA sequence we sought to relate our measured effective charge values to simple rigid rod models of the ssNAs within a continuum electrostatics framework. Although for dsDNA the rigid rod electrostatic model is an excellent approximation of the atomic level reality, its appropriateness for describing the interactions of single stranded nucleic acids is less obvious38,40,62. Specifically, given the high conformational flexibility of the disordered polyelectrolyte chain, it is not clear that the charged groups in the molecule could be reasonably expected to lie on the surface of a rectilinear rod as suggested in the simplest models of polyelectrolyte theory40,72,73. Interestingly however, a rigid rod model of ssDNA has been shown to be sufficient in PB modelling of the ion excess as measured by AES experiments62.

We calculated the effective charge of uniformly charged rigid rods of radii r = 0.05 and 0.4 nm using our continuum electrostatics calculation framework (see “Methods”). At the lower limit the rod radius approximates a line charge, whilst the upper radius limit may be thought to incorporate the excluded volume due to the finite size of the ssNA backbone atoms as well as that of a hydrated cation (Fig. 5A). We have previously shown that accounting for the excluded volume due to ions and water of hydration at the molecular interface in this manner permits us to incorporate the role of finite ion size in a point-ion based PB model34. Previous studies on the generation of molecular surfaces for PB calculations suggested that a value of w≈ 0.2 nm produces an ion accessible surface for an atomistic structure that successfully incorporates the role of the finite size of a hydrated ion34. Here, w is a parameter which represents a thickness to the molecular surface in addition to that represented by the vdW surface. Together with the radius of the backbone atoms of ≈ 0.2 nm we may therefore expect a rod of radius r≈ 0.4 nm to provide a reasonable description of a ssNA molecule within the rod model. We calculated the effective charge, qcalc, for rods of length l=nbb carrying a total charge of qstr=-nb+1e, where b is the axial charge spacing (Fig. 5A,C). We focused on rod lengths corresponding to nb=60 as this represents the most highly charge renormalizing regime in our experiments. An estimate of the renormalized charge of two ATTO dye molecules was added to the effective charge value calculated for a rod representing a label-free 60 base DNA in order to facilitate comparison with experiment (see “Methods” and Fig. S3). We determined qcalc values for rods of various axial spacings b ranging from 0.35 to 0.65 nm. Setting qcalc=qeff in the obtained qcalc vs b relationships then yielded an estimate of the value of b corresponding to each of the two values of rod radius for each molecular species (Fig. 5B).Fig. 5 The rigid rod model for a single stranded nucleic acid. (A) Schematic depiction of a flexible ssNA backbone contour consisting of phosphate sites (pink spheres) occupying the center of a tube of finite diameter that represents the volume excluded to the center of mass of an ion. A further thickness of w≈ 0.2 nm produces an ion accessible surface that successfully incorporates the role of the finite size of a hydrated cation. The equivalent rigid rod is depicted as a grey vertical cylinder with axial inter-charge spacing b. (B) Calculated values of effective charge |qcalc| vs axial inter-charge spacing b for the rigid rod models of radius r= 0.05 and 0.4 nm and length l=nbb. ETe measurements for nb= 60 shown as horizontal coloured bands for sequences poly-dT (blue), -dA (red), mixed DNA sequence (purple) and -rU (orange). Inferred b values are bT≈0.61-0.71 nm, and bA≈0.48-0.59 nm for poly-dT and -dA based on the present range of of r values (downward pointing arrows). Linear fit equations to the data corresponding to r= 0.05 and 0.4 nm in the ranges of interest are |qcalc|=44.2b+13.22 and |qcalc|=37.36b+21.25 respectively. (C) Electrostatic potential distribution obtained by solving the PB equation for a representative oxDNA poly-dT structure (nb= 60) converted into an all-atom model with a width parameter of w= 0.2 nm (left). The polymer is characterised by an end-to-end distance R=22.5 nm and contour length lc= 36 nm. The effective charge of this polymer is matched by an equivalent rod model with dimensions r = 0.4 nm and b = 0.5 nm (right). Electrostatic potential contours are shown at ϕ= − 5, − 2.5, − 1.25, − 0.625 and − 0.3125 kBT/e. The average electrical potential on the surfaces of the poly-dT structure and rod model are − 5.295 and − 6.58 kBT/e respectively.

At the upper limit of rod radius, r= 0.4 nm, we obtained inferred values of b= 0.61 nm for poly-dT and b≈ 0.5 nm for poly-dA and poly-rU (Fig. 5B). These values are comparable with the values of contour length per base, bc, obtained for poly-dT and poly-rU in SAXS experiments on 40 base ssNAs of about 0.6 and 0.5 nm respectively24. Specifically, SAXS measurements report bc= 0.56 nm and 0.49 nm for poly-dT and poly-rU respectively (see Table 1) and this trend is in good agreement with our inferences of axial charge spacing based on the rod model. However, we in fact expect a quantitative disparity between b and bc since the rigid rod model effectively replaces a curvilinear chain representing a flexible polyelectrolyte by a rod with a single rectilinear axis and a projected axial charge spacing, b (Fig. 5A). For a disordered single-stranded polynucleotide structure therefore, b is to be interpreted as the mean charge spacing along the rectilinear axis which is expected to be smaller than bc72–74 (Fig. 5A,C; see Supplementary Information Sect. S6 for further detail). The fact that our present results indicate b≈bc may in fact point to a current overestimate of measured qeff, as encountered and previously discussed in the comparison of the experimental qeff with qcalc from MD models.Table 1 Left: tabulated values of the contour length per base, bc, measured experimentally with techniques such as SAXS (*—Ref.24), AFM (†—Ref.98, ‡—Ref.99), transient electric birefringence (TEB, Ref.100), fluorescence correlation spectroscopy data in combination with a mean field theory (FCS, Ref.75) and from X-ray diffraction of ssDNA-protein complexes (XRD, Ref.26). ETe (and ASAXS in one case, Δ—Ref.43) measures a charge renormalization factor, η (bottom rows); right: η values for nb=60, the most highly charge renormalizing regime, can be mapped onto a value of the axial base spacing, b, using various polyelectrolyte models, such as the rigid rod model discussed in the main text, or by means of a theoretical model such as that in Ref.40 (pink shaded rows).

Calculated values of b inferred from studies of DNA melting data in combination with a charge renormalization model, such as in Refs.74,97, are also presented for comparison.

For the simple line charge description (r= 0.05 nm) in turn, we inferred b values of 0.71, 0.65, 0.63 and 0.59 nm for poly-dT, mixed, -rU & -dA sequences respectively (Fig. 5B). These axial charge spacings are higher than most estimates of bc (see Table 1), but do appear to find agreement with some experimental values obtained for the contour length per base for short ssNA oligomers in ssNA-protein complexes as measured by X-ray diffraction (XRD, Ref.26), and inferred for long ssNAs from fluorescence correlation spectroscopy data (FCS, Ref.75). Thus it appears that the rod model may not only provide a meaningful description of the electrostatics, but in doing so appears to be able to quantitatively capture an indication of the interphosphate spacing along the backbone contour, reflecting the essential impact of charge density on the electrostatics problem for ssNAs. This simplified view of the problem may therefore furnish a relatively straightforward modelling framework within which to interpret and parametrise differences in local structure as reflected in charge spacing in disordered biomolecules in solution.

Conclusion

We have measured effective charge values of single stranded nucleic acids of varying lengths and sequences in solution using molecular electrometry. We find that the magnitude of effective charge of poly-dT oligonucleotides is significantly higher than poly-dA and poly-rU sequences, despite all species carrying identical structural charge. These measurements indicate that conformational and structural differences between ssNAs of varying sequences are strongly reflected in their electrostatic properties. Since ETe measurements are performed under thermal and mechanical equilibrium in solution, and in the absence of external fields and sieving matrices, the measurements are amenable to computational modelling using both coarse grained and atomistic approaches. Comparing experimental data with calculated effective charge values determined using all-atom and coarse-grained simulations of ssNAs suggests that the lower measured effective charge values of poly-dA, as compared to poly-dT, arises from more compact molecular conformations, brought about arguably by enhanced inter-base stacking energies for adenine bases compared to thymine. For poly-rU oligonucleotides, structural modelling using a specialised ssRNA forcefield points towards a shorter average contour length than DNA and increased flexibility of the polymer chain, which in turn seem to underpin the consistently lower magnitude of both measured and calculated effective charge values for poly-rU compared to poly-dT. Our comparison of measured values of effective charge for ssNAs with those calculated for the same species using a variety of molecular models has highlighted some deficiencies of popular all-atom forcefields, thus recapitulating the widely recognised need for improved forcefield descriptions of nucleic acids63. Although our measurements of effective charge for the longer oligomers agree with the general magnitude of the corresponding computed values, we note at present, a consistent overestimate of 10–20% of measured values compared to computed quantities. Given that the present level of inaccuracy of ETe is expected to be around 10%, the present results would still suggest a residual disparity whose origin may become clearer with advances in modelling procedures or indeed further experimentation. We have further found that mapping our measured effective charge values for ssNAs on to that of a uniformly charged rod not only captures the essential impact of charge density on the electrostatic problem but also, importantly, yields estimates of a key structural parameter—the average contour length per base, bc—describing the molecule, which we find are in good agreement with measured values using SAXS and X-ray diffraction.

Finally, ETe offers significant operational advantages over state-of the-art structural methods such as AES and SAXS in that the measurement is performed using a standard wide-field fluorescence microscope on sub-zeptomole amounts of sample, e.g., 100 molecules yields a measurement precision on qeff smaller than 3%, and in a typical measurement duration of 1 min which can be significantly improved upon by measuring a larger number of molecules or for longer durations. Furthermore, the small measurement imprecision (< 1%) readily attainable in ETe permits small differences in charge renormalization between species to be measured with high certainty, which may not hold for ion counting experimental readouts (see Supplementary Information Sect. S7). Importantly, similar to FRET, the shorter measurement durations characteristic of single molecule optical detection experiments, and the access to individual molecular behaviour in real time, offers the prospect of detecting molecular conformation dynamics—albeit on the sub-second timescales—that is not readily accessible via other approaches. ETe therefore combines high measurement precision with several advantages of existing experimental methods directed at reading out molecular conformational properties. Finally, although the ETe measurements in this work were performed in low salt conditions (approximately 1 mM NaCl), further improvements to the methodology and device design should in future enable measurements at higher ionic strengths of 50–100 mM. As shown in Figs. 1 and S2, the depth of the potential well depends on the electrical potential at the midplane of the slit, ϕm=2ϕsexp(-κh), which suggests that reduction of the slit height to the range of h= 10–20 nm will support measurements at ionic strengths of around 50 mM. Divalent or multivalent cations may further be present in the electrolyte noting that they contribute more strongly to the total ionic strength than monovalent salts. Future improvements may also entail label-free imaging of the escape dynamics and molecular property measurements in the solution phase76.

ETe thus offers a facile and rapid measurement approach that opens up the dimension of effective charge as a means to test and validate outputs of theoretical and computational platforms for structural and conformational modelling of charged disordered biomolecules. The ability to detect subtle differences in 3D conformation of nucleic acids has immediate ramifications for the rapid, high throughput, coarse-grained characterization of biomolecular structures and conformation in the solution phase. This combined approach will likely foster structural and conformational inferences and modelling of more complex molecules such as folded RNAs for which structural studies lag significantly behind their protein counterparts77. Beyond examining the conformation of molecular species in the pure state, the approach could also provide a platform to screen and investigate molecular conformational changes induced by the binding of small-molecule ligands and metabolites to proteins and riboswitches, a process which has far reaching implications in gene expression, enzymatic activity and drug design78,79.

Methods

Experimental methods

Characterization of DNA and RNA samples

The short chain nucleic acid fragments and nb= 60 ssDNA fragments were purchased from Microsynth. nb= 30 dsDNA and nb= 60 ssRNA samples were purchased from IBA Lifesciences (Germany) and all nb= 30 ssNA samples were purchased from Integrated DNA Technologies (IDT). All fragments were doubly labelled with ATTO 532 dye molecules, apart from the nb= 5 ssDNA sample which was singly labelled (see Fig. S1 for a chemical diagram). The ssNA samples were examined with 20% polyacrylamide native gel electrophoresis (Novex™ TBE Gels, 20%, ThermoFisher Scientific) (Fig. 3C). Sequences for the short chain DNA fragments are as follows: nb= 12: TAG AAC TAG TGG, nb= 11: TAG AAC TAG TG, nb= 9: TAG AAC TAG, nb= 8: TAG AAC TA nb= 5: TAG AA.

Measuring the effective charge of single stranded nucleic acids (ssNAs) using Escape Time Electrometry (ETe)

Devices for ETe measurements were fabricated to the same specification and using procedures described previously32–34,36,37,53. Landscapes of electrostatic traps were created in nanoslits of height 2h= 70–75 nm and channel width = 5 μm, with nanostructure pocket depths d= 280 nm, and radius = 300 nm (see Fig. 1A). The height of the nanoslits were measured using atomic force microscopy. Nanoslits were loaded with a suspension of the molecular species at a concentration of 100–125 pM. Solution pH and conductivity was measured before and after the measurement using a micro-pH electrode (SevenDirectSD23, Mettler Toledo, UK) and microconductivity meter (Laquatwin, Horiba Scientific, Japan). ssNAs were loaded into an ETe device which consists of parallel nanoslits each containing arrays of electrostatic fluidic traps as shown in Fig. 1.

Wide-field fluorescence microscopy of the escape dynamics of molecules in the trap landscape was performed using excitation from a 532 nm DPSS laser (DMPV-532-3, Del Mar Photonics) that was focused at the back aperture of a 60 ×, NA = 1.2 water immersion objective (Olympus, UK). Time lapse videos were recorded at 10–20 Hz using an exposure time texp= 5 ms, for a measurement time of 3–5 min using a sCMOS camera (Kinetix Photometrix). Extensive details of the ETe measurement technique are provided in previous work, contained in Refs.32,36,37,53. As in previous work, the probability density function of residence times of molecules in the traps, PΔt∝1tescexp-Δttesc, was recorded and used to extract the average molecular escape time, tesc. In general we do not include in the fitting process 1–2 data points at the shortest measured lag times, Δt. This is because escape events recorded at the shortest durations can contain contributions from transiently trapped, weakly charged molecular species such as free dye, degraded material, or free ssDNA strands in the background solution in the case of the dsDNA calibrator molecules. The average escape time is related to the well depth W according to Kramer’s relation tesc∝expWkBT80.

In order to convert measurements of tesc to the depth of the trap, W, we performed Brownian dynamics (BD) simulations of the escape process. Full details of BD simulations are covered extensively in our previous works, particularly in Refs.32,36,37,53. In our BD simulations, we treat molecules as effective spheres of a radius equal to the measured hydrodynamic radius, rH, of the molecule, which was measured using fluorescence correlation spectroscopy. Observation of a large number of escape events (N≈104) allows W to be determined with a statistical uncertainty of about 1%37. Widths of presented Gaussian distributions of tesc therefore arise predominantly from the statistical error on the fitted tesc values which is approximately tesc/N. The total free energy W can be written as the sum of the electrostatic free energy of interaction ΔFel and an entropic term arising from the spatial fluctuations of the molecule ΔFtrans. A substantial contribution to W results from the electrostatic interaction free energy, ΔFel=Fslit-Fpocket, which is the difference in the axial minimum of electrostatic free energies for a molecule residing in the “slit” and “pocket” states, as shown in Figs. 1A and S232,36,37,53. ΔFel may also be calculated for a given molecular structure as described in a subsequent section.

We have previously shown that effective charge values, qeff can be deduced via the relation ΔFel=qeffΔϕmid35. Here, ϕmid is the electric potential at the midplane of the nanoslit in the experimental device and is given by ϕmid≈ϕm≈2ϕs exp(-κh) where ϕs is the effective surface potential at the slit walls, whereas in the pocket region, ϕmid is zero by design. In order to accurately measure the effective charge qeff for a molecule of interest we require knowledge of the surface potential ϕs of the ETe measurement device which is a priori unknown. Therefore, in order to determine ϕs, a calibration measurement was first performed using a molecule of well-known effective charge, namely a 30 bp or 60 bp dsDNA molecule carrying two attached ATTO-532 dyes (|qeff|= 28.2 and 45.6 e respectively, in this study). The effective charge values of the calibrator molecules were determined from previous work. Specifically, Ref.34 showed that the electrostatic properties of the double helix could be well described by a cylinder of radius r= 1.3 nm (including that of a hydrated Na+ cation) and a rise per basepair of b= 0.31 nm. The imprecisions in the measured values of tesc and qeff are ∼ 0.5 and ∼ 1.5% for the 30 and 60 bp calibrator molecules respectively, ∼ 0.5% for 30 base ssDNA, and ∼ 1% for 60 base ssDNA, which are suitable for capturing qualitative trends and differences arising from disparities in conformation between molecular species that lie in the range of ≈ 5–10% in this study.

Simulation methods

Simulation protocol

This section discusses in detail the methods and procedures used to carry out the MD simulations in this work. Simulation input files, trajectories, and code for the analysis of the simulations performed in this study are available on our Figshare repository: 10.6084/m9.figshare.25392787.

Initial nucleic acid structures were generated with the Nuleic Acid Builder (NAB) as part of the AmberTools software81. These were then solvated via the CHARMM-GUI solution builder platform in a truncated octahedral box, leaving at least 1 nm of space between the molecule and the edge of the box82. The net charge of the molecule was neutralised with the corresponding number of monovalent Na+ counterions, and no additional background electrolyte concentration added so as to best approximate the added salt concentration of ≈ 1 mM NaCl in the experiments.

MD simulations were performed with GROMACS software83. The Lennard–Jones and short-range electrostatic interactions were calculated with cutoff values of 1.2 nm and 0.9 nm for simulations employing the CHARMM36 and the AMBER or AMBER variant (DES-Amber/HB-CUFIX) forcefields respectively. The particle mesh Ewald (PME) method was used to evaluate the long-range electrostatic interactions with 3D periodic boundaries and employed a 0.12 nm grid spacing. Constraints were applied with the LINCS algorithm to all bonds involving hydrogen atoms in the system. After energy minimization with the steepest descents algorithm, a short 1 ns equilibration was performed in the NVT ensemble at 100 K using a v-rescale thermostat and a timestep of 2 fs, followed by an equilibration in the NPT ensemble (1 atm, 300 K, 2 fs timestep) using a Berendsen barostat with a coupling constant of 1 ps—all steps with position restraints of 500 kJ mol–1 nm–2 on the heavy atoms of the molecule. In another 20 ns NPT equilibration (1 atm, 300 K, 2 fs timestep) using the Parrinello-Rahman barostat, these restraints were then gradually released before production simulations of at least 1 μs duration (see Table 2). For those models which displayed greater exploration of conformational parameter space we ran longer simulations of up to 3 μs (see Table 2). The first 100 ns of the production simulation trajectories were discarded in the subsequent analyses. We point out that although the required timescales needed for the full simulation convergence of flexible ssNAs may indeed be much larger than those considered here, we found that extending simulations beyond 1 μs did not appreciably alter the conformational (Rg,R) landscapes and therefore the results of this study.Table 2 Summary of the MD simulations performed in this work.

Force field + water model	System	Duration (μs)	
Amberbsc1 + TIP3P	dT30, dA30	1	
Charmm36 + TIP3P	dT30	2	
Charmm36 + TIP3P	dA30, rU30	1	
DES-AmberSF1.0 + TIP4D	dT30, dA30	1	
DES-Amber_RNA + TIP4D	rU30	3	
HB-CUFIX + TIP3P	rU30	2	
DES-AmberSF1.0 + TIP4D	nb= 5–12 (short fragments)	0.5	

oxDNA simulations

We simulated ssNAs at 298 K, in 1 mM salt, using oxDNA version 2 following a protocol outlined in Ref.84. We performed two independent VMMC simulations for each system with 2 × 107 attempted steps per particle. The first 105 moves were discarded in the subsequent analysis. Sequence dependent stacking strengths and hydrogen bond parameters from oxDNA2 were implemented. Representative structures taken from oxDNA trajectories, chosen using the same method as described below for the all-atom MD trajectories, were then converted to all atom PDB format with tacoxDNA and a short energy minimization was performed with the amberbsc1 forcefield85.

Poisson Boltzmann (PB) free energy calculations to calculate theoretical molecular effective charge, qcalc

The calculation of qeff has robust theoretical underpinnings in Poisson-Boltzmann theory for solution phase electrostatics which provides a satisfactory theoretical description of the experiments in the regime of low monovalent salt concentrations35,38,54. Model structures of ssNAs used in PB calculations were obtained from MD simulations as described in the main text. Atomic coordinates of the molecule in PDB format were converted into PQR format (containing atomic charge information) with the MDAnalysis package86. However, for all atomic forcefields considered in this study, a polymer chain with number of bases nb is characterised by a structural charge of qstr=-nb-1 e, due to the non-integer charges of the terminal residues (which for instance contribute ≈ − 0.7 and − 0.3 e respectively in the Amber formalism, rather than a total charge of − 1 e characteristic of the interior bases). Prior to performing a PB calculation on the model structure, we therefore modified the atomic partial charges of the hydrogen termini atoms (named H5T, H3T) by − 1 e each to give a model with charge -(nb+1) e. This slight modification ensures that the model structures of the polymer chain have the same structural charge as the ‘nucleic acid core’ of the molecules in experiments (see Fig. S1 for a chemical diagram). In practice the molecules also carry covalently linked dye molecules that contribute to both the molecule’s total structural charge (qstr =-(nb+3) e) as well its effective charge. The renormalized charge of the effective of the attached ATTO dye molecules, qeff,dye, (the value of which was itself found to depend on nb, see Fig. S3) was added to the calculated molecular effective charge in order to facilitate comparison of qcalc with the experimentally measured qeff values. Atomic charges for the model of an ATTO dye molecule were taken from the AMBER-DYES library87. Note that these calculated values for qcalc,dye= − 0.62 e (for nb = 60) are in excellent agreement with the previously measured values under the same experimental conditions34.

In order to calculate effective charge values for a given molecular structure we consider molecular representations based on an ion accessible surface (IAS) deduced from atomistic structures. Previous work inferring structural parameters for dsDNA and dsRNA using measured effective charge values from ETe in combination with PB molecular electrostatic models pointed towards the presence of a thin region around the nucleic acids that is inaccessible to point ions in solution. The thickness of this layer, w, was found to correlate with values for the hydrated cation radii34,88. In previous work, inclusion of an additional thickness in addition to the van der Waals (vdW) surface, 0.14≤w≤0.3 nm, permitted measured electrostatic free energies to be captured using a PB view based on point-ions in solution that does not explicitly consider finite ion size. Previous work suggested a value of w≈ 0.2 nm could capture the IAS for experiments with Na+ ions in solution88. Therefore, for the molecular electrostatic models considered in this study, we assigned each atom (excluding hydrogen atoms) a radius equal to the vdW radius (taken from the MD forcefield), and added a further 0.2 nm to generate an effective ion accessible surface. We used the Nanoshaper package in order to generate molecular surfaces of atomic coordinates using a probe size of 0.1 nm38,89. These surfaces were subsequently processed with the GMSH package to generate the required 3D mesh for our finite-element calculation90.

We then proceeded to calculate theoretical qeff values of model ssNA structures in our PB interaction free energy framework, as previously outlined extensively in Refs.33–35,38. The governing equations were solved for the cylinder models and model molecular surfaces using the COMSOL Multiphysics and FEniCS finite-element software packages respectively91–93. We solved the non-linear PB equation for the outer region of the molecular surface (i.e. in the electrolyte) for both the ‘slit’ and ‘pocket’ states in the experimental geometry, as shown in Figs. 1A and S2. The principal axis of the molecule was aligned parallel to the slit in the xy plane (i.e., this configuration represents the local minimum in free energy). In the electrolyte region, the electric potential, ϕ, is given by the solution of the PB equation: -∇·ϵϵ0∇ϕ=ρe , where ϵϵ0 is the permittivity of the medium and, ρe=ΣiciNAeziexp(-zieϕkBT) is the ion density in solution. Here ci denotes the bulk concentration of ions of species i and is simply the salt concentration, c0, in solution for monovalent salts used in these experiments.

We uniformly distributed the partial charges of each atom over spherical seeds of radius 0.2 Å which is substantially smaller than typical values of the vdW radii of atoms. It should be noted that although the absolute electrostatic energy of the interior depends on the radius of the spherical seeds, which is an arbitrary parameter, the free energy difference ,ΔF, is independent of this parameter due to a cancellation effect which arises from subtraction of electrostatic free energies in the pocket and slit, ΔFel=Fslit-Fpocket (Fig. S2). In practice, the contribution of the interior region of the atoms to the free energy difference and hence the effective charge is negligible38. Therefore, in our free energy calculation, we only consider the contributions from the solvent region. Accordingly, the free energy functional is given by Fel=∫V[ϵϵ02E.E-2c0NAkBT(coshψ-ψsinhψ-1)]dV where V is the volume of the solvent region, E=-∇ϕ is the electric field and ψ=eϕ/kBT is the non-dimensional electric potential54,94.

For cylinders the PB equation was solved using Neumann boundary conditions that specify constant charge densities at the surface of the cylinder models. A constant charge approximation of ssNAs holds due to the acidic nature of the phosphate groups under the conditions of the ETe experiment. The silica walls were also treated with constant charge boundary conditions, as discussed previously34. For the silica surfaces, we have -ϵ0ϵ∇ϕ.n=σw, where σw=-0.1e/nm2 is a nominal surface charge density and n is the normal vector pointing into the electrolyte.

The following matching boundary conditions were also satisfied at the molecular surface: ϵ0ϵmol∇ϕmol.n=ϵ0ϵ∇ϕ.n and n×∇ϕ-∇ϕmol=0. Here, n is the normal vector pointing into the electrolyte and the above matching conditions imply the continuity of the normal component of the displacement field, D=-ϵ0ϵ∇ϕ, and the tangential component of E across the interface, respectively. We then deduced the effective charge of the molecule from the relationship ΔFel=qcalcϕm where ΔFel=Fslit-Fpocket and ϕm is the electrical potential at the midplane of the slit in the absence of the molecule35,54.

Electrostatic visualization of MD structures

Solving the PB equation determines the electric potential on a 3D finite element mesh inside and outside of the molecular surface which contains the partial charges of atoms. We first labelled each point in the mesh according to whether it resided in either the molecular volume or the electrolyte region. The molecular volume was defined by generating a molecular surface from the atomic PQR file using the Nanoshaper package using a probe size of 0.1 nm as mentioned previously38,89. The electrical potential on the molecular surface can be obtained by interpolation of data from the 3D mesh onto the 2D surface. The process of interpolation and visualization of the data can be efficiently handled by Paraview which provides a wide range of filters for processing large datasets95.

Supplementary Information

Supplementary Information.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-70641-x.

Acknowledgements

The authors thank Stefano Piana-Agostinetti from the David. E. Shaw research group, and Jonathan Doye for discussions. The authors acknowledge the use of the University of Oxford Advanced Research Computing (ARC) facility in carrying out this work. Molecular graphics were generated with Protein Imager96.

Author contributions

RWG and XZ contributed equally to this work. XZ and TJDB performed experiments and analyzed data. RWG performed molecular simulations and molecular modelling. RWG and AB performed calculations. MK designed and supervised the research. RWG and MK wrote the manuscript.

Funding

This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (Grant Agreement No. 724180) and the UKRI Biotechnology and Biological Sciences Research Council (BBSRC) (award reference BB/W017415/1).

Data availability

The datasets generated during and/or analysed during the current study are available in the ‘ssNA electrometry and molecular simulations’ repository at 10.6084/m9.figshare.25392787.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

These authors contributed equally: Rowan Walker-Gibbons and Xin Zhu.
==== Refs
References

1. Watson JD Crick FH Molecular structure of nucleic acids; a structure for deoxyribose nucleic acid Nature 1953 171 737 738 10.1038/171737a0 13054692
Watson, J. D. & Crick, F. H. Molecular structure of nucleic acids; a structure for deoxyribose nucleic acid. Nature 171, 737–738 (1953).13054692 10.1038/171737a0
2. Dickerson RE The anatomy of A-, B-, and Z-DNA Science 1982 216 475 485 10.1126/science.7071593 7071593
Dickerson, R. E. et al. The anatomy of A-, B-, and Z-DNA. Science 216, 475–485 (1982).7071593 10.1126/science.7071593
3. Tjandra N Tate S-I Ono A Kainosho M Bax A The NMR structure of a DNA Dodecamer in an aqueous dilute liquid crystalline phase J. Am. Chem. Soc. 2000 122 6190 6200 10.1021/ja000324n
Tjandra, N., Tate, S.-I., Ono, A., Kainosho, M. & Bax, A. The NMR structure of a DNA Dodecamer in an aqueous dilute liquid crystalline phase. J. Am. Chem. Soc. 122, 6190–6200 (2000).10.1021/ja000324n
4. Steve, S. & Pollack, L. Determining the locations of ions and water around DNA from X-Ray scattering measurements. Biophys. J. 108, 2886–2895 (2015).
5. Singh AK Wen C Cheng S Vinh NQ Long-range DNA-water interactions Biophys. J. 2021 120 4966 4979 10.1016/j.bpj.2021.10.016 34687717
Singh, A. K., Wen, C., Cheng, S. & Vinh, N. Q. Long-range DNA-water interactions. Biophys. J. 120, 4966–4979 (2021).34687717 10.1016/j.bpj.2021.10.016
6. Cozzolino F Iacobucci I Monaco V Monti M Protein–DNA/RNA interactions: An overview of investigation methods in the -omics era J. Proteome Res. 2021 20 3018 3030 10.1021/acs.jproteome.1c00074 33961438
Cozzolino, F., Iacobucci, I., Monaco, V. & Monti, M. Protein–DNA/RNA interactions: An overview of investigation methods in the -omics era. J. Proteome Res. 20, 3018–3030 (2021).33961438 10.1021/acs.jproteome.1c00074
7. Eliezer D Biophysical characterization of intrinsically disordered proteins Curr. Opin. Struct. Biol. 2009 19 23 30 10.1016/j.sbi.2008.12.004 19162471
Eliezer, D. Biophysical characterization of intrinsically disordered proteins. Curr. Opin. Struct. Biol. 19, 23–30 (2009).19162471 10.1016/j.sbi.2008.12.004
8. Evans R Ramisetty S Kulkarni P Weninger K Illuminating intrinsically disordered proteins with integrative structural biology Biomolecules 2023 13 124 10.3390/biom13010124 36671509
Evans, R., Ramisetty, S., Kulkarni, P. & Weninger, K. Illuminating intrinsically disordered proteins with integrative structural biology. Biomolecules 13, 124 (2023).36671509 10.3390/biom13010124
9. Strutt, S. C., Torrez, R. M., Kaya, E., Negrete, O. A. & Doudna, J. A. RNA-dependent RNA targeting by CRISPR-Cas9. eLife 7 (2018).
10. Knott GJ Doudna JA CRISPR-Cas guides the future of genetic engineering Science 2018 361 866 869 10.1126/science.aat5011 30166482
Knott, G. J. & Doudna, J. A. CRISPR-Cas guides the future of genetic engineering. Science 361, 866–869 (2018).30166482 10.1126/science.aat5011
11. Maharana S RNA buffers the phase separation behavior of prion-like RNA binding proteins Science 2018 360 918 921 10.1126/science.aar7366 29650702
Maharana, S. et al. RNA buffers the phase separation behavior of prion-like RNA binding proteins. Science 360, 918–921 (2018).29650702 10.1126/science.aar7366
12. Tom JKA Onuchic PL Deniz AA Short PolyA RNA homopolymers Undergo Mg2+-mediated kinetically arrested condensation J. Phys. Chem. B 2022 126 9715 9725 10.1021/acs.jpcb.2c05935 36378781
Tom, J. K. A., Onuchic, P. L. & Deniz, A. A. Short PolyA RNA homopolymers Undergo Mg2+-mediated kinetically arrested condensation. J. Phys. Chem. B 126, 9715–9725 (2022).36378781 10.1021/acs.jpcb.2c05935
13. Cieniková Z Damberger FF Hall J Allain FHT Maris C Structural and mechanistic insights into Poly(uridine) tract recognition by the hnRNP C RNA recognition motif J. Am. Chem. Soc. 2014 136 14536 14544 10.1021/ja507690d 25216038
Cieniková, Z., Damberger, F. F., Hall, J., Allain, F. H. T. & Maris, C. Structural and mechanistic insights into Poly(uridine) tract recognition by the hnRNP C RNA recognition motif. J. Am. Chem. Soc. 136, 14536–14544 (2014).25216038 10.1021/ja507690d
14. Thurm AR Beren C Duran-Meza AL Knobler CM Gelbart WM RNA homopolymers form higher-curvature virus-like particles than do normal-composition RNAs Biophys. J. 2019 117 1331 1341 10.1016/j.bpj.2019.08.012 31514968
Thurm, A. R., Beren, C., Duran-Meza, A. L., Knobler, C. M. & Gelbart, W. M. RNA homopolymers form higher-curvature virus-like particles than do normal-composition RNAs. Biophys. J. 117, 1331–1341 (2019).31514968 10.1016/j.bpj.2019.08.012
15. Plumridge A Andresen K Pollack L Visualizing disordered single-stranded RNA: Connecting sequence, structure, and electrostatics J. Am. Chem. Soc. 2020 142 109 119 10.1021/jacs.9b04461 31804813
Plumridge, A., Andresen, K. & Pollack, L. Visualizing disordered single-stranded RNA: Connecting sequence, structure, and electrostatics. J. Am. Chem. Soc. 142, 109–119 (2020).31804813 10.1021/jacs.9b04461
16. Statello L Guo C-J Chen L-L Huarte M Gene regulation by long non-coding RNAs and its biological functions Nat. Rev. Mol. Cell Biol. 2021 22 96 118 10.1038/s41580-020-00315-9 33353982
Statello, L., Guo, C.-J., Chen, L.-L. & Huarte, M. Gene regulation by long non-coding RNAs and its biological functions. Nat. Rev. Mol. Cell Biol. 22, 96–118 (2021).33353982 10.1038/s41580-020-00315-9
17. Passmore LA Coller J Roles of mRNA poly(A) tails in regulation of eukaryotic gene expression Nat. Rev. Mol. Cell Biol. 2022 23 93 106 10.1038/s41580-021-00417-y 34594027
Passmore, L. A. & Coller, J. Roles of mRNA poly(A) tails in regulation of eukaryotic gene expression. Nat. Rev. Mol. Cell Biol. 23, 93–106 (2022).34594027 10.1038/s41580-021-00417-y
18. Oksuz O Transcription factors interact with RNA to regulate genes Mol. Cell. 2023 83 2449 2463 10.1016/j.molcel.2023.06.012 37402367
Oksuz, O. et al. Transcription factors interact with RNA to regulate genes. Mol. Cell. 83, 2449–2463 (2023).37402367 10.1016/j.molcel.2023.06.012
19. Han D Single-stranded DNA and RNA origami Science 2017 358 eaao2648 10.1126/science.aao2648 29242318
Han, D. et al. Single-stranded DNA and RNA origami. Science 358, eaao2648 (2017).29242318 10.1126/science.aao2648
20. Rinaldi C Wood MJA Antisense oligonucleotides: The next frontier for treatment of neurological disorders Nat. Rev. Neurol. 2018 14 9 21 10.1038/nrneurol.2017.148 29192260
Rinaldi, C. & Wood, M. J. A. Antisense oligonucleotides: The next frontier for treatment of neurological disorders. Nat. Rev. Neurol. 14, 9–21 (2018).29192260 10.1038/nrneurol.2017.148
21. Oweida TJ Kim HS Donald JM Singh A Yingling YG Assessment of AMBER force fields for simulations of ssDNA J. Chem. Theory Comput. 2021 17 1208 1217 10.1021/acs.jctc.0c00931 33434436
Oweida, T. J., Kim, H. S., Donald, J. M., Singh, A. & Yingling, Y. G. Assessment of AMBER force fields for simulations of ssDNA. J. Chem. Theory Comput. 17, 1208–1217 (2021).33434436 10.1021/acs.jctc.0c00931
22. Scharner J Aznarez I Clinical applications of single-stranded oligonucleotides: Current landscape of approved and in-development therapeutics Mol. Ther. 2021 29 540 554 10.1016/j.ymthe.2020.12.022 33359792
Scharner, J. & Aznarez, I. Clinical applications of single-stranded oligonucleotides: Current landscape of approved and in-development therapeutics. Mol. Ther. 29, 540–554 (2021).33359792 10.1016/j.ymthe.2020.12.022
23. Plumridge A Meisburger SP Pollack L Visualizing single-stranded nucleic acids in solution Nucleic Acids Res. 2017 45 e66 e66 28034955
Plumridge, A., Meisburger, S. P. & Pollack, L. Visualizing single-stranded nucleic acids in solution. Nucleic Acids Res. 45, e66–e66 (2017).28034955
24. Chen H Ionic strength-dependent persistence lengths of single-stranded RNA and DNA Proc. Natl. Acad. Sci. 2012 109 799 804 10.1073/pnas.1119057109 22203973
Chen, H. et al. Ionic strength-dependent persistence lengths of single-stranded RNA and DNA. Proc. Natl. Acad. Sci. 109, 799–804 (2012).22203973 10.1073/pnas.1119057109
25. Tanikawa, J., Nishimura, Y., Hirao, I. & Miura, K. NMR spectroscopic study of single-stranded DNA fragments of d(CGGCGAAAGCCG) and d(CGGCAAAAGCCG). Nucleic Acids Symp. Ser. 47–48 (1991).
26. Murphy MC Rasnik I Cheng W Lohman TM Ha T Probing single-stranded DNA conformational flexibility using fluorescence spectroscopy Biophys. J. 2004 86 2530 2537 10.1016/S0006-3495(04)74308-8 15041689
Murphy, M. C., Rasnik, I., Cheng, W., Lohman, T. M. & Ha, T. Probing single-stranded DNA conformational flexibility using fluorescence spectroscopy. Biophys. J. 86, 2530–2537 (2004).15041689 10.1016/S0006-3495(04)74308-8
27. Isaksson J Acharya S Barman J Cheruku P Chattopadhyaya J Single-stranded adenine-rich DNA and RNA retain structural characteristics of their respective double-stranded conformations and show directional differences in stacking pattern Biochemistry 2004 43 15996 16010 10.1021/bi048221v 15609994
Isaksson, J., Acharya, S., Barman, J., Cheruku, P. & Chattopadhyaya, J. Single-stranded adenine-rich DNA and RNA retain structural characteristics of their respective double-stranded conformations and show directional differences in stacking pattern. Biochemistry 43, 15996–16010 (2004).15609994 10.1021/bi048221v
28. Sim AYL Lipfert J Herschlag D Doniach S Salt dependence of the radius of gyration and flexibility of single-stranded DNA in solution probed by small-angle x-ray scattering Phys. Rev. E 2012 86 1 10.1103/PhysRevE.86.021901
Sim, A. Y. L., Lipfert, J., Herschlag, D. & Doniach, S. Salt dependence of the radius of gyration and flexibility of single-stranded DNA in solution probed by small-angle x-ray scattering. Phys. Rev. E 86, 1 (2012).10.1103/PhysRevE.86.021901
29. Tinland B Pluen A Sturm J Weill G Persistence length of single-stranded DNA Macromolecules 1997 30 5763 5765 10.1021/ma970381+
Tinland, B., Pluen, A., Sturm, J. & Weill, G. Persistence length of single-stranded DNA. Macromolecules 30, 5763–5765 (1997).10.1021/ma970381+
30. Laurence TA Kong X Jäger M Weiss S Probing structural heterogeneities and fluctuations of nucleic acids and denatured proteins Proc. Natl. Acad. Sci. 2005 102 17348 17353 10.1073/pnas.0508584102 16287971
Laurence, T. A., Kong, X., Jäger, M. & Weiss, S. Probing structural heterogeneities and fluctuations of nucleic acids and denatured proteins. Proc. Natl. Acad. Sci. 102, 17348–17353 (2005).16287971 10.1073/pnas.0508584102
31. Stellwagen E Stellwagen NC Electrophoretic mobility of DNA in solutions of high ionic strength Biophys. J. 2020 118 2783 2789 10.1016/j.bpj.2020.02.034 32445623
Stellwagen, E. & Stellwagen, N. C. Electrophoretic mobility of DNA in solutions of high ionic strength. Biophys. J. 118, 2783–2789 (2020).32445623 10.1016/j.bpj.2020.02.034
32. Ruggeri F Single-molecule electrometry Nature Nanotechnol. 2017 12 488 495 10.1038/nnano.2017.26 28288117
Ruggeri, F. et al. Single-molecule electrometry. Nature Nanotechnol. 12, 488–495 (2017).28288117 10.1038/nnano.2017.26
33. Kloes G Far-field electrostatic signatures of macromolecular 3D conformation Nano Lett. 2022 22 7834 7840 10.1021/acs.nanolett.2c02485 36125326
Kloes, G. et al. Far-field electrostatic signatures of macromolecular 3D conformation. Nano Lett. 22, 7834–7840 (2022).36125326 10.1021/acs.nanolett.2c02485
34. Bespalova M Behjatian A Karedla N Walker-Gibbons R Krishnan M Opto-electrostatic determination of nucleic acid double-helix dimensions and the structure of the molecule-solvent interface Macromolecules 2022 55 6200 6210 10.1021/acs.macromol.2c00657 35910310
Bespalova, M., Behjatian, A., Karedla, N., Walker-Gibbons, R. & Krishnan, M. Opto-electrostatic determination of nucleic acid double-helix dimensions and the structure of the molecule-solvent interface. Macromolecules 55, 6200–6210 (2022).35910310 10.1021/acs.macromol.2c00657
35. Krishnan M A simple model for electrical charge in globular macromolecules and linear polyelectrolytes in solution J. Chem. Phys. 2017 146 205101 10.1063/1.4983485 28571334
Krishnan, M. A simple model for electrical charge in globular macromolecules and linear polyelectrolytes in solution. J. Chem. Phys. 146, 205101 (2017).28571334 10.1063/1.4983485
36. Ruggeri F Krishnan M Spectrally resolved single-molecule electrometry J. Chem. Phys. 2017 148 123307 10.1063/1.5008936
Ruggeri, F. & Krishnan, M. Spectrally resolved single-molecule electrometry. J. Chem. Phys. 148, 123307 (2017).10.1063/1.5008936
37. Ruggeri F Krishnan M Entropic trapping of a singly charged molecule in solution Nano Lett. 2018 18 3773 3779 10.1021/acs.nanolett.8b01011 29688720
Ruggeri, F. & Krishnan, M. Entropic trapping of a singly charged molecule in solution. Nano Lett. 18, 3773–3779 (2018).29688720 10.1021/acs.nanolett.8b01011
38. Behjatian A Krishnan M Electrostatic free energies carry structural information on nucleic acid molecules in solution J. Chem. Phys. 2022 156 134201 10.1063/5.0080008 35395894
Behjatian, A. & Krishnan, M. Electrostatic free energies carry structural information on nucleic acid molecules in solution. J. Chem. Phys. 156, 134201 (2022).35395894 10.1063/5.0080008
39. Manning GS Electrostatic free energies of spheres, cylinders, and planes in counterion condensation theory with some applications Macromolecules 2007 40 8071 8081 10.1021/ma071457x
Manning, G. S. Electrostatic free energies of spheres, cylinders, and planes in counterion condensation theory with some applications. Macromolecules 40, 8071–8081 (2007).10.1021/ma071457x
40. Netz RR Orland H Variational charge renormalization in charged systems Eur. Phys. J. E 2003 11 301 311 10.1140/epje/i2002-10159-0 15011050
Netz, R. R. & Orland, H. Variational charge renormalization in charged systems. Eur. Phys. J. E 11, 301–311 (2003).15011050 10.1140/epje/i2002-10159-0
41. Manning GS Limiting laws and counterion condensation in polyelectrolyte solutions I. Colligative properties J. Chem. Phys. 2003 51 924 933 10.1063/1.1672157
Manning, G. S. Limiting laws and counterion condensation in polyelectrolyte solutions I. Colligative properties. J. Chem. Phys. 51, 924–933 (2003).10.1063/1.1672157
42. Aubouy M Trizac E Bocquet LR Effective charge versus bare charge: An analytical estimate for colloids in the infinite dilution limit J. Phys. A Math. General 2003 36 5835 5840 10.1088/0305-4470/36/22/302
Aubouy, M., Trizac, E. & Bocquet, L. R. Effective charge versus bare charge: An analytical estimate for colloids in the infinite dilution limit. J. Phys. A Math. General 36, 5835–5840 (2003).10.1088/0305-4470/36/22/302
43. Meisburger SP Polyelectrolyte properties of single stranded DNA measured using SAXS and single molecule FRET: Beyond the wormlike chain model Biopolymers 2013 99 1032 1045 10.1002/bip.22265 23606337
Meisburger, S. P. et al. Polyelectrolyte properties of single stranded DNA measured using SAXS and single molecule FRET: Beyond the wormlike chain model. Biopolymers 99, 1032–1045 (2013).23606337 10.1002/bip.22265
44. Lamm G Wong L Pack GR Monte Carlo and Poisson-Boltzmann calculations of the fraction of counterions bound to DNA Biopolymers 1994 34 227 237 10.1002/bip.360340209 8142591
Lamm, G., Wong, L. & Pack, G. R. Monte Carlo and Poisson-Boltzmann calculations of the fraction of counterions bound to DNA. Biopolymers 34, 227–237 (1994).8142591 10.1002/bip.360340209
45. Misra VK Honig B On the magnitude of the electrostatic contribution to ligand-DNA interactions Proc. Natl. Acad. Sci. 1995 92 4691 4695 10.1073/pnas.92.10.4691 7753866
Misra, V. K. & Honig, B. On the magnitude of the electrostatic contribution to ligand-DNA interactions. Proc. Natl. Acad. Sci. 92, 4691–4695 (1995).7753866 10.1073/pnas.92.10.4691
46. Fogolari F Brigo A Molinari H The Poisson–Boltzmann equation for biomolecular electrostatics: A tool for structural biology J. Mol. Recogn. 2002 15 377 392 10.1002/jmr.577
Fogolari, F., Brigo, A. & Molinari, H. The Poisson–Boltzmann equation for biomolecular electrostatics: A tool for structural biology. J. Mol. Recogn. 15, 377–392 (2002).10.1002/jmr.577
47. Fleck C Netz RR Von Grünberg HH Poisson–Boltzmann theory for membranes with mobile charged lipids and the pH-dependent interaction of a DNA molecule with a membrane Biophys. J. 2002 82 76 92 10.1016/S0006-3495(02)75375-7 11751297
Fleck, C., Netz, R. R. & Von Grünberg, H. H. Poisson–Boltzmann theory for membranes with mobile charged lipids and the pH-dependent interaction of a DNA molecule with a membrane. Biophys. J. 82, 76–92 (2002).11751297 10.1016/S0006-3495(02)75375-7
48. Kirmizialtin S Alexander E Marcia R The ionic atmosphere around A-RNA: Poisson–Boltzmann and molecular dynamics simulations Biophys. J. 2012 102 829 838 10.1016/j.bpj.2011.12.055 22385854
Kirmizialtin, S., Alexander, E. & Marcia, R. The ionic atmosphere around A-RNA: Poisson–Boltzmann and molecular dynamics simulations. Biophys. J. 102, 829–838 (2012).22385854 10.1016/j.bpj.2011.12.055
49. Bertonati C Honig B Alexov E Poisson–Boltzmann calculations of nonspecific salt effects on protein-protein binding free energies Biophys. J. 2007 92 1891 1899 10.1529/biophysj.106.092122 17208980
Bertonati, C., Honig, B. & Alexov, E. Poisson–Boltzmann calculations of nonspecific salt effects on protein-protein binding free energies. Biophys. J. 92, 1891–1899 (2007).17208980 10.1529/biophysj.106.092122
50. Chu VB Bai Y Lipfert J Herschlag D Doniach S Evaluation of ion binding to dna duplexes using a size-modified Poisson–Boltzmann theory Biophys. J. 2007 93 3202 3209 10.1529/biophysj.106.099168 17604318
Chu, V. B., Bai, Y., Lipfert, J., Herschlag, D. & Doniach, S. Evaluation of ion binding to dna duplexes using a size-modified Poisson–Boltzmann theory. Biophys. J. 93, 3202–3209 (2007).17604318 10.1529/biophysj.106.099168
51. Maffeo C DNA–DNA Interactions in tight supercoils are described by a small effective charge density Phys. Rev. Lett. 2010 105 158101 10.1103/PhysRevLett.105.158101 21230940
Maffeo, C. et al. DNA–DNA Interactions in tight supercoils are described by a small effective charge density. Phys. Rev. Lett. 105, 158101 (2010).21230940 10.1103/PhysRevLett.105.158101
52. Batys P Luukkonen S Sammalkorpi M Ability of the Poisson-Boltzmann equation to capture molecular dynamics predicted ion distribution around polyelectrolytes Phys. Chem. Chem. Phys. 2017 19 24583 24593 10.1039/C7CP02547E 28853454
Batys, P., Luukkonen, S. & Sammalkorpi, M. Ability of the Poisson-Boltzmann equation to capture molecular dynamics predicted ion distribution around polyelectrolytes. Phys. Chem. Chem. Phys. 19, 24583–24593 (2017).28853454 10.1039/C7CP02547E
53. Ruggeri F Krishnan M Lattice diffusion of a single molecule in solution Phys. Rev. E 2017 96 062406 10.1103/PhysRevE.96.062406 29347432
Ruggeri, F. & Krishnan, M. Lattice diffusion of a single molecule in solution. Phys. Rev. E 96, 062406 (2017).29347432 10.1103/PhysRevE.96.062406
54. Krishnan M Electrostatic free energy for a confined nanoscale object in a fluid J. Chem. Phys. 2013 138 114906 10.1063/1.4795087 23534661
Krishnan, M. Electrostatic free energy for a confined nanoscale object in a fluid. J. Chem. Phys. 138, 114906 (2013).23534661 10.1063/1.4795087
55. Kjellander R Mitchell DJ Dressed ion theory for electric double layer structure and interactions; an exact analysis Mol. Phys. 1997 91 173 188 10.1080/002689797171490
Kjellander, R. & Mitchell, D. J. Dressed ion theory for electric double layer structure and interactions; an exact analysis. Mol. Phys. 91, 173–188 (1997).10.1080/002689797171490
56. Ulander J Greberg H Kjellander R Primary and secondary effective charges for electrical double layer systems with asymmetric electrolytes J. Chem. Phys. 2001 115 7144 7160 10.1063/1.1398587
Ulander, J., Greberg, H. & Kjellander, R. Primary and secondary effective charges for electrical double layer systems with asymmetric electrolytes. J. Chem. Phys. 115, 7144–7160 (2001).10.1063/1.1398587
57. Yu B Iwahara J Experimental approaches for investigating ion atmospheres around nucleic acids and proteins Comput. Struct. Biotechnol. J. 2021 19 2279 2285 10.1016/j.csbj.2021.04.033 33995919
Yu, B. & Iwahara, J. Experimental approaches for investigating ion atmospheres around nucleic acids and proteins. Comput. Struct. Biotechnol. J. 19, 2279–2285 (2021).33995919 10.1016/j.csbj.2021.04.033
58. Lipfert J Doniach S Das R Herschlag D Understanding nucleic acid-ion interactions Annu. Rev. Biochem. 2014 83 813 841 10.1146/annurev-biochem-060409-092720 24606136
Lipfert, J., Doniach, S., Das, R. & Herschlag, D. Understanding nucleic acid-ion interactions. Annu. Rev. Biochem. 83, 813–841 (2014).24606136 10.1146/annurev-biochem-060409-092720
59. Gebala M Cation-anion interactions within the nucleic acid ion atmosphere revealed by ion counting J. Am. Chem. Soc. 2015 137 14705 14715 10.1021/jacs.5b08395 26517731
Gebala, M. et al. Cation-anion interactions within the nucleic acid ion atmosphere revealed by ion counting. J. Am. Chem. Soc. 137, 14705–14715 (2015).26517731 10.1021/jacs.5b08395
60. Plumridge A Meisburger SP Andresen K Pollack L The impact of base stacking on the conformations and electrostatics of single-stranded DNA Nucleic Acids Res. 2017 45 3932 3943 10.1093/nar/gkx140 28334825
Plumridge, A., Meisburger, S. P., Andresen, K. & Pollack, L. The impact of base stacking on the conformations and electrostatics of single-stranded DNA. Nucleic Acids Res. 45, 3932–3943 (2017).28334825 10.1093/nar/gkx140
61. Zhao J Cramer SM McGown LB Mechanism of sequence-based separation of single-stranded DNA in capillary zone electrophoresis Electrophoresis 2020 41 705 713 10.1002/elps.201900418 32031267
Zhao, J., Cramer, S. M. & McGown, L. B. Mechanism of sequence-based separation of single-stranded DNA in capillary zone electrophoresis. Electrophoresis 41, 705–713 (2020).32031267 10.1002/elps.201900418
62. Jacobson DR Saleh OA Quantifying the ion atmosphere of unfolded, single-stranded nucleic acids using equilibrium dialysis and single-molecule methods Nucleic Acids Res. 2016 44 3763 3771 10.1093/nar/gkw196 27036864
Jacobson, D. R. & Saleh, O. A. Quantifying the ion atmosphere of unfolded, single-stranded nucleic acids using equilibrium dialysis and single-molecule methods. Nucleic Acids Res. 44, 3763–3771 (2016).27036864 10.1093/nar/gkw196
63. Liebl K Zacharias M The development of nucleic acids force fields: From an unchallenged past to a competitive future Biophys. J. 2022 122 2841 2851 10.1016/j.bpj.2022.12.022 36540025
Liebl, K. & Zacharias, M. The development of nucleic acids force fields: From an unchallenged past to a competitive future. Biophys. J. 122, 2841–2851 (2022).36540025 10.1016/j.bpj.2022.12.022
64. Tucker MR Piana S Tan D Levine MV Shaw DE Development of force field parameters for the simulation of single- and double-stranded DNA molecules and DNA–protein complexes J. Phys. Chem. B 2022 126 4442 4457 10.1021/acs.jpcb.1c10971 35694853
Tucker, M. R., Piana, S., Tan, D., Levine, M. V. & Shaw, D. E. Development of force field parameters for the simulation of single- and double-stranded DNA molecules and DNA–protein complexes. J. Phys. Chem. B 126, 4442–4457 (2022).35694853 10.1021/acs.jpcb.1c10971
65. Tan D Piana S Dirks RM Shaw DE RNA force field with accuracy comparable to state-of-the-art protein force fields Proc. Natl. Acad. Sci. 2018 115 E1346 E1355 10.1073/pnas.1713027115 29378935
Tan, D., Piana, S., Dirks, R. M. & Shaw, D. E. RNA force field with accuracy comparable to state-of-the-art protein force fields. Proc. Natl. Acad. Sci. 115, E1346–E1355 (2018).29378935 10.1073/pnas.1713027115
66. He W Naleem N Kleiman D Kirmizialtin S Refining the RNA force field with small-angle X-ray scattering of helix–junction–helix RNA J. Phys. Chem. Lett. 2022 13 3400 3408 10.1021/acs.jpclett.2c00359 35404614
He, W., Naleem, N., Kleiman, D. & Kirmizialtin, S. Refining the RNA force field with small-angle X-ray scattering of helix–junction–helix RNA. J. Phys. Chem. Lett. 13, 3400–3408 (2022).35404614 10.1021/acs.jpclett.2c00359
67. Snodin BEK Introducing improved structural properties and salt dependence into a coarse-grained model of DNA J. Chem. Phys. 2015 142 234901 10.1063/1.4921957 26093573
Snodin, B. E. K. et al. Introducing improved structural properties and salt dependence into a coarse-grained model of DNA. J. Chem. Phys. 142, 234901 (2015).26093573 10.1063/1.4921957
68. Ke C Humeniuk M Gracz H Marszalek PE Direct measurements of base stacking interactions in DNA by single-molecule atomic-force spectroscopy Phys. Rev. Lett. 2007 99 018302 10.1103/PhysRevLett.99.018302 17678193
Ke, C., Humeniuk, M., Gracz, H. & Marszalek, P. E. Direct measurements of base stacking interactions in DNA by single-molecule atomic-force spectroscopy. Phys. Rev. Lett. 99, 018302 (2007).17678193 10.1103/PhysRevLett.99.018302
69. Dustin A Duggan G Gouil Q Omar A Sequence-dependent elasticity and electrostatics of single-stranded DNA: Signatures of base-stacking Biophys. J. 2014 106 659 666 10.1016/j.bpj.2013.12.018 24507606
Dustin, A., Duggan, G., Gouil, Q. & Omar, A. Sequence-dependent elasticity and electrostatics of single-stranded DNA: Signatures of base-stacking. Biophys. J. 106, 659–666 (2014).24507606 10.1016/j.bpj.2013.12.018
70. Jacobson DR McIntosh DB Stevens MJ Rubinstein M Saleh OA Single-stranded nucleic acid elasticity arises from internal electrostatic tension Proc. Natl. Acad. Sci. 2017 114 5095 5100 10.1073/pnas.1701132114 28461493
Jacobson, D. R., McIntosh, D. B., Stevens, M. J., Rubinstein, M. & Saleh, O. A. Single-stranded nucleic acid elasticity arises from internal electrostatic tension. Proc. Natl. Acad. Sci. 114, 5095–5100 (2017).28461493 10.1073/pnas.1701132114
71. Brown RF Andrews CT Elcock AH Stacking free energies of all DNA and RNA nucleoside pairs and dinucleoside-monophosphates computed using recently revised AMBER parameters and compared with experiment J. Chem. Theory Comput. 2015 11 2315 2328 10.1021/ct501170h 26574427
Brown, R. F., Andrews, C. T. & Elcock, A. H. Stacking free energies of all DNA and RNA nucleoside pairs and dinucleoside-monophosphates computed using recently revised AMBER parameters and compared with experiment. J. Chem. Theory Comput. 11, 2315–2328 (2015).26574427 10.1021/ct501170h
72. Manning GS The molecular theory of polyelectrolyte solutions with applications to the electrostatic properties of polynucleotides Q. Rev. Biophys. 1978 11 179 246 10.1017/S0033583500002031 353876
Manning, G. S. The molecular theory of polyelectrolyte solutions with applications to the electrostatic properties of polynucleotides. Q. Rev. Biophys. 11, 179–246 (1978).353876 10.1017/S0033583500002031
73. Manning, G. S. The application of polyelectrolyte limiting laws to the helix–coil transition of DNA. VI. The numerical value of the axial phosphate spacing for the coil form. Biopolymers 15, 2385–2390 (1976).
74. Record MT Woodbury CP Lohman TM Na+ effects on transitions of DNA and polynucleotides of variable linear charge density Biopolymers 1976 15 893 915 10.1002/bip.1976.360150507 4169
Record, M. T., Woodbury, C. P. & Lohman, T. M. Na+ effects on transitions of DNA and polynucleotides of variable linear charge density. Biopolymers 15, 893–915 (1976).4169 10.1002/bip.1976.360150507
75. Chi Q Wang G Jiang J The persistence length and length per base of single-stranded DNA obtained from fluorescence correlation spectroscopy measurements using mean field theory Phys. A Stat. Mech. Appl. 2013 392 1072 1079 10.1016/j.physa.2012.09.022
Chi, Q., Wang, G. & Jiang, J. The persistence length and length per base of single-stranded DNA obtained from fluorescence correlation spectroscopy measurements using mean field theory. Phys. A Stat. Mech. Appl. 392, 1072–1079 (2013).10.1016/j.physa.2012.09.022
76. Špačková B Label-free nanofluidic scattering microscopy of size and mass of single diffusing molecules and nanoparticles Nat. Methods 2022 19 751 758 10.1038/s41592-022-01491-6 35637303
Špačková, B. et al. Label-free nanofluidic scattering microscopy of size and mass of single diffusing molecules and nanoparticles. Nat. Methods 19, 751–758 (2022).35637303 10.1038/s41592-022-01491-6
77. Townshend RJL Geometric deep learning of RNA structure Science 2021 373 1047 1051 10.1126/science.abe5650 34446608
Townshend, R. J. L. et al. Geometric deep learning of RNA structure. Science 373, 1047–1051 (2021).34446608 10.1126/science.abe5650
78. Panchal V Brenk R Riboswitches as drug targets for antibiotics Antibiotics 2021 10 45 10.3390/antibiotics10010045 33466288
Panchal, V. & Brenk, R. Riboswitches as drug targets for antibiotics. Antibiotics. 10, 45 (2021).33466288 10.3390/antibiotics10010045
79. Soukup JK Soukup GA Riboswitches exert genetic control through metabolite-induced conformational change Curr. Opin. Struct. Biol. 2004 14 344 349 10.1016/j.sbi.2004.04.007 15193315
Soukup, J. K. & Soukup, G. A. Riboswitches exert genetic control through metabolite-induced conformational change. Curr. Opin. Struct. Biol. 14, 344–349 (2004).15193315 10.1016/j.sbi.2004.04.007
80. Kramers HA Brownian motion in a field of force and the diffusion model of chemical reactions Physica 1940 7 284 304 10.1016/S0031-8914(40)90098-2
Kramers, H. A. Brownian motion in a field of force and the diffusion model of chemical reactions. Physica 7, 284–304 (1940).10.1016/S0031-8914(40)90098-2
81. Case DA The Amber biomolecular simulation programs J. Comput. Chem. 2005 26 1668 1688 10.1002/jcc.20290 16200636
Case, D. A. et al. The Amber biomolecular simulation programs. J. Comput. Chem. 26, 1668–1688 (2005).16200636 10.1002/jcc.20290
82. Lee J CHARMM-GUI input generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM simulations using the CHARMM36 additive force field J. Chem. Theory Comput. 2016 12 405 413 10.1021/acs.jctc.5b00935 26631602
Lee, J. et al. CHARMM-GUI input generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM simulations using the CHARMM36 additive force field. J. Chem. Theory Comput. 12, 405–413 (2016).26631602 10.1021/acs.jctc.5b00935
83. Abraham MJ GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers SoftwareX 2015 1–2 19 25 10.1016/j.softx.2015.06.001
Abraham, M. J. et al. GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 1–2, 19–25 (2015).10.1016/j.softx.2015.06.001
84. Sengar, A., Ouldridge, T. E., Henrich, O., Rovigatti, L. & Šulc, P. A Primer on the oxDNA model of DNA: When to use it, how to simulate it and how to interpret the results. Front. Mol. Biosci. 8 (2021).
85. Suma A TacoxDNA: A user-friendly web server for simulations of complex DNA structures, from single strands to origami J. Comput. Chem. 2019 40 2586 2595 10.1002/jcc.26029 31301183
Suma, A. et al. TacoxDNA: A user-friendly web server for simulations of complex DNA structures, from single strands to origami. J. Comput. Chem. 40, 2586–2595 (2019).31301183 10.1002/jcc.26029
86. Michaud-Agrawal N Denning EJ Woolf TB Beckstein O MDAnalysis: A toolkit for the analysis of molecular dynamics simulations J. Comput. Chem. 2011 32 2319 2327 10.1002/jcc.21787 21500218
Michaud-Agrawal, N., Denning, E. J., Woolf, T. B. & Beckstein, O. MDAnalysis: A toolkit for the analysis of molecular dynamics simulations. J. Comput. Chem. 32, 2319–2327 (2011).21500218 10.1002/jcc.21787
87. Graen T Hoefling M Grubmüller H AMBER-DYES: Characterization of charge fluctuations and force field parameterization of fluorescent dyes for molecular dynamics simulations J. Chem. Theory Comput. 2014 10 5505 5512 10.1021/ct500869p 26583233
Graen, T., Hoefling, M. & Grubmüller, H. AMBER-DYES: Characterization of charge fluctuations and force field parameterization of fluorescent dyes for molecular dynamics simulations. J. Chem. Theory Comput. 10, 5505–5512 (2014).26583233 10.1021/ct500869p
88. Pau PCF Berg JO McMillan WG Application of Stokes' law to ions in aqueous solution J. Phys. Chem. 1990 94 2671 2679 10.1021/j100369a080
Pau, P. C. F., Berg, J. O. & McMillan, W. G. Application of Stokes’ law to ions in aqueous solution. J. Phys. Chem. 94, 2671–2679 (1990).10.1021/j100369a080
89. Decherchi S Rocchia W A general and robust ray-casting-based algorithm for triangulating surfaces at the nanoscale PLoS ONE 2013 8 e59744 10.1371/journal.pone.0059744 23577073
Decherchi, S. & Rocchia, W. A general and robust ray-casting-based algorithm for triangulating surfaces at the nanoscale. PLoS ONE 8, e59744 (2013).23577073 10.1371/journal.pone.0059744
90. Geuzaine C Remacle J-F Gmsh: A 3-D finite element mesh generator with built-in pre- and post-processing facilities Int. J. Numer. Methods Eng. 2009 79 1309 1331 10.1002/nme.2579
Geuzaine, C. & Remacle, J.-F. Gmsh: A 3-D finite element mesh generator with built-in pre- and post-processing facilities. Int. J. Numer. Methods Eng. 79, 1309–1331 (2009).10.1002/nme.2579
91. Scroggs MW Dokken JS Richardson CN Wells GN Construction of arbitrary order finite element degree-of-freedom maps on polygonal and polyhedral cell meshes ACM Trans. Math. Softw. 2022 48 1 23 10.1145/3524456
Scroggs, M. W., Dokken, J. S., Richardson, C. N. & Wells, G. N. Construction of arbitrary order finite element degree-of-freedom maps on polygonal and polyhedral cell meshes. ACM Trans. Math. Softw. 48, 1–23 (2022).10.1145/3524456
92. Scroggs MW Baratta IA Richardson CN Wells GN Basix: A runtime finite element basis evaluation library J. Open Sour. Softw. 2022 7 3982 10.21105/joss.03982
Scroggs, M. W., Baratta, I. A., Richardson, C. N. & Wells, G. N. Basix: A runtime finite element basis evaluation library. J. Open Sour. Softw. 7, 3982 (2022).10.21105/joss.03982
93. Alnæs MS Logg A Ølgaard KB Rognes ME Wells GN Unified form language ACM Trans. Math. Softw. 2014 40 1 37 10.1145/2566630
Alnæs, M. S., Logg, A., Ølgaard, K. B., Rognes, M. E. & Wells, G. N. Unified form language. ACM Trans. Math. Softw. 40, 1–37 (2014).10.1145/2566630
94. Theodoor J Overbeek G The role of energy and entropy in the electrical double layer Colloids Surf. 1990 51 61 75 10.1016/0166-6622(90)80132-N
Theodoor, J. & Overbeek, G. The role of energy and entropy in the electrical double layer. Colloids Surf. 51, 61–75 (1990).10.1016/0166-6622(90)80132-N
95. Ahrens, J. P., Geveci, B. & Law, C. C. in The Visualization Handbook. (2005).
96. Tomasello G Armenia I Molla G The Protein Imager: A full-featured online molecular viewer interface with server-side HQ-rendering capabilities Bioinformatics 2020 36 2909 2911 10.1093/bioinformatics/btaa009 31930403
Tomasello, G., Armenia, I. & Molla, G. The Protein Imager: A full-featured online molecular viewer interface with server-side HQ-rendering capabilities. Bioinformatics 36, 2909–2911 (2020).31930403 10.1093/bioinformatics/btaa009
97. Pless RC Ts'o PO Duplex formation of a nonionic oligo(deoxythymidylate) analogue (heptadeoxythymidylyl-(3'-5')-deoxythymidine heptaethyl ester (d-(Tp(Et))7T)) with poly(deoxyadenylate). Evaluation of the electrostatic interaction Biochemistry 1977 16 1239 1250 10.1021/bi00625a033 14676
Pless, R. C. & Ts’o, P. O. Duplex formation of a nonionic oligo(deoxythymidylate) analogue (heptadeoxythymidylyl-(3’-5’)-deoxythymidine heptaethyl ester (d-(Tp(Et))7T)) with poly(deoxyadenylate). Evaluation of the electrostatic interaction. Biochemistry 16, 1239–1250 (1977).14676 10.1021/bi00625a033
98. Hansma HG Revenko I Kim K Laney DE Atomic force microscopy of long and short double-stranded, single-stranded and triple-stranded nucleic acids Nucleic Acids Res. 1996 24 713 720 10.1093/nar/24.4.713 8604315
Hansma, H. G., Revenko, I., Kim, K. & Laney, D. E. Atomic force microscopy of long and short double-stranded, single-stranded and triple-stranded nucleic acids. Nucleic Acids Res. 24, 713–720 (1996).8604315 10.1093/nar/24.4.713
99. Adamcik J Klinov DV Witz G Sekatskii SK Dietler G Observation of single-stranded DNA on mica and highly oriented pyrolytic graphite by atomic force microscopy FEBS Lett. 2006 580 5671 5675 10.1016/j.febslet.2006.09.017 17007844
Adamcik, J., Klinov, D. V., Witz, G., Sekatskii, S. K. & Dietler, G. Observation of single-stranded DNA on mica and highly oriented pyrolytic graphite by atomic force microscopy. FEBS Lett. 580, 5671–5675 (2006).17007844 10.1016/j.febslet.2006.09.017
100. Mills JB Vacano E Hagerman PJ Flexibility of single-stranded DNA: use of gapped duplex helices to determine the persistence lengths of Poly(dT) and Poly(dA)11Edited by B Honig J. Mol. Biol. 1999 285 245 257 10.1006/jmbi.1998.2287 9878403
Mills, J. B., Vacano, E. & Hagerman, P. J. Flexibility of single-stranded DNA: use of gapped duplex helices to determine the persistence lengths of Poly(dT) and Poly(dA). J. Mol. Biol. 285, 245–257 (1999).9878403 10.1006/jmbi.1998.2287
