A coarse-grained model for DNA origami

Models and simulations of structural DNA nanotechnology reveal fundamental principles of self-assembly

DNA is not only a centrally important molecule in biology: the specificity of bonding that allows it to be the primary information storage medium for life has also allowed it to become one of the most promising materials for designing intricate, self-assembling structures at the nanoscale. While the applications of these structures are both broad and highly promising, the self-assembly process itself has attracted interest not only for the practical applications of designing structures with more efficient assembly pathways, but also due to a desire to understand the principles underlying self-assembling systems more generally, of which DNA-based systems provide intriguing and unique examples. Here, we review the fundamental physical principles that underpin the self-assembly process in the field of DNA nanotechnology, with a specific focus on simulation and modelling and what we can learn from them. In particular, we compare and contrast DNA origami and bricks and briefly outline other approaches, with an overview of concepts such as cooperativity, nucleation and hysteresis; we also explain how nucleation barriers can be controlled and why they can be helpful in ensuring error-free assembly. While high-resolution models may be needed to obtain accurate system-specific properties, often very simple coarse-grained models are sufficient to extract the fundamentals of the underlying physics and can enable us to gain deep insight. By combining experimental and simulation approaches to understand the details of the self-assembly process, we can optimise its yields and fidelity, which may in turn facilitate its use in practical applications.

Coarse-grained modeling of DNA-RNA hybrids

We introduce oxNA, a new model for the simulation of DNA-RNA hybrids that is based on two previously developed coarse-grained models-oxDNA and oxRNA. The model naturally reproduces the physical properties of hybrid duplexes, including their structure, persistence length, and force-extension characteristics. By parameterizing the DNA-RNA hydrogen bonding interaction, we fit the model's thermodynamic properties to experimental data using both average-sequence and sequence-dependent parameters. To demonstrate the model's applicability, we provide three examples of its use-calculating the free energy profiles of hybrid strand displacement reactions, studying the resolution of a short R-loop, and simulating RNA-scaffolded wireframe origami.

Application and prospects of nucleic acid nanomaterials in tumor therapy

Cancer poses a great threat to human life, and current cancer treatments, such as radiotherapy, chemotherapy, and surgery, have significant side effects and limitations that hinder their application. Nucleic acid nanomaterials have specific spatial configurations and can be used as nanocarriers to deliver different therapeutic drugs, thereby enabling various biomedical applications, such as biosensors and cancer therapy. In recent decades, a variety of DNA nanostructures have been synthesized, and they have demonstrated remarkable potential in cancer therapy related applications, such as DNA origami structures, tetrahedral framework nucleic acids, and dynamic DNA nanostructures. Importantly, more attention is also being paid to RNA nanostructures, which play an important role in gene therapy. Therefore, this review introduces the developmental history of nucleic acid nanotechnology, summarizes the applications of DNA and RNA nanostructures for tumor treatment, and discusses the development opportunities for nucleic acid nanomaterials in the future.

Computational Modeling of DNA 3D Structures: From Dynamics and Mechanics to Folding

DNA carries the genetic information required for the synthesis of RNA and proteins and plays an important role in many processes of biological development. Understanding the three-dimensional (3D) structures and dynamics of DNA is crucial for understanding their biological functions and guiding the development of novel materials. In this review, we discuss the recent advancements in computer methods for studying DNA 3D structures. This includes molecular dynamics simulations to analyze DNA dynamics, flexibility, and ion binding. We also explore various coarse-grained models used for DNA structure prediction or folding, along with fragment assembly methods for constructing DNA 3D structures. Furthermore, we also discuss the advantages and disadvantages of these methods and highlight their differences.

Ab initio predictions for 3D structure and stability of single- and double-stranded DNAs in ion solutions

The three-dimensional (3D) structure and stability of DNA are essential to understand/control their biological functions and aid the development of novel materials. In this work, we present a coarse-grained (CG) model for DNA based on the RNA CG model proposed by us, to predict 3D structures and stability for both dsDNA and ssDNA from the sequence. Combined with a Monte Carlo simulated annealing algorithm and CG force fields involving the sequence-dependent base-pairing/stacking interactions and an implicit electrostatic potential, the present model successfully folds 20 dsDNAs (≤52nt) and 20 ssDNAs (≤74nt) into the corresponding native-like structures just from their sequences, with an overall mean RMSD of 3.4Å from the experimental structures. For DNAs with various lengths and sequences, the present model can make reliable predictions on stability, e.g., for 27 dsDNAs with/without bulge/internal loops and 24 ssDNAs including pseudoknot, the mean deviation of predicted melting temperatures from the corresponding experimental data is only ~2.0°C. Furthermore, the model also quantificationally predicts the effects of monovalent or divalent ions on the structure stability of ssDNAs/dsDNAs.

To determine 3D structures and quantify stability of single- (ss) and double-stranded (ds) DNAs is essential to unveil the mechanisms of their functions and to further guide the production and development of novel materials. Although many DNA models have been proposed to reproduce the basic structural, mechanical, or thermodynamic properties of dsDNAs based on the secondary structure information or preset constraints, there are very few models can be used to investigate the ssDNA folding or dsDNA assembly from the sequence. Furthermore, due to the polyanionic nature of DNAs, metal ions (e.g., Na⁺ and Mg²⁺) in solutions can play an essential role in DNA folding and dynamics. Nevertheless, ab initio predictions for DNA folding in ion solutions are still an unresolved problem. In this work, we developed a novel coarse-grained model to predict 3D structures and thermodynamic stabilities for both ssDNAs and dsDNAs in monovalent/divalent ion solutions from their sequences. As compared with the extensive experimental data and available existing models, we showed that the present model can successfully fold simple DNAs into their native-like structures, and can also accurately reproduce the effects of sequence and monovalent/divalent ions on structure stability for ssDNAs including pseudoknot and dsDNAs with/without bulge/internal loops.

Predicting genome organisation and function with mechanistic modelling

Fitting-free mechanistic models based on polymer simulations predict chromatin folding in 3D by focussing on the underlying biophysical mechanisms. This class of models has been increasingly used in conjunction with experiments to study the spatial organisation of eukaryotic chromosomes. Feedback from experiments to models leads to successive model refinement and has previously led to the discovery of new principles for genome organisation. Here, we review the basis of mechanistic polymer simulations, explain some of the more recent approaches and the contexts in which they have been useful to explain chromosome biology, and speculate on how they might be used in the future.

Disulfide Dimerization of Neuronal Calcium Sensor-1: Implications for Zinc and Redox Signaling

Neuronal calcium sensor-1 (NCS-1) is a four-EF-hand ubiquitous signaling protein modulating neuronal function and survival, which participates in neurodegeneration and carcinogenesis. NCS-1 recognizes specific sites on cellular membranes and regulates numerous targets, including G-protein coupled receptors and their kinases (GRKs). Here, with the use of cellular models and various biophysical and computational techniques, we demonstrate that NCS-1 is a redox-sensitive protein, which responds to oxidizing conditions by the formation of disulfide dimer (dNCS-1), involving its single, highly conservative cysteine C38. The dimer content is unaffected by the elevation of intracellular calcium levels but increases to 10–30% at high free zinc concentrations (characteristic of oxidative stress), which is accompanied by accumulation of the protein in punctual clusters in the perinuclear area. The formation of dNCS-1 represents a specific Zn²⁺-promoted process, requiring proper folding of the protein and occurring at redox potential values approaching apoptotic levels. The dimer binds Ca²⁺ only in one EF-hand per monomer, thereby representing a unique state, with decreased α-helicity and thermal stability, increased surface hydrophobicity, and markedly improved inhibitory activity against GRK1 due to 20-fold higher affinity towards the enzyme. Furthermore, dNCS-1 can coordinate zinc and, according to molecular modeling, has an asymmetrical structure and increased conformational flexibility of the subunits, which may underlie their enhanced target-binding properties. In HEK293 cells, dNCS-1 can be reduced by the thioredoxin system, otherwise accumulating as protein aggregates, which are degraded by the proteasome. Interestingly, NCS-1 silencing diminishes the susceptibility of Y79 cancer cells to oxidative stress-induced apoptosis, suggesting that NCS-1 may mediate redox-regulated pathways governing cell death/survival in response to oxidative conditions.

Integrated computer-aided engineering and design for DNA assemblies

Recently, DNA has been used to make nanodevices for a myriad of applications across fields including medicine, nanomanufacturing, synthetic biology, biosensing and biophysics. However, current DNA nanodevices rely primarily on geometric design, and it remains challenging to rationally design functional properties such as force-response or actuation behaviour. Here we report an iterative design pipeline for DNA assemblies that integrates computer-aided engineering based on coarse-grained molecular dynamics with a versatile computer-aided design approach that combines top-down automation with bottom-up control over geometry. This intuitive framework allows for rapid construction of large, multicomponent assemblies from three-dimensional models with finer control over the geometrical, mechanical and dynamical properties of the DNA structures in an automated manner. This approach expands the scope of structural complexity and enhances mechanical and dynamic design of DNA assemblies.

Coarse-Grained Simulations of DNA Reveal Angular Dependence of Sticky-End Binding

Annealing between sticky ends of DNA is an intermediate step in ligation. It can also be utilized to program specific binding sites for DNA tile and origami assembly. This reaction is generally understood as a bimolecular reaction dictated by the local concentration of the sticky ends. Its dependence on the relative orientation between the sticky ends, however, is less understood. Here we report on the interactions between DNA sticky ends using the coarse-grained oxDNA model; specifically, we consider how the orientational alignment of the double-stranded DNA (dsDNA) segments affects the time required for the sticky ends to bind, τ_b. We specify the orientation of the dsDNA segments with three parameters: θ, which measures the angle between the helical axes, and ϕ₁ and ϕ₂, which measure rotations of each strand around the helical axis. We find that the binding time depends strongly on both θ and ϕ₂: ∼20-fold change with θ and 10-fold change with ϕ₂. The binding time is the fastest when the helical axes of duplexes are pointing toward each other and the sticky ends protrude from the farthest two points. Our result is relevant for predicting hybridization efficiency of sticky ends that are rotationally restricted.

Rapid Computational Analysis of DNA Origami Assemblies at Near-Atomic Resolution

Structural DNA nanotechnology plays an ever-increasing role in advanced biomolecular applications. Here, we present a computational method to analyze structured DNA assemblies rapidly at near-atomic resolution. Both high computational efficiency and molecular-level accuracy are achieved by developing a multiscale analysis framework. The sequence-dependent relative geometry and mechanical properties of DNA motifs are characterized by the all-atom molecular dynamics simulation and incorporated into the structural finite element model successfully without significant loss of atomic information. The proposed method can predict the three-dimensional shape, equilibrium dynamic properties, and mechanical rigidities of monomeric to hierarchically assembled DNA structures at near-atomic resolution without adjusting any model parameters. The calculation takes less than only 15 min for most origami-scale DNA nanostructures consisting of 7000-8000 base-pairs. Hence, it is expected to be highly utilized in an iterative design-analysis-revision process for structured DNA assemblies.

Measuring Internal Forces in Single-Stranded DNA: Application to a DNA Force Clamp

We present a new method for calculating internal forces in DNA structures using coarse-grained models and demonstrate its utility with the oxDNA model. The instantaneous forces on individual nucleotides are explored and related to model potentials, and using our framework, internal forces are calculated for two simple DNA systems and for a recently published nanoscopic force clamp. Our results highlight some pitfalls associated with conventional methods for estimating internal forces, which are based on elastic polymer models, and emphasize the importance of carefully considering secondary structure and ionic conditions when modeling the elastic behavior of single-stranded DNA. Beyond its relevance to the DNA nanotechnological community, we expect our approach to be broadly applicable to calculations of internal force in a variety of structures-from DNA to protein-and across other coarse-grained simulation models.

MrDNA: a multi-resolution model for predicting the structure and dynamics of DNA systems

Although the field of structural DNA nanotechnology has been advancing with an astonishing pace, de novo design of complex 3D nanostructures and functional devices remains a laborious and time-consuming process. One reason for that is the need for multiple cycles of experimental characterization to elucidate the effect of design choices on the actual shape and function of the self-assembled objects. Here, we demonstrate a multi-resolution simulation framework, mrdna, that, in 30 min or less, can produce an atomistic-resolution structure of a self-assembled DNA nanosystem. We demonstrate fidelity of our mrdna framework through direct comparison of the simulation results with the results of cryo-electron microscopy (cryo-EM) reconstruction of multiple 3D DNA origami objects. Furthermore, we show that our approach can characterize an ensemble of conformations adopted by dynamic DNA nanostructures, the equilibrium structure and dynamics of DNA objects constructed using off-lattice self-assembly principles, i.e. wireframe DNA objects, and to study the properties of DNA objects under a variety of environmental conditions, such as applied electric field. Implemented as an open source Python package, our framework can be extended by the community and integrated with DNA design and molecular graphics tools.

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.

From quantum to subcellular scales: multi-scale simulation approaches and the SIRAH force field

Modern molecular and cellular biology profits from astonishing resolution structural methods, currently even reaching the whole cell level. This is encompassed by the development of computational methods providing a deep view into the structure and dynamics of molecular processes happening at very different scales in time and space. Linking such scales is of paramount importance when aiming at far-reaching biological questions. Computational methods at the interface between classical and coarse-grained resolutions are gaining momentum with several research groups dedicating important efforts to their development and tuning. An overview of such methods is addressed herein, with special emphasis on the SIRAH force field for coarse-grained and multi-scale simulations. Moreover, we provide proof of concept calculations on the implementation of a multi-scale simulation scheme including quantum calculations on a classical fine-grained/coarse-grained representation of double-stranded DNA. This opens the possibility to include the effect of large conformational fluctuations in chromatin segments on, for instance, the reactivity of particular base pairs within the same simulation framework.

Characterization of Hydroxyproline-Containing Hairpin-Like Antimicrobial Peptide EcAMP1-Hyp from Barnyard Grass (Echinochloa crusgalli L.) Seeds: Structural Identification and Comparative Analysis of Antifungal Activity

Herein, we describe a modified form of the antimicrobial hairpin-like peptide EcAMP1, isolated from barnyard grass (E. crusgalli) seeds, which is structurally characterized by a combination of high-pressure liquid chromatography, mass spectrometry, and automated Edman sequencing. This derivate has a single amino acid substitution (Pro19Hyp) in the second α-helical region of the molecule, which is critical for the formation of the hydrophobic core and the secondary structure elements. Comparing the antifungal activity of these two peptides, we found that the modified EcAMP1-Hyp had a significantly weaker activity towards the most-sensitive plant pathogenic fungus Fusarium solani. Molecular dynamics simulations and in vitro binding to the commercial polysaccharides allowed us to conclude that the Pro-19 residue is important for binding to carbohydrates located in the spore cell wall and it chiefly exhibits a fungistatic action representing the hyphal growth inhibition. These data are novel and significant for understanding a role of α-hairpinins in plant immunity.

QM/MM Description of Newly Selected Catalytic Bioscavengers Against Organophosphorus Compounds Revealed Reactivation Stimulus Mediated by Histidine Residue in the Acyl-Binding Loop

Butyrylcholinesterase (BChE) is considered as an efficient stoichiometric antidote against organophosphorus (OP) poisons. Recently we utilized combination of calculations and ultrahigh-throughput screening (uHTS) to select BChE variants capable of catalytic destruction of OP pesticide paraoxon. The purpose of this study was to elucidate the molecular mechanism underlying enzymatic hydrolysis of paraoxon by BChE variants using hybrid quantum mechanical/molecular mechanical (QM/MM) calculations. Detailed analysis of accomplished QM/MM runs revealed that histidine residues introduced into the acyl-binding loop are always located in close proximity with aspartate residue at position 70. Histidine residue acts as general base thus leading to attacking water molecule activation and subsequent SN2 inline hydrolysis resulting in BChE reactivation. This combination resembles canonical catalytic triad found in active centers of various proteases. Carboxyl group activates histidine residue by altering its pK_a, which in turn promotes the activation of water molecule in terms of its nucleophilicity. Observed re-protonation of catalytic serine residue at position 198 from histidine residue at position 438 recovers initial configuration of the enzyme’s active center, facilitating next catalytic cycle. We therefore suggest that utilization of uHTS platform in combination with deciphering of molecular mechanisms by QM/MM calculations may significantly improve our knowledge of enzyme function, propose new strategies for enzyme design and open new horizons in generation of catalytic bioscavengers against OP poisons.

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.