1
|
Paschek D, Busch J, Mock E, Ludwig R, Strate A. Computing the frequency-dependent NMR relaxation of 1H nuclei in liquid water. J Chem Phys 2024; 160:074102. [PMID: 38364003 DOI: 10.1063/5.0191052] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2023] [Accepted: 01/23/2024] [Indexed: 02/18/2024] Open
Abstract
We present a computational framework for reliably determining the frequency-dependent intermolecular and intramolecular nuclear magnetic resonance (NMR) dipole-dipole relaxation rates of spin 1/2 nuclei from Molecular Dynamics (MD) simulations. This approach avoids the alterations caused by the well-known finite-size effects of translational diffusion. Moreover, a procedure is derived to control and correct for effects caused by fixed distance-sampling cutoffs and periodic boundary conditions. By construction, this approach is capable of accurately predicting the correct low-frequency scaling behavior of the intermolecular NMR dipole-dipole relaxation rate and thus allows for the reliable calculation of the frequency-dependent relaxation rate over many orders of magnitude. Our approach is based on the utilization of the theory of Hwang and Freed for the intermolecular dipole-dipole correlation function and its corresponding spectral density [L.-P. Hwang and J. H. Freed, J. Chem. Phys. 63, 4017-4025 (1975)] and its combination with data from MD simulations. The deviations from the Hwang and Freed theory caused by periodic boundary conditions and sampling distance cutoffs are quantified by means of random walker Monte Carlo simulations. An expression based on the Hwang and Freed theory is also suggested for correcting those effects. As a proof of principle, our approach is demonstrated by computing the frequency-dependent intermolecular and intramolecular dipolar NMR relaxation rates of 1H nuclei in liquid water at 273 and 298 K based on the simulations of the TIP4P/2005 model. Our calculations are suggesting that the intermolecular contribution to the 1H NMR relaxation rate of the TIP4P/2005 model in the extreme narrowing limit has previously been substantially underestimated.
Collapse
Affiliation(s)
- Dietmar Paschek
- Institut für Chemie, Abteilung Physikalische und Theoretische Chemie, Universität Rostock, Albert-Einstein-Str. 27, D-18059 Rostock, Germany
| | - Johanna Busch
- Institut für Chemie, Abteilung Physikalische und Theoretische Chemie, Universität Rostock, Albert-Einstein-Str. 27, D-18059 Rostock, Germany
| | - Eduard Mock
- Institut für Chemie, Abteilung Physikalische und Theoretische Chemie, Universität Rostock, Albert-Einstein-Str. 27, D-18059 Rostock, Germany
| | - Ralf Ludwig
- Institut für Chemie, Abteilung Physikalische und Theoretische Chemie, Universität Rostock, Albert-Einstein-Str. 27, D-18059 Rostock, Germany
- Leibniz-Institut für Katalyse an der Universität Rostock e.V., Albert-Einstein-Str. 29a, D-18059 Rostock, Germany
- Department Life, Light & Matter, Universität Rostock, Albert-Einstein-Str. 25, D-18059 Rostock, Germany
| | - Anne Strate
- Institut für Chemie, Abteilung Physikalische und Theoretische Chemie, Universität Rostock, Albert-Einstein-Str. 27, D-18059 Rostock, Germany
| |
Collapse
|
2
|
Voloshin V, Smolin N, Geiger A, Winter R, Medvedev NN. Dynamics of TMAO and urea in the hydration shell of the protein SNase. Phys Chem Chem Phys 2019; 21:19469-19479. [PMID: 31461098 DOI: 10.1039/c9cp03184g] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Using all-atom molecular dynamics simulations of aqueous solutions of the globular protein SNase, the dynamic behavior of water molecules and cosolvents (trimethylamine-N-oxide (TMAO) and urea) in the hydration shell of the protein was studied for different solvent compositions. TMAO is a potent protein-stabilizing osmolyte, whereas urea is known to destabilize proteins. For molecules that are initially located in successive narrow layers at a given distance from the protein, the mean displacements and the distribution of displacements for short time intervals are calculated. For molecules that are initially located in solvation shells of a given thickness around the protein, the characteristic residence times in these shells are determined to characterize the dynamic behavior of the solvent molecules as a function of the distance to the protein. A combined consideration of these characteristics allows to reveal additional features of the dynamics of the cosolvents. It is shown that TMAO molecules leave the nearest vicinity of the protein faster than urea molecules, despite the fact that the mobility of TMAO molecules, measured by their mean displacements, is lower than that of urea. Moreover, we show that the rate of release of TMAO molecules from the hydration shell is lower in ternary (TMAO + urea + H2O) solvent mixtures than in the binary ones. This is consistent with a recent observation that the fraction of TMAO near the protein decreases in the presence of urea. From the analysis of the decay of the number of particles initially located in the region of the first peak of the distribution function of solvent molecules around the protein, we estimated that about 20 water molecules and 6-7 urea molecules stay near the protein for more than 1000 ps.
Collapse
Affiliation(s)
- Vladimir Voloshin
- Institute of Chemical Kinetics and Combustion, SB RAS, 630090 Novosibirsk, Russia.
| | - Nikolai Smolin
- Department of Cell and Molecular Physiology, Loyola University Chicago, Maywood, Illinois 60153, USA
| | - Alfons Geiger
- Physikalische Chemie, Fakultät für Chemie und Chemische Biologie, Technische Universität Dortmund, Otto-Hahn-Straße 4a, 44221 Dortmund, Germany.
| | - Roland Winter
- Physikalische Chemie, Fakultät für Chemie und Chemische Biologie, Technische Universität Dortmund, Otto-Hahn-Straße 4a, 44221 Dortmund, Germany.
| | - Nikolai N Medvedev
- Institute of Chemical Kinetics and Combustion, SB RAS, 630090 Novosibirsk, Russia. and Novosibirsk State University, 630090 Novosibirsk, Russia
| |
Collapse
|
3
|
Chemical shift extremum of 129Xe(aq) reveals details of hydrophobic solvation. Sci Rep 2018; 8:7023. [PMID: 29728689 PMCID: PMC5935698 DOI: 10.1038/s41598-018-25418-4] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2018] [Accepted: 04/20/2018] [Indexed: 11/08/2022] Open
Abstract
The 129Xe chemical shift in an aqueous solution exhibits a non-monotonic temperature dependence, featuring a maximum at 311 K. This is in contrast to most liquids, where the monotonic decrease of the shift follows that of liquid density. In particular, the shift maximum in water occurs at a higher temperature than that of the maximum density. We replicate this behaviour qualitatively via a molecular dynamics simulation and computing the 129Xe chemical shift for snapshots of the simulation trajectory. We also construct a semianalytical model, in which the Xe atom occupies a cavity constituted by a spherical water shell, consisting of an even distribution of solvent molecules. The temperature dependence of the shift is seen to result from a product of the decreasing local water density and an increasing term corresponding to the energetics of the Xe-H2O collisions. The latter moves the chemical shift maximum up in temperature, as compared to the density maximum. In water, the computed temperature of the shift maximum is found to be sensitive to both the details of the binary chemical shift function and the coordination number. This work suggests that, material parameters allowing, the maximum should be exhibited by other liquids, too.
Collapse
|
4
|
Aidas K, Ågren H, Kongsted J, Laaksonen A, Mocci F. A quantum mechanics/molecular dynamics study of electric field gradient fluctuations in the liquid phase. The case of Na+in aqueous solution. Phys Chem Chem Phys 2013; 15:1621-31. [DOI: 10.1039/c2cp41993a] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
|
5
|
Lippens G, Van Belle D, Wodak S, Jeener J. T 1 relaxation time of water from a molecular dynamics simulation. Mol Phys 2006. [DOI: 10.1080/00268979300103151] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
Affiliation(s)
- G. Lippens
- a Unité de Conformation des macromolécules Biologiques , Université Libre de Bruxelles , CP 160/16, P2, Avenue P. Héger, B-1050 , Bruxelles , Belgium
- b Laboratorium voor Genetika , Universiteit Gent , Ledeganckstraat 35, B-9000 , Ghent , Belgium
| | - D. Van Belle
- a Unité de Conformation des macromolécules Biologiques , Université Libre de Bruxelles , CP 160/16, P2, Avenue P. Héger, B-1050 , Bruxelles , Belgium
| | - S.J. Wodak
- a Unité de Conformation des macromolécules Biologiques , Université Libre de Bruxelles , CP 160/16, P2, Avenue P. Héger, B-1050 , Bruxelles , Belgium
| | - J. Jeener
- c Laboratoire de Résonance Magnétique , Université Libre de Bruxelles , CP 232, Campus Plaine, B-1050 , Bruxelles , Belgium
| |
Collapse
|
6
|
Yamaguchi T, Matsuoka T, Koda S. Mode-coupling study on the dynamics of hydrophobic hydration. J Chem Phys 2006; 120:7590-601. [PMID: 15267671 DOI: 10.1063/1.1687319] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The molecular motion of water in water-hydrophobic solute mixtures was investigated by the mode-coupling theory for molecular liquids based on the interaction-site description. When the model Lennard-Jones solute was mixed with water, both the translational and reorientational motions of solvent water become slower, in harmony with various experiments and molecular dynamics simulations. We compared the mechanism of the slowing down with that of the pressure dependence of the molecular motion of neat water [T. Yamaguchi, S.-H. Chong, and F. Hirata, J. Chem. Phys. 119, 1021 (2003)]. We found that the decrease in the solvent mobility caused by the solute can essentially be elucidated by the same mechanism: That is, the fluctuation of the number density of solvent due to the cavity formation by the solute strengthens the friction on the collective polarization through the dielectric friction mechanism: We also employed the solute molecule that is the same as solvent water except for the amount of partial charges, in order to alter the strength of the solute-solvent interaction continuously. The mobility of the solvent water was reduced both by the hydrophobic and strongly hydrophilic solutes, but it was enhanced in the intermediate case. Such a behavior was discussed in connection with the concept of positive and negative hydrations.
Collapse
Affiliation(s)
- T Yamaguchi
- Department of Molecular Design and Engineering, Graduate School of Engineering, Nagoya University, Chikusa, Nagoya, Aichi 464-8603, Japan.
| | | | | |
Collapse
|
7
|
Jameson CJ, Sears DN, Murad S. Molecular dynamics averaging of Xe chemical shifts in liquids. J Chem Phys 2004; 121:9581-92. [PMID: 15538880 DOI: 10.1063/1.1807817] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
The Xe nuclear magnetic resonance chemical shift differences that afford the discrimination between various biological environments are of current interest for biosensor applications and medical diagnostic purposes. In many such environments the Xe signal appears close to that in water. We calculate average Xe chemical shifts (relative to the free Xe atom) in solution in eleven liquids: water, isobutane, perfluoro-isobutane, n-butane, n-pentane, neopentane, perfluoroneopentane, n-hexane, n-octane, n-perfluorooctane, and perfluorooctyl bromide. The latter is a liquid used for intravenous Xe delivery. We calculate quantum mechanically the Xe shielding response in Xe-molecule van der Waals complexes, from which calculations we develop Xe (atomic site) interpolating functions that reproduce the ab initio Xe shielding response in the complex. By assuming additivity, these Xe-site shielding functions can be used to calculate the shielding for any configuration of such molecules around Xe. The averaging over configurations is done via molecular dynamics (MD). The simulations were carried out using a MD technique that one of us had developed previously for the simulation of Henry's constants of gases dissolved in liquids. It is based on separating a gaseous compartment in the MD system from the solvent using a semipermeable membrane that is permeable only to the gas molecules. We reproduce the experimental trends in the Xe chemical shifts in n-alkanes with increasing number of carbons and the large chemical shift difference between Xe in water and in perfluorooctyl bromide. We also reproduce the trend for a given solvent of decreasing Xe chemical shift with increasing temperature. We predict chemical shift differences between Xe in alkanes vs their perfluoro counterparts.
Collapse
Affiliation(s)
- Cynthia J Jameson
- Department of Chemistry, MC-111, University of Illinois at Chicago, Chicago, Illinois 60607-7061, USA
| | | | | |
Collapse
|
8
|
Mocci F, Laaksonen A, Lyubartsev A, Saba G. Molecular Dynamics Investigation of23Na NMR Relaxation in Oligomeric DNA Aqueous Solution. J Phys Chem B 2004. [DOI: 10.1021/jp047744+] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
9
|
Perng BC, Ladanyi BM. A dielectric theory of spin-lattice relaxation for nuclei with electric quadrupole moments. J Chem Phys 1998. [DOI: 10.1063/1.476606] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
10
|
KIRCHNER BARBARA, ERMAKOVA ELENA, STEINEBRUNNER GEROLD, DYSON ANTHONYJ, HUBER HANSPETER. Ab initio calculation of the NMR spin-lattice relaxation time and the diffusion coefficient of 21Ne in liquid and supercritical states. Mol Phys 1998. [DOI: 10.1080/002689798168123] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
11
|
Odelius M, Holz M, Laaksonen A. Quadrupolar Relaxation of 21Ne, 83Kr, and 131Xe Dissolved in Acetonitrile. A Molecular Dynamics Study. J Phys Chem A 1997. [DOI: 10.1021/jp972506i] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Michael Odelius
- Department of Physical Chemistry, Uppsala University, Box 532, S-751 21 Uppsala, Sweden
| | - Manfred Holz
- Institut für Physikalische Chemie, Universität Karlsruhe, D-76128 Karlsruhe, Germany
| | - Aatto Laaksonen
- Division of Physical Chemistry, Arrhenius Laboratory, University of Stockholm, S-106 91 Stockholm, Sweden
| |
Collapse
|
12
|
Kaatze U. The dielectric properties of water in its different states of interaction. J SOLUTION CHEM 1997. [DOI: 10.1007/bf02768829] [Citation(s) in RCA: 128] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
13
|
Affiliation(s)
- Michel Luhmer
- Chimie Organique, CP 165, Université Libre de Bruxelles, Av. F.D. Roosevelt, 50, 1050 Brussels, Belgium
| | - Kristin Bartik
- Chimie Organique, CP 165, Université Libre de Bruxelles, Av. F.D. Roosevelt, 50, 1050 Brussels, Belgium
| |
Collapse
|
14
|
Bagno A, Scorrano G. Site of Ionization of Polyfunctional Bases and Acids. 2. Ab Initio Electric Field Gradients at Nitrogen, Oxygen, Phosphorus, and Sulfur in Neutral and Ionized Forms. ACTA ACUST UNITED AC 1996. [DOI: 10.1021/jp9518421] [Citation(s) in RCA: 22] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Alessandro Bagno
- Centro CNR Meccanismi Reazioni Organiche, Dipartimento di Chimica Organica, Università di Padova, via Marzolo 1, 35131 Padova, Italy
| | - Gianfranco Scorrano
- Centro CNR Meccanismi Reazioni Organiche, Dipartimento di Chimica Organica, Università di Padova, via Marzolo 1, 35131 Padova, Italy
| |
Collapse
|
15
|
Weingärtner H, Haselmeier R, Holz M. Effect of Xenon upon the Dynamical Anomalies of Supercooled Water. A Test of Scaling-Law Behavior. ACTA ACUST UNITED AC 1996. [DOI: 10.1021/jp951424w] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Hermann Weingärtner
- Institut für Physikalische Chemie und Elektrochemie der Universität Karlsruhe, D-76128 Karlsruhe, Germany
| | - Ralf Haselmeier
- Institut für Physikalische Chemie und Elektrochemie der Universität Karlsruhe, D-76128 Karlsruhe, Germany
| | - Manfred Holz
- Institut für Physikalische Chemie und Elektrochemie der Universität Karlsruhe, D-76128 Karlsruhe, Germany
| |
Collapse
|
16
|
Odelius M, Laaksonen A. Molecular dynamics simulations of quadrupolar relaxation of131Xe in carbon tetrachloride, acetonitrile, and methanol. Mol Phys 1994. [DOI: 10.1080/00268979400100364] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
17
|
Sacco A, Weingärtner H, Braun BM, Holz M. Study of the structure-breaking effect in aqueous CsCl solutions based on H2O/D2O isotope effects on transport coefficients and microdynamical properties. ACTA ACUST UNITED AC 1994. [DOI: 10.1039/ft9949000849] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
18
|
|
19
|
|
20
|
|
21
|
Quist PO, Halle B. Curvature defects in a lamellar phase revealed by nuclear-spin-relaxation anisotropy. PHYSICAL REVIEW. E, STATISTICAL PHYSICS, PLASMAS, FLUIDS, AND RELATED INTERDISCIPLINARY TOPICS 1993; 47:3374-3395. [PMID: 9960390 DOI: 10.1103/physreve.47.3374] [Citation(s) in RCA: 18] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
22
|
Luhmer M, van Belle D, Reisse J, Odelius M, Kowalewski J, Laaksonen A. Magnetic relaxation of xenon‐131 dissolved in benzene. A study by molecular dynamics and Monte Carlo simulations. J Chem Phys 1993. [DOI: 10.1063/1.464273] [Citation(s) in RCA: 26] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
23
|
Eggenberger R, Gerber S, Huber H, Searles D, Welker M. Abinitiocalculation of the deuterium quadrupole coupling in liquid water. J Chem Phys 1992. [DOI: 10.1063/1.463749] [Citation(s) in RCA: 75] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
24
|
Weingärtner H, Haselmeier R, Holz M. 129Xe NMR as a new tool for studying gas diffusion in liquids: self-diffusion of xenon in water. Chem Phys Lett 1992. [DOI: 10.1016/0009-2614(92)85569-v] [Citation(s) in RCA: 22] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
25
|
Kaatze U, Pottel R. Dielectric properties of organic solute/water mixtures. Hydrophobic hydration and relaxation. J Mol Liq 1992. [DOI: 10.1016/0167-7322(92)80103-o] [Citation(s) in RCA: 33] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
26
|
|
27
|
83Kr and 131Xe NMR of the noble gases dissolved in nematic liquid crystals. Determination of quadrupole couplings and electric field gradients. ACTA ACUST UNITED AC 1991. [DOI: 10.1016/0022-2364(91)90257-t] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
28
|
|
29
|
Dejaegere A, Luhmer M, Stien ML, Reisse J. Study of NMR relaxation of xenon-131 in quadrupolar solvents. ACTA ACUST UNITED AC 1991. [DOI: 10.1016/0022-2364(91)90198-3] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
30
|
Diehl P, Jokisaari J. Nuclear magnetic relaxation of the 129 Xe and 131 Xe isotopes of xenon gas dissolved in isotropic and anisotropic liquids. ACTA ACUST UNITED AC 1990. [DOI: 10.1016/0022-2364(90)90299-o] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|