Next Article in Journal
Oxidized Low-Density Lipoprotein Promotes In Vitro Calcification
Next Article in Special Issue
LCAO Electronic Structure of Nucleic Acid Bases and Other Heterocycles and Transfer Integrals in B-DNA, Including Structural Variability
Previous Article in Journal
Epitaxial Growth of GaN on Magnetron Sputtered AlN/Hexagonal BN/Sapphire Substrates
Previous Article in Special Issue
Hole Transfer in Open Carbynes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Base-Pairs’ Correlated Oscillation Effects on the Charge Transfer in Double-Helix B-DNA Molecules

Departamento de Física de Materiales, Facultad CC. Físicas, Universidad Complutense de Madrid, E-28040 Madrid, Spain
Materials 2020, 13(22), 5119; https://doi.org/10.3390/ma13225119
Submission received: 26 October 2020 / Revised: 9 November 2020 / Accepted: 9 November 2020 / Published: 13 November 2020
(This article belongs to the Special Issue Computational Modeling and Simulation of Polymers and Biopolymers)

Abstract

:
By introducing a suitable renormalization process, the charge carrier and phonon dynamics of a double-stranded helical DNA molecule are expressed in terms of an effective Hamiltonian describing a linear chain, where the renormalized transfer integrals explicitly depend on the relative orientations of the Watson–Crick base pairs, and the renormalized on-site energies are related to the electronic parameters of consecutive base pairs along the helix axis, as well as to the low-frequency phonons’ dispersion relation. The existence of synchronized collective oscillations enhancing the π - π orbital overlapping among different base pairs is disclosed from the study of the obtained analytical dynamical equations. The role of these phonon-correlated, long-range oscillation effects on the charge transfer properties of double-stranded DNA homopolymers is discussed in terms of the resulting band structure.

1. Introduction

In physiological conditions, the DNA double helix exhibits a full-fledged three-dimensional (3D) geometry, where every two consecutive Watson–Crick base pairs (bps) stand nearly parallel to each other, and they are twisted by a certain angle ( θ 0 36 in equilibrium conditions) around the helix axis. In their pioneering work, Eley and Spivey pointed out that a double-stranded DNA (dsDNA) molecule might behave as a one-dimensional (1D) aromatic crystal displaying a π - π based electrical conductivity along the helical axis [1]. The reasoning behind this proposal was that dsDNA’s nucleobases adenine (A), guanine (G), cytosine (C), and thymine (T) are aromatic compounds whose atomic p z orbitals perpendicular to the plane of the base can form rather delocalized π bonding and π * antibonding molecular orbitals. If the orbital overlap between the bps is strong enough, this could lead to extended electronic states along the helical axis, thereby promoting charge transfer (CT) between consecutive bps in an efficient way over long distances through the aromatic bp stack within the DNA helix. Accordingly, CT depends on the intimate coupling among stacked bases, as determined by their relative separation and twist angle, and hence, any perturbation in that stacking, altering the optimal overlapping face-to-face configuration ( θ 0 = 0 ), will significantly affect DNA charge migration. Consequently, one expects structural fluctuations to be an important factor, influencing charge carriers’ transport through dsDNA molecules [2,3,4]. In fact, at physiological temperatures, the relative orientation of neighboring bases becomes a function of time, thereby modifying their mutual overlapping in an oscillatory way. The motion of bases can either occur in a synchronized manner (normal modes’ propagation) or incoherently. The role of thermal fluctuations on the CT efficiency has been discussed in a number of previous works, where the structural fluctuations of the DNA double helix are described by sampling the initial angular velocities and twist angles from a Boltzmann distribution at a given temperature [5,6,7,8,9,10,11,12,13]. Not surprisingly, it was found that the uncorrelated motion of bps, randomly twisting back and forth around the helix axis due to thermal fluctuations, generally reduces π - π stacking overlap, hence degrading the CT efficiency.
On the contrary, the presence of synchronized, collective twist motions of the Watson–Crick bps in DNA duplexes can efficiently enhance the π π orbital overlapping between non-consecutive bps via a long-range, phonon-correlated tunneling effect [14,15]. In this work, we will focus on coherent charge transport promoted by the coupling between both twist and radial vibration modes and charge motion through duplex DNA, thereby extending previously obtained results [16]. In order to analyze the interplay between the dsDNA low frequency bps’ dynamics and CT efficiency, we will study the coupling between the oscillations of complementary bases along the transversal direction and the twisting motion of each bp as a whole through the helical sugar-phosphate backbone structure, explicitly taking into account its characteristic helical geometry, which has been shown to be very important in biological processes, such as denaturation and transcription [17]. To this end, we will explicitly take into account the stacking interaction, mediated by the orbital overlapping between adjacent bps along the helix, as well as hydrogen bond stretch motions, as described in the Peyrard–Dauxois–Bishop (PDB) model [18,19,20], in the phonon dynamical equations. In our approach, the charge carrier dynamics through a helical dsDNA molecule is expressed in terms of an effective renormalized Hamiltonian describing a diatomic linear chain, where the renormalized transfer integrals explicitly depend on the relative orientations of the Watson–Crick bps and the renormalized on-site energies are related to both the electronic parameters of consecutive codon units along the helix axis, as well as the low-frequency phonon dispersion relation. The corresponding effective hopping terms include both helical and dynamical effects in an intertwined fashion, allowing for a unified treatment of charge-lattice coupled dynamics in a fully analytical way. Thus, we disclose a number of remarkable symmetries of the motion equations themselves, which may be implemented with accurate charge transfer parameters derived from quantum chemistry and molecular dynamics approaches in a straightforward way. Our main conclusion is that a significant improvement of CT can occur in dsDNA via charge-phonon coupling mediated by synchronized helical waves stemming from collective, long-range correlated bps’ oscillation modes.
The paper is organized as follows. In Section 2, we introduce the model Hamiltonian describing the lattice and electronic dynamics. The lattice contribution is expressed in cylindrical coordinates, in order to explicitly take into account the 3D geometry of the double helix DNA molecule. Then, we derive the π π electronic coupling term describing the molecular orbitals’ overlap through the helix axis in terms of these cylindrical coordinates. This term, describing the transfer integral between successive bps, allows one to relate the dynamical behavior to the CT process by means of a suitable tight-binding electronic model. A convenient feature of our adopted approach is that the resulting 3D fishbone model can be properly mapped into a mathematically simpler effective 1D chain model, still retaining much physico-chemical information in the corresponding renormalization parameters. In Section 3, we obtain the linearized canonical equations of motion for twist and radial lattice variables. In doing so, we introduce a number of characteristic frequencies along with their related time scales. For the sake of simplicity, in Section 3.2, we focus on the dynamics of homopolymer dsDNA molecules, showing the presence of collective oscillations in the form of helical waves. The related dispersion relations for the acoustic and optical branches are analytically derived, and the obtained results are compared to some available experimental results. In this section, we also disclose a very interesting relationship between the dynamics of the dsDNA molecule as a whole and that corresponding to its codon building blocks. In Section 4, we solve the Schrödinger equation for the effective 1D Hamiltonian previously introduced, making explicit use of the helical wave solutions in the transfer integral term. In this way, the charge-phonon coupling effect is fully incorporated in the resulting 3D energy spectrum. Finally, the main conclusions of this work are summarized in Section 5.

2. DNA Model Hamiltonian

Two kinds of order coexist in biological DNA, each one related to two separate subsystems in the DNA helix, namely the nucleobase and backbone systems [21]. The informative chemical order determined by the sequence of Watson–Crick bps can be suitably characterized by ab-initio quantum chemistry calculations [22,23], which properly highlight the emergence of molecular orbitals, as is shown in Figure 1b. In order to describe most basic properties of dsDNA molecules, we must consider a model Hamiltonian accounting for different scales of time and space by means of an adequate choice of generalized coordinates including both electronic and dynamic degrees of freedom. According to the Born–Oppenheimer approximation, the lattice and charge dynamics of a dsDNA molecule can be split in terms of the general Hamiltonian H = H l + H e , where H l describes the double strand dynamics, and H e describes the CT across nucleobases, as illustrated in the structural lattice model and the electronic tight-binding model depicted in Figure 1a,c, respectively.

2.1. Lattice Hamiltonian

In our lattice model, we treat each nucleotide (base + sugar + phosphate) as a point mass, helically arranged and mutually connected by means of elastic rods describing: (1) the sugar-phosphate backbone along a given strand and (2) the interstrand H-bonding between complementary bases (see Figure 1a) [24,25,26]. We explicitly take into consideration the mass difference among the four nucleobases, namely, m G = 347.05 , m C = 307.05 , m A = 331.06 , and m T = 322.05 amu, and so, we realize that the mass of each bp as a whole is essentially the same, i.e., M m G + m C m A + m T 653.5 ± 0.5 amu. Adopting the reference frame indicated in Figure 1a, the position of the nth nucleobase can be expressed as x n = r n cos φ n , y n = r n sin φ n , and z n = c φ n , where n labels the considered bp along the dsDNA, r n and φ n are the usual cylindrical coordinates, and c = h 0 / θ 0 , where h 0 0.34 nm is the equilibrium distance between two successive bp planes along the Z axis in the B-DNA form, while θ 0 is the equilibrium relative angular separation between neighboring bps. We note that in this model, the distance between successive bps along the Z direction is proportional to the twist angle. In this way, the helical structure is naturally preserved during the dynamical evolution, in agreement with dispersion relation data reported from inelastic X-ray scattering measurements [27]. Thus, we can express the Euclidean distance between adjacent bases along a given strand as:
d n , n ± 1 = c 2 θ n , n ± 1 2 + R 0 + ρ n ± 1 2 + R 0 + ρ n 2 2 R 0 + ρ n ± 1 R 0 + ρ n cos θ n , n ± 1
where we defined θ n , n ± 1 = ± φ n ± 1 φ n as the relative angle between two neighboring bps, and ρ n = r n R 0 is the radial displacement about the equilibrium position ( R 0 = 1 nm). We will further assume that the location of the dsDNA as a whole remains fixed, so that the center of mass is constant for each bp. Therefore, the radial displacements about the equilibrium position satisfy the relationship ρ ¯ n = λ n ρ n , where λ n m n / m ¯ n , and henceforth, the upper bar denotes the physical magnitudes of the complementary bases. Therefore, d ¯ n , n ± 1 for the opposite strand can be obtained by simply replacing ρ n with ρ ¯ n in Equation (1). In equilibrium, both distances reduce to the value l 0 = d n , n ± 1 e q . = d ¯ n , n ± 1 e q . = h 0 2 + 4 R 0 2 sin 2 θ 0 / 2 0.685 nm, where we adopted θ 0 = π / 5.2 34.6 .
When describing the phonon dynamics of DNA at a molecular scale, one can disregard the inner degrees of freedom of the bases, since we can separate the fast vibrational motions of atoms about their equilibrium positions from the slower motions of molecular groups. Accordingly, we can write the dsDNA molecule lattice Hamiltonian as [16]:
H l = 1 2 M n = 1 N P ρ n 2 λ n + P φ n 2 ξ 2 + λ n ρ n 2 + 4 R 0 ρ n m n M 1 + U H + U S + U B ,
where n runs over the number N of bps, P ρ n and P φ n are the conjugate momenta of the nth bp radial and twist variables, respectively, and ξ = c 2 + R 0 2 1.147 nm is related to the helical geometry of the system, so that ξ θ n , n ± 1 measures the helix arc length providing the shortest path between two points along a helical coil. In the limit of small radial and twist oscillations ( r n R 0 , θ n , n + 1 1 ) , Equation (1) reads d n , n ± 1 = R 0 2 + c 2 θ n , n ± 1 ξ θ n , n ± 1 , so that the Euclidean distance coincides with the helix arc length in this case [14].
The three elastic potential terms in Equation (2) describe the different interactions between the bases within the framework of the PDB model [18,19], namely:
U H = n = 1 N D n e α n 2 1 + λ n ρ n 1 2 ,
represents the radial stretching of the hydrogen bonds connecting complementary bases in the opposite strands of the double helix by means of Morse potentials of depth D n and width α n [19,28]. This potential term includes both the attraction due to the H-bonds forming the bps and the repulsion of the negatively charged phosphates in the backbone of the two strands, which is, in turn, screened by the surrounding solvent water molecules and positively charged counterions. The sequence dependence can be considered by adopting different values for the model parameters D n and α n , accounting for the different number of H-bonds in the G≡C and A=T bps [29]. The potential term:
U S = 1 8 n = 1 N 1 k n , n + 1 S 1 + e b 2 u n , n + 1 + u n , n + 1 2 ,
with u n , n + 1 ± = 1 + λ n ρ n ± 1 + λ n + 1 ρ n + 1 , describes the stacking interaction between adjacent bps, whose role is to inhibit configurations with large relative radial displacements between neighboring pairs. This interaction is characterized by the exponential term that effectively modulates an otherwise harmonic radial oscillation. This term accounts for local constraints in nucleotide motions, measured in terms of the stacking stiffness k n , n + 1 S and the interaction range b parameter, resulting in long-range cooperative elastic effects due to the distortion of hydrogen bonds and the overlap of the π -type orbitals [28]. The description of the radial degree of freedom via the non-linear potentials given by Equations (3) and (4) is more realistic than a purely harmonic approach and has been successful in capturing denaturation, as well as transcription initiation processes in several DNA model chains [19,30,31,32]. Finally, the term:
U B = k B 2 n = 1 N 1 d n , n + 1 l 0 2 + d ¯ n , n + 1 l 0 2 ,
describes the harmonic interaction between neighboring bases along each backbone’s strand.

2.2. The π - π Electronic Coupling

The radial and twist oscillations of bps have a significant impact on the molecular orbitals’ overlap throughout the π -stacking, so that the resulting electronic transfer integrals’ values explicitly depend on the dynamical degrees of freedom ρ n , ρ n ± 1 , and θ n , n ± 1 . The π and π * molecular orbitals are formed by the C, N, and O atomic p z orbitals perpendicular to the bps and pointing along the helical axis, as is illustrated in Figure 2. The p z orbitals from different bps couple by p p σ > 0 and p p π < 0 hybridization, the different signs arising from the respective atomic orbital’s lobe sign. According to the Slater–Koster theory, the transfer matrix element between two p z orbitals on neighboring bps is given by the combination of p p σ and p p π hybridization contributions as V i j = V p p σ sin 2 ζ + V p p π cos 2 ζ , where ζ measures the rise angle between successive bps ( sin ζ = h i j / d i j ), and:
V p p x = η p p x ħ 2 m d i j 2 exp ( d i j / R c )
where η p p π and η p p σ describe the hybridization matrix elements, m is the electron mass, d i j = l i j 2 + h i j 2 is the Euclidean distance between atoms belonging to neighboring nucleobases, and R c describes the exponential tails of the atomic wave functions [33,34].
The π - π coupling generally involves one base (say, Y) of the ( n + 1 )th bp and another base (say, X) of the nth bp, both belonging to the same strand (i.e., 5’-XY-3’), so that the π - π transfer integral between successive bps along the dsDNA helix is then given by:
t n , n + 1 X Y = i = 1 N 1 j = 1 N 2 V i j n , n + 1 c i n + 1 c j n ,
where N 1 and N 2 are the number of p z orbitals in the bps n + 1 and n, respectively, and c j n is the jth linear combination of atomic orbitals (LCAO) coefficient of the π molecular frontier orbital (HOMO or LUMO) of bp n. Making use of Equation (6) in Equation (7), assuming that the mean distances among atoms belonging to different nucleobases can be roughly approximated as d i j d and l i j 2 l n , n ± 1 2 = r n 2 + r n ± 1 2 2 r n r n ± 1 cos θ n , n ± 1 , respectively, we obtain: [14]
t n , n ± 1 X Y ( ρ , θ ) = t 0 X Y 1 η ¯ d 2 R 0 + ρ n 2 + R 0 + ρ n ± 1 2 2 R 0 + ρ n ± 1 R 0 + ρ n cos θ n , n ± 1 ,
where:
t 0 X Y = η p p σ ħ 2 m d 2 exp ( d / R c ) i = 1 N 1 j = 1 N 2 c i n + 1 c j n ,
is the transfer integral corresponding to the optimal face-to-face geometry (i.e., θ n , n ± 1 0 , ρ n = 0 , n ) and η ¯ 1 + | η p p π | / η p p σ = 1 + 2.26 / 5.27 1.429 [35]. In the B-DNA form equilibrium configuration ( θ n , n ± 1 θ 0 , ρ n = 0 , n ), Equation (8) reads:
t n , n ± 1 X Y ( 0 , θ 0 ) = t 0 X Y 1 η ¯ 2 R 0 l 0 sin θ 0 2 2 .
Since η ¯ > 0 , we get t n , n ± 1 X Y ( 0 , θ 0 ) < t 0 X Y , hence indicating that by explicitly considering helical geometry in the equilibrium configuration, the π π base coupling strength is significantly reduced below that corresponding to the optimal face-to-face geometry. If we relax the equilibrium structure, allowing for the propagation of low frequency twist oscillations (acoustic modes), though keeping the radial variable describing H-bonding stretch oscillations fixed (no optical modes), Equation (8) can be approximated as:
t n , n ± 1 X Y ( 0 , θ n , n ± 1 ) t 0 X Y 1 χ θ n , n ± 1 2 ,
for small enough twists ( cos θ n , n ± 1 1 θ n , n ± 1 2 / 2 ), where the dimensionless parameter χ η ¯ ( R 0 / l 0 ) 2 2.92 measures the electron-phonon coupling strength. Despite its approximate nature, Equation (11) reasonably reproduces the main features of the transfer integral versus twist angle dependence derived from detailed quantum-chemistry calculations in the regime of low energies [12,14]. If we now allow for radial oscillations, still keeping within the small twist angle regime, Equation (8) adopts the form:
t n , n ± 1 X Y ( v n , n ± 1 , θ n , n ± 1 ) t 0 X Y 1 χ v n v n ± 1 2 + ( 1 + v n ) ( 1 + v n ± 1 ) θ n , n ± 1 2 ,
where v n ρ n / R 0 . Finally, in the limit v n 1 , Equation (12) can be approximated as:
t n , n ± 1 X Y ( v n , , n ± 1 , θ n , n ± 1 ) t 0 X Y 1 χ v n v n ± 1 2 + θ n , n ± 1 2 ,
where only terms up to the second order are retained.

2.3. Electronic Hamiltonian

In order to obtain a realistic description of the rich dsDNA physico-chemistry, keeping at the same time the convenient mathematical simplicity, we will exploit the three-step renormalization approach sketched in Figure 3.
In the fist step, the Watson–Crick bps present in the triplet codon shown in Figure 3a are renormalized to obtain the tight-binding model depicted in Figure 3b. The renormalized on-site energies and transfer integrals are respectively given by [36,37]: ε i j = t i j = ε j i , which describes the charge carrier hopping from one base to its complementary base within the same bp [38], and:
τ j = t P + ε j t P ( E γ j ) , j = { G , C , A , T } ,
are the effective transfer integrals between the bps and the sugar-phosphate groups, where t P is the glycosidic bond transfer integral, E is the charge carrier energy, γ j measure the sugar-phosphate groups’ on-site energies, and ε j are the on-site energies of the corresponding nucleobases. In general, γ j will depend on the nature of the neighboring base, as well as the presence of water molecules and/or counterions attached to the backbone. Thus, the renormalized model parameters ε i j and τ j entail substantial physicochemical information concerning nucleotide interactions and backbone gating effects [36,39].
We note that the model depicted in Figure 3b can be properly regarded as a 3D generalization of the so-called fishbone model in 2D [40,41]. Accordingly, in the second renormalization step, the sugar-phosphate groups contribution is decimated [36,40,42], so that the original dsDNA molecule is mapped into the equivalent 1D binary lattice shown in Figure 3c, where the renormalized on-site energies labeled α and β correspond to the Watson–Crick complementary bps, and they are explicitly given by ( E γ j ):
ε ˜ n X Y ( E , ρ , θ ) = t n , n ± 1 X Y ( ρ , θ ) + τ X 2 E γ X + τ Y 2 E γ Y t n , n ± 1 X Y ( ρ , θ ) + ϵ n X Y ( E ) ,
where t n , n ± 1 X Y ( ρ , θ ) accounts for the aromatic base stacking between adjacent nucleotides, given by Equations (8)–(13). Therefore, the renormalized on-site potentials ε ˜ n X Y ( E ) explicitly depend on the charge carrier energy, as well as on the angular and radial coordinates describing dsDNA oscillations, thereby enclosing all the relevant physicochemical information of the considered system. In this way, one obtains a realistic description, including 14 electronic model parameters, { ε j , t j , γ j , t G C , t A T }, fully describing CT throughout dsDNA molecules in terms of just two main functions, namely, ϵ n X Y ( E ) and t n , n ± 1 X Y ( ρ , θ ) , in a unified way in terms of the effective 1D Hamiltonian:
H ˜ e 1 D = n = 1 N t n , n + 1 X Y ( ρ , θ ) + ϵ n X Y ( E ) c n c n n = 1 N 1 t n , n + 1 X Y ( ρ , θ ) ( c n + 1 c n + c n c n + 1 ) ,
where c n ( c n ) is the creation (annihilation) operator for a charge at the nth site in the chain. In this way, the Hamiltonian given by Equation (16) provides a realistic treatment of CT mechanisms in dsDNA under physiological conditions, properly taking into account the influence of the dynamical state of the macromolecule on the CT efficiency. It is worth noting that epigenetic processes such as methylation (the addition of a methyl group (-CH 3 ) to one of the bases) will modify the on-site nucleobase energy and its effective mass alike, thereby changing both the mass ratio parameter λ among nucleotides and their γ parameter value. While the role of methylation-related on-site energy changes has been discussed in several recent works [43,44,45], the role of methylation-related dynamical effects in the CT efficiency has not. Due to the presence of both an on-site energy term ( ϵ n X Y ( E ) ) and a transfer integral term ( t n , n ± 1 X Y ( ρ , θ ) ) in the diagonal term of the effective Hamiltonian given by Equation (16), one should expect the possible existence of resonance effects involving both electronic and dynamical physical parameters in an intertwined fashion.

3. Dynamical Equations of Motion

3.1. General Expressions

From the lattice Hamiltonian given by Equation (2), we can straightforwardly obtain the canonical equations of motion:
ξ 2 + P ( ρ n ) φ ¨ n = k B M 1 l 0 d n , n 1 f n 1 1 l 0 d n , n + 1 f n + 1 + 1 l 0 d ¯ n , n 1 f ¯ n 1 1 l 0 d ¯ n , n + 1 f ¯ n + 1 ,
where P ( ρ n ) = λ n ρ n 2 + 4 m n M R 0 ρ n λ n ρ n 2 + 2 R 0 ρ n and:
f n ± 1 ( R 0 2 + R 0 ( ρ n + ρ n ± 1 ) + ρ n ρ n ± 1 ) sin θ n , n ± 1 + c 2 θ n , n ± 1 , f ¯ n ± 1 ( R 0 2 + R 0 ( λ n ρ n + λ n ± 1 ρ n ± 1 ) + λ n λ n ± 1 ρ n ρ n ± 1 ) sin θ n , n ± 1 + c 2 θ n , n ± 1 ,
along with:
ρ ¨ n = D n α n M 1 + λ n λ n ( e α n 2 ( 1 + λ n ) ρ n 1 ) e α n 2 ( 1 + λ n ) ρ n + k S 8 M 1 + λ n λ n u n 1 , n 2 + e b 2 u n , n 1 + 2 + b 2 u n 1 , n + u n + 1 , n 2 + e b 2 u n + 1 , n + 2 + b 2 u n + 1 , n k B M λ n 1 1 l 0 d n , n 1 g n 1 + λ n 1 1 l 0 d n , n + 1 g n + 1 + 1 l 0 d ¯ n , n 1 g ¯ n 1 + 1 l 0 d ¯ n , n + 1 g ¯ n + 1 ,
where:
g n ± 1 R 0 + ρ n ( R 0 + ρ n ± 1 ) cos θ n , n ± 1 , g ¯ n ± 1 R 0 + λ n ρ n ( R 0 + λ n ρ n ± 1 ) cos θ n , n ± 1 ,
and as a first approximation, we assumed all the stacking stiffness parameters to take on the same value (i.e., k n , n + 1 S = k S n ). In obtaining Equations (17) and (18), we neglected non-linear contributions related to the φ ˙ n ρ ˙ n and φ ˙ n 2 terms. Keeping only linear terms in the Taylor series of the functions appearing in the above expressions, we get the linearized equations of motion:
φ ¨ n + ω φ 2 2 φ n φ n 1 φ n + 1 = ω φ 2 l B 1 u n + 1 , n 1 ,
where we reasonably assumed 2 R 0 ρ n ξ 2 and introduced the twist frequency:
ω φ 2 2 k B M f 0 ξ l 0 2 ,
along with the characteristic length l B 2 f 0 / g 0 , with f 0 = c 2 θ 0 + R 0 2 sin θ 0 0.759 nm 2 and g 0 = R 0 1 cos θ 0 0.177 nm, so that l B 8.575 nm (i.e., about 25 bps) and:
ρ ¨ n + ω φ ρ , n 2 ρ n 1 2 ω φ S , n 1 2 ρ n 1 1 2 ω φ S , n + 1 2 ρ n + 1 = a B 2 ( 1 + λ n 1 ) ω φ 2 φ n 1 φ n + 1 ,
where a B g 0 ξ 2 / f 0 0.307 nm is a characteristic length whose value is comparable to the equilibrium bps separation h 0 . We introduced the coupled frequencies:
ω φ ρ , n 2 ( 1 + λ n ) 2 λ n ( ω H , n 2 + ω S 2 ) + 1 + λ n 2 λ n b B 2 ω φ 2 ,
and:
ω φ S , n ± 1 2 ( 1 + λ n ) ( 1 + λ n ± 1 ) λ n ω S 2 1 + λ n λ n ± 1 λ n b B 2 ω φ 2 ,
where b B a B / ξ = g 0 ξ / f 0 0.267 is a dimensionless factor. Therefore, ω φ ρ , n depends on the twist frequency, as well as the radial stretch H-bonding and lateral stacking oscillations of bps, whose characteristic frequencies are ω H , n 2 = D n α n 2 / ( 2 M ) , and ω S 2 = k S / M , respectively. The site label in ω φ ρ , n arises from the presence of the λ n factor, as well as the fact that ω H , n 2 is site dependent, due to the different Morse potential parameters values for G≡C and A=T bps. On the other hand, ω φ S , n ± 1 involves lateral stacking and twist oscillations. In this case, the site label dependence involves all the λ k terms. The set of coupled Equations (19) and (21) describes the dynamics of general dsDNA molecules, where two kinds of bps can be arranged either periodically or aperiodically [21,46,47,48,49].

3.2. Dynamics of Homopolymer dsDNA Macromolecules

For the sake of simplicity, we will consider in this section the homopolymer case (i.e., polyA-polyT or polyG-polyC chains). In this case, the renormalized chain shown in Figure 3c becomes an effective monoatomic lattice (i.e., α β ), where the renormalized on-site potentials depend on both the electron energy E and the phonon wavevector q, due to the presence of the π - π transfer integral in Equation (15). In addition, λ n = λ X Y , D n = D X Y D and α n = α X Y α (see Table 1), so that the frequencies ω H , n ω H , ω φ ρ , n ω φ ρ and ω φ S , n ± 1 ω φ S are no longer site dependent. Hence, Equations (19)–(21) can be rewritten as (henceforth, we will drop the subscript XY in the λ parameters for the sake of clarity):
φ ¨ n + ω φ 2 2 φ n φ n 1 φ n + 1 = A λ ω φ 2 ( ρ n + 1 ρ n 1 ) ,
ρ ¨ n + 1 2 ω φ ρ 2 ( 2 ρ n + ρ n 1 + ρ n + 1 ) 1 2 ω H S 2 ( ρ n 1 + ρ n + 1 ) = B λ ω φ 2 φ n 1 φ n + 1 ,
where A λ ( 1 + λ ) l B 1 and B λ a B ( 1 + λ 1 ) / 2 are constants, and we introduced the frequency ω H S 2 μ ( ω H 2 + 2 ω S 2 ) , describing the coupling between the stacking and H-bond radial stretch oscillations, where μ λ 1 ( 1 + λ ) 2 can take on two values. Making use of the model parameters listed in Table 1, we obtain the values listed in Table 2 for the characteristic frequencies just introduced, along with their related time scales.
From the data listed in Table 2, we see that the time scale of angular motions, determined by the twist frequency, amounts to ∼6 ps, which are an order of magnitude slower than those corresponding to the twist-radial coupled oscillations. The time scale related to bps’ H-bond and stacking motions occupy an intermediate position, whereas coupled oscillations involving stretch and stacking motions are about 2.6 times quicker than the H-bond-mediated stretch oscillations alone. For the sake of comparison, the transition times reported for intrastrand hole transfer in ds-GT n GGG oligonucleotides range from τ = 0.5 ps for n = 1 to τ = 315 ps for n = 4 [63]. Quite interestingly, the electrical response of biological dsDNA chains to light irradiation has been recently investigated in order to engineer a DNA based molecular switch. In these experiments, it was observed that the electrical current turns on when the frequency of the incident light is above the 2 THz threshold [64], a value that coincides with that listed for ν φ ρ in Table 2. On the other hand, it is worth mentioning that the ν H S frequency value listed in Table 2 is smaller than the 2.83 and 3.04 THz frequencies experimentally observed by optical Kerr-effect spectroscopy in ds-GGCGGCCCGCGCGGGCCGCC and ds-ATTATTATTATATTA oligonucleotides, respectively. These oscillations are consistent with delocalized optical phonon modes with a wavelength extending throughout the molecule as a whole, and they are assumed to be related to the H-bonding dynamics between the two strands [65,66,67].
The mathematical structure of Equations (24) and (25) clearly indicates the correlated nature of next-neighboring bps dynamics. It is then convenient to zoom out our perspective and focus our attention on the dynamics of consecutive triplets of bps along the double-helix, which are closely related to the so-called codon units in genomics. To this end, we properly add up the dynamical equations corresponding to the consecutive bps corresponding to sites n 1 , n, and n + 1 , grouping the resulting expression in terms of the collective variables x n 2 φ n φ n 1 φ n + 1 and y n 2 ρ n ρ n 1 ρ n + 1 , to obtain:
x ¨ n + ω φ 2 2 x n x n 1 x n + 1 = A λ ω φ 2 ( y n + 1 y n 1 ) ,
y ¨ n + 1 2 ω φ ρ 2 ( 2 y n + y n 1 + y n + 1 ) 1 2 ω H S 2 ( y n 1 + y n + 1 ) = B λ ω φ 2 x n 1 x n + 1 .
Remarkably enough, we realize that Equations (26) and (27), describing the codon dynamics as a whole, are formally identical to Equations (24) and (25), which describe the motion of their constituent bps, since they are invariant upon the simultaneous variable exchange φ n x n and ρ n y n . This property can be regarded as expressing a renormalization of the dynamical equations when going from the bp local scale to the longer triplet codon scale. Quite interestingly, the very mathematical structure of the codon dynamics prescribed by Equations (26) and (27) guarantees that we will obtain similar dynamical equations by grouping codons in successive triplets of nested codon units recurrently, all the way up to the entire DNA molecule itself. Accordingly, the set of Equations (24) and (25) exhibits a self-similar symmetry upon triplet renormalization operation, so that by solving this fundamental dynamical equation set, we are actually disclosing the main features of the dynamics of the entire dsDNA macromolecule as a whole.
Inspired by previous results [14,15], we look for solutions to Equations (24) and (25) of the form φ n = φ 0 e i ( ω t n q ξ ) and ρ n = ρ 0 e i ( ω t n q ξ ) , describing a helical wave propagating throughout the dsDNA with frequency ω and wave vector q, where φ 0 8 = 0.14 rad and ρ 0 0.05 nm are the twist and radial oscillation amplitudes at ambient temperature, respectively [68]. In so doing, Equations (24) and (25) can be expressed in the matrix form:
2 ω φ 2 [ 1 cos ( q ξ ) ] ω 2 2 i A λ ω φ 2 sin ( q ξ ) 2 i B λ ω φ 2 sin ( q ξ ) ω φ ρ 2 ω φ S 2 cos ( q ξ ) ω 2 φ 0 ρ 0 = 0 0 .
The solution to Equation (28) requires the matrix determinant to identically vanish, thereby leading to a biquadratic equation whose solutions yield the dispersion relations for the acoustic and optical phonon branches given by:
ω ± 2 = 1 2 G ( q ) + H ( q ) ± 1 2 H ( q ) G ( q ) 2 + 16 C λ ω φ 4 sin 2 ( q ξ ) ,
where G ( q ) 4 ω φ 2 sin 2 ( q ξ / 2 ) , H ( q ) ω φ ρ 2 ω φ S 2 cos ( q ξ ) , and C λ A λ B λ . The exact dispersion relations given by Equation (29) are plotted in Figure 4, though they can be very well approximated by the simpler expressions:
ν 2 = 4 ν φ 2 sin 2 q ξ 2 ν + 2 = μ ν H 2 + 2 ν S 2 sin 2 q ξ 2 ,
which extend previously reported results [25,30]. As we see, the acoustic branch is completely determined by twist oscillations, whereas the optical branch depends on both stretch and stacking oscillations. In particular, the q = 0 bandgap, Δ ν ( 0 ) ν + ( 0 ) ν ( 0 ) = μ ν H 1.70 THz (≃7 meV), is fixed by the H-bond frequency value. The maximum bandgap width occurs for q * = π / ξ 2.732 nm 1 , with Δ ν ( q * ) = ν H S 2 ν φ 1.91 THz. The sound velocity obtained from the acoustic dispersion curve is 1.2 km s 1 , a figure smaller than the experimentally reported values ranging from 1.7 to 4.3 km s 1 , depending on the employed technique [69].
An interesting subject that can be addressed within the framework presented in this work refers to the behavior of the correlated oscillations when the DNA molecule is bonded with other small molecules. Generally speaking, the presence of ligands attached onto the sugar-phosphate backbone simultaneously affects the electronic properties and the mass distribution in a local region of the DNA molecule. This twofold effect can be properly accounted for in terms of a change in the γ value parameter, which modifies the electronic bandgap according to Equation (40) and a change in the values of parameters λ and M, thereby modifying the frequency values given by Equations (20)–(23). Thus, any effective bp mass increase leads to a slow down of frequencies ω φ , ω H , and ω S , along with their related coupled frequencies. In addition, the reduction of the twist frequency value makes the acoustic dispersion relation slope decline, so that the sound speed is reduced as well, ultimately leading to a lower thermal conductivity around the place where the small molecule has been bonded. On the other hand, a smaller ω φ value will widen the gap between the acoustic and optical branches (see Figure 4).

4. Charge Transfer through dsDNA Homopolymers

Making use of the helical wave complex expressions for the variables φ n ( t ) and ρ n ( t ) into Equation (13), we get:
t n , n ± 1 X Y ( ρ n , θ n , n ± 1 ) t 0 X Y 1 4 χ A 0 2 sin 2 q ξ 2 ,
where A 0 2 φ 0 2 + ( ρ 0 / R 0 ) 2 0.02 , so that 4 χ A 0 2 0.26 eV and t n , n ± 1 X Y > 0 , q . The resulting transfer integrals become site independent, which is a natural consequence of the synchronized motion mediated by helical waves propagating in the dsDNA molecule. Accordingly, we can write t n , n ± 1 X Y ( ρ n , θ n , n ± 1 ) = t X Y ( q ) t X X ( q ) in the homopolymer case, and the Schrödinger nearest-neighbor tight-binding equation corresponding to the electronic Hamiltonian given by Equation (16) reads:
[ E ϵ X Y ( E ) t X X ( q ) ] ψ n t X X ( q ) ( ψ n + 1 + ψ n 1 ) = 0 ,
where ψ n is the electronic wave function at site n. Equation (32) can be expressed in the matrix form:
ψ n + 1 ψ n = E ϵ X Y ( E ) t X X ( q ) 1 1 1 0 ψ n ψ n 1 M ( E , q ) ψ n ψ n 1 ,
and within the framework of the transfer matrix formalism, the charge carrier dispersion relation of a dsDNA molecule containing N bps is given by (assuming periodic boundary conditions):
cos ( κ N h 0 ) = 1 2 tr M N ( E , q ) 1 2 tr n = N 1 M ( E , q ) = 1 2 tr M N ( E , q ) ,
where κ is the charge carrier wavevector. Since det M ( E , q ) = 1 , we can use the Cayley–Hamilton theorem for unimodular matrices in order to calculate the required power matrix as M N ( E , q ) = U N 1 M ( E , q ) U N 2 I , where I is the identity matrix and U m ( x ) sin [ ( m + 1 ) ϕ ] / sin ϕ , with x cos ϕ tr M / 2 , are Chebyshev polynomials of the second kind, satisfying the recurrence relationship U m + 1 2 x U m + U m 1 = 0 . In this way, we obtain [21]:
M N ( E , q ) = U N 1 2 x 1 1 0 U N 2 1 0 0 1 = U N U N 1 U N 1 U N 2 .
Taking into account the relationship U m U m 2 = 2 cos ( m ϕ ) , the charge carrier dispersion relation can then be expressed as:
E = ϵ X Y ( E ) + t X X ( q ) ( 1 + 2 cos ( κ h 0 ) ) .
According to Equation (15), the ϵ X Y ( E ) function entails detailed information regarding the electronic structure model of the dsDNA molecule. X-ray experiments indicated that counterions condense around the nucleic acid chain in a tightly-bound layer, in agreement with early model calculations [70]. Therefore, a homogeneous charge distribution through the backbone can be assumed as a first approximation, that is γ j γ . In that case, we can write:
ϵ X Y ( E ) = 2 t P 2 + a X Y ( E γ ) + b X Y ( E γ ) 2 ( E γ ) 1 ,
where a X Y ε X + ε Y and b X Y ( ε X 2 + ε Y 2 ) / ( 2 t P 2 ) . Thus, plugging Equations (31) and (37) into Equation (36), we obtain the charge carrier dispersion relation in the explicit polynomial form:
B 2 X Y E 2 + B 1 X Y + F ( κ , q ) E + 2 B 0 X Y γ F ( κ , q ) = 0 ,
where B 2 X Y 2 b X Y 1 , B 1 X Y γ ( 1 4 b X Y ) + 2 a X Y , and B 0 X Y γ ( γ b X Y a X Y ) + t P 2 , and we introduced the auxiliary function:
F ( κ , q ) t 0 X X 1 4 χ A 0 2 sin 2 q ξ 2 [ 1 + 2 cos ( κ h 0 ) ] .
The resulting energy spectrum structure consists of two slightly asymmetric bands, E ± ( κ , q ) , with relatively small widths ( W ± ), separated by a gap, E g ( 0 , q ) , whose value depends on the phonon wavevector q, as is illustrated in Figure 5 and Table 3 for the particular case γ = 0 . By inspecting Figure 5, we can clearly appreciate that the phonon coupling gives rise to a systematic reduction of the bandgap width E g ( 0 , 0 ) as q increases up to the value q * = π / ξ (see Figure 4), thereby enhancing the CT efficiency. Thus, the bandgap relative variation amounts to about 9% ( 6.5 %) for polyG-polyC (polyA-polyT), respectively, as compared to the E g ( 0 , 0 ) value. We note that the bandgap values listed in Table 3 are remarkably smaller than those usually reported in other studies [71,72,73,74,75,76]. Indeed, the dsDNA electronic structure is very sensitive to the precise value of the sugar-phosphate on-site energy. Thus, for the general case γ 0 , the bandgap width can be explicitly expressed as:
E g ( 0 , 0 ) = 1 2 b X Y 1 ( 3 t 0 X X γ + 2 a X Y ) 2 8 ( 2 b X Y 1 ) t P 2 ,
so that a semiconductor-semimetal transition can be promoted by properly tuning the adopted γ value [36].

5. Conclusions

The physical picture inspiring this work relies on the fact that the presence of collective, orchestrated oscillation motions of bps within the DNA double helix structure can efficiently enhance the π π orbital overlapping between bps along the backbone chain, hence promoting charge transfer via a long-range, phonon-correlated tunneling effect involving bps, which are relatively far apart. This property is intimately related to the helical geometry of the nucleobases’ arrangement along the duplex chain, which makes possible that a local hopping process, involving a relatively small number of neighboring nucleotides, ultimately extends over the entire DNA chain as a consequence of the synchronized nature of the resulting helical wave. The possible presence of these helical waves may be relevant in the study of CT properties in dsDNA polymers exhibiting extensive chemically homogeneous regions along their strands, such as those reported in genomic studies of the telomere sequences of certain invertebrates [77] or in tandem repeats [78], involving mononucleotide triplet motifs (AAA, TTT, GGG, or CCC) [79]. Indeed, the existence of DNA-mediated charge migration has been related to the understanding of the damage recognition process, including the presence of lesions and mismatches. Furthermore, since CT dependence can be sensed electrically, it can be exploited for nanotechnological applications through base modifications or DNA-protein binding or with the task of designing nanoscale sensing of genomic mutations, opening new challenges for emerging nanobiotechnologies [80,81,82,83]. In addition, the fundamental dynamical mechanisms reported in this work are expected to also take place in other π π molecular wires, such as G based quadruplexes [84].

Funding

This research received no external funding.

Acknowledgments

I sincerely thank Constantinos Simserides for his continued interest in my research works and Evgeni B. Starikov for illuminating conversations on the physico-chemical aspects of DNA molecules. I thank Victoria Hernández for a critical reading of the manuscript.

Conflicts of Interest

The author declares no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
dsDNAdouble-stranded DNA
bpbase pair
CTcharge transfer
Aadenine
Ccytosine
Gguanine
Tthymine
PDBPeyrard–Dauxois–Bishop
HOMOhighest occupied molecular orbital
LUMOlowest occupied molecular orbital
LCAOlinear combination of atomic orbitals

References

  1. Eley, D.D.; Spivey, D.I. Semiconductivity of organic substances. Part 9. Nucleic Acid in the Dry State. Trans. Faraday Soc. 1962, 58, 411–415. [Google Scholar] [CrossRef]
  2. Starikov, E.B.; Fujita, T.; Wanatabe, H.; Sengoku, Y.; Tanaka, S.; Wenzel, W. Effects of molecular motion on charge transfer/transport through DNA duplexes with and without base pair mismatch. Mol. Simul. 2006, 32, 759–764. [Google Scholar] [CrossRef]
  3. Berlin, Y.A.; Grozema, F.C.; Siebbeles, L.D.A.; Ratner, M.A. Charge transfer in donor-bridge-acceptor systems: Static disorder, dynamic fluctuations, and complex kinetics. J. Phys. Chem. C 2008, 112, 10988–11000. [Google Scholar] [CrossRef]
  4. Starikov, E.B.; Quintilla, A.; Nganou, C.; Lee, K.H.; Cuniberti, G.; Wenzel, W. Single-molecule DNA conductance in water solutions: Role of DNA low-frequency dynamics. Chem. Phys. Lett. 2009, 467, 369–374. [Google Scholar] [CrossRef]
  5. Bruinsma, R.; Grüner, G.; D’Orsogna, M.R.; Rudnick, J. Fluctuation-facilitated charge migration along DNA. Phys. Rev. Lett. 2000, 85, 4393–4396. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Roche, S. Sequence Dependent DNA-Mediated Conduction. Phys. Rev. Lett. 2003, 91, 108101. [Google Scholar] [CrossRef]
  7. Zhu, Y.; Kaun, C.C.; Guo, H. Contact, charging, and disorder effects on charge transport through a model DNA molecule. Phys. Rev. B 2004, 69, 245112. [Google Scholar] [CrossRef]
  8. Ren, W.; Wang, J.; Ma, Z.; Guo, H. Effect of thermal fluctuations of twist angles on charge transport in DNA: A model calculation. Phys. Rev. B 2005, 72, 035456. [Google Scholar] [CrossRef]
  9. Guo, A.M.; Xu, H. Sequence dependent charge transport through DNA molecules: The role of periodicity and long-range correlations. Physica B 2007, 391, 292–298. [Google Scholar] [CrossRef]
  10. Gutiérrez, R.; Mandal, S.; Cuniberti, G. Dissipative effects in the electronic transport through DNA molecular wires. Phys. Rev. B 2005, 71, 235116. [Google Scholar] [CrossRef] [Green Version]
  11. Gutiérrez, R.; Mohapatra, S.; Cohen, H.; Porath, D.; Cuniberti, G. Inelastic quantum transport in a ladder model: Implications for DNA conduction and comparison to experiments on suspended DNA oligomers. Phys. Rev. B 2006, 74, 235105. [Google Scholar] [CrossRef] [Green Version]
  12. Lemaur, V.; da Silva Filho, D.A.; Coropceanu, V.; Lehmann, M.; Geerts, Y.; Piris, J.; Debije, M.G.; van de Craats, A.M.; Senthilkumar, K.; Siebbeles, L.D.A.; et al. Charge transport properties in discotic liquid crystals: A quantum-chemical insight into structure-property relationships. J. Am. Chem. Soc. 2004, 126, 3271–3279. [Google Scholar] [CrossRef] [PubMed]
  13. Gutiérrez, R.; Caetano, R.; B Woiczikowski, P.B.; Kubar, T.; Elstner, M.; Cuniberti, G. Structural fluctuations and quantum transport through DNA molecular wires: A combined molecular dynamics and model Hamiltonian approach. New J. Phys. 2010, 12, 023022. [Google Scholar] [CrossRef]
  14. Maciá, E. Electrical conductance in duplex DNA: Helical effects and low-frequency vibrational coupling. Phys. Rev. B 2007, 76, 245123. [Google Scholar] [CrossRef] [Green Version]
  15. Maciá, E. π-π orbital resonance in twisting douplex DNA: Dynamical phyllotaxis and electronic structure effects. Phys. Rev. B 2009, 80, 125102. [Google Scholar] [CrossRef] [Green Version]
  16. Torrellas, G.; Maciá, E. Twist–radial normal mode analysis in double-stranded DNA chains. Phys. Lett. A 2012, 376, 3407–3410. [Google Scholar] [CrossRef]
  17. Peyrard, M. Nonlinear dynamics and statistical physics of DNA. Nonlinearity 2004, 17, R1–R40. [Google Scholar] [CrossRef]
  18. Peyrard, M.; Bishop, A.R. Statistical mechanics of a nonlinear model for DNA denaturation. Phys. Rev. Lett. 1989, 62, 2755–2758. [Google Scholar] [CrossRef]
  19. Dauxois, T.; Peyrard, M. Entropy-driven transition in a one-dimensional system. Phys. Rev. E 1995, 51, 4027–4040. [Google Scholar] [CrossRef]
  20. Peyrard, M.; Cuesta-López, S.; Angelov, D. Experimental and theoretical studies of sequence effects on the fluctuation and melting of short DNA molecule. J. Phys. Condens. Matter 2009, 21, 034103. [Google Scholar] [CrossRef] [Green Version]
  21. Maciá-Barber, E. Aperiodic Structures in Condensed Matter: Fundamentals and Applications; Taylor and Francis, CRC Press: Boca Raton, FL, USA, 2009; pp. 209–297. [Google Scholar]
  22. Shinwari, W.; Deen, M.J.; Starikov, E.B.; Cuniberti, G. Electrical Conductance in Biological Molecules. Adv. Funct. Mater. 2010, 20, 1865–1883. [Google Scholar] [CrossRef] [Green Version]
  23. Artacho, E.; Machado, M.; Sánchez-Portal, D.; Ordejón, P.; Soler, J.M. Electrons in dry DNA from density functional calculations. Mol. Phys. 2003, 101, 1587–1594. [Google Scholar] [CrossRef] [Green Version]
  24. Barbi, M.; Cocco, S.; Peyrard, M. Helicoidal model for DNA opening. Phys. Lett. A 1999, 253, 358–369. [Google Scholar] [CrossRef]
  25. Cocco, S.; Monasson, R. Statistical mechanics of torque induced denaturation of DNA. Phys. Rev. Lett. 1999, 83, 5178–5181. [Google Scholar] [CrossRef] [Green Version]
  26. Agarwal, J.; Hennig, D. Breather solutions of a nonlinear DNA model including a longitudinal degree of freedom. Physica A 2003, 323, 519–533. [Google Scholar] [CrossRef]
  27. Krisch, M.; Mermet, A.; Grimm, H.; Forsyth, V.T.; Rupprecht, A. Phonon dispersion of oriented DNA by inelastic x-ray scattering. Phys. Rev. E 2006, 73, 061909. [Google Scholar] [CrossRef] [Green Version]
  28. Dauxois, T.; Peyrard, M.; Bishop, A.R. Entropy-driven DNA denaturation. Phys. Rev. E 1993, 47, R44–R47. [Google Scholar] [CrossRef]
  29. Zhu, J.X.; Rasmussen, K.O.; Balatsky, A.V.; Bishop, A.R. Local electronic structure in the Peyrard–Bishop–Holstein model. J. Phys. Condens. Matter 2007, 19, 136203. [Google Scholar] [CrossRef]
  30. Cocco, S.; Monasson, R. Theoretical study of collective modes in DNA at ambient temperature. J. Chem. Phys. 2000, 112, 10017–10033. [Google Scholar] [CrossRef] [Green Version]
  31. Campa, A.; Giansanti, A. Experimental tests of the Peyrard-Bishop model applied to the melting of very short DNA chains. Phys. Rev. E 1998, 58, 3585–3588. [Google Scholar] [CrossRef] [Green Version]
  32. Cule, D.; Hwa, T. Denaturation of Heterogeneous DNA. Phys. Rev. Lett. 1997, 79, 2375–2378. [Google Scholar] [CrossRef] [Green Version]
  33. Endres, R.G.; Cox, D.L.; Singh, R.R.P. Colloquium: The quest for high-conductance DNA. Rev. Mod. Phys. 2004, 76, 195–214. [Google Scholar] [CrossRef]
  34. Hawke, L.G.D.; Kalosakas, G.; Simserides, C. Empirical LCAO parameters for π molecular orbitals in planar organic molecules. Mol. Phys. 2009, 107, 1755–1771. [Google Scholar] [CrossRef]
  35. Endres, R.G.; Cox, D.L.; Singh, R.R.P. Electronic properties of DNA: Structural and chemical influence on the quest for high conductance and charge transfer. arXiv 2002, arXiv:cond-mat/0201404. [Google Scholar]
  36. Maciá, E.; Roche, S. Backbone-induced effects in the charge transport efficiency of synthetic DNA molecules. Nanotechnology 2006, 17, 3002–3007. [Google Scholar] [CrossRef]
  37. De Almeida, M.L.; Ourique, G.S.; Fulco, U.L.; Albuquerque, E.L.; de Moura, F.A.B.F.; Lyra, M.L. Charge transport properties of a twisted DNA molecule: A renormalization approach. Chem. Phys. 2016, 478, 48–54. [Google Scholar] [CrossRef]
  38. Hawke, L.G.D.; Kalosakas, G.; Simserides, C. Electronic parameters for charge transfer along DNA. Eur. Phys. J. E 2010, 32, 291–305. [Google Scholar] [CrossRef]
  39. Maciá, E. Electronic structure and transport properties of double-stranded Fibonacci DNA. Phys. Rev. B 2006, 74, 245105. [Google Scholar] [CrossRef] [Green Version]
  40. Cuniberti, G.; Craco, L.; Porath, D.; Dekker, C. Backbone-induced semiconducting behavior in short DNA wires. Phys. Rev. B 2002, 65, 241314. [Google Scholar] [CrossRef] [Green Version]
  41. Cuniberti, G.; Maciá, E.; Rodríguez, A.; Römer, R.A. Charge Migration in DNA: Physics, Chemistry and Biology Perspectives; Chakraborty, T., Ed.; Springer: Berlin, Germany, 2007. [Google Scholar]
  42. Sarmento, R.G.; Fulco, U.L.; Albuquerque, E.L.; Caetano, E.W.S.; Freire, V.N. A renormalization approach to describe charge transport in quasiperiodic dangling backbone ladder (DBL)-DNA molecules. Phys. Lett. A 2011, 375, 3993–3996. [Google Scholar] [CrossRef] [Green Version]
  43. De Moura, F.A.B.F.; Lyra, M.L.; de Almeida, M.L.; Ourique, G.S.; Fulco, U.L.; Albuquerque, E.L. Methylation effect on the ohmic resistance of a poly-GC DNA-like chain. Phys. Lett. A 2016, 380, 3559–3563. [Google Scholar] [CrossRef]
  44. De Almeida, M.L.; Oliveira, J.I.N.; Lima Neto, J.X.; Gomes, C.E.M.; Fulco, U.L.; Albuquerque, E.L.; Freire, V.N.; Caetano, E.W.S.; de Moura, F.A.B.F.; Lyra, M.L. Electronic transport in methylated fragments of DNA. Appl. Phys. Lett. 2015, 107, 203701. [Google Scholar] [CrossRef]
  45. Qi, J.; Govind, N.; Anantram, M.P. The role of cytosine methylation on charge transport through a DNA strand. J. Phys. Lett. 2015, 143, 094306. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Maciá, E. Charge transfer in DNA: Effective Hamiltonian approaches. Z. Kristallogr. 2009, 224, 91–95. [Google Scholar] [CrossRef] [Green Version]
  47. Mantela, M.; Lambropoulos, K.; Theodorakou, M.; Simserides, C. Quasi-Periodic and Fractal Polymers: Energy Structure and Carrier Transfer. Materials 2019, 12, 2177. [Google Scholar] [CrossRef] [Green Version]
  48. Lambropoulos, K.; Simserides, C. Periodic, quasiperiodic, fractal, Kolakoski, and random binary polymers: Energy structure and carrier transport. Phys. Rev. E 2019, 99, 032415. [Google Scholar] [CrossRef] [Green Version]
  49. Lambropoulos, K.; Simserides, C. Tight-Binding Modeling of Nucleic Acid Sequences: Interplay between Various Types of Order or Disorder and Charge Transport. Symmetry 2019, 11, 968. [Google Scholar] [CrossRef] [Green Version]
  50. Lee, O.; Jeon, J.H.; Sung, W. How double-stranded DNA breathing enhances its flexibility and instability on short length scales. Phys. Rev. E 2010, 81, 021906. [Google Scholar] [CrossRef] [Green Version]
  51. Michoel, T.; Van de Peer, Y. Helicoidal transfer matrix model for inhomogeneous DNA melting. Phys. Rev. E 2006, 73, 011908. [Google Scholar] [CrossRef] [Green Version]
  52. Okonogi, T.M.; Alley, S.C.; Harwood, E.A.; Hopkins, P.B.; BHRobinson, B.H. Phosphate backbone neutralization increases duplex DNA flexibility. Proc. Natl. Acad. Sci. USA 2002, 99, 4156–4160. [Google Scholar] [CrossRef] [Green Version]
  53. Duduială, C.I.; Wattis, J.A.D.; Dryden, I.L.; Laughton, C.A. Nonlinear breathing modes at a defect site in DNA. Phys. Rev. E 2009, 80, 061906. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Zdravković, S.; Satarixcx, M.V. Single-molecule unzippering experiments on DNA and Peyrard-Bishop-Dauxois model. Phys. Rev. E 2006, 73, 021905. [Google Scholar] [CrossRef] [PubMed]
  55. Ghorbani, M.; Rafiee, F.M. Geometrical correlations in the nucleosomal DNA conformation and the role of the covalent bonds rigidity. Nucleic Acids Res. 2011, 39, 1220–1230. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Barbi, M.; Lepri, S.; Peyrard, M.; Theodorakopoulos, N. Thermal denaturation of a helicoidal DNA model. Phys. Rev. E 2003, 68, 061909. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Voityuk, A.A. Fluctuation of the electronic coupling in DNA: Multistate versus two-state model. Chem. Phys. Lett. 2007, 439, 162–165. [Google Scholar] [CrossRef]
  58. Mehrez, H.; Anantram, M.P. Interbase electronic coupling for transport through DNA. Phys. Rev. B 2005, 71, 115405. [Google Scholar] [CrossRef] [Green Version]
  59. Peluso, A.; Caruso, T.; Landi, A.; Capobianco, A. The dynamics of hole transfer in DNA. Molecules 2019, 24, 4044. [Google Scholar] [CrossRef] [Green Version]
  60. Leal, M.R.; Weber, G. Sharp DNA denaturation in a helicoidal mesoscopic model. Chem. Phys. Lett. 2020, 755, 137781. [Google Scholar] [CrossRef]
  61. Blancafort, L.; Voityuk, A.A. CASSCF/CAS-PT2 Study of Hole Transfer in Stacked DNA Nucleobases. J. Phys. Chem. A 2006, 110, 6426–6432. [Google Scholar] [CrossRef]
  62. Iguchi, K. π-Electrons in a single strand of DNA: A phenomenological approach. Int. J. Mod. Phys. B 2004, 18, 1845–1910. [Google Scholar] [CrossRef] [Green Version]
  63. Landi, A.; Borrelli, R.; Capobianco, A.; Peluso, A. Transient and enduring electronic resonances drive coherent long distance charge transport in molecular wires. J. Phys. Chem. Lett. 2019, 10, 1845–1851. [Google Scholar] [CrossRef] [PubMed]
  64. Behnia, S.; Fathizadeh, S.; Javanshour, E.; Nemati, F. Light-driven modulation of electrical current through DNA sequences: Engineering of a molecular optical switch. J. Phys. Chem.B 2020, 124, 3261–3270. [Google Scholar] [CrossRef] [PubMed]
  65. González-Jiménez, M.; Ramakrishnan, G.; Harwood, T.; Lapthorn, A.J.; Kelly, S.M.; Ellis, E.M.; Wynne, K. LESS Observation of coherent delocalized phonon-like modes in DNA under physiological conditions. Nat. Commun. 2016, 7, 11799. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Hithell, G.; González-Jiménez, M.; Greetham, G.M.; Donaldson, P.M.; Towrie, M.; Parker, A.W.; Burley, G.A.; Wynne, K.; Hunt, N.T. Ultrafast 2D-IR and optical Kerr effect spectroscopy reveal the impact of duplex melting on the structural dynamics of DNA. Phys. Chem. Chem. Phys. 2017, 19, 10333–10342. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Paez, C.J.; Schulz, P.A. Delocalization of vibrational normal modes in double chains: Application to DNA systems. Microelectron. J. 2008, 39, 1222–1223. [Google Scholar] [CrossRef]
  68. Young, M.A.; Ravishanker, G.; Beveridge, D.L. A 5-nanosecond molecular dynamics trajectory for B-DNA: Analysis of structure, motions, and solvation. Biophys. J. 1997, 73, 2313–2336. [Google Scholar] [CrossRef] [Green Version]
  69. Paciaroni, A.; Comez, L.; Longo, M.; Sebastiani, F.; Bianchi, F.; Orecchini, A.; Zanatta, M.; Verbeni, R.; Bosak, A.; Sacchetti, F.; et al. Terahertz collective dynamics of DNA as affected by hydration and counterions. J. Mol. Liq. 2020, 318, 113956. [Google Scholar] [CrossRef]
  70. Andresen, K.; Das, R.; Park, H.Y.; Smith, H.; Kwok, L.W.; Lamb, J.S.; Kirkl, E.J.; Herschlag, D.; Finkelstein, K.D.; Pollack, L. Spatial Distribution of Competing Ions around DNA in Solution. Phys. Rev. Lett. 2004, 93, 248103. [Google Scholar] [CrossRef]
  71. Iguchi, K. The Huckel theoretical calculation for the electronic structure of DNA. In DNA Engineering: Properties and Applications; Mizoguch, K., Sakamoto, H., Eds.; Pan Stanford: New York, NY, USA, 2016. [Google Scholar]
  72. Lakhno, V.D.; Sultanov, V.B. Band structure of the spectra of Hamiltonians of regular polynucleotide duplexes. Theor. Math. Phys. 2013, 176, 1194–1206. [Google Scholar] [CrossRef]
  73. Yamada, H.; Iguchi, K. Some effective tight-binding models for electrons in DNA conduction: A review. Adv. Cond. Mat. Phys. 2010, 380710. [Google Scholar] [CrossRef] [Green Version]
  74. Hennig, D.; Starikov, E.B.; Archilla, J.F.R.; Palmero, F. Charge Transport in Poly(dG)–Poly(dC) and Poly(dA)–Poly(dT) DNA Polymers. J. Biol. Phys. 2004, 30, 227–238. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Zhang, W.; Govorov, A.O.; Ulloa, S.E. Polarons with a twist. Phys. Rev. B 2002, 66, 060303(R). [Google Scholar] [CrossRef] [Green Version]
  76. Iguchi, K. Semiconductivity and band gap of a double strand of DNA. J. Phys. Soc. Jpn. 2001, 70, 593–597. [Google Scholar] [CrossRef]
  77. Available online: http://telomerase.asu.edu/sequences_telomere.html (accessed on 20 August 2020).
  78. A Tandem Repeat Is a Sequence of Two or More DNA Base Pairs That Is Repeated in Such a Way that the Repeats Lie Adjacent to Each Other on the Chromosome. Tandem Repeats are Generally Associated with Non-coding DNA. In Some Instances, the Number of Times the DNA Sequence Is Repeated Is Variable. Such Variable Tandem Repeats Are Used in DNA Fingerprinting Procedures. Available online: https://www.genome.gov/genetics-glossary/Tandem-Repeat (accessed on 20 August 2020).
  79. Lopes Cardoso, I.; Marques, V. Trinucleotide repeat diseases—Antecipation diseases. J. Clin. Gen. Genom. 2018, 1, 4–9. [Google Scholar]
  80. Beeker, D.; Adhikary, A.; Sevilla, M.D. Charge Migration in DNA. Perspectives from Physics, Chemistry and Biology; Chakraborty, T., Ed.; Springer: Berlin, Germany, 2007; p. 139. [Google Scholar]
  81. Braun, E.; Eichen, Y.; Sivan, U.; Ben-Yoseph, G. DNA-templated assembly and electrode attachment of a conducting silver wire. Nature 1998, 391, 775–779. [Google Scholar] [CrossRef]
  82. Treadway, C.; Hill, M.G.; Barton, J.K. Charge Transport Through a Molecular Pi-stack: Double Helical DNA. Chem. Phys. 2002, 281, 409–428. [Google Scholar] [CrossRef]
  83. Zwang, T.J.; Tse, E.C.M.; Barton, J.K. Sensing DNA through DNA Charge Transport. ACS Chem. Biol. 2018, 13, 1799–1809. [Google Scholar] [CrossRef]
  84. Di Felice, R.; Calzolari, A.; Zhang, H. Towards metalated DNA-based structures. Nanotechnology 2004, 15, 1256–1263. [Google Scholar] [CrossRef]
Figure 1. (a) Diagram of the dsDNA lattice model showing the harmonic bonds between adjacent bps through the sugar-phosphate backbone (zig-zag lines) and between complementary bases due to H-bonds (solid lines). (b) Full-atom dsDNA model with surfaces of constant charge density for the states corresponding to the π -like lowest unoccupied molecular orbital (LUMO, in red) of the C bases and the highest occupied molecular orbital (HOMO, in blue) of the G bases of a polyG-polyC molecule in the dry conditions A form (Courtesy of Emilio Artacho) [23]. (c) Diagram of the dsDNA electronic model showing the π - π channel transfer integrals t n , n ± 1 X Y (cylindrical rods), the interstrand transfer integrals between complementary bases t X Y (perpendicular bars), the glucosidic transfer integrals t P between nucleobases (spheres) and sugar-phosphate groups (cubes), and the on-site energies of these groups γ j .
Figure 1. (a) Diagram of the dsDNA lattice model showing the harmonic bonds between adjacent bps through the sugar-phosphate backbone (zig-zag lines) and between complementary bases due to H-bonds (solid lines). (b) Full-atom dsDNA model with surfaces of constant charge density for the states corresponding to the π -like lowest unoccupied molecular orbital (LUMO, in red) of the C bases and the highest occupied molecular orbital (HOMO, in blue) of the G bases of a polyG-polyC molecule in the dry conditions A form (Courtesy of Emilio Artacho) [23]. (c) Diagram of the dsDNA electronic model showing the π - π channel transfer integrals t n , n ± 1 X Y (cylindrical rods), the interstrand transfer integrals between complementary bases t X Y (perpendicular bars), the glucosidic transfer integrals t P between nucleobases (spheres) and sugar-phosphate groups (cubes), and the on-site energies of these groups γ j .
Materials 13 05119 g001
Figure 2. The p z atomic orbitals overlapping between neighboring nucleobases are described in terms of their Euclidean distance d i j , which is determined by the relative twist and radial variables θ n , n + 1 , r n + 1 , and r n . Reprinted figure with permission from Maciá, E. Physical Review B, 76, 245123, 2007. Copyright (2007) by the American Physical Society.
Figure 2. The p z atomic orbitals overlapping between neighboring nucleobases are described in terms of their Euclidean distance d i j , which is determined by the relative twist and radial variables θ n , n + 1 , r n + 1 , and r n . Reprinted figure with permission from Maciá, E. Physical Review B, 76, 245123, 2007. Copyright (2007) by the American Physical Society.
Materials 13 05119 g002
Figure 3. Sketch illustrating the three-step renormalization process mapping the 3D dsDNA electronic model shown in (a) into the 3D fishbone model shown in (b), then into the effective 1D diatomic and polyatomic lattice models displayed in (c) and (d), respectively, the latter corresponding to a chain made of codon triplet building blocks.
Figure 3. Sketch illustrating the three-step renormalization process mapping the 3D dsDNA electronic model shown in (a) into the 3D fishbone model shown in (b), then into the effective 1D diatomic and polyatomic lattice models displayed in (c) and (d), respectively, the latter corresponding to a chain made of codon triplet building blocks.
Materials 13 05119 g003
Figure 4. Acoustic and optical phonon branches for polyG-polyC and polyA-polyT dsDNA homopolymers (their ν ( q ) and ν + ( q ) curves respectively coincide up to the third decimal place), where q is measured in nm 1 . We used the characteristic and coupled frequency values listed in Table 2.
Figure 4. Acoustic and optical phonon branches for polyG-polyC and polyA-polyT dsDNA homopolymers (their ν ( q ) and ν + ( q ) curves respectively coincide up to the third decimal place), where q is measured in nm 1 . We used the characteristic and coupled frequency values listed in Table 2.
Materials 13 05119 g004
Figure 5. Electronic band structure for a polyG-polyC homopolymer. We used the electronic model parameter values listed in Table 1 and the choice γ = 0 eV. The wavevectors k and q are measured in nm 1 , and q * = π / ξ 2.732 nm 1 . The band structure of a polyA-polyT homopolymer is very similar (see Table 3) with the center of mass of the bands shifted upwards in energy by about 0.02 eV.
Figure 5. Electronic band structure for a polyG-polyC homopolymer. We used the electronic model parameter values listed in Table 1 and the choice γ = 0 eV. The wavevectors k and q are measured in nm 1 , and q * = π / ξ 2.732 nm 1 . The band structure of a polyA-polyT homopolymer is very similar (see Table 3) with the center of mass of the bands shifted upwards in energy by about 0.02 eV.
Materials 13 05119 g005
Table 1. Geometrical, dynamical, lattice model, and electronic model parameters adopted in the dsDNA homopolymer model studied in this work. The same effective Morse potential is used to describe H-bonding in both GC and AT bps. The spring constant k B is difficult to estimate, and different possible values, ranging from 0.04 to 0.5 eV Å 2 , have been reported in the literature [50,51,52,53,54,55].
Table 1. Geometrical, dynamical, lattice model, and electronic model parameters adopted in the dsDNA homopolymer model studied in this work. The same effective Morse potential is used to describe H-bonding in both GC and AT bps. The spring constant k B is difficult to estimate, and different possible values, ranging from 0.04 to 0.5 eV Å 2 , have been reported in the literature [50,51,52,53,54,55].
GeometricalDynamicalLatticeElectronic (eV)Electronic (eV)
θ 0 = π / 5.2 rad M = 653.5 amu b = 0.5 Å 1 [25,56] t G C = 0.01 [38,57,58] ε G = 7.8 8.2 [39,59]
h 0 = 0.34 nm λ G C = 1.130 α = 5 Å 1 [24,56,60] t A T = 0.02 [38] ε A = 8.2 [59]
l 0 = 0.68 nm λ A T = 1.028 D = 0.15 eV [30,56,60] t 0 G G = 0.08 [38,59,61] ε C = 8.9 [59]
R 0 = 1.00 nm μ G C = 4.015 k S = 0.7 eV Å 2 [30] t 0 A A = 0.09 [38,58,59] ε T = 9.0 9.1 [39,59]
ξ = 1.15 nm μ A T = 4.001 k B = 0.04 eV Å 2 [56] t P = 1.5 [62] 0 γ 12 [14,62]
Table 2. Characteristic frequencies and scale times in the lattice dynamics of dsDNA polyG-polyC homopolymers. The  ν φ ρ , ν φ S , and  ν H S values for polyA-polyT homopolymers are slightly smaller by just about 4 GHz.
Table 2. Characteristic frequencies and scale times in the lattice dynamics of dsDNA polyG-polyC homopolymers. The  ν φ ρ , ν φ S , and  ν H S values for polyA-polyT homopolymers are slightly smaller by just about 4 GHz.
Oscillation ω k ( 10 12 rad s 1 ) ν k (THz) τ k (ps)
Twist ω φ 1.07 0.17 5.87
Stacking ω S 3.26 0.52 1.93
H-Bonding ω H 5.34 0.85 1.18
Twist-Stacking ω φ S 6.52 1.04 0.96
Twist-Radial ω φ ρ 12.53 2.00 0.50
Stretch-Stacking ω H S 14.13 2.25 0.45
Table 3. PolyG-polyC and polyA-polyT band structure properties (measured in meV). We adopted the values γ = 0 , a G C = 16.6 eV, a A T = 17.4 eV, b G C = 30.9 , and b A T = 33.6 .
Table 3. PolyG-polyC and polyA-polyT band structure properties (measured in meV). We adopted the values γ = 0 , a G C = 16.6 eV, a A T = 17.4 eV, b G C = 30.9 , and b A T = 33.6 .
Homopolymer W + W E g ( 0 , 0 ) E g ( 0 , q * )
PolyG-PolyC 29.7 24.4 80.3 73.0
PolyA-PolyT 22.0 16.6 93.5 87.4
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Maciá, E. Base-Pairs’ Correlated Oscillation Effects on the Charge Transfer in Double-Helix B-DNA Molecules. Materials 2020, 13, 5119. https://doi.org/10.3390/ma13225119

AMA Style

Maciá E. Base-Pairs’ Correlated Oscillation Effects on the Charge Transfer in Double-Helix B-DNA Molecules. Materials. 2020; 13(22):5119. https://doi.org/10.3390/ma13225119

Chicago/Turabian Style

Maciá, Enrique. 2020. "Base-Pairs’ Correlated Oscillation Effects on the Charge Transfer in Double-Helix B-DNA Molecules" Materials 13, no. 22: 5119. https://doi.org/10.3390/ma13225119

APA Style

Maciá, E. (2020). Base-Pairs’ Correlated Oscillation Effects on the Charge Transfer in Double-Helix B-DNA Molecules. Materials, 13(22), 5119. https://doi.org/10.3390/ma13225119

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop