Calculation of binary magnetic properties and potential energy curve in xenon dimer: second virial coefficient of (129)Xe nuclear shielding

Quantum chemical calculations of the nuclear shielding tensor, the nuclear quadrupole coupling tensor, and the spin-rotation tensor are reported for the Xe dimer using ab initio quantum chemical methods. The binary chemical shift delta, the anisotropy of the shielding tensor Delta sigma, the nuclear quadrupole coupling tensor component along the internuclear axis chi( parallel ), and the spin-rotation constant C( perpendicular ) are presented as a function of internuclear distance. The basis set superposition error is approximately corrected for by using the counterpoise correction (CP) method. Electron correlation effects are systematically studied via the Hartree-Fock, complete active space self-consistent field, second-order M?ller-Plesset many-body perturbation, and coupled-cluster singles and doubles (CCSD) theories, the last one without and with noniterative triples, at the nonrelativistic all-electron level. We also report a high-quality theoretical interatomic potential for the Xe dimer, gained using the relativistic effective potential/core polarization potential scheme. These calculations used valence basis set of cc-pVQZ quality supplemented with a set of midbond functions. The second virial coefficient of Xe nuclear shielding, which is probably the experimentally best-characterized intermolecular interaction effect in nuclear magnetic resonance spectroscopy, is computed as a function of temperature, and compared to experiment and earlier theoretical results. The best results for the second virial coefficient, obtained using the CCSD(CP) binary chemical shift curve and either our best theoretical potential or the empirical potentials from the literature, are in good agreement with experiment. Zero-point vibrational corrections of delta, Delta sigma, chi (parallel), and C (perpendicular) in the nu=0, J=0 rovibrational ground state of the xenon dimer are also reported.     

Evaluation and optimization of coherence transfer in high molecular weight systems

For very large proteins in the highest magnetic fields, the large chemical shift anisotropy (CSA) of carbonyl carbon deteriorates coherence transfer efficiency in experiments designed for unambiguous sequential backbone assignment. In this communication, coherence throughput of several TROSY experiments is evaluated. Two new experiments, MP-HNCA and HN(CO)CANH, are also introduced as attractive alternatives for sequential assignment purposes of large proteins with correlation time over 50 ns. Their theoretical coherence transfer efficiencies for the interresidual (13)C(alpha) correlations are significantly better than in recently introduced MP-CT-HNCA and sequential HNCA experiments. The improvement with the new experiments is observed already on 60.8 kDa homodimer of protein Cel6A at 800 (1)H MHz.     

Optimized pathway selection in intraresidual triple-resonance experiments

An optimized intraresidual pulse sequence element with better sensitivity and suppression of sequential cross peaks is presented. Concatenation of three magnetization transfer delays allows their independent setting, in accordance with the relaxation properties of the individual spins, without concomitantly prolonging the pulse sequence. Additionally, implementations of the scheme to HNCA, HNCACB, and the TROSY based triple-resonance experiments are proposed. The feasibility of the new element was verified by recording HNCA and HNCACB on the small 8.6 kDa protein ubiquitin. The corresponding HNCA-TROSY experiment was tested on a larger protein, the 30.4 kDa Cel6A from the thermophilic soil bacterium Thermobifida fusca at 800 (1)H MHz.     

Measurement of Homonuclear J-Couplings from Spin-State Selective Double-/Zero-Quantum Two-Dimensional NMR Spectra

¹H-detected two-dimensional double-/zero-quantum experiments are described for measurement of homonuclear ²J_HH-couplings of NH₂ or CH₂ groups in proteins. These experiments utilize multiple-quantum coherence for determination of the size and the absolute sign of the geminal scalar and dipolar couplings in the presence of broad lines. Spectra are simplified by gradient selection and spin-state selective filters.     

Measurement of homonuclear (2)J-couplings from spin-state selective double-/zero-quantum two-dimensional NMR spectra.

(1)H-detected two-dimensional double-/zero-quantum experiments are described for measurement of homonuclear (2)J(HH)-couplings of NH(2) or CH(2) groups in proteins. These experiments utilize multiple-quantum coherence for determination of the size and the absolute sign of the geminal scalar and dipolar couplings in the presence of broad lines. Spectra are simplified by gradient selection and spin-state selective filters.     

Toward calculations of the 129Xe chemical shift in Xe@C60 at experimental conditions: relativity, correlation, and dynamics

We calculate the 129Xe chemical shift in endohedral Xe@C60 with systematic inclusion of the contributing physical effects to model the real experimental conditions. These are relativistic effects, electron correlation, the temperature-dependent dynamics, and solvent effects. The ultimate task is to obtain the right result for the right reason and to develop a physically justified methodological model for calculations and simulations of endohedral Xe fullerenes and other confined Xe systems. We use the smaller Xe...C6H6 model to calibrate density functional theory approaches against accurate correlated wave function methods. Relativistic effects as well as the coupling of relativity and electron correlation are evaluated using the leading-order Breit-Pauli perturbation theory. The dynamic effects are treated in two ways. In the first approximation, quantum dynamics of the Xe atom in a rigid cage takes advantage of the centrosymmetric potential for Xe within the thermally accessible distance range from the center of the cage. This reduces the problem of obtaining the solution of a diatomic rovibrational problem. In the second approach, first-principles classical molecular dynamics on the density functional potential energy hypersurface is used to produce the dynamical trajectory for the whole system, including the dynamic cage. Snapshots from the trajectory are used for calculations of the dynamic contribution to the absorption 129Xe chemical shift. The calculated nonrelativistic Xe shift is found to be highly sensitive to the optimized molecular structure and to the choice of the exchange-correlation functional. Relativistic and dynamic effects are significant and represent each about 10% of the nonrelativistic static shift at the minimum structure. While the role of the Xe dynamics inside of the rigid cage is negligible, the cage dynamics turns out to be responsible for most of the dynamical correction to the 129Xe shift. Solvent effects evaluated with a polarized continuum model are found to be very small.     

Nuclear magnetic resonance chemical shifts and quadrupole couplings for different hydrogen-bonding cases occurring in liquid water: a computational study

Nuclear magnetic resonance (NMR) parameters are determined theoretically for the oxygen and hydrogen/deuterium nuclei of differently hydrogen-bonded water molecules in liquid water at 300 K. The parameters are the chemical shift, the shielding anisotropy, the asymmetry parameter of shielding, the nuclear quadrupole coupling constant, and the asymmetry parameter of the nuclear quadrupole coupling. We sample instantaneous configurations from a Car-Parrinello molecular dynamics simulation and feed nuclear coordinates into a quantum chemical program for the calculation of NMR parameters using density-functional theory with the three-parameter hybrid exchange-correlation (B3LYP) functional. In the subsequent analysis, molecules are divided into groups according to the number of hydrogen bonds they possess, and the full average NMR tensors are calculated separately for each group. The classification of the hydrogen-bonding cases is performed using a simple distance-based criterion. The analysis reveals in detail how the NMR tensors evolve as the environment changes gradually from gas to liquid upon increasing the number of hydrogen bonds to the molecule of interest. Liquid-state distributions of the instantaneous values of the NMR properties show a wide range of values for each hydrogen-bonding species with significant overlap between the different cases. Our study shows how local changes in the environment, along with classical thermal averaging, affect the NMR parameters in liquid water. For example, a broken or alternatively extra hydrogen bond induces major changes in the NMR tensors, and the effect is more pronounced for hydrogen or deuterium than for oxygen. The data sheds light on the usefulness of NMR experiments in investigating the local coordination of liquid water.     

SESAME-HSQC for simultaneous measurement of NH and CH scalar and residual dipolar couplings

We present a novel pulse sequence, SESAME-HSQC, for the simultaneous measurement of several NH and CH scalar and residual dipolar couplings in double labeled proteins. The proposed Spin-statE Selective All Multiplicity Edited (SESAME)-HSQC combines gradient selected and sensitivity enhanced (15)N- and constant-time (13)C-HSQC experiments with the recently introduced spin-state selective method (Nolis et al., J. Magn. Reson. 180 (2006) 39-50) for measuring couplings simultaneously at amide and aliphatic regions. Excellent resolution and high sensitivity is warranted by removing all coupling interactions during the indirectly detected t(1) period, and by employing pulsed field gradients for coherence selection and utilizing coherence order selective spin-state selection. The scalar and residual dipolar couplings can be readily measured from a two-dimensional (15)N/(13)C-HSQC spectrum without additional spectral crowding. SESAME-HSQC can be used for epitope mapping by observing chemical shift changes in both amide and aliphatic regions. Simultaneously, potential conversion in protein conformation can be probed by analyzing changes in residual dipolar couplings induced by ligand binding. The pulse sequence is experimentally verified with a sample of (15)N/(13)C enriched human ubiquitin. The internuclear vector directions determined from the residual dipolar couplings are found to be in excellent correlation with those predicted from ubiquitin's refined solution structure.     

Calculations of nuclear quadrupole coupling in noble gas-noble metal fluorides: interplay of relativistic and electron correlation effects

The nuclear quadrupole coupling constants (NQCCs) of noble gas and noble metal nuclei in the recently found noble gas-noble metal fluorides (NgMF, where Ng=Ar,Kr,Xe and M=Cu,Ag,Au) are obtained theoretically by high-level ab initio calculations, where both relativistic and electron correlation effects are included, and compared to experimental results. Fully relativistic four-component Dirac-Hartree-Fock (DHF) calculations are carried out at the basis set limit for electric field gradient that couples with the electric quadrupole moment of the nucleus, and uncorrelated relativistic effects are extracted by comparing DHF results to nonrelativistic (NR) HF calculations. Electron correlation effects are investigated both at fully relativistic second-order Moller-Plesset (DMP2) and at NR MP2 levels of theory, as well as at the NR coupled-cluster singles and doubles with perturbational triples [CCSD(T)] level. The validity of the approximation where relativistic effects, on the one hand, and nonrelativistically obtained correlation effects, on the other hand, are evaluated separately and assumed to be additive, is investigated by comparison with the DMP2 results. Inclusion of relativistic effects is shown to be necessary for obtaining the correct NQCC trends as the nucleus of interest and/or its neighbors become heavier. Electron correlation treatment is needed for approaching quantitative agreement with the experimental NQCCs. The assumption of additive electron correlation and relativistic effects, corresponding to the NR correlation treatment added on top of relativistic DHF data, gives qualitatively correct noble gas NQCCs. For noble metal NQCCs, correlation treatment at the relativistic level of theory is mandatory for reaching agreement with experimental results. Current work also confirms the experimental trends of NQCCs, which have been taken as an indication of nearly covalent interaction between noble gas and noble metal in the heaviest present systems, especially in XeAuF.