1
|
Muniz MI, Bustos AH, Slott S, Astakhova K, Weber G. Cation valence dependence of hydrogen bond and stacking potentials in DNA mesoscopic models. Biophys Chem 2023; 294:106949. [PMID: 36706510 DOI: 10.1016/j.bpc.2022.106949] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2022] [Revised: 12/23/2022] [Accepted: 12/25/2022] [Indexed: 12/30/2022]
Abstract
Monovalent and divalent cations play a crucial role in living cells and for molecular techniques such as PCR. Here we evaluate DNA melting temperatures in magnesium (Mg2+) and magnesium‑potassium (Mg2++ K+) buffers with a mesoscopic model that allows us to estimate hydrogen bonds and stacking interaction potentials. The Mg2+ and Mg2++ K+ results are compared to previous calculations for sodium ions (Na+), in terms of equivalent sodium concentration and ionic strength. Morse potentials, related to hydrogen bonding, were found to be essentially constant and unaffected by cation conditions. However, for stacking interactions we find a clear dependence with ionic strength and cation valence. The highest ionic strength variations, for both hydrogen bonds and stacking interactions, was found at the sequence terminals. This suggests that end-to-end interactions in DNA will be strongly dependent on cation valence and ionic strength.
Collapse
Affiliation(s)
- Maria Izabel Muniz
- Departamento de Física, Universidade Federal de Minas Gerais, 31270-901 Belo Horizonte, MG, Brazil
| | - Adrian H Bustos
- Department of Chemistry, Technical University of Denmark, 206-207 Kemitorvet, 2800 Kongens Lyngby, Denmark
| | - Sofie Slott
- Department of Chemistry, Technical University of Denmark, 206-207 Kemitorvet, 2800 Kongens Lyngby, Denmark
| | - Kira Astakhova
- Department of Chemistry, Technical University of Denmark, 206-207 Kemitorvet, 2800 Kongens Lyngby, Denmark
| | - Gerald Weber
- Departamento de Física, Universidade Federal de Minas Gerais, 31270-901 Belo Horizonte, MG, Brazil.
| |
Collapse
|
2
|
Mesoscopic model confirms strong base pair metal mediated bonding for T-Hg 2+-T and weaker for C-Ag +-C. Chem Phys Lett 2022. [DOI: 10.1016/j.cplett.2022.139847] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
3
|
Ferreira I, Slott S, Astakhova K, Weber G. Complete Mesoscopic Parameterization of Single LNA Modifications in DNA Applied to Oncogene Probe Design. J Chem Inf Model 2021; 61:3615-3624. [PMID: 34251211 DOI: 10.1021/acs.jcim.1c00470] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The use of mesoscopic models to describe the thermodynamic properties of locked nucleic acid (LNA)-modified nucleotides can provide useful insights into their properties, such as hydrogen-bonding and stacking interactions. In addition, the mesoscopic parameters can be used to optimize LNA insertion in probes, to achieve accurate melting temperature predictions, and to obtain duplex opening profiles at the base-pair level. Here, we applied this type of model to parameterize a large set of melting temperatures for LNA-modified sequences, from published sources, covering all possible nearest-neighbor configurations. We have found a very large increase in Morse potentials, which indicates very strong hydrogen bonding as the main cause of improved LNA thermodynamic stability. LNA-modified adenine-thymine (AT) was found to have similar hydrogen bonding to unmodified cytosine-guanine (CG) base pairs, while for LNA CG, we found exceptionally large hydrogen bonding. In contrast, stacking interactions, which were thought to be behind the stability of LNA, were similar to unmodified DNA in most cases. We applied the new LNA parameters to the design of BRAF, KRAS, and EGFR oncogene variants by testing all possible LNA modifications. Selected sequences were then synthesized and had their hybridization temperatures measured, achieving a prediction accuracy within 1 °C. We performed a detailed base-pair opening analysis to discuss specific aspects of these probe hybridizations that may be relevant for probe design.
Collapse
Affiliation(s)
- Izabela Ferreira
- Departamento de Física, Universidade Federal de Minas Gerais, 31270-901 Belo Horizonte, MG, Brazil.,Programa Interunidades de Pós-Graduação em Bioinformática, Universidade Federal de Minas Gerais, 31270-901 Belo Horizonte, MG, Brazil
| | - Sofie Slott
- Department of Chemistry, Technical University of Denmark, Kemitorvet, Bygning 207, 2800 Kgs. Lyngby, Denmark
| | - Kira Astakhova
- Department of Chemistry, Technical University of Denmark, Kemitorvet, Bygning 207, 2800 Kgs. Lyngby, Denmark
| | - Gerald Weber
- Departamento de Física, Universidade Federal de Minas Gerais, 31270-901 Belo Horizonte, MG, Brazil
| |
Collapse
|
4
|
|
5
|
Optical and theoretical study of strand recognition by nucleic acid probes. Commun Chem 2020; 3:111. [PMID: 36703315 PMCID: PMC9814704 DOI: 10.1038/s42004-020-00362-5] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2020] [Accepted: 07/23/2020] [Indexed: 01/29/2023] Open
Abstract
Detection of nucleic acids is crucial to the study of their basic properties and consequently to applying this knowledge to the determination of pathologies such as cancer. In this work, our goal is to determine new trends for creating diagnostic tools for cancer driver mutations. Herein, we study a library of natural and modified oligonucleotide duplexes by a combination of optical and theoretical methods. We report a profound effect of additives on the duplexes, including nucleic acids as an active crowder. Unpredictably and inconsistent with DNA+LNA/RNA duplexes, locked nucleic acids contribute poorly to mismatch discrimination in the DNA+LNA/DNA duplexes. We develop a theoretical framework that explains poor mismatch discrimination in KRAS oncogene. We implement our findings in a bead-bait genotyping assay to detect mutated human cancer RNA. The performance of rationally designed probes in this assay is superior to the LNA-primer polymerase chain reaction, and it agrees with sequencing data.
Collapse
|
6
|
Oliveira LM, Long AS, Brown T, Fox KR, Weber G. Melting temperature measurement and mesoscopic evaluation of single, double and triple DNA mismatches. Chem Sci 2020; 11:8273-8287. [PMID: 34094181 PMCID: PMC8163305 DOI: 10.1039/d0sc01700k] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023] Open
Abstract
Unlike the canonical base pairs AT and GC, the molecular properties of mismatches such as hydrogen bonding and stacking interactions are strongly dependent on the identity of the neighbouring base pairs. As a result, due to the sheer number of possible combinations of mismatches and flanking base pairs, only a fraction of these have been studied in varying experiments or theoretical models. Here, we report on the melting temperature measurement and mesoscopic analysis of contiguous DNA mismatches in nearest-neighbours and next-nearest neighbour contexts. A total of 4032 different mismatch combinations, including single, double and triple mismatches were covered. These were compared with 64 sequences containing all combinations of canonical base pairs in the same location under the same conditions. For a substantial number of single mismatch configurations, 15%, the measured melting temperatures were higher than the least stable AT base pair. The mesoscopic calculation, using the Peyrard-Bishop model, was performed on the set of 4096 sequences, and resulted in estimates of on-site and nearest-neighbour interactions that can be correlated to hydrogen bonding and base stacking. Our results confirm many of the known properties of mismatches, including the peculiar sheared stacking of tandem GA mismatches. More intriguingly, it also reveals that a number of mismatches present strong hydrogen bonding when flanked on both sites by other mismatches. To highlight the applicability of our results, we discuss a number of practical situations such as enzyme binding affinities, thymine DNA glycosylase repair activity, and trinucleotide repeat expansions.
Collapse
Affiliation(s)
- Luciana M Oliveira
- Departamento de Física, Universidade Federal de Minas Gerais 31270-901 Belo Horizonte MG Brazil +55 31 3409 5600 +55 31 3409 6616
| | - Adam S Long
- School of Biological Sciences, University of Southampton Life Sciences Building 85 Southampton SO17 1BJ UK
| | - Tom Brown
- Department of Chemistry, University of Oxford Oxford UK
| | - Keith R Fox
- School of Biological Sciences, University of Southampton Life Sciences Building 85 Southampton SO17 1BJ UK
| | - Gerald Weber
- Departamento de Física, Universidade Federal de Minas Gerais 31270-901 Belo Horizonte MG Brazil +55 31 3409 5600 +55 31 3409 6616
| |
Collapse
|
7
|
Muniz MI, Lackey HH, Heemstra JM, Weber G. DNA/TNA mesoscopic modeling of melting temperatures suggests weaker hydrogen bonding of CG than in DNA/RNA. Chem Phys Lett 2020. [DOI: 10.1016/j.cplett.2020.137413] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
|
8
|
de Oliveira Martins E, Basílio Barbosa V, Weber G. DNA/RNA hybrid mesoscopic model shows strong stability dependence with deoxypyrimidine content and stacking interactions similar to RNA/RNA. Chem Phys Lett 2019. [DOI: 10.1016/j.cplett.2018.11.015] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
|
9
|
Mesoscopic modelling of Cy3 and Cy5 dyes attached to DNA duplexes. Biophys Chem 2017; 230:62-67. [DOI: 10.1016/j.bpc.2017.08.007] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2017] [Revised: 08/21/2017] [Accepted: 08/27/2017] [Indexed: 11/19/2022]
|
10
|
de Oliveira Martins E, Weber G. An asymmetric mesoscopic model for single bulges in RNA. J Chem Phys 2017; 147:155102. [PMID: 29055303 DOI: 10.1063/1.5006948] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
Simple one-dimensional DNA or RNA mesoscopic models are of interest for their computational efficiency while retaining the key elements of the molecular interactions. However, they only deal with perfectly formed DNA or RNA double helices and consider the intra-strand interactions to be the same on both strands. This makes it difficult to describe highly asymmetric structures such as bulges and loops and, for instance, prevents the application of mesoscopic models to determine RNA secondary structures. Here we derived the conditions for the Peyrard-Bishop mesoscopic model to overcome these limitations and applied it to the calculation of single bulges, the smallest and simplest of these asymmetric structures. We found that these theoretical conditions can indeed be applied to any situation where stacking asymmetry needs to be considered. The full set of parameters for group I RNA bulges was determined from experimental melting temperatures using an optimization procedure, and we also calculated average opening profiles for several RNA sequences. We found that guanosine bulges show the strongest perturbation on their neighboring base pairs, considerably reducing the on-site interactions of their neighboring base pairs.
Collapse
Affiliation(s)
- Erik de Oliveira Martins
- Departamento de Física, Universidade Federal de Minas Gerais, 31270-901 Belo Horizonte, MG, Brazil
| | - Gerald Weber
- Departamento de Física, Universidade Federal de Minas Gerais, 31270-901 Belo Horizonte, MG, Brazil
| |
Collapse
|
11
|
Amarante TD, Weber G. Evaluating Hydrogen Bonds and Base Stacking of Single, Tandem and Terminal GU Mismatches in RNA with a Mesoscopic Model. J Chem Inf Model 2015; 56:101-9. [DOI: 10.1021/acs.jcim.5b00571] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Tauanne D. Amarante
- Departamento de Física, Universidade Federal de Minas Gerais, 31270-901 Belo
Horizonte-MG, Brazil
| | - Gerald Weber
- Departamento de Física, Universidade Federal de Minas Gerais, 31270-901 Belo
Horizonte-MG, Brazil
| |
Collapse
|
12
|
Ferreira I, Amarante TD, Weber G. DNA terminal base pairs have weaker hydrogen bonds especially for AT under low salt concentration. J Chem Phys 2015; 143:175101. [DOI: 10.1063/1.4934783] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
Affiliation(s)
- Izabela Ferreira
- Departamento de Física, Universidade Federal de Minas Gerais, 31270-901 Belo Horizonte-MG, Brazil
| | - Tauanne D. Amarante
- Departamento de Física, Universidade Federal de Minas Gerais, 31270-901 Belo Horizonte-MG, Brazil
| | - Gerald Weber
- Departamento de Física, Universidade Federal de Minas Gerais, 31270-901 Belo Horizonte-MG, Brazil
| |
Collapse
|
13
|
Maximiano RV, Weber G. Deoxyinosine mismatch parameters calculated with a mesoscopic model result in uniform hydrogen bonding and strongly variable stacking interactions. Chem Phys Lett 2015. [DOI: 10.1016/j.cplett.2015.04.045] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
14
|
Weber G. Optimization method for obtaining nearest-neighbour DNA entropies and enthalpies directly from melting temperatures. Bioinformatics 2014; 31:871-7. [DOI: 10.1093/bioinformatics/btu751] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|