Double-stranded RNA under force and torque: similarities to and striking differences from double-stranded DNA

Multiplateau Force-Extension Curves of Long Double-Stranded DNA Molecules

When highly stretched, double-stranded DNA exhibits a plateau region in its force-extension curve. Using a bead-spring coarse-grained dynamic model based on a nonconvex potential, we predict that a long double-stranded DNA fragment made of several consecutive segments with substantially different plateau force values for each segment will exhibit multiple distinct plateau regions in the force-extension curve under physiologically relevant solvent conditions. For example, a long composite double-stranded (ds) DNA fragment consisting of two equal-length segments characterized by two different plateau force values, such as the poly(dA-dT)-poly(dG-dC) fragment, is predicted to exhibit two distinct plateau regions in its force-extension curve; a long composite dsDNA fragment consisting of three segments having three different plateau force values is predicted to have three distinct plateau regions. The formation of mixed states of slightly and highly stretched DNA, coexisting with macroscopically distinct phases of uniformly stretched DNA is also predicted. When one of the segments overstretches, the extensions of the segments can differ drastically. For example, for the poly(dA-dT)-poly(dG-dC) composite fragment, in the middle of the first plateau, 96.7% of the total extension of the fragment (relative to L x /L 0 ≈ 1.0) comes from the poly(dA-dT) segment, while only 3.3% of it comes from the poly(dG-dC) segment. The order of the segments has little effect on the force-extension curve or the distribution of conformational states. We speculate that the distinct structural states of stretched double-stranded DNA may have functional importance. For example, these states may modulate, in a sequence-dependent manner, the rate of double-stranded DNA processing by key cellular machines.

Effect of protein binding on the twist-stretch coupling of double-stranded RNA

The elasticities of RNAs are generally essential for their biological functions, and RNAs often become functional when interacting with their binding proteins. However, the effects of binding proteins on the elasticities of double-stranded (ds) RNAs, such as twist-stretch coupling, still remain little understood. Here, our extensive all-atom molecular dynamics simulations show that the twist-stretch coupling of dsRNAs can be reversed from positive to negative by their binding proteins. Our analyses revealed that such a reversing effect of binding proteins is attributed to the protein anchoring across the major groove of dsRNAs, which alters the dominating deformation pathway from a major-groove-mediated one to a helical-radius-mediated one through two base-pair parameters of slide and inclination. Meanwhile, the anchoring effect from binding proteins on dsRNAs is further ascribed to the strong electrostatic attractions between dsRNAs and the positively charged binding domain of the proteins.

Stretching and twisting of double-stranded RNA under forces: Unwinding mechanism and base-pair dependent elasticity

We used all-atom molecular dynamics simulations to investigate the mechanical response of double-stranded RNA (dsRNA) by applying various forces. We used the helical rise and helical twist, as well as a newly defined helical diameter, to characterize the stretching and twisting of dsRNA. The results indicate that dsRNA unwinds when stretched, accompanied by a linear increase in helical rise and helical diameter. Then, we utilized the normal modes, which are linear combinations of helical modes, to elucidate the underlying mechanism of dsRNA unwinding from an energetic perspective. On the other hand, we employed a stiffness matrix based on a rigid base pair model to examine the base-pair dependence of twist elasticity for dsRNA, as well as stretch elasticities with respect to the helical rise and helical diameter. The results show that the force induces variations in the local elasticities and their couplings of dsRNA, which are closely related to the distributions of base pairs. The mean stretch and twist elasticities can be considered as constants within the measurement uncertainties; however, their couplings demonstrate a slight linear dependency on applied force.

RNA aptamer-mediated RNA nanotechnology for potential treatment of cardiopulmonary diseases

Ribonucleic acid (RNA) aptamers are single-stranded RNAs that bind to target proteins or other molecules with high specificity and affinity, modulating biological functions through distinct mechanisms. These aptamers can act n as antagonists to block pathological interactions, agonists to activate signaling pathways, or delivery vehicles for therapeutic cargos such as siRNAs and miRNAs. The advances in RNA nanotechnology further enhances the versatility of RNA aptamers, offering scalable platforms for engineering. In this review, we have summarized recent developments in RNA aptamer-mediated RNA nanotechnology and provide an overview of its potential in treating cardiovascular and respiratory disorders, including atherosclerosis, acute coronary syndromes, heart failure, lung cancer, pulmonary hypertension, asthma, chronic obstructive pulmonary disease (COPD), acute lung injury, viral respiratory infections, and pulmonary fibrosis. By integrating aptamer technologies with innovative delivery systems, RNA aptamers hold the potential to revolutionize the treatment landscape for cardiopulmonary diseases.

Multi-plateau force-extension curves of long double-stranded DNA molecules

When highly stretched, double-stranded DNA exhibits a plateau region in its force-extension curve. Using a bead-spring coarse-grained dynamic model based on a non-convex potential, we predict that a long double-stranded DNA fragment made of several consecutive segments with substantially different plateau force values for each segment will exhibit multiple distinct plateau regions in the force-extension curve under physiologically relevant solvent conditions. For example, a long composite double-stranded (ds) DNA fragment consisting of two equal-length segments characterized by two different plateau force values, such as the poly(dA-dT)-poly(dG-dC) fragment, is predicted to exhibit two distinct plateau regions in its force-extension curve; a long composite dsDNA fragment consisting of three segments having three different plateau force values is predicted to have three distinct plateau regions. The formation of mixed states of slightly and highly stretched DNA, co-existing with macroscopically distinct phases of uniformly stretched DNA is also predicted. When one of the segments overstretches, the extensions of the segments can differ drastically. For example, for the poly(dA-dT)-poly(dG-dC) composite fragment, in the middle of the first plateau, 96.7 % of the total extension of the fragment (relative toL x / L 0 ≈ 1.0 ) comes from the poly(dA-dT) segment, while only 3.3 % of it comes from the poly(dG-dC) segment. The order of the segments has little effect on the force-extension curve or the distribution of conformational states. We speculate that the distinct structural states of stretched double-stranded DNA may have functional importance. For example, these states may modulate, in a sequence-dependent manner, the rate of double-stranded DNA processing by key cellular machines.

Temperature-induced swelling and unwinding of double-stranded DNA

We utilized all-atom molecular dynamics simulations to investigate the temperature-induced swelling and unwinding of double-stranded DNA (dsDNA). We adopted three helical parameters, specifically helical twist, helical rise, and diameter, to quantitatively describe the deformations and elastic properties associated with swelling and unwinding processes within an orthogonal cylindrical coordinate system. The results indicate that as temperature increases, dsDNA experiences a weak swelling accompanied by unwinding. This is associated with a slight increase in helical rise, while the helical diameter almost remains unchanged and the helical twist decreases. We evaluated all potential pathways for unwinding and elucidated that twist-diameter coupling drives the unwinding from an entropy perspective. On the other hand, we employed the rigid base pair model to examine the swelling and unwinding elasticities, with a focus on the stiffnesses of twist and diameter. The results suggest that the temperature induces variations in the local twist and diameter elasticities, as well as their couplings of dsDNA, which are closely related to the distance between the base pairs, attributed to its thermal fluctuations and correlations. The global twist elasticity reduces as the temperature rises; nonetheless, the global diameter elasticity and the twist-diameter coupling can be considered as constants, which indicate independence from the increasing temperature.

A hammerhead ribozyme selects mechanically stable conformations for catalysis against viral RNA

Ribozymes, widely found in prokaryotes and eukaryotes, target nucleic acids and can be engineered as biotechnical tools or for gene regulation or immune therapy. Among them, hammerhead is the smallest and best characterized ribozyme. However, the structure and biochemical data of ribozymes have been disagreed on, making the understanding of its catalysis mechanism a longstanding issue. Particularly, the role of conformational dynamics in ribozyme catalysis remains elusive. Here, we use single-molecule magnetic tweezers to reveal a concerted catalysis mechanism of mechanical conformational selection for a mini hammerhead ribozyme against a viral RNA sequence from the SARS-CoV-2. We identify a conformational set containing five mechanical conformers of the mini ribozyme, where magnesium ions select the active one. Our results are supported by molecular dynamics simulations. Our understanding of the RNA catalytic mechanism will be beneficial for ribozyme's biotechnological applications and as potential therapeutics against RNA viruses.

The impact of sequence periodicity on DNA mechanics: investigating the origin of A-tract's curvature

Periodic sequences in phase with DNA helical shape are prevalent in genomes due to their capacity to modulate DNA elasticity on a global scale. However, how this occurs is not well understood. We use all-atom molecular dynamics simulations on 40 bp DNA fragments to assess the effect of periodicity on bending, twisting, and stretch elasticity. We observe that DNA static curvature is the mechanical parameter most influenced by periodicity, with A-tract sequences having the greatest effect. A-tracts generate global curvature by bending in distinct directions (minor groove and backbones) that complement the bending of the rest of DNA, which predominantly is towards the major groove. Even if A-tracts are rigid at the local scale, these small bends integrate with the greater bends from the sequences between, producing an amplifying effect. As a result, our findings support a 'delocalized bend' model in which the A-tract operates as an 'adaptable mechanical part'. By understanding how global curvature emerges from local fluctuations, we reconcile previous contradictory theories and open an avenue for manipulating DNA mechanics through sequence design.

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

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

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

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

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

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

Insights into elastic properties of coarse-grained DNA models: q-stiffness of cgDNA vs cgDNA

Coarse-grained models have emerged as valuable tools to simulate long DNA molecules while maintaining computational efficiency. These models aim at preserving interactions among coarse-grained variables in a manner that mirrors the underlying atomistic description. We explore here a method for testing coarse-grained vs all-atom models using stiffness matrices in Fourier space (q-stiffnesses), which are particularly suited to probe DNA elasticity at different length scales. We focus on a class of coarse-grained rigid base DNA models known as cgDNA and its most recent version, cgDNA+. Our analysis shows that while cgDNA+ closely follows the q-stiffnesses of the all-atom model, the original cgDNA shows some deviations for twist and bending variables, which are rather strong in the q → 0 (long length scale) limit. The consequence is that while both cgDNA and cgDNA+ give a suitable description of local elastic behavior, the former misses some effects that manifest themselves at longer length scales. In particular, cgDNA performs poorly on twist stiffness, with a value much lower than expected for long DNA molecules. Conversely, the all-atom and cgDNA+ twist are strongly length scale dependent: DNA is torsionally soft at a few base pair distances but becomes more rigid at distances of a few dozen base pairs. Our analysis shows that the bending persistence length in all-atom and cgDNA+ is somewhat overestimated.

The origin of different bending stiffness between double-stranded RNA and DNA revealed by magnetic tweezers and simulations

The subtle differences in the chemical structures of double-stranded (ds) RNA and DNA lead to significant variations in their biological roles and medical implications, largely due to their distinct biophysical properties, such as bending stiffness. Although it is well known that A-form dsRNA is stiffer than B-form dsDNA under physiological salt conditions, the underlying cause of this difference remains unclear. In this study, we employ high-precision magnetic-tweezer experiments along with molecular dynamics simulations and reveal that the relative bending stiffness between dsRNA and dsDNA is primarily determined by the structure- and salt-concentration-dependent ion distribution around their helical structures. At near-physiological salt conditions, dsRNA shows a sparser ion distribution surrounding its phosphate groups compared to dsDNA, causing its greater stiffness. However, at very high monovalent salt concentrations, phosphate groups in both dsRNA and dsDNA become fully neutralized by excess ions, resulting in a similar intrinsic bending persistence length of approximately 39 nm. This similarity in intrinsic bending stiffness of dsRNA and dsDNA is coupled to the analogous fluctuations in their total groove widths and further coupled to the similar fluctuation of base-pair inclination, despite their distinct A-form and B-form helical structures.

Systematic Comparison of Atomistic Force Fields for the Mechanical Properties of Double-Stranded DNA

The response of double-stranded DNA to external mechanical stress plays a central role in its interactions with the protein machinery in the cell. Modern atomistic force fields have been shown to provide highly accurate predictions for the fine structural features of the duplex. In contrast, and despite their pivotal function, less attention has been devoted to the accuracy of the prediction of the elastic parameters. Several reports have addressed the flexibility of double-stranded DNA via all-atom molecular dynamics, yet the collected information is insufficient to have a clear understanding of the relative performance of the various force fields. In this work, we fill this gap by performing a systematic study in which several systems, characterized by different sequence contexts, are simulated with the most popular force fields within the AMBER family, bcs1 and OL15, as well as with CHARMM36. Analysis of our results, together with their comparison with previous work focused on bsc0, allows us to unveil the differences in the predicted rigidity between the newest force fields and suggests a roadmap to test their performance against experiments. In the case of the stretch modulus, we reconcile these differences, showing that a single mapping between sequence-dependent conformation and elasticity via the crookedness parameter captures simultaneously the results of all force fields, supporting the key role of crookedness in the mechanical response of double-stranded DNA.

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

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

Temperature-dependent elasticity of DNA, RNA, and hybrid double helices

Nucleic acid double helices in their DNA, RNA, and DNA-RNA hybrid form play a fundamental role in biology and are main building blocks of artificial nanostructures, but how their properties depend on temperature remains poorly understood. Here, we report thermal dependence of dynamic bending persistence length, twist rigidity, stretch modulus, and twist-stretch coupling for DNA, RNA, and hybrid duplexes between 7°C and 47°C. The results are based on all-atom molecular dynamics simulations using different force field parameterizations. We first demonstrate that unrestrained molecular dynamics can reproduce experimentally known mechanical properties of the duplexes at room temperature. Beyond experimentally known features, we also infer the twist rigidity and twist-stretch coupling of the hybrid duplex. As for the temperature dependence, we found that increasing temperature softens all the duplexes with respect to bending, twisting, and stretching. The relative decrease of the stretch moduli is 0.003-0.004/°C, similar for all the duplex variants despite their very different stretching stiffness, whereas RNA twist stiffness decreases by 0.003/°C, and smaller values are found for the other elastic moduli. The twist-stretch couplings are nearly unaffected by temperature. The stretching, bending, and twisting stiffness all include an important entropic component. Relation of our results to the two-state model of DNA flexibility is discussed. Our work provides temperature-dependent elasticity of nucleic acid duplexes at the microsecond scale relevant for initial stages of protein binding.

Engineering ssRNA tile filaments for (dis)assembly and membrane binding

Cytoskeletal protein filaments such as actin and microtubules confer mechanical support to cells and facilitate many cellular functions such as motility and division. Recent years have witnessed the development of a variety of molecular scaffolds that mimic such filaments. Indeed, filaments that are programmable and compatible with biological systems may prove useful in studying or substituting such proteins. Here, we explore the use of ssRNA tiles to build and modify filaments in vitro. We engineer a number of functionalities that are crucial to the function of natural proteins filaments into the ssRNA tiles, including the abilities to assemble or disassemble filaments, to tune the filament stiffness, to induce membrane binding, and to bind proteins. This work paves the way for building dynamic cytoskeleton-mimicking systems made out of rationally designed ssRNA tiles that can be transcribed in natural or synthetic cells.

Liquid crystalline inverted lipid phases encapsulating siRNA enhance lipid nanoparticle mediated transfection

Efficient cytosolic delivery of RNA molecules remains a formidable barrier for RNA therapeutic strategies. Lipid nanoparticles (LNPs) serve as state-of-the-art carriers that can deliver RNA molecules intracellularly, as exemplified by the recent implementation of several vaccines against SARS-CoV-2. Using a bottom-up rational design approach, we assemble LNPs that contain programmable lipid phases encapsulating small interfering RNA (siRNA). A combination of cryogenic transmission electron microscopy, cryogenic electron tomography and small-angle X-ray scattering reveals that we can form inverse hexagonal structures, which are present in a liquid crystalline nature within the LNP core. Comparison with lamellar LNPs reveals that the presence of inverse hexagonal phases enhances the intracellular silencing efficiency over lamellar structures. We then demonstrate that lamellar LNPs exhibit an in situ transition from a lamellar to inverse hexagonal phase upon interaction with anionic membranes, whereas LNPs containing pre-programmed liquid crystalline hexagonal phases bypass this transition for a more efficient one-step delivery mechanism, explaining the increased silencing effect. This rational design of LNPs with defined lipid structures aids in the understanding of the nano-bio interface and adds substantial value for LNP design, optimization and use.

Temperature-Dependent Twist of Double-Stranded RNA Probed by Magnetic Tweezer Experiments and Molecular Dynamics Simulations

RNA plays critical roles in the transmission and regulation of genetic information and is increasingly used in biomedical and biotechnological applications. Functional RNAs contain extended double-stranded regions, and the structure of double-stranded RNA (dsRNA) has been revealed at high resolution. However, the dependence of the properties of the RNA double helix on environmental effects, notably temperature, is still poorly understood. Here, we use single-molecule magnetic tweezer measurements to determine the dependence of the dsRNA twist on temperature. We find that dsRNA unwinds with increasing temperature, even more than DNA, with ΔTwRNA = -14.4 ± 0.7°/(°C·kbp), compared to ΔTwDNA = -11.0 ± 1.2°/(°C·kbp). All-atom molecular dynamics (MD) simulations using a range of nucleic acid force fields, ion parameters, and water models correctly predict that dsRNA unwinds with rising temperature but significantly underestimate the magnitude of the effect. These MD data, together with additional MD simulations involving DNA and DNA-RNA hybrid duplexes, reveal a linear correlation between the twist temperature decrease and the helical rise, in line with DNA but at variance with RNA experimental data. We speculate that this discrepancy might be caused by some unknown bias in the RNA force fields tested or by as yet undiscovered transient alternative structures in the RNA duplex. Our results provide a baseline to model more complex RNA assemblies and to test and develop new parametrizations for RNA simulations. They may also inspire physical models of the temperature-dependent dsRNA structure.

Optical Tweezers to Study Viruses

A virus is a complex molecular machine that propagates by channeling its genetic information from cell to cell. Unlike macroscopic engines, it operates in a nanoscopic world under continuous thermal agitation. Viruses have developed efficient passive and active strategies to pack and release nucleic acids. Some aspects of the dynamic behavior of viruses and their substrates can be studied using structural and biochemical techniques. By the turn of the millennium, physical techniques have been applied to dynamic studies of viruses in which their intrinsic mechanical activity can be measured directly. Optical tweezers are a technology that can be used to measure the force, torque, and strain produced by molecular motors, as a function of time and at the single-molecule level. Thanks to this technique, some bacteriophages are now known to be powerful nanomachines; they exert force in the piconewton range and their motors work in a highly coordinated fashion for packaging the viral nucleic acid genome. Nucleic acids, whose elasticity and condensation behavior are inherently coupled to the viral packaging mechanisms, virion assembly, and virion-cell interactions are also amenable to examination with optical tweezers. In this chapter, we provide a comprehensive analysis of this laser-based tool, its combination with imaging methods, and its application to the study of viruses and viral molecules.

Twist-stretch relations in nucleic acids

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

Mechanical codes of chemical-scale specificity in DNA motifs

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

A microfluidic electrophoretic dual dynamic staining method for the identification and relative quantitation of dsRNA contaminants in mRNA vaccines

mRNA vaccines (i.e., COVID-19 vaccine) offer various advantages over traditional vaccines in preventing and reducing disease and shortening the time between pathogen discovery and vaccine creation. Production of mRNA vaccines results in several nucleic acid and enzymatic by-products, most of which can be detected and removed; however, double-stranded RNA (dsRNA) contaminants pose a particular challenge. Current purification and detection platforms for dsRNA vary in effectiveness, with problems in scalability for mass mRNA vaccine production. Effectively detecting dsRNA is crucial in ensuring the safety and efficacy of the vaccines, as these strands can cause autoimmune reactions with length-symptom dependency and enhance mRNA degradation. We present a new microfluidics method to rapidly identify and quantify dsRNA fragments in mRNA samples. Our innovation exploits the differences in the dynamic staining behavior between mRNA and dsRNA molecules to detect dsRNA contaminants in a high throughput approach. The limit of detection of the system for dsRNA was estimated to be between 17.7-76.6 pg μL-1 with a maximum loading capacity of mRNA of 12.99 ng μL-1. Based on these estimated values, our method allows for the detection of dsRNA contaminants present in percentages as low as 0.14-0.59% compared to the total mRNA concentration. Here, we discuss the molecular mechanism of the dynamic staining behavior of dsRNA and mRNA for two different stains. We believe our method will accelerate the mRNA vaccine development from initial development to quality control workflows.

The organization of double-stranded RNA in the chikungunya virus replication organelle

Alphaviruses are mosquito-borne, positive-sense single-stranded RNA viruses. Amongst the alphaviruses, chikungunya virus is notable as a large source of human illness, especially in tropical and subtropical regions. When they invade a cell, alphaviruses generate dedicated organelles for viral genome replication, so-called spherules. Spherules form as outward-facing buds at the plasma membrane, and it has recently been shown that the thin membrane neck that connects this membrane bud with the cytoplasm is guarded by a two-megadalton protein complex that contains all the enzymatic functions necessary for RNA replication. The lumen of the spherules contains a single copy of the negative-strand template RNA, present in a duplex with newly synthesized positive-sense RNA. Less is known about the organization of this double-stranded RNA as compared to the protein components of the spherule. Here, we analyzed cryo-electron tomograms of chikungunya virus spherules in terms of the organization of the double-stranded RNA replication intermediate. We find that the double-stranded RNA has a shortened apparent persistence length as compared to unconstrained double-stranded RNA. Around half of the genome is present in either of five conformations identified by subtomogram classification, each representing a relatively straight segment of ~25-32 nm. Finally, the RNA occupies the spherule lumen at a homogeneous density, but has a preferred orientation to be perpendicular to a vector pointing from the membrane neck towards the spherule center. Taken together, this analysis lays another piece of the puzzle of the highly coordinated alphavirus genome replication.

High-yield ligation-free assembly of DNA constructs with nucleosome positioning sequence repeats for single-molecule manipulation assays

Force and torque spectroscopy have provided unprecedented insights into the mechanical properties, conformational transitions, and dynamics of DNA and DNA-protein complexes, notably nucleosomes. Reliable single-molecule manipulation measurements require, however, specific and stable attachment chemistries to tether the molecules of interest. Here, we present a functionalization strategy for DNA that enables high-yield production of constructs for torsionally constrained and very stable attachment. The method is based on two subsequent PCRs: first ∼380 bp long DNA strands are generated that contain multiple labels, which are used as "megaprimers" in a second PCR to generate ∼kbp long double-stranded DNA constructs with multiple labels at the respective ends. To achieve high-force stability, we use dibenzocyclooctyne-based click chemistry for covalent attachment to the surface and biotin-streptavidin coupling to the bead. The resulting tethers are torsionally constrained and extremely stable under load, with an average lifetime of 70 ± 3 h at 45 pN. The high yield of the approach enables nucleosome reconstitution by salt dialysis on the functionalized DNA, and we demonstrate proof-of-concept measurements on nucleosome assembly statistics and inner turn unwrapping under force. We anticipate that our approach will facilitate a range of studies of DNA interactions and nucleoprotein complexes under forces and torques.

Revolving hexameric ATPases as asymmetric motors to translocate double-stranded DNA genome along one strand

DsDNA translocation through nanoscale pores is generally accomplished by ATPase biomotors. The discovery of the revolving dsDNA translocation mechanism, as opposed to rotation, in bacteriophage phi29 elucidated how ATPase motors move dsDNA. Revolution-driven, hexameric dsDNA motors have been reported in herpesvirus, bacterial FtsK, Streptomyces TraB, and T7 phage. This review explores the common relationship between their structure and mechanisms. Commonalities include moving along the 5'→3' strand, inchworm sequential action leading to an asymmetrical structure, channel chirality, channel size, and 3-step channel gating for controlling motion direction. The revolving mechanism and contact with one of the dsDNA strands addresses the historic controversy of dsDNA packaging using nicked, gapped, hybrid, or chemically modified DNA. These controversies surrounding dsDNA packaging activity using modified materials can be answered by whether the modification was introduced into the 3'→5' or 5'→3' strand. Perspectives concerning solutions to the controversy of motor structure and stoichiometry are also discussed.

Filament assembly underpins the double-stranded DNA specificity of AIM2-like receptors

Upon sensing cytosolic- and/or viral double-stranded (ds)DNA, absent-in-melanoma-2 (AIM2)-like-receptors (ALRs) assemble into filamentous signaling platforms to initiate inflammatory responses. The versatile yet critical roles of ALRs in host innate defense are increasingly appreciated; however, the mechanisms by which AIM2 and its related IFI16 specifically recognize dsDNA over other nucleic acids remain poorly understood (i.e. single-stranded (ss)DNA, dsRNA, ssRNA and DNA:RNA hybrid). Here, we find that although AIM2 can interact with various nucleic acids, it preferentially binds to and assembles filaments faster on dsDNA in a duplex length-dependent manner. Moreover, AIM2 oligomers assembled on nucleic acids other than dsDNA not only display less ordered filamentous structures, but also fail to induce the polymerization of downstream ASC. Likewise, although showing broader nucleic acid selectivity than AIM2, IFI16 binds to and oligomerizes most readily on dsDNA in a duplex length-dependent manner. Nevertheless, IFI16 fails to form filaments on single-stranded nucleic acids and does not accelerate the polymerization of ASC regardless of bound nucleic acids. Together, we reveal that filament assembly is integral to nucleic acid distinction by ALRs.

Force-dependent elasticity of nucleic acids

The functioning of double-stranded (ds) nucleic acids (NAs) in cellular processes is strongly mediated by their elastic response. These processes involve proteins that interact with dsDNA or dsRNA and distort their structures. The perturbation of the elasticity of NAs arising from these deformations is not properly considered by most theoretical frameworks. In this work, we introduce a novel method to assess the impact of mechanical stress on the elastic response of dsDNA and dsRNA through the analysis of the fluctuations of the double helix. Application of this approach to atomistic simulations reveals qualitative differences in the force dependence of the mechanical properties of dsDNA with respect to those of dsRNA, which we relate to structural features of these molecules by means of physically-sound minimalistic models.

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

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

Nucleic acid nanostructures for in vivo applications: The influence of morphology on biological fate

The development of programmable biomaterials for use in nanofabrication represents a major advance for the future of biomedicine and diagnostics. Recent advances in structural nanotechnology using nucleic acids have resulted in dramatic progress in our understanding of nucleic acid-based nanostructures (NANs) for use in biological applications. As the NANs become more architecturally and functionally diverse to accommodate introduction into living systems, there is a need to understand how critical design features can be controlled to impart desired performance in vivo. In this review, we survey the range of nucleic acid materials utilized as structural building blocks (DNA, RNA, and xenonucleic acids), the diversity of geometries for nanofabrication, and the strategies to functionalize these complexes. We include an assessment of the available and emerging characterization tools used to evaluate the physical, mechanical, physiochemical, and biological properties of NANs in vitro. Finally, the current understanding of the obstacles encountered along the in vivo journey is contextualized to demonstrate how morphological features of NANs influence their biological fates. We envision that this summary will aid researchers in the designing novel NAN morphologies, guide characterization efforts, and design of experiments and spark interdisciplinary collaborations to fuel advancements in programmable platforms for biological applications.

Maturation of siRNA by strand separation: Steered molecular dynamics study

RNA interference, particularly siRNA induced gene silencing is becoming an important avenue of modern therapeutics. The siRNA is delivered to the cells as short double helical RNA which becomes single stranded for forming the RISC complex. Significant experimental evidence is available for most of the steps except the process of the separation of the two strands. We have attempted to understand the pathway for double stranded siRNA (dsRNA) to single stranded (ssRNA) molecules using steered molecular dynamics simulations. As the process is completely unexplored we have applied force from all possible directions restraining all possible residues to convert dsRNA to ssRNA. We found pulling one strand along the helical axis direction restraining the far end of the other strand demands excessive force for ssRNA formation. Pulling a central residue of one strand, in a direction perpendicular to the helix axis, while keeping the base paired residue fixed requires intermediate force for strand separation. Moreover, we found that in this process the force requirement is quite high for the first bubble formation (nucleation energy) and the bubble propagation energies are quite small. We believe the success rate of the design of siRNA sequences for gene silencing may increase if this mechanistic knowledge is utilized for such a design process.Communicated by Ramaswamy H. Sarma.

Biophysical Modeling of SARS-CoV-2 Assembly: Genome Condensation and Budding

The COVID-19 pandemic caused by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has spurred unprecedented and concerted worldwide research to curtail and eradicate this pathogen. SARS-CoV-2 has four structural proteins: Envelope (E), Membrane (M), Nucleocapsid (N), and Spike (S), which self-assemble along with its RNA into the infectious virus by budding from intracellular lipid membranes. In this paper, we develop a model to explore the mechanisms of RNA condensation by structural proteins, protein oligomerization and cellular membrane-protein interactions that control the budding process and the ultimate virus structure. Using molecular dynamics simulations, we have deciphered how the positively charged N proteins interact and condense the very long genomic RNA resulting in its packaging by a lipid envelope decorated with structural proteins inside a host cell. Furthermore, considering the length of RNA and the size of the virus, we find that the intrinsic curvature of M proteins is essential for virus budding. While most current research has focused on the S protein, which is responsible for viral entry, and it has been motivated by the need to develop efficacious vaccines, the development of resistance through mutations in this crucial protein makes it essential to elucidate the details of the viral life cycle to identify other drug targets for future therapy. Our simulations will provide insight into the viral life cycle through the assembly of viral particles de novo and potentially identify therapeutic targets for future drug development.

Nucleic Acid Thermodynamics Derived from Mechanical Unzipping Experiments

Force-spectroscopy techniques have led to significant progress in studying the physicochemical properties of biomolecules that are not accessible in bulk assays. The application of piconewton forces with laser optical tweezers to single nucleic acids has permitted the characterization of molecular thermodynamics and kinetics with unprecedented accuracy. Some examples are the hybridization reaction between complementary strands in DNA and the folding of secondary, tertiary, and other heterogeneous structures, such as intermediate and misfolded states in RNA. Here we review the results obtained in our lab on deriving the nearest-neighbor free energy parameters in DNA and RNA duplexes from mechanical unzipping experiments. Remarkable nonequilibrium effects are also observed, such as the large irreversibility of RNA unzipping and the formation of non-specific secondary structures in single-stranded DNA. These features originate from forming stem-loop structures along the single strands of the nucleic acid. The recently introduced barrier energy landscape model quantifies kinetic trapping effects due to stem-loops being applicable to both RNA and DNA. The barrier energy landscape model contains the essential features to explain the many behaviors observed in heterogeneous nucleic-acid folding.

Mechanical properties of nucleic acids and the non-local twistable wormlike chain model

Mechanical properties of nucleic acids play an important role in many biological processes that often involve physical deformations of these molecules. At sufficiently long length scales (say, above ∼20-30 base pairs), the mechanics of DNA and RNA double helices is described by a homogeneous Twistable Wormlike Chain (TWLC), a semiflexible polymer model characterized by twist and bending stiffnesses. At shorter scales, this model breaks down for two reasons: the elastic properties become sequence-dependent and the mechanical deformations at distal sites get coupled. We discuss in this paper the origin of the latter effect using the framework of a non-local Twistable Wormlike Chain (nlTWLC). We show, by comparing all-atom simulations data for DNA and RNA double helices, that the non-local couplings are of very similar nature in these two molecules: couplings between distal sites are strong for tilt and twist degrees of freedom and weak for roll. We introduce and analyze a simple double-stranded polymer model that clarifies the origin of this universal distal couplings behavior. In this model, referred to as the ladder model, a nlTWLC description emerges from the coarsening of local (atomic) degrees of freedom into angular variables that describe the twist and bending of the molecule. Different from its local counterpart, the nlTWLC is characterized by a length-scale-dependent elasticity. Our analysis predicts that nucleic acids are mechanically softer at the scale of a few base pairs and are asymptotically stiffer at longer length scales, a behavior that matches experimental data.

Twisting DNA by salt

The structure and properties of DNA depend on the environment, in particular the ion atmosphere. Here, we investigate how DNA twist -one of the central properties of DNA- changes with concentration and identity of the surrounding ions. To resolve how cations influence the twist, we combine single-molecule magnetic tweezer experiments and extensive all-atom molecular dynamics simulations. Two interconnected trends are observed for monovalent alkali and divalent alkaline earth cations. First, DNA twist increases monotonously with increasing concentration for all ions investigated. Second, for a given salt concentration, DNA twist strongly depends on cation identity. At 100 mM concentration, DNA twist increases as Na+ < K+ < Rb+ < Ba2+ < Li+ ≈ Cs+ < Sr2+ < Mg2+ < Ca2+. Our molecular dynamics simulations reveal that preferential binding of the cations to the DNA backbone or the nucleobases has opposing effects on DNA twist and provides the microscopic explanation of the observed ion specificity. However, the simulations also reveal shortcomings of existing force field parameters for Cs+ and Sr2+. The comprehensive view gained from our combined approach provides a foundation for understanding and predicting cation-induced structural changes both in nature and in DNA nanotechnology.

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.

Accurate Sequence-Dependent Coarse-Grained Model for Conformational and Elastic Properties of Double-Stranded DNA

We introduce MADna, a sequence-dependent coarse-grained model of double-stranded DNA (dsDNA), where each nucleotide is described by three beads localized at the sugar, at the base moiety, and at the phosphate group, respectively. The sequence dependence is included by considering a step-dependent parametrization of the bonded interactions, which are tuned in order to reproduce the values of key observables obtained from exhaustive atomistic simulations from the literature. The predictions of the model are benchmarked against an independent set of all-atom simulations, showing that it captures with high fidelity the sequence dependence of conformational and elastic features beyond the single step considered in its formulation. A remarkably good agreement with experiments is found for both sequence-averaged and sequence-dependent conformational and elastic features, including the stretching and torsion moduli, the twist–stretch and twist–bend couplings, the persistence length, and the helical pitch. Overall, for the inspected quantities, the model has a precision comparable to atomistic simulations, hence providing a reliable coarse-grained description for the rationalization of single-molecule experiments and the study of cellular processes involving dsDNA. Owing to the simplicity of its formulation, MADna can be straightforwardly included in common simulation engines. Particularly, an implementation of the model in LAMMPS is made available on an online repository to ease its usage within the DNA research community.

Nanoscale structures and mechanics of peptide nucleic acids

Peptide nucleic acids (PNAs) are charge-neutral polyamide oligomers having extremely favorable thermal stability and high affinity to cell membranes when coupled with cationic cell-penetrating peptides (CPPs), as well as the encouraging antisense and antigene activity in cell-free systems. The study of the mechanical properties of short PNA molecules is rare both in experiments and theoretical calculations. Here, we studied the microscopic structures and elastic properties; namely, persistence length, stretch modulus, twist-stretch coupling, and structural crookedness of double-stranded PNA (dsPNA) and their hybrid derivatives using all-atom MD simulation and compared them with those of double-stranded DNA (dsDNA) and double-stranded RNA (dsRNA). The stretch modulus of the dsPNA is found to be ∼160 pN, an order of magnitude lower than that of dsDNA and smaller than dsRNA, respectively. Similarly, the persistence length of dsPNA is found to be ∼35 nm, significantly smaller than those of dsDNA and dsRNA. The PNA-DNA and PNA-RNA hybrid duplexes have elastic properties lying between that of dsPNA and dsDNA/dsRNA. We argue that the neutral backbones of the PNA make it less stiff than dsDNA and dsRNA molecules. Measurement of structural crookedness and principal component analysis additionally support the bending flexibility of dsPNA. Detailed analysis of the helical-rise coupled to helical-twist indicates that the PNA-DNA hybrid over-winds like dsDNA, while PNA-PNA and PNA-RNA unwind like dsRNA upon stretching. Because of the highly flexible nature of PNA, it can bind other biomolecules by adopting a wide range of conformations and is believed to be crucial for future nanobiotechnology research studies.

Diameter Dependent Melting and Softening of dsDNA Under Cylindrical Confinement

Carbon nanotubes (CNTs) are considered promising candidates for biomolecular confinement, including DNA encapsulation for gene delivery. Threshold values of diameters have been reported for double-stranded DNA (dsDNA) encapsulation inside CNTs. We have performed all-atom molecular dynamics (MD) simulations of dsDNAs confined inside single-walled CNTs (SWCNTs) at the physiologically relevant temperature of 300 K. We found that the dsDNA can be confined without being denatured only when the diameter of the SWCNT exceeds a threshold value. Below this threshold diameter, the dsDNA gets denatured and melts even at the temperature of 300 K. Our simulations using SWCNTs with chirality indices (20,20) to (30,30) at 300 K found the critical diameter to be 3.25 nm (corresponding to (24,24) chirality). Analyses of the hydrogen bonds (H-bonds), Van der Walls (VdW) energy, and other inter-base interactions show drastic reduction in the number of H-bonds, VdW energy, and electrostatic energies between the bases of dsDNA when it is confined in narrower SWCNTs (up to diameter of 3.12 nm). On the other hand, the higher interaction energy between the dsDNA and the SWCNT surface in narrower SWCNTs assists in the melting of the dsDNA. Electrostatic mapping and hydration status analyses show that the dsDNA is not adequately hydrated and the counter ion distribution is not uniform below the critical diameter of the SWCNT. As properly hydrated counter ions provide stability to the dsDNA, we infer that the inappropriate hydration of counter ions and their non-uniform distribution around the dsDNA cause the melting of the dsDNA inside SWCNTs of diameter below the critical value of 3.25 nm. For confined dsDNAs that do not get denatured, we computed their elastic properties. The persistence length of dsDNA was found to increase by a factor of about two and the torsional stiffness by a factor of 1.5 for confinement inside SWCNTs of diameters up to 3.79 nm, the stretch modulus also following nearly the same trend. Interestingly, for higher diameters of SWCNT, 3.79 nm and above, the dsDNA becomes more flexible, demonstrating that the mechanical properties of the dsDNA under cylindrical confinement depend non-monotonically on the confinement diameter.

Twist-diameter coupling drives DNA twist changes with salt and temperature

DNA deformations upon environmental changes, e.g., salt and temperature, play crucial roles in many biological processes and material applications. Here, our magnetic tweezers experiments observed that the increase in NaCl, KCl, or RbCl concentration leads to substantial DNA overwinding. Our simulations and theoretical calculation quantitatively explain the salt-induced twist change through the mechanism: More salt enhances the screening of interstrand electrostatic repulsion and hence reduces DNA diameter, which is transduced to twist increase through twist-diameter coupling. We determined that the coupling constant is 4.5 ± 0.8 k_BT/(degrees∙nm) for one base pair. The coupling comes from the restraint of the contour length of DNA backbone. On the basis of this coupling constant and diameter-dependent DNA conformational entropy, we predict the temperature dependence of DNA twist Δω_bp/ΔT ≈ -0.01 degree/°C, which agrees with our and previous experimental results. Our analysis suggests that twist-diameter coupling is a common driving force for salt- and temperature-induced DNA twist changes.

Multivalent Cations Reverse the Twist-Stretch Coupling of RNA

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

Modeling multiple duplex DNA attachments in a force-extension experiment

Optical tweezers-based DNA stretching often relies on tethering a single end-activated DNA molecule between optically manipulated end-binding beads. Measurement success can depend on DNA concentration. At lower DNA concentrations tethering is less common, and many trials may be required to observe a single-molecule stretch. At higher DNA concentrations tethering is more common; however, the resulting force-extensions observed are more complex and may vary from measurement to measurement. Typically these more complex results are attributed to the formation of multiple tethers between the beads; however, to date there does not appear to have been a critical examination of this hypothesis or the potential usefulness of such data. Here we examine stretches at a higher DNA concentration and use analysis and simulation to show how the more complex force-extensions observed can be understood in terms of multiple DNA attachments.

Nanoscale structures and mechanics of peptide nucleic acids

Peptide nucleic acids are charge-neutral polyamide oligomers with extremely flexible backbones that have a strong affinity for hybridization with complementary DNA or RNA, as well as encouraging antisense and antigene activity in cell-free systems.

Cargo Release from Nonenveloped Viruses and Virus-like Nanoparticles: Capsid Rupture or Pore Formation

Virus-like nanoparticles are protein shells similar to wild-type viruses, and both aim to deliver their content into a cell. Unfortunately, the release mechanism of their cargo/genome remains elusive. Pores on the symmetry axes were proposed to enable the slow release of the viral genome. In contrast, cryo-EM images showed that capsids of nonenveloped RNA viruses can crack open and rapidly release the genome. We combined in vitro cryo-EM observations of the genome release of three viruses with coarse-grained simulations of generic virus-like nanoparticles to investigate the cargo/genome release pathways. Simulations provided details on both slow and rapid release pathways, including the success rates of individual releases. Moreover, the simulated structures from the rapid release pathway were in agreement with the experiment. Slow release occurred when interactions between capsid subunits were long-ranged, and the cargo/genome was noncompact. In contrast, rapid release was preferred when the interaction range was short and/or the cargo/genome was compact. These findings indicate a design strategy of virus-like nanoparticles for drug delivery.

RNA Secondary Structures Regulate Adsorption of Fragments onto Flat Substrates

RNA is a functionally rich molecule with multilevel, hierarchical structures whose role in the adsorption to molecular substrates is only beginning to be elucidated. Here, we introduce a multiscale simulation approach that combines a tractable coarse-grained RNA structural model with an interaction potential of a structureless flat adsorbing substrate. Within this approach, we study the specific role of stem-hairpin and multibranch RNA secondary structure motifs on its adsorption phenomenology. Our findings identify a dual regime of adsorption for short RNA fragments with and without the secondary structure and underline the adsorption efficiency in both cases as a function of the surface interaction strength. The observed behavior results from an interplay between the number of contacts formed at the surface and the conformational entropy of the RNA molecule. The adsorption phenomenology of RNA seems to persist also for much longer RNAs as qualitatively observed by comparing the trends of our simulations with a theoretical approach based on an ideal semiflexible polymer chain.

Tumuc1: A New Accurate DNA Force Field Consistent with High-Level Quantum Chemistry

An accurate molecular mechanics force field forms the basis of Molecular Dynamics simulations to obtain a realistic view of the structure and dynamics of biomolecules such as DNA. Although frequently updated to improve agreement with available experimental data, DNA force fields still rely in part on parameters introduced more than 20 years ago. We have developed an entirely new DNA force field, Tumuc1, derived from quantum mechanical calculations to obtain a consistent set of bonded parameters and partial atomic charges. The performance of the force field was extensively tested on a variety of DNA molecules. It excels in accuracy of B-DNA simulations but also performs very well on other types of DNA structures and structure formation processes such as hairpin folding, duplex formation, and dynamics of DNA-protein complexes. It can complement existing force fields in order to provide an increasingly accurate description of the structure and dynamics of DNA during simulation studies.

The structure of nontypeable Haemophilus influenzae SapA in a closed conformation reveals a constricted ligand-binding cavity and a novel RNA binding motif

Nontypeable Haemophilus influenzae (NTHi) is a significant pathogen in respiratory disease and otitis media. Important for NTHi survival, colonization and persistence in vivo is the Sap (sensitivity to antimicrobial peptides) ABC transporter system. Current models propose a direct role for Sap in heme and antimicrobial peptide (AMP) transport. Here, the crystal structure of SapA, the periplasmic component of Sap, in a closed, ligand bound conformation, is presented. Phylogenetic and cavity volume analysis predicts that the small, hydrophobic SapA central ligand binding cavity is most likely occupied by a hydrophobic di- or tri- peptide. The cavity is of insufficient volume to accommodate heme or folded AMPs. Crystal structures of SapA have identified surface interactions with heme and dsRNA. Heme binds SapA weakly (Kd 282 μM) through a surface exposed histidine, while the dsRNA is coordinated via residues which constitute part of a conserved motif (estimated Kd 4.4 μM). The RNA affinity falls within the range observed for characterized RNA/protein complexes. Overall, we describe in molecular-detail the interactions of SapA with heme and dsRNA and propose a role for SapA in the transport of di- or tri-peptides.