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

Influence of temperature on bend, twist and twist-bend coupling of dsDNA

The temperature-dependent bend and twist elasticities of dsDNA, as well as their couplings, were explored through all-atom molecular dynamics simulations. Three rotational parameters, tilt, roll, and twist, were employed to assess the bend and twist elasticities through their stiffness matrix. Our analysis indicates that the bend and twist stiffnesses decrease as the temperature rises, primarily owing to entropic influences stemming from thermodynamic fluctuations. Furthermore, the couplings between these rotational parameters also exhibit a decline with increasing temperature, although the roll-twist coupling displays greater strength than the tilt-roll and tilt-twist couplings, attributed to its more robust correction component. We elucidated the influence of temperature on bend and twist elasticities based on the comparisons between various models and existing data.

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.

The Response of the Replication Apparatus to Leading Template Strand Blocks

Duplication of the genome requires the replication apparatus to overcome a variety of impediments, including covalent DNA adducts, the most challenging of which is on the leading template strand. Replisomes consist of two functional units, a helicase to unwind DNA and polymerases to synthesize it. The helicase is a multi-protein complex that encircles the leading template strand and makes the first contact with a leading strand adduct. The size of the channel in the helicase would appear to preclude transit by large adducts such as DNA: protein complexes (DPC). Here we discuss some of the extensively studied pathways that support replication restart after replisome encounters with leading template strand adducts. We also call attention to recent work that highlights the tolerance of the helicase for adducts ostensibly too large to pass through the central channel.

Permeability of antioxidants through a lipid bilayer model with coarse-grained simulations

Oxidative stress caused by pollution and lifestyle changes causes an excess of free radicals that react chemically with cell constituents leading to irreversible damage. There are molecules known as antioxidants that reduce the levels of free radicals. Some pigments of fruits and vegetables known as anthocyanins have antioxidant properties. Their interaction with the cell membrane becomes a crucial step in studying these substances. In this research, molecular dynamics simulations, particularly, coarse-grained molecular dynamics (CGMD) were used. This technique aims to replace functional groups with corresponding beads that represent their level of polarity and affinities to other chemical groups. Also, umbrella sampling was carried out to obtain free energy profiles that describe well the orientation and location of antioxidants in a membrane considering Trolox, Cyanidin, Gallic Acid, and Resveratrol molecules to study the structural effects they cause on it. It was concluded in this study that an antioxidant when crossing the membrane does not cause either damage to the structural properties or the loss of packing and stratification of phospholipids. it was also observed that the most reactive part of the molecules could easily approach area A prone to lipid oxidation, which can describe the antioxidant capacity of these molecules.Communicated by Ramaswamy H. Sarma.

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.

Mechanistic insight into the structure, thermodynamics and dynamics of equilibrium gels of multi-armed DNA nanostars

The unique sequence specificity rule of DNA makes it an ideal molecular building block for constructing periodic arrays and devices with nanoscale accuracy and precision. Here, we present the self-assembly of DNA nanostars having three, four and five arms into a gel phase using a simplistic coarse-grained bead-spring model developed by Z. Xing, C. Ness, D. Frenkel and E. Eiser (Macromolecules, 2019, 52, 504-512). Our simulations show that the DNA nanostars form a thermodynamically stable fully bonded gel phase from an unstructured liquid phase with the lowering of temperature. We characterize the phase transition by calculating several structural features such as the radial distribution function and structure factor. The thermodynamics of gelation is quantified by the potential energy and translational pair-entropy of the system. The phase transition from an arrested gel phase to an unstructured liquid phase has been modelled using a two-state theoretical model. We find that this transition is enthalpy driven, and loss of configuration and translational entropy is counterpoised by enthalpic interaction of the DNA sticky-ends, which gives rise to a gel phase at low temperature. The absolute rotational and translational entropy of the systems, measured using a two-phase thermodynamic model, also substantiates the gel transition. The slowing down of the dynamics upon approaching the transition temperature from a high temperature demonstrates the phase transition to a gel phase. A detailed numerical simulation study of the morphology, dynamics and thermodynamics of DNA gelation can provide guidance for future experiments, is easily extensible to other polymeric systems, and is expected to help in understanding the physics of self-assembly.

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.

Nucleic acid nanostructures for in vivo applications: The influence of morphology on biological fate

The development of programmable biomaterials for use in nanofabrication represents a major advance for the future of biomedicine and diagnostics. Recent advances in structural nanotechnology using nucleic acids have resulted in dramatic progress in our understanding of nucleic acid-based nanostructures (NANs) for use in biological applications. As the NANs become more architecturally and functionally diverse to accommodate introduction into living systems, there is a need to understand how critical design features can be controlled to impart desired performance in vivo. In this review, we survey the range of nucleic acid materials utilized as structural building blocks (DNA, RNA, and xenonucleic acids), the diversity of geometries for nanofabrication, and the strategies to functionalize these complexes. We include an assessment of the available and emerging characterization tools used to evaluate the physical, mechanical, physiochemical, and biological properties of NANs in vitro. Finally, the current understanding of the obstacles encountered along the in vivo journey is contextualized to demonstrate how morphological features of NANs influence their biological fates. We envision that this summary will aid researchers in the designing novel NAN morphologies, guide characterization efforts, and design of experiments and spark interdisciplinary collaborations to fuel advancements in programmable platforms for biological applications.

Nanoscale structures and mechanics of peptide nucleic acids

Peptide nucleic acids (PNAs) are charge-neutral polyamide oligomers having extremely favorable thermal stability and high affinity to cell membranes when coupled with cationic cell-penetrating peptides (CPPs), as well as the encouraging antisense and antigene activity in cell-free systems. The study of the mechanical properties of short PNA molecules is rare both in experiments and theoretical calculations. Here, we studied the microscopic structures and elastic properties; namely, persistence length, stretch modulus, twist-stretch coupling, and structural crookedness of double-stranded PNA (dsPNA) and their hybrid derivatives using all-atom MD simulation and compared them with those of double-stranded DNA (dsDNA) and double-stranded RNA (dsRNA). The stretch modulus of the dsPNA is found to be ∼160 pN, an order of magnitude lower than that of dsDNA and smaller than dsRNA, respectively. Similarly, the persistence length of dsPNA is found to be ∼35 nm, significantly smaller than those of dsDNA and dsRNA. The PNA-DNA and PNA-RNA hybrid duplexes have elastic properties lying between that of dsPNA and dsDNA/dsRNA. We argue that the neutral backbones of the PNA make it less stiff than dsDNA and dsRNA molecules. Measurement of structural crookedness and principal component analysis additionally support the bending flexibility of dsPNA. Detailed analysis of the helical-rise coupled to helical-twist indicates that the PNA-DNA hybrid over-winds like dsDNA, while PNA-PNA and PNA-RNA unwind like dsRNA upon stretching. Because of the highly flexible nature of PNA, it can bind other biomolecules by adopting a wide range of conformations and is believed to be crucial for future nanobiotechnology research studies.

Diameter Dependent Melting and Softening of dsDNA Under Cylindrical Confinement

Carbon nanotubes (CNTs) are considered promising candidates for biomolecular confinement, including DNA encapsulation for gene delivery. Threshold values of diameters have been reported for double-stranded DNA (dsDNA) encapsulation inside CNTs. We have performed all-atom molecular dynamics (MD) simulations of dsDNAs confined inside single-walled CNTs (SWCNTs) at the physiologically relevant temperature of 300 K. We found that the dsDNA can be confined without being denatured only when the diameter of the SWCNT exceeds a threshold value. Below this threshold diameter, the dsDNA gets denatured and melts even at the temperature of 300 K. Our simulations using SWCNTs with chirality indices (20,20) to (30,30) at 300 K found the critical diameter to be 3.25 nm (corresponding to (24,24) chirality). Analyses of the hydrogen bonds (H-bonds), Van der Walls (VdW) energy, and other inter-base interactions show drastic reduction in the number of H-bonds, VdW energy, and electrostatic energies between the bases of dsDNA when it is confined in narrower SWCNTs (up to diameter of 3.12 nm). On the other hand, the higher interaction energy between the dsDNA and the SWCNT surface in narrower SWCNTs assists in the melting of the dsDNA. Electrostatic mapping and hydration status analyses show that the dsDNA is not adequately hydrated and the counter ion distribution is not uniform below the critical diameter of the SWCNT. As properly hydrated counter ions provide stability to the dsDNA, we infer that the inappropriate hydration of counter ions and their non-uniform distribution around the dsDNA cause the melting of the dsDNA inside SWCNTs of diameter below the critical value of 3.25 nm. For confined dsDNAs that do not get denatured, we computed their elastic properties. The persistence length of dsDNA was found to increase by a factor of about two and the torsional stiffness by a factor of 1.5 for confinement inside SWCNTs of diameters up to 3.79 nm, the stretch modulus also following nearly the same trend. Interestingly, for higher diameters of SWCNT, 3.79 nm and above, the dsDNA becomes more flexible, demonstrating that the mechanical properties of the dsDNA under cylindrical confinement depend non-monotonically on the confinement diameter.

Characterizing and Harnessing the Mechanical Properties of Short Single-Stranded DNA in Structured Assemblies

Precise engineering of DNA structures is of growing interest to solve challenging problems in biomolecular applications and beyond. The introduction of single-stranded DNA (ssDNA) into the DNA structure can play a pivotal role in providing high controllability of critical structural features. Herein, we present a computational model of ssDNA with structural applications to harness its characteristics. The nonlinear properties of nucleotide gaps are systematically characterized to construct a structural model of the ssDNA across length scales with the incorporation of a finite element framework. The proposed method shows the programmability of structural bending, twisting, and persistence length by implementing the ssDNA in various DNA structures with experimental validation. Our results have significant implications for DNA nanotechnology in expanding the boundary of design and analysis of structural shape and stiffness.