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.

All-Atom Molecular Dynamics Simulations of Membrane-Spanning DNA Origami Nanopores

Building on the recent technological advances, all-atom molecular dynamics (MD) simulations have become an indispensable tool to study the molecular behavior at nanoscale. Molecular simulations have been used to characterize the structure, dynamics, and mechanical and electrical properties of DNA origami objects. In this chapter we describe a method to build all-atom model of lipid-spanning DNA origami nanopores and perform molecular dynamics simulations in explicit electrolyte solutions.

Non-linear Hamiltonian models for DNA

Nucleic acids' physical properties have been investigated by theoretical methods based both on fully atomistic representations and on coarse-grained models, e.g., the worm-like-chain, taken from polymer physics. In this review article, I discuss an intermediate (mesoscopic) approach and show how to build a three-dimensional Hamiltonian model which accounts for the main interactions responsible for the stability of the helical molecules. While the 3D mesoscopic model yields a sufficiently detailed description of the helix at the level of the base pair, it also allows one to predict the thermodynamical and structural properties of molecules in solution. Relying on the idea that the base pair fluctuations can be conceived as trajectories, I have built over the past years a computational method based on the time-dependent path integral formalism to derive the partition function. While the main features of the method are presented, I focus here in particular on a newly developed statistical method to set the maximum amplitude of the base pair fluctuations, a key parameter of the theory. Some applications to the calculation of DNA flexibility properties are discussed together with the available experimental data.

Mechanical properties of DNA and DNA nanostructures: comparison of atomistic, Martini and oxDNA models

The flexibility and stiffness of small DNA molecules play a fundamental role ranging from several biophysical processes to nano-technological applications. Here, we estimate the mechanical properties of short double-stranded DNA (dsDNA) with lengths ranging from 12 base-pairs (bp) to 56 bp, paranemic crossover (PX) DNA and hexagonal DNA nanotubes (DNTs) using two widely used coarse-grained models - Martini and oxDNA. To calculate the persistence length (Lp) and the stretch modulus (γ) of the dsDNA, we incorporate the worm-like chain and elastic rod model, while for the DNTs, we implement our previously developed theoretical framework. We compare and contrast all of the results with previously reported all-atom molecular dynamics (MD) simulations and experimental results. The mechanical properties of dsDNA (Lp ∼ 50 nm, γ ∼ 800-1500 pN), PX DNA (γ ∼ 1600-2000 pN) and DNTs (Lp ∼ 1-10 μm, γ ∼ 6000-8000 pN) estimated using the Martini soft elastic network and oxDNA are in very good agreement with the all-atom MD and experimental values, while the stiff elastic network Martini reproduces values of Lp and γ which are an order of magnitude higher. The high flexibility of small dsDNA is also depicted in our calculations. However, Martini models proved inadequate to capture the salt concentration effects on the mechanical properties with increasing salt molarity. oxDNA captures the salt concentration effect on the small dsDNA mechanics. But it is found to be ineffective for reproducing the salt-dependent mechanical properties of DNTs. Also, unlike Martini, the time evolved PX DNA and DNT structures from the oxDNA models are comparable to the all-atom MD simulated structures. Our findings provide a route to study the mechanical properties of DNA and DNA based nanostructures with increased time and length scales and has a remarkable implication in the context of DNA nanotechnology.

Chiral Systems Made from DNA

The very chemical structure of DNA that enables biological heredity and evolution has non-trivial implications for the self-organization of DNA molecules into larger assemblies and provides limitless opportunities for building functional nanostructures. This progress report discusses the natural organization of DNA into chiral structures and recent advances in creating synthetic chiral systems using DNA as a building material. How nucleic acid chirality naturally comes into play in a diverse array of situations is considered first, at length scales ranging from an individual nucleotide to entire chromosomes. Thereafter, chiral liquid crystal phases formed by dense DNA mixtures are discussed, including the ongoing efforts to understand their origins. The report then summarizes recent efforts directed toward building chiral structures, and other structures of complex topology, using the principle of DNA self-assembly. Discussed last are existing and proposed functional man-made nanostructures designed to either probe or harness DNA's chirality, from plasmonics and spintronics to biosensing.

Computing the Elastic Mechanical Properties of Rodlike DNA Nanostructures

To study the elastic properties of rodlike DNA nanostructures, we perform long simulations of these structures using the oxDNA coarse-grained model. By analyzing the fluctuations in these trajectories, we obtain estimates of the bend and twist persistence lengths and the underlying bend and twist elastic moduli and couplings between them. Only on length scales beyond those associated with the spacings between the interhelix crossovers do the bending fluctuations behave like those of a wormlike chain. The obtained bending persistence lengths are much larger than that for double-stranded DNA and increase nonlinearly with the number of helices, whereas the twist moduli increase approximately linearly. To within the numerical error in our data, the twist-bend coupling constants are of order zero. That the bending persistence lengths that we obtain are generally somewhat higher than in experiment probably reflects both that the simulated origamis have no assembly defects and that the oxDNA extensional modulus for double-stranded DNA is too large.

First-passage probability: a test for DNA Hamiltonian parameters

A method is developed to chose the set of input parameters for DNA mesoscopic Hamiltonian models.

Atomic structures of RNA nanotubes and their comparison with DNA nanotubes

We present a computational framework to model RNA based nanostructures and study their microscopic structures. We model hexagonal nanotubes made of 6 dsRNA (RNTs) connected by double crossover (DX) at different positions. Using several hundred nano-second (ns) long all-atom molecular dynamics simulations, we study the atomic structure, conformational change and elastic properties of RNTs in the presence of explicit water and ions. Based on several structural quantities such as root mean square deviation (RMSD) and root mean square fluctuation (RMSF), we find that the RNTs are almost as stable as DNA nanotubes (DNTs). Although the central portion of the RNTs maintain its cylindrical shape, both the terminal regions open up to give rise to a gating like behavior which can play a crucial role in drug delivery. From the bending angle distribution, we observe that the RNTs are more flexible than DNTs. The calculated persistence length of the RNTs is in the micron range which is an order of magnitude higher than that of a single dsRNA. The stretch modulus of the RNTs from the contour length distribution is in the range of 4-7 nN depending on the sequence. The calculated persistence length and stretch modulus are in the same range of values as in the case of DNTs. To understand the structural properties of RNTs at the individual base-pair level we have also calculated all the helicoidal parameters and analyzed the relative flexibility and rigidity of RNTs having a different sequence. These findings emphasized the fascinating properties of RNTs which will expedite further theoretical and experimental studies in this field.

Biomimetic DNA Nanotubes: Nanoscale Channel Design and Applications

Biomacromolecular nanotubes play important physiological roles in transmembrane ion/molecule channeling, intracellular transport, and inter-cellular communications. While genetically encoded protein nanotubes are prevalent in vivo, the in vitro construction of biomimetic DNA nanotubes has attracted intense interest with the rise of structural DNA nanotechnology. The abiotic use of DNA assembly provides a powerful bottom-up approach for the rational construction of complex materials with arbitrary size and shape at the nanoscale. More specifically, a typical DNA nanotube can be assembled either with parallel-aligned DNA duplexes or by closing DNA tile lattices. These artificial DNA nanotubes can be tailored and site-specifically modified to realize biomimetic functions including ionic or molecular channeling, bioreactors, drug delivery, and biomolecular sensing. In this Minireview, we aim to summarize recent advances in design strategies, including the characterization and applications of biomimetic DNA nanotubes.

Ionic liquids make DNA rigid

Persistence length of double-stranded DNA (dsDNA) is known to decrease with an increase in ionic concentration of the solution. In contrast to this, here we show that the persistence length of dsDNA increases dramatically as a function of ionic liquid (IL) concentration. Using all atom explicit solvent molecular dynamics simulations and theoretical models, we present, for the first time, a systematic study to determine the mechanical properties of dsDNA in various hydrated ILs at different concentrations. We find that dsDNA in 50 wt % ILs have lower persistence length and stretch modulus in comparison to 80 wt % ILs. We further observe that both the persistence length and stretch modulus of dsDNA increase as we increase the concentration of ILs. The present trend of the stretch modulus and persistence length of dsDNA with IL concentration supports the predictions of the macroscopic elastic theory, in contrast to the behavior exhibited by dsDNA in monovalent salt. Our study further suggests the preferable ILs that can be used for maintaining DNA stability during long-term storage.

Charge Neutralization Drives the Shape Reconfiguration of DNA Nanotubes

Reconfiguration of membrane protein channels for gated transport is highly regulated under physiological conditions. However, a mechanistic understanding of such channels remains challenging owing to the difficulty in probing subtle gating-associated structural changes. Herein, we show that charge neutralization can drive the shape reconfiguration of a biomimetic 6-helix bundle DNA nanotube (6HB). Specifically, 6HB adopts a compact state when its charge is neutralized by Mg2+ ; whereas Na+ switches it to the expanded state, as revealed by MD simulations, small-angle X-ray scattering (SAXS), and FRET characterization. Furthermore, partial neutralization of the DNA backbone charges by chemical modification renders 6HB compact and insensitive to ions, suggesting an interplay between electrostatic and hydrophobic forces in the channels. This system provides a platform for understanding the structure-function relationship of biological channels and designing rules for the shape control of DNA nanostructures in biomedical applications.

Structure and electrical properties of DNA nanotubes embedded in lipid bilayer membranes

Engineering the synthetic nanopores through lipid bilayer membrane to access the interior of a cell is a long persisting challenge in biotechnology. Here, we demonstrate the stability and dynamics of a tile-based 6-helix DNA nanotube (DNT) embedded in POPC lipid bilayer using the analysis of 0.2 μs long equilibrium MD simulation trajectories. We observe that the head groups of the lipid molecules close to the lumen cooperatively tilt towards the hydrophilic sugar-phosphate backbone of DNA and form a toroidal structure around the patch of DNT protruding in the membrane. Further, we explore the effect of ionic concentrations to the in-solution structure and stability of the lipid-DNT complex. Transmembrane ionic current measurements for the constant electric field MD simulation provide the I-V characteristics of the water filled DNT lumen in lipid membrane. With increasing salt concentrations, the measured values of transmembrane ionic conductance of the porous DNT lumen vary from 4.3 to 20.6 nS. Simulations of the DNTs with ssDNA and dsDNA overhangs at the mouth of the pore show gating effect with remarkable difference in the transmembrane ionic conductivities for open and close state nanopores.

Conformational Dynamics of Mechanically Compliant DNA Nanostructures from Coarse-Grained Molecular Dynamics Simulations

Structural DNA nanotechnology, the assembly of rigid 3D structures of complex yet precise geometries, has recently been used to design dynamic, mechanically compliant nanostructures with tunable equilibrium conformations and conformational distributions. Here we use coarse-grained molecular dynamics simulations to provide insights into the conformational dynamics of a set of mechanically compliant DNA nanostructures-DNA hinges that use single-stranded DNA "springs" to tune the equilibrium conformation of a layered double-stranded DNA "joint" connecting two stiff "arms" constructed from DNA helix bundles. The simulations reproduce the experimentally measured equilibrium angles between hinge arms for a range of hinge designs. The hinges are found to be structurally stable, except for some fraying of the open ends of the DNA helices comprising the hinge arms and some loss of base-pairing interactions in the joint regions coinciding with the crossover junctions, especially in hinges designed to exhibit a small bending angle that exhibit large local stresses resulting in strong kinks in their joints. Principal component analysis reveals that while the hinge dynamics are dominated by bending motion, some twisting and sliding of hinge arms relative to each other also exists. Forced deformation of the hinges reveals distinct bending mechanisms for hinges with short, inextensible springs versus those with longer, more extensible springs. Lastly, we introduce an approach for rapidly predicting equilibrium hinge angles from individual force-deformation behaviors of its single- and double-stranded DNA components. Taken together, these results demonstrate that coarse-grained modeling is a promising approach for designing, predicting, and studying the dynamics of compliant DNA nanostructures, where conformational fluctuations become important, multiple deformation mechanisms exist, and continuum approaches may not yield accurate properties.

Probing the structure and in silico stability of cargo loaded DNA icosahedra using MD simulations

Platonic solids such as polyhedra based on DNA have been deployed for multifarious applications such as RNAi delivery, biological targeting and bioimaging. All of these applications hinge on the capability of DNA polyhedra for molecular display with high spatial precision. Therefore high resolution structural models of such polyhedra are critical to widen their applications in both materials and biology. Here, we present an atomistic model of a well-characterized DNA icosahedron, with demonstrated versatile functionalities in biological systems. We study the structure and dynamics of this DNA icosahedron using fully atomistic molecular dynamics (MD) simulation in explicit water and ions. The major modes of internal motion have been identified using principal component analysis. We provide a quantitative estimate of the radius of gyration (R_g), solvent accessible surface area (SASA) and volume of the icosahedron which is essential to estimate its maximal cargo carrying capacity. Importantly, our simulation of gold nanoparticles (AuNPs) encapsulated within DNA icosahedra revealed enhanced stability of the AuNP loaded DNA icosahedra compared to empty icosahedra. This is consistent with the experimental results that show high yields of cargo-encapsulated DNA icosahedra that have led to its diverse applications for precision targeting. These studies reveal that the stabilizing interactions between the cargo and the DNA scaffold powerfully position DNA polyhedra as targetable nanocapsules for payload delivery. These insights can be exploited for precise molecular display for diverse biological applications.

Studying the adsorption of DNA nanostructures on graphene in the aqueous phase using molecular dynamic simulations

DNA nanostructures can undergo large structural fluctuations and deviate from their intended configurations. In this work, two model DNA nanostructures (i.e., Nan and Kai) were designed based on the shape of the two Chinese characters of the name of Nankai University, and additional single-stranded DNA fragments were added to interact with graphene. During four 50-ns molecular dynamic simulations in aqueous solution, the DNA nanostructures adsorbed onto graphene demonstrated more stable conformations with lower root mean square deviations and smaller coordinate changes in the z-axis direction than the DNA nanostructures that were not adsorbed onto graphene. The interaction analyses and energetic calculations show that π-π interactions between single-stranded DNA and graphene are necessary for adsorption of the DNA nanostructures. Overall, this work examined the interactions between DNA and graphene at a large spatial scale with the hope that it provides a new strategy to stabilize DNA nanostructures.

Quantum dot-loaded monofunctionalized DNA icosahedra for single-particle tracking of endocytic pathways

Functionalization of quantum dots (QDs) with a single biomolecular tag using traditional approaches in bulk solution has met with limited success. DNA polyhedra consist of an internal void bounded by a well-defined three-dimensional structured surface. The void can house cargo and the surface can be functionalized with stoichiometric and spatial precision. Here, we show that monofunctionalized QDs can be realized by encapsulating QDs inside DNA icosahedra and functionalizing the DNA shell with an endocytic ligand. We deployed the DNA-encapsulated QDs for real-time imaging of three different endocytic ligands-folic acid, galectin-3 (Gal3) and the Shiga toxin B-subunit (STxB). Single-particle tracking of Gal3- or STxB-functionalized QD-loaded DNA icosahedra allows us to monitor compartmental dynamics along endocytic pathways. These DNA-encapsulated QDs, which bear a unique stoichiometry of endocytic ligands, represent a new class of molecular probes for quantitative imaging of endocytic receptor dynamics.

Nanoscale Structure and Elasticity of Pillared DNA Nanotubes

We present an atomistic model of pillared DNA nanotubes (DNTs) and their elastic properties which will facilitate further studies of these nanotubes in several important nanotechnological and biological applications. In particular, we introduce a computational design to create an atomistic model of a 6-helix DNT (6HB) along with its two variants, 6HB flanked symmetrically with two double helical DNA pillars (6HB+2) and 6HB flanked symmetrically by three double helical DNA pillars (6HB+3). Analysis of 200 ns all-atom simulation trajectories in the presence of explicit water and ions shows that these structures are stable and well behaved in all three geometries. Hydrogen bonding is well maintained for all variants of 6HB DNTs. From the equilibrium bending angle distribution, we calculate the persistence lengths of these tubes. The measured persistence lengths of these nanotubes are ∼10 μm, which is 2 orders of magnitude larger than that of dsDNA. We also find a gradual increase of persistence length with an increasing number of pillars, in quantitative agreement with previous experimental findings. To have a quantitative understanding of the stretch modulus of these tubes, we carried out nonequilibrium steered molecular dynamics (SMD). The linear part of the force-extension plot gives a stretch modulus in the range 6500 pN for 6HB without pillars, which increases to 11 000 pN for tubes with three pillars. The values of the stretch modulus calculated using contour length distribution obtained from equilibrium MD simulations are similar to those obtained from nonequilibrium SMD simulations. The addition of pillars makes these DNTs very rigid.

Comparative simulative analysis of single and double stranded truncated octahedral DNA nanocages

Spacefill view of double (DSL) and single (SSL) stranded linkers DNA cages. The blue atoms represent the shared cages scaffold, while the yellow atoms show the single stranded DNA oligonucleotides shaping the double stranded linkers of the DSL cage.

Gating-like Motions and Wall Porosity in a DNA Nanopore Scaffold Revealed by Molecular Simulations

Recently developed synthetic membrane pores composed of folded DNA enrich the current range of natural and engineered protein pores and of nonbiogenic channels. Here we report all-atom molecular dynamics simulations of a DNA nanotube (DNT) pore scaffold to gain fundamental insight into its atomic structure, dynamics, and interactions with ions and water. Our multiple simulations of models of DNTs that are composed of a six-duplex bundle lead to a coherent description. The central tube lumen adopts a cylindrical shape while the mouth regions at the two DNT openings undergo gating-like motions which provide a possible molecular explanation of a lower conductance state observed in our previous experimental study on a membrane-spanning version of the DNT (ACS Nano 2015, 9, 1117-26). Similarly, the central nanotube lumen is filled with water and ions characterized by bulk diffusion coefficients while the gating regions exhibit temporal fluctuations in their aqueous volume. We furthermore observe that the porous nature of the walls allows lateral leakage of ions and water. This study will benefit rational design of DNA nanopores of enhanced stability of relevance for sensing applications, of nanodevices with tunable gating properties that mimic gated ion channels, or of nanopores featuring defined permeation behavior.