1
|
Vanness BC, Linz TH. Multiplexed miRNA and Protein Analysis Using Digital Quantitative PCR in Microwell Arrays. Anal Chem 2024; 96:1371-1379. [PMID: 38183281 PMCID: PMC11168192 DOI: 10.1021/acs.analchem.3c05213] [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: 01/08/2024]
Abstract
Proteins and microRNAs (miRNAs) act in tandem within biological pathways to regulate cellular functions, and their misregulation has been correlated to numerous diseases. Because of their interconnectedness, both miRNAs and proteins must be evaluated together to obtain accurate insights into the molecular pathways of pathogenesis. However, few analytical techniques can measure both classes of biomolecules in parallel from a single biological sample. Here, microfluidic digital quantitative PCR (dqPCR) was developed to simultaneously quantify miRNA and protein targets in a multiplexed assay using a single detection chemistry. This streamlined analysis was achieved by integrating base-stacking PCR and immuno-PCR in a microfluidic array platform. Analyses of let-7a (miRNA) and IL-6 (protein) were first optimized separately to identify thermocycling and capture conditions amenable to both biomolecules. Singleplex dqPCR studies exhibited the expected digital signals and quantification cycles for both analytes over a range of concentrations. Multiplexed analyses were then conducted to quantify both let-7a and IL-6 with high sensitivity (LODs ∼ 3 fM) over a broad dynamic range (5-5000 fM) using only standard PCR reagents. This multiplexed dqPCR was then translated to the analysis of HEK293 cell lysate, where endogenous let-7a and IL-6 were measured simultaneously without sample purification or pretreatment. Collectively, these studies demonstrate that the integration of BS-PCR and immuno-PCR achieves a sensitive and streamlined approach for multiplexed analyses of miRNAs and proteins, which will enable researchers to gain better insights into disease pathogenesis in future applications.
Collapse
Affiliation(s)
- Brice C. Vanness
- Department of Chemistry, Wayne State University, Detroit, MI 48202
| | - Thomas H. Linz
- Department of Chemistry, Wayne State University, Detroit, MI 48202
| |
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
|
Three-Way DNA Junction as an End Label for DNA in Atomic Force Microscopy Studies. Int J Mol Sci 2022; 23:ijms231911404. [PMID: 36232705 PMCID: PMC9569629 DOI: 10.3390/ijms231911404] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2022] [Revised: 09/17/2022] [Accepted: 09/23/2022] [Indexed: 11/25/2022] Open
Abstract
Atomic Force Microscopy (AFM) is widely used for topographic imaging of DNA and protein-DNA complexes in ambient conditions with nanometer resolution. In AFM studies of protein-DNA complexes, identifying the protein’s location on the DNA substrate is one of the major goals. Such studies require distinguishing between the DNA ends, which can be accomplished by end-specific labeling of the DNA substrate. We selected as labels three-way DNA junctions (3WJ) assembled from synthetic DNA oligonucleotides with two arms of 39–40 bp each. The third arm has a three-nucleotide overhang, GCT, which is paired with the sticky end of the DNA substrate generated by the SapI enzyme. Ligation of the 3WJ results in the formation of a Y-type structure at the end of the linear DNA mole cule, which is routinely identified in the AFM images. The yield of labeling is 69%. The relative orientation of arms in the Y-end varies, such dynamics were directly visualized with time-lapse AFM studies using high-speed AFM (HS-AFM). This labeling approach was applied to the characterization of the nucleosome arrays assembled on different DNA templates. HS-AFM experiments revealed a high dynamic of nucleosomes resulting in a spontaneous unraveling followed by disassembly of nucleosomes.
Collapse
|
6
|
Dutkiewicz Z. Computational methods for calculation of protein-ligand binding affinities in structure-based drug design. PHYSICAL SCIENCES REVIEWS 2022. [DOI: 10.1515/psr-2020-0034] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
Abstract
Drug design is an expensive and time-consuming process. Any method that allows reducing the time the costs of the drug development project can have great practical value for the pharmaceutical industry. In structure-based drug design, affinity prediction methods are of great importance. The majority of methods used to predict binding free energy in protein-ligand complexes use molecular mechanics methods. However, many limitations of these methods in describing interactions exist. An attempt to go beyond these limits is the application of quantum-mechanical description for all or only part of the analyzed system. However, the extensive use of quantum mechanical (QM) approaches in drug discovery is still a demanding challenge. This chapter briefly reviews selected methods used to calculate protein-ligand binding affinity applied in virtual screening (VS), rescoring of docked poses, and lead optimization stage, including QM methods based on molecular simulations.
Collapse
Affiliation(s)
- Zbigniew Dutkiewicz
- Department of Chemical Technology of Drugs , Poznan University of Medical Sciences , ul. Grunwaldzka 6 , 60-780 Poznań , Poznan , 60-780, Poland
| |
Collapse
|
7
|
Szczuka A, Horton J, Evans KJ, DiPietri VT, Sivey JD, Wigginton KR. Chloride Enhances DNA Reactivity with Chlorine under Conditions Relevant to Water Treatment. ENVIRONMENTAL SCIENCE & TECHNOLOGY 2022; 56:13347-13356. [PMID: 36027047 PMCID: PMC9494735 DOI: 10.1021/acs.est.2c03267] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/10/2022] [Revised: 08/09/2022] [Accepted: 08/09/2022] [Indexed: 06/15/2023]
Abstract
Free available chlorine (FAC) is widely used to inactivate viruses by oxidizing viral components, including genomes. It is commonly assumed that hypochlorous acid (HOCl) is the chlorinating agent responsible for virus inactivation; however, recent studies have underscored that minor constituents of FAC existing in equilibrium with HOCl, such as molecular chlorine (Cl2), can influence FAC reactivity toward select organic compounds. This study measures the FAC reaction kinetics with dsDNA and ssDNA extracted from representative bacteriophages (T3 and ϕX174) in samples augmented with chloride. Herein, chloride enhances FAC reactivity toward dsDNA and, to a lesser extent, toward ssDNA, especially at pH < 7.5. The enhanced reactivity can be attributed to the formation of Cl2. Second-order rate constants were determined for reactions of ssDNA and dsDNA with HOCl and Cl2. DNA chlorination kinetics followed the reactivity-selectivity principle, where the more-reactive nucleophilic species (ssDNA, ∼100× more reactive than dsDNA) reacted less selectively with electrophilic FAC species. The addition of chloride was also shown to enhance the inactivation of bacteriophage T3 (dsDNA genome) by FAC but did not enhance the inactivation of bacteriophage ϕX174 (ssDNA genome). Overall, the results suggest that Cl2 is an important chlorinating agent of nucleic acids and viruses.
Collapse
Affiliation(s)
- Aleksandra Szczuka
- Department
of Civil and Environmental Engineering, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Jordon Horton
- Department
of Civil and Environmental Engineering, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Kelsey J. Evans
- Department
of Chemistry, Towson University, Towson, Maryland 21252, United States
| | - Vincent T. DiPietri
- Department
of Chemistry, Towson University, Towson, Maryland 21252, United States
| | - John D. Sivey
- Department
of Chemistry, Towson University, Towson, Maryland 21252, United States
| | - Krista R. Wigginton
- Department
of Civil and Environmental Engineering, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
8
|
Vahid H, Scacchi A, Yang X, Ala-Nissila T, Sammalkorpi M. Modified Poisson–Boltzmann theory for polyelectrolytes in monovalent salt solutions with finite-size ions. J Chem Phys 2022; 156:214906. [DOI: 10.1063/5.0092273] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a soft-potential-enhanced Poisson–Boltzmann (SPB) theory to efficiently capture ion distributions and electrostatic potential around rodlike charged macromolecules. The SPB model is calibrated with a coarse-grained particle-based model for polyelectrolytes (PEs) in monovalent salt solutions as well as compared to a full atomistic molecular dynamics simulation with the explicit solvent. We demonstrate that our modification enables the SPB theory to accurately predict monovalent ion distributions around a rodlike PE in a wide range of ion and charge distribution conditions in the weak-coupling regime. These include excess salt concentrations up to 1M and ion sizes ranging from small ions, such as Na+ or Cl−, to softer and larger ions with a size comparable to the PE diameter. The work provides a simple way to implement an enhancement that effectively captures the influence of ion size and species into the PB theory in the context of PEs in aqueous salt solutions.
Collapse
Affiliation(s)
- Hossein Vahid
- Department of Chemistry and Materials Science, Aalto University, P.O. Box 16100, FI-00076 Aalto, Finland
- Department of Applied Physics, Aalto University, P.O. Box 11000, FI-00076 Aalto, Finland
| | - Alberto Scacchi
- Department of Chemistry and Materials Science, Aalto University, P.O. Box 16100, FI-00076 Aalto, Finland
- Department of Applied Physics, Aalto University, P.O. Box 11000, FI-00076 Aalto, Finland
| | - Xiang Yang
- Department of Applied Physics, Aalto University, P.O. Box 11000, FI-00076 Aalto, Finland
- Quantum Technology Finland Center of Excellence, Department of Applied Physics, Aalto University, P.O. Box 11000, FI-00076 Aalto, Finland
| | - Tapio Ala-Nissila
- Quantum Technology Finland Center of Excellence, Department of Applied Physics, Aalto University, P.O. Box 11000, FI-00076 Aalto, Finland
- Interdisciplinary Centre for Mathematical Modelling and Department of Mathematical Sciences, Loughborough University, Loughborough, Leicestershire LE11 3TU, United Kingdom
| | - Maria Sammalkorpi
- Department of Chemistry and Materials Science, Aalto University, P.O. Box 16100, FI-00076 Aalto, Finland
- Department of Bioproducts and Biosystems, Aalto University, P.O. Box 16100, FI-00076 Aalto, Finland
- Academy of Finland Center of Excellence in Life-Inspired Hybrid Materials (LIBER), Aalto University, P.O. Box 16100, FI-00076 Aalto, Finland
| |
Collapse
|
9
|
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
|
10
|
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
|
11
|
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
|
12
|
Kolesnikov ES, Gushchin IY, Zhilyaev PA, Onufriev AV. Similarities and Differences between Na + and K + Distributions around DNA Obtained with Three Popular Water Models. J Chem Theory Comput 2021; 17:7246-7259. [PMID: 34633813 DOI: 10.1021/acs.jctc.1c00332] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We have compared distributions of sodium and potassium ions around double-stranded DNA, simulated using fixed charge SPC/E, TIP3P, and OPC water models and the Joung/Cheatham (J/C) ion parameter set, as well as the Li/Merz HFE 6-12 (L/M HFE) ion parameters for OPC water. In all the simulations, the ion distributions are in qualitative agreement with Manning's condensation theory and the Debye-Hückel theory, where expected. In agreement with experiment, binding affinity of monovalent ions to DNA does not depend on ion type in every solvent model. However, behavior of deeply bound ions, including ions bound to specific sites, depends strongly on the solvent model. In particular, the number of potassium ions in the minor groove of AT-tracts differs at least 3-fold between the solvent models tested. The number of sodium ions associated with the DNA agrees quantitatively with the experiment for the OPC water model, followed closely by TIP3P+J/C; the largest deviation from the experiment, ∼10%, is seen for SPC/E+J/C. On the other hand, SPC/E+J/C model is most consistent (67%) with the experimental potassium binding sites, followed by OPC+J/C (60%), TIP3P+J/C (53%), and OPC+L/M HFE (27%). The use of NBFIX correction with TIP3P+J/C improves its consistency with the experiment. In summary, the choice of the solvent model matters little for simulating the diffuse atmosphere of sodium and potassium ions around DNA, but ion distributions become increasingly sensitive to the solvent model near the helical axis. We offer an explanation for these trends. There is no single gold standard solvent model, although OPC water with J/C ions or TIP3P with J/C + NBFIX may offer an imperfect compromise for practical simulations of ionic atmospheres around DNA.
Collapse
Affiliation(s)
- Egor S Kolesnikov
- Research Center for Molecular Mechanisms of Aging and Age-Related Diseases, Moscow Institute of Physics and Technology, Dolgoprudny 141700, Russia
| | - Ivan Yu Gushchin
- Research Center for Molecular Mechanisms of Aging and Age-Related Diseases, Moscow Institute of Physics and Technology, Dolgoprudny 141700, Russia
| | - Petr A Zhilyaev
- Center for Design, Manufacturing and Materials, Skolkovo Institute of Science and Technology, Bolshoy Boulevard 30, bld. 1, Moscow 121205, Russia
| | - Alexey V Onufriev
- Department of Computer Science, Virginia Tech, Blacksburg 24061-0131, United States.,Department of Physics, Virginia Tech, Blacksburg 24061-0131, United States.,Center for Soft Matter and Biological Physics, Virginia Tech, Blacksburg 24061-0131, United States
| |
Collapse
|
13
|
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
|
14
|
Largy E, König A, Ghosh A, Ghosh D, Benabou S, Rosu F, Gabelica V. Mass Spectrometry of Nucleic Acid Noncovalent Complexes. Chem Rev 2021; 122:7720-7839. [PMID: 34587741 DOI: 10.1021/acs.chemrev.1c00386] [Citation(s) in RCA: 34] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Nucleic acids have been among the first targets for antitumor drugs and antibiotics. With the unveiling of new biological roles in regulation of gene expression, specific DNA and RNA structures have become very attractive targets, especially when the corresponding proteins are undruggable. Biophysical assays to assess target structure as well as ligand binding stoichiometry, affinity, specificity, and binding modes are part of the drug development process. Mass spectrometry offers unique advantages as a biophysical method owing to its ability to distinguish each stoichiometry present in a mixture. In addition, advanced mass spectrometry approaches (reactive probing, fragmentation techniques, ion mobility spectrometry, ion spectroscopy) provide more detailed information on the complexes. Here, we review the fundamentals of mass spectrometry and all its particularities when studying noncovalent nucleic acid structures, and then review what has been learned thanks to mass spectrometry on nucleic acid structures, self-assemblies (e.g., duplexes or G-quadruplexes), and their complexes with ligands.
Collapse
Affiliation(s)
- Eric Largy
- Univ. Bordeaux, CNRS, INSERM, ARNA, UMR 5320, U1212, IECB, F-33600 Pessac, France
| | - Alexander König
- Univ. Bordeaux, CNRS, INSERM, ARNA, UMR 5320, U1212, IECB, F-33600 Pessac, France
| | - Anirban Ghosh
- Univ. Bordeaux, CNRS, INSERM, ARNA, UMR 5320, U1212, IECB, F-33600 Pessac, France
| | - Debasmita Ghosh
- Univ. Bordeaux, CNRS, INSERM, ARNA, UMR 5320, U1212, IECB, F-33600 Pessac, France
| | - Sanae Benabou
- Univ. Bordeaux, CNRS, INSERM, ARNA, UMR 5320, U1212, IECB, F-33600 Pessac, France
| | - Frédéric Rosu
- Univ. Bordeaux, CNRS, INSERM, IECB, UMS 3033, F-33600 Pessac, France
| | - Valérie Gabelica
- Univ. Bordeaux, CNRS, INSERM, ARNA, UMR 5320, U1212, IECB, F-33600 Pessac, France
| |
Collapse
|
15
|
Complex Conformational Dynamics of the Heart Failure-Associated Pre-miRNA-377 Hairpin Revealed by Single-Molecule Optical Tweezers. Int J Mol Sci 2021; 22:ijms22169008. [PMID: 34445712 PMCID: PMC8396532 DOI: 10.3390/ijms22169008] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2021] [Revised: 08/13/2021] [Accepted: 08/18/2021] [Indexed: 11/17/2022] Open
Abstract
Pre-miRNA-377 is a hairpin-shaped regulatory RNA associated with heart failure. Here, we use single-molecule optical tweezers to unzip pre-miRNA-377 and study its stability and dynamics. We show that magnesium ions have a strong stabilizing effect, and that sodium ions stabilize the hairpin more than potassium ions. The hairpin unfolds in a single step, regardless of buffer composition. Interestingly, hairpin folding occurs either in a single step (type 1) or through the formation of intermediates, in multiple steps (type 2) or gradually (type 3). Type 3 occurs only in the presence of both sodium and magnesium, while type 1 and 2 take place in all buffers, with type 1 being the most prevalent. By reducing the size of the native hairpin loop from fourteen to four nucleotides, we demonstrate that the folding heterogeneity originates from the large size of the hairpin loop. Further, while efficient pre-miRNA-377 binders are lacking, we demonstrate that the recently developed C2 ligand displays bimodal activity: it enhances the mechanical stability of the pre-miRNA-377 hairpin and perturbs its folding. The knowledge regarding pre-miRNA stability and dynamics that we provide is important in understanding its regulatory function and how it can be modulated to achieve a therapeutic effect, e.g., in heart failure treatment.
Collapse
|
16
|
Hadži S, Lah J. Origin of heat capacity increment in DNA folding: The hydration effect. Biochim Biophys Acta Gen Subj 2020; 1865:129774. [PMID: 33164852 DOI: 10.1016/j.bbagen.2020.129774] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2020] [Revised: 10/09/2020] [Accepted: 10/20/2020] [Indexed: 02/03/2023]
Abstract
BACKGROUND Understanding DNA folding thermodynamics is crucial for prediction of DNA thermal stability. It is now well established that DNA folding is accompanied by a decrease of the heat capacity ∆cp, F, however its molecular origin is not understood. In analogy to protein folding it has been assumed that this is due to dehydration of DNA constituents, however no evidence exists to support this conclusion. METHODS Here we analyze partial molar heat capacity of nucleic bases and nucleosides in aqueous solutions obtained from calorimetric experiments and calculate the hydration heat capacity contribution ∆cphyd. RESULTS We present hydration heat capacity contributions of DNA constituents and show that they correlate with the solvent accessible surface area. The average contribution for nucleic base dehydration is +0.56 J mol-1 K-1 Å-2 and can be used to estimate the ∆cp, F contribution for DNA folding. CONCLUSIONS We show that dehydration is one of the major sources contributing to the observed ∆cp, F increment in DNA folding. Other possible sources contributing to the overall ∆cp, F should be significant but appear to compensate each other to high degree. The calculated ∆cphyd for duplexes and noncanonical DNA structures agree excellently with the overall experimental ∆cp, F values. By contrast, empirical parametrizations developed for proteins result in poor ∆cphyd predictions and should not be applied to DNA folding. GENERAL SIGNIFICANCE Heat capacity is one of the main thermodynamic quantities that strongly affects thermal stability of macromolecules. At the molecular level the heat capacity in DNA folding stems from removal of water from nucleobases.
Collapse
Affiliation(s)
- S Hadži
- Department of Physical Chemistry, Faculty of Chemistry and Chemical Technology, University of Ljubljana, 1000 Ljubljana, Slovenia.
| | - J Lah
- Department of Physical Chemistry, Faculty of Chemistry and Chemical Technology, University of Ljubljana, 1000 Ljubljana, Slovenia.
| |
Collapse
|
17
|
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
|
18
|
Nepal S, Kumar P. Growth, Cell Division, and Gene Expression of Escherichia coli at Elevated Concentrations of Magnesium Sulfate: Implications for Habitability of Europa and Mars. Microorganisms 2020; 8:E637. [PMID: 32349403 PMCID: PMC7285182 DOI: 10.3390/microorganisms8050637] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2020] [Revised: 04/15/2020] [Accepted: 04/21/2020] [Indexed: 01/20/2023] Open
Abstract
We perform quantitative studies of the growth, death, and gene expression of Escherichia coli in a wide range of magnesium sulfate (MgSO 4 ) concentrations (0-2.5 M). Elevated concentration of MgSO 4 causes the inhibition of cell growth, leading to an increase in the population doubling time. We find that cells exhibit three distinct morphological phenotypes-(i) normal, (ii) filamentous, and (iii) small cells at 1 . 25 M MgSO 4 . Filamentous cells arise due to the lack of cell division, while the small cells arise due to the partial plasmolysis of the cells. We further find that cell death starts for salt concentrations >1.25 M and increases with an increasing concentration of MgSO 4 . For salt concentrations ≥1.66 M, the growth of cells stops and all the cells become smaller than the control cells, suggesting the plasmolysis of the population. Cells grown at salt concentration up to 2 . 07 M are reversible in both the growth rate and morphology upon the removal of the salt stress. The time scale of reversibility increases with increasing salt concentration. Finally, we investigate the expression of an osmotically inducible gene (osmC), genes involved in magnesium transport (corA), sulfate transport (cysP), and osmotically driven transport of water (aqpZ). We find that a high concentration of magnesium sulfate leads to the upregulation of cysP and osmC.
Collapse
Affiliation(s)
- Sudip Nepal
- Department of Physics, University of Arkansas, Fayetteville, AR 72701, USA;
- Microelectronics and Photonics Graduate Program, University of Arkansas, Fayetteville, AR 72701, USA
| | - Pradeep Kumar
- Department of Physics, University of Arkansas, Fayetteville, AR 72701, USA;
- Microelectronics and Photonics Graduate Program, University of Arkansas, Fayetteville, AR 72701, USA
| |
Collapse
|
19
|
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
|
20
|
Galindo-Murillo R, Cheatham TE. Lessons learned in atomistic simulation of double-stranded DNA: Solvation and salt concerns [Article v1.0]. ACTA ACUST UNITED AC 2019; 1. [PMID: 33073182 DOI: 10.33011/livecoms.1.2.9974] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Nucleic acids are highly charged macromolecules sensitive to their surroundings of water, salt, and other biomolecules. Molecular dynamics simulations with accurate biomolecular force fields provide a detailed atomistic view into DNA and RNA that has been useful to study the structure and dynamics of these molecules and their biological relevance. In this work we study the Drew-Dickerson dodecamer duplex with the sequence d(GCGCAATTGCGC)2 in three different salt concentrations and using different monvalent salt types to detect possible structural influence. Overall, the DNA shows no major structural changes regardless of amount or type of monovalent ions used. Our results show that only at very high salt conditions (5M) is a small structural effect observed in the DNA duplex, which mainly consist of narrowing of the grooves due to increased residence of ions. We also present the importance of sampling time to achieve a converged ensemble, which is of major relevance in any simulation to avoid biased or non-meaningful results.
Collapse
Affiliation(s)
- Rodrigo Galindo-Murillo
- Department of Medicinal Chemistry, L. S. Skaggs Pharmacy Institute, University of Utah, Salt Lake City, UT 84112
| | - Thomas E Cheatham
- Department of Medicinal Chemistry, L. S. Skaggs Pharmacy Institute, University of Utah, Salt Lake City, UT 84112
| |
Collapse
|
21
|
Denesyuk NA, Hori N, Thirumalai D. Molecular Simulations of Ion Effects on the Thermodynamics of RNA Folding. J Phys Chem B 2018; 122:11860-11867. [PMID: 30468380 DOI: 10.1021/acs.jpcb.8b08142] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
How ions affect RNA folding thermodynamics and kinetics is an important but a vexing problem that remains unsolved. Experiments have shown that the free-energy change, Δ G( c), of RNA upon folding varies with the salt concentration ( c) as, Δ G( c) = k c ln c + const, where the coefficient k c is proportional to the difference in the ion preferential coefficient, ΔΓ. We performed simulations of a coarse-grained model, by modeling electrostatic interactions implicitly and with explicit representation of ions, to elucidate the molecular underpinnings of the relationship between Δ G and ΔΓ. The simulations quantitatively reproduce the heat capacity for a pseudoknot, thus validating the model. We show that Δ G( c), calculated directly from ΔΓ, varies linearly with ln c ( c < 0.2 M), for a hairpin and the pseudoknot, demonstrating a molecular link between the two quantities. Explicit ion simulations also show the linear dependence of Δ G( c) on ln c at all c with k c = 2 kB T, except that Δ G( c) values are shifted by ∼2 kcal/mol higher than experiments. The discrepancy is due to an underestimation of Γ for both the folded and unfolded states while giving accurate values for ΔΓ. The predictions for the salt dependence of ΔΓ are amenable to test using single-molecule pulling experiments. The framework provided here can be used to obtain accurate thermodynamics for other RNA molecules as well.
Collapse
Affiliation(s)
- Natalia A Denesyuk
- Department of Chemistry and Biochemistry and Biophysics Program, Institute for Physical Science and Technology , University of Maryland , College Park , Maryland 20742 , United States
| | - Naoto Hori
- Department of Chemistry , University of Texas at Austin , Austin , Texas 78712 , United States
| | - D Thirumalai
- Department of Chemistry , University of Texas at Austin , Austin , Texas 78712 , United States
| |
Collapse
|
22
|
The 3-D structure of VNG0258H/RosR - A haloarchaeal DNA-binding protein in its ionic shell. J Struct Biol 2018; 204:191-198. [PMID: 30110657 DOI: 10.1016/j.jsb.2018.08.008] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2018] [Revised: 08/05/2018] [Accepted: 08/09/2018] [Indexed: 11/21/2022]
Abstract
Protein-DNA interactions are highly dependent on salt concentration. To gain insight into how such interactions are maintained in the highly saline cytoplasm of halophilic archaea, we determined the 3-D structure of VNG0258H/RosR, the first haloarchaeal DNA-binding protein from the extreme halophilic archaeon Halobactrium salinarum. It is a dimeric winged-helix-turn-helix (wHTH) protein with unique features due to adaptation to the halophilic environment. As ions are major players in DNA binding processes, particularly in halophilic environments, we investigated the solution structure of the ionic envelope and located anions in the first shell around the protein in the crystal using anomalous scattering. Anions that were found to be tightly bound to residues in the positively charged DNA-binding site would probably be released upon DNA binding and will thus make significant contribution to the driving force of the binding process. Unexpectedly, ions were also found in a buried internal cavity connected to the external medium by a tunnel. Our structure lays a solid groundwork for future structural, computational and biochemical studies on complexes of the protein with cognate DNA sequences, with implications to protein-DNA interactions in hyper-saline environments.
Collapse
|
23
|
Ryde U, Söderhjelm P. Ligand-Binding Affinity Estimates Supported by Quantum-Mechanical Methods. Chem Rev 2016; 116:5520-66. [DOI: 10.1021/acs.chemrev.5b00630] [Citation(s) in RCA: 175] [Impact Index Per Article: 21.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ulf Ryde
- Department of Theoretical
Chemistry and ‡Department of Biophysical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Pär Söderhjelm
- Department of Theoretical
Chemistry and ‡Department of Biophysical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| |
Collapse
|
24
|
Gebala M, Giambasu GM, Lipfert J, Bisaria N, Bonilla S, Li G, York DM, Herschlag D. Cation-Anion Interactions within the Nucleic Acid Ion Atmosphere Revealed by Ion Counting. J Am Chem Soc 2015; 137:14705-15. [PMID: 26517731 PMCID: PMC4739826 DOI: 10.1021/jacs.5b08395] [Citation(s) in RCA: 58] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The ion atmosphere is a critical structural, dynamic, and energetic component of nucleic acids that profoundly affects their interactions with proteins and ligands. Experimental methods that "count" the number of ions thermodynamically associated with the ion atmosphere allow dissection of energetic properties of the ion atmosphere, and thus provide direct comparison to theoretical results. Previous experiments have focused primarily on the cations that are attracted to nucleic acid polyanions, but have also showed that anions are excluded from the ion atmosphere. Herein, we have systematically explored the properties of anion exclusion, testing the zeroth-order model that anions of different identity are equally excluded due to electrostatic repulsion. Using a series of monovalent salts, we find, surprisingly, that the extent of anion exclusion and cation inclusion significantly depends on salt identity. The differences are prominent at higher concentrations and mirror trends in mean activity coefficients of the electrolyte solutions. Salts with lower activity coefficients exhibit greater accumulation of both cations and anions within the ion atmosphere, strongly suggesting that cation-anion correlation effects are present in the ion atmosphere and need to be accounted for to understand electrostatic interactions of nucleic acids. To test whether the effects of cation-anion correlations extend to nucleic acid kinetics and thermodynamics, we followed the folding of P4-P6, a domain of the Tetrahymena group I ribozyme, via single-molecule fluorescence resonance energy transfer in solutions with different salts. Solutions of identical concentration but lower activity gave slower and less favorable folding. Our results reveal hitherto unknown properties of the ion atmosphere and suggest possible roles of oriented ion pairs or anion-bridged cations in the ion atmosphere for electrolyte solutions of salts with reduced activity. Consideration of these new results leads to a reevaluation of the strengths and limitations of Poisson-Boltzmann theory and highlights the need for next-generation atomic-level models of the ion atmosphere.
Collapse
Affiliation(s)
- Magdalena Gebala
- Department of Biochemistry, Stanford University, Stanford, California 94305, United States
| | - George M. Giambasu
- BioMaPS Institute for Quantitative Biology and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, United States
| | - Jan Lipfert
- Department of Physics, Nanosystems Initiative Munich, and Center for Nanoscience, Ludwig Maximilian University of Munich, 80799 Munich, Germany
| | - Namita Bisaria
- Department of Biochemistry, Stanford University, Stanford, California 94305, United States
| | - Steve Bonilla
- Department of Chemical Engineering, Stanford University, Stanford, California 94305, United States
| | - Guangchao Li
- School of Earth, Energy and Environment Sciences, Stanford University, Stanford, California 94305, United States
| | - Darrin M. York
- BioMaPS Institute for Quantitative Biology and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, United States
| | - Daniel Herschlag
- Department of Biochemistry, Stanford University, Stanford, California 94305, United States
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| |
Collapse
|
25
|
Siddique JA, Sharma S, Naqvi S, Ahmad A, Khatoon A, Mohd-Setapar SH. Apparent Molal Volume and Compressibility of Glucose and Maltose at Different Temperatures in Lysozyme Solution. ARABIAN JOURNAL FOR SCIENCE AND ENGINEERING 2015. [DOI: 10.1007/s13369-015-1818-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
26
|
Giambaşu GM, Gebala MK, Panteva MT, Luchko T, Case DA, York DM. Competitive interaction of monovalent cations with DNA from 3D-RISM. Nucleic Acids Res 2015; 43:8405-15. [PMID: 26304542 PMCID: PMC4787805 DOI: 10.1093/nar/gkv830] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2015] [Revised: 08/05/2015] [Accepted: 08/06/2015] [Indexed: 12/15/2022] Open
Abstract
The composition of the ion atmosphere surrounding nucleic acids affects their folding, condensation and binding to other molecules. It is thus of fundamental importance to gain predictive insight into the formation of the ion atmosphere and thermodynamic consequences when varying ionic conditions. An early step toward this goal is to benchmark computational models against quantitative experimental measurements. Herein, we test the ability of the three dimensional reference interaction site model (3D-RISM) to reproduce preferential interaction parameters determined from ion counting (IC) experiments for mixed alkali chlorides and dsDNA. Calculations agree well with experiment with slight deviations for salt concentrations >200 mM and capture the observed trend where the extent of cation accumulation around the DNA varies inversely with its ionic size. Ion distributions indicate that the smaller, more competitive cations accumulate to a greater extent near the phosphoryl groups, penetrating deeper into the grooves. In accord with experiment, calculated IC profiles do not vary with sequence, although the predicted ion distributions in the grooves are sequence and ion size dependent. Calculations on other nucleic acid conformations predict that the variation in linear charge density has a minor effect on the extent of cation competition.
Collapse
Affiliation(s)
- George M Giambaşu
- BioMaPS Institute for Quantitative Biology and Department of Chemistry and Chemical Biology, Rutgers University 174 Frelinghuysen Road, Piscataway, NJ 08854, USA
| | - Magdalena K Gebala
- Department of Biochemistry, Stanford University, Stanford, CA 94305, USA
| | - Maria T Panteva
- BioMaPS Institute for Quantitative Biology and Department of Chemistry and Chemical Biology, Rutgers University 174 Frelinghuysen Road, Piscataway, NJ 08854, USA
| | - Tyler Luchko
- Department of Physics & Astronomy, California State University, Northridge, CA 91330, USA
| | - David A Case
- BioMaPS Institute for Quantitative Biology and Department of Chemistry and Chemical Biology, Rutgers University 174 Frelinghuysen Road, Piscataway, NJ 08854, USA
| | - Darrin M York
- BioMaPS Institute for Quantitative Biology and Department of Chemistry and Chemical Biology, Rutgers University 174 Frelinghuysen Road, Piscataway, NJ 08854, USA
| |
Collapse
|
27
|
Abstract
Ions surround nucleic acids in what is referred to as an ion atmosphere. As a result, the folding and dynamics of RNA and DNA and their complexes with proteins and with each other cannot be understood without a reasonably sophisticated appreciation of these ions' electrostatic interactions. However, the underlying behavior of the ion atmosphere follows physical rules that are distinct from the rules of site binding that biochemists are most familiar and comfortable with. The main goal of this review is to familiarize nucleic acid experimentalists with the physical concepts that underlie nucleic acid-ion interactions. Throughout, we provide practical strategies for interpreting and analyzing nucleic acid experiments that avoid pitfalls from oversimplified or incorrect models. We briefly review the status of theories that predict or simulate nucleic acid-ion interactions and experiments that test these theories. Finally, we describe opportunities for going beyond phenomenological fits to a next-generation, truly predictive understanding of nucleic acid-ion interactions.
Collapse
Affiliation(s)
- Jan Lipfert
- Department of Bionanoscience, Kavli Institute of Nanoscience, Delft University of Technology, 2628 CJ Delft, Netherlands;
| | | | | | | |
Collapse
|
28
|
Xiao S, Zhu H, Wang L, Liang H. DNA conformational flexibility study using phosphate backbone neutralization model. SOFT MATTER 2014; 10:1045-1055. [PMID: 24983118 DOI: 10.1039/c3sm52345d] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/03/2023]
Abstract
Due to the critical role of DNA in the processes of the cell cycle, the structural and physicochemical properties of DNA have long been of concern. In the present work, the effect of interplay between the DNA duplex and metal ions in solution on the DNA structure and conformational flexibility is studied by comparing the structure and dynamic conformational behavior of a duplex in a normal form and its “null isomer” using molecular dynamics methods. It was found that the phosphate neutralization changes the cation atmosphere around the DNA duplex greatly, increases the major groove width, decreases the minor groove width, and reduces the global bending direction preference. We also noted that the probability of BI phosphate linkages increases significantly because of the charge reduction in the backbone phosphate groups. More importantly, we found that the electrostatic effect on the DNA conformational flexibility is dependent on the sequence; that is, the phosphate backbone neutralization induces the global dynamic bending to be less flexible for GC-rich sequences but more flexible for AT-rich sequences.
Collapse
|
29
|
Li PTX. Analysis of diffuse K+ and Mg2+ ion binding to a two-base-pair kissing complex by single-molecule mechanical unfolding. Biochemistry 2013; 52:4991-5001. [PMID: 23842027 DOI: 10.1021/bi400646x] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The folding and stability of RNA tertiary interactions depend critically on cationic conditions. It is usually difficult, however, to isolate such effects on tertiary interactions from those on the entire RNA. By manipulating conformations of single RNA molecules using optical tweezers, we distinguished individual steps of breaking and forming of a two-base-pair kissing interaction from those of secondary folding. The binding of metal ions to the small tertiary structure appeared to be saturable with an apparent Kd of 160 mM for K(+) and 1.5 mM for Mg(2+). The kissing formation was estimated to be associated with binding of ~2-3 diffuse K(+) or Mg(2+) ions. At their saturated binding, Mg(2+) provided ~3 kcal/mol more stabilizing energy to the structure than K(+). Furthermore, the cations change the unkissing forces significantly more than the kissing ones. For example, the presence of Mg(2+) ions increased the average unkissing force from 21 pN to 44 pN, surprisingly high for breaking merely two base pairs; in contrast, the mean kissing force was changed by only 4.5 pN. Interestingly, the differential salt effects on the transition forces were not caused by different changes in the height of the kinetic barriers but were instead attributed to how different molecular structures respond to the applied force. Our results showed the importance of diffuse cation binding to the stability of tertiary interaction and demonstrated the utility of mechanical unfolding in studying tertiary interactions.
Collapse
Affiliation(s)
- Pan T X Li
- Department of Biological Sciences and The RNA Institute, University at Albany, SUNY , Albany, New York 12222, United States
| |
Collapse
|
30
|
El Jarroudi M, Brillard A. Asymptotic analysis of the Poisson-Boltzmann equation in biological membrane channels. Math Biosci 2013; 243:46-56. [PMID: 23429183 DOI: 10.1016/j.mbs.2013.01.011] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2012] [Revised: 01/23/2013] [Accepted: 01/25/2013] [Indexed: 11/29/2022]
Abstract
The Poisson-Boltzmann equation has been increasingly used for the description of biomolecular systems in order to derive their electrostatic properties. We here consider a domain consisting of two living cells which communicate through a system of proteins which assemble at specific membrane areas building microchannels called gap junctions. We describe the asymptotic behavior of the solution of the Poisson-Boltzmann equation posed in this domain. Using Γ-convergence tools, we derive some electrostatic properties of the biological membrane with respect to a vanishing parameter which is simultaneously associated to the membrane thinness, to the diameter of the gap junction microchannels and to the Debye length parameter which characterizes the spatial scale electrostatic interactions between particles within the gap junctions.
Collapse
Affiliation(s)
- Mustapha El Jarroudi
- Université Abdelmalek Essadi, FST Tanger, Département de Mathématiques, BP 416 Tanger, Morocco
| | | |
Collapse
|
31
|
Anthony PC, Sim AY, Chu VB, Doniach S, Block SM, Herschlag D. Electrostatics of nucleic acid folding under conformational constraint. J Am Chem Soc 2012; 134:4607-14. [PMID: 22369617 PMCID: PMC3303965 DOI: 10.1021/ja208466h] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
RNA folding is enabled by interactions between the nucleic acid and its ion atmosphere, the mobile sheath of aqueous ions that surrounds and stabilizes it. Understanding the ion atmosphere requires the interplay of experiment and theory. However, even an apparently simple experiment to probe the ion atmosphere, measuring the dependence of DNA duplex stability upon ion concentration and identity, suffers from substantial complexity, because the unfolded ensemble contains many conformational states that are difficult to treat accurately with theory. To minimize this limitation, we measured the unfolding equilibrium of a DNA hairpin using a single-molecule optical trapping assay, in which the unfolded state is constrained to a limited set of elongated conformations. The unfolding free energy increased linearly with the logarithm of monovalent cation concentration for several cations, such that smaller cations tended to favor the folded state. Mg(2+) stabilized the hairpin much more effectively at low concentrations than did any of the monovalent cations. Poisson-Boltzmann theory captured trends in hairpin stability measured for the monovalent cation titrations with reasonable accuracy, but failed to do so for the Mg(2+) titrations. This finding is consistent with previous work, suggesting that Poisson-Boltzmann and other mean-field theories fail for higher valency cations where ion-ion correlation effects may become significant. The high-resolution data herein, because of the straightforward nature of both the folded and the unfolded states, should serve as benchmarks for the development of more accurate electrostatic theories that will be needed for a more quantitative and predictive understanding of nucleic acid folding.
Collapse
Affiliation(s)
| | - Adelene Y.L. Sim
- Department of Applied Physics, Stanford University, Stanford, CA 94305
| | - Vincent B. Chu
- Department of Applied Physics, Stanford University, Stanford, CA 94305
| | - Sebastian Doniach
- Department of Applied Physics, Stanford University, Stanford, CA 94305
- Department of Physics, Stanford University, Stanford, CA 94305
| | - Steven M. Block
- Department of Applied Physics, Stanford University, Stanford, CA 94305
- Department of Biology, Stanford University, Stanford, CA 94305
| | - Daniel Herschlag
- Department of Biochemistry, Stanford University, Stanford, CA 94305
| |
Collapse
|
32
|
Kirmizialtin S, Pabit S, Meisburger S, Pollack L, Elber R. RNA and its ionic cloud: solution scattering experiments and atomically detailed simulations. Biophys J 2012; 102:819-28. [PMID: 22385853 PMCID: PMC3283807 DOI: 10.1016/j.bpj.2012.01.013] [Citation(s) in RCA: 84] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2011] [Revised: 12/22/2011] [Accepted: 01/12/2012] [Indexed: 01/12/2023] Open
Abstract
RNA molecules play critical roles in many cellular processes. Traditionally viewed as genetic messengers, RNA molecules were recently discovered to have diverse functions related to gene regulation and expression. RNA also has great potential as a therapeutic and a tool for further investigation of gene regulation. Metal ions are an integral part of RNA structure and should be considered in any experimental or theoretical study of RNA. Here, we report a multidisciplinary approach that combines anomalous small-angle x-ray scattering and molecular-dynamics (MD) simulations with explicit solvent and ions around RNA. From experiment and simulation results, we find excellent agreement in the number and distribution of excess monovalent and divalent ions around a short RNA duplex. Although similar agreement can be obtained from a continuum description of the solvent and mobile ions (by solving the Poisson-Boltzmann equation and accounting for finite ion size), the use of MD is easily extended to flexible RNA systems with thermal fluctuations. Therefore, we also model a short RNA pseudoknot and find good agreement between the MD results and the experimentally derived solution structures. Surprisingly, both deviate from crystal structure predictions. These favorable comparisons of experiment and simulations encourage work on RNA in all-atom dynamic models.
Collapse
Affiliation(s)
- Serdal Kirmizialtin
- Department of Chemistry and Biochemistry, The University of Texas at Austin, Austin, Texas
- Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas
| | - Suzette A. Pabit
- School of Applied and Engineering Physics, Cornell University, Ithaca, New York
| | - Steve P. Meisburger
- School of Applied and Engineering Physics, Cornell University, Ithaca, New York
| | - Lois Pollack
- School of Applied and Engineering Physics, Cornell University, Ithaca, New York
| | - Ron Elber
- Department of Chemistry and Biochemistry, The University of Texas at Austin, Austin, Texas
- Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas
| |
Collapse
|
33
|
Kirmizialtin S, Silalahi A, Elber R, Fenley M. The ionic atmosphere around A-RNA: Poisson-Boltzmann and molecular dynamics simulations. Biophys J 2012; 102:829-38. [PMID: 22385854 PMCID: PMC3283773 DOI: 10.1016/j.bpj.2011.12.055] [Citation(s) in RCA: 49] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2011] [Revised: 12/22/2011] [Accepted: 12/28/2011] [Indexed: 10/28/2022] Open
Abstract
The distributions of different cations around A-RNA are computed by Poisson-Boltzmann (PB) equation and replica exchange molecular dynamics (MD). Both the nonlinear PB and size-modified PB theories are considered. The number of ions bound to A-RNA, which can be measured experimentally, is well reproduced in all methods. On the other hand, the radial ion distribution profiles show differences between MD and PB. We showed that PB results are sensitive to ion size and functional form of the solvent dielectric region but not the solvent dielectric boundary definition. Size-modified PB agrees with replica exchange molecular dynamics much better than nonlinear PB when the ion sizes are chosen from atomistic simulations. The distribution of ions 14 Å away from the RNA central axis are reasonably well reproduced by size-modified PB for all ion types with a uniform solvent dielectric model and a sharp dielectric boundary between solvent and RNA. However, this model does not agree with MD for shorter distances from the A-RNA. A distance-dependent solvent dielectric function proposed by another research group improves the agreement for sodium and strontium ions, even for shorter distances from the A-RNA. However, Mg(2+) distributions are still at significant variances for shorter distances.
Collapse
Affiliation(s)
- Serdal Kirmizialtin
- Department of Chemistry and Biochemistry, Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas
| | - Alexander R.J. Silalahi
- Department of Physics and the Institute of Molecular Biophysics, Florida State University, Tallahassee, Florida
| | - Ron Elber
- Department of Chemistry and Biochemistry, Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas
| | - Marcia O. Fenley
- Department of Physics and the Institute of Molecular Biophysics, Florida State University, Tallahassee, Florida
| |
Collapse
|
34
|
Li W, Nordenskiöld L, Mu Y. Sequence-Specific Mg2+–DNA Interactions: A Molecular Dynamics Simulation Study. J Phys Chem B 2011; 115:14713-20. [DOI: 10.1021/jp2052568] [Citation(s) in RCA: 50] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Weifeng Li
- School of Physical and Mathematical Sciences, 21 Nanyang Link, and ‡School of Biological Sciences, 60 Nanyang Drive, Nanyang Technological University, Singapore
| | - Lars Nordenskiöld
- School of Physical and Mathematical Sciences, 21 Nanyang Link, and ‡School of Biological Sciences, 60 Nanyang Drive, Nanyang Technological University, Singapore
| | - Yuguang Mu
- School of Physical and Mathematical Sciences, 21 Nanyang Link, and ‡School of Biological Sciences, 60 Nanyang Drive, Nanyang Technological University, Singapore
| |
Collapse
|
35
|
Ryde U, Mata RA, Grimme S. Does DFT-D estimate accurate energies for the binding of ligands to metal complexes? Dalton Trans 2011; 40:11176-83. [DOI: 10.1039/c1dt10867k] [Citation(s) in RCA: 73] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
36
|
Pabit SA, Meisburger SP, Li L, Blose JM, Jones CD, Pollack L. Counting ions around DNA with anomalous small-angle X-ray scattering. J Am Chem Soc 2010; 132:16334-6. [PMID: 21047071 DOI: 10.1021/ja107259y] [Citation(s) in RCA: 73] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The majority of charge-compensating ions around nucleic acids form a diffuse counterion "cloud" that is not amenable to investigation by traditional methods that rely on rigid structural interactions. Although various techniques have been employed to characterize the ion atmosphere around nucleic acids, only anomalous small-angle X-ray scattering (ASAXS) provides information about the spatial distribution of ions. Here we present an experimentally straightforward extension of ASAXS that can be used to count the number of ions around nucleic acids.
Collapse
Affiliation(s)
- Suzette A Pabit
- School of Applied and Engineering Physics, Cornell University, Ithaca, New York 14853, United States
| | | | | | | | | | | |
Collapse
|
37
|
Lee DJ, Wynveen A, Kornyshev AA, Leikin S. Undulations enhance the effect of helical structure on DNA interactions. J Phys Chem B 2010; 114:11668-80. [PMID: 20718454 PMCID: PMC2937169 DOI: 10.1021/jp104552u] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
During the past decade, theory and experiments have provided clear evidence that specific helical patterns of charged groups and adsorbed (condensed) counterions on the DNA surface are responsible for many important features of DNA-DNA interactions in hydrated aggregates. The effects of helical structure on DNA-DNA interactions result from a preferential juxtaposition of the negatively charged sugar phosphate backbone with counterions bound within the grooves of the opposing molecule. Analysis of X-ray diffraction experiments confirmed the mutual alignment of parallel molecules in hydrated aggregates required for such juxtaposition. However, it remained unclear how this alignment and molecular interactions might be affected by intrinsic and thermal fluctuations, which cause structural deviations away from an ideal double helical conformation. We previously argued that the torsional flexibility of DNA allows the molecules to adapt their structure to accommodate a more electrostatically favorable alignment between molecules, partially compensating disruptive fluctuation effects. In the present work, we develop a more comprehensive theory, incorporating also stretching and bending fluctuations of DNA. We found the effects of stretching to be qualitatively and quantitatively similar to those of twisting fluctuations. However, this theory predicts more dramatic and surprising effects of bending. Undulations of DNA in hydrated aggregates strongly amplify rather than weaken the helical structure effects. They enhance the structural adaptation, leading to better alignment of neighboring molecules and pushing the geometry of the DNA backbone closer to that of an ideal helix. These predictions are supported by a quantitative comparison of the calculated and measured osmotic pressures in DNA.
Collapse
Affiliation(s)
- D. J. Lee
- To whom correspondence should be addressed. (D.J.L.) . (A.W.) . (A.A.K.) . (S.L.) Tel: 1-301-594-8314; FAX: 1-301-402-0292;
| | - A. Wynveen
- To whom correspondence should be addressed. (D.J.L.) . (A.W.) . (A.A.K.) . (S.L.) Tel: 1-301-594-8314; FAX: 1-301-402-0292;
| | - A. A Kornyshev
- To whom correspondence should be addressed. (D.J.L.) . (A.W.) . (A.A.K.) . (S.L.) Tel: 1-301-594-8314; FAX: 1-301-402-0292;
| | - S. Leikin
- To whom correspondence should be addressed. (D.J.L.) . (A.W.) . (A.A.K.) . (S.L.) Tel: 1-301-594-8314; FAX: 1-301-402-0292;
| |
Collapse
|
38
|
Christian EL, Anderson VE, Carey PR, Harris ME. A quantitative Raman spectroscopic signal for metal-phosphodiester interactions in solution. Biochemistry 2010; 49:2869-79. [PMID: 20180599 DOI: 10.1021/bi901866u] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Accurate identification and quantification of metal ion-phosphodiester interactions are essential for understanding the role of metal ions as determinants of three-dimensional folding of large RNAs and as cofactors in the active sites of both RNA and protein phosphodiesterases. Accomplishing this goal is difficult due to the dynamic and complex mixture of direct and indirect interactions formed with nucleic acids and other phosphodiesters in solution. To address this issue, Raman spectroscopy has been used to measure changes in bond vibrational energies due to metal interactions. However, the contributions of inner-sphere, H-bonding, and electrostatic interactions to the Raman spectrum of phosphoryl oxygens have not been analyzed quantitatively. Here, we report that all three forms of metal ion interaction result in attenuation of the Raman signal for the symmetric vibration of the nonbridging phosphate oxygens (nu(s)PO(2)(-)), while only inner-sphere coordination gives rise to an apparent shift of nu(s)PO(2)(-) to higher wavenumbers (nu(s)PO(2)(-)M) in solution. Formation of nu(s)PO(2)(-)M is shown to be both dependent on metal ion identity and an accurate measure of site-specific metal ion binding. In addition, the spectroscopic parameter reflecting the energetic difference between nu(s)PO(2)(-) and nu(s)PO(2)(-)M (DeltanuM) is largely insensitive to changes in phosphodiester structure but strongly dependent on the absolute electronegativity and hardness of the interacting metal ion. Together, these studies provide strong experimental support for the use of nu(s)PO(2)(-)M and DeltanuM as general spectroscopic features for the quantitative analysis of metal binding affinity and the identification of metal ions associated with phosphodiesters in solution.
Collapse
Affiliation(s)
- Eric L Christian
- Center for RNA Molecular Biology, Case Western Reserve University School of Medicine,Cleveland, Ohio 44106, USA.
| | | | | | | |
Collapse
|
39
|
Effect of temperature on volumetric and isentropic compressibility of glycine in aqueous KCl/NaCl solutions. J Mol Liq 2010. [DOI: 10.1016/j.molliq.2010.01.003] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
40
|
Temiz AN, Benos PV, Camacho CJ. Electrostatic hot spot on DNA-binding domains mediates phosphate desolvation and the pre-organization of specificity determinant side chains. Nucleic Acids Res 2010; 38:2134-44. [PMID: 20047959 PMCID: PMC2853105 DOI: 10.1093/nar/gkp1132] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
A major obstacle towards elucidating the molecular basis of transcriptional regulation is the lack of a detailed understanding of the interplay between non-specific and specific protein–DNA interactions. Based on molecular dynamics simulations of C2H2 zinc fingers (ZFs) and engrailed homeodomain transcription factors (TFs), we show that each of the studied DNA-binding domains has a set of highly constrained side chains in preset configurations ready to form hydrogen bonds with the DNA backbone. Interestingly, those domains that bury their recognition helix into the major groove are found to have an electrostatic hot spot for Cl− ions located on the same binding cavity as the most buried DNA phosphate. The spot is characterized by three protein hydrogen bond donors, often including two basic side chains. If bound, Cl− ions, likely mimicking phosphates, steer side chains that end up forming specific contacts with bases into bound-like conformations. These findings are consistent with a multi-step DNA-binding mechanism in which a pre-organized set of TF side chains assist in the desolvation of phosphates into well defined sites, prompting the re-organization of specificity determining side chains into conformations suitable for the recognition of their cognate sequence.
Collapse
Affiliation(s)
- Alpay N Temiz
- Department of Computational Biology, School of Medicine, University of Pittsburgh, Pittsburgh, PA 15260, USA
| | | | | |
Collapse
|
41
|
Kongsted J, Söderhjelm P, Ryde U. How accurate are continuum solvation models for drug-like molecules? J Comput Aided Mol Des 2009; 23:395-409. [DOI: 10.1007/s10822-009-9271-6] [Citation(s) in RCA: 63] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2009] [Accepted: 04/15/2009] [Indexed: 12/01/2022]
|
42
|
Pabit SA, Qiu X, Lamb JS, Li L, Meisburger SP, Pollack L. Both helix topology and counterion distribution contribute to the more effective charge screening in dsRNA compared with dsDNA. Nucleic Acids Res 2009; 37:3887-96. [PMID: 19395592 PMCID: PMC2709557 DOI: 10.1093/nar/gkp257] [Citation(s) in RCA: 67] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
The recent discovery of the RNA interference mechanism emphasizes the biological importance of short, isolated, double-stranded (ds) RNA helices and calls for a complete understanding of the biophysical properties of dsRNA. However, most previous studies of the electrostatics of nucleic acid duplexes have focused on DNA. Here, we present a comparative investigation of electrostatic effects in RNA and DNA. Using resonant (anomalous) and non-resonant small-angle X-ray scattering, we characterized the charge screening efficiency and counterion distribution around short (25 bp) dsDNA and RNA molecules of comparable sequence. Consistent with theoretical predictions, we find counterion mediated screening to be more efficient for dsRNA than dsDNA. Furthermore, the topology of the RNA A-form helix alters the spatial distribution of counterions relative to B-form DNA. The experimental results reported here agree well with ion-size-corrected non-linear Poisson-Boltzmann calculations. We propose that differences in electrostatic properties aid in selective recognition of different types of short nucleic acid helices by target binding partners.
Collapse
Affiliation(s)
- Suzette A Pabit
- School of Applied and Engineering Physics, Cornell University, Ithaca, NY 14853, USA
| | | | | | | | | | | |
Collapse
|
43
|
Greenfeld M, Herschlag D. Probing nucleic acid-ion interactions with buffer exchange-atomic emission spectroscopy. Methods Enzymol 2009; 469:375-89. [PMID: 20946799 DOI: 10.1016/s0076-6879(09)69018-2] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
The ion atmosphere of nucleic acids directly affects measured biochemical and biophysical properties. However, study of the ion atmosphere is difficult due to its diffuse and dynamic nature. Standard techniques available have significant limitations in sensitivity, specificity, and directness of the assays. Buffer exchange-atomic emission spectroscopy (BE-AES) was developed to overcome many of the limitations of previously available techniques. This technique can provide a complete accounting of all ions constituting the ionic atmosphere of a nucleic acid at thermodynamic equilibrium. Although initially developed for the study of the ion atmosphere of nucleic acids, BE-AES has also been applied to study site-bound ions in RNA and protein.
Collapse
Affiliation(s)
- Max Greenfeld
- Department of Chemical Engineering and Biochemistry, Stanford University, Stanford, California, USA
| | | |
Collapse
|
44
|
Abstract
Understanding how RNA folds and what causes it to unfold has become more important as knowledge of the diverse functions of RNA has increased. Here we review the contributions of single-molecule experiments to providing answers to questions such as: How much energy is required to unfold a secondary or tertiary structure? How fast is the process? How do helicases unwind double helices? Are the unwinding activities of RNA-dependent RNA polymerases and of ribosomes different from other helicases? We discuss the use of optical tweezers to monitor the unfolding activities of helicases, polymerases, and ribosomes, and to apply force to unfold RNAs directly. We also review the applications of fluorescence and fluorescence resonance energy transfer to measure RNA dynamics.
Collapse
Affiliation(s)
- Pan T X Li
- Department of Biological Sciences, University at Albany, State University of New York, Albany, NY 12222, USA.
| | | | | |
Collapse
|
45
|
Bai Y, Chu VB, Lipfert J, Pande VS, Herschlag D, Doniach S. Critical assessment of nucleic acid electrostatics via experimental and computational investigation of an unfolded state ensemble. J Am Chem Soc 2008; 130:12334-41. [PMID: 18722445 DOI: 10.1021/ja800854u] [Citation(s) in RCA: 62] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Electrostatic forces, acting between helices and modulated by the presence of the ion atmosphere, are key determinants in the energetic balance that governs RNA folding. Previous studies have employed Poisson-Boltzmann (PB) theory to compute the energetic contribution of these forces in RNA folding. However, the complex interaction of these electrostatic forces with RNA features such as tertiary contact formation, specific ion-binding, and complex interhelical junctions present in prior studies precluded a rigorous evaluation of PB theory, especially in physiologically important Mg(2+) solutions. To critically assess PB theory, we developed a model system that isolates these electrostatic forces. The model system, composed of two DNA duplexes tethered by a polyethylene glycol junction, is an analog for the unfolded state of canonical helix-junction-helix motifs found in virtually all structured RNAs. This model system lacks the complicating features that have precluded a critical assessment of PB in prior studies, ensuring that interhelical electrostatic forces dominate the behavior of the system. The system's simplicity allows PB predictions to be directly compared with small-angle X-ray scattering experiments over a range of monovalent and divalent ion concentrations. These comparisons indicate that PB is a reasonable description of the underlying electrostatic energies for monovalent ions, but large deviations are observed for divalent ions. The validation of PB for monovalent solutions allows analysis of the change in the conformational ensemble of this simple motif as salt concentration is changed. Addition of ions allows the motif to sample more compact microstates, increasing its conformational entropy. The increase of conformational entropy presents an additional barrier to folding by stabilizing the unfolded state. Neglecting this effect will adversely impact the accuracy of folding analyses and models.
Collapse
Affiliation(s)
- Yu Bai
- Department of Biochemistry, Stanford University, California 94305, USA
| | | | | | | | | | | |
Collapse
|
46
|
Ozkirimli E, Yadav SS, Miller WT, Post CB. An electrostatic network and long-range regulation of Src kinases. Protein Sci 2008; 17:1871-80. [PMID: 18687871 DOI: 10.1110/ps.037457.108] [Citation(s) in RCA: 47] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Abstract
The regulatory mechanism of Src tyrosine kinases includes conformational activation by a change in the catalytic domain tertiary structure and in domain-domain contacts between the catalytic domain and the SH2/SH3 regulatory domains. The kinase is activated when tyrosine phosphorylation occurs on the activation loop, but without phosphorylation of the C-terminal tail. Activation also occurs by allostery when contacts between the catalytic domain (CD) and the regulatory SH3 and SH2 domains are released as a result of exogenous protein binding. The aim of this work is to examine the proposed role of an electrostatic network in the conformational transition and to elucidate the molecular mechanism for long-range, allosteric conformational activation by using a combination of experimental enzyme kinetics and nonequilibrium molecular dynamics simulations. Salt dependence of the induction phase is observed in kinetic assays and supports the role of an electrostatic network in the transition. In addition, simulations provide evidence that allosteric activation involves a concerted motion coupling highly conserved residues, and spanning several nanometers from the catalytic site to the regulatory domain interface to communicate between the CD and the regulatory domains.
Collapse
Affiliation(s)
- Elif Ozkirimli
- 1Medicinal Chemistry and Molecular Pharmacology Department, Markey Center for Structural Biology and Purdue Cancer Center, Purdue University, West Lafayette, Indiana 47907-2091, USA
| | | | | | | |
Collapse
|
47
|
Wittayanarakul K, Hannongbua S, Feig M. Accurate prediction of protonation state as a prerequisite for reliable MM-PB(GB)SA binding free energy calculations of HIV-1 protease inhibitors. J Comput Chem 2008; 29:673-85. [PMID: 17849388 DOI: 10.1002/jcc.20821] [Citation(s) in RCA: 57] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
Binding free energies were calculated for the inhibitors lopinavir, ritonavir, saquinavir, indinavir, amprenavir, and nelfinavir bound to HIV-1 protease. An MMPB/SA-type analysis was applied to conformational samples from 3 ns explicit solvent molecular dynamics simulations of the enzyme-inhibitor complexes. Binding affinities and the sampled conformations of the inhibitor and enzyme were compared between different HIV-1 protease protonation states to find the most likely protonation state of the enzyme in the complex with each of the inhibitors. The resulting set of protonation states leads to good agreement between calculated and experimental binding affinities. Results from the MMPB/SA analysis are compared with an explicit/implicit hybrid scheme and with MMGB/SA methods. It is found that the inclusion of explicit water molecules may offer a slight advantage in reproducing absolute binding free energies while the use of the Generalized Born approximation significantly affects the accuracy of the calculated binding affinities.
Collapse
|
48
|
Bai Y, Greenfeld M, Travers K, Chu VB, Lipfert J, Doniach S, Herschlag D. Quantitative and comprehensive decomposition of the ion atmosphere around nucleic acids. J Am Chem Soc 2007; 129:14981-8. [PMID: 17990882 PMCID: PMC3167487 DOI: 10.1021/ja075020g] [Citation(s) in RCA: 215] [Impact Index Per Article: 12.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The ion atmosphere around nucleic acids critically affects biological and physical processes such as chromosome packing, RNA folding, and molecular recognition. However, the dynamic nature of the ion atmosphere renders it difficult to characterize. The basic thermodynamic description of this atmosphere, a full accounting of the type and number of associated ions, has remained elusive. Here we provide the first complete accounting of the ion atmosphere, using buffer equilibration and atomic emission spectroscopy (BE-AES) to accurately quantitate the cation association and anion depletion. We have examined the influence of ion size and charge on ion occupancy around simple, well-defined DNA molecules. The relative affinity of monovalent and divalent cations correlates inversely with their size. Divalent cations associate preferentially over monovalent cations; e.g., with Na+ in 4-fold excess of Mg2+ (20 vs 5 mM), the ion atmosphere nevertheless has 3-fold more Mg2+ than Na+. Further, the dicationic polyamine putrescine2+ does not compete effectively for association relative to divalent metal ions, presumably because of its lower charge density. These and other BE-AES results can be used to evaluate and guide the improvement of electrostatic treatments. As a first step, we compare the BE-AES results to predictions from the widely used nonlinear Poisson Boltzmann (NLPB) theory and assess the applicability and precision of this theory. In the future, BE-AES in conjunction with improved theoretical models, can be applied to complex binding and folding equilibria of nucleic acids and their complexes, to parse the electrostatic contribution from the overall thermodynamics of important biological processes.
Collapse
Affiliation(s)
| | | | | | | | | | | | - Daniel Herschlag
- Corresponding Author: Department of Biochemistry, Beckman Center B471, Stanford University, CA 94305, USA. Tel. 650 723-9442, FAX. 650 723-6783;
| |
Collapse
|
49
|
Bugs MR, Cornélio ML. Analysis of the Ethidium Bromide Bound to DNA by Photoacoustic and FTIR Spectroscopy¶. Photochem Photobiol 2007. [DOI: 10.1562/0031-8655(2001)0740512aotebb2.0.co2] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
50
|
Abascal JLF, Domercq M, Montoro JCG. Computer Simulation of the Ionic Atmosphere around Z-DNA. J Phys Chem B 2006; 110:25080-90. [PMID: 17149933 DOI: 10.1021/jp064199z] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
We describe a coarse-grained model for Z-DNA that mimics the DNA shape with a relatively small number of repulsive interaction sites. In addition, negative charges are placed at the phosphate positions. The ionic atmosphere around this grooved Z-DNA model is then investigated with Monte Carlo simulation. Cylindrically averaged concentration profiles as well as the spatial distribution of ions have been calculated. The results are compared to those for other DNA models differing in the repulsive core. This allows the examination of the effect of the DNA shape in the ionic distribution. It is seen that the penetrability of the ions to the DNA groove plays an important role in the ionic distribution. The results are also compared with those reported for B-DNA. In both conformers the ions are structured in alternating layers of positive and negative charge. In Z-DNA the layers are more or less concentric to the molecular axis. Besides, no coions enter into the single groove of this conformer. On the contrary, the alternating layers of B-DNA are also structured along the axial coordinate with some coions penetrating into the major groove. In both cases we have found five preferred locations of the counterions and two for the coions. The concentration of counterions reaches its absolute maximum at the narrow Z-DNA groove and at the minor groove of B-DNA, the value of the maximum being higher in the Z conformer.
Collapse
Affiliation(s)
- J L F Abascal
- Departamento de Química-Física, Facultad de Ciencias Químicas, Universidad Complutense, 28040 Madrid, Spain
| | | | | |
Collapse
|