Introducing improved structural properties and salt dependence into a coarse-grained model of DNA

Programming rigidity into size-defined wireframe DNA nanotubes

Nanotubes built from DNA hold promise for several biological and materials applications, due to their high aspect ratio and encapsulation potential. A particularly appealing goal is to control the size, shape, and dynamic behaviour of DNA nanotubes with minimal design alteration, as nanostructures of varying morphologies and lengths have been shown to exhibit distinct cellular uptake, encapsulation behaviour, and in vivo biodistribution. Herein, we report a systematic investigation, combining experimental and computational design, to modulate the length, flexibility, and longitudinal patterns of wireframe DNA nanotubes. Subtle design changes govern the structure and properties of our nanotubes, which are built from a custom-made, long, and size-defined template strand to which DNA rungs and linkers are attached. Unlike DNA origami, these custom-made strands possess regions with repeating sequences at strategic locations, thereby reducing the number of strands necessary for assembly. Through strand displacement, the nanotubes can be reversibly altered between extended and collapsed morphologies. These design concepts enable fine-tuning of the nanotube stiffness and may pave the way for the development of designer nanotubes for a variety of applications, including the study of cellular internalization, biodistribution, and uptake mechanisms for structures of varied shapes and sizes.

Design, molecular dynamics simulation, and investigation of the mechanical behavior of DNA origami nanotubes with auxetic and honeycomb structures

Numerous applications of DNA origami nanotubes for load-bearing purposes necessitate the improvement of properties and mechanical behavior of these types of structures, as well as the use of innovative structures such as metamaterials. To this end, the present study aims to investigate the design, molecular dynamics (MD) simulation, and mechanical behavior of DNA origami nanotube structures consisting of honeycomb and re-entrant auxetic cross-sections. The results revealed both structures kept their structural stability. In addition, DNA origami based-nanotube with auxetic cross-section exhibits negative Poisson's ratio (NPR) under tensile loading. Furthermore, MD simulation results demonstrated that the values of stiffness, specific stiffness, energy absorption, and specific energy absorption in the structure with an auxetic cross-section are higher than that of a honeycomb cross-section, similar to their behavior in macro-scale structures. The finding of this study is to propose re-entrant auxetic structure as the next generation of DNA origami nanotubes. In addition, it can be utilized to aid scientists with the design and fabrication of novel auxetic DNA origami structures.Communicated by Ramaswamy H. Sarma.

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.

DNA Droplets: Intelligent, Dynamic Fluid

Breathtaking advances in DNA nanotechnology have established DNA as a promising biomaterial for the fabrication of programmable higher-order nano/microstructures. In the context of developing artificial cells and tissues, DNA droplets have emerged as a powerful platform for creating intelligent, dynamic cell-like machinery. DNA droplets are a microscale membrane-free coacervate of DNA formed through phase separation. This new type of DNA system couples dynamic fluid-like property with long-established DNA programmability. This hybrid nature offers an advantageous route to facile and robust control over the structures, functions, and behaviors of DNA droplets. This review begins by describing programmable DNA condensation, commenting on the physical properties and fabrication strategies of DNA hydrogels and droplets. By presenting an overview of the development pathways leading to DNA droplets, it is shown that DNA technology has evolved from static, rigid systems to soft, dynamic systems. Next, the basic characteristics of DNA droplets are described as intelligent, dynamic fluid by showcasing the latest examples highlighting their distinctive features related to sequence-specific interactions and programmable mechanical properties. Finally, this review discusses the potential and challenges of numerical modeling able to connect a robust link between individual sequences and macroscopic mechanical properties of DNA droplets.

Interplay of the mechanical and structural properties of DNA nanostructures determines their electrostatic interactions with lipid membranes

Nucleic acids and lipids function in close proximity in biological processes, as well as in nanoengineered constructs for therapeutic applications. As both molecules carry a rich charge profile, and frequently coexist in complex ionic solutions, the electrostatics surely play a pivotal role in interactions between them. Here we discuss how each component of a DNA/ion/lipid system determines its electrostatic attachment. We examine membrane binding of a library of DNA molecules varying from nanoengineered DNA origami through plasmids to short DNA domains, demonstrating the interplay between the molecular structure of the nucleic acid and the phase of lipid bilayers. Furthermore, the magnitude of DNA/lipid interactions is tuned by varying the concentration of magnesium ions in the physiologically relevant range. Notably, we observe that the structural and mechanical properties of DNA are critical in determining its attachment to lipid bilayers and demonstrate that binding is correlated positively with the size, and negatively with the flexibility of the nucleic acid. The findings are utilized in a proof-of-concept comparison of membrane interactions of two DNA origami designs - potential nanotherapeutic platforms - showing how the results can have a direct impact on the choice of DNA geometry for biotechnological applications.

Nonequilibrium Thermodynamics of DNA Nanopore Unzipping

Using theory and simulations, we carried out a first systematic characterization of DNA unzipping via nanopore translocation. Starting from partially unzipped states, we found three dynamical regimes depending on the applied force f: (i) heterogeneous DNA retraction and rezipping (f<17  pN), (ii) normal (17  pN<f<60  pN), and (iii) anomalous (f>60  pN) drift-diffusive behavior. We show that the normal drift-diffusion regime can be effectively modeled as a one-dimensional stochastic process in a tilted periodic potential. We use the theory of stochastic processes to recover the potential from nonequilibrium unzipping trajectories and show that it corresponds to the free-energy landscape for single-base-pair unzipping. Applying this general approach to other single-molecule systems with periodic potentials ought to yield detailed free-energy landscapes from out-of-equilibrium trajectories.

Non-Invasive, Reliable, and Fast Quantification of DNA Loading on Gold Nanoparticles by a One-Step Optical Measurement

An exquisite, versatile, and reproducible quantification of DNA loading on gold nanoparticles (Au NPs) has long been pursued because this loading influences the analytical, therapeutic, and self-assembly behaviors of DNA-Au NPs. Nevertheless, the existing methods used thus far rely solely on the invasive detachment and subsequent spectroscopic quantification of DNA, which are error-prone and highly dependent on trained personnel. Here, we present a non-invasive optical framework that can determine the number of DNA strands on Au NPs by versatile one-step measurement of the visible absorption spectra of DNA-Au NP solutions without any invasive modifications or downstream processes. Using effective medium theory in conjunction with electromagnetic numerical calculation, the change in DNA loading density, resulting from varying the ion concentration, Au NP size, DNA strand length, and surrounding temperature, can be tracked in situ merely by the one-step measurement of visible absorption spectra, which is otherwise impossible to achieve. Moreover, the simplicity and robustness of this method promote reproducible DNA loading quantification regardless of experimental adeptness, which is in stark contrast with existing invasive and multistep methods. Overall, the optical framework outlined in this work can contribute to democratizing research on DNA-Au NPs and facilitating their rapid adoption in transformative applications.

DNA origami tubes with reconfigurable cross-sections

Structural DNA nanotechnology has enabled the design and construction of complex nanoscale structures with precise geometry and programmable dynamic and mechanical properties. Recent efforts have led to major advances in the capacity to actuate shape changes of DNA origami devices and incorporate DNA origami into larger assemblies, which open the prospect of using DNA to design shape-morphing assemblies as components of micro-scale reconfigurable or sensing materials. Indeed, a few studies have constructed higher order assemblies with reconfigurable devices; however, these demonstrations have utilized structures with relatively simple motion, primarily hinges that open and close. To advance the shape changing capabilities of DNA origami assemblies, we developed a multi-component DNA origami 6-bar mechanism that can be reconfigured into various shapes and can be incorporated into larger assemblies while maintaining capabilities for a variety of shape transformations. We demonstrate the folding of the 6-bar mechanism into four different shapes and demonstrate multiple transitions between these shapes. We also studied the shape preferences of the 6-bar mechanism in competitive folding reactions to gain insight into the relative free energies of the shapes. Furthermore, we polymerized the 6-bar mechanism into tubes with various cross-sections, defined by the shape of the individual mechanism, and we demonstrate the ability to change the shape of the tube cross-section. This expansion of current single-device reconfiguration to higher order scales provides a foundation for nano to micron scale DNA nanotechnology applications such as biosensing or materials with tunable properties.

The oxDNA Coarse-Grained Model as a Tool to Simulate DNA Origami

This chapter introduces how to run molecular dynamics simulations for DNA origami using the oxDNA coarse-grained model.

Electrochemical response of surface-attached redox DNA governed by low activation energy electron transfer kinetics

We demonstrate, using high scan rate cyclic voltammetry and molecular dynamics simulations, that the electrochemical response of electrode-attached redox DNA is governed by low reorganization energy electron transfer kinetics.

Automated design of 3D DNA origami with non-rasterized 2D curvature

Improving the precision and function of encapsulating three-dimensional (3D) DNA nanostructures via curved geometries could have transformative impacts on areas such as molecular transport, drug delivery, and nanofabrication. However, the addition of non-rasterized curvature escalates design complexity without algorithmic regularity, and these challenges have limited the ad hoc development and usage of previously unknown shapes. In this work, we develop and automate the application of a set of previously unknown design principles that now includes a multilayer design for closed and curved DNA nanostructures to resolve past obstacles in shape selection, yield, mechanical rigidity, and accessibility. We design, analyze, and experimentally demonstrate a set of diverse 3D curved nanoarchitectures, showing planar asymmetry and examining partial multilayer designs. Our automated design tool implements a combined algorithmic and numerical approximation strategy for scaffold routing and crossover placement, which may enable wider applications of general DNA nanostructure design for nonregular or oblique shapes.

The wending rhombus: Self-assembling 3D DNA crystals

In this perspective, we provide a summary of recent developments in self-assembling three-dimensional (3D) DNA crystals. Starting from the inception of this subfield, we describe the various advancements in structure that have led to an increase in the diversity of macromolecular crystal motifs formed through self-assembly, and we further comment on the future directions of the field, which exploit noncanonical base pairing interactions beyond Watson-Crick. We then survey the current applications of self-assembling 3D DNA crystals in reversibly active nanodevices and materials engineering and provide an outlook on the direction researchers are taking these structures. Finally, we compare 3D DNA crystals with DNA origami and suggest how these distinct subfields might work together to enhance biomolecule structure solution, nanotechnological motifs, and their applications.

Energy landscapes of rotary DNA origami devices determined by fluorescence particle tracking

Biomolecular nanomechanical devices are of great interest as tools for the processing and manipulation of molecules, thereby mimicking the function of nature's enzymes. DNA nanotechnology provides the capability to build molecular analogs of mechanical machine elements such as joints and hinges via sequence-programmable self-assembly, which are otherwise known from traditional mechanical engineering. Relative to their size, these molecular machine elements typically do not reach the same relative precision and reproducibility that we know from their macroscopic counterparts; however, as they are scaled down to molecular sizes, physical effects typically not considered by mechanical engineers such as Brownian motion, intramolecular forces, and the molecular roughness of the devices begin to dominate their behavior. In order to investigate the effect of different design choices on the roughness of the mechanical energy landscapes of DNA nanodevices in greater detail, we here study an exemplary DNA origami-based structure, a modularly designed rotor-stator arrangement, which resembles a rotatable nanorobotic arm. Using fluorescence tracking microscopy, we follow the motion of individual rotors and record their corresponding energy landscapes. We then utilize the modular construction of the device to exchange its constituent parts individually and systematically test the effect of different design variants on the movement patterns. This allows us to identify the design parameters that most strongly affect the shape of the energy landscapes of the systems. Taking into account these insights, we are able to create devices with significantly flatter energy landscapes, which translates to mechanical nanodevices with improved performance and behaviors more closely resembling those of their macroscopic counterparts.

Computer Aided Development of Nucleic Acid Applications in Nanotechnologies

Utilization of nucleic acids (NAs) in nanotechnologies and nanotechnology-related applications is a growing field with broad application potential, ranging from biosensing up to targeted cell delivery. Computer simulations are useful techniques that can aid design and speed up development in this field. This review focuses on computer simulations of hybrid nanomaterials composed of NAs and other components. Current state-of-the-art molecular dynamics simulations, empirical force fields (FFs), and coarse-grained approaches for the description of deoxyribonucleic acid and ribonucleic acid are critically discussed. Challenges in combining biomacromolecular and nanomaterial FFs are emphasized. Recent applications of simulations for modeling NAs and their interactions with nano- and biomaterials are overviewed in the fields of sensing applications, targeted delivery, and NA templated materials. Future perspectives of development are also highlighted.

A bistable and reconfigurable molecular system with encodable bonds

Molecular systems with ability to controllably transform between different conformations play pivotal roles in regulating biochemical functions. Here, we report the design of a bistable DNA origami four-way junction (DOJ) molecular system that adopts two distinct stable conformations with controllable reconfigurability by using conformation-controlled base stacking. Exquisite control over DOJ's conformation and transformation is realized by programming the stacking bonds (quasi-blunt-ends) within the junction to induce prescribed coaxial stacking of neighboring junction arms. A specific DOJ conformation may be achieved by encoding the stacking bonds with binary stacking sequences based on thermodynamic calculations. Dynamic transformations of DOJ between various conformations are achieved by using specific environmental and molecular stimulations to reprogram the stacking codes. This work provides a useful platform for constructing self-assembled DNA nanostructures and nanomachines and insights for future design of artificial molecular systems with increasing complexity and reconfigurability.

Molecular states and spin crossover of hemin studied by DNA origami enabled single-molecule surface-enhanced Raman scattering

The study of biologically relevant molecules and their interaction with external stimuli on a single molecular scale is of high importance due to the availability of distributed rather than averaged information. Surface enhanced Raman scattering (SERS) provides direct chemical information, but is rather challenging on the single molecule (SM) level, where it is often assumed to require a direct contact of analyte molecules with the metal surface. Here, we detect and investigate the molecular states of single hemin by SM-SERS. A DNA aptamer based G-quadruplex mediated recognition of hemin directs its placement in the SERS hot-spot of a DNA Origami Nanofork Antenna (DONA). The configuration of the DONA structure allows the molecule to be trapped at the plasmonic hot-spot preferentially in no-contact configuration with the metal surface. Owing to high field enhancement at the plasmonic hot spot, the detection of a single folded G-quadruplex becomes possible. For the first time, we present a systematic study by SM-SERS where most hemin molecule adopt a high spin and oxidation state (III) that showed state crossover to low spin upon strong-field-ligand binding. The present study therefore, provides a platform for studying biologically relevant molecules and their properties at SM sensitivity along with demonstrating a conceptual advancement towards successful monitoring of single molecular chemical interaction using DNA aptamers.

Mechanical deformation behaviors and structural properties of ligated DNA crystals

DNA self-assembly has emerged as a powerful strategy for constructing complex nanostructures. While the mechanics of individual DNA strands have been studied extensively, the deformation behaviors and structural properties of self-assembled architectures are not well understood. This is partly due to the small dimensions and limited experimental methods available. DNA crystals are macroscopic crystalline structures assembled from nanoscale motifs via sticky-end association. The large DNA constructs may thus be an ideal platform to study structural mechanics. Here, we investigate the fundamental mechanical properties and behaviors of ligated DNA crystals made of tensegrity triangular motifs. We perform coarse-grained molecular dynamics simulations and confirm the results with nanoindentation experiments using atomic force microscopy. We observe various deformation modes, including untension, linear elasticity, duplex dissociation, and single-stranded component stretch. We find that the mechanical properties of a DNA architecture are correlated with those of its components. However, the structure shows complex behaviors which may not be predicted by components alone and the architectural design must be considered.

Structure and Dynamics of dsDNA in Cell-like Environments

Deoxyribonucleic acid (DNA) is a fundamental biomolecule for correct cellular functioning and regulation of biological processes. DNA's structure is dynamic and has the ability to adopt a variety of structural conformations in addition to its most widely known double-stranded DNA (dsDNA) helix structure. Stability and structural dynamics of dsDNA play an important role in molecular biology. In vivo, DNA molecules are folded in a tightly confined space, such as a cell chamber or a channel, and are highly dense in solution; their conformational properties are restricted, which affects their thermodynamics and mechanical properties. There are also many technical medical purposes for which DNA is placed in a confined space, such as gene therapy, DNA encapsulation, DNA mapping, etc. Physiological conditions and the nature of confined spaces have a significant influence on the opening or denaturation of DNA base pairs. In this review, we summarize the progress of research on the stability and dynamics of dsDNA in cell-like environments and discuss current challenges and future directions. We include studies on various thermal and mechanical properties of dsDNA in ionic solutions, molecular crowded environments, and confined spaces. By providing a better understanding of melting and unzipping of dsDNA in different environments, this review provides valuable guidelines for predicting DNA thermodynamic quantities and for designing DNA/RNA nanostructures.

Predicting accurate ab initio DNA electron densities with equivariant neural networks

One of the fundamental limitations of accurately modeling biomolecules like DNA is the inability to perform quantum chemistry calculations on large molecular structures. We present a machine learning model based on an equivariant Euclidean neural network framework to obtain accurate ab initio electron densities for arbitrary DNA structures that are much too large for conventional quantum methods. The model is trained on representative B-DNA basepair steps that capture both base pairing and base stacking interactions. The model produces accurate electron densities for arbitrary B-DNA structures with typical errors of less than 1%. Crucially, the error does not increase with system size, which suggests that the model can extrapolate to large DNA structures with negligible loss of accuracy. The model also generalizes reasonably to other DNA structural motifs such as the A- and Z-DNA forms, despite being trained on only B-DNA configurations. The model is used to calculate electron densities of several large-scale DNA structures, and we show that the computational scaling for this model is essentially linear. We also show that this machine learning electron density model can be used to calculate accurate electrostatic potentials for DNA. These electrostatic potentials produce more accurate results compared with classical force fields and do not show the usual deficiencies at short range.

Coarse-Grained Simulations for the Characterization and Optimization of Hybrid Protein-DNA Nanostructures

We present here the combination of experimental and computational modeling tools for the design and characterization of protein-DNA hybrid nanostructures. Our work incorporates several features in the design of these nanostructures: (1) modeling of the protein-DNA linker identity and length; (2) optimizing the design of protein-DNA cages to account for mechanical stresses; (3) probing the incorporation efficiency of protein-DNA conjugates into DNA nanostructures. The modeling tools were experimentally validated using structural characterization methods like cryo-TEM and AFM. Our method can be used for fitting low-resolution electron density maps when structural insights cannot be deciphered from experiments, as well as enable in-silico validation of nanostructured systems before their experimental realization. These tools will facilitate the design of complex hybrid protein-DNA nanostructures that seamlessly integrate the two different biomolecules.

Heterogeneous migration routes of DNA triplet repeat slip-outs

It is unclear how the length of a repetitive DNA tract determines the onset and progression of repeat expansion diseases, but the dynamics of secondary DNA structures formed by repeat sequences are believed to play an important role. It was recently shown that three-way DNA junctions containing slip-out hairpins of CAG or CTG repeats and contiguous triplet repeats in the adjacent duplex displayed single-molecule FRET (smFRET) dynamics that were ascribed to both local conformational motions and longer-range branch migration. Here we explore these so-called "mobile" slip-out structures through a detailed kinetic analysis of smFRET trajectories and coarse-grained modeling. Despite the apparent structural simplicity, with six FRET states resolvable, most smFRET states displayed biexponential dwell-time distributions, attributed to structural heterogeneity and overlapping FRET states. Coarse-grained modeling for a (GAC)₁₀ repeat slip-out included trajectories that corresponded to a complete round of branch migration; the structured free energy landscape between slippage events supports the dynamical complexity observed by smFRET. A hairpin slip-out with 40 CAG repeats, which is above the repeat length required for disease in several triplet repeat disorders, displayed smFRET dwell times that were on average double those of 3WJs with 10 repeats. The rate of secondary-structure rearrangement via branch migration, relative to particular DNA processing pathways, may be an important factor in the expansion of triplet repeat expansion diseases.

Mutual effects between single-stranded DNA conformation and Na+-Mg2+ ion competition in mixed salt solutions

The ion-dependence of single-stranded DNA (ssDNA) conformational changes has attracted growing attention because of its biological and technological importance. Although single-species ion effects have been extensively explored, it is challenging to study the ssDNA conformational properties under mixed monovalent/divalent ion conditions due to the complications of ssDNA flexibility and ion-ion competition. In this study, we apply Langevin dynamics simulations to investigate mixed Na⁺/Mg²⁺ ion-dependent ssDNA conformations. The ssDNA structure is described using a coarse-grained model, in which the phosphate, base, and sugar of each nucleotide are represented by three different beads. A novel improvement in our simulation model is that mixed-salt-related electrostatic interactions are computed via combining Manning counterion condensation (MCC) theory with the Monte Carlo tightly bound ion (MCTBI) model. Based on this MCC-MCTBI combination, we report new empirical functions to describe the ion-concentration-dependent and ssDNA conformation/structure-dependent electrostatic effects. The calculation results relating to the ion binding properties and the simulation results relating to the ssDNA conformational properties are validated against experimental results. In addition, our simulation results suggest a quantitative relationship between the ssDNA conformation and Na⁺-Mg²⁺ competition; this in turn reveals their mutual impact in the ion atmosphere.

Redox-labelled electrochemical aptasensors with nanosupported cancer cells

The transfer of redox-labelled bioelectrochemical sensors from proteins to cells is not straightforward because of the cell downward force issue on the surface of the sensors. In this paper, 20-nm-thick nanopillars are introduced to overcome this issue, in a well-controlled manner. We show on both molecular dynamics simulations and experiments that suspending cells a few nanometers above an electrode surface enables redox-labelled tethered DNA aptamer probes to move freely, while remaining at an interaction distance from a target membrane protein, i. e. epithelial cell adhesion molecule (EpCAM), which is typically overexpressed in cancer cells. By this nanopillar configuration, the interaction of aptamer with cancer cells is clearly observable, with 13 cells as the lower limit of detection. Nanoconfinement induced by the gap between the electrode surface and the cell membrane appears to improve the limit of detection and to lower the melting temperature of DNA aptamer hairpins, offering an additional degree of freedom to optimize molecular recognition mechanisms. This novel nanosupported electrochemical DNA cell sensor scheme including Brownian-fluctuating redox species opens new opportunities for the design of all-electrical sensors using redox-labelled probes.

Selective Integrin α5β1 Targeting through Spatially Constrained Multivalent DNA-Based Nanoparticles

Targeting cells specifically based on receptor expression levels remains an area of active research to date. Selective binding of receptors cannot be achieved by increasing the individual binding strength, as this does not account for differing distributions of receptor density across healthy and diseased cells. Engaging receptors above a threshold concentration would be desirable in devising selective diagnostics. Integrins are prime target candidates as they are readily available on the cell surface and have been reported to be overexpressed in diseases. Insights into their spatial organization would therefore be advantageous to design selective targeting agents. Here, we investigated the effect of activation method on integrin α₅β₁ clustering by immunofluorescence and modeled the global neighbor distances with input from an immuno-staining assay and image processing of microscopy images. This data was used to engineer spatially-controlled DNA-scaffolded bivalent ligands, which we used to compare trends in spatial-selective binding observed across HUVEC, CHO and HeLa in resting versus activated conditions in confocal microscopy images. For HUVEC and CHO, the data demonstrated an improved selectivity and localisation of binding for smaller spacings ~7 nm and ~24 nm, in good agreement with the model. A deviation from the mode predictions for HeLa was observed, indicative of a clustered, instead of homogeneous, integrin organization. Our findings demonstrate how low-technology imaging methods can guide the design of spatially controlled ligands to selectively differentiate between cell type and integrin activation state.

Stability and Existence of Noncanonical I-motif DNA Structures in Computer Simulations Based on Atomistic and Coarse-Grained Force Fields

Cytosine-rich DNA sequences are able to fold into noncanonical structures, in which semi-protonated cytosine pairs develop extra hydrogen bonds, and these bonds are responsible for the overall stability of a structure called the i-motif. The i-motif can be formed in many regions of the genome, but the most representative is the telomeric region in which the CCCTAA sequences are repeated thousands of times. The ability to reverse folding/unfolding in response to pH change makes the above sequence and i-motif very promising components of nanomachines, extended DNA structures, and drug carriers. Molecular dynamics analysis of such structures is highly beneficial due to direct insights into the microscopic structure of the considered systems. We show that Amber force fields for DNA predict the stability of the i-motif over a long timescale; however, these force fields are not able to predict folding of the cytosine-rich sequences into the i-motif. The reason is the kinetic partitioning of the folding process, which makes the transitions between various intermediates too time-consuming in atomistic force field representation. Application of coarse-grained force fields usually highly accelerates complex structural transitions. We, however, found that three of the most popular coarse-grained force fields for DNA (oxDNA, 3SPN, and Martini) were not able to predict the stability of the i-motif structure. Obviously, they were not able to accelerate the folding of unfolded states into an i-motif. This observation must be strongly highlighted, and the need to develop suitable extensions of coarse-grained force fields for DNA is pointed out. However, it will take a great deal of effort to successfully solve these problems.

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.

Probing the Mechanical Properties of DNA Nanostructures with Metadynamics

Molecular dynamics simulations are often used to provide feedback in the design workflow of DNA nanostructures. However, even with coarse-grained models, the convergence of distributions from unbiased simulation is slow, limiting applications to equilibrium structural properties. Given the increasing interest in dynamic, reconfigurable, and deformable devices, methods that enable efficient quantification of large ranges of motion, conformational transitions, and mechanical deformation are critically needed. Metadynamics is an automated biasing technique that enables the rapid acquisition of molecular conformational distributions by flattening free energy landscapes. Here we leveraged this approach to sample the free energy landscapes of DNA nanostructures whose unbiased dynamics are nonergodic, including bistable Holliday junctions and part of a bistable DNA origami structure. Taking a DNA origami-compliant joint as a case study, we further demonstrate that metadynamics can predict the mechanical response of a full DNA origami device to an applied force, showing good agreement with experiments. Our results exemplify the efficient computation of free energy landscapes and force response in DNA nanodevices, which could be applied for rapid feedback in iterative design workflows and generally facilitate the integration of simulation and experiments. Metadynamics will be particularly useful to guide the design of dynamic devices for nanorobotics, biosensing, or nanomanufacturing applications.

Probing the Mechanical Properties of DNA Nanostructures with Metadynamics

Molecular dynamics simulations are often used to provide feedback in the design workflow of DNA nanostructures. However, even with coarse-grained models, the convergence of distributions from unbiased simulation is slow, limiting applications to equilibrium structural properties. Given the increasing interest in dynamic, reconfigurable, and deformable devices, methods that enable efficient quantification of large ranges of motion, conformational transitions, and mechanical deformation are critically needed. Metadynamics is an automated biasing technique that enables the rapid acquisition of molecular conformational distributions by flattening free energy landscapes. Here we leveraged this approach to sample the free energy landscapes of DNA nanostructures whose unbiased dynamics are nonergodic, including bistable Holliday junctions and part of a bistable DNA origami structure. Taking a DNA origami-compliant joint as a case study, we further demonstrate that metadynamics can predict the mechanical response of a full DNA origami device to an applied force, showing good agreement with experiments. Our results exemplify the efficient computation of free energy landscapes and force response in DNA nanodevices, which could be applied for rapid feedback in iterative design workflows and generally facilitate the integration of simulation and experiments. Metadynamics will be particularly useful to guide the design of dynamic devices for nanorobotics, biosensing, or nanomanufacturing applications.

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.

Design and simulation of DNA, RNA and hybrid protein-nucleic acid nanostructures with oxView

Molecular simulation has become an integral part of the DNA/RNA nanotechnology research pipeline. In particular, understanding the dynamics of structures and single-molecule events has improved the precision of nanoscaffolds and diagnostic tools. Here we present oxView, a design tool for visualization, design, editing and analysis of simulations of DNA, RNA and nucleic acid-protein nanostructures. oxView provides an accessible software platform for designing novel structures, tweaking existing designs, preparing them for simulation in the oxDNA/RNA molecular simulation engine and creating visualizations of simulation results. In several examples, we present procedures for using the tool, including its advanced features that couple the design capabilities with a coarse-grained simulation engine and scripting interface that can programmatically edit structures and facilitate design of complex structures from multiple substructures. These procedures provide a practical basis from which researchers, including experimentalists with limited computational experience, can integrate simulation and 3D visualization into their existing research programs.

Salt-induced conformational switching of a flat rectangular DNA origami structure

A rectangular DNA origami structure is one of the most studied and often used motif for applications in DNA nanotechnology. Here, we present two assays to study structural changes in DNA nanostructures and reveal a reversible rolling-up of the rectangular DNA origami structure induced by bivalent cations such as magnesium or calcium. First, we applied one-color and two-color superresolution DNA-PAINT with protruding strands along the long edges of the DNA origami rectangle. At increasing salt concentration, a single line instead of two lines is observed as a first indicator of rolling-up. Two-color measurements also revealed different conformations with parallel and angled edges. Second, we placed a gold nanoparticle and a dye molecule at different positions on the DNA origami structure. Distance dependent fluorescence quenching by the nanoparticle reports on dynamic transitions as well as it provides evidence that the rolling-up occurs preferentially along the diagonal of the DNA origami rectangle. The results will be helpful to test DNA structural models and the assays presented will be useful to study further structural transitions in DNA nanotechnology.

Dataset of normalized probability distributions of virtual bond lengths, bond angles, and dihedral angles for the coarse-grained single-stranded DNA structures

The utility of the coarse-grained (CG) single-stranded DNA (ssDNA) model can drastically reduce the compute time for simulating the ssDNA dynamics. The model-matched CG potentials and the inherent potential constants can be derived by coarse-graining the experimentally measured ssDNA structures. A useful and widespread treatment of the CG model is to use three different pseudo-atoms P, S, and B to represent the atomic groups of phosphate, sugar, and base, respectively, in each nucleotide of the ssDNA structures. The three pseudo-atoms generate nine types of the structural parameters to characterize the unstructured ssDNA conformations, including three (virtual) bond lengths (P-S, S-B, and S-P) between two neighbouring beads, four bond angles (P-S-P, S-P-S, P-S-B, and B-S-P) between three adjacent bonds, and two dihedral angles (P-S-P-S and S-P-S-P) between three successive bonds. This paper mainly presents the data of normalized probability distributions of the bond lengths, bond angles, and dihedral angles for the CG ssDNAs.

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.

Environment-Dependent Stability and Mechanical Properties of DNA Origami Six-Helix Bundles with Different Crossover Spacings

The internal design of DNA nanostructures defines how they behave in different environmental conditions, such as endonuclease-rich or low-Mg²⁺ solutions. Notably, the inter-helical crossovers that form the core of such DNA objects have a major impact on their mechanical properties and stability. Importantly, crossover design can be used to optimize DNA nanostructures for target applications, especially when developing them for biomedical environments. To elucidate this, two otherwise identical DNA origami designs are presented that have a different number of staple crossovers between neighboring helices, spaced at 42- and 21- basepair (bp) intervals, respectively. The behavior of these structures is then compared in various buffer conditions, as well as when they are exposed to enzymatic digestion by DNase I. The results show that an increased number of crossovers significantly improves the nuclease resistance of the DNA origami by making it less accessible to digestion enzymes but simultaneously lowers its stability under Mg²⁺ -free conditions by reducing the malleability of the structures. Therefore, these results represent an important step toward rational, application-specific DNA nanostructure design.

OxDNA to Study Species Interactions

Molecular ecology uses molecular genetic data to answer traditional ecological questions in biogeography and biodiversity, among others. Several ecological principles, such as the niche hypothesis and the competitive exclusions, are based on the fact that species compete for resources. More in generally, it is now recognized that species interactions play a crucial role in determining the coexistence and abundance of species. However, experimentally controllable platforms, which allow us to study and measure competitions among species, are rare and difficult to implement. In this work, we suggest exploiting a Molecular Dynamics coarse-grained model to study interactions among single strands of DNA, representing individuals of different species, which compete for binding to other oligomers considered as resources. In particular, the well-established knowledge of DNA–DNA interactions at the nanoscale allows us to test the hypothesis that the maximum consecutive overlap between pairs of oligomers measure the species’ competitive advantages. However, we suggest that a more complex structure also plays a role in the ability of the species to successfully bind to the target resource oligomer. We complement the simulations with experiments on populations of DNA strands which qualitatively confirm our hypotheses. These tools constitute a promising starting point for further developments concerning the study of controlled, DNA-based, artificial ecosystems.

A localized DNA finite-state machine with temporal resolution

The identity and timing of environmental stimulus play a pivotal role in living organisms in programming their signaling networks and developing specific phenotypes. The ability to unveil history-dependent signals will advance our understanding of temporally regulated biological processes. Here, we have developed a two-input, five-state DNA finite-state machine (FSM) to sense and record the temporally ordered inputs. The spatial organization of the processing units on DNA origami enables facile modulation of the energy landscape of DNA strand displacement reactions, allowing precise control of the reactions along predefined paths for different input orders. The use of spatial constraints brings about a simple, modular design for the FSM with a minimum set of orthogonal components and confers minimized leaky reactions and fast kinetics. The FSM demonstrates the capability of sensing the temporal orders of two microRNAs, highlighting its potential for temporally resolved biosensing and smart therapeutics.

Predicting the Free-Form Shape of Structured DNA Assemblies from Their Lattice-Based Design Blueprint

Structured DNA assemblies have been designed primarily on a three-dimensional lattice because it is easy to arrange and cross-link the helices there. However, when we design free-form structures including wireframes and topologically closed circular objects on a lattice, artificially stretched bonds connecting bases are inevitably and arbitrarily formed. They often lead to nonconvergence or convergence to a wrong configuration in computational analysis to predict the equilibrium shape of the structure when started from its lattice-based configuration, which hinders the design process of free-form structures. Here, we present a computational procedure enabling the shape prediction of free-form structures from their lattice-based design blueprint without any convergence issue. It automatically partitions the structure into substructures and relocates them into a new configuration. When the analysis for calculating the equilibrium shape begins from this configuration, no convergence issue occurs because substructures and stretched bonds connecting them do not overlap and intertwine each other during analysis. Using the proposed approach, we could obtain the free-form shape of a comprehensive set of wireframe and circular structures accurately and quickly. We further demonstrated that it also facilitated a design of wireframe structures with nonstraight edges.

The interplay of supercoiling and thymine dimers in DNA

Thymine dimers are a major mutagenic photoproduct induced by UV radiation. While they have been the subject of extensive theoretical and experimental investigations, questions of how DNA supercoiling affects local defect properties, or, conversely, how the presence of such defects changes global supercoiled structure, are largely unexplored. Here, we introduce a model of thymine dimers in the oxDNA forcefield, parametrized by comparison to melting experiments and structural measurements of the thymine dimer induced bend angle. We performed extensive molecular dynamics simulations of double-stranded DNA as a function of external twist and force. Compared to undamaged DNA, the presence of a thymine dimer lowers the supercoiling densities at which plectonemes and bubbles occur. For biologically relevant supercoiling densities and forces, thymine dimers can preferentially segregate to the tips of the plectonemes, where they enhance the probability of a localized tip-bubble. This mechanism increases the probability of highly bent and denatured states at the thymine dimer site, which may facilitate repair enzyme binding. Thymine dimer-induced tip-bubbles also pin plectonemes, which may help repair enzymes to locate damage. We hypothesize that the interplay of supercoiling and local defects plays an important role for a wider set of DNA damage repair systems.

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.

Multi-resolution simulation of DNA transport through large synthetic nanostructures

Modeling and simulation has become an invaluable partner in development of nanopore sensing systems. The key advantage of the nanopore sensing method - the ability to rapidly detect individual biomolecules as a transient reduction of the ionic current flowing through the nanopore - is also its key deficiency, as the current signal itself rarely provides direct information about the chemical structure of the biomolecule. Complementing experimental calibration of the nanopore sensor readout, coarse-grained and all-atom molecular dynamics simulations have been used extensively to characterize the nanopore translocation process and to connect the microscopic events taking place inside the nanopore to the experimentally measured ionic current blockades. Traditional coarse-grained simulations, however, lack the precision needed to predict ionic current blockades with atomic resolution whereas traditional all-atom simulations are limited by the length and time scales amenable to the method. Here, we describe a multi-resolution framework for modeling electric field-driven passage of DNA molecules and nanostructures through to-scale models of synthetic nanopore systems. We illustrate the method by simulating translocation of double-stranded DNA through a solid-state nanopore and a micron-scale slit, capture and translocation of single-stranded DNA in a double nanopore system, and modeling ionic current readout from a DNA origami nanostructure passage through a nanocapillary. We expect our multi-resolution simulation framework to aid development of the nanopore field by providing accurate, to-scale modeling capability to research laboratories that do not have access to leadership supercomputer facilities.

Computational investigation of the impact of core sequence on immobile DNA four-way junction structure and dynamics

Immobile four-way junctions (4WJs) are core structural motifs employed in the design of programmed DNA assemblies. Understanding the impact of sequence on their equilibrium structure and flexibility is important to informing the design of complex DNA architectures. While core junction sequence is known to impact the preferences for the two possible isomeric states that junctions reside in, previous investigations have not quantified these preferences based on molecular-level interactions. Here, we use all-atom molecular dynamics simulations to investigate base-pair level structure and dynamics of four-way junctions, using the canonical Seeman J1 junction as a reference. Comparison of J1 with equivalent single-crossover topologies and isolated nicked duplexes reveal conformational impact of the double-crossover motif. We additionally contrast J1 with a second junction core sequence termed J24, with equal thermodynamic preference for each isomeric configuration. Analyses of the base-pair degrees of freedom for each system, free energy calculations, and reduced-coordinate sampling of the 4WJ isomers reveal the significant impact base sequence has on local structure, isomer bias, and global junction dynamics.

Nanobase.org: a repository for DNA and RNA nanostructures

We introduce a new online database of nucleic acid nanostructures for the field of DNA and RNA nanotechnology. The database implements an upload interface, searching and database browsing. Each deposited nanostructures includes an image of the nanostructure, design file, an optional 3D view, and additional metadata such as experimental data, protocol or literature reference. The database accepts nanostructures in any preferred format used by the uploader for the nanostructure design. We further provide a set of conversion tools that encourage design file conversion into common formats (oxDNA and PDB) that can be used for setting up simulations, interactive editing or 3D visualization. The aim of the repository is to provide to the DNA/RNA nanotechnology community a resource for sharing their designs for further reuse in other systems and also to function as an archive of the designs that have been achieved in the field so far. Nanobase.org is available at https://nanobase.org/.

Unified Nanotechnology Format: One Way to Store Them All

The domains of DNA and RNA nanotechnology are steadily gaining in popularity while proving their value with various successful results, including biosensing robots and drug delivery cages. Nowadays, the nanotechnology design pipeline usually relies on computer-based design (CAD) approaches to design and simulate the desired structure before the wet lab assembly. To aid with these tasks, various software tools exist and are often used in conjunction. However, their interoperability is hindered by a lack of a common file format that is fully descriptive of the many design paradigms. Therefore, in this paper, we propose a Unified Nanotechnology Format (UNF) designed specifically for the biomimetic nanotechnology field. UNF allows storage of both design and simulation data in a single file, including free-form and lattice-based DNA structures. By defining a logical and versatile format, we hope it will become a widely accepted and used file format for the nucleic acid nanotechnology community, facilitating the future work of researchers and software developers. Together with the format description and publicly available documentation, we provide a set of converters from existing file formats to simplify the transition. Finally, we present several use cases visualizing example structures stored in UNF, showcasing the various types of data UNF can handle.

Massively Parallelized Molecular Force Manipulation with On-Demand Thermal and Optical Control

In single-molecule force spectroscopy (SMFS), a tethered molecule is stretched using a specialized instrument to study how macromolecules extend under force. One problem in SMFS is the serial and slow nature of the measurements, performed one molecule at a time. To address this long-standing challenge, we report on the origami polymer force clamp (OPFC) which enables parallelized manipulation of the mechanical forces experienced by molecules without the need for dedicated SMFS instruments or surface tethering. The OPFC positions target molecules between a rigid nanoscale DNA origami beam and a responsive polymer particle that shrinks on demand. As a proof-of-concept, we record the steady state and time-resolved mechanical unfolding dynamics of DNA hairpins using the fluorescence signal from ensembles of molecules and confirm our conclusion using modeling.

DNA sequence and methylation prescribe the inside-out conformational dynamics and bending energetics of DNA minicircles

Eukaryotic genome and methylome encode DNA fragments' propensity to form nucleosome particles. Although the mechanical properties of DNA possibly orchestrate such encoding, the definite link between 'omics' and DNA energetics has remained elusive. Here, we bridge the divide by examining the sequence-dependent energetics of highly bent DNA. Molecular dynamics simulations of 42 intact DNA minicircles reveal that each DNA minicircle undergoes inside-out conformational transitions with the most likely configuration uniquely prescribed by the nucleotide sequence and methylation of DNA. The minicircles' local geometry consists of straight segments connected by sharp bends compressing the DNA's inward-facing major groove. Such an uneven distribution of the bending stress favors minimum free energy configurations that avoid stiff base pair sequences at inward-facing major grooves. Analysis of the minicircles' inside-out free energy landscapes yields a discrete worm-like chain model of bent DNA energetics that accurately account for its nucleotide sequence and methylation. Experimentally measuring the dependence of the DNA looping time on the DNA sequence validates the model. When applied to a nucleosome-like DNA configuration, the model quantitatively reproduces yeast and human genomes' nucleosome occupancy. Further analyses of the genome-wide chromatin structure data suggest that DNA bending energetics is a fundamental determinant of genome architecture.

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.