1
|
Sengupta S, Bhattacharjee SM, Mishra G. Large bubble drives circular DNA melting. Phys Chem Chem Phys 2024. [PMID: 39027987 DOI: 10.1039/d4cp01536c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/20/2024]
Abstract
We investigate the melting transition of non-supercoiled circular DNA of different lengths, employing Brownian dynamics simulations. In the absence of supercoiling, we find that melting of circular DNA is driven by a large bubble, which agrees with the previous predictions of circular DNA melting in the presence of supercoiling. By analyzing sector-wise changes in average base-pair distance, our study reveals that the melting behavior of circular DNA closely resembles that of linear DNA. Additionally, we find a marked difference in the thermal stability of circular DNA over linear DNA at very short length scales, an effect that diminishes as the length of circular DNA increases. The stability of smaller circular DNA is linked to the occurrence of transient small bubbles, characterized by a lower probability of growth.
Collapse
Affiliation(s)
- Souradeep Sengupta
- Department of Physics, Ashoka University, Sonipat, Haryana - 131029, India.
| | | | - Garima Mishra
- Department of Physics, Ashoka University, Sonipat, Haryana - 131029, India.
| |
Collapse
|
2
|
Zoli M. Twist-stretch relations in nucleic acids. EUROPEAN BIOPHYSICS JOURNAL : EBJ 2023; 52:641-650. [PMID: 37357224 DOI: 10.1007/s00249-023-01669-6] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/04/2023] [Revised: 05/15/2023] [Accepted: 06/10/2023] [Indexed: 06/27/2023]
Abstract
Nucleic acids are highly deformable helical molecules constantly stretched, twisted and bent in their biological functioning. Single molecule experiments have shown that double stranded (ds)-RNA and standard ds-DNA have opposite twist-stretch patterns and stretching properties when overwound under a constant applied load. The key structural features of the A-form RNA and B-form DNA helices are here incorporated in a three-dimensional mesoscopic Hamiltonian model which accounts for the radial, bending and twisting fluctuations of the base pairs. Using path integral techniques which sum over the ensemble of the base pair fluctuations, I compute the average helical repeat of the molecules as a function of the load. The obtained twist-stretch relations and stretching properties, for short A- and B-helical fragments, are consistent with the opposite behaviors observed in kilo-base long molecules.
Collapse
Affiliation(s)
- Marco Zoli
- School of Science and Technology, University of Camerino, 62032, Camerino, Italy.
| |
Collapse
|
3
|
Kabir A, Bhattarai M, Rasmussen KØ, Shehu A, Usheva A, Bishop AR, Alexandrov B. Examining DNA breathing with pyDNA-EPBD. Bioinformatics 2023; 39:btad699. [PMID: 37991847 PMCID: PMC10681863 DOI: 10.1093/bioinformatics/btad699] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2023] [Revised: 10/23/2023] [Accepted: 11/21/2023] [Indexed: 11/24/2023] Open
Abstract
MOTIVATION The two strands of the DNA double helix locally and spontaneously separate and recombine in living cells due to the inherent thermal DNA motion. This dynamics results in transient openings in the double helix and is referred to as "DNA breathing" or "DNA bubbles." The propensity to form local transient openings is important in a wide range of biological processes, such as transcription, replication, and transcription factors binding. However, the modeling and computer simulation of these phenomena, have remained a challenge due to the complex interplay of numerous factors, such as, temperature, salt content, DNA sequence, hydrogen bonding, base stacking, and others. RESULTS We present pyDNA-EPBD, a parallel software implementation of the Extended Peyrard-Bishop-Dauxois (EPBD) nonlinear DNA model that allows us to describe some features of DNA dynamics in detail. The pyDNA-EPBD generates genomic scale profiles of average base-pair openings, base flipping probability, DNA bubble probability, and calculations of the characteristically dynamic length indicating the number of base pairs statistically significantly affected by a single point mutation using the Markov Chain Monte Carlo algorithm. AVAILABILITY AND IMPLEMENTATION pyDNA-EPBD is supported across most operating systems and is freely available at https://github.com/lanl/pyDNA_EPBD. Extensive documentation can be found at https://lanl.github.io/pyDNA_EPBD/.
Collapse
Affiliation(s)
- Anowarul Kabir
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87544, United States
- Department of Computer Science, George Mason University, Fairfax, VA 22030, United States
| | - Manish Bhattarai
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87544, United States
| | - Kim Ø Rasmussen
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87544, United States
| | - Amarda Shehu
- Department of Computer Science, George Mason University, Fairfax, VA 22030, United States
| | - Anny Usheva
- Department of Surgery, Brown University, Providence, RI 02912, United States
| | - Alan R Bishop
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87544, United States
| | - Boian Alexandrov
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87544, United States
| |
Collapse
|
4
|
Genthon A, Dvirnas A, Ambjörnsson T. Equilibrium melting probabilities of a DNA molecule with a defect: An exact solution of the Poland-Scheraga model. J Chem Phys 2023; 159:145102. [PMID: 37815110 DOI: 10.1063/5.0168915] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2023] [Accepted: 09/20/2023] [Indexed: 10/11/2023] Open
Abstract
In this study we derive analytically the equilibrium melting probabilities for basepairs of a DNA molecule with a defect site. We assume that the defect is characterized by a change in the Watson-Crick basepair energy of the defect basepair, and in the associated two stacking energies for the defect, as compared to the remaining parts of the DNA. The defect site could, for instance, occur due to DNA basepair mismatching, cross-linking, or by the chemical modifications when attaching fluorescent labels, such as fluorescent-quencher pairs, to DNA. Our exact solution of the Poland-Scheraga model for DNA melting provides the probability that the labeled basepair, and its neighbors, are open at different temperatures. Our work is of direct importance, for instance, for studies where fluorophore-quencher pairs are used for studying single basepair fluctuations of designed DNA molecules.
Collapse
Affiliation(s)
- Arthur Genthon
- Max Planck Institute for the Physics of Complex Systems, 01187 Dresden, Germany
| | - Albertas Dvirnas
- Computational Biology and Biological Physics, Centre for Environmental and Climate Science, Lund University, SE-223 62 Lund, Sweden
| | - Tobias Ambjörnsson
- Computational Biology and Biological Physics, Centre for Environmental and Climate Science, Lund University, SE-223 62 Lund, Sweden
| |
Collapse
|
5
|
Basov A, Dorohova A, Malyshko V, Moiseev A, Svidlov A, Bezhenar M, Nechipurenko Y, Dzhimak S. Influence of a Single Deuterium Substitution for Protium on the Frequency Generation of Different-Size Bubbles in IFNA17. Int J Mol Sci 2023; 24:12137. [PMID: 37569512 PMCID: PMC10418495 DOI: 10.3390/ijms241512137] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2023] [Revised: 07/21/2023] [Accepted: 07/23/2023] [Indexed: 08/13/2023] Open
Abstract
The influence of a single 2H/1H replacement on the frequency generation of different-size bubbles in the human interferon alpha-17 gene (IFNA17) under various energies was studied by a developed algorithm and mathematical modeling without simplifications or averaging. This new approach showed the efficacy of researching DNA bubbles and open states both when all hydrogen bonds in nitrogenous base pairs are protium and after an 2H-substitution. After a single deuterium substitution under specific energies, it was demonstrated that the non-coding region of IFNA17 had a more significant regulatory role in bubble generation in the whole gene than the promoter had. It was revealed that a single deuterium substitution for protium has an influence on the frequency generation of DNA bubbles, which also depends on their size and is always higher for the smaller bubbles under the largest number of the studied energies. Wherein, compared to the natural condition under the same critical value of energy, the bigger raises of the bubble frequency occurrence (maximums) were found for 11-30 base pair (bp) bubbles (higher by 319%), 2-4 bp bubbles (higher by 300%), and 31 bp and over ones (higher by 220%); whereas the most significant reductions of the indicators (minimums) were observed for 11-30 bp bubbles (lower by 43%) and bubbles size over 30 bp (lower by 82%). In this study, we also analyzed the impact of several circumstances on the AT/GC ratio in the formation of DNA bubbles, both under natural conditions and after a single hydrogen isotope exchange. Moreover, based on the obtained data, substantial positive and inverse correlations were revealed between the AT/GC ratio and some factors (energy values, size of DNA bubbles). So, this modeling and variant of the modified algorithm, adapted for researching DNA bubbles, can be useful to study the regulation of replication and transcription in the genes under different isotopic substitutions in the nucleobases.
Collapse
Affiliation(s)
- Alexandr Basov
- Department of Fundamental and Clinical Biochemistry, Kuban State Medical University, Krasnodar 350063, Russia; (A.B.); (V.M.)
- Department of Radiophysics and Nanotechnology, Kuban State University, Krasnodar 350040, Russia; (A.D.); (A.S.); (S.D.)
| | - Anna Dorohova
- Department of Radiophysics and Nanotechnology, Kuban State University, Krasnodar 350040, Russia; (A.D.); (A.S.); (S.D.)
- Laboratory of Problems of Stable Isotope Spreading in Living Systems, Federal Research Center of the Southern Scientific Center of the Russian Academy of Sciences, Rostov-on-Don 344006, Russia
| | - Vadim Malyshko
- Department of Fundamental and Clinical Biochemistry, Kuban State Medical University, Krasnodar 350063, Russia; (A.B.); (V.M.)
- Laboratory of Problems of Stable Isotope Spreading in Living Systems, Federal Research Center of the Southern Scientific Center of the Russian Academy of Sciences, Rostov-on-Don 344006, Russia
| | - Arkadii Moiseev
- Scientific Department, Kuban State Agrarian University, Krasnodar 350004, Russia;
| | - Alexandr Svidlov
- Department of Radiophysics and Nanotechnology, Kuban State University, Krasnodar 350040, Russia; (A.D.); (A.S.); (S.D.)
- Laboratory of Problems of Stable Isotope Spreading in Living Systems, Federal Research Center of the Southern Scientific Center of the Russian Academy of Sciences, Rostov-on-Don 344006, Russia
| | - Maria Bezhenar
- Department of Function Theory, Kuban State University, Krasnodar 350040, Russia;
| | - Yury Nechipurenko
- Engelhardt Institute of Molecular Biology, Russian Academy of Sciences, Moscow 119991, Russia
| | - Stepan Dzhimak
- Department of Radiophysics and Nanotechnology, Kuban State University, Krasnodar 350040, Russia; (A.D.); (A.S.); (S.D.)
- Laboratory of Problems of Stable Isotope Spreading in Living Systems, Federal Research Center of the Southern Scientific Center of the Russian Academy of Sciences, Rostov-on-Don 344006, Russia
| |
Collapse
|
6
|
Bubble Relaxation Dynamics in Homopolymer DNA Sequences. Molecules 2023; 28:molecules28031041. [PMID: 36770707 PMCID: PMC9920605 DOI: 10.3390/molecules28031041] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2022] [Revised: 01/11/2023] [Accepted: 01/15/2023] [Indexed: 01/22/2023] Open
Abstract
Understanding the inherent timescales of large bubbles in DNA is critical to a thorough comprehension of its physicochemical characteristics, as well as their potential role on helix opening and biological function. In this work, we employ the coarse-grained Peyrard-Bishop-Dauxois model of DNA to study relaxation dynamics of large bubbles in homopolymer DNA, using simulations up to the microsecond time scale. By studying energy autocorrelation functions of relatively large bubbles inserted into thermalised DNA molecules, we extract characteristic relaxation times from the equilibration process for both adenine-thymine (AT) and guanine-cytosine (GC) homopolymers. Bubbles of different amplitudes and widths are investigated through extensive statistics and appropriate fittings of their relaxation. Characteristic relaxation times increase with bubble amplitude and width. We show that, within the model, relaxation times are two orders of magnitude longer in GC sequences than in AT sequences. Overall, our results confirm that large bubbles leave a lasting impact on the molecule's dynamics, for times between 0.5-500 ns depending on the homopolymer type and bubble shape, thus clearly affecting long-time evolutions of the molecule.
Collapse
|
7
|
Zoli M. Non-linear Hamiltonian models for DNA. EUROPEAN BIOPHYSICS JOURNAL : EBJ 2022; 51:431-447. [PMID: 35976412 DOI: 10.1007/s00249-022-01614-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/16/2022] [Revised: 07/23/2022] [Accepted: 07/25/2022] [Indexed: 06/15/2023]
Abstract
Nucleic acids' physical properties have been investigated by theoretical methods based both on fully atomistic representations and on coarse-grained models, e.g., the worm-like-chain, taken from polymer physics. In this review article, I discuss an intermediate (mesoscopic) approach and show how to build a three-dimensional Hamiltonian model which accounts for the main interactions responsible for the stability of the helical molecules. While the 3D mesoscopic model yields a sufficiently detailed description of the helix at the level of the base pair, it also allows one to predict the thermodynamical and structural properties of molecules in solution. Relying on the idea that the base pair fluctuations can be conceived as trajectories, I have built over the past years a computational method based on the time-dependent path integral formalism to derive the partition function. While the main features of the method are presented, I focus here in particular on a newly developed statistical method to set the maximum amplitude of the base pair fluctuations, a key parameter of the theory. Some applications to the calculation of DNA flexibility properties are discussed together with the available experimental data.
Collapse
Affiliation(s)
- Marco Zoli
- School of Science and Technology, University of Camerino, 62032, Camerino, Italy.
| |
Collapse
|
8
|
Abstract
We study dsDNA (double strand DNA) melting in detail within varying strip-like confinement in a two-dimensional lattice model. The interplay between reduced configurational entropy and attractive base-pairing energy results in a non-monotonic melting profile of DNA. Structural transitions associated with confined DNA melting reveal a stretched or extended state for very strong confinement. By using the exact enumeration method, we investigate the emergence of a local denatured zone e.g. bubbles during DNA melting. The survival time of a single bubble within varying strip width is studied from the Fokker-Planck formalism by considering the bubble size as a reaction co-ordinate. We show that a simple lattice model can capture the sequence heterogeneity effect on DNA melting and bubble dynamics within the strip. Different time scales of bubble zipping for different DNA sequences are found, which may have potential applications in denaturation mapping.
Collapse
Affiliation(s)
- Dibyajyoti Mohanta
- Department of Physics, IIT (BHU), Varanasi 221005, India.
- The Institute of Mathematical Sciences, C.I.T Campus, Taramani, Chennai 600113, India
- Homi Bhaba National Institute, Training School Complex, Anushakti Nagar, Mumbai 400094, India
| |
Collapse
|
9
|
Hillebrand M, Kalosakas G, Bishop AR, Skokos C. Bubble lifetimes in DNA gene promoters and their mutations affecting transcription. J Chem Phys 2021; 155:095101. [PMID: 34496591 DOI: 10.1063/5.0060335] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023] Open
Abstract
Relative lifetimes of inherent double stranded DNA openings with lengths up to ten base pairs are presented for different gene promoters and corresponding mutants that either increase or decrease transcriptional activity in the framework of the Peyrard-Bishop-Dauxois model. Extensive microcanonical simulations are used with energies corresponding to physiological temperature. The bubble lifetime profiles along the DNA sequences demonstrate a significant reduction of the average lifetime at the mutation sites when the mutated promoter decreases transcription, while a corresponding enhancement of the bubble lifetime is observed in the case of mutations leading to increased transcription. The relative difference in bubble lifetimes between the mutated and wild type promoters at the position of mutation varies from 20% to more than 30% as the bubble length decreases.
Collapse
Affiliation(s)
- M Hillebrand
- Nonlinear Dynamics and Chaos Group, Department of Mathematics and Applied Mathematics, University of Cape Town, Rondebosch 7701, South Africa
| | - G Kalosakas
- Department of Materials Science, University of Patras, GR-26504 Rio, Greece
| | - A R Bishop
- Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Ch Skokos
- Nonlinear Dynamics and Chaos Group, Department of Mathematics and Applied Mathematics, University of Cape Town, Rondebosch 7701, South Africa
| |
Collapse
|
10
|
Abstract
A statistical method is developed to estimate the maximum amplitude of the base pair fluctuations in a three dimensional mesoscopic model for nucleic acids. The base pair thermal vibrations around the helix diameter are viewed as a Brownian motion for a particle embedded in a stable helical structure. The probability to return to the initial position is computed, as a function of time, by integrating over the particle paths consistent with the physical properties of the model potential. The zero time condition for the first-passage probability defines the constraint to select the integral cutoff for various macroscopic helical conformations, obtained by tuning the twist, bending, and slide motion between adjacent base pairs along the molecule stack. Applying the method to a short homogeneous chain at room temperature, we obtain meaningful estimates for the maximum fluctuations in the twist conformation with ∼10.5 base pairs per helix turn, typical of double stranded DNA helices. Untwisting the double helix, the base pair fluctuations broaden and the integral cutoff increases. The cutoff is found to increase also in the presence of a sliding motion, which shortens the helix contour length, a situation peculiar of dsRNA molecules.
Collapse
Affiliation(s)
- Marco Zoli
- School of Science and Technology, University of Camerino, I-62032 Camerino, Italy
| |
Collapse
|