Polyelectrolyte properties of single stranded DNA measured using SAXS and single-molecule FRET: Beyond the wormlike chain model

Sensing the structural and conformational properties of single-stranded nucleic acids using electrometry and molecular simulations

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.

Nanosecond chain dynamics of single-stranded nucleic acids

The conformational dynamics of single-stranded nucleic acids are fundamental for nucleic acid folding and function. However, their elementary chain dynamics have been difficult to resolve experimentally. Here we employ a combination of single-molecule Förster resonance energy transfer, nanosecond fluorescence correlation spectroscopy, and nanophotonic enhancement to determine the conformational ensembles and rapid chain dynamics of short single-stranded nucleic acids in solution. To interpret the experimental results in terms of end-to-end distance dynamics, we utilize the hierarchical chain growth approach, simple polymer models, and refinement with Bayesian inference to generate structural ensembles that closely align with the experimental data. The resulting chain reconfiguration times are exceedingly rapid, in the 10-ns range. Solvent viscosity-dependent measurements indicate that these dynamics of single-stranded nucleic acids exhibit negligible internal friction and are thus dominated by solvent friction. Our results provide a detailed view of the conformational distributions and rapid dynamics of single-stranded nucleic acids.

Na+-Mg2+ ion effects on conformation and translocation dynamics of single-stranded RNA: Cooperation and competition

The nanopore-based translocation of a single-stranded RNA (ssRNA) in mixed salt solution has garnered increasing interest for its biological and technological significance. However, it is challenging to comprehensively understand the effects of the mixed ion species on the translocation dynamics due to their cooperation and competition, which can be directly reflected by the ion screening and neutralizing effects, respectively. In this study, Langevin dynamics simulation is employed to investigate the properties of ssRNA conformation and translocation in mixed Na+-Mg2+ ion environments. Simulation results reveal that the ion screening effect dominates the change in the ssRNA conformational size, the ion neutralizing effect controls the capture rate of the ssRNA by the nanopore, and both of them take charge of the different changes in translocation time of the ssRNA under various mixed ion environments. Under high Na+ ion concentration, as Mg2+ concentration increases, the ion neutralizing effect strengthens, weakening the driving force inside the nanopore, leading to longer translocation time. Conversely, at low Na+ concentration, an increase in Mg2+ concentration enhances the ion screening effect, aiding in faster translocation. Furthermore, these simulation results will be explained by quantitative analysis, advancing a deeper understanding of the complicated effects of the mixed Na+-Mg2+ ions.

Competition between Stacking and Divalent Cation-Mediated Electrostatic Interactions Determines the Conformations of Short DNA Sequences

Interplay between divalent cations (Mg2+ and Ca2+) and single-stranded DNA (ssDNA) and double-stranded DNA (dsDNA), as well as stacking interactions, is important in nucleosome stability and phase separation in nucleic acids. Quantitative techniques accounting for ion-DNA interactions are needed to obtain insights into these and related problems. Toward this end, we created a sequence-dependent computational TIS-ION model that explicitly accounts for monovalent and divalent ions. Simulations of the rigid 24 base-pair (bp) dsDNA and flexible ssDNA sequences, dT30 and dA30, with varying amounts of the divalent cations show that the calculated excess number of ions around the dsDNA and ssDNA agree quantitatively with ion-counting experiments. Using an ensemble of all-atom structures generated from coarse-grained simulations, we calculated the small-angle X-ray scattering profiles, which are in excellent agreement with experiments. Although ion-counting experiments mask the differences between Mg2+ and Ca2+, we find that Mg2+ binds to the minor grooves and phosphate groups, whereas Ca2+ binds specifically to the minor groove. Both Mg2+ and Ca2+ exhibit a tendency to bind to the minor groove of DNA as opposed to the major groove. The dA30 conformations are dominated by stacking interactions, resulting in structures with considerable helical order. The near cancellation of the favorable stacking and unfavorable electrostatic interactions leads to dT30 populating an ensemble of heterogeneous conformations. The successful applications of the TIS-ION model are poised to confront many problems in DNA biophysics.

Can Hindered Transport Models for Rigid Spheres Predict the Rejection of Single Stranded DNA from Porous Membranes?

In this paper, predictions from a theoretical model describing the rejection of a rigid spherical solute from porous membranes are compared to experimental results for a single stranded DNA (ssDNA) with 60 thymine nucleotides. Experiments were conducted with different pore size track-etched membranes at different transmembrane pressures and different NaCl concentrations. The model includes both hydrodynamic and electrostatic solute-pore wall interactions; predictions were made using different size parameters for the ssDNA (radius of gyration, hydrodynamic radius, and root mean square end-to-end distance). At low transmembrane pressures, experimental results are in good agreement with rejection predictions made using the hard sphere model for the ssDNA when the solute size is described using its root mean square end-to-end distance. When the ssDNA size is characterized using the radius of gyration or the hydrodynamic radius, the hard sphere model under-predicts rejection. Not surprisingly, the model overestimates ssDNA rejection at conditions where flow induced elongation of the DNA is expected. The results from this study are encouraging because they mean that a relatively simple hindered transport model can be used to estimate the rejection of a small DNA from porous membranes.

Molecular Structure of Single-Stranded DNA on the ZnS Surface of Quantum Dots

DNA-based nanoparticle assemblies have emerged as leading candidates in the development of bioimaging materials, photonic devices, and computing materials. Here, we combine atomistic simulations and experiments to characterize the wrapping mechanism of chimeric single-stranded DNA (ssDNA) on CdSe-ZnS (core-shell) quantum dots (QDs) at different ratios of the phosphorothioate (PS) modification of the bases. We use an implicit solvent, all-atom ssDNA model to match the experimentally calculated ssDNA conformation at low salt concentrations. Through simulation, we find that 3-mercaptopropionic acid (MPA) induces electrostatic repulsion and O-(2-mercaptoethyl)-Ó-methyl-hexa (ethylene glycol) (mPEG) induces steric exclusion, and both reduce the binding affinity of ssDNA. In both simulation and experiment, we find that ssDNA is closer to the QD surface when the QD size is larger. The effect of the PS-base ratio on the conformation of ssDNA is also elaborated in this work. We found through MD simulations, and confirmed by transmission electron microscopy, that the maximum valence numbers are 1, 2, and 3 on QDs of 6, 9, and 14 nm in diameter, respectively. We conclude that the maximum ssDNA valence number is linearly related to the QD size, n ∝ R, and justify this finding through an electrostatic repulsion mechanism.

Assessment of AMBER Force Fields for Simulations of ssDNA

Single-stranded DNA (ssDNA) plays an important role in biological processes and is used in DNA nanotechnology and other novel applications. Many important research questions can be addressed with molecular simulations of ssDNA molecules; however, no dedicated force field for ssDNA has been developed, and there is limited experimental information about ssDNA structures. This study assesses the accuracy and applicability of existing Amber force fields for all-atom simulations of ssDNA, such as ff99, bsc0, bsc1, and OL15, in implicit and explicit solvents via comparison to available experimental data, such as Forster resonance energy transfer and small angle X-ray scattering. We observed that some force fields agree better with experiments than others mainly due to the difference in parameterization of the propensity for hydrogen bonding and base stacking. Overall, the Amber ff99 force field in the IGB5 or IGB8 implicit solvent and the bsc1 force field in the explicit TIP3P solvent had the best agreement with experiment.

Structure-based screening of binding affinities via small-angle X-ray scattering

Protein-protein and protein-ligand interactions often involve conformational changes or structural rearrangements that can be quantified by solution small-angle X-ray scattering (SAXS). These scattering intensity measurements reveal structural details of the bound complex, the number of species involved and, additionally, the strength of interactions if carried out as a titration. Although a core part of structural biology workflows, SAXS-based titrations are not commonly used in drug discovery contexts. This is because prior knowledge of expected sample requirements, throughput and prediction accuracy is needed to develop reliable ligand screens. This study presents the use of the histidine-binding protein (26 kDa) and other periplasmic binding proteins to benchmark ligand screen performance. Sample concentrations and exposure times were varied across multiple screening trials at four beamlines to investigate the accuracy and precision of affinity prediction. The volatility ratio between titrated scattering curves and a common apo reference is found to most reliably capture the extent of structural and population changes. This obviates the need to explicitly model scattering intensities of bound complexes, which can be strongly ligand-dependent. Where the dissociation constant is within 102 of the protein concentration and the total exposure times exceed 20 s, the titration protocol presented at 0.5 mg ml-1 yields affinities comparable to isothermal titration calorimetry measurements. Estimated throughput ranges between 20 and 100 ligand titrations per day at current synchrotron beamlines, with the limiting step imposed by sample handling and cleaning procedures.

Mechanism of sequence-based separation of single-stranded DNA in capillary zone electrophoresis

Separation of DNA by length using CGE is a mature field. Separation of DNA by sequence, in contrast, is a more difficult problem. Existing techniques generally rely upon changes in intrinsic or induced differences in conformation. Previous work in our group showed that sets of ssDNA of the same length differing in sequence by as little as a single base could be separated by CZE using simple buffers at high ionic strength. Here, we explore the basis of the separation using circular dichroism spectroscopy, fluorescence anisotropy, and small angle X-ray scattering. The results reveal sequence-dependent differences among the same length strands, but the trends in the differences are not correlated to the migration order of the strands in the CZE separation. They also indicate that the separation is based on intrinsic differences among the strands that do not change with increasing ionic strength; rather, increasing ionic strength has a greater effect on electroosmotic mobility in the normal direction than on electrophoretic mobility of the strands in the reverse direction. This increases the migration time of the strands in the normal direction, allowing more time for the same-length strands to be teased apart based on very small differences in the intrinsic properties of the strands of different sequence. Regression analysis was used to model the intrinsic differences among DNA strands in order to gain insight into the relationship between mobility and sequence that underlies the separation.

Visualizing Disordered Single-Stranded RNA: Connecting Sequence, Structure, and Electrostatics

Disordered homopolymeric regions of single-stranded RNA, such as U or A tracts, are found within functional RNAs where they play distinct roles in defining molecular structure and facilitating recognition by partners. Despite this prominence, details of conformational and biophysical properties of these regions have not yet been resolved. We apply a number of experimental techniques to investigate the conformations of these biologically important motifs and provide quantitative measurements of their ion atmospheres. Single strands of RNA display pronounced sequence-dependent conformations that relate to the unique ion atmospheres each attracts. Chains of rU bases are relatively unstructured under all conditions, while chains of rA bases display distinct ordering through stacking or clustering motifs, depending on the composition of the surrounding solution. These dramatic structural differences are consistent with the measured disparity in ion composition and atmospheres around each homopolymer, revealing a complex interplay of base, ion, and single-strand ordering. The unique structural and ionic signatures of homopolymer ssRNAs explains their role(s) in folding structured RNAs and may explain their distinct recognition by protein partners.

Structure, stability and specificity of the binding of ssDNA and ssRNA with proteins

Recognition of single-stranded DNA (ssDNA) or single-stranded RNA (ssRNA) is important for many fundamental cellular functions. A variety of single-stranded DNA-binding proteins (ssDBPs) and single-stranded RNA-binding proteins (ssRBPs) have evolved that bind ssDNA and ssRNA, respectively, with varying degree of affinities and specificities to form complexes. Structural studies of these complexes provide key insights into their recognition mechanism. However, computational modeling of the specific recognition process and to predict the structure of the complex is challenging, primarily due to the heterogeneity of their binding energy landscape and the greater flexibility of ssDNA or ssRNA compared with double-stranded nucleic acids. Consequently, considerably fewer computational studies have explored interactions between proteins and single-stranded nucleic acids compared with protein interactions with double-stranded nucleic acids. Here, we report a newly developed energy-based coarse-grained model to predict the structure of ssDNA–ssDBP and ssRNA–ssRBP complexes and to assess their sequence-specific interactions and stabilities. We tuned two factors that can modulate specific recognition: base–aromatic stacking strength and the flexibility of the single-stranded nucleic acid. The model was successfully applied to predict the binding conformations of 12 distinct ssDBP and ssRBP structures with their cognate ssDNA and ssRNA partners having various sequences. Estimated binding energies agreed well with the corresponding experimental binding affinities. Bound conformations from the simulation showed a funnel-shaped binding energy distribution where the native-like conformations corresponded to the energy minima. The various ssDNA–protein and ssRNA–protein complexes differed in the balance of electrostatic and aromatic energies. The lower affinity of the ssRNA–ssRBP complexes compared with the ssDNA–ssDBP complexes stems from lower flexibility of ssRNA compared to ssDNA, which results in higher rate constants for the dissociation of the complex (k_off) for complexes involving the former.

Quantifying bimolecular self-assembly is pivotal to understanding cellular function. In recent years, a large progress has been made in understanding the structure and biophysics of protein-protein interactions. Particularly, various computational tools are available for predicting these structures and to estimate their stability and the driving forces of their formation. The understating of the interactions between proteins and nucleic acids, however, is still limited, presumably due to the involvement of non-specific interactions as well as the high conformational plasticity that may demand an induced-fit mechanism. In particular, the interactions between proteins and single-stranded nucleic acids (i.e., single-stranded DNA and RNA) is very challenging due to their high flexibility. Furthermore, the interface between proteins and single-stranded nucleic acids is often chemically more heterogeneous than the interface between proteins and double-stranded DNA. In this study, we developed a coarse-grained computational model to predict the structure of complexes between proteins and single-stranded nucleic acids. The model was applied to estimate binding affinities and the estimated binding energies agreed well with the corresponding experimental binding affinities. The kinetics of association as well as the specificity of the complexes between proteins and ssDNA are different than those with ssRNA, mostly due to differences in their conformational flexibility.

Dispersion Correction Alleviates Dye Stacking of Single-Stranded DNA and RNA in Simulations of Single-Molecule Fluorescence Experiments

We combine single-molecule Förster resonance energy transfer (single-molecule FRET) experiments with extensive all-atom molecular dynamics (MD) simulations (>100 μs) to characterize the conformational ensembles of single-stranded (ss) DNA and RNA in solution. From MD simulations with explicit dyes attached to single-stranded nucleic acids via flexible linkers, we calculate FRET efficiencies and fluorescence anisotropy decays. We find that dispersion-corrected water models alleviate the problem of overly abundant interactions between fluorescent dyes and the aromatic ring systems of nucleobases. To model dye motions in a computationally efficient and conformationally exhaustive manner, we introduce a dye-conformer library, built from simulations of dinucleotides with covalently attached dye molecules. We use this library to calculate FRET efficiencies for dT₁₉, dA₁₉, and rA₁₉ simulated without explicit labels over a wide range of salt concentrations. For end-labeled homopolymeric pyrimidine ssDNA, MD simulations with the parmBSC1 force field capture the overall trend in salt-dependence of single-molecule FRET based distance measurements. For homopolymeric purine ssRNA and ssDNA, the DESRES and parmBSC1 force fields, respectively, provide useful starting points, even though our comparison also identifies clear deviations from experiment.

The role of small-angle scattering in structure-based screening applications

In many biomolecular interactions, changes in the assembly states and structural conformations of participants can act as a complementary reporter of binding to functional and thermodynamic assays. This structural information is captured by a number of structural biology and biophysical techniques that are viable either as primary screens in small-scale applications or as secondary screens to complement higher throughput methods. In particular, small-angle X-ray scattering (SAXS) reports the average distance distribution between all atoms after orientational averaging. Such information is important when for example investigating conformational changes involved in inhibitory and regulatory mechanisms where binding events do not necessarily cause functional changes. Thus, we summarise here the current and prospective capabilities of SAXS-based screening in the context of other methods that yield structural information. Broad guidelines are also provided to assist readers in preparing screening protocols that are tailored to available X-ray sources.

Enzymatic Synthesis of Nucleobase-Modified Single-Stranded DNA Offers Tunable Resistance to Nuclease Degradation

We synthesized long, nucleobase-modified, single-stranded DNA (ssDNA) using terminal deoxynucleotidyl transferase (TdT) enzymatic polymerization. Specifically, we investigated the effect of unnatural nucleobase size and incorporation density on ssDNA resistance to exo- and endonuclease degradation. We discovered that increasing the size and density of unnatural nucleobases enhances ssDNA resistance to degradation in the presence of exonuclease I, DNase I, and human serum. We also studied the mechanism of this resistance enhancement using molecular dynamics simulations. Our results show that the presence of unnatural nucleobases in ssDNA decreases local chain flexibility and hampers nuclease access to the ssDNA backbone, which hinders nuclease binding to ssDNA and slows its degradation. Our discoveries suggest that incorporating nucleobase-modified nucleotides into ssDNA, using enzymatic polymerization, is an easy and efficient strategy to prolong and tune the half-life of DNA-based materials in nucleases-containing environments.

The impact of base stacking on the conformations and electrostatics of single-stranded DNA

Single-stranded DNA (ssDNA) is notable for its interactions with ssDNA binding proteins (SSBs) during fundamentally important biological processes including DNA repair and replication. Previous work has begun to characterize the conformational and electrostatic properties of ssDNA in association with SSBs. However, the conformational distributions of free ssDNA have been difficult to determine. To capture the vast array of ssDNA conformations in solution, we pair small angle X-ray scattering with novel ensemble fitting methods, obtaining key parameters such as the size, shape and stacking character of strands with different sequences. Complementary ion counting measurements using inductively coupled plasma atomic emission spectroscopy are employed to determine the composition of the ion atmosphere at physiological ionic strength. Applying this combined approach to poly dA and poly dT, we find that the global properties of these sequences are very similar, despite having vastly different propensities for single-stranded helical stacking. These results suggest that a relatively simple mechanism for the binding of ssDNA to non-specific SSBs may be at play, which explains the disparity in binding affinities observed for these systems.

Counting the ions surrounding nucleic acids

Nucleic acids are strongly negatively charged, and thus electrostatic interactions—screened by ions in solution—play an important role in governing their ability to fold and participate in biomolecular interactions. The negative charge creates a region, known as the ion atmosphere, in which cation and anion concentrations are perturbed from their bulk values. Ion counting experiments quantify the ion atmosphere by measuring the preferential ion interaction coefficient: the net total number of excess ions above, or below, the number expected due to the bulk concentration. The results of such studies provide important constraints on theories, which typically predict the full three-dimensional distribution of the screening cloud. This article reviews the state of nucleic acid ion counting measurements and critically analyzes their ability to test both analytical and simulation-based models.

Martini Coarse-Grained Force Field: Extension to RNA

RNA has an important role not only as the messenger of genetic information but also as a regulator of gene expression. Given its central role in cell biology, there is significant interest in studying the structural and dynamic behavior of RNA in relation to other biomolecules. Coarse-grain molecular dynamics simulations are a key tool to that end. Here, we have extended the coarse-grain Martini force field to include RNA after our recent extension to DNA. In the same way DNA was modeled, the tertiary structure of RNA is constrained using an elastic network. This model, therefore, is not designed for applications involving RNA folding but rather offers a stable RNA structure for studying RNA interactions with other (bio)molecules. The RNA model is compatible with all other Martini models and opens the way to large-scale explicit-solvent molecular dynamics simulations of complex systems involving RNA.

Visualizing single-stranded nucleic acids in solution

Single-stranded nucleic acids (ssNAs) are ubiquitous in many key cellular functions. Their flexibility limits both the number of high-resolution structures available, leaving only a small number of protein-ssNA crystal structures, while forcing solution investigations to report ensemble averages. A description of the conformational distributions of ssNAs is essential to more fully characterize biologically relevant interactions. We combine small angle X-ray scattering (SAXS) with ensemble-optimization methods (EOM) to dynamically build and refine sets of ssNA structures. By constructing candidate chains in representative dinucleotide steps and refining the models against SAXS data, a broad array of structures can be obtained to match varying solution conditions and strand sequences. In addition to the distribution of large scale structural parameters, this approach reveals, for the first time, intricate details of the phosphate backbone and underlying strand conformations. Such information on unperturbed strands will critically inform a detailed understanding of an array of problems including protein-ssNA binding, RNA folding and the polymer nature of NAs. In addition, this scheme, which couples EOM selection with an iteratively refining pool to give confidence in the underlying structures, is likely extendable to the study of other flexible systems.

Conformational Changes and Flexibility of DNA Devices Observed by Small-Angle X-ray Scattering

Self-assembled DNA origami nanostructures enable the creation of precisely defined shapes at the molecular scale. Dynamic DNA devices that are capable of switching between defined conformations could afford completely novel functionalities for diagnostic, therapeutic, or engineering applications. Developing such objects benefits strongly from experimental feedback about conformational changes and 3D structures, ideally in solution, free of potential biases from surface attachment or labeling. Here, we demonstrate that small-angle X-ray scattering (SAXS) can quantitatively resolve the conformational changes of a DNA origami two-state switch device as a function of the ionic strength of the solution. In addition, we show how SAXS data allow for refinement of the predicted idealized three-dimensional structure of the DNA object using a normal mode approach based on an elastic network model. The results reveal deviations from the idealized design geometries that are otherwise difficult to resolve. Our results establish SAXS as a powerful tool to investigate conformational changes and solution structures of DNA origami and we anticipate our methodology to be broadly applicable to increasingly complex DNA and RNA devices.

Quantifying the ion atmosphere of unfolded, single-stranded nucleic acids using equilibrium dialysis and single-molecule methods

To form secondary structure, nucleic acids (NAs) must overcome electrostatic strand–strand repulsion, which is moderated by the surrounding atmosphere of screening ions. The free energy of NA folding therefore depends on the interactions of this ion atmosphere with both the folded and unfolded states. We quantify such interactions using the preferential ion interaction coefficient or ion excess: the number of ions present near the NA in excess of the bulk concentration. The ion excess of the folded, double-helical state has been extensively studied; however, much less is known about the salt-dependent ion excess of the unfolded, single-stranded state. We measure this quantity using three complementary approaches: a direct approach of Donnan equilibrium dialysis read out by atomic emission spectroscopy and two indirect approaches involving either single-molecule force spectroscopy or existing thermal denaturation data. The results of these three approaches, each involving an independent experimental technique, are in good agreement. Even though the single-stranded NAs are flexible polymers that are expected to adopt random-coil configurations, we find that their ion atmosphere is quantitatively described by rod-like models that neglect large-scale conformational freedom, an effect that we explain in terms of the competition between the relevant structural and electrostatic length scales.

Determining the Locations of Ions and Water around DNA from X-Ray Scattering Measurements

Nucleic acids carry a negative charge, attracting salt ions and water. Interactions with these components of the solvent drive DNA to condense, RNA to fold, and proteins to bind. To understand these biological processes, knowledge of solvent structure around the nucleic acids is critical. Yet, because they are often disordered, ions and water evade detection by x-ray crystallography and other high-resolution methods. Small-angle x-ray scattering (SAXS) is uniquely sensitive to the spatial correlations between solutes and the surrounding solvent. Thus, SAXS provides an experimental constraint to guide or test emerging solvation theories. However, the interpretation of SAXS profiles is nontrivial because of the difficulty in separating the scattering signals of each component: the macromolecule, ions, and hydration water. Here, we demonstrate methods for robustly deconvoluting these signals, facilitating a more straightforward comparison with theory. Using SAXS data collected on an absolute intensity scale for short DNA duplexes in solution with Na(+), K(+), Rb(+), or Cs(+) counterions, we mathematically decompose the scattering profiles into components (DNA, water, and ions) and validate the decomposition using anomalous scattering measurements. In addition, we generate a library of physically motivated ion atmosphere models and rank them by agreement with the scattering data. The best-fit models have relatively compact ion atmospheres when compared to predictions from the mean-field Poisson-Boltzmann theory of electrostatics. Thus, the x-ray scattering methods presented here provide a valuable measurement of the global structure of the ion atmosphere that can be used to test electrostatics theories that go beyond the mean-field approximation.

Integrating molecular dynamics simulations with chemical probing experiments using SHAPE-FIT

Integration and calibration of molecular dynamics simulations with experimental data remain a challenging endeavor. We have developed a novel method to integrate chemical probing experiments with molecular simulations of RNA molecules by using a native structure-based model. Selective 2'-hydroxyl acylation by primer extension (SHAPE) characterizes the mobility of each residue in the RNA. Our method, SHAPE-FIT, automatically optimizes the potential parameters of the force field according to measured reactivities from SHAPE. The optimized parameter set allows simulations of dynamics highly consistent with SHAPE probing experiments. Such atomistic simulations, thoroughly grounded in experiment, can open a new window on RNA structure-function relations.

RNA-guided assembly of Rev-RRE nuclear export complexes

HIV replication requires nuclear export of unspliced and singly spliced viral transcripts. Although a unique RNA structure has been proposed for the Rev-response element (RRE) responsible for viral mRNA export, how it recruits multiple HIV Rev proteins to form an export complex has been unclear. We show here that initial binding of Rev to the RRE triggers RNA tertiary structural changes, enabling further Rev binding and the rapid formation of a viral export complex. Analysis of the Rev-RRE assembly pathway using SHAPE-Seq and small-angle X-ray scattering (SAXS) reveals two major steps of Rev-RRE complex formation, beginning with rapid Rev binding to a pre-organized region presenting multiple Rev binding sites. This step induces long-range remodeling of the RNA to expose a cryptic Rev binding site, enabling rapid assembly of additional Rev proteins into the RNA export complex. This kinetic pathway may help maintain the balance between viral replication and maturation.DOI: http://dx.doi.org/10.7554/eLife.03656.001.

Many-body effect in ion binding to RNA

Ion-mediated electrostatic interactions play an important role in RNA folding stability. For a RNA in a solution with higher Mg(2+) ion concentration, more counterions in the solution can bind to the RNA, causing a strong many-body coupling between the bound ions. The many-body effect can change the effective potential of mean force between the tightly bound ions. This effect tends to dampen ion binding and lower RNA folding stability. Neglecting the many-body effect leads to a systematic error (over-estimation) of RNA folding stability at high Mg(2+) ion concentrations. Using the tightly bound ion model combined with a conformational ensemble model, we investigate the influence of the many-body effect on the ion-dependent RNA folding stability. Comparisons with the experimental data indicate that including the many-body effect led to much improved predictions for RNA folding stability at high Mg(2+) ion concentrations. The results suggest that the many-body effect can be important for RNA folding in high concentrations of multivalent ions. Further investigation showed that the many-body effect can influence the spatial distribution of the tightly bound ions and the effect is more pronounced for compact RNA structures and structures prone to the formation of local clustering of ions.

Multi-enzyme complexes on DNA scaffolds capable of substrate channelling with an artificial swinging arm

Swinging arms are a key functional component of multistep catalytic transformations in many naturally occurring multi-enzyme complexes. This arm is typically a prosthetic chemical group that is covalently attached to the enzyme complex via a flexible linker, allowing the direct transfer of substrate molecules between multiple active sites within the complex. Mimicking this method of substrate channelling outside the cellular environment requires precise control over the spatial parameters of the individual components within the assembled complex. DNA nanostructures can be used to organize functional molecules with nanoscale precision and can also provide nanomechanical control. Until now, protein-DNA assemblies have been used to organize cascades of enzymatic reactions by controlling the relative distance and orientation of enzymatic components or by facilitating the interface between enzymes/cofactors and electrode surfaces. Here, we show that a DNA nanostructure can be used to create a multi-enzyme complex in which an artificial swinging arm facilitates hydride transfer between two coupled dehydrogenases. By exploiting the programmability of DNA nanostructures, key parameters including position, stoichiometry and inter-enzyme distance can be manipulated for optimal activity.