Multi-resolution simulation of DNA transport through large synthetic nanostructures

Ionic flow through partially blocked nanopores

Employing atomistic molecular dynamics simulations, we investigate the ionic conductivity mechanisms in a partially blocked nanopore containing a centrally positioned spherical constriction, exploring the effects of pore diameter, surface charge, and blockage size. Our results show that ionic mobilities are significantly influenced by the polarity of the surface charge and the size of the pore gap. Particularly, we observe ion-specific effects for K+ and Cl- ions based on their size and charge, especially in sub-nanometer pore gaps. Furthermore, we find that the current flow in partially blocked nanopores sensitively depends on the surface charges, consistent with the calculated free energy profiles. The percentage of the current drop is found to be correlated to the volume of the spherical constriction with the effects more pronounced when the sizes of the spherical blockage and nanopore are comparable.

Nucleosomal DNA unwinding pathway through canonical and non-canonical histone disassembly

The nucleosome including H2A.B, a mammalian-specific H2A variant, plays pivotal roles in spermatogenesis, embryogenesis, and oncogenesis, indicating unique involvement in transcriptional regulation distinct from canonical H2A nucleosomes. Despite its significance, the exact regulatory mechanism remains elusive. This study utilized solid-state nanopores to investigate DNA unwinding dynamics, applying local force between DNA and histones. Comparative analysis of canonical H2A and H2A.B nucleosomes demonstrated that the H2A.B variant required a lower voltage for complete DNA unwinding. Furthermore, synchronization analysis and molecular dynamics simulations indicate that the H2A.B variant rapidly unwinds DNA, causing the H2A-H2B dimer to dissociate from DNA immediately upon disassembly of the histone octamer. In contrast, canonical H2A nucleosomes unwind DNA at a slower rate, suggesting that the H2A-H2B dimer undergoes a state of stacking at the pore. These findings suggest that nucleosomal DNA in the H2A.B nucleosomes undergoes a DNA unwinding process involving histone octamer disassembly distinct from that of canonical H2A nucleosomes, enabling smoother unwinding. The integrated approach of MD simulations and nanopore measurements is expected to evolve into a versatile tool for studying molecular interactions, not only within nucleosomes but also through the forced dissociation of DNA-protein complexes.

Exploring ssDNA translocation through α-hemolysin using coarse-grained steered molecular dynamics

Protein nanopores have proven to be effective for single-molecule studies, particularly for single-stranded DNA (ssDNA) translocation. Previous experiments demonstrated their ability to distinguish differences in purine and pyrimidine bases and in the orientation of the ssDNA molecule inside nanopores. Unfortunately, the microscopic details of ssDNA translocation over experimental time scales, which are not accessible through all-atom molecular dynamics (MD), have yet to be examined. However, coarse-grained (CG) MD simulations enable systems to be simulated over longer characteristic times closer to experiments than all-atom MD. This paper studies ssDNA translocation through α-hemolysin nanopores exploiting steered MD using the MARTINI CG force field. The impacts of the sequence length, orientation inside the nanopore and DNA charges on translocation dynamics as well as the conformational dynamics of ssDNA during the translocation are explored. Our results highlight the efficacy of CG molecular dynamics in capturing the experimental properties of ssDNA translocation, including a wide distribution in translocation times per base. In particular, the phosphate charges of the DNA molecule are crucial in the translocation dynamics and impact the translocation rate. Additionally, the influence of the ssDNA molecule orientation on the translocation rate is explained by the conformational differences of ssDNA inside the nanopore during its translocation. Our study emphasizes the significance of obtaining sufficient statistics via CG MD, which can elucidate the great variety of translocation processes.

A lumen-tunable triangular DNA nanopore for molecular sensing and cross-membrane transport

Synthetic membrane nanopores made of DNA are promising systems to sense and control molecular transport in biosensing, sequencing, and synthetic cells. Lumen-tunable nanopore like the natural ion channels and systematically increasing the lumen size have become long-standing desires in developing nanopores. Here, we design a triangular DNA nanopore with a large tunable lumen. It allows in-situ transition from expanded state to contracted state without changing its stable triangular shape, and vice versa, in which specific DNA bindings as stimuli mechanically pinch and release the three corners of the triangular frame. Transmission electron microscopy images and molecular dynamics simulations illustrate the stable architectures and the high shape retention. Single-channel current recordings and fluorescence influx studies demonstrate the low-noise repeatable readouts and the controllable cross-membrane macromolecular transport. We envision that the proposed DNA nanopores could offer powerful tools in molecular sensing, drug delivery, and the creation of synthetic cells.

Nanopore Translocation Reveals Electrophoretic Force on Noncanonical RNA:DNA Double Helix

Electrophoretic transport plays a pivotal role in advancing sensing technologies. So far, systematic studies have focused on the translocation of canonical B-form or A-form nucleic acids, while direct RNA analysis is emerging as the new frontier for nanopore sensing and sequencing. Here, we compare the less-explored dynamics of noncanonical RNA:DNA hybrids in electrophoretic transport to the well-researched transport of B-form DNA. Using DNA/RNA nanotechnology and solid-state nanopores, the translocation of RNA:DNA (RD) and DNA:DNA (DD) duplexes was examined. Notably, RD duplexes were found to translocate through nanopores faster than DD duplexes, despite containing the same number of base pairs. Our experiments reveal that RD duplexes present a noncanonical helix, with distinct transport properties from B-form DD molecules. We find that RD and DD molecules, with the same contour length, move with comparable velocity through nanopores. We examined the physical characteristics of both duplex forms using atomic force microscopy, atomistic molecular dynamics simulations, agarose gel electrophoresis, and dynamic light scattering measurements. With the help of coarse-grained and molecular dynamics simulations, we find the effective force per unit length applied by the electric field to a fragment of RD or DD duplex in nanopores with various geometries or shapes to be approximately the same. Our results shed light on the significance of helical form in nucleic acid translocation, with implications for RNA sensing, sequencing, and the molecular understanding of electrophoretic transport.

Nanopore translocation reveals electrophoretic force on non-canonical RNA:DNA double helix

Electrophoretic transport plays a pivotal role in advancing sensing technologies. So far, systematic studies have focused on translocation of canonical B-form or A-form nucleic acids, while direct RNA analysis is emerging as the new frontier for nanopore sensing and sequencing. Here, we compare the less-explored dynamics of non-canonical RNA:DNA hybrids in electrophoretic transport with the well-researched transport of B-form DNA. Using DNA/RNA nanotechnology and solid-state nanopores, the translocation of RNA:DNA (RD) and DNA:DNA (DD) duplexes was examined. Notably, RD duplexes were found to translocate through nanopores faster than DD duplexes, despite containing the same number of base pairs. Our experiments reveal that RD duplexes present a non-canonical helix with distinct transport properties from B-form DD molecules. We find RD and DD molecules with the same contour length move with comparable velocity through nanopores. We examined the physical characteristics of both duplex forms using atomic force microscopy, atomistic molecular dynamics simulations, agarose gel electrophoresis, and dynamic light scattering measurements. With the help of coarse-grained and molecular dynamics simulations, we find the effective force per unit length applied by the electric field to a fragment of RD or DD duplex in nanopores with various geometries or shapes to be approximately the same within experimental errors. Our results shed light on the significance of helical form in nucleic acid translocation, with implications for RNA sensing, sequencing, and molecular understanding of electrophoretic transport.

Engineering Biological Nanopore Approaches toward Protein Sequencing

Biotechnological innovations have vastly improved the capacity to perform large-scale protein studies, while the methods we have for identifying and quantifying individual proteins are still inadequate to perform protein sequencing at the single-molecule level. Nanopore-inspired systems devoted to understanding how single molecules behave have been extensively developed for applications in genome sequencing. These nanopore systems are emerging as prominent tools for protein identification, detection, and analysis, suggesting realistic prospects for novel protein sequencing. This review summarizes recent advances in biological nanopore sensors toward protein sequencing, from the identification of individual amino acids to the controlled translocation of peptides and proteins, with attention focused on device and algorithm development and the delineation of molecular mechanisms with the aid of simulations. Specifically, the review aims to offer recommendations for the advancement of nanopore-based protein sequencing from an engineering perspective, highlighting the need for collaborative efforts across multiple disciplines. These efforts should include chemical conjugation, protein engineering, molecular simulation, machine-learning-assisted identification, and electronic device fabrication to enable practical implementation in real-world scenarios.

Fast prediction of antibiotic permeability through membrane channels using Brownian dynamics

The efficient permeation across the Gram-negative bacterial membrane is an important step in the overall process of antibacterial action of a molecule and the one that has posed a significant hurdle on the way toward approved antibiotics. Predicting the permeability for a large library of molecules and assessing the effect of different molecular transformations on permeation rates of a given molecule is critical to the development of effective antibiotics. We present a computational approach for obtaining estimates of molecular permeability through a porin channel in a matter of hours using a Brownian dynamics approach. The fast sampling using a temperature acceleration scheme enables the approximate estimation of permeability using the inhomogeneous solubility diffusion model. Although the method is a significant approximation to similar all-atom approaches tested previously, we show that the present approach predicts permeabilities that correlate fairly well with the respective experimental permeation rates from liposome swelling experiments and accumulation rates from antibiotic accumulation assays, and is significantly, i.e., about 14 times, faster compared with a previously reported approach. The possible applications of the scheme in high-throughput screening for fast permeators are discussed.

The development of nucleic acids force fields: From an unchallenged past to a competitive future

Molecular dynamics simulations have strongly matured as a method to study biomolecular processes. Their validity, however, is determined by the accuracy of the underlying force fields that describe the forces between all atoms. In this article, we review the development of nucleic acids force fields. We describe the early attempts in the 1990s and emphasize their strong influence on recent force fields. State-of-the-art force fields still use the same Lennard-Jones parameters derived 25 years ago in spite of the fact that these parameters were in general not fitted for nucleic acids. In addition, electrostatic parameters also are deprecated, which may explain some of the current force field deficiencies. We compare different force fields for various systems and discuss new tests of the recently developed Tumuc1 force field. The OL-force fields and Tumuc1 are arguably the best force fields to describe the DNA double helix. However, no force field is flawless. In particular, the description of sugar-puckering remains a problem for nucleic acids force fields. Future refinements are required, so we review methods for force field refinement and give an outlook to the future of force fields.

Super-Resolution Detection of DNA Nanostructures Using a Nanopore

High-resolution analysis of biomolecules has brought unprecedented insights into fundamental biological processes and dramatically advanced biosensing. Notwithstanding the ongoing resolution revolution in electron microscopy and optical imaging, only a few methods are presently available for high-resolution analysis of unlabeled single molecules in their native states. Here, label-free electrical sensing of structured single molecules with a spatial resolution down to single-digit nanometers is demonstrated. Using a narrow solid-state nanopore, the passage of a series of nanostructures attached to a freely translocating DNA molecule is detected, resolving individual nanostructures placed as close as 6 nm apart and with a surface-to-surface gap distance of only 2 nm. Such super-resolution ability is attributed to the nanostructure-induced enhancement of the electric field at the tip of the nanopore. This work demonstrates a general approach to improving the resolution of single-molecule nanopore sensing and presents a critical advance towards label-free, high-resolution DNA sequence mapping, and digital information storage independent of molecular motors.

Emerging Approaches to DNA Data Storage: Challenges and Prospects

With the total amount of worldwide data skyrocketing, the global data storage demand is predicted to grow to 1.75 × 10¹⁴ GB by 2025. Traditional storage methods have difficulties keeping pace given that current storage media have a maximum density of 10³ GB/mm³. As such, data production will far exceed the capacity of currently available storage methods. The costs of maintaining and transferring data, as well as the limited lifespans and significant data losses associated with current technologies also demand advanced solutions for information storage. Nature offers a powerful alternative through the storage of information that defines living organisms in unique orders of four bases (A, T, C, G) located in molecules called deoxyribonucleic acid (DNA). DNA molecules as information carriers have many advantages over traditional storage media. Their high storage density, potentially low maintenance cost, ease of synthesis, and chemical modification make them an ideal alternative for information storage. To this end, rapid progress has been made over the past decade by exploiting user-defined DNA materials to encode information. In this review, we discuss the most recent advances of DNA-based data storage with a major focus on the challenges that remain in this promising field, including the current intrinsic low speed in data writing and reading and the high cost per byte stored. Alternatively, data storage relying on DNA nanostructures (as opposed to DNA sequence) as well as on other combinations of nanomaterials and biomolecules are proposed with promising technological and economic advantages. In summarizing the advances that have been made and underlining the challenges that remain, we provide a roadmap for the ongoing research in this rapidly growing field, which will enable the development of technological solutions to the global demand for superior storage methodologies.

How capture affects polymer translocation in a solitary nanopore

DNA capture with high fidelity is an essential part of nanopore translocation. We report several important aspects of the capture process and subsequent translocation of a model DNA polymer through a solid-state nanopore in the presence of an extended electric field using the Brownian dynamics simulation that enables us to record statistics of the conformations at every stage of the translocation process. By releasing the equilibrated DNAs from different equipotentials, we observe that the capture time distribution depends on the initial starting point and follows a Poisson process. The field gradient elongates the DNA on its way toward the nanopore and favors a successful translocation even after multiple failed threading attempts. Even in the limit of an extremely narrow pore, a fully flexible chain has a finite probability of hairpin-loop capture, while this probability decreases for a stiffer chain and promotes single file translocation. Our in silico studies identify and differentiate characteristic distributions of the mean first passage time due to single file translocation from those due to translocation of different types of folds and provide direct evidence of the interpretation of the experimentally observed folds [M. Gershow and J. A. Golovchenko, Nat. Nanotechnol. 2, 775 (2007) and Mihovilovic et al., Phys. Rev. Lett. 110, 028102 (2013)] in a solitary nanopore. Finally, we show a new finding-that a charged tag attached at the 5' end of the DNA enhances both the multi-scan rate and the uni-directional translocation (5' → 3') probability that would benefit the genomic barcoding and sequencing experiments.

Atomistic Simulation of Molecules Interacting with Biological Nanopores: From Current Understanding to Future Directions

Biological nanopores have been at the focus of numerous studies due to their role in many biological processes as well as their (prospective) technological applications. Among many other topics, recent studies on nanopores have addressed two key areas: antibiotic permeation through bacterial channels and sensing of analytes. Although the two areas are quite far apart in terms of their objectives, in both cases atomistic simulations attempt to understand the solute dynamics and the solute-protein interactions within the channel lumen. While decades of studies on various channels have culminated in an improved understanding of the key molecular factors and led to practical applications in some cases, successful utilization is limited. In this Perspective we summarize recent progress in understanding key issues in molecular simulations of antibiotic translocation and in the development of nanopore sensors. Moreover, we comment on possible advancements in computational algorithms that can potentially resolve some of the issues.

Discrimination of RNA fiber structures using solid-state nanopores

RNA fibers are a class of biomaterials that can be assembled using HIV-like kissing loop interactions. Because of the programmability of molecular design and low immunorecognition, these structures present an interesting opportunity to solve problems in nanobiotechnology and synthetic biology. However, the experimental tools to fully characterize and discriminate among different fiber structures in solution are limited. Herein, we utilize solid-state nanopore experiments and Brownian dynamics simulations to characterize and distinguish several RNA fiber structures that differ in their degrees of branching. We found that, regardless of the electrolyte type and concentration, fiber structures that have more branches produce longer and deeper ionic current blockades in comparison to the unbranched fibers. Experiments carried out at temperatures ranging from 20-60 °C revealed almost identical distributions of current blockade amplitudes, suggesting that the kissing loop interactions in fibers are resistant to heating within this range.