1
|
Arteaga S, Dolenz BJ, Znosko BM. Competitive Influence of Alkali Metals in the Ion Atmosphere on Nucleic Acid Duplex Stability. ACS OMEGA 2024; 9:1287-1297. [PMID: 38222622 PMCID: PMC10785066 DOI: 10.1021/acsomega.3c07563] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/29/2023] [Revised: 12/04/2023] [Accepted: 12/11/2023] [Indexed: 01/16/2024]
Abstract
The nonspecific atmosphere around nucleic acids, often termed the ion atmosphere, encompasses a collection of weak ion-nucleic acid interactions. Although nonspecific, the ion atmosphere has been shown to influence nucleic acid folding and structural stability. Studies investigating the composition of the ion atmosphere have shown competitive occupancy of the atmosphere between metal ions in the same solution. Many studies have investigated single ion effects on nucleic acid secondary structure stability; however, no comprehensive studies have investigated how the competitive occupancy of mixed ions in the ion atmosphere influences nucleic acid secondary structure stability. Here, six oligonucleotides were optically melted in buffers containing molar quantities, or mixtures, of either XCl (X = Li, K, Rb, or Cs) or NaCl. A correction factor was developed to better predict RNA duplex stability in solutions containing mixed XCl/NaCl. For solutions containing a 1:1 mixture of XCl/NaCl, one alkali metal chloride contributed more to duplex stability than the other. Overall, there was a 54% improvement in predictive capabilities with the correction factor compared with the standard 1.0 M NaCl nearest-neighbor models. This correction factor can be used in models to better predict RNA secondary structure in solutions containing mixed XCl/NaCl.
Collapse
Affiliation(s)
- Sebastian
J. Arteaga
- Department of Chemistry, Saint Louis University, Saint
Louis, Missouri 63103, United States
| | - Bruce J. Dolenz
- Department of Chemistry, Saint Louis University, Saint
Louis, Missouri 63103, United States
| | - Brent M. Znosko
- Department of Chemistry, Saint Louis University, Saint
Louis, Missouri 63103, United States
| |
Collapse
|
2
|
Templeton C, Hamilton I, Russell R, Elber R. Impact of Ion-Mixing Entropy on Orientational Preferences of DNA Helices: FRET Measurements and Computer Simulations. J Phys Chem B 2023; 127:8796-8808. [PMID: 37815452 DOI: 10.1021/acs.jpcb.3c04354] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/11/2023]
Abstract
Biological processes require DNA and RNA helices to pack together in specific interhelical orientations. While electrostatic repulsion between backbone charges is expected to be maximized when helices are in parallel alignment, such orientations are commonplace in nature. To better understand how the repulsion is overcome, we used experimental and computational approaches to investigate how the orientational preferences of DNA helices depend on the concentration and valence of mobile cations. We used Förster resonance energy transfer (FRET) to probe the relative orientations of two 24-bp helices held together via a freely rotating PEG linker. At low cation concentrations, the helices preferred more "cross"-like orientations over those closer to parallel, and this preference was reduced with increasing salt concentrations. The results were in good quantitative agreement with Poisson-Boltzmann (PB) calculations for monovalent salt (Na+). However, PB underestimated the ability of mixtures of monovalent and divalent ions (Mg2+) to reduce the conformational preference. As a complementary approach, we performed all-atom molecular dynamics (MD) simulations and found better agreement with the experimental results. While MD and PB predict similar electrostatic forces, MD predicts a greater accumulation of Mg2+ in the ion atmosphere surrounding the DNA. Mg2+ occupancy is predicted to be greater in conformations close to the parallel orientation than in conformations close to the crossed orientation, enabling a greater release of Na+ ions and providing an entropic gain (one bound ion for two released). MD predicts an entropy gain larger than that of PB because of the increased Mg2+ occupancy. The entropy changes have a negligible effect at low Mg2+ concentrations because the free energies are dominated by electrostatic repulsion. However, as the Mg2+ concentration increases, charge screening is more effective and the mixing entropy produces readily detectable changes in packing preferences. Our results underline the importance of mixing entropy of counterions in nucleic acid interactions and provide a new understanding on the impact of a mixed ion atmosphere on the packing of DNA helices.
Collapse
Affiliation(s)
- Clark Templeton
- Department of Chemical Engineering, University of Texas at Austin, Austin, Texas 78712, United States
- Department of Physics, FU Berlin, 14195 Berlin, Germany
| | - Ian Hamilton
- Department of Molecular Biosciences, University of Texas at Austin, Austin, Texas 78712, United States
| | - Rick Russell
- Department of Molecular Biosciences, University of Texas at Austin, Austin, Texas 78712, United States
| | - Ron Elber
- Institute for Computational Engineering and Science, Department of Chemistry, University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
3
|
Kolesnikov ES, Gushchin IY, Zhilyaev PA, Onufriev AV. Why Na+ has higher propensity than K+ to condense DNA in a crowded environment. J Chem Phys 2023; 159:145103. [PMID: 37815107 DOI: 10.1063/5.0159341] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2023] [Accepted: 08/22/2023] [Indexed: 10/11/2023] Open
Abstract
Experimentally, in the presence of the crowding agent polyethylene glycol (PEG), sodium ions compact double-stranded DNA more readily than potassium ions. Here, we have used molecular dynamics simulations and the "ion binding shells model" of DNA condensation to provide an explanation for the observed variations in condensation of short DNA duplexes in solutions containing different monovalent cations and PEG; several predictions are made. According to the model we use, externally bound ions contribute the most to the ion-induced aggregation of DNA duplexes. The simulations reveal that for two adjacent DNA duplexes, the number of externally bound Na+ ions is larger than the number of K+ ions over a wide range of chloride concentrations in the presence of PEG, providing a qualitative explanation for the higher propensity of sodium ions to compact DNA under crowded conditions. The qualitative picture is confirmed by an estimate of the corresponding free energy of DNA aggregation that is at least 0.2kBT per base pair more favorable in solution with NaCl than with KCl at the same ion concentration. The estimated attraction free energy of DNA duplexes in the presence of Na+ depends noticeably on the DNA sequence; we predict that AT-rich DNA duplexes are more readily condensed than GC-rich ones in the presence of Na+. Counter-intuitively, the addition of a small amount of a crowding agent with high affinity for the specific condensing ion may lead to the weakening of the ion-mediated DNA-DNA attraction, shifting the equilibrium away from the DNA condensed phase.
Collapse
Affiliation(s)
- Egor S Kolesnikov
- Research Center for Molecular Mechanisms of Aging and Age-Related Diseases, Moscow Institute of Physics and Technology, 141700 Dolgoprudny, Russia
- Department of Physics, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Ivan Yu Gushchin
- Research Center for Molecular Mechanisms of Aging and Age-Related Diseases, Moscow Institute of Physics and Technology, 141700 Dolgoprudny, Russia
| | - Petr A Zhilyaev
- The Center for Materials Technologies, Skolkovo Institute of Science and Technology, Bolshoy Boulevard 30, bld. 1, Moscow 121205, Russia
| | - Alexey V Onufriev
- Department of Physics, Virginia Tech, Blacksburg, Virginia 24061, USA
- Department of Computer Science, Virginia Tech, 2160C Torgersen Hall, Blacksburg, Virginia 24061, USA
- Center for Soft Matter and Biological Physics, Virginia Tech, Blacksburg, Virginia 24061, USA
| |
Collapse
|
4
|
Kührová P, Mlýnský V, Otyepka M, Šponer J, Banáš P. Sensitivity of the RNA Structure to Ion Conditions as Probed by Molecular Dynamics Simulations of Common Canonical RNA Duplexes. J Chem Inf Model 2023; 63:2133-2146. [PMID: 36989143 PMCID: PMC10091408 DOI: 10.1021/acs.jcim.2c01438] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2022] [Indexed: 03/30/2023]
Abstract
RNA molecules play a key role in countless biochemical processes. RNA interactions, which are of highly diverse nature, are determined by the fact that RNA is a highly negatively charged polyelectrolyte, which leads to intimate interactions with an ion atmosphere. Although RNA molecules are formally single-stranded, canonical (Watson-Crick) duplexes are key components of folded RNAs. A double-stranded (ds) RNA is also important for the design of RNA-based nanostructures and assemblies. Despite the fact that the description of canonical dsRNA is considered the least problematic part of RNA modeling, the imperfect shape and flexibility of dsRNA can lead to imbalances in the simulations of larger RNAs and RNA-containing assemblies. We present a comprehensive set of molecular dynamics (MD) simulations of four canonical A-RNA duplexes. Our focus was directed toward the characterization of the influence of varying ion concentrations and of the size of the solvation box. We compared several water models and four RNA force fields. The simulations showed that the A-RNA shape was most sensitive to the RNA force field, with some force fields leading to a reduced inclination of the A-RNA duplexes. The ions and water models played a minor role. The effect of the box size was negligible, and even boxes with a small fraction of the bulk solvent outside the RNA hydration sphere were sufficient for the simulation of the dsRNA.
Collapse
Affiliation(s)
- Petra Kührová
- Regional
Centre of Advanced Technologies and Materials, Czech Advanced Technology
and Research Institute (CATRIN), Palacký
University Olomouc, Šlechtitelů 27, 779 00 Olomouc, Czech Republic
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 00 Brno, Czech Republic
| | - Vojtěch Mlýnský
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 00 Brno, Czech Republic
| | - Michal Otyepka
- Regional
Centre of Advanced Technologies and Materials, Czech Advanced Technology
and Research Institute (CATRIN), Palacký
University Olomouc, Šlechtitelů 27, 779 00 Olomouc, Czech Republic
- IT4Innovations, VSB − Technical University of Ostrava, 17. listopadu 2172/15, 708 00 Ostrava, Poruba, Czech Republic
| | - Jiří Šponer
- Regional
Centre of Advanced Technologies and Materials, Czech Advanced Technology
and Research Institute (CATRIN), Palacký
University Olomouc, Šlechtitelů 27, 779 00 Olomouc, Czech Republic
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 00 Brno, Czech Republic
| | - Pavel Banáš
- Regional
Centre of Advanced Technologies and Materials, Czech Advanced Technology
and Research Institute (CATRIN), Palacký
University Olomouc, Šlechtitelů 27, 779 00 Olomouc, Czech Republic
| |
Collapse
|
5
|
Kim SE, Hong SC. Two Opposing Effects of Monovalent Cations on the Stability of i-Motif Structure. J Phys Chem B 2023; 127:1932-1939. [PMID: 36811958 DOI: 10.1021/acs.jpcb.2c07069] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/24/2023]
Abstract
At acidic pH, cytosine-rich single-stranded DNA can be folded into a tetraplex structure called i-motif (iM). In recent studies, the effect of monovalent cations on the stability of iM structure has been addressed, but a consensus about the issue has not been reached yet. Thus, we investigated the effects of various factors on the stability of iM structure using fluorescence resonance energy transfer (FRET)-based analysis for three types of iM derived from human telomere sequences. We confirmed that the protonated cytosine-cytosine (C:C+) base pair is destabilized as the concentration of monovalent cations (Li+, Na+, K+) increases and that Li+ has the greatest tendency of destabilization. Intriguingly, monovalent cations would play an ambivalent role in iM formation by making single-stranded DNA flexible and pliant for an iM structure. In particular, we found that Li+ has a notably greater flexibilizing effect than Na+ and K+. All taken together, we conclude that the stability of iM structure is controlled by the subtle balance of the two counteractive effects of monovalent cations: electrostatic screening and disruption of cytosine base pairing.
Collapse
Affiliation(s)
- Sung Eun Kim
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science, Seoul, Korea, Department of Physics, Korea University, Seoul 02841, Korea
| | - Seok-Cheol Hong
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science, Seoul, Korea, Department of Physics, Korea University, Seoul 02841, Korea
| |
Collapse
|
6
|
Nishio T, Masaoka T, Yoshikawa Y, Sadakane K, Kenmotsu T, Schiessel H, Yoshikawa K. Markedly Different Effects of Monovalent Cations on the Efficiency of Gene Expression. Adv Biol (Weinh) 2023; 7:e2200164. [PMID: 36328593 DOI: 10.1002/adbi.202200164] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2022] [Revised: 09/20/2022] [Indexed: 11/06/2022]
Abstract
The effect of monovalent cations on a cell-free transcription-translation (TX-TL) system is examined using a luciferase assay. It is found that the potency for all ions analyzed here is in the order Rb+ > K+ > Cs+ > Na+ ≈ Li+ > (CH3 )4 N+ , where Rb+ is most efficient at promoting TX-TL and the ions of Li+ , Na+ , and (CH3 )4 N+ exhibit an inhibitory effect. Similar promotion/inhibition effects are observed for cell-free TL alone with an mRNA template.
Collapse
Affiliation(s)
- Takashi Nishio
- Faculty of Life and Medical Sciences, Doshisha University, Kyoto, 610-0394, Japan
- Cluster of Excellence Physics of Life, Technical University of Dresden, 01307, Dresden, Germany
| | - Tomoya Masaoka
- Faculty of Life and Medical Sciences, Doshisha University, Kyoto, 610-0394, Japan
| | - Yuko Yoshikawa
- Faculty of Life and Medical Sciences, Doshisha University, Kyoto, 610-0394, Japan
| | - Koichiro Sadakane
- Faculty of Life and Medical Sciences, Doshisha University, Kyoto, 610-0394, Japan
| | - Takahiro Kenmotsu
- Faculty of Life and Medical Sciences, Doshisha University, Kyoto, 610-0394, Japan
| | - Helmut Schiessel
- Cluster of Excellence Physics of Life, Technical University of Dresden, 01307, Dresden, Germany
| | - Kenichi Yoshikawa
- Faculty of Life and Medical Sciences, Doshisha University, Kyoto, 610-0394, Japan
- Center for Integrative Medicine and Physics Institute for Advanced Study, Kyoto University, Kyoto, 606-8501, Japan
| |
Collapse
|
7
|
Chau C, Marcuccio F, Soulias D, Edwards MA, Tuplin A, Radford SE, Hewitt E, Actis P. Probing RNA Conformations Using a Polymer-Electrolyte Solid-State Nanopore. ACS NANO 2022; 16:20075-20085. [PMID: 36279181 PMCID: PMC9798860 DOI: 10.1021/acsnano.2c08312] [Citation(s) in RCA: 21] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/10/2023]
Abstract
Nanopore systems have emerged as a leading platform for the analysis of biomolecular complexes with single-molecule resolution. The conformation of biomolecules, such as RNA, is highly dependent on the electrolyte composition, but solid-state nanopore systems often require high salt concentration to operate, precluding analysis of macromolecular conformations under physiologically relevant conditions. Here, we report the implementation of a polymer-electrolyte solid-state nanopore system based on alkali metal halide salts dissolved in 50% w/v poly(ethylene) glycol (PEG) to augment the performance of our system. We show that polymer-electrolyte bath governs the translocation dynamics of the analyte which correlates with the physical properties of the salt used in the bath. This allowed us to identify CsBr as the optimal salt to complement PEG to generate the largest signal enhancement. Harnessing the effects of the polymer-electrolyte, we probed the conformations of the Chikungunya virus (CHIKV) RNA genome fragments under physiologically relevant conditions. Our system was able to fingerprint CHIKV RNA fragments ranging from ∼300 to ∼2000 nt length and subsequently distinguish conformations between the co-transcriptionally folded and the natively refolded ∼2000 nt CHIKV RNA. We envision that the polymer-electrolyte solid-state nanopore system will further enable structural and conformational analyses of individual biomolecules under physiologically relevant conditions.
Collapse
Affiliation(s)
- Chalmers Chau
- School
of Molecular and Cellular Biology and Astbury Centre for Structural
Molecular Biology, University of Leeds, Leeds LS2 9JT, U.K.
- School
of Electronic and Electrical Engineering and Pollard Institute, University of Leeds, Leeds LS2 9JT, U.K.
- Bragg
Centre for Materials Research, University
of Leeds, Leeds LS2 9JT, U.K.
| | - Fabio Marcuccio
- School
of Electronic and Electrical Engineering and Pollard Institute, University of Leeds, Leeds LS2 9JT, U.K.
- Bragg
Centre for Materials Research, University
of Leeds, Leeds LS2 9JT, U.K.
| | - Dimitrios Soulias
- School
of Electronic and Electrical Engineering and Pollard Institute, University of Leeds, Leeds LS2 9JT, U.K.
- Bragg
Centre for Materials Research, University
of Leeds, Leeds LS2 9JT, U.K.
| | - Martin Andrew Edwards
- Department
of Chemistry & Biochemistry, University
of Arkansas, Fayetteville, Arkansas 72701, United States
| | - Andrew Tuplin
- School
of Molecular and Cellular Biology and Astbury Centre for Structural
Molecular Biology, University of Leeds, Leeds LS2 9JT, U.K.
| | - Sheena E. Radford
- School
of Molecular and Cellular Biology and Astbury Centre for Structural
Molecular Biology, University of Leeds, Leeds LS2 9JT, U.K.
| | - Eric Hewitt
- School
of Molecular and Cellular Biology and Astbury Centre for Structural
Molecular Biology, University of Leeds, Leeds LS2 9JT, U.K.
| | - Paolo Actis
- School
of Electronic and Electrical Engineering and Pollard Institute, University of Leeds, Leeds LS2 9JT, U.K.
- Bragg
Centre for Materials Research, University
of Leeds, Leeds LS2 9JT, U.K.
| |
Collapse
|
8
|
Cruz-León S, Schwierz N. RNA Captures More Cations than DNA: Insights from Molecular Dynamics Simulations. J Phys Chem B 2022; 126:8646-8654. [PMID: 36260822 PMCID: PMC9639116 DOI: 10.1021/acs.jpcb.2c04488] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
The distribution of cations around nucleic acids is essential for a broad variety of processes ranging from DNA condensation and RNA folding to the detection of biomolecules in biosensors. Predicting the exact distribution of ions remains challenging since the distribution and, hence, a broad variety of nucleic acid properties depend on the salt concentration, the valency of the ions, and the ion type. Despite the importance, a general theory to quantify ion-specific effects for highly charged biomolecules is still lacking. Moreover, recent experiments reveal that despite their similar building blocks, DNA and RNA duplexes can react differently to the same ionic conditions. The aim of our current work is to provide a comprehensive set of molecular dynamics simulations using more than 180 μs of simulation time. For the mono- and divalent cations Li+, Na+, K+, Cs+, Ca2+, Sr2+, and Ba2+, the simulations allow us to reveal the ion-specific distributions and binding patterns for DNA and RNA duplexes. The microscopic insights from the simulations display the origin of ion-specificity and shed light on the question of why DNA and RNA show opposing behavior in the same ionic conditions. Finally, the detailed binding patterns from the simulations reveal why RNA can capture more cations than DNA.
Collapse
Affiliation(s)
- Sergio Cruz-León
- Department
of Theoretical Biophysics, Max Planck Institute
of Biophysics, Max-von-Laue-Str. 3, 60438Frankfurt am Main, Germany
| | - Nadine Schwierz
- Department
of Theoretical Biophysics, Max Planck Institute
of Biophysics, Max-von-Laue-Str. 3, 60438Frankfurt am Main, Germany,Institute
of Physics, University of Augsburg, Universitätsstraße 1, 86159Augsburg, Germany,E-mail:
| |
Collapse
|
9
|
Yu B, Bien KG, Wang T, Iwahara J. Diffusion NMR-based comparison of electrostatic influences of DNA on various monovalent cations. Biophys J 2022; 121:3562-3570. [PMID: 35754184 PMCID: PMC9515368 DOI: 10.1016/j.bpj.2022.06.022] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2022] [Revised: 06/03/2022] [Accepted: 06/22/2022] [Indexed: 11/18/2022] Open
Abstract
Counterions are important constituents for the structure and function of nucleic acids. Using 7Li and 133Cs nuclear magnetic resonance (NMR) spectroscopy, we investigated how ionic radii affect the behavior of counterions around DNA through diffusion measurements of Li+ and Cs+ ions around a 15-bp DNA duplex. Together with our previous data on 23Na+ and 15NH4+ ions around the same DNA under the same conditions, we were able to compare the dynamics of four different monovalent ions around DNA. From the apparent diffusion coefficients at varied concentrations of DNA, we determined the diffusion coefficients of these cations inside and outside the ion atmosphere around DNA (Db and Df, respectively). We also analyzed ionic competition with K+ ions for the ion atmosphere and assessed the relative affinities of these cations for DNA. Interestingly, all cations (i.e., Li+, Na+, NH4+, and Cs+) analyzed by diffusion NMR spectroscopy exhibited nearly identical Db/Df ratios despite the differences in their ionic radii, relative affinities, and diffusion coefficients. These results, along with the theoretical relationship between diffusion and entropy, suggest that the entropy change due to the release of counterions from the ion atmosphere around DNA is also similar regardless of the monovalent ion types. These findings and the experimental diffusion data on the monovalent ions are useful for examination of computational models for electrostatic interactions or ion solvation.
Collapse
Affiliation(s)
- Binhan Yu
- Department of Biochemistry & Molecular Biology, Sealy Center for Structural Biology & Molecular Biophysics, University of Texas Medical Branch, Galveston, Texas
| | - Karina G Bien
- Department of Biochemistry & Molecular Biology, Sealy Center for Structural Biology & Molecular Biophysics, University of Texas Medical Branch, Galveston, Texas
| | - Tianzhi Wang
- Department of Biochemistry & Molecular Biology, Sealy Center for Structural Biology & Molecular Biophysics, University of Texas Medical Branch, Galveston, Texas
| | - Junji Iwahara
- Department of Biochemistry & Molecular Biology, Sealy Center for Structural Biology & Molecular Biophysics, University of Texas Medical Branch, Galveston, Texas.
| |
Collapse
|
10
|
Ekesan Ş, McCarthy E, Case DA, York DM. RNA Electrostatics: How Ribozymes Engineer Active Sites to Enable Catalysis. J Phys Chem B 2022; 126:5982-5990. [PMID: 35862934 PMCID: PMC9496635 DOI: 10.1021/acs.jpcb.2c03727] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Electrostatic interactions are fundamental to RNA structure and function, and intimately influenced by solvation and the ion atmosphere. RNA enzymes, or ribozymes, are catalytic RNAs that are able to enhance reaction rates over a million-fold, despite having only a limited repertoire of building blocks and available set of chemical functional groups. Ribozyme active sites usually occur at junctions where negatively charged helices come together, and in many cases leverage this strained electrostatic environment to recruit metal ions in solution that can assist in catalysis. Similar strategies have been implicated in related artificially engineered DNA enzymes. Herein, we apply Poisson-Boltzmann, 3D-RISM, and molecular simulations to study a set of metal-dependent small self-cleaving ribozymes (hammerhead, pistol, and Varkud satellite) as well as an artificially engineered DNAzyme (8-17) to examine electrostatic features and their relation to the recruitment of monovalent and divalent metal ions important for activity. We examine several fundamental roles for these ions that include: (1) structural integrity of the catalytically active state, (2) pKa tuning of residues involved in acid-base catalysis, and (3) direct electrostatic stabilization of the transition state via Lewis acid catalysis. Taken together, these examples demonstrate how RNA electrostatics orchestrates the site-specific and territorial binding of metal ions to play important roles in catalysis.
Collapse
Affiliation(s)
- Şölen Ekesan
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | - Erika McCarthy
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | - David A. Case
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | - Darrin M. York
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| |
Collapse
|
11
|
Henning-Knechtel A, Thirumalai D, Kirmizialtin S. Differences in ion-RNA binding modes due to charge density variations explain the stability of RNA in monovalent salts. SCIENCE ADVANCES 2022; 8:eabo1190. [PMID: 35857829 PMCID: PMC9299541 DOI: 10.1126/sciadv.abo1190] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/08/2023]
Abstract
The stability of RNA increases as the charge density of the alkali metal cations increases. The molecular mechanism for this phenomenon remains elusive. To fill this gap, we performed all-atom molecular dynamics pulling simulations of HIV-1 trans-activation response RNA. We first established that the free energy landscape obtained in the simulations is in excellent agreement with the single-molecule optical tweezer experiments. The origin of the stronger stability in sodium compared to potassium is found to be due to the differences in the charge density-related binding modes. The smaller hydrated sodium ion preferentially binds to the highly charged phosphates that have high surface area. In contrast, the larger potassium ions interact with the major grooves. As a result, more cations condense around phosphate groups in the case of sodium ions, leading to the reduction of electrostatic repulsion. Because the proposed mechanism is generic, we predict that the same conclusions are valid for divalent alkaline earth metal cations.
Collapse
Affiliation(s)
| | - D. Thirumalai
- Department of Chemistry, University of Texas, Austin, TX 78712, USA
- Corresponding author. (D.T.); (S.K.)
| | - Serdal Kirmizialtin
- Chemistry Program, Math and Sciences, New York University Abu Dhabi, Abu Dhabi, UAE
- Corresponding author. (D.T.); (S.K.)
| |
Collapse
|
12
|
Cruz-León S, Vanderlinden W, Müller P, Forster T, Staudt G, Lin YY, Lipfert J, Schwierz N. Twisting DNA by salt. Nucleic Acids Res 2022; 50:5726-5738. [PMID: 35640616 PMCID: PMC9177979 DOI: 10.1093/nar/gkac445] [Citation(s) in RCA: 24] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2021] [Revised: 05/06/2022] [Accepted: 05/10/2022] [Indexed: 12/11/2022] Open
Abstract
The structure and properties of DNA depend on the environment, in particular the ion atmosphere. Here, we investigate how DNA twist -one of the central properties of DNA- changes with concentration and identity of the surrounding ions. To resolve how cations influence the twist, we combine single-molecule magnetic tweezer experiments and extensive all-atom molecular dynamics simulations. Two interconnected trends are observed for monovalent alkali and divalent alkaline earth cations. First, DNA twist increases monotonously with increasing concentration for all ions investigated. Second, for a given salt concentration, DNA twist strongly depends on cation identity. At 100 mM concentration, DNA twist increases as Na+ < K+ < Rb+ < Ba2+ < Li+ ≈ Cs+ < Sr2+ < Mg2+ < Ca2+. Our molecular dynamics simulations reveal that preferential binding of the cations to the DNA backbone or the nucleobases has opposing effects on DNA twist and provides the microscopic explanation of the observed ion specificity. However, the simulations also reveal shortcomings of existing force field parameters for Cs+ and Sr2+. The comprehensive view gained from our combined approach provides a foundation for understanding and predicting cation-induced structural changes both in nature and in DNA nanotechnology.
Collapse
Affiliation(s)
- Sergio Cruz-León
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics, Max-von-Laue-Str. 3, 60438 Frankfurt am Main, Germany
| | - Willem Vanderlinden
- Department of Physics and Center for Nanoscience (CeNS), LMU Munich, Amalienstr. 54, 80799 Munich, Germany
| | - Peter Müller
- Department of Physics and Center for Nanoscience (CeNS), LMU Munich, Amalienstr. 54, 80799 Munich, Germany
| | - Tobias Forster
- Department of Physics and Center for Nanoscience (CeNS), LMU Munich, Amalienstr. 54, 80799 Munich, Germany
| | - Georgina Staudt
- Department of Physics and Center for Nanoscience (CeNS), LMU Munich, Amalienstr. 54, 80799 Munich, Germany
| | - Yi-Yun Lin
- Department of Physics and Center for Nanoscience (CeNS), LMU Munich, Amalienstr. 54, 80799 Munich, Germany
| | - Jan Lipfert
- Department of Physics and Center for Nanoscience (CeNS), LMU Munich, Amalienstr. 54, 80799 Munich, Germany
| | - Nadine Schwierz
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics, Max-von-Laue-Str. 3, 60438 Frankfurt am Main, Germany
| |
Collapse
|
13
|
Zhou Y, Jiang Y, Chen SJ. RNA-ligand molecular docking: advances and challenges. WILEY INTERDISCIPLINARY REVIEWS. COMPUTATIONAL MOLECULAR SCIENCE 2022; 12:e1571. [PMID: 37293430 PMCID: PMC10250017 DOI: 10.1002/wcms.1571] [Citation(s) in RCA: 23] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/26/2021] [Accepted: 07/20/2021] [Indexed: 12/16/2022]
Abstract
With rapid advances in computer algorithms and hardware, fast and accurate virtual screening has led to a drastic acceleration in selecting potent small molecules as drug candidates. Computational modeling of RNA-small molecule interactions has become an indispensable tool for RNA-targeted drug discovery. The current models for RNA-ligand binding have mainly focused on the docking-and-scoring method. Accurate docking and scoring should tackle four crucial problems: (1) conformational flexibility of ligand, (2) conformational flexibility of RNA, (3) efficient sampling of binding sites and binding poses, and (4) accurate scoring of different binding modes. Moreover, compared with the problem of protein-ligand docking, predicting ligand binding to RNA, a negatively charged polymer, is further complicated by additional effects such as metal ion effects. Thermodynamic models based on physics-based and knowledge-based scoring functions have shown highly encouraging success in predicting ligand binding poses and binding affinities. Recently, kinetic models for ligand binding have further suggested that including dissociation kinetics (residence time) in ligand docking would result in improved performance in estimating in vivo drug efficacy. More recently, the rise of deep-learning approaches has led to new tools for predicting RNA-small molecule binding. In this review, we present an overview of the recently developed computational methods for RNA-ligand docking and their advantages and disadvantages.
Collapse
Affiliation(s)
- Yuanzhe Zhou
- Department of Physics and Astronomy, Department of Biochemistry, Institute of Data Sciences and Informatics, University of Missouri, Columbia, MO 65211-7010, USA
| | - Yangwei Jiang
- Department of Physics and Astronomy, Department of Biochemistry, Institute of Data Sciences and Informatics, University of Missouri, Columbia, MO 65211-7010, USA
| | - Shi-Jie Chen
- Department of Physics and Astronomy, Department of Biochemistry, Institute of Data Sciences and Informatics, University of Missouri, Columbia, MO 65211-7010, USA
| |
Collapse
|
14
|
Hamilton I, Gebala M, Herschlag D, Russell R. Direct Measurement of Interhelical DNA Repulsion and Attraction by Quantitative Cross-Linking. J Am Chem Soc 2022; 144:1718-1728. [PMID: 35073489 PMCID: PMC8815069 DOI: 10.1021/jacs.1c11122] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2021] [Indexed: 12/30/2022]
Abstract
To better understand the forces that mediate nucleic acid compaction in biology, we developed the disulfide cross-linking approach xHEED (X-linking of Helices to measure Electrostatic Effects at Distance) to measure the distance-dependent encounter frequency of two DNA helices in solution. Using xHEED, we determined the distance that the electrostatic potential extends from DNA helices, the dependence of this distance on ionic conditions, and the magnitude of repulsion when two helices approach one another. Across all conditions tested, the potential falls to that of the bulk solution within 15 Å of the major groove surface. For separations of ∼30 Å, we measured a repulsion of 1.8 kcal/mol in low monovalent ion concentration (30 mM Na+), with higher Na+ concentrations ameliorating this repulsion, and 2 M Na+ or 100 mM Mg2+ eliminating it. Strikingly, we found full screening at very low Co3+ concentrations and net attraction at higher concentrations, without the higher-order DNA condensation that typically complicates studies of helical attraction. Our measurements define the relevant distances for electrostatic interactions of nucleic-acid helices in biology and introduce a new method to propel further understanding of how these forces impact biological processes.
Collapse
Affiliation(s)
- Ian Hamilton
- Department
of Molecular Biosciences, University of
Texas at Austin, Austin, Texas 78712, United States
| | - Magdalena Gebala
- Department
of Biochemistry, Stanford University, Stanford California 94305, United States
| | - Daniel Herschlag
- Department
of Biochemistry, Stanford University, Stanford California 94305, United States
| | - Rick Russell
- Department
of Molecular Biosciences, University of
Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
15
|
Jambrec D, Gebala M. DNA Electrostatics: From Theory to Application. ChemElectroChem 2022. [DOI: 10.1002/celc.202101415] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Daliborka Jambrec
- Analytische Chemie – Elektroanalytik & Sensorik Ruhr-Universität Bochum Universitätsstr. 150 D-44780 Bochum Germany
| | - Magdalena Gebala
- Department of Biochemistry Stanford University Stanford 94305, CA USA
| |
Collapse
|
16
|
Yu B, Bien KG, Pletka CC, Iwahara J. Dynamics of Cations around DNA and Protein as Revealed by 23Na Diffusion NMR Spectroscopy. Anal Chem 2022; 94:2444-2452. [PMID: 35080384 PMCID: PMC8829827 DOI: 10.1021/acs.analchem.1c04197] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Counterions are vital for the structure and function of biomolecules. However, the behavior of counterions remains elusive due to the difficulty in characterizing mobile ions. Here, we demonstrate that the dynamics of cations around biological macromolecules can be revealed by 23Na diffusion nuclear magnetic resonance (NMR) spectroscopy. NMR probe hardware capable of generating strong magnetic field gradients enables 23Na NMR-based diffusion measurements for Na+ ions in solutions of biological macromolecules and their complexes. The dynamic properties of Na+ ions interacting with the macromolecules can be investigated using apparent 23Na diffusion coefficients measured under various conditions. Our diffusion data clearly show that Na+ ions retain high mobility within the ion atmosphere around DNA. The 23Na diffusion NMR method also permits direct observation of the release of Na+ ions from nucleic acids upon protein-nucleic acid association. The entropy change due to the ion release can be estimated from the diffusion data.
Collapse
Affiliation(s)
- Binhan Yu
- Department of Biochemistry and Molecular Biology, Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, Galveston, Texas 77555-1068 United States
| | - Karina G Bien
- Department of Biochemistry and Molecular Biology, Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, Galveston, Texas 77555-1068 United States
| | - Channing C Pletka
- Department of Biochemistry and Molecular Biology, Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, Galveston, Texas 77555-1068 United States
| | - Junji Iwahara
- Department of Biochemistry and Molecular Biology, Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, Galveston, Texas 77555-1068 United States
| |
Collapse
|
17
|
Kanapeckaitė A, Burokienė N, Mažeikienė A, Cottrell GS, Widera D. Biophysics is reshaping our perception of the epigenome: from DNA-level to high-throughput studies. BIOPHYSICAL REPORTS 2021; 1:100028. [PMID: 36425454 PMCID: PMC9680810 DOI: 10.1016/j.bpr.2021.100028] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/26/2021] [Accepted: 09/24/2021] [Indexed: 06/16/2023]
Abstract
Epigenetic research holds great promise to advance our understanding of biomarkers and regulatory processes in health and disease. An increasing number of new approaches, ranging from molecular to biophysical analyses, enable identifying epigenetic changes on the level of a single gene or the whole epigenome. The aim of this review is to highlight how the field is shifting from completely molecular-biology-driven solutions to multidisciplinary strategies including more reliance on biophysical analysis tools. Biophysics not only offers technical advancements in imaging or structure analysis but also helps to explore regulatory interactions. New computational methods are also being developed to meet the demand of growing data volumes and their processing. Therefore, it is important to capture these new directions in epigenetics from a biophysical perspective and discuss current challenges as well as multiple applications of biophysical methods and tools. Specifically, we gradually introduce different biophysical research methods by first considering the DNA-level information and eventually higher-order chromatin structures. Moreover, we aim to highlight that the incorporation of bioinformatics, machine learning, and artificial intelligence into biophysical analysis allows gaining new insights into complex epigenetic processes. The gained understanding has already proven useful in translational and clinical research providing better patient stratification options or new therapeutic insights. Together, this offers a better readiness to transform bench-top experiments into industrial high-throughput applications with a possibility to employ developed methods in clinical practice and diagnostics.
Collapse
Affiliation(s)
- Austė Kanapeckaitė
- Algorithm379, Laisvės g. 7, LT 12007, Vilnius, Lithuania
- Reading School of Pharmacy, Whiteknights, Reading, UK, RG6 6UB
| | - Neringa Burokienė
- Clinics of Internal Diseases, Family Medicine and Oncology, Institute of Clinical Medicine, Faculty of Medicine, Vilnius University, M. K. Čiurlionio str. 21/27, LT-03101 Vilnius, Lithuania
| | - Asta Mažeikienė
- Department of Physiology, Biochemistry, Microbiology and Laboratory Medicine, Institute of Biomedical Sciences, Faculty of Medicine, M. K. Čiurlionio str. 21/27, LT-03101 Vilnius, Lithuania
| | | | - Darius Widera
- Reading School of Pharmacy, Whiteknights, Reading, UK, RG6 6UB
| |
Collapse
|
18
|
Comparison of the ionic effects of Ca2+ and Mg2+ on nucleic acids in liquids. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.117781] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
19
|
Ma CY, Pezzotti S, Schwaab G, Gebala M, Herschlag D, Havenith M. Cation enrichment in the ion atmosphere is promoted by local hydration of DNA. Phys Chem Chem Phys 2021; 23:23203-23213. [PMID: 34622888 PMCID: PMC8797164 DOI: 10.1039/d1cp01963e] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Electrostatic interactions are central to the structure and function of nucleic acids, including their folding, condensation, and interaction with proteins and other charged molecules. These interactions are profoundly affected by ions surrounding nucleic acids, the constituents of the so-called ion atmosphere. Here, we report precise Fourier Transform-Terahertz/Far-Infrared (FT-THz/FIR) measurements in the frequency range 30-500 cm-1 for a 24-bp DNA solvated in a series of alkali halide (NaCl, NaF, KCl, CsCl, and CsF) electrolyte solutions which are sensitive to changes in the ion atmosphere. Cation excess in the ion atmosphere is detected experimentally by observation of cation modes of Na+, K+, and Cs+ in the frequency range between 70-90 cm-1. Based on MD simulations, we propose that the magnitude of cation excess (which is salt specific) depends on the ability of the electrolyte to perturb the water network at the DNA interface: In the NaF atmosphere, the ions reduce the strength of interactions between water and the DNA more than in case of a NaCl electrolyte. Here, we explicitly take into account the solvent contribution to the chemical potential in the ion atmosphere: A decrease in the number of bound water molecules in the hydration layer of DNA is correlated with enhanced density fluctuations, which decrease the free energy cost of ion-hydration, thus promoting further ion accumulation within the DNA atmosphere. We propose that taking into account the local solvation is crucial for understanding the ion atmosphere.
Collapse
Affiliation(s)
- Chun Yu Ma
- Department of Physical Chemistry II, Ruhr-University Bochum, 44780 Bochum, Germany.
| | - Simone Pezzotti
- Department of Physical Chemistry II, Ruhr-University Bochum, 44780 Bochum, Germany.
| | - Gerhard Schwaab
- Department of Physical Chemistry II, Ruhr-University Bochum, 44780 Bochum, Germany.
| | - Magdalena Gebala
- Department of Biochemistry, Stanford University, Stanford, California 94305, USA
| | - Daniel Herschlag
- Department of Biochemistry, Stanford University, Stanford, California 94305, USA
| | - Martina Havenith
- Department of Physical Chemistry II, Ruhr-University Bochum, 44780 Bochum, Germany.
| |
Collapse
|
20
|
Boychuk BTA, Rebecca Jeong YE, Wetmore SD. Assessment of the Accuracy of DFT-Predicted Li +-Nucleic Acid Binding Energies. J Chem Theory Comput 2021; 17:5392-5408. [PMID: 34339194 DOI: 10.1021/acs.jctc.1c00401] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Understanding how lithium interacts with complex biosystems is crucial for uncovering the roles of this alkali metal in biology and designing extraction techniques for battery production and environmental remediation. In this light, fundamental information about Li+ binding to nucleic acids is required. Herein, a new database of Li+-nucleic acid interactions is presented that contains CCSD(T)/CBS benchmark energies for all nucleobase and phosphate binding locations. Furthermore, the performance of 54 DFT functionals in combination with three triple-zeta (TZ) basis sets (6-311+G(3df,2p), aug-cc-pVTZ, and def2-TZVPP) is tested. The results identify a range of functionals across different families (B2-PLYP, PBE-QIDH, ωB97, ωB97X-D, MN15, B3PW91, B97-2, TPSS, BP86-D3(BJ), and PBE) that can accurately describe coordinated Li+-nucleic acid interactions, with the average mean percent error (AMPE) across binding positions and basis sets being below 2%. Nevertheless, only three functionals tested (B2-PLYP, PBE-QIDH, and ωB97X-D) preserve this accuracy for metal cation-π interactions, suggesting that caution is warranted when choosing a functional to describe a diverse range of Li+-nucleic acid complexes. Removal of counterpoise corrections has very little impact on the reliability of most functionals, while the effect of empirical dispersion corrections varies depending on the functional choice and interaction type. While increasing the basis set to quadruple-zeta quality had little impact on the AMPE, the accuracy of double-zeta basis sets varies with family. Importantly, DFT methods reproduce the CCSD(T)/CBS trend in the preferred binding position for a given nucleic acid component and the global trend across components (phosphate ≫ G > C ≫ A ∼ T = U), as well as the geometries of the metal-nucleic acid complexes. The overall top performing functional is PBE-QIDH, which results in deviations from CCSD(T)/CBS values as small as ∼0.1 kcal/mol for nucleobase contacts and ∼1 kcal/mol for phosphate interactions. The most accurate DFT methods identified in the present work are recommended for future investigations of lithium interactions in larger nucleic acid systems to provide insights into the biological roles of this metal and the design of novel biosensing strategies.
Collapse
Affiliation(s)
- Briana T A Boychuk
- Department of Chemistry and Biochemistry, University of Lethbridge, 4401 University Drive West, Lethbridge, Alberta T1K 3M4, Canada
| | - Ye Eun Rebecca Jeong
- Department of Chemistry and Biochemistry, University of Lethbridge, 4401 University Drive West, Lethbridge, Alberta T1K 3M4, Canada
| | - Stacey D Wetmore
- Department of Chemistry and Biochemistry, University of Lethbridge, 4401 University Drive West, Lethbridge, Alberta T1K 3M4, Canada
| |
Collapse
|
21
|
Abstract
Electrostatic potentials computed from three-dimensional structures of biomolecules by solving the Poisson-Boltzmann equation are widely used in molecular biophysics, structural biology, and medicinal chemistry. Despite the approximate nature of the Poisson-Boltzmann theory, validation of the computed electrostatic potentials around biological macromolecules is rare and methodologically limited. Here, we present a unique and powerful NMR method that allows for straightforward and extensive comparison with electrostatic models for biomolecules and their complexes. This method utilizes paramagnetic relaxation enhancement arising from analogous cationic and anionic cosolutes whose spatial distributions around biological macromolecules reflect electrostatic potentials. We demonstrate that this NMR method enables de novo determination of near-surface electrostatic potentials for individual protein residues without using any structural information. We applied the method to ubiquitin and the Antp homeodomain-DNA complex. The experimental data agreed well with predictions from the Poisson-Boltzmann theory. Thus, our experimental results clearly support the validity of the theory for these systems. However, our experimental study also illuminates certain weaknesses of the Poisson-Boltzmann theory. For example, we found that the theory predicts stronger dependence of near-surface electrostatic potentials on ionic strength than observed in the experiments. Our data also suggest that conformational flexibility or structural uncertainties may cause large errors in theoretical predictions of electrostatic potentials, particularly for highly charged systems. This NMR-based method permits extensive assessment of near-surface electrostatic potentials for various regions around biological macromolecules and thereby may facilitate improvement of the computational approaches for electrostatic potentials.
Collapse
|
22
|
Yu B, Iwahara J. Experimental approaches for investigating ion atmospheres around nucleic acids and proteins. Comput Struct Biotechnol J 2021; 19:2279-2285. [PMID: 33995919 PMCID: PMC8102144 DOI: 10.1016/j.csbj.2021.04.033] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2021] [Revised: 04/14/2021] [Accepted: 04/14/2021] [Indexed: 01/26/2023] Open
Abstract
Ionic interactions are crucial to biological functions of DNA, RNA, and proteins. Experimental research on how ions behave around biological macromolecules has lagged behind corresponding theoretical and computational research. In the 21st century, quantitative experimental approaches for investigating ionic interactions of biomolecules have become available and greatly facilitated examinations of theoretical electrostatic models. These approaches utilize anomalous small-angle X-ray scattering, atomic emission spectroscopy, mass spectrometry, or nuclear magnetic resonance (NMR) spectroscopy. We provide an overview on the experimental methodologies that can quantify and characterize ions within the ion atmospheres around nucleic acids, proteins, and their complexes.
Collapse
Affiliation(s)
- Binhan Yu
- Department of Biochemistry and Molecular Biology, Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, Galveston, TX 77555-1068, USA
| | - Junji Iwahara
- Department of Biochemistry and Molecular Biology, Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, Galveston, TX 77555-1068, USA
| |
Collapse
|
23
|
Abstract
![]()
Molecular association of proteins with nucleic
acids is required
for many biological processes essential to life. Electrostatic interactions
via ion pairs (salt bridges) of nucleic acid phosphates and protein
side chains are crucial for proteins to bind to DNA or RNA. Counterions
around the macromolecules are also key constituents for the thermodynamics
of protein–nucleic acid association. Until recently, there
had been only a limited amount of experiment-based information about
how ions and ionic moieties behave in biological macromolecular processes.
In the past decade, there has been significant progress in quantitative
experimental research on ionic interactions with nucleic acids and
their complexes with proteins. The highly negatively charged surfaces
of DNA and RNA electrostatically attract and condense cations, creating
a zone called the ion atmosphere. Recent experimental studies were
able to examine and validate theoretical models on ions and their
mobility and interactions with macromolecules. The ionic interactions
are highly dynamic. The counterions rapidly diffuse within the ion
atmosphere. Some of the ions are released from the ion atmosphere
when proteins bind to nucleic acids, balancing the charge via intermolecular
ion pairs of positively charged side chains and negatively charged
backbone phosphates. Previously, the release of counterions had been
implicated indirectly by the salt-concentration dependence of the
association constant. Recently, direct detection of counterion
release by NMR spectroscopy
has become possible and enabled more accurate and quantitative analysis
of the counterion release and its entropic impact on the thermodynamics
of protein–nucleic acid association. Recent studies also revealed
the dynamic nature of ion pairs of protein side chains and nucleic
acid phosphates. These ion pairs undergo transitions between two major
states. In one of the major states, the cation and the anion are in
direct contact and form hydrogen bonds. In the other major state,
the cation and the anion are separated by water. Transitions between
these states rapidly occur on a picosecond to nanosecond time scale.
When proteins interact with nucleic acids, interfacial arginine (Arg)
and lysine (Lys) side chains exhibit considerably different behaviors.
Arg side chains show a higher propensity to form rigid contacts with
nucleotide bases, whereas Lys side chains tend to be more mobile at
the molecular interfaces. The dynamic ionic interactions may facilitate
adaptive molecular recognition and play both thermodynamic and kinetic
roles in protein–nucleic acid interactions.
Collapse
Affiliation(s)
- Binhan Yu
- Department of Biochemistry & Molecular Biology, Sealy Center for Structural Biology & Molecular Biophysics, University of Texas Medical Branch, Galveston, Texas 77555-1068, United States
| | - B. Montgomery Pettitt
- Department of Biochemistry & Molecular Biology, Sealy Center for Structural Biology & Molecular Biophysics, University of Texas Medical Branch, Galveston, Texas 77555-1068, United States
| | - Junji Iwahara
- Department of Biochemistry & Molecular Biology, Sealy Center for Structural Biology & Molecular Biophysics, University of Texas Medical Branch, Galveston, Texas 77555-1068, United States
| |
Collapse
|
24
|
Cruz-León S, Schwierz N. Hofmeister Series for Metal-Cation-RNA Interactions: The Interplay of Binding Affinity and Exchange Kinetics. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2020; 36:5979-5989. [PMID: 32366101 PMCID: PMC7304902 DOI: 10.1021/acs.langmuir.0c00851] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
A large variety of physicochemical properties involving RNA depends on the type of metal cation present in solution. In order to gain microscopic insight into the origin of these ion specific effects, we apply molecular dynamics simulations to describe the interactions of metal cations and RNA. For the three most common ion binding sites on RNA, we calculate the binding affinities and exchange rates of eight different mono- and divalent metal cations. Our results reveal that binding sites involving phosphate groups preferentially bind metal cations with high charge density (such as Mg2+) in inner-sphere conformations while binding sites involving N7 or O6 atoms preferentially bind cations with low charge density (such as K+). The binding affinity therefore follows a direct Hofmeister series at the backbone but is reversed at the nucleobases leading to a high selectivity of ion binding sites on RNA. In addition, the exchange rates for cation binding cover almost 5 orders of magnitude, leading to a vastly different time scale for the lifetimes of contact pairs. Taken together, the site-specific binding affinities and the specific lifetime of contact pairs provide the microscopic explanation of ion specific effects observed in a wide variety of macroscopic RNA properties. Finally, combining the results from atomistic simulations with extended Poisson-Boltzmann theory allows us to predict the distribution of metal cations around double-stranded RNA at finite concentrations and to reproduce the results of ion counting experiments with good accuracy.
Collapse
|
25
|
Sung HL, Nesbitt DJ. Single-molecule kinetic studies of DNA hybridization under extreme pressures. Phys Chem Chem Phys 2020; 22:23491-23501. [DOI: 10.1039/d0cp04035e] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Pressure-responsive dynamics of DNA hairpin hybridization/dehybridization is directly visualized at the single molecule level.
Collapse
Affiliation(s)
- Hsuan-Lei Sung
- JILA
- National Institute of Standards and Technology and University of Colorado
- Boulder
- USA
- Department of Chemistry
| | - David J. Nesbitt
- JILA
- National Institute of Standards and Technology and University of Colorado
- Boulder
- USA
- Department of Chemistry
| |
Collapse
|
26
|
Plumridge A, Andresen K, Pollack L. Visualizing Disordered Single-Stranded RNA: Connecting Sequence, Structure, and Electrostatics. J Am Chem Soc 2019; 142:109-119. [PMID: 31804813 DOI: 10.1021/jacs.9b04461] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Disordered homopolymeric regions of single-stranded RNA, such as U or A tracts, are found within functional RNAs where they play distinct roles in defining molecular structure and facilitating recognition by partners. Despite this prominence, details of conformational and biophysical properties of these regions have not yet been resolved. We apply a number of experimental techniques to investigate the conformations of these biologically important motifs and provide quantitative measurements of their ion atmospheres. Single strands of RNA display pronounced sequence-dependent conformations that relate to the unique ion atmospheres each attracts. Chains of rU bases are relatively unstructured under all conditions, while chains of rA bases display distinct ordering through stacking or clustering motifs, depending on the composition of the surrounding solution. These dramatic structural differences are consistent with the measured disparity in ion composition and atmospheres around each homopolymer, revealing a complex interplay of base, ion, and single-strand ordering. The unique structural and ionic signatures of homopolymer ssRNAs explains their role(s) in folding structured RNAs and may explain their distinct recognition by protein partners.
Collapse
Affiliation(s)
- Alex Plumridge
- School of Applied and Engineering Physics , Cornell University , Ithaca , New York 14853 , United States
| | - Kurt Andresen
- Department of Physics , Gettysburg College , Gettysburg , Pennsylvania 17325 , United States
| | - Lois Pollack
- School of Applied and Engineering Physics , Cornell University , Ithaca , New York 14853 , United States
| |
Collapse
|
27
|
Sun LZ, Zhou Y, Chen SJ. Predicting Monovalent Ion Correlation Effects in Nucleic Acids. ACS OMEGA 2019; 4:13435-13446. [PMID: 31460472 PMCID: PMC6705202 DOI: 10.1021/acsomega.9b01689] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/09/2019] [Accepted: 07/18/2019] [Indexed: 05/14/2023]
Abstract
Ion correlation and fluctuation can play a potentially significant role in metal ion-nucleic acid interactions. Previous studies have focused on the effects for multivalent cations. However, the correlation and fluctuation effects can be important also for monovalent cations around the nucleic acid surface. Here, we report a model, gMCTBI, that can explicitly treat discrete distributions of both monovalent and multivalent cations and can account for the correlation and fluctuation effects for the cations in the solution. The gMCTBI model enables investigation of the global ion binding properties as well as the detailed discrete distributions of the bound ions. Accounting for the ion correlation effect for monovalent ions can lead to more accurate predictions, especially in a mixed monovalent and multivalent salt solution, for the number and location of the bound ions. Furthermore, although the monovalent ion-mediated correlation does not show a significant effect on the number of bound ions, the correlation may enhance the accumulation of monovalent ions near the nucleic acid surface and hence affect the ion distribution. The study further reveals novel ion correlation-induced effects in the competition between the different cations around nucleic acids.
Collapse
Affiliation(s)
- Li-Zhen Sun
- Department
of Applied Physics, Zhejiang University
of Technology, Hangzhou 310023, China
- Department
of Physics, Department of Biochemistry, and Informatics Institute, University of Missouri, Columbia, Missouri 65211, United States
| | - Yuanzhe Zhou
- Department
of Physics, Department of Biochemistry, and Informatics Institute, University of Missouri, Columbia, Missouri 65211, United States
| | - Shi-Jie Chen
- Department
of Physics, Department of Biochemistry, and Informatics Institute, University of Missouri, Columbia, Missouri 65211, United States
- E-mail:
| |
Collapse
|
28
|
Quantitative Studies of an RNA Duplex Electrostatics by Ion Counting. Biophys J 2019; 117:1116-1124. [PMID: 31466697 DOI: 10.1016/j.bpj.2019.08.007] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2019] [Revised: 08/01/2019] [Accepted: 08/05/2019] [Indexed: 01/22/2023] Open
Abstract
RNAs are one of the most charged polyelectrolytes in nature, and understanding their electrostatics is fundamental to their structure and biological functions. An effective way to characterize the electrostatic field generated by nucleic acids is to quantify interactions between nucleic acids and ions that surround the molecules. These ions form a loosely associated cloud referred to as an ion atmosphere. Although theoretical and computational studies can describe the ion atmosphere around RNAs, benchmarks are needed to guide the development of these approaches, and experiments to date that read out RNA-ion interactions are limited. Here, we present ion counting studies to quantify the number of ions surrounding well-defined model systems of RNA and DNA duplexes. We observe that the RNA duplex attracts more cations and expels fewer anions compared to the DNA duplex, and the RNA duplex interacts significantly stronger with the divalent cation Mg2+, despite their identical total charge. These experimental results suggest that the RNA duplex generates a stronger electrostatic field than DNA, as is predicted based on the structural differences between their helices. Theoretical calculations using a nonlinear Poisson-Boltzmann equation give excellent agreement with experiments for monovalent ions but underestimate Mg2+-DNA and Mg2+-RNA interactions by 20%. These studies provide needed stringent benchmarks to use against other all-atom theoretical models of RNA-ion interactions, interactions that likely must be accurately accounted for in structural, dynamic, and energetic terms to confidently model RNA structure, interactions, and function.
Collapse
|
29
|
Gebala M, Johnson SL, Narlikar GJ, Herschlag D. Ion counting demonstrates a high electrostatic field generated by the nucleosome. eLife 2019; 8:e44993. [PMID: 31184587 PMCID: PMC6584128 DOI: 10.7554/elife.44993] [Citation(s) in RCA: 36] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2019] [Accepted: 06/08/2019] [Indexed: 01/04/2023] Open
Abstract
In eukaryotes, a first step towards the nuclear DNA compaction process is the formation of a nucleosome, which is comprised of negatively charged DNA wrapped around a positively charged histone protein octamer. Often, it is assumed that the complexation of the DNA into the nucleosome completely attenuates the DNA charge and hence the electrostatic field generated by the molecule. In contrast, theoretical and computational studies suggest that the nucleosome retains a strong, negative electrostatic field. Despite their fundamental implications for chromatin organization and function, these opposing views of nucleosome electrostatics have not been experimentally tested. Herein, we directly measure nucleosome electrostatics and find that while nucleosome formation reduces the complex charge by half, the nucleosome nevertheless maintains a strong negative electrostatic field. Our studies highlight the importance of considering the polyelectrolyte nature of the nucleosome and its impact on processes ranging from factor binding to DNA compaction.
Collapse
Affiliation(s)
- Magdalena Gebala
- Department of BiochemistryStanford UniversityStanfordUnited States
| | - Stephanie L Johnson
- Department of Biochemistry and BiophysicsUniversity of California, San FranciscoSan FranciscoUnited States
| | - Geeta J Narlikar
- Department of Biochemistry and BiophysicsUniversity of California, San FranciscoSan FranciscoUnited States
| | - Dan Herschlag
- Department of BiochemistryStanford UniversityStanfordUnited States
- Department of ChemistryStanford UniversityStanfordUnited States
- ChEM-H InstituteStanford UniversityStanfordUnited States
| |
Collapse
|
30
|
Giambasu GM, Case DA, York DM. Predicting Site-Binding Modes of Ions and Water to Nucleic Acids Using Molecular Solvation Theory. J Am Chem Soc 2019; 141:2435-2445. [PMID: 30632365 PMCID: PMC6574206 DOI: 10.1021/jacs.8b11474] [Citation(s) in RCA: 41] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Site binding of ions and water shapes nucleic acids folding, dynamics, and biological function, complementing the more diffuse, nonspecific "territorial" ion binding. Unlike territorial binding, prediction of site-specific binding to nucleic acids remains an unsolved challenge in computational biophysics. This work presents a new toolset based on the 3D-RISM molecular solvation theory and topological analysis that predicts cation and water site binding to nucleic acids. 3D-RISM is shown to accurately capture alkali cations and water binding to the central channel, transversal loops, and grooves of the Oxytricha nova's telomeres' G-quadruplex ( Oxy-GQ), in agreement with high-resolution crystallographic data. To improve the computed cation occupancy along the Oxy-GQ central channel, it was necessary to refine and validate new cation-oxygen parameters using structural and thermodynamic data available for crown ethers and ion channels. This single set of parameters that describes both localized and delocalized binding to various biological systems is used to gain insight into cation occupancy along the Oxy-GQ channel under various salt conditions. The paper concludes with prospects for extending the method to predict divalent cation binding to nucleic acids. This work advances the forefront of theoretical methods able to provide predictive insight into ion atmosphere effects on nucleic acids function.
Collapse
Affiliation(s)
- George M. Giambasu
- Department of Chemistry and Chemical Biology, Rutgers, The State University of New Jersey, New Brunswick, New Jersey 08901, United States
| | - David A. Case
- Department of Chemistry and Chemical Biology, Rutgers, The State University of New Jersey, New Brunswick, New Jersey 08901, United States
| | - Darrin M. York
- Department of Chemistry and Chemical Biology, Rutgers, The State University of New Jersey, New Brunswick, New Jersey 08901, United States
- Laboratory for Biomolecular Simulation Research, Rutgers, The State University of New Jersey, New Brunswick, New Jersey 08901, United States
- Center for Integrative Proteomics Research, Rutgers, The State University of New Jersey, New Brunswick, New Jersey 08901, United States
| |
Collapse
|
31
|
Guilbaud S, Salomé L, Destainville N, Manghi M, Tardin C. Dependence of DNA Persistence Length on Ionic Strength and Ion Type. PHYSICAL REVIEW LETTERS 2019; 122:028102. [PMID: 30720315 DOI: 10.1103/physrevlett.122.028102] [Citation(s) in RCA: 53] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/28/2018] [Indexed: 06/09/2023]
Abstract
Even though the persistence length L_{P} of double-stranded DNA plays a pivotal role in cell biology and nanotechnologies, its dependence on ionic strength I lacks a consensual description. Using a high-throughput single-molecule technique and statistical physics modeling, we measure L_{P} in the presence of monovalent (Li^{+}, Na^{+}, K^{+}) and divalent (Mg^{2+}, Ca^{2+}) metallic and alkyl ammonium ions, over a large range 0.5 mM≤I≤5 M. We show that linear Debye-Hückel-type theories do not describe even part of these data. By contrast, the Netz-Orland and Trizac-Shen formulas, two approximate theories including nonlinear electrostatic effects and the finite DNA radius, fit our data with divalent and monovalent ions, respectively, over the whole I range. Furthermore, the metallic ion type does not influence L_{P}(I), in contrast to alkyl ammonium monovalent ions at high I.
Collapse
Affiliation(s)
- Sébastien Guilbaud
- Institut de Pharmacologie et de Biologie Structurale, Université de Toulouse, CNRS, UPS, 31 077 Toulouse, France
| | - Laurence Salomé
- Institut de Pharmacologie et de Biologie Structurale, Université de Toulouse, CNRS, UPS, 31 077 Toulouse, France
| | - Nicolas Destainville
- Laboratoire de Physique Théorique (IRSAMC), Université de Toulouse, CNRS, UPS, 31 062 Toulouse, France
| | - Manoel Manghi
- Laboratoire de Physique Théorique (IRSAMC), Université de Toulouse, CNRS, UPS, 31 062 Toulouse, France
| | - Catherine Tardin
- Institut de Pharmacologie et de Biologie Structurale, Université de Toulouse, CNRS, UPS, 31 077 Toulouse, France
| |
Collapse
|
32
|
Jahan M, Uline MJ. Quantifying Mg 2+ Binding to ssDNA Oligomers: A Self-Consistent Field Theory Study at Varying Ionic Strengths and Grafting Densities. Polymers (Basel) 2018; 10:polym10121403. [PMID: 30961328 PMCID: PMC6401855 DOI: 10.3390/polym10121403] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2018] [Revised: 12/12/2018] [Accepted: 12/14/2018] [Indexed: 02/07/2023] Open
Abstract
The performance of aptamer-based biosensors is crucially impacted by their interactions with physiological metal ions, which can alter their structures and chemical properties. Therefore, elucidating the nature of these interactions carries the utmost importance in the robust design of highly efficient biosensors. We investigated Mg2+ binding to varying sequences of polymers to capture the effects of ionic strength and grafting density on ion binding and molecular reorganization of the polymer layer. The polymers are modeled as ssDNA aptamers using a self-consistent field theory, which accounts for non-covalent ion binding by integrating experimentally-derived binding constants. Our model captures the typical polyelectrolyte behavior of chain collapse with increased ionic strength for the ssDNA chains at low grafting density and exhibits the well-known re-entrant phenomena of stretched chains with increased ionic strength at high grafting density. The binding results suggest that electrostatic attraction between the monomers and Mg2+ plays the dominant role in defining the ion cloud around the ssDNA chains and generates a nearly-uniform ion distribution along the chains containing varying monomer sequences. These findings are in qualitative agreement with recent experimental results for Mg2+ binding to surface-bound ssDNA.
Collapse
Affiliation(s)
- Merina Jahan
- Department of Chemical Engineering, University of South Carolina, Columbia, SC 29208, USA.
| | - Mark J Uline
- Department of Chemical Engineering, University of South Carolina, Columbia, SC 29208, USA.
- Biomedical Engineering Program, University of South Carolina, Columbia, SC 29208, USA.
| |
Collapse
|
33
|
Yildirim A, Brenner N, Sutherland R, Feig M. Role of protein interactions in stabilizing canonical DNA features in simulations of DNA in crowded environments. BMC BIOPHYSICS 2018; 11:8. [PMID: 30555686 PMCID: PMC6286541 DOI: 10.1186/s13628-018-0048-y] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/08/2018] [Accepted: 11/21/2018] [Indexed: 11/28/2022]
Abstract
Background Cellular environments are highly crowded with biological macromolecules resulting in frequent non-specific interactions. While the effect of such crowding on protein structure and dynamics has been studied extensively, very little is known how cellular crowding affects the conformational sampling of nucleic acids. Results The effect of protein crowding on the conformational preferences of DNA (deoxyribonucleic acid) is described from fully atomistic molecular dynamics simulations of systems containing a DNA dodecamer surrounded by protein crowders. From the simulations, it was found that DNA structures prefer to stay in B-like conformations in the presence of the crowders. The preference for B-like conformations results from non-specific interactions of crowder proteins with the DNA sugar-phosphate backbone. Moreover, the simulations suggest that the crowder interactions narrow the conformational sampling to canonical regions of the conformational space. Conclusions The overall conclusion is that crowding effects may stabilize the canonical features of DNA that are most important for biological function. The results are complementary to a previous study of DNA in reduced dielectric environments where reduced dielectric environments alone led to a conformational shift towards A-DNA. Such a shift was not observed here suggested that the reduced dielectric response of cellular environments is counteracted by non-specific interactions with protein crowders under in vivo conditions. Electronic supplementary material The online version of this article (10.1186/s13628-018-0048-y) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Asli Yildirim
- 1Department of Chemistry, Michigan State University, East Lansing, MI 48824 USA
| | - Nathalie Brenner
- 2Faculty of Mathematics and Natural Sciences, University of Düsseldorf, 40225 Düsseldorf, Germany.,3Department of Biochemistry & Molecular Biology, Michigan State University, 603 Wilson Road, Room BCH 218, East Lansing, MI 48824 USA
| | - Robert Sutherland
- 3Department of Biochemistry & Molecular Biology, Michigan State University, 603 Wilson Road, Room BCH 218, East Lansing, MI 48824 USA
| | - Michael Feig
- 3Department of Biochemistry & Molecular Biology, Michigan State University, 603 Wilson Road, Room BCH 218, East Lansing, MI 48824 USA
| |
Collapse
|
34
|
Roh JH, Kilburn D, Behrouzi R, Sung W, Briber RM, Woodson SA. Effects of Preferential Counterion Interactions on the Specificity of RNA Folding. J Phys Chem Lett 2018; 9:5726-5732. [PMID: 30211556 PMCID: PMC6351067 DOI: 10.1021/acs.jpclett.8b02086] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
The real-time search for native RNA structure is essential for the operation of regulatory RNAs. We previously reported that a fraction of the Azoarcus ribozyme achieves a compact structure in less than a millisecond. To scrutinize the forces that drive initial folding steps, we used time-resolved SAXS to compare the folding dynamics of this ribozyme in thermodynamically isostable concentrations of different counterions. The results show that the size of the fast-folding population increases with the number of available counterions and correlates with the flexibility of initial RNA structures. Within 1 ms of folding, Mg2+ exhibits a smaller preferential interaction coefficient per charge, ΔΓ+/ Z, than Na+ or [Co(NH3)6]3+. The lower ΔΓ+/ Z corresponds to a smaller yield of folded RNA, although Mg2+ stabilizes native RNA more efficiently than other ions at equilibrium. These results suggest that strong Mg2+-RNA interactions impede the search for globally native structure during early folding stages.
Collapse
Affiliation(s)
- Joon Ho Roh
- Department of Materials Science and Engineering , University of Maryland , College Park , Maryland 20742 , United States
- T. C. Jenkins Department of Biophysics , Johns Hopkins University , Baltimore , Maryland 21218 , United States
| | - Duncan Kilburn
- T. C. Jenkins Department of Biophysics , Johns Hopkins University , Baltimore , Maryland 21218 , United States
| | - Reza Behrouzi
- Cell Biology , Harvard Medical School , Boston , Massachusetts 02115 , United States
| | - Wokyung Sung
- Department of Physics , Pohang University of Science and Technology , Pohang 37673 , Republic of Korea
| | - R M Briber
- Department of Materials Science and Engineering , University of Maryland , College Park , Maryland 20742 , United States
| | - Sarah A Woodson
- T. C. Jenkins Department of Biophysics , Johns Hopkins University , Baltimore , Maryland 21218 , United States
| |
Collapse
|
35
|
Ferapontova EE. DNA Electrochemistry and Electrochemical Sensors for Nucleic Acids. ANNUAL REVIEW OF ANALYTICAL CHEMISTRY (PALO ALTO, CALIF.) 2018; 11:197-218. [PMID: 29894229 DOI: 10.1146/annurev-anchem-061417-125811] [Citation(s) in RCA: 101] [Impact Index Per Article: 16.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/24/2023]
Abstract
Sensitive, specific, and fast analysis of nucleic acids (NAs) is strongly needed in medicine, environmental science, biodefence, and agriculture for the study of bacterial contamination of food and beverages and genetically modified organisms. Electrochemistry offers accurate, simple, inexpensive, and robust tools for the development of such analytical platforms that can successfully compete with other approaches for NA detection. Here, electrode reactions of DNA, basic principles of electrochemical NA analysis, and their relevance for practical applications are reviewed and critically discussed.
Collapse
Affiliation(s)
- Elena E Ferapontova
- Interdisciplinary Nanoscience Center (iNANO), Aarhus University, DK-8000 Aarhus C, Denmark;
| |
Collapse
|
36
|
Bonilla S, Limouse C, Bisaria N, Gebala M, Mabuchi H, Herschlag D. Single-Molecule Fluorescence Reveals Commonalities and Distinctions among Natural and in Vitro-Selected RNA Tertiary Motifs in a Multistep Folding Pathway. J Am Chem Soc 2017; 139:18576-18589. [PMID: 29185740 PMCID: PMC5748328 DOI: 10.1021/jacs.7b08870] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
![]()
Decades
of study of the RNA folding problem have revealed that
diverse and complex structured RNAs are built from a common set of
recurring structural motifs, leading to the perspective that a generalizable
model of RNA folding may be developed from understanding of the folding
properties of individual structural motifs. We used single-molecule
fluorescence to dissect the kinetic and thermodynamic properties of
a set of variants of a common tertiary structural motif, the tetraloop/tetraloop-receptor
(TL/TLR). Our results revealed a multistep TL/TLR folding pathway
in which preorganization of the ubiquitous AA-platform submotif precedes
the formation of the docking transition state and tertiary A-minor
hydrogen bond interactions form after the docking transition state.
Differences in ion dependences between TL/TLR variants indicated the
occurrence of sequence-dependent conformational rearrangements prior
to and after the formation of the docking transition state. Nevertheless,
varying the junction connecting the TL/TLR produced a common kinetic
and ionic effect for all variants, suggesting that the global conformational
search and compaction electrostatics are energetically independent
from the formation of the tertiary motif contacts. We also found that in vitro-selected variants, despite their similar stability
at high Mg2+ concentrations, are considerably less stable
than natural variants under near-physiological ionic conditions, and
the occurrence of the TL/TLR sequence variants in Nature correlates
with their thermodynamic stability in isolation. Overall, our findings
are consistent with modular but complex energetic properties of RNA
structural motifs and will aid in the eventual quantitative description
of RNA folding from its secondary and tertiary structural elements.
Collapse
Affiliation(s)
- Steve Bonilla
- Department of Chemical Engineering, ‡Department of Applied Physics, §Department of Biochemistry, ∥Department of Chemistry, ⊥Stanford ChEM-H, Stanford University , Stanford, California 94305, United States
| | - Charles Limouse
- Department of Chemical Engineering, ‡Department of Applied Physics, §Department of Biochemistry, ∥Department of Chemistry, ⊥Stanford ChEM-H, Stanford University , Stanford, California 94305, United States
| | - Namita Bisaria
- Department of Chemical Engineering, ‡Department of Applied Physics, §Department of Biochemistry, ∥Department of Chemistry, ⊥Stanford ChEM-H, Stanford University , Stanford, California 94305, United States
| | - Magdalena Gebala
- Department of Chemical Engineering, ‡Department of Applied Physics, §Department of Biochemistry, ∥Department of Chemistry, ⊥Stanford ChEM-H, Stanford University , Stanford, California 94305, United States
| | - Hideo Mabuchi
- Department of Chemical Engineering, ‡Department of Applied Physics, §Department of Biochemistry, ∥Department of Chemistry, ⊥Stanford ChEM-H, Stanford University , Stanford, California 94305, United States
| | - Daniel Herschlag
- Department of Chemical Engineering, ‡Department of Applied Physics, §Department of Biochemistry, ∥Department of Chemistry, ⊥Stanford ChEM-H, Stanford University , Stanford, California 94305, United States
| |
Collapse
|
37
|
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
|
38
|
Quantitative tests of a reconstitution model for RNA folding thermodynamics and kinetics. Proc Natl Acad Sci U S A 2017; 114:E7688-E7696. [PMID: 28839094 DOI: 10.1073/pnas.1703507114] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
Decades of study of the architecture and function of structured RNAs have led to the perspective that RNA tertiary structure is modular, made of locally stable domains that retain their structure across RNAs. We formalize a hypothesis inspired by this modularity-that RNA folding thermodynamics and kinetics can be quantitatively predicted from separable energetic contributions of the individual components of a complex RNA. This reconstitution hypothesis considers RNA tertiary folding in terms of ΔGalign, the probability of aligning tertiary contact partners, and ΔGtert, the favorable energetic contribution from the formation of tertiary contacts in an aligned state. This hypothesis predicts that changes in the alignment of tertiary contacts from different connecting helices and junctions (ΔGHJH) or from changes in the electrostatic environment (ΔG+/-) will not affect the energetic perturbation from a mutation in a tertiary contact (ΔΔGtert). Consistent with these predictions, single-molecule FRET measurements of folding of model RNAs revealed constant ΔΔGtert values for mutations in a tertiary contact embedded in different structural contexts and under different electrostatic conditions. The kinetic effects of these mutations provide further support for modular behavior of RNA elements and suggest that tertiary mutations may be used to identify rate-limiting steps and dissect folding and assembly pathways for complex RNAs. Overall, our model and results are foundational for a predictive understanding of RNA folding that will allow manipulation of RNA folding thermodynamics and kinetics. Conversely, the approaches herein can identify cases where an independent, additive model cannot be applied and so require additional investigation.
Collapse
|
39
|
Sun LZ, Kranawetter C, Heng X, Chen SJ. Predicting Ion Effects in an RNA Conformational Equilibrium. J Phys Chem B 2017; 121:8026-8036. [PMID: 28780864 DOI: 10.1021/acs.jpcb.7b03873] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
We develop a partial charge-based tightly bound ion (PCTBI) model for the ion effects in RNA folding. On the basis of the Monte Carlo tightly bound ion (MCTBI) approach, the model can account for ion fluctuation and correlation effects, and can predict the ion distribution around the RNA. Furthermore, unlike the previous coarse-grained RNA charge models, where negative charges are placed on the phosphates only, the current new model considers the detailed all-atom partial charge distribution on the RNA. Thus, the model not only keeps the advantage of the MCTBI model, but also has the potential to provide important detailed information unattainable by the previous MCTBI models. For example, the model predicts the reduction in ion binding upon protein binding and ion-induced conformational switches. For hepatitis C virus genomic RNA, the model predicts a Mg2+-induced stabilization of a kissing motif for a cis-acting regulatory element in the genomic RNA. Extensive theory-experiment comparisons support the reliability of the theoretical predictions. Therefore, the model may serve as a robust starting point for further development of an accurate method for ion effects in an RNA conformational equilibrium and RNA-cofactor interactions.
Collapse
Affiliation(s)
- Li-Zhen Sun
- Department of Physics, Department of Biochemistry, and Informatics Institute and ‡Department of Biochemistry, University of Missouri , Columbia, Missouri 65211, United States
| | - Clayton Kranawetter
- Department of Physics, Department of Biochemistry, and Informatics Institute and ‡Department of Biochemistry, University of Missouri , Columbia, Missouri 65211, United States
| | - Xiao Heng
- Department of Physics, Department of Biochemistry, and Informatics Institute and ‡Department of Biochemistry, University of Missouri , Columbia, Missouri 65211, United States
| | - Shi-Jie Chen
- Department of Physics, Department of Biochemistry, and Informatics Institute and ‡Department of Biochemistry, University of Missouri , Columbia, Missouri 65211, United States
| |
Collapse
|
40
|
Allred BE, Gebala M, Herschlag D. Determination of Ion Atmosphere Effects on the Nucleic Acid Electrostatic Potential and Ligand Association Using AH +·C Wobble Formation in Double-Stranded DNA. J Am Chem Soc 2017; 139:7540-7548. [PMID: 28489947 PMCID: PMC5466006 DOI: 10.1021/jacs.7b01830] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
![]()
The high charge density of nucleic
acids and resulting ion atmosphere
profoundly influence the conformational landscape of RNA and DNA and
their association with small molecules and proteins. Electrostatic
theories have been applied to quantitatively model the electrostatic
potential surrounding nucleic acids and the effects of the surrounding
ion atmosphere, but experimental measures of the potential and tests
of these models have often been complicated by conformational changes
and multisite binding equilibria, among other factors. We sought a
simple system to further test the basic predictions from electrostatics
theory and to measure the energetic consequences of the nucleic acid
electrostatic field. We turned to a DNA system developed by Bevilacqua
and co-workers that involves a proton as a ligand whose binding is
accompanied by formation of an internal AH+·C wobble
pair [Siegfried, N. A., et al. Biochemistry, 2010, 49, 3225]. Consistent with predictions
from polyelectrolyte models, we observed logarithmic dependences of
proton affinity versus salt concentration of −0.96 ± 0.03
and −0.52 ± 0.01 with monovalent and divalent cations,
respectively, and these results help clarify prior results that appeared
to conflict with these fundamental models. Strikingly, quantitation
of the ion atmosphere content indicates that divalent cations are
preferentially lost over monovalent cations upon A·C protonation,
providing experimental indication of the preferential localization
of more highly charged cations to the inner shell of the ion atmosphere.
The internal AH+·C wobble system further allowed us
to parse energetic contributions and extract estimates for the electrostatic
potential at the position of protonation. The results give a potential
near the DNA surface at 20 mM Mg2+ that is much less substantial
than at 20 mM K+ (−120 mV vs −210 mV). These
values and difference are similar to predictions from theory, and
the potential is substantially reduced at higher salt, also as predicted;
however, even at 1 M K+ the potential remains substantial,
counter to common assumptions. The A·C protonation module allows
extraction of new properties of the ion atmosphere and provides an
electrostatic meter that will allow local electrostatic potential
and energetics to be measured within nucleic acids and their complexes
with proteins.
Collapse
Affiliation(s)
- Benjamin E Allred
- Department of Biochemistry, Stanford University , Stanford, California 94305, United States
| | - Magdalena Gebala
- Department of Biochemistry, Stanford University , Stanford, California 94305, United States
| | - Daniel Herschlag
- Department of Biochemistry, Stanford University , Stanford, California 94305, United States.,Department of Chemistry, Stanford University , Stanford, California 94305, United States.,ChEM-H Institute, Stanford University , Stanford, California 94305, United States
| |
Collapse
|
41
|
Andrews CT, Campbell BA, Elcock AH. Direct Comparison of Amino Acid and Salt Interactions with Double-Stranded and Single-Stranded DNA from Explicit-Solvent Molecular Dynamics Simulations. J Chem Theory Comput 2017; 13:1794-1811. [PMID: 28288277 DOI: 10.1021/acs.jctc.6b00883] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Given the ubiquitous nature of protein-DNA interactions, it is important to understand the interaction thermodynamics of individual amino acid side chains for DNA. One way to assess these preferences is to perform molecular dynamics (MD) simulations. Here we report MD simulations of 20 amino acid side chain analogs interacting simultaneously with both a 70-base-pair double-stranded DNA and with a 70-nucleotide single-stranded DNA. The relative preferences of the amino acid side chains for dsDNA and ssDNA match well with values deduced from crystallographic analyses of protein-DNA complexes. The estimated apparent free energies of interaction for ssDNA, on the other hand, correlate well with previous simulation values reported for interactions with isolated nucleobases, and with experimental values reported for interactions with guanosine. Comparisons of the interactions with dsDNA and ssDNA indicate that, with the exception of the positively charged side chains, all types of amino acid side chain interact more favorably with ssDNA, with intercalation of aromatic and aliphatic side chains being especially notable. Analysis of the data on a base-by-base basis indicates that positively charged side chains, as well as sodium ions, preferentially bind to cytosine in ssDNA, and that negatively charged side chains, and chloride ions, preferentially bind to guanine in ssDNA. These latter observations provide a novel explanation for the lower salt dependence of DNA duplex stability in GC-rich sequences relative to AT-rich sequences.
Collapse
Affiliation(s)
- Casey T Andrews
- Department of Biochemistry, University of Iowa , Iowa City, Iowa 52242, United States
| | - Brady A Campbell
- Department of Biochemistry, University of Iowa , Iowa City, Iowa 52242, United States
| | - Adrian H Elcock
- Department of Biochemistry, University of Iowa , Iowa City, Iowa 52242, United States
| |
Collapse
|
42
|
Kékedy-Nagy L, Ferapontova EE, Brand I. Submolecular Structure and Orientation of Oligonucleotide Duplexes Tethered to Gold Electrodes Probed by Infrared Reflection Absorption Spectroscopy: Effect of the Electrode Potentials. J Phys Chem B 2017; 121:1552-1565. [PMID: 28177253 DOI: 10.1021/acs.jpcb.6b12363] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
Unique electronic and ligand recognition properties of the DNA double helix provide basis for DNA applications in biomolecular electronic and biosensor devices. However, the relation between the structure of DNA at electrified interfaces and its electronic properties is still not well understood. Here, potential-driven changes in the submolecular structure of DNA double helices composed of either adenine-thymine (dAdT)25 or cytosine-guanine (dGdC)20 base pairs tethered to the gold electrodes are for the first time analyzed by in situ polarization modulation infrared reflection absorption spectroscopy (PM IRRAS) performed under the electrochemical control. It is shown that the conformation of the DNA duplexes tethered to gold electrodes via the C6 alkanethiol linker strongly depends on the nucleic acid sequence composition. The tilt of purine and pyrimidine rings of the complementary base pairs (dAdT and dGdC) depends on the potential applied to the electrode. By contrast, neither the conformation nor orientation of the ionic in character phosphate-sugar backbone is affected by the electrode potentials. At potentials more positive than the potential of zero charge (pzc), a gradual tilting of the double helix is observed. In this tilted orientation, the planes of the complementary purine and pyrimidine rings lie ideally parallel to each other. These potentials do not affect the integral stability of the DNA double helix at the charged interface. At potentials more negative than the pzc, DNA helices adopt a vertical to the gold surface orientation. Tilt of the purine and pyrimidine rings depends on the composition of the double helix. In monolayers composed of (dAdT)25 molecules the rings of the complementary base pairs lie parallel to each other. By contrast, the tilt of purine and pyrimidine rings in (dGdC)20 helices depends on the potential applied to the electrode. Such potential-induced mobility of the complementary base pairs can destabilize the helix structure at a submolecular level. These pioneer results on the potential-driven changes in the submolecular structure of double stranded DNA adsorbed on conductive supports contribute to further understanding of the potential-driven sequence-specific electronic properties of surface-tethered oligonucleotides.
Collapse
Affiliation(s)
- László Kékedy-Nagy
- Interdisciplinary Nanoscience Center (iNANO) and Center for DNA Nanotechnology (CDNA), Science and Technology, Aarhus University , Gustav Wieds Vej 14, DK-8000 Aarhus-C, Denmark
| | - Elena E Ferapontova
- Interdisciplinary Nanoscience Center (iNANO) and Center for DNA Nanotechnology (CDNA), Science and Technology, Aarhus University , Gustav Wieds Vej 14, DK-8000 Aarhus-C, Denmark
| | - Izabella Brand
- Department of Chemistry, University of Oldenburg , 26111 Oldenburg, Germany
| |
Collapse
|