
A small circular DNA - plasmid pBR322, is considered as a complex dynamic system where nonlinear conformational perturbations which are often named open states or kinks, can arise and propagate. To describe the internal dynamics of the plasmid we use mathematical model consisting of two coupled sine-Gordon equations that in the average field approximation are transformed to two sine-Gordon independent equations with renormalized coefficients. The first equation describes angular oscillations of nitrous bases of the main chain. The second equation describes angular oscillations of nitrous bases in the complementary chain. As a result, two types of kink-like solutions have been obtained. One type kinks were the solutions of the first equation, and the other kinks were the solutions of the second equation. We calculated the main characteristics of the kink motion, including the time dependences of the kink velocity, coordinates, and total energy. These calculations were performed at the initial velocity equal to 1881 m/s which was chosen to avoid reflections from energy barriers corresponding to CDS-1 and CDS-2. The movement of the kinks was investigated by the method of the double energy profile. The maximum complete set of the DNA dynamic parameters was used to calculate the double profile. To calculate the velocity, energy and trajectory of the kinks, the block method was used. The results obtained made it possible to explain in which region of the plasmid the formation of a transcription bubble is most likely, as well to understand in which direction the bubble will move and the transcription process will go.
Citation: Ludmila V. Yakushevich, Larisa A. Krasnobaeva. Double energy profile of pBR322 plasmid[J]. AIMS Biophysics, 2021, 8(2): 221-232. doi: 10.3934/biophy.2021016
[1] | Larisa A. Krasnobaeva, Ludmila V. Yakushevich . DNA kinks behavior in the potential pit-trap. AIMS Biophysics, 2022, 9(2): 130-146. doi: 10.3934/biophy.2022012 |
[2] | Larisa A. Krasnobaeva, Ludmila V. Yakushevich . On the dimensionless model of the transcription bubble dynamics. AIMS Biophysics, 2023, 10(2): 205-219. doi: 10.3934/biophy.2023013 |
[3] | Audrey Lejart, Gilles Salbert, Sébastien Huet . Cytosine hydroxymethylation by TET enzymes: From the control of gene expression to the regulation of DNA repair mechanisms, and back. AIMS Biophysics, 2018, 5(3): 182-193. doi: 10.3934/biophy.2018.3.182 |
[4] | Roman Marks . Bubble mediated polymerization of RNA and DNA. AIMS Biophysics, 2022, 9(2): 96-107. doi: 10.3934/biophy.2022009 |
[5] | Michael J. Hamilton, Matthew D. Young, Silvia Sauer, Ernest Martinez . The interplay of long non-coding RNAs and MYC in cancer. AIMS Biophysics, 2015, 2(4): 794-809. doi: 10.3934/biophy.2015.4.794 |
[6] | Yuri M. Moshkin . Chromatin—a global buffer for eukaryotic gene control. AIMS Biophysics, 2015, 2(4): 531-554. doi: 10.3934/biophy.2015.4.531 |
[7] | Arsha Moorthy, Alireza Sarvestani, Whitney Massock, Chamaree de Silva . Twisting and extension: Application of magnetic tweezers to DNA studies. AIMS Biophysics, 2023, 10(3): 317-346. doi: 10.3934/biophy.2023020 |
[8] | Davood Norouzi, Ataur Katebi, Feng Cui, Victor B. Zhurkin . Topological diversity of chromatin fibers: Interplay between nucleosome repeat length, DNA linking number and the level of transcription. AIMS Biophysics, 2015, 2(4): 613-629. doi: 10.3934/biophy.2015.4.613 |
[9] | Théo Lebeaupin, Hafida Sellou, Gyula Timinszky, Sébastien Huet . Chromatin dynamics at DNA breaks: what, how and why?. AIMS Biophysics, 2015, 2(4): 458-475. doi: 10.3934/biophy.2015.4.458 |
[10] | Sarada Ghosh, Guruprasad Samanta, Manuel De la Sen . Feature selection and classification approaches in gene expression of breast cancer. AIMS Biophysics, 2021, 8(4): 372-384. doi: 10.3934/biophy.2021029 |
A small circular DNA - plasmid pBR322, is considered as a complex dynamic system where nonlinear conformational perturbations which are often named open states or kinks, can arise and propagate. To describe the internal dynamics of the plasmid we use mathematical model consisting of two coupled sine-Gordon equations that in the average field approximation are transformed to two sine-Gordon independent equations with renormalized coefficients. The first equation describes angular oscillations of nitrous bases of the main chain. The second equation describes angular oscillations of nitrous bases in the complementary chain. As a result, two types of kink-like solutions have been obtained. One type kinks were the solutions of the first equation, and the other kinks were the solutions of the second equation. We calculated the main characteristics of the kink motion, including the time dependences of the kink velocity, coordinates, and total energy. These calculations were performed at the initial velocity equal to 1881 m/s which was chosen to avoid reflections from energy barriers corresponding to CDS-1 and CDS-2. The movement of the kinks was investigated by the method of the double energy profile. The maximum complete set of the DNA dynamic parameters was used to calculate the double profile. To calculate the velocity, energy and trajectory of the kinks, the block method was used. The results obtained made it possible to explain in which region of the plasmid the formation of a transcription bubble is most likely, as well to understand in which direction the bubble will move and the transcription process will go.
The construction of energy profiles and their analysis is a new method for studying the physical and functional properties of plasmids. This method came to DNA science from theoretical physics and nonlinear mathematics. In the frameworks of the method, the DNA molecule is considered as a complex dynamic system in which nonlinear conformational perturbations-kinks, can arise and propagate [1]–[4].
Since the work [5], kinks have been actively used to simulate the DNA open states which are small (only a few base pairs long) unwound DNA regions inside which hydrogen bonds between complementary bases are broken. The behavior of the kinks is in many ways similar to that of ordinary particles moving in the DNA potential field. Having calculated the energy profile of the potential field, one can completely solve the problem of the kink dynamics including calculation of time dependence of the kink velocity and construction of the kink trajectory.
A particular example of the open state is the transcription bubble (Figure 1) which is formed at the initial stage of the transcription process as a result of the interaction of RNA polymerase with DNA promoter region [6],[7]. One of the first attempts to calculate the DNA potential field in which the transcriptional bubble moves was made by Salerno [8]. However, the absence at that time of a complete set of dynamic parameters of DNA, as well as the computer power required for such calculations, forced the author to confine himself to only a small fragment of a DNA molecule, which did not allow to obtain reliable results on the kink dynamics and to relate them with functional properties of the molecule.
In this work, we solve this problem on the basis of new data and new approaches to the study of the dynamics of open states in inhomogeneous DNA. In the calculations we use the full set of dynamic parameters of DNA, which was first obtained in [9] and later refined in [10]. The block method and the method of constructing the kinks trajectories are also applied [11],[12]. In contrast to work [12], where only one plasmid profile was calculated, here we use the method of the double profile, which makes it possible to compare the dynamic characteristics of kinks activated in the main and complementary DNA strands. Such a comparison makes it possible to explain in which region of the plasmid the formation of a transcription bubble is most likely, as well to understand in which direction the bubble will move and the transcription process will go.
Plasmid pBR322 [13] is a small circular DNA the sequence of which consists of 4361 bases [14] (Appendix). This sequence contains three coding regions:
CDS-1, encoding tetracycline resistance protein,
CDS-2, encoding ROP protein,
CDS-3, encoding beta-lactamase,
and three promoters: P1, P2 and P3 marked in Appendix by grey. Coordinates of the CDS regions are (86..1276), (1915..2106) and (3293..4153). Coordinates of promoters are (27-33), (43-49) and (4188-4194). The double-stranded scheme of the plasmid is shown in Figure 2a. A linear version of the scheme obtained by cutting at the point Q, is presented in Figure 2b.
Region number | Region coordinates | NA | NT | NG | NC | N |
1 | 1..85 | 22/29 | 29/22 | 17/17 | 17/17 | 85 |
2 (CDS-1) | 86..1276 | 190/268 | 268/190 | 353/381 | 381/353 | 1192 |
3 | 1277..1914 | 142/145 | 145/142 | 160/191 | 191/160 | 638 |
4 (CDS-2) | 1915..2106 | 54/37 | 37/54 | 48/53 | 53/48 | 192 |
5 | 2107..3292 | 294/268 | 268/294 | 319/305 | 305/319 | 1186 |
6 (CDS-3) | 3293..4153 | 213/223 | 223/213 | 203/222 | 222/203 | 861 |
7 | 4154-4361 | 67/64 | 64/67 | 34/42 | 42/34 | 207 |
It is convenient to renumber all regions including CDS regions and the regions between them, starting from the point S. A total of 7 regions turned out. Data on their composition are presented in Table 1.
However, the splitting of the sequence into 7 regions does not reflect the circular nature of the plasmid. To take this into account, we combined the 1st and 7th regions into a single region with number (7 + 1). As a result, the total number of regions decreased by one and became equal to six. In what follows, these six regions will be referred to as blocks.
To describe angular oscillations of nitrous bases in a separate block of the DNA molecule consisting of N bp, let us use the improved Englander model [5], which is a system of 2N coupled nonlinear differential equations [15]:
Here
where
After the averaging procedure (3), equations (1)–(2) are transformed to the form:
In the continuum approximation, these equations take the following form:
In the approximation of average field [17], two coupled equations (6) and (7) are transformed into two independent equations:
where a is the distance between the nearest base pairs. In this approximation, the movement of kinks in the main and complementary chains can be considered separately. The validity of this approximation for DNA was shown in [18], where it was demonstrated numerically that taking into account the interaction between kinks activated in the main and complementary chains leads only to a slight perturbation of the kink shape.
One-soliton solutions of equations (8)–(9), kinks, have then the form :
Here
Equation for the kink velocity was obtained by McLaughlin and Scott [19]. In the case of plasmid pBR322, this equation has the form:
The analytical solution of equation (1) found in [20], has the form :
where
The solution for the kink coordinate
where z0 is the coordinate of the kink at the initial moment of time t0.
The total kink energy is determined by the following formula [21]:
where
In the case of small velocities (
where (
To construct the energy profile, we first calculated the averaged values of the main dynamic parameters of the plasmid:
and also calculated the values of the rest energies of the kinks activated in each of these regions:
Here
Number of the region | Coordinate of the region | ||||
1+7 | (1..85)+(4154..4361) | 6.13/6.28 | 1.94/1.97 | 2.06/2.10 | 1.60/1.63 |
2 (CDS-1) | 86..1276 | 6.05/6.33 | 1.90/1.97 | 2.26/2.32 | 1.66/1.71 |
3 | 1277..1914 | 6.09/6.30 | 1.92/1.96 | 2.21/2.26 | 1.65/1.68 |
4 (CDS-2) | 1915..2106 | 6.26/6.13 | 1.97/1.92 | 2.23/2.19 | 1.67/1.64 |
5 | 2107..3292 | 6.25/6.14 | 1.95/1.93 | 2.22/2.20 | 1.67/1.65 |
6 (CDS-3) | 3293..4153 | 6.14/6.25 | 1.93/1.95 | 2.17/2.20 | 1.64/1.66 |
With the help of the data from the last column of Table 2 we constructed the double energy profile of the potential field of the plasmid (Figure 3). The energy was calculated not at every point, but for each range into which we divided the plasmid. The coordinates of the bases are plotted in Figure 3 along the horizontal in units of base pairs (bp), 1 bp = 3.4×10–10 m. The values of the rest energy of kinks are plotted along the vertical.
Figure 3 shows that the deepest well in the energy profile of both the main (solid line) and complementary (dashed line) sequences corresponds to the (7 + 1)-th region. This means that the probability of kink activation in this region is the highest. This result is in good agreement with the experimental data summarized in [14]. Indeed. the transcription process usually begins with the binding of RNA polymerase to the promoter region of the DNA sequence, and this event is accompanied by local unwinding double-stranded DNA. According to [14] there are three promoters in the (7 + 1)-th region: the P2 promoter located in the main chain, and the P1 and P3 promoters located in the complementary chain. That saturation of the (7 + 1)-th region with promoters indicates that just in this region with high probability the formation of local unwound regions - kinks, takes place.
It can be seen also that in the (7 + 1)-th region, the rest energy of the kinks activated in the main strand is less than the rest energy of the kinks activated in the complementary strand. This means that the probability of the kink activation in the main strand containing P2 promoter is greater than in the complementary strand containing P1 and P3 promoters. This allows us to suggest that the P1 and P3 promoters are much weaker than the P2 promoter.
Using formulas (2)-(4) and the data of Table 2, we plotted the time dependence of the velocity, coordinates and total energy of the plasmid pBR322 kinks (Figures 4–6).
Figure 4 shows that the velocity of the kinks activated in both strands decreases with time, which is explained by the effects of dissipation. The velocity of the kink activated in the complementary strand, decreases faster. The zigzags observed on both curves are explained by abrupt changes in the velocity of the kinks when passing the boundaries between the regions.
Figure 5 shows that the coordinate of the kinks grows smoothly and reaches a constant value, which means that the kink movement stops. Kink activated in the main strand stops at the point 2332 bp, and kink activated in the complementary strand stops at 2558 bp.
From Figure 6 it can be seen that at the same value of the initial kink velocity (1881 m/s), the values of the initial total energy of the kinks activated in different strands differ. For kinks activated in the main strand, the initial total energy is 10.42×10−18 (J), and for kinks activated in the complementary strand, the initial total energy is 12.25×10−18 (J). Later both energy curves decreases and reaches the values of the rest energies.
To construct trajectories of kinks in the 3D space, we solved numerically a system of coupled differential equations for the velocity and coordinates of the kinks:
The results are presented in the form of trajectories of the kink motion in the space {v, z, t} (Figure 7).
In the present work, we investigated the dynamics of the nonlinear conformation distortions – kinks, activated in the plasmid pBR322. The kinks were considered as quasi particles moving in the potential field of the plasmid. The key moment of the method used is the construction of the double energetic profile the analysis of which gives valuable information about the kink behavior.
To imitate the internal dynamics of the plasmid we used mathematical model consisting of two coupled sine-Gordon equations that in the average field approximation were transformed to two sine-Gordon independent equations with renormalized coefficients. The first equation described angular oscillations of nitrous bases of the main chain in the complementary chain. The second equation described angular oscillations of nitrous bases in the complementary chain. As a result, we obtained two types of kink-like solutions: one type are kinks which were the solutions of the first equation, and the other are kinks which were the solutions of the second equation.
To take into account the circular nature of DNA, we combined the end regions into a single block. Using calculated values of the rest energies of the kinks we constructed two energy profiles: the first for the main sequence, the second for the complementary one.
Analysis of these profiles showed that the deepest well in the energy profiles corresponds to the (7 + 1)-th region. This means that the probability of kink activation in this region is the highest.
Because the kink rest energy in the (7 + 1)-th region of the main sequence is less than in that of the complementary sequence we concluded that the P2 promoter located in the (7 + 1)-th region of the main sequence is stronger than the P1 and P3 promoters located in the (7 + 1)-th region of the complementary sequence. Therefore, it can be assumed that the transcription bubble formed in the (7 + 1)-th region is most likely to start moving in the direction of the CDS-1 region.
We calculated the main characteristics of the kink motion, including the time dependences of the kink velocity, coordinates, and total energy. These calculations were performed at the initial velocity equal to 1881 m/s which was chosen to avoid reflections from energy barriers corresponding to CDS-1 and CDS-2. Despite the fact that here we have considered only one value of the initial kink velocity the proposed approach is applicable to any values of the initial kink velocity. It is necessary only to keep in mind one important condition: the initial values should be less than the sound velocity. This condition arises from the mathematical requirements for the existence and stability of the solution in the form of kink.
It should be noted, also, that these results were obtained under three additional constraints. First, we used a simplified model that took into account angular oscillations of the nitrous bases in only one of the two DNA strands. Second, only one type of internal DNA motions-angular oscillations of the nitrous bases, was taken into accountю. Third, the model used did not take into account the DNA torque. In spite of all these limitations, we can expect that the results of the investigations obtained by simple and convenient method correctly convey the basic laws of the kink dynamic behavior in plasmid pBR322 and can be applied to another DNA sequences.
[1] | Yakushevich LV (2006) Nonlinear physics of DNA Weinheim: Wiley. |
[2] | Dauxois T, Peyrard M (2006) Physics of Solitons Cambridge: Cambridge University Press. |
[3] | Shapovalov AV, Krasnobaeva LA (2009) Solitons of the sine-Gordon equation Tomsk: TGU, Available from: http://vital.lib.tsu.ru/vital/access/manager/Repository/vtls:000392302. |
[4] | Ivancevic VG, Ivancevic TT (2013) Sine-Gordon solitons, kinks and breathers as physical models of nonlinear excitations in living cellular structures. J Geom Symmetry Phys 31: 1-56. |
[5] |
Englander SW, Kallenbach NR, Heeger AJ, et al. (1980) Nature of the open state in long polynucleotide double helices: possibility of soliton excitations. Proc Natl Acad Sci USA 77: 7222-7226. doi: 10.1073/pnas.77.12.7222
![]() |
[6] |
Forth S, Sheinin MY, Inman J, et al. (2013) Torque measurement at the single-molecule level. Annu Rev Biophys 42: 583-604. doi: 10.1146/annurev-biophys-083012-130412
![]() |
[7] | Severin ES (2020) Biochemistry Moscow: GEOTAR–Media, Available from: https://expose.gpntbsib.ru/expose/vnp-4f32781c/book/Е2019-3257830172062. |
[8] |
Salerno M (1995) Nonlinear dynamics of plasmid pBR322 promoters. Nonlinear Excitations in Biomolecules New York: Springer, 147-153. doi: 10.1007/978-3-662-08994-1_10
![]() |
[9] | Yakushevich LV, Krasnobaeva LA, Shapovalov AV, et al. (2005) One-and two-soliton solutions of the sine-Gordon equation as applied to DNA. Biophysics 50: 450-455. |
[10] |
Yakushevich LV, Krasnobaeva LA (2016) Forced oscillations of DNA bases. Biophysics 61: 241-250. doi: 10.1134/S000635091602024X
![]() |
[11] |
Grinevich AA, Ryasik AA, Yakushevich LV (2015) Trajectories of DNA bubbles. Chaos, Soliton Fract 75: 62-75. doi: 10.1016/j.chaos.2015.02.009
![]() |
[12] |
Yakushevich LV, Krasnobaeva LA (2019) Plasmid pBR322 and nonlinear conformational distortions (kinks). Math Biol Bioinform 14: 327-339. doi: 10.17537/2019.14.327
![]() |
[13] |
Watson N (1988) A new revision of the sequence of plasmid pBR322. Gene 70: 399-403. doi: 10.1016/0378-1119(88)90212-0
![]() |
[14] | The pBR322 DNA sequence Available from: http://www.ncbi.nlm.nih.gov/nuccore/J01749.1. |
[15] |
Krasnobaeva LA, Yakushevich LV (2015) Rotational dynamics of bases in the gene coding interferon alpha 17 (IFNA17). J Bioinf Comput Biol 13: 1540002. doi: 10.1142/S0219720015400028
![]() |
[16] |
Yakushevich LV, Krasnobaeva LA (2008) A new approach to studies of non-linear dynamics of kinks activated in inhomogeneous polynucleotide chains. Int J Nonlin Mech 43: 1074-1081. doi: 10.1016/j.ijnonlinmec.2008.05.003
![]() |
[17] |
Yakushevich L (2017) On the mechanical analogue of DNA. J Biol Phys 43: 113-125. doi: 10.1007/s10867-016-9437-0
![]() |
[18] |
Yakushevich LV, Savin AV, Manevitch LI (2002) Nonlinear dynamics of topological solitons in DNA. Phys Rev E 66: 016614. doi: 10.1103/PhysRevE.66.016614
![]() |
[19] |
McLaughlin DW, Scott AC (1978) Perturbation analysis of fuxon dynamics. Phys Rev A 18: 1652-1658. doi: 10.1103/PhysRevA.18.1652
![]() |
[20] |
Yakushevich LV, Krasnobaeva LA (2007) Effects of dissipation and external fields on the dynamics of conformational distortions in DNA. Biophysics 52: 237-243. doi: 10.1134/S0006350907020066
![]() |
[21] |
Yakushevich LV, Ryasik AA (2012) Dynamical characteristics of DNA kinks and antikinks. Comput Res Model 4: 209-217. doi: 10.20537/2076-7633-2012-4-1-209-217
![]() |
![]() |
![]() |
1. | L. Krasnobaeva, L. Yakushevich, FEATURES OF THE POTENTIAL FIELD STRUCTURE OF THE PPF1 PLASMID AND THEIR INFLUENCE ON THE CHARACTER OF MOTION OF NONLINEAR CONFORMATIONAL PERTURBATIONS – KINKS, 2022, 7, 2499-9962, 99, 10.29039/rusjbpc.2022.0490 |
Region number | Region coordinates | NA | NT | NG | NC | N |
1 | 1..85 | 22/29 | 29/22 | 17/17 | 17/17 | 85 |
2 (CDS-1) | 86..1276 | 190/268 | 268/190 | 353/381 | 381/353 | 1192 |
3 | 1277..1914 | 142/145 | 145/142 | 160/191 | 191/160 | 638 |
4 (CDS-2) | 1915..2106 | 54/37 | 37/54 | 48/53 | 53/48 | 192 |
5 | 2107..3292 | 294/268 | 268/294 | 319/305 | 305/319 | 1186 |
6 (CDS-3) | 3293..4153 | 213/223 | 223/213 | 203/222 | 222/203 | 861 |
7 | 4154-4361 | 67/64 | 64/67 | 34/42 | 42/34 | 207 |
Number of the region | Coordinate of the region | ||||
1+7 | (1..85)+(4154..4361) | 6.13/6.28 | 1.94/1.97 | 2.06/2.10 | 1.60/1.63 |
2 (CDS-1) | 86..1276 | 6.05/6.33 | 1.90/1.97 | 2.26/2.32 | 1.66/1.71 |
3 | 1277..1914 | 6.09/6.30 | 1.92/1.96 | 2.21/2.26 | 1.65/1.68 |
4 (CDS-2) | 1915..2106 | 6.26/6.13 | 1.97/1.92 | 2.23/2.19 | 1.67/1.64 |
5 | 2107..3292 | 6.25/6.14 | 1.95/1.93 | 2.22/2.20 | 1.67/1.65 |
6 (CDS-3) | 3293..4153 | 6.14/6.25 | 1.93/1.95 | 2.17/2.20 | 1.64/1.66 |
Region number | Region coordinates | NA | NT | NG | NC | N |
1 | 1..85 | 22/29 | 29/22 | 17/17 | 17/17 | 85 |
2 (CDS-1) | 86..1276 | 190/268 | 268/190 | 353/381 | 381/353 | 1192 |
3 | 1277..1914 | 142/145 | 145/142 | 160/191 | 191/160 | 638 |
4 (CDS-2) | 1915..2106 | 54/37 | 37/54 | 48/53 | 53/48 | 192 |
5 | 2107..3292 | 294/268 | 268/294 | 319/305 | 305/319 | 1186 |
6 (CDS-3) | 3293..4153 | 213/223 | 223/213 | 203/222 | 222/203 | 861 |
7 | 4154-4361 | 67/64 | 64/67 | 34/42 | 42/34 | 207 |
Number of the region | Coordinate of the region | ||||
1+7 | (1..85)+(4154..4361) | 6.13/6.28 | 1.94/1.97 | 2.06/2.10 | 1.60/1.63 |
2 (CDS-1) | 86..1276 | 6.05/6.33 | 1.90/1.97 | 2.26/2.32 | 1.66/1.71 |
3 | 1277..1914 | 6.09/6.30 | 1.92/1.96 | 2.21/2.26 | 1.65/1.68 |
4 (CDS-2) | 1915..2106 | 6.26/6.13 | 1.97/1.92 | 2.23/2.19 | 1.67/1.64 |
5 | 2107..3292 | 6.25/6.14 | 1.95/1.93 | 2.22/2.20 | 1.67/1.65 |
6 (CDS-3) | 3293..4153 | 6.14/6.25 | 1.93/1.95 | 2.17/2.20 | 1.64/1.66 |