Overstretching Double-Stranded RNA, Double-Stranded DNA, and RNA-DNA Duplexes

A predictive model for the thermomechanical melting transition of double stranded DNA

By extending the classical Peyrard-Bishop model, we are able to obtain a fully analytical description for the mechanical response of DNA under stretching at variable values of temperature, number of base pairs and intrachains and interchains bonds stiffness. In order to compare elasticity and temperature effects, we first analyze the system in the zero temperature mechanical limit, important to describe several experimental effects including possible hysteresis. We then analyze temperature effects in the framework of equilibrium Statistical Mechanics. In particular, we obtain an analytical expression for the temperature-dependent melting force and unzipping assigned displacement in the thermodynamical limit, also depending on the relative stability of intra vs. inter molecular bonds. Such results coincide with the purely mechanical model in the limit of zero temperature and with the denaturation temperature that we obtain with the classical transfer integral method. Based on our analytical results, we obtain explicitly phase diagrams and cooperativity parameters, where also discreteness effect can be accounted for. The obtained results are successfully applied in reproducing the thermomechanical experimental melting of DNA and the response of DNA hairpins. Due to the generality of the model, exemplified in the proposed analysis of both overstretching and unzipping experiments, we argue that the proposed approach can be extended to other thermomechanically induced molecular melting phenomena. STATEMENT OF SIGNIFICANCE: We obtain a fully analytical description of the complex wiggly energy landscape of two stranded macromolecules under unzipping loading. Based on Equilibrium Statistical Mechanics, we describe the combined thermomechanical effects and the melting transition of double stranded molecules such as nucleic acids. This is proved by quantitatively predicting the experimental behavior of both melting of DNA and DNA hairpins opening. While analytical results have been previously attained under special conditions on the relative stiffness of the covalent vs. non-covalent bonds of the base pairs, our model is completely general in this respect, thus representing a tool in the perspective of the design at the molecular scale. We show that the obtained model can be fully inscribed in the theory of phase transitions giving a new interpretation of the thermomechanical behavior of double stranded molecules.

Structural domains of SARS-CoV-2 nucleocapsid protein coordinate to compact long nucleic acid substrates

The SARS-CoV-2 nucleocapsid (N) protein performs several functions including binding, compacting, and packaging the ∼30 kb viral genome into the viral particle. N protein consists of two ordered domains, with the N terminal domain (NTD) primarily associated with RNA binding and the C terminal domain (CTD) primarily associated with dimerization/oligomerization, and three intrinsically disordered regions, an N-arm, a C-tail, and a linker that connects the NTD and CTD. We utilize an optical tweezers system to isolate a long single-stranded nucleic acid substrate to measure directly the binding and packaging function of N protein at a single molecule level in real time. We find that N protein binds the nucleic acid substrate with high affinity before oligomerizing and forming a highly compact structure. By comparing the activities of truncated protein variants missing the NTD, CTD, and/or linker, we attribute specific steps in this process to the structural domains of N protein, with the NTD driving initial binding to the substrate and ensuring high localized protein density that triggers interprotein interactions mediated by the CTD, which forms a compact and stable protein-nucleic acid complex suitable for packaging into the virion.

Stretching of Long Double-Stranded DNA and RNA Described by the Same Approach

We propose an approach to help interpret polymer force-extension curves that exhibit plateau regimes. When coupled to a bead-spring dynamic model, the approach accurately reproduces a variety of experimental force-extension curves of long double-stranded DNA and RNA, including torsionally constrained and unconstrained DNA and negatively supercoiled DNA. A key feature of the model is a specific nonconvex energy function of the spring. We provide an algorithm to obtain the five required parameters of the model from experimental force-extension curves. The applicability of the approach to the force-extension curves of double-stranded (ds) DNA of variable GC content as well as to a DNA/RNA hybrid structure is explored and confirmed. We use the approach to explain counterintuitive sequence-dependent trends and make predictions. In the plateau region of the force-extension curves, our molecular dynamics simulations show that the polymer separates into a mix of weakly and strongly stretched states without forming macroscopically distinct phases. The distribution of these states is predicted to depend on the sequence.

The nucleotide addition cycle of the SARS-CoV-2 polymerase

Coronaviruses have evolved elaborate multisubunit machines to replicate and transcribe their genomes. Central to these machines are the RNA-dependent RNA polymerase subunit (nsp12) and its intimately associated cofactors (nsp7 and nsp8). We use a high-throughput magnetic-tweezers approach to develop a mechanochemical description of this core polymerase. The core polymerase exists in at least three catalytically distinct conformations, one being kinetically consistent with incorporation of incorrect nucleotides. We provide evidence that the RNA-dependent RNA polymerase (RdRp) uses a thermal ratchet instead of a power stroke to transition from the pre- to post-translocated state. Ultra-stable magnetic tweezers enable the direct observation of coronavirus polymerase deep and long-lived backtracking that is strongly stimulated by secondary structures in the template. The framework we present here elucidates one of the most important structure-dynamics-function relationships in human health today and will form the grounds for understanding the regulation of this complex.

Formation of frameshift-stimulating RNA pseudoknots is facilitated by remodeling of their folding intermediates

Programmed –1 ribosomal frameshifting is an essential regulation mechanism of translation in viruses and bacteria. It is stimulated by mRNA structures inside the coding region. As the structure is unfolded repeatedly by consecutive translating ribosomes, whether it can refold properly each time is important in performing its function. By using single-molecule approaches and molecular dynamics simulations, we found that a frameshift-stimulating RNA pseudoknot folds sequentially through its upstream stem S1 and downstream stem S2. In this pathway, S2 folds from the downstream side and tends to be trapped in intermediates. By masking the last few nucleotides to mimic their gradual emergence from translating ribosomes, S2 can be directed to fold from the upstream region. The results show that the intermediates are greatly suppressed, suggesting that mRNA refolding may be modulated by ribosomes. Moreover, masking the first few nucleotides of S1 favors the folding from S2 and yields native pseudoknots, which are stable enough to retrieve the masked nucleotides. We hypothesize that translating ribosomes can remodel an intermediate mRNA structure into a stable conformation, which may in turn stimulate backward slippage of the ribosome. This supports an interactive model of ribosomal frameshifting and gives an insightful account addressing previous experimental observations.

The nucleotide addition cycle of the SARS-CoV-2 polymerase

Coronaviruses have evolved elaborate multisubunit machines to replicate and transcribe their genomes. Central to these machines are the RNA-dependent RNA polymerase subunit (nsp12) and its intimately associated cofactors (nsp7 and nsp8). We have used a high-throughput magnetic-tweezers approach to develop a mechanochemical description of this core polymerase. The core polymerase exists in at least three catalytically distinct conformations, one being kinetically consistent with incorporation of incorrect nucleotides. We provide the first evidence that an RdRp uses a thermal ratchet instead of a power stroke to transition from the pre- to post-translocated state. Ultra-stable magnetic tweezers enables the direct observation of coronavirus polymerase deep and long-lived backtrack that are strongly stimulated by secondary structure in the template. The framework presented here elucidates one of the most important structure-dynamics-function relationships in human health today, and will form the grounds for understanding the regulation of this complex.

Detection of genetic variation and base modifications at base-pair resolution on both DNA and RNA

Accurate decoding of nucleic acid variation is critical to understand the complexity and regulation of genome function. Here we use a single-molecule magnetic tweezer (MT) platform to identify sequence variation and map a range of important epigenetic base modifications with high sensitivity, specificity, and precision in the same single molecules of DNA or RNA. We have also developed a highly specific amplification-free CRISPR-Cas enrichment strategy to isolate genomic regions from native DNA. We demonstrate enrichment of DNA from both E. coli and the FMR1 5'UTR coming from cells derived from a Fragile X carrier. From these kilobase-length enriched molecules we could characterize the differential levels of adenine and cytosine base modifications on E. coli, and the repeat expansion length and methylation status of FMR1. Together these results demonstrate that our platform can detect a variety of genetic, epigenetic, and base modification changes concomitantly within the same single molecules.

GC-Content Dependence of Elastic and Overstretching Properties of DNA:RNA Hybrid Duplexes

DNA:RNA hybrid duplex plays important roles in various biological processes. Both its structural stability and interactions with proteins are highly sequence dependent. In this study, we utilize homebuilt optical tweezers to investigate how GC contents in the sequence influence the structural and mechanical properties of DNA:RNA hybrid by measuring its contour length, elasticities, and overstretching dynamics. Our results support that the DNA:RNA hybrid adopts a conformation between the A- and B-form helix, and the GC content does not affect its structural and elastic parameters obviously when varying from 40 to 60% before the overstretching transition of DNA:RNA hybrid occurs. In the overstretching transition, however, our study unravels significant heterogeneity and strong sequence dependence, suggesting the presence of a highly dynamic competition between the two processes, namely the S-form duplex formation (nonhysteretic) and the unpeeling (hysteretic). Analyzing the components left in DNA:RNA hybrid after the overstretching transition suggests that the RNA strand is more easily unpeeled than the DNA strand, whereas an increase in the GC content from 40 to 60% can significantly reduce unpeeling. Large hysteresis is observed between the stretching and relaxation processes, which is also quantitatively correlated with the percentage of unpeeling in the DNA:RNA duplex. Increasing in both the salt concentration and GC content can effectively reduce the hysteresis with the latter being more significant. Together, our study reveals that the mechanical properties of DNA:RNA hybrid duplexes are significantly different from double-stranded DNA and double-stranded RNA, and its overstretching behavior is highly sequence dependent. These results should be taken into account in the future studies on DNA:RNA-hybrid-related functional structures and motor proteins.