1
|
Mu ZC, Tan YL, Zhang BG, Liu J, Shi YZ. Ab initio predictions for 3D structure and stability of single- and double-stranded DNAs in ion solutions. PLoS Comput Biol 2022; 18:e1010501. [PMID: 36260618 PMCID: PMC9621594 DOI: 10.1371/journal.pcbi.1010501] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2022] [Revised: 10/31/2022] [Accepted: 09/27/2022] [Indexed: 11/07/2022] Open
Abstract
The three-dimensional (3D) structure and stability of DNA are essential to understand/control their biological functions and aid the development of novel materials. In this work, we present a coarse-grained (CG) model for DNA based on the RNA CG model proposed by us, to predict 3D structures and stability for both dsDNA and ssDNA from the sequence. Combined with a Monte Carlo simulated annealing algorithm and CG force fields involving the sequence-dependent base-pairing/stacking interactions and an implicit electrostatic potential, the present model successfully folds 20 dsDNAs (≤52nt) and 20 ssDNAs (≤74nt) into the corresponding native-like structures just from their sequences, with an overall mean RMSD of 3.4Å from the experimental structures. For DNAs with various lengths and sequences, the present model can make reliable predictions on stability, e.g., for 27 dsDNAs with/without bulge/internal loops and 24 ssDNAs including pseudoknot, the mean deviation of predicted melting temperatures from the corresponding experimental data is only ~2.0°C. Furthermore, the model also quantificationally predicts the effects of monovalent or divalent ions on the structure stability of ssDNAs/dsDNAs. To determine 3D structures and quantify stability of single- (ss) and double-stranded (ds) DNAs is essential to unveil the mechanisms of their functions and to further guide the production and development of novel materials. Although many DNA models have been proposed to reproduce the basic structural, mechanical, or thermodynamic properties of dsDNAs based on the secondary structure information or preset constraints, there are very few models can be used to investigate the ssDNA folding or dsDNA assembly from the sequence. Furthermore, due to the polyanionic nature of DNAs, metal ions (e.g., Na+ and Mg2+) in solutions can play an essential role in DNA folding and dynamics. Nevertheless, ab initio predictions for DNA folding in ion solutions are still an unresolved problem. In this work, we developed a novel coarse-grained model to predict 3D structures and thermodynamic stabilities for both ssDNAs and dsDNAs in monovalent/divalent ion solutions from their sequences. As compared with the extensive experimental data and available existing models, we showed that the present model can successfully fold simple DNAs into their native-like structures, and can also accurately reproduce the effects of sequence and monovalent/divalent ions on structure stability for ssDNAs including pseudoknot and dsDNAs with/without bulge/internal loops.
Collapse
Affiliation(s)
- Zi-Chun Mu
- Research Center of Nonlinear Science, School of Mathematical & Physical Sciences, Wuhan Textile University, Wuhan, China
- School of Computer Science and Artificial Intelligence, Wuhan Textile University, Wuhan, China
| | - Ya-Lan Tan
- Research Center of Nonlinear Science, School of Mathematical & Physical Sciences, Wuhan Textile University, Wuhan, China
| | - Ben-Gong Zhang
- Research Center of Nonlinear Science, School of Mathematical & Physical Sciences, Wuhan Textile University, Wuhan, China
| | - Jie Liu
- Research Center of Nonlinear Science, School of Mathematical & Physical Sciences, Wuhan Textile University, Wuhan, China
| | - Ya-Zhou Shi
- Research Center of Nonlinear Science, School of Mathematical & Physical Sciences, Wuhan Textile University, Wuhan, China
- * E-mail:
| |
Collapse
|
2
|
Satpathi S, Endoh T, Sugimoto N. Applicability of the nearest-neighbour model for pseudoknot RNAs. Chem Commun (Camb) 2022; 58:5952-5955. [PMID: 35451430 DOI: 10.1039/d1cc07094k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The validity of the nearest-neighbour (NN) model was verified in an RNA pseudoknot (PK) structure. The thermodynamic parameters of the second hairpin stem (S2) region, which separates the PK from a hairpin structure, were monitored using CD and UV melting. Different PKs with identical NN base pairs in the S2 region exhibited similar thermodynamic parameters, highlighting the validity of the NN model in this RNA tertiary structure motif.
Collapse
Affiliation(s)
- Sagar Satpathi
- Frontier Institute for Biomolecular Engineering Research (FIBER), Konan University, 7-1-20 Minatojima-minamimachi, Kobe, 650-0047, Japan
| | - Tamaki Endoh
- Frontier Institute for Biomolecular Engineering Research (FIBER), Konan University, 7-1-20 Minatojima-minamimachi, Kobe, 650-0047, Japan
| | - Naoki Sugimoto
- Frontier Institute for Biomolecular Engineering Research (FIBER), Konan University, 7-1-20 Minatojima-minamimachi, Kobe, 650-0047, Japan.,Graduate School of Frontiers of Innovative Research in Science and Technology (FIRST), Konan University, 7-1-20 Minatojima-minamimachi, Kobe, 650-0047, Japan.
| |
Collapse
|
3
|
Malousi A, Andreou AZ, Georgiou E, Tzimagiorgis G, Kovatsi L, Kouidou S. Age-dependent methylation in epigenetic clock CpGs is associated with G-quadruplex, co-transcriptionally formed RNA structures and tentative splice sites. Epigenetics 2018; 13:808-821. [PMID: 30270726 PMCID: PMC6224212 DOI: 10.1080/15592294.2018.1514232] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
Horvath's epigenetic clock consists of 353 CpGs whose methylation levels can accurately predict the age of individuals. Using bioinformatics analysis, we investigated the conformation, energy characteristics and presence of tentative splice sites of the sequences surrounding the epigenetic clock CpGs, in relation to the median methylation changes in different ages, the presence of CpG islands and their position in genes. Common characteristics in the 100 nt sequences surrounding the epigenetic clock CpGs are G-quadruplexes and/or tentative splice site motifs. Median methylation increases significantly in sequences which adopt less stable structures during transcription. Methylation is higher when CpGs overlap with G-quadruplexes than when they precede them. Median methylation in epigenetic clock CpGs is higher in sequences expressed as single products rather than in multiple products and those containing single donors and multiple acceptors. Age-related methylation variation is significant in sequences without G-quadruplexes, particularly those producing low stability nascent RNA and those with splice sites. CpGs in sequences close to transcription start sites and those which are possibly never expressed (hypothetical proteins) undergo similar extent of age-related median methylation decrease and increase. Preservation of methylation is observed in CpG islands without G-quadruplexes, contrary to CpGs far from CpG islands (open sea). Sequences containing G-quadruplexes and RNA pseudoknots, determining the recognition by H3K27 histone methyltransferase, are hypomethylated. The presented structural DNA and co-transcriptional RNA analysis of epigenetic clock sequences, foreshadows the association of age-related methylation changes with the principle biological processes of DNA and histone methylation, splicing and chromatin silencing.
Collapse
Affiliation(s)
- Andigoni Malousi
- a Laboratory of Biological Chemistry , Medical School, Aristotle University of Thessaloniki , Thessaloniki , Greece
| | | | - Elisavet Georgiou
- a Laboratory of Biological Chemistry , Medical School, Aristotle University of Thessaloniki , Thessaloniki , Greece
| | - Georgios Tzimagiorgis
- a Laboratory of Biological Chemistry , Medical School, Aristotle University of Thessaloniki , Thessaloniki , Greece
| | - Leda Kovatsi
- c Laboratory of Forensic Medicine & Toxicology , Medical School, Aristotle University of Thessaloniki , Thessaloniki , Greece
| | - Sofia Kouidou
- a Laboratory of Biological Chemistry , Medical School, Aristotle University of Thessaloniki , Thessaloniki , Greece
| |
Collapse
|
4
|
Reiling-Steffensmeier C, Marky LA. Structural Insight into the Unbound State of the DNA Analogue of the PreQ 1 Riboswitch: A Thermodynamic Approach. Biochemistry 2017; 56:6231-6239. [PMID: 29076719 DOI: 10.1021/acs.biochem.7b00596] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The preQ1 riboswitch aptamer domain is very dynamic in its unbound state with the ability to form multiple structures: a hairpin, kissing hairpins, and pseudoknot-like structure. The aim of this study is to determine whether the DNA analogue (PreQ1) is able to form structures similar to that of the reported RNA aptamer. Using a thermodynamic approach, we report on structural determination using differential scanning calorimetry under different salt conditions. Further analysis of the primary sequence allowed us to design modified molecules to determine what potential structures are forming in this single-stranded DNA analogue. We found, in a 16 mM Na+ solution, PreQ1 has three transitions with TM values of 14.8, 19.4, and 26.2 °C and a total ΔH of -44.7 kcal/mol. With the increase in salt concentration to 116 mM, there are TM values of 22.3, 28.7, and 38.9 °C and a ΔH of -69.1 kcal/mol, while at 216 mM, the three transitions have TM values of 24.4, 31.6, and 42.9 °C with a total ΔH of -71.5 kcal/mol. Therefore, the increase in enthalpy is due to the formation of additional base-pair stacks. The modified molecules, which would inhibit pseudoknot formation, kissing hairpins, and internal loop interactions, were fully characterized and compared to the native DNA analogue. The analysis of the enthalpy and differential binding of counterions allows us to conclude this single-stranded DNA analogue under physiological conditions is not forming a pseudoknot-like structure. Instead, two potential structures, Compact-Hairpin and Kissing-Complex, are more likely and could be in equilibrium.
Collapse
Affiliation(s)
- Calliste Reiling-Steffensmeier
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Nebraska Medical Center , 986025 Nebraska Medical Center, Omaha, Nebraska 68198-6025, United States
| | - Luis A Marky
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Nebraska Medical Center , 986025 Nebraska Medical Center, Omaha, Nebraska 68198-6025, United States
| |
Collapse
|
5
|
Lee HT, Carr CE, Khutsishvili I, Marky LA. Effect of Loop Length and Sequence on the Stability of DNA Pyrimidine Triplexes with TAT Base Triplets. J Phys Chem B 2017; 121:9175-9184. [PMID: 28875701 DOI: 10.1021/acs.jpcb.7b07591] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
We report the thermodynamic contributions of loop length and loop sequence to the overall stability of DNA intramolecular pyrimidine triplexes. Two sets of triplexes were designed: in the first set, the C5 loop closing the triplex stem was replaced with 5'-CTnC loops (n = 1-5), whereas in the second set, both the duplex and triplex loops were replaced with a 5'-GCAA or 5'-AACG tetraloop. For the triplexes with a 5'-CTnC loop, the triplex with five bases in the loop has the highest stability relative to the control. A loop length lower than five compromises the strength of the base-pair stacks without decreasing the thermal stability, leading to a decreased enthalpy, whereas an increase in the loop length leads to a decreased enthalpy and a higher entropic penalty. The incorporation of the GCAA loop yielded more stable triplexes, whereas the incorporation of AACG in the triplex loop yielded a less stable triplex due to an unfavorable enthalpy term. Thus, addition of the GCAA tetraloop can cause an increase in the thermodynamics of the triplex without affecting the sequence or melting behavior and may result in an additional layer of genetic regulation.
Collapse
Affiliation(s)
- Hui-Ting Lee
- Department of Pharmaceutical Sciences, University of Nebraska Medical Center , 986025 Nebraska Medical Center, Omaha, Nebraska 68198-6025, United States
| | - Carolyn E Carr
- Department of Pharmaceutical Sciences, University of Nebraska Medical Center , 986025 Nebraska Medical Center, Omaha, Nebraska 68198-6025, United States
| | - Irine Khutsishvili
- Department of Pharmaceutical Sciences, University of Nebraska Medical Center , 986025 Nebraska Medical Center, Omaha, Nebraska 68198-6025, United States
| | - Luis A Marky
- Department of Pharmaceutical Sciences, University of Nebraska Medical Center , 986025 Nebraska Medical Center, Omaha, Nebraska 68198-6025, United States
| |
Collapse
|
6
|
Jacobson DR, Saleh OA. Counting the ions surrounding nucleic acids. Nucleic Acids Res 2017; 45:1596-1605. [PMID: 28034959 PMCID: PMC5389524 DOI: 10.1093/nar/gkw1305] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2016] [Accepted: 12/21/2016] [Indexed: 01/29/2023] Open
Abstract
Nucleic acids are strongly negatively charged, and thus electrostatic interactions—screened by ions in solution—play an important role in governing their ability to fold and participate in biomolecular interactions. The negative charge creates a region, known as the ion atmosphere, in which cation and anion concentrations are perturbed from their bulk values. Ion counting experiments quantify the ion atmosphere by measuring the preferential ion interaction coefficient: the net total number of excess ions above, or below, the number expected due to the bulk concentration. The results of such studies provide important constraints on theories, which typically predict the full three-dimensional distribution of the screening cloud. This article reviews the state of nucleic acid ion counting measurements and critically analyzes their ability to test both analytical and simulation-based models.
Collapse
Affiliation(s)
- David R Jacobson
- Department of Physics, University of California, Santa Barbara, CA 93106, USA
| | - Omar A Saleh
- Materials Department and BMSE Program, University of California, Santa Barbara, CA 93106, USA
| |
Collapse
|
7
|
Q Nguyen KK, Gomez YK, Bakhom M, Radcliffe A, La P, Rochelle D, Lee JW, Sorin EJ. Ensemble simulations: folding, unfolding and misfolding of a high-efficiency frameshifting RNA pseudoknot. Nucleic Acids Res 2017; 45:4893-4904. [PMID: 28115636 PMCID: PMC5416846 DOI: 10.1093/nar/gkx012] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2016] [Accepted: 01/11/2017] [Indexed: 12/11/2022] Open
Abstract
Massive all-atom molecular dynamics simulations were conducted across a distributed computing network to study the folding, unfolding, misfolding and conformational plasticity of the high-efficiency frameshifting double mutant of the 26 nt potato leaf roll virus RNA pseudoknot. Our robust sampling, which included over 40 starting structures spanning the spectrum from the extended unfolded state to the native fold, yielded nearly 120 μs of cumulative sampling time. Conformational microstate transitions on the 1.0 ns to 10.0 μs timescales were observed, with post-equilibration sampling providing detailed representations of the conformational free energy landscape and the complex folding mechanism inherent to the pseudoknot motif. Herein, we identify and characterize two alternative native structures, three intermediate states, and numerous misfolded states, the latter of which have not previously been characterized via atomistic simulation techniques. While in line with previous thermodynamics-based models of a general RNA folding mechanism, our observations indicate that stem-strand-sequence-separation may serve as an alternative predictor of the order of stem formation during pseudoknot folding. Our results contradict a model of frameshifting based on structural rigidity and resistance to mechanical unfolding, and instead strongly support more recent studies in which conformational plasticity is identified as a determining factor in frameshifting efficiency.
Collapse
Affiliation(s)
- Khai K Q Nguyen
- Department of Chemistry & Biochemistry, California State University Long Beach, Long Beach, CA 90840, USA.,Department of Computer Engineering & Computer Science, California State University Long Beach, Long Beach, CA 90840, USA
| | - Yessica K Gomez
- Department of Chemistry & Biochemistry, California State University Long Beach, Long Beach, CA 90840, USA.,Department of Physics & Astronomy, California State University Long Beach, Long Beach, CA 90840, USA
| | - Mona Bakhom
- Department of Chemistry & Biochemistry, California State University Long Beach, Long Beach, CA 90840, USA
| | - Amethyst Radcliffe
- Department of Physics & Astronomy, California State University Long Beach, Long Beach, CA 90840, USA
| | - Phuc La
- Department of Chemistry & Biochemistry, California State University Long Beach, Long Beach, CA 90840, USA
| | - Dakota Rochelle
- Department of Chemistry & Biochemistry, California State University Long Beach, Long Beach, CA 90840, USA
| | - Ji Won Lee
- Department of Chemistry & Biochemistry, California State University Long Beach, Long Beach, CA 90840, USA
| | - Eric J Sorin
- Department of Chemistry & Biochemistry, California State University Long Beach, Long Beach, CA 90840, USA
| |
Collapse
|
8
|
Jacobson DR, Saleh OA. Quantifying the ion atmosphere of unfolded, single-stranded nucleic acids using equilibrium dialysis and single-molecule methods. Nucleic Acids Res 2016; 44:3763-71. [PMID: 27036864 PMCID: PMC4856996 DOI: 10.1093/nar/gkw196] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2016] [Accepted: 03/14/2016] [Indexed: 12/31/2022] Open
Abstract
To form secondary structure, nucleic acids (NAs) must overcome electrostatic strand–strand repulsion, which is moderated by the surrounding atmosphere of screening ions. The free energy of NA folding therefore depends on the interactions of this ion atmosphere with both the folded and unfolded states. We quantify such interactions using the preferential ion interaction coefficient or ion excess: the number of ions present near the NA in excess of the bulk concentration. The ion excess of the folded, double-helical state has been extensively studied; however, much less is known about the salt-dependent ion excess of the unfolded, single-stranded state. We measure this quantity using three complementary approaches: a direct approach of Donnan equilibrium dialysis read out by atomic emission spectroscopy and two indirect approaches involving either single-molecule force spectroscopy or existing thermal denaturation data. The results of these three approaches, each involving an independent experimental technique, are in good agreement. Even though the single-stranded NAs are flexible polymers that are expected to adopt random-coil configurations, we find that their ion atmosphere is quantitatively described by rod-like models that neglect large-scale conformational freedom, an effect that we explain in terms of the competition between the relevant structural and electrostatic length scales.
Collapse
Affiliation(s)
- David R Jacobson
- Department of Physics, University of California, Santa Barbara, CA 93106, USA
| | - Omar A Saleh
- Materials Department and BMSE Program, University of California, Santa Barbara, CA, USA
| |
Collapse
|
9
|
The Complementarity of the Loop to the Stem in DNA Pseudoknots Gives Rise to Local TAT Base-Triplets. Methods Enzymol 2016. [PMID: 26794363 DOI: 10.1016/bs.mie.2015.07.029] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/15/2023]
Abstract
Pseudoknots belong to an RNA structural motif that has significant roles in the biological function of RNA. An example is ribosomal frameshifting; in this mechanism, the formation of a local triplex changes the reading frame that allows for differences in the translation of mRNAs. In this work, we have used a combination of temperature-dependent UV spectroscopy and differential scanning calorimetry (DSC) to determine the unfolding thermodynamics of a set of DNA pseudoknots with the following sequence: d(TCTCTTnAAAAAAAAGAGAT5TTTTTTT), where "Tn" is a thymine loop with n=5 (PsK-5), 7 (PsK-7), 9 (PsK-9), or 11 (PsK-11). All four oligonucleotides form intramolecular pseudoknots, and the increase in the length of this loop yielded more stable pseudoknots due to higher transition temperatures and higher unfolding enthalpies. This indicates formation of one and three TAT/TAT stacks in PsK-9 and PsK-11, respectively. We have flipped one AT for a TA base pair in the core stem of these pseudoknots, preventing in this way the formation of these base-triplet stacks. The DSC curves of these pseudoknots yielded lower unfolding enthalpies, confirming the formation of a local triplex in PsK-9 and PsK-11. Furthermore, we have investigated the reaction of PsK-5 and PsK-9 with their partially complementary strands: directly by isothermal titration calorimetry and indirectly by creating a Hess cycle with the DSC data. Relative to the PsK-5 reaction, PsK-9 reacts with its complementary strand with less favorable free energy and enthalpy contributions; this indicates PsK-9 is more stable and more compact due to the formation of a local triplex.
Collapse
|
10
|
Application of differential scanning calorimetry to measure the differential binding of ions, water and protons in the unfolding of DNA molecules. Biochim Biophys Acta Gen Subj 2015; 1860:990-998. [PMID: 26450631 DOI: 10.1016/j.bbagen.2015.10.002] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2015] [Revised: 09/28/2015] [Accepted: 10/01/2015] [Indexed: 11/24/2022]
Abstract
BACKGROUND The overall stability of DNA molecules globally depends on base-pair stacking, base-pairing, polyelectrolyte effect and hydration contributions. In order to understand how they carry out their biological roles, it is essential to have a complete physical description of how the folding of nucleic acids takes place, including their ion and water binding. SCOPE OF REVIEW To investigate the role of ions, water and protons in the stability and melting behavior of DNA structures, we report here an experimental approach i.e., mainly differential scanning calorimetry (DSC), to determine linking numbers: the differential binding of ions (Δnion), water (ΔnW) and protons (ΔnH(+)) in the helix-coil transition of DNA molecules. GENERAL SIGNIFICANCE We use DSC and temperature-dependent UV spectroscopic techniques to measure the differential binding of ions, water, and protons for the unfolding of a variety of DNA molecules: salmon testes DNA (ST-DNA), one dodecamer, one undecamer and one decamer duplexes, nine hairpin loops, and two triplexes. These methods can be applied to any conformational transition of a biomolecule. MAJOR CONCLUSIONS We determined complete thermodynamic profiles, including all three linking numbers, for the unfolding of each molecule. The favorable folding of a DNA helix results from a favorable enthalpy-unfavorable entropy compensation. DSC thermograms and UV melts as a function of salt, osmolyte and proton concentrations yielded releases of ions and water. Therefore, the favorable folding of each DNA molecule results from the formation of base-pair stacks and uptake of both counterions and water molecules. In addition, the triplex with C(+)GC base triplets yielded an uptake of protons. Furthermore, the folding of a DNA duplex is accompanied by a lower uptake of ions and a similar uptake of four water molecules as the DNA helix gets shorter. In addition, the oligomer duplexes and hairpin thermodynamic data suggest ion and water binding depends on the DNA sequence rather than DNA composition.
Collapse
|