Characterizing DNA Star-Tile-Based Nanostructures Using a Coarse-Grained Model

Comparative Simulative Analysis and Design of Single-Chain Self-Assembled Protein Cages

Coiled-coil protein origami (CCPO) is a modular strategy for the de novo design of polypeptide nanostructures. It represents a type of modular design based on pairwise-interacting coiled-coil (CC) units with a single-chain protein programmed to fold into a polyhedral cage. However, the mechanisms underlying the self-assembly of the protein tetrahedron are still not fully understood. In the present study, 18 CCPO cages with three different topologies were modeled in silico. Then, molecular dynamics simulations and CC parameters were calculated to characterize the dynamic properties of protein tetrahedral cages at both the local and global levels. Furthermore, a deformed CC unit was redesigned, and the stability of the new cage was significantly improved.

Mechanisms of Biomolecular Self-Assembly Investigated Through In Situ Observations of Structures and Dynamics

Biomolecular self-assembly of hierarchical materials is a precise and adaptable bottom-up approach to synthesizing across scales with considerable energy, health, environment, sustainability, and information technology applications. To achieve desired functions in biomaterials, it is essential to directly observe assembly dynamics and structural evolutions that reflect the underlying energy landscape and the assembly mechanism. This review will summarize the current understanding of biomolecular assembly mechanisms based on in situ characterization and discuss the broader significance and achievements of newly gained insights. In addition, we will also introduce how emerging deep learning/machine learning-based approaches, multiparametric characterization, and high-throughput methods can boost the development of biomolecular self-assembly. The objective of this review is to accelerate the development of in situ characterization approaches for biomolecular self-assembly and to inspire the next generation of biomimetic materials.

Structure and dynamics of an archetypal DNA nanoarchitecture revealed via cryo-EM and molecular dynamics simulations

DNA can be folded into rationally designed, unique, and functional materials. To fully realise the potential of these DNA materials, a fundamental understanding of their structure and dynamics is necessary, both in simple solvents as well as more complex and diverse anisotropic environments. Here we analyse an archetypal six-duplex DNA nanoarchitecture with single-particle cryo-electron microscopy and molecular dynamics simulations in solvents of tunable ionic strength and within the anisotropic environment of biological membranes. Outside lipid bilayers, the six-duplex bundle lacks the designed symmetrical barrel-type architecture. Rather, duplexes are arranged in non-hexagonal fashion and are disorted to form a wider, less elongated structure. Insertion into lipid membranes, however, restores the anticipated barrel shape due to lateral duplex compression by the bilayer. The salt concentration has a drastic impact on the stability of the inserted barrel-shaped DNA nanopore given the tunable electrostatic repulsion between the negatively charged duplexes. By synergistically combining experiments and simulations, we increase fundamental understanding into the environment-dependent structural dynamics of a widely used nanoarchitecture. This insight will pave the way for future engineering and biosensing applications.

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.

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

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

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.

Simulative Analysis of a Family of DNA Tetrahedrons Produced by Changing the Twisting Number of Each Double Helix

In this paper, three tetrahedral nanocages, composed of six DNA double helix edges with all having the twist number 1, 2 or 3, have been characterized using classical molecular dynamics simulation to measure the specific structural and conformational features produced by only changing the twisting number of each double helix. The simulation result indicates that three tetrahedral cages are relatively stable and are maintained along the entire trajectory. Each double helix is more inclined to behave as a whole in the 2TD and 3TD cages than in the 1TD cage according to the cross-correlation maps for three nanocages, and also their local motions are more easily induced by the conformational variability of the thymidine linkers due to the increased flexibility of each helix. Hence, the double helices become the important factors on the structural stability of total cages with the DNA twisting number, and also give the signification contributions to the sizes of these cages conferring the larger spaces of the 2TD and 3TD cages than the 1TD cage. Our result provides an insight into which roles the double helix edges play in assembling DNA polyhedron, and also contribute to improving the loading capacity of DNA tetrahedron in drug delivery.

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

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

The Art of Designing DNA Nanostructures with CAD Software

Since the arrival of DNA nanotechnology nearly 40 years ago, the field has progressed from its beginnings of envisioning rather simple DNA structures having a branched, multi-strand architecture into creating beautifully complex structures comprising hundreds or even thousands of unique strands, with the possibility to exactly control the positions down to the molecular level. While the earliest construction methodologies, such as simple Holliday junctions or tiles, could reasonably be designed on pen and paper in a short amount of time, the advent of complex techniques, such as DNA origami or DNA bricks, require software to reduce the time required and propensity for human error within the design process. Where available, readily accessible design software catalyzes our ability to bring techniques to researchers in diverse fields and it has helped to speed the penetration of methods, such as DNA origami, into a wide range of applications from biomedicine to photonics. Here, we review the historical and current state of CAD software to enable a variety of methods that are fundamental to using structural DNA technology. Beginning with the first tools for predicting sequence-based secondary structure of nucleotides, we trace the development and significance of different software packages to the current state-of-the-art, with a particular focus on programs that are open source.

Valency and Binding Affinity Variations Can Regulate the Multilayered Organization of Protein Condensates with Many Components

Biomolecular condensates, which assemble via the process of liquid-liquid phase separation (LLPS), are multicomponent compartments found ubiquitously inside cells. Experiments and simulations have shown that biomolecular condensates with many components can exhibit multilayered organizations. Using a minimal coarse-grained model for interacting multivalent proteins, we investigate the thermodynamic parameters governing the formation of multilayered condensates through changes in protein valency and binding affinity. We focus on multicomponent condensates formed by scaffold proteins (high-valency proteins that can phase separate on their own via homotypic interactions) and clients (proteins recruited to condensates via heterotypic scaffold-client interactions). We demonstrate that higher valency species are sequestered to the center of the multicomponent condensates, while lower valency proteins cluster towards the condensate interface. Such multilayered condensate architecture maximizes the density of LLPS-stabilizing molecular interactions, while simultaneously reducing the surface tension of the condensates. In addition, multilayered condensates exhibit rapid exchanges of low valency proteins in and out, while keeping higher valency proteins-the key biomolecules involved in condensate nucleation-mostly within. We also demonstrate how modulating the binding affinities among the different proteins in a multicomponent condensate can significantly transform its multilayered structure, and even trigger fission of a condensate into multiple droplets with different compositions.

Self-Limiting Polymerization of DNA Origami Subunits with Strain Accumulation

Biology demonstrates how a near infinite array of complex systems and structures at many scales can originate from the self-assembly of component parts on the nanoscale. But to fully exploit the benefits of self-assembly for nanotechnology, a crucial challenge remains: How do we rationally encode well-defined global architectures in subunits that are much smaller than their assemblies? Strain accumulation via geometric frustration is one mechanism that has been used to explain the self-assembly of global architectures in diverse and complex systems a posteriori. Here we take the next step and use strain accumulation as a rational design principle to control the length distributions of self-assembling polymers. We use the DNA origami method to design and synthesize a molecular subunit known as the PolyBrick, which perturbs its shape in response to local interactions via flexible allosteric blocking domains. These perturbations accumulate at the ends of polymers during growth, until the deformation becomes incompatible with further extension. We demonstrate that the key thermodynamic factors for controlling length distributions are the intersubunit binding free energy and the fundamental strain free energy, both which can be rationally encoded in a PolyBrick subunit. While passive polymerization yields geometrical distributions, which have the highest statistical length uncertainty for a given mean, the PolyBrick yields polymers that approach Gaussian length distributions whose variance is entirely determined by the strain free energy. We also show how strain accumulation can in principle yield length distributions that become tighter with increasing subunit affinity and approach distributions with uniform polymer lengths. Finally, coarse-grained molecular dynamics and Monte Carlo simulations delineate and quantify the dominant forces influencing strain accumulation in a molecular system. This study constitutes a fundamental investigation of the use of strain accumulation as a rational design principle in molecular self-assembly.

Understanding DNA interactions in crowded environments with a coarse-grained model

Nucleic acid interactions under crowded environments are of great importance for biological processes and nanotechnology. However, the kinetics and thermodynamics of nucleic acid interactions in a crowded environment remain poorly understood. We use a coarse-grained model of DNA to study the kinetics and thermodynamics of DNA duplex and hairpin formation in crowded environments. We find that crowders can increase the melting temperature of both an 8-mer DNA duplex and a hairpin with a stem of 6-nt depending on the excluded volume fraction of crowders in solution and the crowder size. The crowding induced stability originates from the entropic effect caused by the crowding particles in the system. Additionally, we study the hybridization kinetics of DNA duplex formation and the formation of hairpin stems, finding that the reaction rate k_on is increased by the crowding effect, while k_off is changed only moderately. The increase in k_on mostly comes from increasing the probability of reaching a transition state with one base pair formed. A DNA strand displacement reaction in a crowded environment is also studied with the model and we find that rate of toehold association is increased, with possible applications to speeding up strand displacement cascades in nucleic acid nanotechnology.

Computing the Elastic Mechanical Properties of Rodlike DNA Nanostructures

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

An electrochemical biosensor exploiting binding-induced changes in electron transfer of electrode-attached DNA origami to detect hundred nanometer-scale targets

The specific detection in clinical samples of analytes with dimensions in the tens to hundreds of nanometers, such as viruses and large proteins, would improve disease diagnosis. Detection of these "mesoscale" analytes (as opposed to their nanoscale components), however, is challenging as it requires the simultaneous binding of multiple recognition sites often spaced over tens of nanometers. In response, we have adapted DNA origami, with its unparalleled customizability to precisely display multiple target-binding sites over the relevant length scale, to an electrochemical biosensor platform. Our proof-of-concept employs triangular origami covalently attached to a gold electrode and functionalized with redox reporters. Electrochemical interrogation of this platform successfully monitors mesoscale, target-binding-induced changes in electron transfer in a manner consistent with coarse-grained molecular dynamics simulations. Our approach enables the specific detection of analytes displaying recognition sites that are separated by ∼40 nm, a spacing significantly greater than that achieved in similar sensor architectures employing either antibodies or aptamers.

Coarse-grained modelling of the structural properties of DNA origami

We use the oxDNA coarse-grained model to provide a detailed characterization of the fundamental structural properties of DNA origami, focussing on archetypal 2D and 3D origami. The model reproduces well the characteristic pattern of helix bending in a 2D origami, showing that it stems from the intrinsic tendency of anti-parallel four-way junctions to splay apart, a tendency that is enhanced both by less screened electrostatic interactions and by increased thermal motion. We also compare to the structure of a 3D origami whose structure has been determined by cryo-electron microscopy. The oxDNA average structure has a root-mean-square deviation from the experimental structure of 8.4 Å, which is of the order of the experimental resolution. These results illustrate that the oxDNA model is capable of providing detailed and accurate insights into the structure of DNA origami, and has the potential to be used to routinely pre-screen putative origami designs and to investigate the molecular mechanisms that regulate the properties of DNA origami.

DNA bipedal motor walking dynamics: an experimental and theoretical study of the dependency on step size

We present a detailed coarse-grained computer simulation and single molecule fluorescence study of the walking dynamics and mechanism of a DNA bipedal motor striding on a DNA origami. In particular, we study the dependency of the walking efficiency and stepping kinetics on step size. The simulations accurately capture and explain three different experimental observations. These include a description of the maximum possible step size, a decrease in the walking efficiency over short distances and a dependency of the efficiency on the walking direction with respect to the origami track. The former two observations were not expected and are non-trivial. Based on this study, we suggest three design modifications to improve future DNA walkers. Our study demonstrates the ability of the oxDNA model to resolve the dynamics of complex DNA machines, and its usefulness as an engineering tool for the design of DNA machines that operate in the three spatial dimensions.

Uncertainty quantification of a DNA origami mechanism using a coarse-grained model and kinematic variance analysis

Significant advances have been made towards the design, fabrication, and actuation of dynamic DNA nanorobots including the development of DNA origami mechanisms. These DNA origami mechanisms integrate relatively stiff links made of bundles of double-stranded DNA and relatively flexible joints made of single-stranded DNA to mimic the design of macroscopic machines and robots. Despite reproducing the complex configurations of macroscopic machines, these DNA origami mechanisms exhibit significant deviations from their intended motion behavior since nanoscale mechanisms are subject to significant thermal fluctuations that lead to variations in the geometry of the underlying DNA origami components. Understanding these fluctuations is critical to assess and improve the performance of DNA origami mechanisms and to enable precise nanoscale robotic functions. Here, we report a hybrid computational framework combining coarse-grained modeling with kinematic variance analysis to predict uncertainties in the motion pathway of a multi-component DNA origami mechanism. Coarse-grained modeling was used to evaluate the variation in geometry of individual components due to thermal fluctuations. This variation was incorporated in kinematic analyses to predict the motion pathway uncertainty of the entire mechanism, which agreed well with experimental characterization of motion. We further demonstrated the ability to predict the probability density of DNA origami mechanism conformations based on analysis of mechanical properties of individual joints. This integration of computational analysis, modeling tools, and experimental methods establish the foundation to predict and manage motion uncertainties of general DNA origami mechanisms to guide the design of DNA-based nanoscale machines and robots.

Force-Induced Unravelling of DNA Origami

The mechanical properties of DNA nanostructures are of widespread interest as applications that exploit their stability under constant or intermittent external forces become increasingly common. We explore the force response of DNA origami in comprehensive detail by combining AFM single molecule force spectroscopy experiments with simulations using oxDNA, a coarse-grained model of DNA at the nucleotide level, to study the unravelling of an iconic origami system: the Rothemund tile. We contrast the force-induced melting of the tile with simulations of an origami 10-helix bundle. Finally, we simulate a recently proposed origami biosensor, whose function takes advantage of origami behavior under tension. We observe characteristic stick-slip unfolding dynamics in our force-extension curves for both the Rothemund tile and the helix bundle and reasonable agreement with experimentally observed rupture forces for these systems. Our results highlight the effect of design on force response: we observe regular, modular unfolding for the Rothemund tile that contrasts with strain-softening of the 10-helix bundle which leads to catastropic failure under monotonically increasing force. Further, unravelling occurs straightforwardly from the scaffold ends inward for the Rothemund tile, while the helix bundle unfolds more nonlinearly. The detailed visualization of the yielding events provided by simulation allows preferred pathways through the complex unfolding free-energy landscape to be mapped, as a key factor in determining relative barrier heights is the extensional release per base pair broken. We shed light on two important questions: how stable DNA nanostructures are under external forces and what design principles can be applied to enhance stability.

Characterizing the Motion of Jointed DNA Nanostructures Using a Coarse-Grained Model

As detailed structural characterizations of large complex DNA nanostructures are hard to obtain experimentally, particularly if they have substantial flexibility, coarse-grained modeling can potentially provide an important complementary role. Such modeling can provide a detailed view of both the average structure and the structural fluctuations, as well as providing insight into how the nanostructure's design determines its structural properties. Here, we present a case study of jointed DNA nanostructures using the oxDNA model. In particular, we consider archetypal hinge and sliding joints, as well as more complex structures involving a number of such coupled joints. Our results highlight how the nature of the motion in these structures can sensitively depend on the precise details of the joints. Furthermore, the generally good agreement with experiments illustrates the power of this approach and suggests the use of such modeling to prescreen the properties of putative designs.

Irregular model DNA particles self-assemble into a regular structure

DNA nanoparticles with three-fold coordination have been observed to self-assemble in experiment into a network equivalent to the hexagonal (6.6.6) tiling, and a network equivalent to the 4.8.8 Archimedean tiling. Both networks are built from a single type of vertex. Here we use analytic theory and equilibrium and dynamic simulation to show that a model particle, whose rotational properties lie between those of the vertices of the 6.6.6 and 4.8.8 networks, can self-assemble into a network built from three types of vertex. Important in forming this network is the ability of the particle to rotate when bound, thereby allowing the formation of more than one type of binding motif. The network in question is equivalent to a false tiling, a periodic structure built from irregular polygons, and possesses 40 particles in its unit cell. The emergence of this complex structure, whose symmetry properties are not obviously related to those of its constituent particles, highlights the potential for creating new structures from simple variants of existing nanoparticles.

Structure and conformational dynamics of scaffolded DNA origami nanoparticles

Synthetic DNA is a highly programmable nanoscale material that can be designed to self-assemble into 3D structures that are fully determined by underlying Watson–Crick base pairing. The double crossover (DX) design motif has demonstrated versatility in synthesizing arbitrary DNA nanoparticles on the 5–100 nm scale for diverse applications in biotechnology. Prior computational investigations of these assemblies include all-atom and coarse-grained modeling, but modeling their conformational dynamics remains challenging due to their long relaxation times and associated computational cost. We apply all-atom molecular dynamics and coarse-grained finite element modeling to DX-based nanoparticles to elucidate their fine-scale and global conformational structure and dynamics. We use our coarse-grained model with a set of secondary structural motifs to predict the equilibrium solution structures of 45 DX-based DNA origami nanoparticles including a tetrahedron, octahedron, icosahedron, cuboctahedron and reinforced cube. Coarse-grained models are compared with 3D cryo-electron microscopy density maps for these five DNA nanoparticles and with all-atom molecular dynamics simulations for the tetrahedron and octahedron. Our results elucidate non-intuitive atomic-level structural details of DX-based DNA nanoparticles, and offer a general framework for efficient computational prediction of global and local structural and mechanical properties of DX-based assemblies that are inaccessible to all-atom based models alone.

The Need for Restructuring the Disordered Science of Amorphous Drug Formulations

The alarming numbers of poorly soluble discovery compounds have centered the efforts towards finding strategies to improve the solubility. One of the attractive approaches to enhance solubility is via amorphization despite the stability issue associated with it. Although the number of amorphous-based research reports has increased tremendously after year 2000, little is known on the current research practice in designing amorphous formulation and how it has changed after the concept of solid dispersion was first introduced decades ago. In this review we try to answer the following questions: What model compounds and excipients have been used in amorphous-based research? How were these two components selected and prepared? What methods have been used to assess the performance of amorphous formulation? What methodology have evolved and/or been standardized since amorphous-based formulation was first introduced and to what extent have we embraced on new methods? Is the extent of research mirrored in the number of marketed amorphous drug products? We have summarized the history and evolution of amorphous formulation and discuss the current status of amorphous formulation-related research practice. We also explore the potential uses of old experimental methods and how they can be used in tandem with computational tools in designing amorphous formulation more efficiently than the traditional trial-and-error approach.

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

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

Intrinsic Dynamics Analysis of a DNA Octahedron by Elastic Network Model

DNA is a fundamental component of living systems where it plays a crucial role at both functional and structural level. The programmable properties of DNA make it an interesting building block for the construction of nanostructures. However, molecular mechanisms for the arrangement of these well-defined DNA assemblies are not fully understood. In this paper, the intrinsic dynamics of a DNA octahedron has been investigated by using two types of Elastic Network Models (ENMs). The application of ENMs to DNA nanocages include the analysis of the intrinsic flexibilities of DNA double-helices and hinge sites through the calculation of the square fluctuations, as well as the intrinsic collective dynamics in terms of cross-collective map calculation coupled with global motions analysis. The dynamics profiles derived from ENMs have then been evaluated and compared with previous classical molecular dynamics simulation trajectories. The results presented here revealed that ENMs can provide useful insights into the intrinsic dynamics of large DNA nanocages and represent a useful tool in the field of structural DNA nanotechnology.