{'fraction_non_alphanumeric': 0.021953754221875812, 'fraction_numerical': 0.014289425824889582, 'mean_word_length': 4.8067131812181785, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 2, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Structure and electrical properties of DNA nanotubes embedded in lipid bilayer membranes

Abstract Engineering the synthetic nanopores through lipid bilayer membrane to access the interior of a cell is a long persisting challenge in biotechnology. Here, we demonstrate the stability and dynamics of a tile-based 6-helix DNA nanotube (DNT) embedded in POPC lipid bilayer using the analysis of 0.2 μs long equilibrium MD simulation trajectories. We observe that the head groups of the lipid molecules close to the lumen cooperatively tilt towards the hydrophilic sugar-phosphate backbone of DNA and form a toroidal structure around the patch of DNT protruding in the membrane. Further, we explore the effect of ionic concentrations to the in-solution structure and stability of the lipid-DNT complex. Transmembrane ionic current measurements for the constant electric field MD simulation provide the I-V characteristics of the water filled DNT lumen in lipid membrane. With increasing salt concentrations, the measured values of transmembrane ionic conductance of the porous DNT lumen vary from 4.3 to 20.6 nS. Simulations of the DNTs with ssDNA and dsDNA overhangs at the mouth of the pore show gating effect with remarkable difference in the transmembrane ionic conductivities for open and close state nanopores.


INTRODUCTION
Lipid bilayer membranes define the boundaries of a biological cell and regulate various kinds of cellular transport across it. The assembly of lipid molecules in aqueous medium is governed by the amphiphilic interactions where the polar head groups face outside and shield the hydrophobic tails from water to form a bilayer. Various peptides self-assemble to form pores which allow access of the interior of cell through these impermeable bilayers. These pores can be further engineered for various biophysical applications and synthetic biology (1)(2)(3)(4). Over past three decades, structural DNA nanotechnology has emerged as a convenient approach to construct nanostructures of arbitrary shape with sub-nanometer precision (5)(6)(7)(8). Several experimental groups have recently shown that self-assembled DNA nanostructures can mimic the naturally occurring nanopores in lipid bilayer membrane (9)(10)(11)(12). The diameter of these DNA Nanopores (DNPs) varies from 1.6 to 4 nm (13)(14)(15). Due to highly functionalized and automated chemical properties along with the conformational polymorphism, DNPs are easier to customize as compared to the conventional protein nanopores. In order to compensate the energy cost of pore formation in membrane, the charged sugar phosphate backbone of membrane spanning region of DNA is modified by covalently conjugating the hydrophobic lipid anchors like ethyl phosphorothioate (10), cholesterol (11,14,16), streptavidin (15), porphyrin (17) etc. DNPs in membrane have been characterized through advanced imaging techniques like confocal fluorescent microscopy, atomic force microscopy (AFM), aragose gel electrophoresis and single ion channel current recordings. These findings have injected fresh aspirations in field of synthetic nanopores. Recent experiments have shown that these nanopores are capable of the selective transport of molecular cargo across the membrane (15,16). The versatile self-assembled DNPs can be further useful for diverse biotechnological applications like targeted drug delivery, biosensing, genome sequencing and tools for biophysical studies etc. (18). Seifert et al. showed a voltage dependent gating of 6-helix DNP in lipid bilayer membranes associated with two states of conductance (19). With the advancements in atomistic and coarse-grained representations of DNA, MD simulations have become a handy tool to predict the properties of self-assembled DNA nanostructures (20)(21)(22)(23)(24). Yoo and Aksimentiev have studied several types of lipid spanning DNPs using atomistic and coarse grained MD simulations which have revealed many inherent aspects of these DNPs like ionic conductance, mechanical gating, electro-osmatic pressure effects etc (25,26). The assembly of DNA and lipid is very rare in nature due to their contrasting features and interactions. It is imperative to study the interaction between DNA and lipid membranes for their rational design (27). Aiming to understand the microscopic structure and the molecular interaction governing the selfassembly of DNPs in ambient conditions, herein we present an all atom molecular dynamics (MD) simulation study of DNT embedded in lipid bilayer membrane. We have created an atomistic model of a tile based 6-helix DNT (6HB) with an inner diameter of 2 nm, embedded in POPC lipid bilayer membrane. This DNT has been experimentally synthesized by Wang et al. (28) and we have investigated the nanoscale structure and elasticity of the same using the MD simulations (29). As it has been observed that the nature of counterions and salt concentration strongly affects the structure and stability of nucleic acids (30) as well as the behavior of lipid bilayer membranes (31), we also investigate the stability of these DNTs in lipid membrane as a function of monovalent (Na + and K + ) salt concentrations. Moreover, we measure the transmembrane ionic current and the Ohmic conductance of the channel from the MD simulations at constant electric field. Based on the analysis of multiple MD simulation trajectories, we propose a novel mechanism of the reorientation of lipid head groups to form a toroidal structure around DNT lumen in lipid bilayer membrane. This mechanism provides stability to the structure of DNPs in lipid bilayer membrane by reducing the energy barrier of insertion. We find that the structure of the DNT is most stable in 1M NaCl salt solution. We expect that our simulation study with atomistic resolution will provide important information about the physical properties of DNPs in lipid membrane and will guide future design of DNPs in this burgeoning field which is still in its infancy from various aspects.

Computational design of DNT in POPC lipid bilayer membrane
There are two widely known approaches to generate DNTs, namely tiles based self-assembly (32) and DNA origami technique. It has been known that due to the simplicity of the construction methods, higher yields and lower molecular weight, DNTs assembled from crossover tiles are easy to customize further to assemble as membranes nanopore (13). These DNTs are essentially the extension of DX (33) and TX (34) crossover molecules into a tubular geometry. We constructed an atomistic model of 57 base pair (bp) long DNT using a custom-built NAB programme. The NAB code operates with 'base first' strategy, then creating the sugar-phosphate backbone and connecting the phosphodiester bond according to the design. Six dsDNA, kept at the vertices of a hexagon of length 2 nm, were interconnected with adjacent dsDNA at two crossovers points on each side. These crossovers are Holliday like junctions, where the strands of DNA switch its double helical domain (35). The distance between two such crossovers is 7 bp or integer multiples of 7 bp. The design of the crossovers and sequence of nucleotides of DNT are taken from the experimental work by Wang et al. (28) Previously, we have extensively investigated the in situ structure, stability and elasticity of this DNT using molecular simulation (29,36). 1-Palmitoyl 2oleoyl-sn-glycero3-phosphocholine (POPC) is mono-cis unsaturated zwitterionic phospholipid with two asymmetric hydrocarbon chains, P and O, attached to phosphatialdylcholines (PCs). We generated the initial configuration of the POPC lipid bilayer membrane using CHARMM-GUI membrane builder (37). Next, we created a pore in the membrane by removing all the lipid molecules within a radius of 3 nm from the center of the membrane as shown in top panel of Supplementary Figure S1A in supporting information (SI). The DNT was inserted in the membrane pore by aligning the helical axis of DNT and bilayer normal in such a way that the center of mass of the membrane and DNT matches as shown in Figure 1A and bottom panel of Supplementary Figure S1A in SI. The method of generating the nanopore and insertion of DNT is similar to the recent study on cytolysin-A (cly-A) pore in POPC lipid bilayer membrane by Mandal et al. (38). We solvated the DNT embedded lipid bilayer membrane in a rectangular water box using TIP3P water model (39). Here, we ensured 20Å water shell from the outermost atom of DNA along the axis of DNT (z axis). Monovalent counterions (Na + /K + ) were added to neutralize the negative charge of the DNT. Further, we added additional Na + /K + and Clions to create the systems with four different molar concentrations i.e. 0.5 M NaCl, 1 M NaCl and 0.5 M KCl, 1 M KCl. The leap module creates a Coulombic grid of 1Å around the solute and places the ions at lowest electrostatic potential. During the process of ion placement, some overlapping water molecules were replaced by ions. Supplementary Figure  S1B in SI shows the side view of the initial built structure of DNT embedded in lipid membrane immersed into explicit salt solution. Supplementary Table S1 summarizes the atomistic details of all the five simulated systems.

Simulation methodology
All atom MD simulation have been performed using AM-BER14 simulation package (40). We used ff99 (41,42) force field parameters with parmbsc0 (43) and the recent ε dihedral refinements (44) to describe the bonded and nonbonded interactions for DNA. Lipid14 force field parameters were used to describe the interaction involving the POPC lipid molecules (45). We used Joung-Cheatham ion parameters optimized for TIP3P water model along with the Lorentz-Berthelot combination rule to describe their crucial non-bonded interaction in the system (46). These force field parameters are extensively validated and widely used in biomolecular simulations (47). The built structures were subjected to a series of steepest descent and conjugated gradient energy minimization steps to remove any bad contacts in the system while slowly removing the harmonic restraints from the solute. The systems were gradually heated upto 300 K in 105 ps MD simulation with 1 fs time step. While heating, the lipid and DNA atoms were harmonically restrained with 20 kcal/molÅ −2 force constant. Subsequently, we equilibrated the system for 5 ns at 1 atm pressure and 300 K temperature to acquire the correct density. Finally, 205 ns production MD simulation were carried out in NPT ensemble using 2 fs integration time step and periodic boundary condition. The translation motion of the center of mass of the system was removed after every 1000 time steps. We used Langevin thermostat with 2 ps −1 collision frequency and anisotropic pressure coupling with Berendsen barostat to maintain the constant temperature and pressure. We chose a new random number at the interval of 1 ns to avoid the synchronization artifacts in Langevin dynamics. Particle Mesh Ewald (PME) method integrated with AMBER simulation package was implemented to calculate the short range part of the Coulombic interaction as well as the long range electrostatic interactions between the various components of the system (48). We used 10Å cutoff to calculate the non-bonded interaction. All the bonds involving hydrogen were restrained using SHAKE algorithm (49). Supplementary video S1 shows the system trajectory (excluding water) for the 200 ns equilibrium MD simulation in the explicit 1 M NaCl electrolyte solution. To measure the ionic current as a function of transmembrane voltage difference, we performed MD simulation at constant electric field using NAMD code with AMBER force field parameters (50). We have extensively used CPPTRAJ (51) and VMD (52) programming suites for various analysis and graphics presented in this study. We have previously characterized the in silico stability of various DNA based nanostructures using the similar protocol (53)(54)(55).

Tilt angle of lipids chains: lipid DNT interaction
Hydrophilic nature of DNA backbone makes the insertion of the DNT into the hydrophobic core of lipid bilayer membrane energetically unfavorable. Experimentally, researchers have overcome this problem by modifying the DNA backbone with hydrophobic lipid based anchors to compensate the energy loss. Closely following the snapshots of the simulation, we see that for the lipid molecules close to the pore, the lipid head groups cooperatively tilt towards the DNT and shield the tail molecules by forming a toroidal shape around the water filled DNT channel within first few ns of the simulation as shown in schematically on the top panel of Figure 1. This conformation is maintained throughout the 205 ns long MD simulation. Figure 1 shows the instantaneous snapshots of the lipid head group atoms (Nitrogen and Phosphorus) at the beginning of the simulation, after 15 ns and at the end of 205 ns equilibrium MD simulation. Similar type of toroidal lipid structure was also  Here, the segment vector is the line joining n -1 to n + 1 carbon atoms in the alkyl chain. Similar scheme is routinely used to calculate the tilt angle of free standing lipid bilayer membrane in MD simulations. Further, we classified the lipid molecules into two groups, near and far, depending on the distance of the corresponding head group atoms from the center of the pore. We used a 3 nm distance cutoff from the center of the nanopore to define the near and far lipid as shown in Figure 2B where the near lipid have been highlighted in blue color. We calculate the tilt angle for all the lipid molecules and average them to get a mean tilt angle. Figure 2C and D shows the time evolution of average tilt angles for P and O chains of the POPC respectively. These plots describe that the near lipids distinctly tilts toward the central pore by large measure giving rise to higher tilt angle whereas the far lipids largely remains in the initial conformation similar to the free standing lipid bilayer membrane. Similar behavior was found in all the simulations at different molar concentrations. Figure 2E and F represent the normalized probability distribution of the tilt angles of far and near lipid molecules averaged over the last 10 ns of the 200 ns long MD simulation for P and O chains respectively. For both the chains, the most probable tilt angle for near lipids is significantly larger as compared to the far lipids. The values of average tilt angle from the last 10 ns equilibrium MD simulation have been summarized in the Table 1. We have encountered similar reorientation of lipid molecules in our previous electroporation study on POPC and DMPC lipid bilayers (56). This phenomenon has also been observed by Gopfrich et al. in their combined experimental and computational studies through the simulation snapshots of the local density of lipid head groups on ion channels made from porphyrin tagged single duplex DNA (25). The lipid reorientation mechanism is very much likely to play an important role in stability of the amphipathic DNT nanopore in lipid bilayer membrane. A similar mechanism was suggested by Khalid et al. for DNA translocation through the lipid bilayer membrane (57). To see any possible diffusion of DNT across the membrane, we performed additional 0.5 s long simulation for 0 M system. As discussed in appendix S2 of SI, the analysis suggests that the DNT structure is preserved in membrane over longer time scale as well.

Stability and conformational evolution
The structure and stability of biomolecules are strongly dependent on the nature of counterions and molarity of salt solution (58). We have compared the structure and conformational stability of DNT embedded in lipid membrane varying the molarity of the solution both for NaCl and KCl buffer solution. Figure 3 shows the instantaneous snapshots of system after 205 ns equilibrium MD simulations at various salt concentrations. These snapshots suggest that tubular structure of DNT is better maintained at higher molarity. We calculated the Root-Mean-Square-Deviation (RMSD) of DNT with respect to the initial energy minimized structure. Figure 4A compares the time evolution of the RMSD for DNT at different salt concentrations both for NaCl and KCl solution. The higher salt concentration gives additional stability to the system. We find that the RMSD for the DNT structure is least for the 1 M NaCl concentration followed by 1 M KCl, 0.5 M NaCl, 0.5 M KCl and 0 M respectively. We have also calculated the DNT radius profile at various salt concentrations and the radius profile is shown in Figure 4B. To calculate the radius profile, we divided the whole length of DNT into small segments of 1Å along the helical axis of DNT, next we find the center of each segment and calculate its distance from all the atoms in that segment. The average distance or radius of each segment from last 10 ns along with the standard deviation has been plotted against the segment position in Figure 4B.  (59). We find that the higher ionic concentration helps to keep the individual helical domains of DNT intact; hence they better maintain the tubular structure. At no salt (0 M case), system does not have enough ions to screen the electrostatic repulsion between helical domains and as a result the individual helices connected via the Holliday junction, opens up and bend towards the lipid bilayer membrane. The higher ions concentration also adds to the stability of DNT by providing more flexibility to the Holliday like crossover junctions of DNT.

Ion density, distribution and dynamics
It is known from the literature that depending on the size of cation, they possess different binding affinities and different binding sites to DNA (60-62). The van der Waal radii of Na + , K + and Clions used in our simulation are 1.37, 1.70 and 2.51Å respectively. With a smaller radius, Na + ions penetrate deeper into the grooves of DNA, condensate the structure more as compared to K + ions. This process gives rise to the more stable DNT in lipid bilayer membrane in Na + . We observe that the equilibrium structure of DNA nanopores (DNPs) significantly varies for different salt concentrations. To explore the distribution of ions along the bilayer normal we have computed the average number density of cations and anions from the last 10 ns of 0.2 s long equilibrium MD simulation as shown in Supplementary Figure  S3A and S3B respectively of appendix S3 in SI. For all the simulated systems we see a dip in the ion density at the central region containing the impermeable membrane patch except the charge neutral case (0M system). For this system, we have fewer cations compared to other systems. Large fraction of these ions enters inside the DNT lumen to neutralize the charge of DNT backbone and give rise to a relatively flatter ion density profile. We have also calculated the diffusion coefficient of ions from the slope of mean square displacement (MSD) curves shown in Supplementary Figure S3C and S3D. Table 1 gives the values of the diffusion coefficient of various ions in different systems. The details of ion distribution and MSD calculations are given in appendix S3 of SI.

Fluidity of lipid bilayer membrane
To quantify the changes in the structure of membrane due to the insertion of DNT, we have computed the time evolution of the area per head group (APH) available for the lipid molecules in each leaflet as shown in Figure 4C. APH decreases initially (up to 50 ns) and then becomes flatter with small fluctuations around a constant value for rest of the simulation. Table 1 summarizes the values of APH av-eraged over last 10 ns of the simulation. It is observed that the system in NaCl solution shows significantly lesser APH as compared to KCl. We found that the NaCl salt solutions make the system more stable and compact. Due to the electrostatic interactions between zwitterionic lipid membrane and ions, the self-diffusion of the lipid molecules decreases with the ionic strength of the solution. Hence, 1M NaCl solution has lowest APH as compared to the all other simulated systems. Next, we examine the normalized electron Nucleic Acids Research, 2018, Vol. 46, No. 5 2239 density along the bilayer normal averaged over last 10 ns as shown in Figure 4D. This curve reflects various regions with different charge density in the system and shows the compact structure of the DNT embedded in the lipid membrane. The central region with two peaks and one cusp at the middle corresponds to the bilayer thickness with DNT embedded. The peaks are slightly higher in 1 M salt solution reflecting the compact bilayer. This region is followed by the sudden decrease in electron density which originates at the terminal of the lipid bilayer membrane. After this region, bare DNT expands in both directions. The thickness of the lipid membrane is grossly constant for the equilibrium structure in all salt concentration.

Ionic current and I-V characteristics
To calculate the transmembrane ionic current, we performed MD simulations at five different values of transmembrane voltage: 10, 20, 50, 100 and 200 mV. We took the equilibrated structure after 100 ns long MD simulation and apply the electric field in the direction of bilayer normal aligned along z axis. We performed the simulations with the electric field for 50 ns for all the cases (five different salt concentrations and each for five different values of transmembrane voltage). In order to maintain the structural integrity of the porin under the effect of electric field, the terminal atoms of DNT and heavy atoms of membrane were constrained to their equilibrium structure using a harmonic force constant of 0.1 kcal/molÅ −2 . This restrain also prevents escaping of the DNT from the bilayer. The current density and hence the current of cylindrical nanopore can be estimated by Equation (1) as follows where N is total number of charge careers (ions) with drift velocity v d , z i (t + t) and z i (t) are the z coordinate of the ith ion at the interval of current sampling frequency (2 ps), L z is the thickness of the membrane which is taken to be 5 nm and q i is the charge of the ion. We observed the steady flow of ions across the membrane which sets up a constant current for a given transmembrane potential. We calculated the ionic current for cations and anions separately and also calculate the total current. It is found that cations contribute to most of the total ionic current. Because of the electrostatic repulsion of the negatively charged DNT, Clions do not prefer to enter inside the lumen giving rise to lesser contribution to current. We performed the running average of the total current values over 1 ns block. Supplementary Figure  S4A-E shows the ionic current as a function of the simulation time at five different molarities at various transmembrane voltage differences. Further, we averaged the ionic current values for last 10 ns of the simulation and plotted the average transmembrane current as a function of the voltage difference. Figure 5A shows the I-V characteristics of DNT lumen embedded in POPC lipid membrane at various salt concentrations. I-V curve shows that the channel is following the Ohm's law reasonably well. The solid lines show the linear fit for the I-V data points. The slope of the fit gives the Ohmic conductivity of the pore. The increasing molarity provides more charge careers, leading to the higher ionic conductance values. Table 1 summarizes the values of Ohmic conductance for the DNT lumen which vary from 4.31 to 20.64 nS depending on the salt concentration. Ionic conductance increases with increasing salt concentrations. The systems with KCl show larger conductance as compared to the system with NaCl due to the lower binding affinities K+ ions to DNA. This is also consistent with the ions diffusion analysis presented in previous section. It is important to note here that the ionic conductivities of the simulated systems are an order of magnitude higher compared to experimentally reported values from various groups as summarized in Table 2. Conductance of 20.64 nS is by far the highest conductance reported for 6-helix DNPs in experiments. In fact, previously published simulation studies by Aksimentiev et al. also found higher ionic conductance for simulated DNPs as compared to the experimentally measured values (14,26). The simulation results emphasizes that there is a scope of improvement in the experimental understanding of the ionic conductance through membrane protruding DNPs. We believe that the equilibrium orientation of the DNPs achieved through their selfassembly in experimental synthesis could be different than what is in general presumed and implemented in simulation studies. In simulation, we start with an ideal configuration of DNT embedded inside the hydrophobic core of the bilayer and perfectly aligned to bilayer normal. This could be different than what is synthesized by the self-assembly of nucleotides in the experiments. Further, a detail and careful study would be required to assess the fidelity of DNPs in lipid bilayer membranes and address the anomaly of lower ionic conductance measured in experiments. To get a theoretical estimate of this conductance we treat the DNT as a nanopore for a constant electrolyte transport. The conductance of such a cylindrical nanopore with channel length l, pore diameter d can be modeled using the following equation. where is the bulk conductance of electrolyte solution Here the second term in the denominator is the correction due to the hemispherical geometry of the nanopore at both the ends, which is also referred as the access resistance (63). This simple model is known to correctly fit a wide range of experimental data from synthetic membrane nanopores (64). Using the typical values of conductance for the bulk solution of 1M NaCl (8.5 S/m) and 1M KCl (10.86 S/m) at normal temperature (65) in Equation (2), we obtain a theoretical conductance of 1.27 and 1.62 nS respectively. In this theoretical framework, we approximated the DNTs with fix radius of 2 nm. But from the radius profile as shown in Figure  4B, we know that the DNTs swell during the course of simulation to attain the equilibrium structure. As the porous nature of DNT has been ignored in the theoretical estimation, this could be also a reason behind the lower conductance vales.  Figure S5A and B shows the snapshots of initial built and equilibrated structures of NPC and NPO after 100 ns MD simulation respectively. The analysis of the MD simulation trajectories shows that these transmembrane nanopores are stable in ambient conditions. In order to see the gating behavior due to the molecular design of NPC and NPO, we have computed the transmembrane ionic current at various transmembrane voltage differences as shown in Supplementary Figure S5A and B. From the figure we find that, for a given value of the transmembrane voltage difference, NPO always has higher ionic current as compared to NPC. From the average ionic current values shown in Supplementary Figure  S5A and B, we have plotted the I-V characteristics curve for NPO, NPC and NP in Figure 5B which shows Ohmic behavior of the simulated ionic currents. From the linear fitting of the I-V curve, we extract the ionic conductivity of these porous channels which comes out to be 27.53, 19.27 and 20.64 nS for NPO, NPC and NP respectively. In case of the NPC, we have single stranded loops over the native 6-helix DNT. These flexible single strands (each having a length of 14 bp) come close to each other during the MD equilibration to adopt a stable in-solution structure due to solvent effects. In the process, they block the passage of ions and act as a lock to central ionic channel. On the other hand, in case of NPO, by the virtue of its design, these single stranded oligomers are hybridized with complimentary strands. Hence the oligomers stay in their helical domain and allow the uninterrupted passage of electrolyte solution in transmembrane channel resulting higher current. The top view of equilibrium snapshots of the simulated structures shows that the central channel is partially blocked in NPC due to the conformational fluctuation of the overhangs.

CONCLUSIONS
In this study, we have demonstrated the stability and dynamics of a tile-based 6-helix DNT embedded in POPC lipid bilayer using the analysis of 0.2 s long equilibrium MD simulation trajectories. The lipid molecules in vicinity of the DNT reorient themselves to form a toroidal structure around it. We propose that this mechanism reduces the free energy barrier for the formation of a porous lumen Nucleic Acids Research, 2018, Vol. 46, No. 5 2241 of DNT in lipid bilayer membrane. The spontaneous formation of a toroidal pore cost energy penalty to the system. In our previous work we have shown that the activation barrier for a toroidal pore formation in POPC bilayer is ca. 25-27 kJ/mol under an electric field strength of 0.3-0.8 V/nm (56). At zero biasing potential, this could be substantially larger. As reported by Gopfrich et al., a toroidal pore of 2 nm radius costs an energy penalty of ca. 100-200 kJ/mol without any cholesterol anchor which reduces to -1000 kJ/mol for a pore of 6nm radius with 20 cholesterol anchors (14). Hence, to embed DNT porin in the bilayer, the energy barrier can be reduced by the addition of cholesterol anchors to DNA. Currently, there is no study on the microscopic understanding of the mechanism of DNT insertion in the bilayer. All the existing simulation studies report structure and conductivities of DNA porin which is preinserted in the bilayer. A thorough study is required to understand the mechanism of insertion and the self-assembly process of DNT in bilayer. We also show that higher salt concentration makes the DNT more compact and helps to stabilize them better in the lipid environment. Further, we observe the higher stability of the DNT in presence of Na + as compared to K + counterions due to their better condensation. The ionic current measurement under constant electric field allows us to calculate the I-V characteristics of these pores. The system displays nearly Ohmic behavior over a wide range of transmembrane voltage and ionic concentration. The conductance of the DNT varies from ca. 4 nS to ca. 21 nS depending on the ionic concentration. A conductance of ca. 21 nS is by far the highest conductance reported in literature for a DNP with radius of 2 nm. The MD simulations with open and close pore models of DNT capture the molecular gating of these ion channels. The results from this all atom MD simulation study will help to understand the phenomenology of DNPs in lipid bilayer membranes which will lead to their rational design for evolutionary biophysical applications.