Abstract
Brownian particles confined in a system with periodic and asymmetric potential can be transported in a specific direction along the potential by repetitively switching the potential on and off. Here, we propose a DNA-based Brownian ratchet for directional transport of positively charged nanoparticles in which nanoparticle delivery follows the path dictated by a single, long, double-stranded DNA. We performed Brownian dynamics simulations to prove its realization using coarse-grained models. A periodic and asymmetric potential for nanoparticle binding is constructed along a single, long, double-stranded DNA molecule by a novel strategy that uses variation in sequence-dependent DNA flexibility. Directional and processive motion of nanoparticles is achieved by changing salt concentration repetitively over several cycles to switch the asymmetric potential on and off. This work suggests that double-stranded DNA molecules with elaborately designed flexibility variation can be used as a molecule-scale guide for spatial and dynamic control of nanoparticles for future applications.
INTRODUCTION
The migration of particles in solution is controlled in a specific direction, conventionally, by applying external fields including magnetic, electric, thermal, and concentration fields that bias the particle motion. In recent years, more sophisticated strategies based on a Brownian ratchet process have been developed for directional transport and sorting of Brownian particles in the absence of any external bias (1–14). In a type of Brownian ratchet mechanism called on-off ratchet, for example, Brownian particles confined in a system with a periodic and asymmetric potential are transported in a specific direction along the potential by repetitively switching the potential on and off (1–7, 13). In earlier Brownian ratchets, asymmetry in potential was obtained typically by using optical (6, 10, 13) or dielectrophoretic (5, 8, 11) forces. However, they are insufficient to manipulate particles of nanometer scale. Therefore, their use was limited to the transport of micrometer-sized particles, with only a few exceptions (7, 14). In addition, particle motion in other directions perpendicular to the potential is limited by the scale of particle confinement, typically micrometers, making it difficult to have precise control of particle localization in a three-dimensional space.
Here, we propose a DNA-based on-off Brownian ratchet for directional transport of positively charged nanoparticles (NPs) and prove the realization by performing Brownian dynamics (BD) simulations of coarse-grained models of a single, long, double-stranded DNA (dsDNA) molecule and a positively charged NP of 4 nm in diameter. An asymmetric potential for DNA-NP binding is constructed along a dsDNA fragment containing 384 base pairs corresponding to 130 nm in length such that the local flexibility of the base pair sequence is increased gradually along the fragment at physiological salt concentration as shown in Fig. 1 (A and B). The fragment is repeated so that the asymmetric potential is repeated periodically along a single, long, dsDNA molecule with several hundreds or thousands of base pairs, as shown in Fig. 2A. Directional and processive motion of an NP is achieved in a specific direction along the dsDNA by changing the salt concentration repetitively over several cycles (Fig. 2C) to switch the asymmetric potential on and off (Fig. 2, A and B). Since an NP does not dissociate from dsDNA at the concentrations under investigation, the conformation of a single dsDNA molecule dictates the path along which NPs are processively delivered, resulting in a precise control of NP localization from one end to the other along a dsDNA. The transport distance can reach several hundreds or thousands of nanometers along a dsDNA by synthesizing the dsDNA with repeated sequences of flexibility gradient and anchoring its two ends at two distant substrates (15).
(A) Model of dsDNA represented by a semiflexible chain of monomers, each of which models a group of 6 DNA base pairs. A dsDNA fragment with flexibility gradient is modeled by a chain of 64 monomers. The local persistence length (
(A and B) Long dsDNA containing four repeat fragments with the flexibility gradient, together with its DNA-NP binding potential energy along the entire dsDNA molecule at solution salt concentrations of Con = 0.15 M and Coff = 0.81M, respectively. (C) Solution salt concentrations are changed repeatedly between Con = 0.15 M and Coff = 0.81 M to turn the asymmetric DNA-NP binding potential on and off with a time period of τon and τoff.
The novelty of this DNA-based Brownian ratchet is in the use of variation in sequence-dependent dsDNA flexibility in addition to the construction of a periodic and asymmetric potential along a single dsDNA molecule. With increasing capability for precise DNA synthesis with specific base pair sequences (16), there are rapid developments in the field of DNA nanotechnology, including two- and three-dimensional structural assemblies of DNA molecules (17, 18), DNA-templated formation and control of nanoscale structures (15, 19), and construction of programmable DNA machines (20–22). In most of these developments, sequence-specific hybridization of complementary DNA strands has been used as an essential tool for connecting and exchanging different DNA strands. However, sequence-dependent variation of dsDNA flexibility has rarely been exploited. Although dsDNA is very rigid with an average persistence length of lp ~ 50 nm (a measure of rigidity or flexibility), it has been shown that local persistence length (
Particle delivery along a designated path proposed in this work is already familiar to biologists who study biological motors such as kinesin and myosin, the well-known examples of microscopic machinery walking along a cytoskeleton (26). Inspired by biological motors, artificial molecular machines that can move along a predefined track (27, 28) have recently been developed. It has been shown that an interlocking macrocycle can be transported directionally through a cyclic track (29), while DNA walkers (30–33) and small synthetic molecules (34, 35) can walk processively along their tracks, suggesting the possibility of achieving molecule-scale cargo transport along a designated path. However, it is nontrivial to achieve directional and processive motion in the design of the track and walker in these artificial machines. In addition, realization has only been verified along short molecular tracks. In contrast, the design of the flexibility-based DNA ratchet presented in this work is much simpler (by the synthesis of dsDNA containing the base pairs with varying degrees of
RESULTS
A dsDNA fragment with flexibility gradient
A dsDNA fragment with a flexibility gradient is modeled by a semiflexible polymer chain with 64 monomers and by mixing two types of monomers with different local persistence lengths of 40 and 51 nm, as shown in Fig. 1A. A single dsDNA monomer with a diameter of 2 nm and a charge of −12 coarsely models a group of 6 base pairs. These monomers with
Although the effect of flexibility variation on the dsDNA structure itself is not very notable, it becomes significant when negatively charged dsDNA fragment interacts with a highly charged cationic NP. At a physiological salt concentration of 0.15 M, an NP is wrapped around by 16 to 17 dsDNA monomers, as shown in the inset of Fig. 1B. This requires a large degree of dsDNA bending. Therefore, the free energy of DNA-NP binding is significantly influenced by dsDNA bending of sequences with different
Probability of the final location of an NP on a dsDNA fragment with flexibility gradient at Con = 0.15 M after BD simulations for a duration of τon = 104τBD, obtained from 200 independent simulations. The monomer number in the x axis represents the identity of the DNA monomer in the middle among those wrapping around an NP at a salt concentration of 0.15 M. Initial location of an NP was random on a dsDNA fragment at the beginning of the 200 simulations.
Unlike DNA-NP binding at 0.15 M (with 16 to 17 dsDNA monomers wrapping around an NP), electrostatic interaction between dsDNA and an NP is significantly weakened at an increased salt concentration of 0.81 M because of the screening effect of electrostatic interactions, and an NP is bound only with five to six dsDNA monomers. Less DNA binding leads to less dependence of dsDNA bending energy on flexibility, as shown in Fig. 1B, at a salt concentration of 0.81 M. In this case, the NP diffuses almost randomly in either direction on the dsDNA fragment because of negligible dependence of bending energy on the flexibility of dsDNA fragment.
Directional transport of an NP along a long dsDNA with repetitive sequences of flexibility gradient
At a salt concentration of 0.15 M, an NP bound to a dsDNA fragment with flexibility gradient slides directionally toward a more flexible region of the fragment subject to asymmetric bending energy. However, the distance of directional NP transport is limited by synthetic availability of the dsDNA fragment with a flexibility gradient in a desired length that can be prepared by combining DNA base pairs with different
Figure 2 (A and B) shows a long dsDNA with four repeat fragments of flexibility gradient at salt concentrations of 0.15 and 0.81 M. The binding energy between the dsDNA and an NP at a salt concentration of 0.15 M is asymmetric in each repeat fragment but periodic along the entire dsDNA, as shown in Fig. 2A. On the other hand, the variation of DNA-NP binding energy along the dsDNA becomes negligible at 0.81 M, as shown in Fig. 2B. Since asymmetry in DNA-NP binding energy is on and off at salt concentrations of 0.15 and 0.81 M, salt concentrations are hereafter referred to as Con and Coff, respectively. It is mentioned that the simulations are performed for a single repeat fragment. The presence of several repeat fragments is mimicked by using periodic boundary condition (PBC) connecting one end of the single fragment to the other end of its periodically imposed replica. This PBC is typically used in molecular simulations to mimic the macroscopic system by simulations of a nanoscale system.
Directional and processive motion of an NP bound to a long dsDNA is achieved by changing the salt concentration between Con and Coff, as shown in Fig. 2C. When an NP is bound to a specific repeat fragment at Con, the NP is wrapped around by 16 to 17 dsDNA monomers and directed to the minimum of asymmetric DNA-NP binding energy within the repeat fragment during BD simulations for a duration of τON. As shown in Fig. 3, the simulation duration of τON = 104τBD is enough for localized distribution of NP near more flexible regions within the repeat fragment with flexibility gradient. When the salt concentration is suddenly increased to Coff, however, the wrapping of the NP by dsDNA is released, and the NP is bound only by five to six dsDNA monomers. The NP then diffuses randomly in either direction on dsDNA, regardless of flexibility variation, during BD simulations for a duration of τoff. Since the starting location of the diffusing NP (determined by the minimum of the asymmetric binding energy at Con) is biased asymmetrically toward the more flexible side of the repeat fragment, the random diffusion of the NP after τoff at Coff results in a Gaussian distribution of the NP location spanning largely on the original repeat fragment and partially on the next repeat fragment close to the minimum of the original repeat fragment. Therefore, a non-negligible fraction of NPs initially bound to a more flexible region of the original repeat fragment can be found on the next repeat fragment after τoff. When the salt concentration is reduced again to Con, those that remained on the original repeat fragment during random diffusion are subjected to asymmetric binding energy along the original repeat fragment and directed again to the minimum of the binding energy. However, those that moved onto the next repeat fragment are subjected to the asymmetric binding energy of the new repeat fragment and directed to its minimum, resulting in a net forward movement by the length of the repeat fragment.
This is the mechanism of the on-off Brownian ratchet (1–4) that has been applied to several experimental systems (5–7, 13). The success of the on-off Brownian ratchet depends on the time interval τoff, during which the asymmetric potential is turned off and a Brownian particle diffuses randomly in either direction. If τoff is too short, most Brownian particles remain within the boundary of the original asymmetric potential. Then, they are located to its minimum when the potential is turned on. Therefore, the net movement to the forward direction is negligible. If τoff is too large, some Brownian particles can diffuse to the asymmetric potentials in the backward direction as well as to those in the forward direction. As a result, the net movement to the forward is canceled out by that to the backward. Appropriate choice of τoff can lead to efficient directional transport of the Brownian particles.
Therefore, we investigated the effect of different τoff on random diffusion of an NP at Coff, as shown in Fig. 4. For an NP initially located at the most flexible region of an original dsDNA repeat fragment (more specifically, we consider an NP on the 60th dsDNA monomer of an original 64-monomer repeat fragment), BD simulations were performed at Coff for durations between τoff = 5 × 102τBD and 104τBD. The probability data in Fig. 4 were obtained from 300 independent simulations at each τoff. Figure 4A shows the probability of the final location of an NP after each simulation duration of τoff = 5 × 102τBD, 103τBD, and 5 × 103τBD. In the figure, dsDNA monomers in the original repeat fragment are numbered from 1 to 64 with increasing flexibility (light orange–colored shades in the figure). The numbers of the dsDNA monomers in other repeat fragments in the forward direction increase sequentially from 65, whereas those in the backward direction decrease from 0 down to negative numbers. Figure 4B shows the probability of an NP in the original repeat fragment to jump to other repeat fragments. Separation of the final repeat fragment from the original repeat fragment is denoted by ΔNjump. Positive values of ΔNjump imply net NP displacement in the forward direction, whereas negative values imply net displacement in the backward direction.
(A) Probability of final location of an NP initially bound to the 60th monomer (most flexible region of the first repeat fragment) after each simulation duration of τoff = 5 × 102τBD, 103τBD, and 5 × 103τBD, calculated with a bin size of 16. (B) Probability of NP jump from the original repeat fragment to the other repeat fragments, ΔNjump, after τoff. ΔNjump is positive for forward direction and negative for backward direction. Light orange–colored shades indicate the original repeat fragment. These histogram data in the figures were obtained from 300 independent simulations at each τoff. (C) Probability of net directional motion of an NP in the forward direction (Pf − Pb) estimated by the difference between probabilities for forward (Pf) and backward (Pb) diffusion of an NP after τoff.
In Fig. 4, simulations performed for the duration of τoff = 5 × 102τBD resulted in nonzero probability data spanning only on the original repeat fragment (ΔNjump = 0) and the next repeat fragment in the forward direction (ΔNjump = 1). Probabilities of an NP for ΔNjump = 0 and 1 are 0.81 and 0.19, respectively, implying that when the asymmetric potential is turned on at Con, 19% of NPs move forward to the minimum of the next repeat fragment, while the other 81% of NPs remain on the original repeat fragment at the potential minimum. On the other hand, after a longer simulation duration of τoff ≥ 103τBD, not only there were more NPs diffusing in the forward direction (ΔNjump ≥ 1) but also a fraction of NPs were also found in the repeat fragment in the backward direction (ΔNjump ≤ −1). More specifically, at τoff = 5 × 103τBD, 32% of NPs were found on the repeat fragments in the forward direction (ΔNjump = 1, 2, and 3), and 20% of NPs were found in the backward direction (ΔNjump = −1 and −2). In such cases, the net NP displacement in the forward direction is canceled to some degree by that in the backward direction. Thus, the Brownian ratchet process is not very efficient. It is noted that the probability data in Fig. 4A were not very symmetric as expected for a Brownian particle performing random diffusion. This implies that the asymmetric potential in the DNA-NP binding energy is not completely off even at Coff. For simplicity, we will continue to say that the asymmetric potential is off at Coff, neglecting the weak asymmetry at this salt concentration.
Simple prediction of the efficiency of the Brownian ratchet process can be made on the basis of probabilities of NP displacement in forward and backward directions at Coff. Figure 4C shows the difference between forward and backward probabilities, Pf − Pb, at different τoff calculated on the basis of probability data shown in Fig. 4B. Pf is the sum of probabilities in the forward direction with ΔNjump ≥ 1, and Pb is the sum of those in the backward direction with ΔNjump ≤ −1. The larger Pf − Pb is, the more chance of NPs moving in the forward direction at Con. As shown in Fig. 4C, it was found that there was an optimal duration, which was at τoff = 103τBD, for efficient Brownian ratchet process with a maximum Pf − Pb of 0.22.
To prove the realization of this flexibility-based DNA Brownian ratchet process for directional NP transport, we performed BD simulations repetitively at Con and Coff, with a fixed duration of τon = 104τBD and different durations of τoff. A total of 40 cycles of on-off simulations were performed at Con and Coff for τoff = 5 × 102τBD and 103τBD, and 20 cycles of on-off simulations were performed for τoff = 5 × 103τBD and 104τBD. For each τoff, we performed 20 independent sets of these cyclic simulations to obtain the average net displacement.
Figure 5A shows the trajectories of 20 independent simulations for τoff = 103τBD. The x axis represents the simulation time over cycles of BD simulations at Con and Coff, while the y axis represents ΔNjump, which is the repeat fragment to which an NP belongs at each time relative to the starting repeat fragment. In the figure, flat regions are observed, with occasional transitions between different values of ΔNjump. Flat regions with constant ΔNjump were obtained during simulations at Con, while regions with sharp transitions were during simulations at Coff. Flat regions are often prolonged since the random diffusion at Coff does not lead to NP jump to another repeat fragment.
(A) Net displacement of an NP with time starting from an original repeat fragment and arriving at a repeat fragment separated by ΔNjump fragments, presented for all 20 independent simulation sets with 40 repeated changes between salt concentrations of Con and Coff with simulation durations of τon = 104τBD and τoff = 103τBD. (B) Net displacement averaged over the 20 independent simulation sets for different simulation durations of τoff = 5 × 102τBD, 103τBD, 5 × 103τBD, and 104τBD. In each simulation set, the asymmetric potential is turned on and off repeatedly 20 times for τoff = 5 × 103τBD and 104τBD and 40 times for τoff = 5 × 102τBD and 103τBD. Positive values for ΔNjump indicates net displacement in the forward direction.
Figure 5B shows the average of 20 independent simulations at each τoff. As predicted above on the basis of Pf − Pb, the net NP transport was the most efficient for τoff = 103τBD. After 40 cycles of simulations at Con and Coff, an NP is transported through three DNA repeat fragments away from the original repeat fragment. This shows that the flexibility-based DNA Brownian ratchet proposed in this work can be realized by repetitively changing salt concentrations between Con and Coff.
Stability of DNA-NP binding and NP diffusion under external force field
For successful realization of this proposed Brownian ratchet, periodic and repetitive change in the solution salt concentration is critical. One possibility lies in the use of a microfluidic device in which two solutions with different salt concentrations of Con and Coff are alternatively injected (38, 39). Although these alternating-injection microfluidic devices were developed to enhance the mixing efficiency in microfluidic channels, it is possible that the alternating injection of two salt solutions with different concentrations can be used for the periodic and repetitive change in the solution salt concentration required for the Brownian ratchet device proposed in this work. However, the flow of injected solutions in a microfluidic device exerts hydrodynamic drag forces on DNA and an NP, possibly influencing the stability of DNA-NP binding and diffusion of an NP on DNA. Therefore, we performed additional simulations under uniform force field that coarsely mimics hydrodynamic drag forces exerted by the laminar flow in a microfluidic channel. As a representative model system, the binding complex of an NP and a dsDNA with four repeat fragments (containing 256 dsDNA monomers) was simulated under constant, uniform force field in the subpiconewton range (see fig. S1). Two ends of the dsDNA are anchored at a separation of 230σ less than its fully stretched contour length, 256σ, to allow for DNA wrapping around an NP at the salt concentration of 0.15 M. Constant forces are applied to all DNA monomers in the direction perpendicular to the DNA. A twice stronger force is applied to an NP in the same direction because the hydrodynamic drag force increases with the size of an object.
Simulation results show that the wrapping structure of DNA around an NP at the salt concentration of Con = 0.15 M remains stable up to the force of 0.20 pN applied to an NP, as shown in Fig. 6A. At applied forces exceeding 0.20 pN, the DNA wrapping around an NP is disrupted (Fig. 6A), weakening the asymmetry in DNA-NP binding because of the relaxation of sharply bent DNA monomers from the wrapping structure. In the range of weak external forces (≤ 0.20 pN), an NP wrapped around by DNA at Con is kept localized at the minimum of the asymmetric potential in each DNA fragment with flexibility gradient during simulation duration of τON = 104τBD. At the salt concentration of Coff = 0.81 M, the binding of an NP on DNA remains stable at applied forces of 0.60 pN or less, whereas at stronger applied forces, an NP comes off from DNA. In the range of weak external forces (≤ 0.60 pN), an NP at Coff performs random Brownian motion on DNA because the asymmetric potential is off, as in the case of no external force. These results set the upper bound of external forces of 0.20 pN under which the proposed DNA-based Brownian ratchet can work. Under external forces exceeding 0.20 pN, the directional motion of an NP at Con along the flexibility gradient is ineffective because of the disruption of the wrapping structure of DNA-NP binding, while the random Brownian motion of an NP at Coff is less affected by the external force of 0.60 pN or less. Figure 6B presents the probability of the final location of an NP at the applied force of 0.20 pN after simulation duration of τoff = 103τBD at Coff, and the near symmetry in the probability distribution around a starting location, set to zero, proves the random Brownian motion of an NP at Coff (also in fig. S2). In summary, the flexibility-biased localization of an NP at Con and the random Brownian motion at Coff are observed even in the presence of external uniform forces of 0.20 pN or less, ensuring the success of the proposed Brownian ratchet device under weak external force.
(A) Snapshots of an NP with DNA at Con = 0.15 M subject to external forces, Fd, of 0.04, 0.20, and 0.40 pN. (B) Probability of final location of an NP at Coff = 0.81 M after BD simulations of τoff = 103τBD under external forces Fd = 0. 04 and 0.20 pN. An NP is initially bound to the most flexible region of a flexibility gradient, and this initial location is set to the monomer number of zero in the figures. The probability distributions in (B) are obtained from 180 independent simulations.
DISCUSSION
We proposed a DNA-based Brownian ratchet system for directional transport of a positively charged NP and proved its realization by BD simulations. Outstanding advantages of this system over earlier Brownian ratchet devices are that nanometer-scale Brownian particles can be transported in a specific direction by using strong electrostatic interactions between NPs and dsDNA and that the NP delivery follows the path along a single, long, dsDNA molecule, adding spatial precision for NP transport. Although we performed simulations only for a single NP transported along dsDNA, several NPs could be collected on the long dsDNA and processively transported toward the end of the dsDNA. Once NPs are collected at one end, they can be released by increasing salt concentration further above Coff. Such highly localized collection of NPs can be useful in many applications.
Although several molecule-scale walkers have been developed in recent years (27, 28, 30–35), the design of walker and track is nontrivial for processive and directional movements. On the contrary, the construction of the DNA Brownian ratchet proposed in this work is simpler. It only requires the design of a track with proper flexibility variation. Directional and processive transport of NPs can be achieved only by changing salt concentrations without needing alternative fuel solutions with different DNA strands. Last, the transport distance for NPs can be easily extended by duplicating the same DNA sequence with a flexibility gradient. Functionalization of NPs may extend the utility of the present DNA Brownian ratchet for transport of other functional groups attached onto NPs.
The sequence-dependent flexibility of a dsDNA has rarely been used as a tool for manipulating DNA-incorporated nanostructures. In our Brownian ratchet system, the asymmetric potential in each dsDNA fragment is built by aligning DNA base pairs with gradual flexibility variation. Although variation of dsDNA flexibility is not significant in structure determination of the dsDNA itself, the effect is augmented when it binds with highly charged cationic NPs, which requires a large degree of dsDNA bending. Therefore, one can find the utility of this flexibility-based NP control in the development of nanocomposites containing dsDNA and NPs.
Last, repetitive change in the solution salt concentration may induce hydrodynamic drag forces on DNA and NP by the flow during alternate injection of two solutions with different concentrations. It was shown from BD simulations under applied forces that the binding complex of DNA and NP remains stable, and the diffusion of an NP on DNA is not significantly influenced in the subpiconewton range of applied forces. Therefore, the realization of the proposed DNA-based Brownian ratchet relies on the development of a microfluidic device in which two solutions with different concentrations are injected alternately and the resultant hydrodynamic drag forces on nanometer-scale objects are in the subpiconewton range.
MATERIALS AND METHODS
Models of a dsDNA and an NP
As a model of a long dsDNA molecule, we used a semiflexible polymer chain. Each monomer had a diameter of 2 nm and a charge of −12, coarsely mimicking a group of 6 DNA base pairs. The charge per base pair was −2 at physiological pH, and a rise per base pair was 0.34 nm in B-DNA. An NP was modeled by a spherical particle with a diameter of 4 nm and a charge of +64. The size and charge of an NP were chosen such that the NP was wrapped around 1.6 to 1.7 turns by dsDNA at a physiological salt concentration of 0.15 M (25). In this case, conformational fluctuations of the DNA-NP binding complex enabled random sliding along dsDNA with uniform flexibility (25, 36).
The bonded interaction between a pair of DNA monomers was modeled by a combination of finite extension nonlinear elastic (FENE) potential energy
Nonbonded interactions were modeled by a combination of excluded volume interactions and electrostatic interactions. Excluded volume interactions between a pair of nonbonded DNA monomers and between an NP and DNA monomers were described by repulsive LJ potential energy identical to Eq. 2 above, but with a shift of r0 to describe size difference between an NP and DNA monomers
Here, r0 was 0 for a pair of DNA monomers, but 0.5σ for a pair of a DNA monomer and an NP because the size of an NP was twice that of DNA monomers. On the other hand, the linearized Debye-Hückel potential energy in Eq. 5 described electrostatic interactions between a pair of DNA monomers and between an NP and DNA monomers
Use of PBC for a dsDNA with repetitive flexibility gradient
In this work, the construction of a periodic and asymmetric potential along a long dsDNA molecule was critical for directional and processive transport of an NP bound to the dsDNA. For this purpose, a long dsDNA molecule was composed of several repeating fragments of 64 monomers (corresponding to 384 base pairs) as shown in Fig. 2A, each of which had a flexibility gradient. The flexibility gradient was constructed by varying the fraction of monomers with
Technically, however, only a single repeat fragment was explicitly considered in BD simulations with the use of PBC, in which one end of the repeat fragment was connected to the other end of its periodically imposed replica. Therefore, when an NP bound to one end of the repeat fragment moved to the other end across the simulation box, it was interpreted as a jump from the repeat fragment to the next nearest fragment represented as its periodic replica. The use of PBC to mimic a series of the repeat fragments with flexibility gradient is reasonable because the repeat fragment of 64 monomers had a contour length of 130 nm, much longer than the persistence length of our dsDNA model (between 40 and 51 nm in case of mixed monomers with different
BD simulation method
We performed BD simulations with a simulation package, GROMACS v. 4.5.4. BD simulations can simulate stochastic displacements of an NP and DNA monomers through a position Langevin equation, as shown below in Eq. 6, instead of explicitly incorporating solvent molecules into the system. Electrostatic screening by the presence of monovalent ions was implicitly taken into account by using Debye-Hückel potential energy in Eq. 5 instead of their explicit presence. Hydrodynamic interactions between DNA monomers and between a DNA monomer and an NP were not considered in this work. Hydrodynamic interactions are known to influence the diffusion of NPs quantitatively, particularly under strongly confined environment (40). However, the qualitative conclusion in this work that the Brownian ratchet device can be realized by repetitively combining the flexibility-biased NP diffusion at Con and the random, unbiased diffusion at Coff does not seem largely influenced by the quantitative change in NP diffusion.
Given a position of a particle at time t, the updated position at the next time step of t + Δt is determined by
Here,
Main results shown in Fig. 5 were obtained by cycles of BD simulations at Con = 0.15M and Coff = 0.81 M. Each simulation at Con = 0.15 M was performed for a fixed duration of 108 time steps, or τon = 104τBD. On the other hand, simulations at Coff = 0.81 M were performed for durations between 5 × 106(τoff = 5 × 102τBD) and 108 time steps (τoff = 104τBD), as mentioned in relevant parts of the main text. Statistical properties reported in this work were calculated by averaging several tens of independent simulations. This was mentioned in the main text when we discussed the data presented in Figs. 3 to 5.
The BD simulations under external forces were performed by adding an extra force to
SUPPLEMENTARY MATERIALS
Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/5/4/eaav4943/DC1
Fig. S1. Simulation model systems of DNA-NP binding at Con = 0.15 M under external force.
Fig. S2. Probability distribution of an NP at Coff = 0.81 M after Brownian motion for a duration of τoff = 103τBD under external force.
This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.
REFERENCES AND NOTES
- Copyright © 2019 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).