# Continuous Approximation for Interaction Energy Transfer of DNA through Lipid Bilayers

Authors: Mansoor H. Alshehri
Journal: Discrete Dynamics in Nature and Society
 (2018)

---

## Abstract

In this study the interaction energies for single-stranded DNA and double-stranded DNA molecules with a lipid bilayer are investigated. The 6-12 Lennard-Jones potential and continuous approximation are used to derive analytical expressions for these interaction energies. Assuming that there is a circular gap in the lipid bilayer, we determine the relationship of the molecular interaction energy, including the circular gap radius and the perpendicular distance of the single-stranded DNA and double-stranded DNA molecules from the gap. For both single-stranded and double-stranded DNA molecules, the relationship between the minimum energy location and the hole radiusb is calculated; in the case of the double-stranded DNA molecule, we assume that the helical phase angle ϕ is equal to π. By minimizing the total interaction energies, the results demonstrate that the single-stranded DNA and double-stranded DNA molecules move through a lipid bilayer when the gap radius b> 10 Å and b> 13.8 Å, respectively. The results present in this project can be leveraged to understand the interactions between cell-penetrating peptides and biomembranes, which may improve gene and drug delivery.


---

## Body

## 1. Introduction

Recently, nanomaterials have promised an unexpected growth in research and applications in many different areas, due to their nanoscales size, unique features, and distinct properties. Watson and Crick discovered the structure of DNA in 1953, one century after the discovery of the existence of DNA [1]. Owing to the self-assembling, unique physicochemical properties, and the geometric structures of DNA molecules, they have attracted attention as a promising material that has many potential applications in biotechnology and biomedicine such as stochastic biosensors, for the controlled translocation of proteins and as drug carriers across membranes [2–5]. Lipid/DNA complexes hold great promise for future medical applications such as transfect genetic material to the cell core and as novel treatment for various inherited diseases and cancers [6–9]. There are a number of research studies about the interaction behavior of the translocation of molecules into cells, which can help to explain the toxicity of nanoparticles. Maingi et al. use molecular dynamics (MD) simulations to investigate the interactions of a simple DNA nanopore with a lipid bilayer. They demonstrate the close packing of lipids around the stably inserted DNA pore and its cation selectivity [10]. He et al. determined the interaction mechanism between polyarginine peptides and asymmetric membranes by performing a coarse-grained molecular dynamics (CGMD) simulation. Their results show that peptides could be moved through a lipid bilayer by inducing a hydrophilic pore formation in the asymmetric membrane [11, 12]. Khalid et al. utilize (CGMD) simulation to investigate the transfer potential of the DNA through a DPPC/DMTAP bilayer. They also find a high energy barrier to DNA insertion into the bilayer hydrophobic core of the bilayer [13]. In addition, a number of simulations have provided good information about the interactions of DNA and lipids, such as [14–16], where the computational requirements of such a problem did not allow the exploration of the transmigration of the DNA inside the bilayer. Baowan et al. use continuous approximation together with the Lennard-Jones potential function to calculate the interaction energy of silica and carbon nanoparticles with a lipid bilayer [17–19]. Alshehri determined the interactions of boron nitride nanotubes as they moved through lipid bilayer membranes using the same technique [20]. Van der Waals energy has played an important role; hence both electrostatic interactions and van der Waals were considered as a part of the computation. In this study, classical applied mathematics is used to calculate the most important interactions of the system and provide the theoretical results which can pave the way toward further improving this area of study. Using continuous approximation and the Lennard-Jones potential function, the molecular interaction energies between ssDNA and dsDNA molecules and the dipalmitoylphosphatidylcholine bilayer (DPPC) are calculated. The penetration behavior of the DNA molecules through an assumed circular gap in the lipid bilayer is investigated. The model formulations for ssDNA and dsDNA molecules and the lipid bilayer are presented in the following section. The mathematical derivations for the ssDNA and dsDNA interaction with the lipid are detailed in Section 3. Moreover, in Section 4 the numerical results are presented, and finally in Section 5 a brief summary of the project is presented.


## 2. Modelling Approach

We determine the molecular interatomic energy between two molecules by employing continuous approximation and the Lennard-Jones potential. We assume that the 6-12 Lennard-Jones potential is expressed as(1)Pρ=-Aρ6+Bρ12,where ρ is the distance between two well-defined distinct molecules, A=4εσ6 and B=4εσ12 denote the attractive and the repulsive constants, respectively, and σ is the interatomic distance when the potential is zero and ε is the energy well depth. Using a continuous approach, atoms are assumed to be evenly distributed over the molecules and the average atomic densities provided are averaged over its surface, which means that the summation over all atoms involved is replaced by surface integrals. The magnitude of total energy arising from molecular interatomic energy is then obtained mathematically by performing a surface integral over the surface of each molecule, given by(2)E=η1η2∫S1∫S2PρdS1dS2,where η1 and η2 denote mean volume densities or the mean surface densities of the atoms on each molecule, and P(ρ) is the potential function for two unbonded atoms. For convenience, we could introduce the integral In as(3)In=∫S1∫S2ρ-2ndS1dS2,n=3,6.Thus, the total energy can given by(4)E=η1η2-AI3+BI6.At this point, we investigate the energy behavior of the molecular interaction energy of a lipid bilayer and a DNA molecule, by determining the penetration behavior of the ssDNA and dsDNA molecules through an assumed circular gap of radius b in the bilayer. The ssDNA molecule is assumed to be a single helicoid, with a ruled surface having a helix as its boundary. This presumption follows the experimental studies [25–27]. A typical point on the surface of an ssDNA molecule is specified by the coordinates:(5)Hθ,t=x,y,z=rtcosθ,rtsinθ,cθ2π,where c=34 Å and r=5.8 Å are the unit cell length and the radius of the ssDNA helix, respectively, and the two parametric variables θ and t are such that -π<θ<π and 0<t<1. In addition, the dsDNA molecule is modelled as a surface with a double helicoid geometry, and the parametric equation is given by(6)GΘ,t=x,y,z=R2cosΘ+cosΘ-ϕ+tcosΘ-cosΘ-ϕ,R2sinΘ+sinΘ-ϕ+tsinΘ-sinΘ-ϕ,cΘ2π,where ϕ is the helical phase angle parameter, R=10 Å is the radius of the dsDNA helix, c=34 Å is the unit cell length of the dsDNA helix, and the parametric variables Θ and t are such that -π<Θ<π and -1<t<1. Moreover, with reference to Figure 1 the bilayer of dipalmitoylphosphatidylcholine (DPPC) has been selected as a part of this experiment, each lipid molecule is composed of 12 particles, which are shown in the MARTINI force field as a head group consisting of choline (Q0) and phosphate (Qa) groups, an intermediate layer consisting of a glycerol group (Na), and a carbon tail group (C1) [23]. We note that the phase behavior of complexes and the considerable diversity in structure stability are depending on the choice of the uncharged lipid and that the chemical type of divalent ion can impact on the stability of the complexes which involves specific interactions [28]. From the work of Bawan et el. [17], the spacing between the two layers of lipids is assumed to be 3.36 Å, and the initial value for the area per lipid is assumed to be 65 Å2, which is taken from the work of Berger et al. [29]. The head and tail groups could be modelled as infinite boxes of thicknesses 4 and 15, respectively, following the work of Bawan et el. [18].Figure 1
Schematic of the lipid bilayer.Using the rectangular Cartesian coordinate system(x,y,z), the six layers of the lipid bilayer are assumed to have the coordinates (acosϑ,acosϑ,z) and (acosϑ,acosϑ,z-37.36) for the upper and lower head groups, respectively. By assuming that the surface of the upper head group is located on the xy-plane, z=0, the coordinates of a typical point on this plane can be given by (acosϑ,acosϑ,0), where b<a<∞ and 0<ϑ<2π. In addition, the upper and lower intermediate layers have coordinates (acosϑ,acosϑ,z-4) and (acosϑ,acosϑ,z-37.36), respectively, as well as the upper and lower tail groups (acosϑ,acosϑ,z-4) and (acosϑ,acosϑ,z-22.36), respectively.


## 3. Mathematical Derivation

In this study, the techniques of volume and surface integration are used to compute the total interaction energy of the ssDNA and dsDNA molecules with a lipid bilayer. The head group and the tail group are modelled as rectangular boxes, and the intermediate group is assumed to be a flat plane. Moreover, the bilayer is assumed to stay in an equilibrium structure while interacting with the DNA molecule; thus the interaction between the two layers could be ignored, and the total interaction energy of the system is assumed to contain(i)
The interactions between the volumes of the upper and lower head groups and the surfaces of the ssDNA and dsDNA molecules(ii)
The interactions between the surfaces of the upper and lower intermediate layers and the surfaces of the ssDNA and dsDNA molecules(iii)
The interactions between the volumes of the upper and lower tail groups and the surfaces of the ssDNA and dsDNA moleculesIn the following subsections, the mathematical derivations of the interactions of this system are detailed.


### 3.1. Interaction Energy of ssDNA Penetrating a Lipid Bilayer Hole

A unit cell of the ssDNA molecule of radiusr is examined at a distance Z above a gap of radius b>r in the lipid bilayer, where Z is the perpendicular distance from the gap of the uppermost surface of the lipid bilayer to the center of the ssDNA molecule, as shown in Figure 2. Thus, a typical surface element of the ssDNA molecule has the coordinates given by rtcosθ,rtsinθ,cθ/2π+Z, and without loss of generality the coordinates of a typical point on the ssDNA molecule can be given by (α,0,Z+cθ/2π), where α=rt. Firstly, the interaction energy between an ssDNA molecule and a flat plane, which has the coordinates (acosϑ,acosϑ,0), is calculated. Therefore, the distance of the surface of an ssDNA to the plane is given by(7)ρ2=α-acosϑ2+-acosϑ2+Z+cθ2π2=α2-2αcosϑ+a2+Z+cθ2π2=α-a2+4aαsin2ϑ2+Z+cθ2π2.and from (3) the integral In can be written as(8)Inp=∫-ππ∫-ππ∫01∫b∞aα-a2+4aαsin2ϑ2+Z+cθ2π2-n×1+4r2π2c2t21/2dadtdθdϑ,where p is referred to in the plane interaction. By setting ξ=α-a2+Z+cθ/2π2 and β=4aα and defining the integral Kn as(9)Kn=∫-ππξ+βsin2ϑ2-ndϑ,upon making the substitution u=sin2(ϑ/2)⇒dϑ=u-1/21-u-1/2du, the integral Kn becomes(10)Kn=2∫01ξ+βu-nu-1/21-u-1/2du=2ξ-n∫011+βξu-nu-1/21-u-1/2du=2πξ-nFn,12;1;-βξ,where F(a∗,b∗;c∗;z∗) is the standard hypergeometric function [30]. Moreover, by using a Pfaff transformation [30], Kn becomes(11)Kn=2πξ-nξξ+β1/2F1-n,12;1;βξ+β=2πξ-nξξ+β1/2∑k=0n-11-nk1/2kk!2βξ+βk=2πrt-a2+Z+cθ2π2-n+1/2∑k=0n-14artk1-nk1/2kk!2rt+a2+Z+cθ/2π2k.Therefore, there are three remaining integrals for Inp:(12)Inp=2π∑k=0n-14rk1-nk1/2kk!2∫-ππ∫01∫b∞ak+1tkrt-a2+Z+cθ/2π2n-1/2rt+a2+Z+cθ/2π2k×1+4r2π2c2t21/2dadtdθ.The above integrals become extremely complex; therefore we leverage a standard integration package, such as MAPLE, to evaluate them numerically. Thus, the total interaction energy between the ssDNA molecule and the intermediate layers could be given by(13)EssDNAp=ηIηd-AI3p+BI6p.Furthermore, the interaction between the ssDNA molecule and a box which is assumed to have a thickness of l, where -l<z<0, can be identified by following similar steps undertaken with the term (Z+cθ/2π) replaced by (Z+cθ/2π-z); thus the integral In which is given in (3) becomes(14)Inb=2π∑k=0n-14rk1-nk1/2kk!2∫-l0∫-ππ∫01∫b∞ak+1tkrt-a2+Z+cθ/2π-z2n-1/2rt+a2+Z+cθ/2π-z2k×1+4r2π2c2t21/2dadtdθdz,where b denotes the box interaction, and the thickness of the lipid bilayer head group is l=4 Å, and l=15 for that of the tail group. Thus, the total interaction energy between the ssDNA molecule and the head or the tail group is given by(15)EssDNAb=ηiηd-AI3b+BI6b,where i can be either H or T.Figure 2
An ssDNA molecule penetrating a lipid bilayer gap.


### 3.2. Interaction Energy of dsDNA Penetrating a Lipid Bilayer Gap

In this subsection, the energy behavior of a unit of the dsDNA molecule is considered. With reference to Figure3, the dsDNA molecule is assumed to be located on the z-axis at a distance Z above the gap in the lipid bilayer. A typical point on the dsDNA has the following coordinates:(16)R2cosΘ+cosΘ-ϕ+tcosΘ-cosΘ-ϕ,R2sinΘ+sinΘ-ϕ+tsinΘ-sinΘ-ϕ,cΘ2π+Z.Here, the helical phase angle ϕ is assumed to be equal to π(ϕ=π), which is close to the physical value of the helical angle which is equal to 12π/17(ϕ=12π/17). We comment that this assumption can lead to a substantial simplification in the analytical expression for the interactions. Thus, the parametric equation of a typical surface element of the dsDNA molecule interacting with the lipid bilayer is given by RtcosΘ,RtsinΘ,cΘ/2π+Z. As before, without loss of generality, the coordinates of a typical point on the dsDNA molecule can be given by (λ,0,Z+cΘ/2π), where λ=Rt. Thus, the distance between two typical points on the surface of a dsDNA and the plane surface is given by(17)ρ2=λ-a2+4aλsin2ϑ2+Z+cΘ2π2,and the total interaction energy between the dsDNA molecule and the intermediate layers can be given by(18)EdsDNAP=ηIηD-AI3P+BI6P,where(19)InP=2π∑k=0n-14Rk1-nk1/2kk!2∫-ππ∫-11∫b∞ak+1tkRt-a2+Z+cΘ/2π2n-1/2Rt+a2+Z+cΘ/2π2k×1+4R2π2c2t21/2dadtdΘ.Moreover, the total energy of the interaction of a dsDNA molecule with the head or the tail group is given by(20)EdsDNAB=ηiηD-AI3B+BI6B,where i can be either H or T, and(21)InB=2π∑k=0n-14Rk1-nk1/2kk!2∫-l0∫-ππ∫-11∫b∞ak+1tkRt-a2+Z+cΘ/2π-z2n-1/2Rt+a2+Z+cΘ/2π-z2k×1+4R2π2c2t21/2dadtdΘdz.Figure 3
A dsDNA molecule penetrating a hole lipid bilayer.


## 3.1. Interaction Energy of ssDNA Penetrating a Lipid Bilayer Hole

A unit cell of the ssDNA molecule of radiusr is examined at a distance Z above a gap of radius b>r in the lipid bilayer, where Z is the perpendicular distance from the gap of the uppermost surface of the lipid bilayer to the center of the ssDNA molecule, as shown in Figure 2. Thus, a typical surface element of the ssDNA molecule has the coordinates given by rtcosθ,rtsinθ,cθ/2π+Z, and without loss of generality the coordinates of a typical point on the ssDNA molecule can be given by (α,0,Z+cθ/2π), where α=rt. Firstly, the interaction energy between an ssDNA molecule and a flat plane, which has the coordinates (acosϑ,acosϑ,0), is calculated. Therefore, the distance of the surface of an ssDNA to the plane is given by(7)ρ2=α-acosϑ2+-acosϑ2+Z+cθ2π2=α2-2αcosϑ+a2+Z+cθ2π2=α-a2+4aαsin2ϑ2+Z+cθ2π2.and from (3) the integral In can be written as(8)Inp=∫-ππ∫-ππ∫01∫b∞aα-a2+4aαsin2ϑ2+Z+cθ2π2-n×1+4r2π2c2t21/2dadtdθdϑ,where p is referred to in the plane interaction. By setting ξ=α-a2+Z+cθ/2π2 and β=4aα and defining the integral Kn as(9)Kn=∫-ππξ+βsin2ϑ2-ndϑ,upon making the substitution u=sin2(ϑ/2)⇒dϑ=u-1/21-u-1/2du, the integral Kn becomes(10)Kn=2∫01ξ+βu-nu-1/21-u-1/2du=2ξ-n∫011+βξu-nu-1/21-u-1/2du=2πξ-nFn,12;1;-βξ,where F(a∗,b∗;c∗;z∗) is the standard hypergeometric function [30]. Moreover, by using a Pfaff transformation [30], Kn becomes(11)Kn=2πξ-nξξ+β1/2F1-n,12;1;βξ+β=2πξ-nξξ+β1/2∑k=0n-11-nk1/2kk!2βξ+βk=2πrt-a2+Z+cθ2π2-n+1/2∑k=0n-14artk1-nk1/2kk!2rt+a2+Z+cθ/2π2k.Therefore, there are three remaining integrals for Inp:(12)Inp=2π∑k=0n-14rk1-nk1/2kk!2∫-ππ∫01∫b∞ak+1tkrt-a2+Z+cθ/2π2n-1/2rt+a2+Z+cθ/2π2k×1+4r2π2c2t21/2dadtdθ.The above integrals become extremely complex; therefore we leverage a standard integration package, such as MAPLE, to evaluate them numerically. Thus, the total interaction energy between the ssDNA molecule and the intermediate layers could be given by(13)EssDNAp=ηIηd-AI3p+BI6p.Furthermore, the interaction between the ssDNA molecule and a box which is assumed to have a thickness of l, where -l<z<0, can be identified by following similar steps undertaken with the term (Z+cθ/2π) replaced by (Z+cθ/2π-z); thus the integral In which is given in (3) becomes(14)Inb=2π∑k=0n-14rk1-nk1/2kk!2∫-l0∫-ππ∫01∫b∞ak+1tkrt-a2+Z+cθ/2π-z2n-1/2rt+a2+Z+cθ/2π-z2k×1+4r2π2c2t21/2dadtdθdz,where b denotes the box interaction, and the thickness of the lipid bilayer head group is l=4 Å, and l=15 for that of the tail group. Thus, the total interaction energy between the ssDNA molecule and the head or the tail group is given by(15)EssDNAb=ηiηd-AI3b+BI6b,where i can be either H or T.Figure 2
An ssDNA molecule penetrating a lipid bilayer gap.


## 3.2. Interaction Energy of dsDNA Penetrating a Lipid Bilayer Gap

In this subsection, the energy behavior of a unit of the dsDNA molecule is considered. With reference to Figure3, the dsDNA molecule is assumed to be located on the z-axis at a distance Z above the gap in the lipid bilayer. A typical point on the dsDNA has the following coordinates:(16)R2cosΘ+cosΘ-ϕ+tcosΘ-cosΘ-ϕ,R2sinΘ+sinΘ-ϕ+tsinΘ-sinΘ-ϕ,cΘ2π+Z.Here, the helical phase angle ϕ is assumed to be equal to π(ϕ=π), which is close to the physical value of the helical angle which is equal to 12π/17(ϕ=12π/17). We comment that this assumption can lead to a substantial simplification in the analytical expression for the interactions. Thus, the parametric equation of a typical surface element of the dsDNA molecule interacting with the lipid bilayer is given by RtcosΘ,RtsinΘ,cΘ/2π+Z. As before, without loss of generality, the coordinates of a typical point on the dsDNA molecule can be given by (λ,0,Z+cΘ/2π), where λ=Rt. Thus, the distance between two typical points on the surface of a dsDNA and the plane surface is given by(17)ρ2=λ-a2+4aλsin2ϑ2+Z+cΘ2π2,and the total interaction energy between the dsDNA molecule and the intermediate layers can be given by(18)EdsDNAP=ηIηD-AI3P+BI6P,where(19)InP=2π∑k=0n-14Rk1-nk1/2kk!2∫-ππ∫-11∫b∞ak+1tkRt-a2+Z+cΘ/2π2n-1/2Rt+a2+Z+cΘ/2π2k×1+4R2π2c2t21/2dadtdΘ.Moreover, the total energy of the interaction of a dsDNA molecule with the head or the tail group is given by(20)EdsDNAB=ηiηD-AI3B+BI6B,where i can be either H or T, and(21)InB=2π∑k=0n-14Rk1-nk1/2kk!2∫-l0∫-ππ∫-11∫b∞ak+1tkRt-a2+Z+cΘ/2π-z2n-1/2Rt+a2+Z+cΘ/2π-z2k×1+4R2π2c2t21/2dadtdΘdz.Figure 3
A dsDNA molecule penetrating a hole lipid bilayer.


## 4. Numerical Results and Discussions

As a part of the experiment, the interaction energies for the ssDNA and dsDNA molecules and the lipid bilayer are calculated and for the dsDNA we assume that the helical phase angleϕ is equal to π. To obtain the numerical results of the above systems, we use the algebraic computer package MAPLE together with the parameter values given in Tables 1 and 2. Firstly, we show the equilibrium positions with distance δ for the ssDNA and dsDNA molecules assuming that the DNA molecules are interacting with a lipid bilayer without a gap. Graphically in Figure 4 the relationship between the interaction energy and the distances δ is shown, where δ is referred to as the closest distance from the upper lipid head group to the DNA molecule and (δ=Z-17). The results demonstrate that for both the ssDNA and dsDNA molecules the equilibrium distances δ are ≈3.60 Å above the bilayer. In addition, due to the difference in the number of atoms between the ssDNA and dsDNA molecules, it can be seen that the value of the energy of the ssDNA molecule is less than the energy value of the dsDNA molecule. In addition, Figure 5 shows the relationship between the total interaction energy and the perpendicular distance Z for various values of the gap radius in the lipid bilayer. Results indicate that when the values of the gap radius of the lipid bilayer are bigger than ≈10.5 Å, the DNA molecule can penetrate the gap. In addition, the minimum energy occurs when the radius of the gap is equal to 10.5 Å. Moreover, Figure 6 shows that the dsDNA molecule can be moved through the gap radius when the value of the radius is larger than ≈14 Å. The optimal value of the radius for a dsDNA molecule which is provided with the minimum energy to penetrate the lipid bilayer is approximately 14.1 Å. Furthermore, for both the ssDNA and dsDNA molecules when the radius of the gap becomes larger, the perpendicular distance Z to the upper head group of the lipid bilayer becomes smaller. In addition, the results demonstrate that the DNA molecule requires a large amount of energy to move inside the bilayer. These results are in agreement with work of Khalid et al. [13].Table 1
Numerical values for constants used in this paper (∗ denotes data from [21], ∗∗ denotes data from [22], and ∗∗∗ denotes data from [18]).
Constant Value ∗Radius of dsDNA R = 10 Å ∗Length of DNA c = 34 Å ∗ ∗Radius of ssDNA r = 5.4  Å ∗ ∗ ∗ Atomic density of the head group η H = 0.0077  Å-3 ∗ ∗ ∗ Atomic density of the intermediate group η I = 0.0308  Å-2 ∗ ∗ ∗ Atomic density of the tail group η T = 0.0082  Å-3 Atomic density of ssDNA η d = 0.80  Å-2 Atomic density of dsDNA η D = 0.84  Å-2Table 2
Lennard-Jones parameters for DPPC lipid interacting with DNA [23, 24].
Interaction ε (kcal/mol) σ (Å) A (Å6 kcal/mol) B (Å12 kcal/mol) Head group 0.2718 4.7996 13540.92566 201188130.2 Intermediate layer 0.2567 4.04596 4989.875741 26282244.42 Tail group 0.2396 4.04596 4657.867035 24533516.72Figure 4
Energy profiles for ssDNA and dsDNA interacting with a lipid bilayer without a gap.Figure 5
Total energy for ssDNA molecule interacting with different values of the gap radiusb with respect to the perpendicular distance Z.
(a)
Gap radiusb=0, 2, 4, 6, 8, and 10 Å (b)
Gap radiusb=10.5-15 ÅFigure 6
Total energy for dsDNA molecule interacting with different values of the gap radiusb with respect to the perpendicular distance Z.
(a)
Gap radiusb=0, 2, 4, 6, 8, and 10 Å (b)
Gap radiusb=11, 12, and 13 Å (c)
Gap radiusb = 14-20 Å


## 5. Conclusion

This study introduces a mathematical modelling approach which is used to adopt the Lennard-Jones potential together with continuous approximation to investigate the interaction energy of single-stranded and double-stranded DNA molecules penetrating into a DPPC lipid bilayer. The numerical results were obtained by employing the algebraic computer package MAPLE. In the absence of a gap in the lipid, the distance from the edge of the ssDNA and dsDNA molecules to the upper head group is obtained where the minimum energy occurs. The results indicate that for both types of molecules the spacing is≈3.60 Å. Moreover, the energy behavior for both the ssDNA and dsDNA molecules is measured while penetrating a presumably circular gap on the infinite plane of a lipid bilayer. Results show that an ssDNA molecule is moved inside the bilayer when the gap radius is larger than 10 Å. In addition, a dsDNA molecule can be located inside the bilayer when the radius of the bilayer gap is larger than 13.8 Å. In addition, the results demonstrate that the potential energy required for both the ssDNA and dsDNA molecules to penetrate the lipid bilayer is high. These results compare favourably with other methods, such as the experiments and (MD) simulations, conducted by [10, 13]. Our study can provide insights into the molecular interactions between a DNA molecule and a lipid bilayer. Therefore, it can also provide direction for targeted drug and gene delivery. In addition, this study can provide additional clarification as to the interaction of DNA with lipids and of the mechanisms by which DNA might move through cells.



---

*Source: 7017416-2018-11-01.xml*