1
|
Deshmukh SH, Nachaki EO, Kuroda DG. Uncovering the binding nature of thiocyanate in contact ion pairs with lithium ions. J Chem Phys 2024; 161:034507. [PMID: 39017430 DOI: 10.1063/5.0216491] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2024] [Accepted: 06/28/2024] [Indexed: 07/18/2024] Open
Abstract
Ion pair formation is a fundamental molecular process that occurs in a wide variety of systems, including electrolytes, biological systems, and materials. In solution, the thiocyanate (SCN-) anion interacts with cations to form contact ion pairs (CIPs). Due to its ambidentate nature, thiocyanate can bind through either its sulfur or nitrogen atoms, depending on the solvent. This study focuses on the binding nature of thiocyanate with lithium ions as a function of the solvents using FTIR, 2D infrared spectroscopy (2DIR) spectroscopies, and theoretical calculations. The study reveals that the SCN- binding mode (S or N end) in CIPs can be identified through 2DIR spectroscopy but not by linear IR spectroscopy. Linear IR spectroscopy shows that the CN stretch frequencies are too close to one another to separate N- and S-bound CIPs. Moreover, the IR spectrum shows that the S-C stretch presents different frequencies for the salt in different solvents, but it is related to the anion speciation rather than to its binding mode. A similar trend is observed for the anion bend. 2DIR spectra show different dynamics for N-bound and S-bound thiocyanate. In particular, the frequency-frequency correlation function (FFCF) dynamics extracted from the 2DIR spectra have a single picosecond exponential decay for N-bound thiocyanate and a biexponential decay for S-bound thiocyanate, consistent with the binding mode of the anion. Finally, it is also observed that the binding mode also affects the line shape parameters, probably due to the different molecular mechanisms of the FFCF for N- and S-bound CIPs.
Collapse
Affiliation(s)
- Samadhan H Deshmukh
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - Ernest O Nachaki
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - Daniel G Kuroda
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| |
Collapse
|
2
|
Yang J, Cong Y, Li Y, Li H. Machine Learning Approach Based on a Range-Corrected Deep Potential Model for Efficient Vibrational Frequency Computation. J Chem Theory Comput 2023; 19:6366-6374. [PMID: 37652890 DOI: 10.1021/acs.jctc.3c00386] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/02/2023]
Abstract
As an ensemble average result, vibrational spectrum simulation can be time-consuming with high accuracy methods. We present a machine learning approach based on the range-corrected deep potential (DPRc) model to improve the computing efficiency. The DPRc method divides the system into "probe region" and "solvent region"; "solvent-solvent" interactions are not counted in the neural network. We applied the approach to two systems: formic acid C═O stretching and MeCN C≡N stretching vibrational frequency shifts in water. All data sets were prepared using the quantum vibration perturbation approach. Effects of different region divisions, one-body correction, cut range, and training data size were tested. The model with a single-molecule "probe region" showed stable accuracy; it ran roughly 10 times faster than regular deep potential and reduced the training time by about four. The approach is efficient, easy to apply, and extendable to calculating various spectra.
Collapse
Affiliation(s)
- Jitai Yang
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, 2519 Jiefang Road, Changchun 130023, P. R. China
| | - Yang Cong
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, 2519 Jiefang Road, Changchun 130023, P. R. China
| | - You Li
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, 2519 Jiefang Road, Changchun 130023, P. R. China
| | - Hui Li
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, 2519 Jiefang Road, Changchun 130023, P. R. China
| |
Collapse
|
3
|
Rapid vibrational relaxation dynamics and strong Solute-Solvent interactions of the nitrite anion. Chem Phys Lett 2023. [DOI: 10.1016/j.cplett.2022.140265] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
|
4
|
Kundu A, Mamatkulov SI, Brünig FN, Bonthuis DJ, Netz RR, Elsaesser T, Fingerhut BP. Short-Range Cooperative Slow-down of Water Solvation Dynamics Around SO 4 2--Mg 2+ Ion Pairs. ACS PHYSICAL CHEMISTRY AU 2022; 2:506-514. [PMID: 36465835 PMCID: PMC9706802 DOI: 10.1021/acsphyschemau.2c00034] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/01/2022] [Revised: 09/15/2022] [Accepted: 09/15/2022] [Indexed: 06/17/2023]
Abstract
The presence of ions affects the structure and dynamics of water on a multitude of length and time scales. In this context, pairs of Mg2+ and SO4 2- ions in water constitute a prototypical system for which conflicting pictures of hydration geometries and dynamics have been reported. Key issues are the molecular pair and solvation shell geometries, the spatial range of electric interactions, and their impact on solvation dynamics. Here, we introduce asymmetric SO4 2- stretching vibrations as new and most specific local probes of solvation dynamics that allow to access ion hydration dynamics at the dilute concentration (0.2 M) of a native electrolyte environment. Highly sensitive heterodyne 2D-IR spectroscopy in the fingerprint region of the SO4 2- ions around 1100 cm-1 reveals a specific slow-down of solvation dynamics for hydrated MgSO4 and for Na2SO4 in the presence of Mg2+ ions, which manifests as a retardation of spectral diffusion compared to aqueous Na2SO4 solutions in the absence of Mg2+ ions. Extensive molecular dynamics and density functional theory QM/MM simulations provide a microscopic view of the observed ultrafast dephasing and hydration dynamics. They suggest a molecular picture where the slow-down of hydration dynamics arises from the structural peculiarities of solvent-shared SO4 2--Mg2+ ion pairs.
Collapse
Affiliation(s)
- Achintya Kundu
- Max-Born-Institut
Für Nichtlineare Optik und Kurzzeitspektroskopie, Berlin12489, Germany
| | - Shavkat I. Mamatkulov
- Institute
of Material Sciences of Uzbekistan Academy of Sciences, Tashkent100084, Uzbekistan
| | | | - Douwe Jan Bonthuis
- Institute
of Theoretical and Computational Physics, Graz University of Technology, Graz8010, Austria
| | - Roland R. Netz
- Fachbereich
Physik, Freie Universität Berlin, Berlin14195, Germany
| | - Thomas Elsaesser
- Max-Born-Institut
Für Nichtlineare Optik und Kurzzeitspektroskopie, Berlin12489, Germany
| | - Benjamin P. Fingerhut
- Max-Born-Institut
Für Nichtlineare Optik und Kurzzeitspektroskopie, Berlin12489, Germany
| |
Collapse
|
5
|
Nir O, Oren Y, Musie, Atshba W, Chandra A, Geller Y, Chaudhary M, Monat L, Singh P, Zevenhoven R. Reactive transport in membrane separation modeling: a perspective. Chem Eng Res Des 2022. [DOI: 10.1016/j.cherd.2022.09.054] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
6
|
Ion-pairing equilibria and kinetics of dimethyl phosphate: A model for counter-ion binding to the phosphate backbone of nucleic acids. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.119868] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
7
|
Puyo-Fourtine J, Juillé M, Hénin J, Clavaguéra C, Duboué-Dijon E. Consistent Picture of Phosphate-Divalent Cation Binding from Models with Implicit and Explicit Electronic Polarization. J Phys Chem B 2022; 126:4022-4034. [PMID: 35608554 DOI: 10.1021/acs.jpcb.2c01158] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The binding of divalent cations to the ubiquitous phosphate group is essential for a number of key biological processes, such as DNA compaction, RNA folding, or interactions of some proteins with membranes. Yet, probing their binding sites, modes, and associated binding free energy is a challenge for both experiments and simulations. In simulations, standard force fields strongly overestimate the interaction between phosphate groups and divalent cations. Here, we examine how different strategies to include electronic polarization effects in force fields─implicitly, through the use of scaled charges or pair-specific Lennard-Jones parameters, or explicitly, with the polarizable force fields Drude and AMOEBA─capture the interactions of a model phosphate compound, dimethyl phosphate, with calcium and magnesium divalent cations. We show that both implicit and explicit approaches, when carefully parameterized, are successful in capturing the overall binding free energy and that common trends emerge from the comparison of different simulation approaches. Overall, the binding is very moderate, slightly weaker for Ca2+ than Mg2+, and the solvent-shared ion pair is slightly more stable than the contact monodentate ion pair. The bidentate ion pair is higher in energy (or even fully unstable for Mg2+). Our results thus suggest practical ways to capture the divalent cations with biomolecular phosphate groups in complex biochemical systems. In particular, the computational efficiency of implicit models makes them ideally suited for large-scale simulations of biological assemblies, with improved accuracy compared to state-of-the-art fixed-charge force fields.
Collapse
Affiliation(s)
- Julie Puyo-Fourtine
- CNRS, Université Paris Cité, UPR9080, Laboratoire de Biochimie Théorique, 13 Rue Pierre et Marie Curie, 75005 Paris, France.,Institut de Biologie Physico-Chimique - Fondation Edmond de Rothschild, PSL Research University, 75005 Paris, France
| | - Marie Juillé
- CNRS, Université Paris Cité, UPR9080, Laboratoire de Biochimie Théorique, 13 Rue Pierre et Marie Curie, 75005 Paris, France.,Institut de Biologie Physico-Chimique - Fondation Edmond de Rothschild, PSL Research University, 75005 Paris, France
| | - Jérôme Hénin
- CNRS, Université Paris Cité, UPR9080, Laboratoire de Biochimie Théorique, 13 Rue Pierre et Marie Curie, 75005 Paris, France.,Institut de Biologie Physico-Chimique - Fondation Edmond de Rothschild, PSL Research University, 75005 Paris, France
| | - Carine Clavaguéra
- Université Paris-Saclay, CNRS, Institut de Chimie Physique, UMR8000, 91405 Orsay, France
| | - Elise Duboué-Dijon
- CNRS, Université Paris Cité, UPR9080, Laboratoire de Biochimie Théorique, 13 Rue Pierre et Marie Curie, 75005 Paris, France.,Institut de Biologie Physico-Chimique - Fondation Edmond de Rothschild, PSL Research University, 75005 Paris, France
| |
Collapse
|
8
|
Price DA, Wedamulla P, Hill TD, Loth TM, Moran SD. The polarization dependence of 2D IR cross-peaks distinguishes parallel-stranded and antiparallel-stranded DNA G-quadruplexes. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2022; 267:120596. [PMID: 34801392 DOI: 10.1016/j.saa.2021.120596] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/20/2021] [Revised: 11/01/2021] [Accepted: 11/04/2021] [Indexed: 06/13/2023]
Abstract
Guanine-rich nucleic acid sequences have a tendency to form four-stranded non-canonical motifs known as G-quadruplexes. These motifs may adopt a wide range of structures characterized by size, strand orientation, guanine base conformation, and fold topology. Using three K+-bound model systems, we show that vibrational coupling between guanine C6 = O and ring modes varies between parallel-stranded and antiparallel-stranded G-quadruplexes, and that such structures can be distinguished by comparison of the polarization dependences of cross-peaks in their two-dimensional infrared (2D IR) spectra. Combined with previously defined vibrational frequency trends, this analysis reveals key features of a 30-nucleotide unimolecular variant of the Bcl-2 proximal promoter that are consistent with its reported structure. This study shows that 2D IR spectroscopy is a convenient method for analyzing G-quadruplex structures that can be applied to complex sequences where traditional high-resolution methods are limited by solubility and disorder.
Collapse
Affiliation(s)
- David A Price
- School of Chemical and Biomolecular Sciences, Southern Illinois University Carbondale, 1245 Lincoln Drive MC 4409, Carbondale, IL 62901, United States
| | - Poornima Wedamulla
- School of Chemical and Biomolecular Sciences, Southern Illinois University Carbondale, 1245 Lincoln Drive MC 4409, Carbondale, IL 62901, United States
| | - Tayler D Hill
- School of Chemical and Biomolecular Sciences, Southern Illinois University Carbondale, 1245 Lincoln Drive MC 4409, Carbondale, IL 62901, United States
| | - Taylor M Loth
- School of Chemical and Biomolecular Sciences, Southern Illinois University Carbondale, 1245 Lincoln Drive MC 4409, Carbondale, IL 62901, United States
| | - Sean D Moran
- School of Chemical and Biomolecular Sciences, Southern Illinois University Carbondale, 1245 Lincoln Drive MC 4409, Carbondale, IL 62901, United States.
| |
Collapse
|
9
|
Cong Y, Zhai Y, Yang J, Grofe A, Gao J, Li H. Quantum vibration perturbation approach with polyatomic probe in simulating infrared spectra. Phys Chem Chem Phys 2021; 24:1174-1182. [PMID: 34932049 DOI: 10.1039/d1cp04490g] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
The quantitative prediction of vibrational spectra of chromophore molecules in solution is challenging and numerous methods have been developed. In this work, we present a quantum vibration perturbation (QVP) approach, which is a procedure that combines molecular quantum vibration and molecular dynamics with perturbation theory. In this framework, an initial Newtonian molecular dynamics simulation is performed, followed by a substitution process to embed molecular quantum vibrational wave functions into the trajectory. The instantaneous vibrational frequency shift at each time step is calculated using the Rayleigh-Schrödinger perturbation theory, where the perturbation operator is the difference in the vibrational potential between the reference chromophore and the perturbed chromophore in the environment. Semi-classical statistical mechanics is employed to obtain the spectral lineshape function. We validated our method using HCOOH·nH2O (n = 1-2) clusters and HCOOH aqueous solution as examples. The QVP method can be employed for rapid prediction of the vibrational spectrum of a specific mode in solution.
Collapse
Affiliation(s)
- Yang Cong
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, 2519 Jiefang Road, Changchun 130023, People's Republic of China.
| | - Yu Zhai
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, 2519 Jiefang Road, Changchun 130023, People's Republic of China.
| | - Jitai Yang
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, 2519 Jiefang Road, Changchun 130023, People's Republic of China.
| | - Adam Grofe
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, 2519 Jiefang Road, Changchun 130023, People's Republic of China.
| | - Jiali Gao
- Department of Chemistry and Supercomputing Institute, University of Minnesota, 207 Pleasant Street, SE, Minneapolis, MN 55455, USA. .,Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China
| | - Hui Li
- Institute of Theoretical Chemistry, College of Chemistry, Jilin University, 2519 Jiefang Road, Changchun 130023, People's Republic of China.
| |
Collapse
|
10
|
Fingerhut BP, Schauss J, Kundu A, Elsaesser T. Contact pairs of RNA with magnesium ions-electrostatics beyond the Poisson-Boltzmann equation. Biophys J 2021; 120:5322-5332. [PMID: 34715079 PMCID: PMC8715182 DOI: 10.1016/j.bpj.2021.10.029] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2020] [Revised: 08/24/2021] [Accepted: 10/22/2021] [Indexed: 11/18/2022] Open
Abstract
The electrostatic interaction of RNA with its aqueous environment is most relevant for defining macromolecular structure and biological function. The attractive interaction of phosphate groups in the RNA backbone with ions in the water environment leads to the accumulation of positively charged ions in the first few hydration layers around RNA. Electrostatics of this ion atmosphere and the resulting ion concentration profiles have been described by solutions of the nonlinear Poisson-Boltzmann equation and atomistic molecular dynamics (MD) simulations. Much less is known on contact pairs of RNA phosphate groups with ions at the RNA surface, regarding their abundance, molecular geometry, and role in defining RNA structure. Here, we present a combined theoretical and experimental study of interactions of a short RNA duplex with magnesium (Mg2+) ions. MD simulations covering a microsecond time range give detailed hydration geometries as well as electrostatics and spatial arrangements of phosphate-Mg2+ pairs, including both pairs in direct contact and separated by a single water layer. The theoretical predictions are benchmarked by linear infrared absorption and nonlinear two-dimensional infrared spectra of the asymmetric phosphate stretch vibration which probes both local interaction geometries and electric fields. Contact pairs of phosphate groups and Mg2+ ions are identified via their impact on the vibrational frequency position and line shape. A quantitative analysis of infrared spectra for a range of Mg2+-excess concentrations and comparison with fluorescence titration measurements shows that on average 20-30% of the Mg2+ ions interacting with the RNA duplex form contact pairs. The experimental and MD results are in good agreement. In contrast, calculations based on the nonlinear Poisson-Boltzmann equation fail in describing the ion arrangement, molecular electrostatic potential, and local electric field strengths correctly. Our results underline the importance of local electric field mapping and molecular-level simulations to correctly account for the electrostatics at the RNA-water interface.
Collapse
|
11
|
Fingerhut BP. The mutual interactions of RNA, counterions and water - quantifying the electrostatics at the phosphate-water interface. Chem Commun (Camb) 2021; 57:12880-12897. [PMID: 34816825 PMCID: PMC8640580 DOI: 10.1039/d1cc05367a] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2021] [Accepted: 11/01/2021] [Indexed: 11/25/2022]
Abstract
The structure and dynamics of polyanionic biomolecules, like RNA, are decisively determined by their electric interactions with the water molecules and the counterions in the environment. The solvation dynamics of the biomolecules involves a subtle balance of non-covalent and many-body interactions with structural fluctuations due to thermal motion occurring in a femto- to subnanosecond time range. This complex fluctuating many particle scenario is crucial in defining the properties of biological interfaces with far reaching significance for the folding of RNA structures and for facilitating RNA-protein interactions. Given the inherent complexity, suited model systems, carefully calibrated and benchmarked by experiments, are required to quantify the relevant interactions of RNA with the aqueous environment. In this feature article we summarize our recent progress in the understanding of the electrostatics at the biological interface of double stranded RNA (dsRNA) and transfer RNA (tRNA). Dimethyl phosphate (DMP) is introduced as a viable and rigorously accessible model system allowing the interaction strength with water molecules and counterions, their relevant fluctuation timescales and the spatial reach of interactions to be established. We find strong (up to ≈90 MV cm-1) interfacial electric fields with fluctuations extending up to ≈20 THz and demonstrate how the asymmetric stretching vibration νAS(PO2)- of the polarizable phosphate group can serve as the most sensitive probe for interfacial interactions, establishing a rigorous link between simulations and experiment. The approach allows for the direct interfacial observation of interactions of biologically relevant Mg2+ counterions with phosphate groups in contact pair geometries via the rise of a new absorption band imposed by exchange repulsion interactions at short interatomic distances. The systematic extension to RNA provides microscopic insights into the changes of the hydration structure that accompany the temperature induced melting of the dsRNA double helix and quantify the ionic interactions in the folded tRNA. The results show that pairs of negatively charged phosphate groups and Mg2+ ions represent a key structural feature of RNA embedded in water. They highlight the importance of binding motifs made of contact pairs in the electrostatic stabilization of RNA structures that have a strong impact on the surface potential and enable the fine tuning of the local electrostatic properties which are expected to be relevant for mediating the interactions between biomolecules.
Collapse
|
12
|
Gresh N, Perahia D. Multimolecular complexes of the phosphodiester anion with Zn(II) or Mg(II) and water molecules-Preliminary validations of a polarizable potential by ab initio quantum chemistry. J Comput Chem 2021; 42:1430-1446. [PMID: 34101861 DOI: 10.1002/jcc.26555] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2021] [Revised: 04/26/2021] [Accepted: 04/27/2021] [Indexed: 11/06/2022]
Abstract
Dimethyl phosphate (DMP- ) is a model for the phosphodiester backbone of DNA, RNA, and phospholipids. It is central for the binding of divalent cations and water along the backbone of nucleic acids. Significant polarization and charge-transfer contributions and nonadditivity come into play in the multimolecular complexes organized around phosphate. Prior to large-scale molecular dynamics (MD) with advanced polarizable potentials, it is essential to evaluate how well the values and trends of intermolecular interaction energies (ΔE) from ab initio quantum chemistry (QC) and their individual contributions are reproduced in a diversity of such complexes. These differ by the starting binding modes of a divalent cation, Zn(II), namely direct, bi- or mono-dentate to anionic and/or ester oxygens, versus through-water binding. We present first the results from automated refinements of the individual contributions of the SIBFA potential with respect to their QC counterparts using a Zn(II) or a water probe. This is followed by validations on eight relaxed multimolecular complexes of DMP- with Zn(II) or Mg(II) and seven waters, then on sixteen complexes of DMP- with Zn(II) and eight waters in arrangements extracted from MD or energy-minimization on a droplet of sixty-four waters. This monitors the compared evolutions of SIBFA and QC ΔE and their individual contributions in the competing arrangements. Some waters, bridging Zn(II) and DMP- , were found to have exceptionally large dipole moments, of up to 3.8 Debye. The perspectives of extension to a flexible phosphodiester backbone are discussed in the context of the SIBFA potential for DNA and RNA.
Collapse
Affiliation(s)
- Nohad Gresh
- Laboratoire de Chimie Théorique, UMR 7616 CNRS, Sorbonne Université, Paris, France
| | - David Perahia
- Laboratoire de Biologie et Pharmacologie Appliquées, UMR 8113 CNRS, Ecole Normale Supérieure Paris-Saclay, Gif-sur-Yvette, France
| |
Collapse
|
13
|
Elsaesser T, Schauss J, Kundu A, Fingerhut BP. Phosphate Vibrations Probe Electric Fields in Hydrated Biomolecules: Spectroscopy, Dynamics, and Interactions. J Phys Chem B 2021; 125:3899-3908. [PMID: 33834783 PMCID: PMC8154594 DOI: 10.1021/acs.jpcb.1c01502] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
Electric interactions
have a strong impact on the structure and
dynamics of biomolecules in their native water environment. Given
the variety of water arrangements in hydration shells and the femto-
to subnanosecond time range of structural fluctuations, there is a
strong quest for sensitive noninvasive probes of local electric fields.
The stretching vibrations of phosphate groups, in particular the asymmetric
(PO2)− stretching vibration νAS(PO2)−, allow for a quantitative
mapping of dynamic electric fields in aqueous environments via a field-induced
redshift of their transition frequencies and concomitant changes of
vibrational line shapes. We present a systematic study of νAS(PO2)− excitations in molecular
systems of increasing complexity, including dimethyl phosphate (DMP),
short DNA and RNA duplex structures, and transfer RNA (tRNA) in water.
A combination of linear infrared absorption, two-dimensional infrared
(2D-IR) spectroscopy, and molecular dynamics (MD) simulations gives
quantitative insight in electric-field tuning rates of vibrational
frequencies, electric field and fluctuation amplitudes, and molecular
interaction geometries. Beyond neat water environments, the formation
of contact ion pairs of phosphate groups with Mg2+ ions
is demonstrated via frequency upshifts of the νAS(PO2)− vibration, resulting in a distinct
vibrational band. The frequency positions of contact geometries are
determined by an interplay of attractive electric and repulsive exchange
interactions.
Collapse
Affiliation(s)
- Thomas Elsaesser
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, Berlin 12489, Germany
| | - Jakob Schauss
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, Berlin 12489, Germany
| | - Achintya Kundu
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, Berlin 12489, Germany
| | - Benjamin P Fingerhut
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, Berlin 12489, Germany
| |
Collapse
|
14
|
Schauss J, Kundu A, Fingerhut BP, Elsaesser T. Magnesium Contact Ions Stabilize the Tertiary Structure of Transfer RNA: Electrostatics Mapped by Two-Dimensional Infrared Spectra and Theoretical Simulations. J Phys Chem B 2021; 125:740-747. [PMID: 33284610 PMCID: PMC7848891 DOI: 10.1021/acs.jpcb.0c08966] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
![]()
Ions interacting with hydrated RNA
play a central role in defining
its secondary and tertiary structure. While spatial arrangements of
ions, water molecules, and phosphate groups have been inferred from
X-ray studies, the role of electrostatic and other noncovalent interactions
in stabilizing compact folded RNA structures is not fully understood
at the molecular level. Here, we demonstrate that contact ion pairs
of magnesium (Mg2+) and phosphate groups embedded in local
water shells stabilize the tertiary equilibrium structure of transfer
RNA (tRNA). Employing dialyzed tRNAPhe from yeast and tRNA
from Escherichia coli, we follow the
population of Mg2+ sites close to phosphate groups of the
ribose-phosphodiester backbone step by step, combining linear and
nonlinear infrared spectroscopy of phosphate vibrations with molecular
dynamics simulations and ab initio vibrational frequency calculations.
The formation of up to six Mg2+/phosphate contact pairs
per tRNA and local field-induced reorientations of water molecules
balance the phosphate–phosphate repulsion in nonhelical parts
of tRNA, thus stabilizing the folded structure electrostatically.
Such geometries display limited sub-picosecond fluctuations in the
arrangement of water molecules and ion residence times longer than
1 μs. At higher Mg2+ excess, the number of contact
ion pairs per tRNA saturates around 6 and weakly interacting ions
prevail. Our results suggest a predominance of contact ion pairs over
long-range coupling of the ion atmosphere and the biomolecule in defining
and stabilizing the tertiary structure of tRNA.
Collapse
Affiliation(s)
- Jakob Schauss
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, Berlin 12489, Germany
| | - Achintya Kundu
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, Berlin 12489, Germany
| | - Benjamin P Fingerhut
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, Berlin 12489, Germany
| | - Thomas Elsaesser
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie, Berlin 12489, Germany
| |
Collapse
|
15
|
Baiz CR, Błasiak B, Bredenbeck J, Cho M, Choi JH, Corcelli SA, Dijkstra AG, Feng CJ, Garrett-Roe S, Ge NH, Hanson-Heine MWD, Hirst JD, Jansen TLC, Kwac K, Kubarych KJ, Londergan CH, Maekawa H, Reppert M, Saito S, Roy S, Skinner JL, Stock G, Straub JE, Thielges MC, Tominaga K, Tokmakoff A, Torii H, Wang L, Webb LJ, Zanni MT. Vibrational Spectroscopic Map, Vibrational Spectroscopy, and Intermolecular Interaction. Chem Rev 2020; 120:7152-7218. [PMID: 32598850 PMCID: PMC7710120 DOI: 10.1021/acs.chemrev.9b00813] [Citation(s) in RCA: 171] [Impact Index Per Article: 42.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Vibrational spectroscopy is an essential tool in chemical analyses, biological assays, and studies of functional materials. Over the past decade, various coherent nonlinear vibrational spectroscopic techniques have been developed and enabled researchers to study time-correlations of the fluctuating frequencies that are directly related to solute-solvent dynamics, dynamical changes in molecular conformations and local electrostatic environments, chemical and biochemical reactions, protein structural dynamics and functions, characteristic processes of functional materials, and so on. In order to gain incisive and quantitative information on the local electrostatic environment, molecular conformation, protein structure and interprotein contacts, ligand binding kinetics, and electric and optical properties of functional materials, a variety of vibrational probes have been developed and site-specifically incorporated into molecular, biological, and material systems for time-resolved vibrational spectroscopic investigation. However, still, an all-encompassing theory that describes the vibrational solvatochromism, electrochromism, and dynamic fluctuation of vibrational frequencies has not been completely established mainly due to the intrinsic complexity of intermolecular interactions in condensed phases. In particular, the amount of data obtained from the linear and nonlinear vibrational spectroscopic experiments has been rapidly increasing, but the lack of a quantitative method to interpret these measurements has been one major obstacle in broadening the applications of these methods. Among various theoretical models, one of the most successful approaches is a semiempirical model generally referred to as the vibrational spectroscopic map that is based on a rigorous theory of intermolecular interactions. Recently, genetic algorithm, neural network, and machine learning approaches have been applied to the development of vibrational solvatochromism theory. In this review, we provide comprehensive descriptions of the theoretical foundation and various examples showing its extraordinary successes in the interpretations of experimental observations. In addition, a brief introduction to a newly created repository Web site (http://frequencymap.org) for vibrational spectroscopic maps is presented. We anticipate that a combination of the vibrational frequency map approach and state-of-the-art multidimensional vibrational spectroscopy will be one of the most fruitful ways to study the structure and dynamics of chemical, biological, and functional molecular systems in the future.
Collapse
Affiliation(s)
- Carlos R. Baiz
- Department of Chemistry, University of Texas at Austin, Austin, TX 78712, U.S.A
| | - Bartosz Błasiak
- Department of Physical and Quantum Chemistry, Wroclaw University of Science and Technology, Wybrzeże Wyspiańskiego 27, 50-370 Wrocław, Poland
| | - Jens Bredenbeck
- Johann Wolfgang Goethe-University, Institute of Biophysics, Max-von-Laue-Str. 1, 60438, Frankfurt am Main, Germany
| | - Minhaeng Cho
- Center for Molecular Spectroscopy and Dynamics, Seoul 02841, Republic of Korea
- Department of Chemistry, Korea University, Seoul 02841, Republic of Korea
| | - Jun-Ho Choi
- Department of Chemistry, Gwangju Institute of Science and Technology, Gwangju 61005, Republic of Korea
| | - Steven A. Corcelli
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, IN 46556, U.S.A
| | - Arend G. Dijkstra
- School of Chemistry and School of Physics and Astronomy, University of Leeds, Leeds LS2 9JT, U.K
| | - Chi-Jui Feng
- Department of Chemistry, James Franck Institute and Institute for Biophysical Dynamics, University of Chicago, Chicago, IL 60637, U.S.A
| | - Sean Garrett-Roe
- Department of Chemistry, University of Pittsburgh, Pittsburgh, PA 15260, U.S.A
| | - Nien-Hui Ge
- Department of Chemistry, University of California at Irvine, Irvine, CA 92697-2025, U.S.A
| | - Magnus W. D. Hanson-Heine
- School of Chemistry, University of Nottingham, Nottingham, University Park, Nottingham, NG7 2RD, U.K
| | - Jonathan D. Hirst
- School of Chemistry, University of Nottingham, Nottingham, University Park, Nottingham, NG7 2RD, U.K
| | - Thomas L. C. Jansen
- University of Groningen, Zernike Institute for Advanced Materials, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - Kijeong Kwac
- Center for Molecular Spectroscopy and Dynamics, Seoul 02841, Republic of Korea
| | - Kevin J. Kubarych
- Department of Chemistry, University of Michigan, 930 N. University Ave., Ann Arbor, MI 48109, U.S.A
| | - Casey H. Londergan
- Department of Chemistry, Haverford College, Haverford, Pennsylvania 19041, U.S.A
| | - Hiroaki Maekawa
- Department of Chemistry, University of California at Irvine, Irvine, CA 92697-2025, U.S.A
| | - Mike Reppert
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, Ontario M5S 3H6, Canada
| | - Shinji Saito
- Department of Theoretical and Computational Molecular Science, Institute for Molecular Science, Myodaiji, Okazaki, 444-8585, Japan
| | - Santanu Roy
- Chemical Sciences Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6110, U.S.A
| | - James L. Skinner
- Institute for Molecular Engineering, University of Chicago, Chicago, IL 60637, U.S.A
| | - Gerhard Stock
- Biomolecular Dynamics, Institute of Physics, Albert Ludwigs University, 79104 Freiburg, Germany
| | - John E. Straub
- Department of Chemistry, Boston University, Boston, MA 02215, U.S.A
| | - Megan C. Thielges
- Department of Chemistry, Indiana University, 800 East Kirkwood, Bloomington, Indiana 47405, U.S.A
| | - Keisuke Tominaga
- Molecular Photoscience Research Center, Kobe University, Nada, Kobe 657-0013, Japan
| | - Andrei Tokmakoff
- Department of Chemistry, James Franck Institute and Institute for Biophysical Dynamics, University of Chicago, Chicago, IL 60637, U.S.A
| | - Hajime Torii
- Department of Applied Chemistry and Biochemical Engineering, Faculty of Engineering, and Department of Optoelectronics and Nanostructure Science, Graduate School of Science and Technology, Shizuoka University, 3-5-1 Johoku, Naka-Ku, Hamamatsu 432-8561, Japan
| | - Lu Wang
- Department of Chemistry and Chemical Biology, Institute for Quantitative Biomedicine, Rutgers University, 174 Frelinghuysen Road, Piscataway, NJ 08854, U.S.A
| | - Lauren J. Webb
- Department of Chemistry, The University of Texas at Austin, 105 E. 24th Street, STOP A5300, Austin, Texas 78712, U.S.A
| | - Martin T. Zanni
- Department of Chemistry, University of Wisconsin-Madison, Madison, Wisconsin 53706-1396, U.S.A
| |
Collapse
|
16
|
Aqueous Contact Ion Pairs of Phosphate Groups with Na+, Ca2+ and Mg2+ – Structural Discrimination by Femtosecond Infrared Spectroscopy and Molecular Dynamics Simulations. Z PHYS CHEM 2020. [DOI: 10.1515/zpch-2020-1614] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Abstract
The extent of contact and solvent shared ion pairs of phosphate groups with Na+, Ca2+ and Mg2+ ions in aqueous environment and their relevance for the stability of polyanionic DNA and RNA structures is highly debated. Employing the asymmetric phosphate stretching vibration of dimethyl phosphate (DMP), a model system of the sugar-phosphate backbone of DNA and RNA, we present linear infrared, femtosecond infrared pump-probe and absorptive 2D-IR spectra that report on contact ion pair formation via the presence of blue shifted spectral signatures. Compared to the linear infrared spectra, the nonlinear spectra reveal contact ion pairs with increased sensitivity because the spectra accentuate differences in peak frequency, transition dipole moment strength, and excited state lifetime. The experimental results are corroborated by long time scale MD simulations, benchmarked by density functional simulations on phosphate-ion-water clusters. The microscopic interpretation reveals subtle structural differences of ion pairs formed by the phosphate group and the ions Na+, Ca2+ and Mg2+. Intricate properties of the solvation shell around the phosphate group and the ion are essential to explain the experimental observations. The present work addresses a challenging to probe topic with the help of a model system and establishes new experimental data of contact ion pair formation, thereby underlining the potential of nonlinear 2D-IR spectroscopy as an analytical probe of phosphate-ion interactions in complex biological systems.
Collapse
|