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

Encoding signal propagation on topology-programmed DNA origami

Biological systems often rely on topological transformation to reconfigure connectivity between nodes to guide the flux of molecular information. Here we develop a topology-programmed DNA origami system that encodes signal propagation at the nanoscale, analogous to topologically efficient information processing in cellular systems. We present a systematic molecular implementation of topological operations involving 'glue-cut' processes that can prompt global conformational change of DNA origami structures, with demonstrated major topological properties including genus, number of boundary components and orientability. By spatially arranging reactive DNA hairpins, we demonstrate signal propagation across transmission paths of varying lengths and orientations, and curvatures on the curved surfaces of three-dimensional origamis. These DNA origamis can also form dynamic scaffolds for regulating the spatial and temporal signal propagations whereby topological transformations spontaneously alter the location of nodes and boundary of signal propagation network. We anticipate that our strategy for topological operations will provide a general route to manufacture dynamic DNA origami nanostructures capable of performing global structural transformations under programmable control.

Simulating the Lyotropic Phase Behavior of a Partially Self-Complementary DNA Tetramer

DNA oligomers in solution have been found to develop liquid crystal phases via a hierarchical process that involves Watson-Crick base pairing, supramolecular assembly into columns of duplexes, and long-range ordering. The multiscale nature of this phenomenon makes it difficult to quantitatively describe and assess the importance of the various contributions, particularly for very short strands. We performed molecular dynamics simulations based on the coarse-grained oxDNA model, aiming to depict all of the assembly processes involved and the phase behavior of solutions of the DNA GCCG tetramers. We find good quantitative matching to experimental data at both levels of molecular association (thermal melting) and collective ordering (phase diagram). We characterize the isotropic state and the low-density nematic and high-density columnar liquid crystal phases in terms of molecular order, size of aggregates, and structure, together with their effects on diffusivity processes. We observe a cooperative aggregation mechanism in which the formation of dimers is less thermodynamically favored than the formation of longer aggregates.

ANNaMo: Coarse-grained modeling for folding and assembly of RNA and DNA systems

The folding of RNA and DNA strands plays crucial roles in biological systems and bionanotechnology. However, studying these processes with high-resolution numerical models is beyond current computational capabilities due to the timescales and system sizes involved. In this article, we present a new coarse-grained model for investigating the folding dynamics of nucleic acids. Our model represents three nucleotides with a patchy particle and is parameterized using well-established nearest-neighbor models. Thanks to the reduction of degrees of freedom and to a bond-swapping mechanism, our model allows for simulations at timescales and length scales that are currently inaccessible to more detailed models. To validate the performance of our model, we conducted extensive simulations of various systems: We examined the thermodynamics of DNA hairpins, capturing their stability and structural transitions, the folding of an MMTV pseudoknot, which is a complex RNA structure involved in viral replication, and also explored the folding of an RNA tile containing a k-type pseudoknot. Finally, we evaluated the performance of the new model in reproducing the melting temperatures of oligomers and the dependence on the toehold length of the displacement rate in toehold-mediated displacement processes, a key reaction used in molecular computing. All in all, the successful reproduction of experimental data and favorable comparisons with existing coarse-grained models validate the effectiveness of the new model.

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

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

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.

Bending Unwinds DNA

In biology, DNA is often tightly bent to small radii. Solely based on the groove asymmetry, a 30-year-old theoretical paper predicted that such bending should unwind DNA, but this effect has not been directly experimentally quantified so far. We developed a ligation-based assay with nicked DNA circles of variable length, thereby decoupling the twist-dependent ligation efficiency from the large bending strain which dominates conventional circularization assays. We demonstrate that tightly bent DNA indeed unwinds to over 11 base pairs/turn, exactly as predicted. Our discovery requires reassessing the molecular mechanisms and energetics of all processes where DNA is tightly bent or relaxed again, including DNA packaging, gene regulation and expression.

A Novel Fractional Brownian Dynamics Method for Simulating the Dynamics of Confined Bottle-Brush Polymers in Viscoelastic Solution

In crowded fluids, polymer segments can exhibit anomalous subdiffusion due to the viscoelasticity of the surrounding environment. Previous single-particle tracking experiments revealed that such anomalous diffusion in complex fluids (e.g., in bacterial cytoplasm) can be described by fractional Brownian motion (fBm). To investigate how the viscoelastic media affects the diffusive behaviors of polymer segments without resolving single crowders, we developed a novel fractional Brownian dynamics method to simulate the dynamics of polymers under confinement. In this work, instead of using Gaussian random numbers ("white Gaussian noise") to model the Brownian force as in the standard Brownian dynamics simulations, we introduce fractional Gaussian noise (fGn) in our homemade fractional Brownian dynamics simulation code to investigate the anomalous diffusion of polymer segments by using a simple "bottle-brush"-type polymer model. The experimental results of the velocity autocorrelation function and the exponent that characterizes the subdiffusion of the confined polymer segments can be reproduced by this simple polymer model in combination with fractional Gaussian noise (fGn), which mimics the viscoelastic media.

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.

Ionic Liquid-Induced Assembly of DNA at Air-Water Interface

DNA nanotechnology is the future of many products in the pharmaceutical and cosmetic industries. Self-assembly of this negatively charged biopolymer at surfaces and interfaces is an essential step to elaborate its field of applications. In this study, the ionic liquid (IL) monolayer-assisted self-assembly of DNA macromolecules at the air-water interface has been closely monitored by employing various quantitative techniques, namely, surface pressure-area (π-A) isotherms, surface potential, interfacial rheology, and X-ray reflectivity (XRR). The π-A isotherms reveal that the IL 1,3-didecyl 3-methyl imidazolium chloride induces DNA self-assembly at the interface, leading to a thick viscoelastic film. The interfacial rheology exhibits a notable rise in the viscoelastic modulus as the surface pressure increases. The values of storage and loss moduli measured as a function of strain frequency suggest a relaxation frequency that depends on the length of the macromolecule. The XRR measurements indicate a considerable increase in DNA layer thickness at the elevated surface pressures depending on the number of base pairs of the DNA. The results are considered in terms of the electrostatic and hydrophobic interactions, allowing a quantitative conclusion about the arrangement of DNA strands underneath the monolayer of the ILs at the air-water interface.

Chemo-mechanical forces modulate the topology dynamics of mesoscale DNA assemblies

The intrinsic complexity of many mesoscale (10-100 nm) cellular machineries makes it challenging to elucidate their topological arrangement and transition dynamics. Here, we exploit DNA origami nanospring as a model system to demonstrate that tens of piconewton linear force can modulate higher-order conformation dynamics of mesoscale molecular assemblies. By switching between two chemical structures (i.e., duplex and tetraplex DNA) in the junctions of adjacent origami modules, the corresponding stretching or compressing chemo-mechanical stress reversibly flips the backbone orientations of the DNA nanosprings. Both coarse-grained molecular dynamics simulations and atomic force microscopy measurements reveal that such a backbone conformational switch does not alter the right-handed chirality of the nanospring helix. This result suggests that mesoscale helical handedness may be governed by the torque, rather than the achiral orientation, of nanospring backbones. It offers a topology-based caging/uncaging concept to present chemicals in response to environmental cues in solution.

Asymmetric Periodic Boundary Conditions for All-Atom Molecular Dynamics and Coarse-Grained Simulations of Nucleic Acids

Periodic boundary conditions are commonly applied in molecular dynamics simulations in the microcanonical (NVE), canonical (NVT), and isothermal-isobaric (NpT) ensembles. In their simplest application, a biological system of interest is placed in the middle of a solvation box, which is chosen 'sufficiently large' to minimize any numerical artifacts associated with the periodic boundary conditions. This practical approach brings limitations to the size of biological systems that can be simulated. Here, we study simulations of effectively infinitely long nucleic acids, which are solvated in the directions perpendicular to the polymer chain, while periodic boundary conditions are also applied along the polymer chain. We study the effects of these asymmetric periodic boundary conditions (APBC) on the simulated results, including the mechanical properties of biopolymers and the properties of the surrounding solvent. To get some further insights into the advantages of using the APBC, a coarse-grained worm-like chain model is first studied, illustrating how the persistence length can be extracted from the local properties of the polymer chain, which are less affected by the APBC than some global averages. This is followed by all-atom molecular dynamics simulations of DNA in ionic solutions, where we use the APBC to investigate sequence-dependent properties of DNA molecules and properties of the surrounding solvent.

Kinetics of radiation-induced DNA double-strand breaks through coarse-grained simulations

Double-strand breaks (DSBs), i.e., the covalent cut of the DNA backbone over both strands, are a detrimental outcome of cell irradiation, bearing chromosomal aberrations and leading to cell apoptosis. In the early stages of the evolution of a DSB, the disruption of the residual interactions between the DNA moieties drives the fracture of the helical layout; in spite of its biological significance, the details of this process are still largely uncertain. Here, we address the mechanical rupture of DNA by DSBs via coarse-grained molecular dynamics simulations: the setup involves a 3855-bp DNA filament and diverse DSB motifs, i.e., within a range of distances between strand breaks (or DSB distance). By employing a coarse-grained model of DNA, we access the molecular details and characteristic timescales of the rupturing process. A sequence-nonspecific, linear correlation is observed between the DSB distance and the internal energy contribution to the disruption of the residual (Watson-Crick and stacking) contacts between DNA moieties, which is seemingly driven by an abrupt, cooperative process. Moreover, we infer an exponential dependence of the characteristic rupture times on the DSB distances, which we associate to an Arrhenius-like law of thermally-activated processes. This work lays the foundations of a detailed, mechanistic assessment of DSBs in silico as a benchmark to both numerical simulations and data from single-molecule experiments.

Dynamics of chromosome organization in a minimal bacterial cell

Computational models of cells cannot be considered complete unless they include the most fundamental process of life, the replication and inheritance of genetic material. By creating a computational framework to model systems of replicating bacterial chromosomes as polymers at 10 bp resolution with Brownian dynamics, we investigate changes in chromosome organization during replication and extend the applicability of an existing whole-cell model (WCM) for a genetically minimal bacterium, JCVI-syn3A, to the entire cell-cycle. To achieve cell-scale chromosome structures that are realistic, we model the chromosome as a self-avoiding homopolymer with bending and torsional stiffnesses that capture the essential mechanical properties of dsDNA in Syn3A. In addition, the conformations of the circular DNA must avoid overlapping with ribosomes identitied in cryo-electron tomograms. While Syn3A lacks the complex regulatory systems known to orchestrate chromosome segregation in other bacteria, its minimized genome retains essential loop-extruding structural maintenance of chromosomes (SMC) protein complexes (SMC-scpAB) and topoisomerases. Through implementing the effects of these proteins in our simulations of replicating chromosomes, we find that they alone are sufficient for simultaneous chromosome segregation across all generations within nested theta structures. This supports previous studies suggesting loop-extrusion serves as a near-universal mechanism for chromosome organization within bacterial and eukaryotic cells. Furthermore, we analyze ribosome diffusion under the influence of the chromosome and calculate in silico chromosome contact maps that capture inter-daughter interactions. Finally, we present a methodology to map the polymer model of the chromosome to a Martini coarse-grained representation to prepare molecular dynamics models of entire Syn3A cells, which serves as an ultimate means of validation for cell states predicted by the WCM.

Reverse engineering DNA origami nanostructure designs from raw scaffold and staple sequence lists

Designs for scaffolded DNA origami nanostructures are commonly and minimally published as the list of DNA staple and scaffold sequences required. In nearly all cases, high-level editable design files (e.g. caDNAno) which generated the low-level sequences are not made available. This de facto 'raw sequence' exchange format allows published origami designs to be re-attempted in the laboratory by other groups, but effectively stops designs from being significantly modified or re-purposed for new future applications. To make the raw sequence exchange format more accessible to further design and engineering, in this work we propose the first algorithmic solution to the inverse problem of converting staple/scaffold sequences back to a 'guide schematic' resembling the original origami schematic. The guide schematic can be used to aid the manual re-input of an origami into a CAD tool like caDNAno, hence recovering a high-level editable design file. Creation of a guide schematic can also be used to double check that a list of staple strand sequences does not have errors and indeed does assemble into a desired origami nanostructure prior to costly laboratory experimentation. We tested our reverse algorithm on 36 diverse origami designs from the literature and found that 29 origamis (81 %) had a good quality guide schematic recovered from raw sequences. Our software is made available at https://revnano.readthedocs.io.

Anomalous Diffusion of Polyelectrolyte Segments on Supported Charged Lipid Bilayers

This work provides mesoscale models for the anomalous diffusion of a polymer chain on a heterogeneous surface with rearranging randomly distributed adsorption sites. Both the "bead-spring" model and oxDNA model were simulated on supported lipid bilayer membranes with various molar fractions of charged lipids, using Brownian dynamics method. Our simulation results demonstrate that "bead-spring" chains exhibit sub-diffusion on charged lipid bilayers which agrees with previous experimental observations for short-time dynamics of DNA segments on membranes. In addition, the non-Gaussian diffusive behaviors of DNA segments have not been observed in our simulations. However, a simulated 17 base pairs double stranded DNA, using oxDNA model, performs normal diffusion on supported cationic lipid bilayers. Due to the number of positively charged lipids attracted by short DNA is small, the energy landscape that the short DNA experiences during diffusion is not as heterogeneous as that experienced by long DNA chains, which results in normal diffusion rather than sub-diffusion for short DNA.

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.

Programmed Self-Assembly of DNA Nanosheets with Discrete Single-Molecule Thickness and Interfacial Mechanics: Design, Simulation, and Characterization

DNA is programmed to hierarchically self-assemble into superstructures spanning from nanometer to micrometer scales. Here, we demonstrate DNA nanosheets assembled out of a rationally designed flexible DNA unit (F-unit), whose shape resembles a Feynman diagram. F-units were designed to self-assemble in two dimensions and to display a high DNA density of hydrophobic moieties. oxDNA simulations confirmed the planarity of the F-unit. DNA nanosheets with a thickness of a single DNA duplex layer and with large coverage (at least 30 μm × 30 μm) were assembled from the liquid phase at the solid/liquid interface, as unambiguously evidenced by atomic force microscopy imaging. Interestingly, single-layer nanodiscs formed in solution at low DNA concentrations. DNA nanosheet superstructures were further assembled at liquid/liquid interfaces, as demonstrated by the fluorescence of a double-stranded DNA intercalator. Moreover, the interfacial mechanical properties of the nanosheet superstructures were measured as a response to temperature changes, demonstrating the control of interfacial shear mechanics based on DNA nanostructure engineering. The rational design of the F-unit, along with the presented results, provide an avenue toward the controlled assembly of reconfigurable/responsive nanosheets and membranes at liquid/liquid interfaces, to be potentially used in the characterization of biomechanical processes and materials transport.

Self-Assembled Artificial DNA Nanocompartments and Their Bioapplications

Compartmentalization is the strategy evolved by nature to control reactions in space and time. The ability to emulate this strategy through synthetic compartmentalization systems has rapidly evolved in the past years, accompanied by an increasing understanding of the effects of spatial confinement on the thermodynamic and kinetic properties of the guest molecules. DNA nanotechnology has played a pivotal role in this scientific endeavor and is still one of the most promising approaches for the construction of nanocompartments with programmable structural features and nanometer-scaled addressability. In this review, the design approaches, bioapplications, and theoretical frameworks of self-assembled DNA nanocompartments are surveyed. From DNA polyhedral cages to virus-like capsules, the construction principles of such intriguing architectures are illustrated. Various applications of DNA nanocompartments, including their use for programmable enzyme scaffolding, single-molecule studies, biosensing, and as artificial nanofactories, ending with an ample description of DNA nanocages for biomedical purposes, are then reported. Finally, the theoretical hypotheses that make DNA nanocompartments, and nanosystems in general, a topic of great interest in modern science, are described and the progresses that have been done until now in the comprehension of the peculiar phenomena that occur within nanosized environments are summarized.

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.

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.

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.

Integrating CRISPR/Cas systems with programmable DNA nanostructures for delivery and beyond

Precise genome editing with CRISPR/Cas paves the way for many biochemical, biotechnological, and medical applications, and consequently, it may enable treatment of already known and still-to-be-found genetic diseases. Meanwhile, another rapidly emerging field—structural DNA nanotechnology—provides a customizable and modular platform for accurate positioning of nanoscopic materials, for e.g., biomedical uses. This addressability has just recently been applied in conjunction with the newly developed gene engineering tools to enable impactful, programmable nanotechnological applications. As of yet, self-assembled DNA nanostructures have been mainly employed to enhance and direct the delivery of CRISPR/Cas, but lately the groundwork has also been laid out for other intriguing and complex functions. These recent advances will be described in this perspective.

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.