Computing the Elastic Mechanical Properties of Rodlike DNA Nanostructures

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.

A lumen-tunable triangular DNA nanopore for molecular sensing and cross-membrane transport

Synthetic membrane nanopores made of DNA are promising systems to sense and control molecular transport in biosensing, sequencing, and synthetic cells. Lumen-tunable nanopore like the natural ion channels and systematically increasing the lumen size have become long-standing desires in developing nanopores. Here, we design a triangular DNA nanopore with a large tunable lumen. It allows in-situ transition from expanded state to contracted state without changing its stable triangular shape, and vice versa, in which specific DNA bindings as stimuli mechanically pinch and release the three corners of the triangular frame. Transmission electron microscopy images and molecular dynamics simulations illustrate the stable architectures and the high shape retention. Single-channel current recordings and fluorescence influx studies demonstrate the low-noise repeatable readouts and the controllable cross-membrane macromolecular transport. We envision that the proposed DNA nanopores could offer powerful tools in molecular sensing, drug delivery, and the creation of synthetic cells.

Prediction and Control in DNA Nanotechnology

DNA nanotechnology is a rapidly developing field that uses DNA as a building material for nanoscale structures. Key to the field's development has been the ability to accurately describe the behavior of DNA nanostructures using simulations and other modeling techniques. In this Review, we present various aspects of prediction and control in DNA nanotechnology, including the various scales of molecular simulation, statistical mechanics, kinetic modeling, continuum mechanics, and other prediction methods. We also address the current uses of artificial intelligence and machine learning in DNA nanotechnology. We discuss how experiments and modeling are synergistically combined to provide control over device behavior, allowing scientists to design molecular structures and dynamic devices with confidence that they will function as intended. Finally, we identify processes and scenarios where DNA nanotechnology lacks sufficient prediction ability and suggest possible solutions to these weak areas.

Correlated Hybrid DNA Structures Explored by the oxDNA Model

Thermodynamically, perfect DNA hybridization can be formed between probes and their corresponding targets due to the favorable energy. However, this is not the case dynamically. Here, we use molecular dynamics (MD) simulations based on the oxDNA model to investigate the process of DNA microarray hybridization. In general, correlated hybrid DNA structures are formed, including one probe associated with several targets as well as one target hybrid with multiple probes leading to the target-mediated hybridization. The formation of these two types of correlated structures largely depends on the surface coverage of the DNA microarray. Moreover, DNA sequence, DNA length, and spacer length have an impact on the structural formation. Our findings shed light on the dynamics of DNA hybridization, which is important for the application of DNA microarray.

Conformational dynamics and mechanical properties of biomimetic RNA, DNA, and RNA-DNA hybrid nanotubes: an atomistic molecular dynamics study

With the nanotechnology boom, artificially designed nucleic acid nanotubes have aroused interest due to their practical applications in nanorobotics, vaccine design, membrane channels, drug delivery, and force sensing. In this paper, computational study was performed to investigate the structural dynamics and mechanical properties of RNA nanotubes (RNTs), DNA nanotubes (DNTs), and RNA-DNA hybrid nanotubes (RDHNTs). So far, the structural and mechanical properties of RDHNTs have not been examined in experiments or theoretical calculations, and there is limited knowledge regarding these properties for RNTs. Here, the simulations were carried out using the equilibrium molecular dynamics (MD) and steered molecular dynamics (SMD) approaches. Using in-house scripting, we modeled hexagonal nanotubes composed of six double-stranded molecules connected by four-way Holliday junctions. Classical MD analyses were performed on the collected trajectory data to investigate structural properties. Analyses of the microscopic structural parameters of RDHNT indicated a structural transition from the A-form to a conformation between the A- and B-forms, which may be attributable to the increased rigidity of RNA scaffolds compared to DNA staples. Comprehensive research on the elastic mechanical properties was also conducted based on spontaneous thermal fluctuations of nanotubes and employing the equipartition theorem. The Young's modulus of RDHNT (E = 165 MPa) and RNT (E = 144 MPa) was found to be almost the same and nearly half of that found for DNT (E = 325 MPa). Furthermore, the results showed that RNT was more resistant to bending, torsional, and volumetric deformations than DNT and RDHNT. We also used non-equilibrium SMD simulations to acquire comprehensive knowledge of the mechanical response of nanotubes to tensile stress.

Two Different Ways of Stress Release in Supercoiled DNA Minicircles under DNA Nick

It is generally believed that DNA nick is an effective way to release stress in supercoiled DNA, resulting from the twisting motion that individual strands rotate around the axis of the DNA helix. Here, we use MD simulations based on the oxDNA model to investigate the relaxation of 336 bp supercoiled minicircular DNA under DNA nick. Our simulations show that stress release, characterized by the abrupt decrease in linking number, may be induced by two types of DNA motion depending on the nick position. Except for the twisting motion, there is a writhing motion, that is, double strands collectively rotating with one plectoneme removal, which may occur in the process of DNA relaxation with the nick position in the loop region. Moreover, the writhing motion is more likely to occur in the DNA with relatively high hardness, such as C-G pairs. Our simulation results uncover the relationship between structural transformation, stress release, and DNA motion during the dynamic process under DNA nick, indicating the influence of nick position on the relaxation of the supercoiled DNA.

Tuning curved DNA origami structures through mechanical design and chemical adducts

The bending and twisting of DNA origami structures are important features for controlling the physical properties of DNA nanodevices. It has not been fully explored yet how to finely tune the bending and twisting of curved DNA structures. Traditional tuning of the curved DNA structures was limited to controlling the in-plane-bending angle through varying the numbers of base pairs of deletions and insertions. Here, we developed two tuning strategies of curved DNA origami structures fromin silicoandin vitroaspects.In silico, the out-of-plane bending and twisting angles of curved structures were introduced, and were tuned through varying the patterns of base pair deletions and insertions.In vitro, a chemical adduct (ethidium bromide) was applied to dynamically tune a curved spiral. The 3D structural conformations, like chirality, of the curved DNA structures were finely tuned through these two strategies. The simulation and TEM results demonstrated that the patterns of base pair insertions and deletions and chemical adducts could effectively tune the bending and twisting of curved DNA origami structures. These strategies expand the programmable accuracy of curved DNA origami structures and have potential in building efficient dynamic functional nanodevices.

The Free-Energy Landscape of a Mechanically Bistable DNA Origami

Molecular simulations using coarse-grained models allow the structure, dynamics and mechanics of DNA origamis to be comprehensively characterized. Here, we focus on the free-energy landscape of a jointed DNA origami that has been designed to exhibit two mechanically stable states and for which a bistable landscape has been inferred from ensembles of structures visualized by electron microscopy. Surprisingly, simulations using the oxDNA model predict that the defect-free origami has a single free-energy minimum. The expected second state is not stable because the hinge joints do not simply allow free angular motion but instead lead to increasing free-energetic penalties as the joint angles relevant to the second state are approached. This raises interesting questions about the cause of this difference between simulations and experiment, such as how assembly defects might affect the ensemble of structures observed experimentally.

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.

Characterizing the free-energy landscapes of DNA origamis

We show how coarse-grained modelling combined with umbrella sampling using distance-based order parameters can be applied to compute the free-energy landscapes associated with mechanical deformations of large DNA nanostructures. We illustrate this approach for the strong bending of DNA nanotubes and the potentially bistable landscape of twisted DNA origami sheets. The homogeneous bending of the DNA nanotubes is well described by the worm-like chain model; for more extreme bending the nanotubes reversibly buckle with the bending deformations localized at one or two "kinks". For a twisted one-layer DNA origami, the twist is coupled to the bending of the sheet giving rise to a free-energy landscape that has two nearly-degenerate minima that have opposite curvatures. By contrast, for a two-layer origami, the increased stiffness with respect to bending leads to a landscape with a single free-energy minimum that has a saddle-like geometry. The ability to compute such landscapes is likely to be particularly useful for DNA mechanotechnology and for understanding stress accumulation during the self-assembly of origamis into higher-order structures.

Design Approaches and Computational Tools for DNA Nanostructures

Designing a structure in nanoscale with desired shape and properties has been enabled by structural DNA nanotechnology. Design strategies in this research field have evolved to interpret various aspects of increasingly more complex nanoscale assembly and to realize molecular-level functionality by exploring static to dynamic characteristics of the target structure. Computational tools have naturally been of significant interest as they are essential to achieve a fine control over both shape and physicochemical properties of the structure. Here, we review the basic design principles of structural DNA nanotechnology together with its computational analysis and design tools.

A Primer on the oxDNA Model of DNA: When to Use it, How to Simulate it and How to Interpret the Results

The oxDNA model of Deoxyribonucleic acid has been applied widely to systems in biology, biophysics and nanotechnology. It is currently available via two independent open source packages. Here we present a set of clearly documented exemplar simulations that simultaneously provide both an introduction to simulating the model, and a review of the model's fundamental properties. We outline how simulation results can be interpreted in terms of-and feed into our understanding of-less detailed models that operate at larger length scales, and provide guidance on whether simulating a system with oxDNA is worthwhile.

Length-scale-dependent elasticity in DNA from coarse-grained and all-atom models

We investigate the influence of nonlocal couplings on the torsional and bending elasticities of DNA. Such couplings have been observed in the past by several simulation studies. Here, we use a description of DNA conformations based on the variables tilt, roll, and twist. Our analysis of both coarse-grained (oxDNA) and all-atom models indicates that these share strikingly similar features: there are strong off-site couplings for tilt-tilt and twist-twist, while they are much weaker in the roll-roll case. By developing an analytical framework to estimate bending and torsional persistence lengths in nonlocal DNA models, we show how off-site interactions generate a length-scale-dependent elasticity. Based on the simulation-generated elasticity data, the theory predicts a significant length-scale-dependent effect on torsional fluctuations but only a modest effect on bending fluctuations. These results are in agreement with experiments probing DNA mechanics from single base pair to kilobase pair scales.