Dependence Of Chemical Shifts Research Articles

Isotopomeric Elucidation of the Mechanism of Temperature Sensitivity in 59Co NMR Molecular Thermometers.

Understanding the mechanisms governing temperature-dependent magnetic resonance properties is essential for enabling thermometry via magnetic resonance imaging. Herein we harness a new molecular design strategy for thermometry─that of effective mass engineering via deuteration in the first coordination shell─to reveal the mechanistic origin of 59Co chemical shift thermometry. Exposure of [Co(en)3]3+ (1; en = ethylenediamine) and [Co(diNOsar)]3+ (2; diNOsar = dinitrosarcophagine) to mixtures of H2O and D2O produces distributions of [Co(en)3]3+-dn (n = 0-12) and [Co(diNOsar)]3+-dn (n = 0-6) isotopomers all resolvable by 59Co NMR. Variable-temperature 59Co NMR analyses reveal a temperature dependence of the 59Co chemical shift, Δδ/ΔT, on deuteration of the N-donor atoms. For 1, deuteration amplifies Δδ/ΔT by 0.07 ppm/°C. Increasing degrees of deuteration yield an opposing influence on 2, diminishing Δδ/ΔT by -0.07 ppm/°C. Solution-phase Raman spectroscopy in the low-frequency 200-600 cm-1 regime reveals a red shift of Raman-active Co-N6 vibrational modes by deuteration. Analysis of the normal vibrational modes shows that Raman modes produce the largest variation in 59Co δ. Finally, partition function analysis of the Raman-active modes shows that increased populations of Raman modes predict greater Δδ/ΔT, representing new experimental insight into the thermometry mechanism.

Analysis of NMR Adsorption Isotherms of Zeolite ZSM-5: Adsorption Profiles Derived from the Pressure and Temperature Dependences of 129Xe NMR Chemical Shift and Signal Intensity.

129Xe NMR spectroscopy of nanomaterials, such as zeolites, can provide valuable information on the nanostructure and physicochemical properties of adsorption. In the present study the pressure and temperature dependences of the 129Xe NMR chemical shift and the signal intensity were investigated in detail with a zeolite ZSM-5. The pressure dependence of the signal intensity at constant temperature was analyzed based on the Langmuir and Dubinin-Radushkevich (D-R) models, from which the thermodynamic parameters and energetic profiles of adsorption were obtained together with information concerning the nanospace size. From this isotherm analysis the coverage, θ, was calculated and used for isotherm analysis of the chemical shift. The θ dependence of the chemical shift was successfully fitted by an exponential function, and the results were discussed in relation to the chemical shift at zero coverage, that at full coverage and the curvature of the exponential function. The chemical shift data reported with the zeolites NaA and KA, where separated signals were observed for the different number of encapsulated Xe atoms in the α cage, were analyzed and discussed collectively.

Temperature-Dependent Solid-State NMR Proton Chemical-Shift Values and Hydrogen Bonding.

Temperature-dependent NMR experiments are often complicated by rather long magnetic-field equilibration times, for example, occurring upon a change of sample temperature. We demonstrate that the fast temporal stabilization of a magnetic field can be achieved by actively stabilizing the temperature of the magnet bore, which allows quantification of the weak temperature dependence of a proton chemical shift, which can be diagnostic for the presence of hydrogen bonds. Hydrogen bonding plays a central role in molecular recognition events from both fields, chemistry and biology. Their direct detection by standard structure-determination techniques, such as X-ray crystallography or cryo-electron microscopy, remains challenging due to the difficulties of approaching the required resolution, on the order of 1 Å. We, herein, explore a spectroscopic approach using solid-state NMR to identify protons engaged in hydrogen bonds and explore the measurement of proton chemical-shift temperature coefficients. Using the examples of a phosphorylated amino acid and the protein ubiquitin, we show that fast magic-angle spinning (MAS) experiments at 100 kHz yield sufficient resolution in proton-detected spectra to quantify the rather small chemical-shift changes upon temperature variations.

Temperature Dependence of 1 H Paramagnetic Chemical Shifts in Actinide Complexes, Beyond Bleaney's Theory: The AnVI O2 2+ -Dipicolinic Acid Complexes (An=Np, Pu) as an Example.

Actinide +VI complexes ( = , and ) with dipicolinic acid derivatives were synthesized and characterized by powder XRD, SQUID magnetometry and NMR spectroscopy. In addition, and complexes were described by first principles CAS based and two-component spin-restricted DFT methods. The analysis of the 1 H paramagnetic NMR chemical shifts for all protons of the ligands according to the X-rays structures shows that the Fermi contact contribution is negligible in agreement with spin density determined by unrestricted DFT. The magnetic susceptibility tensor is determined by combining SQUID, pNMR shifts and Evans' method. The SO-RASPT2 results fit well the experimental magnetic susceptibility and pNMR chemical shifts. The role of the counterions in the solid phase is pointed out; their presence impacts the magnetic properties of the complex. The temperature dependence of the pNMR chemical shifts has a strong contribution, contrarily to Bleaney's theory for lanthanide complexes. The fitting of the temperature dependence of the pNMR chemical shifts and SQUID magnetic susceptibility by a two-Kramers-doublet model for the complex and a non-Kramers-doublet model for the complex allows for the experimental evaluation of energy gaps and magnetic moments of the paramagnetic center.

Thulium complex with DOPC as 1H NMR temperature sensor

The goal of the research was to show the limits within which Tm3+ can be used as a linear response probe in lipid membrane. Temperature dependence of lanthanide-induced chemical shifts (LISs) were studied using the example of DOPC phospholipid by means of NMR spectroscopy. It was demonstrated that these dependences have thermal feature, which could be connected with the changes of the structure of lanthanide-phospholipid complex (Tm-DOPC complex). The analysis of paramagnetic properties such as the temperature sensitivity of LIS indicates the prospect of using kinetically unstable lanthanide complexes as shifting temperature-sensitive probes for determining local temperature using MRI technologies.

Thermal and Nonthermal Microwave Effects of Ethanol and Hexane-Mixed Solution as Revealed by In Situ Microwave Irradiation Nuclear Magnetic Resonance Spectroscopy and Molecular Dynamics Simulation.

Microwave heating is widely used to accelerate the organic synthesis reaction. However, the role of the nonthermal microwave effect in the chemical reaction has not yet been well characterized. The microwave heating processes of an ethanol-hexane mixed solution were investigated using in situ microwave irradiation nuclear magnetic resonance spectroscopy and molecular dynamics (MD) simulation. The temperature of the solution under microwave irradiation was estimated from the temperature dependence of the 1H chemical shifts (chemical shift calibrated (CSC)-temperature). The CSC-temperature increased to 58 °C for CH2 and CH3 protons, while it increased to 42 °C for OH protons during microwave irradiation. The CSC-temperature of CH2 and CH3 protons reflects the bulk temperature of solution by the thermal microwave effect. The lower CSC-temperature of the OH proton can be attributed to a nonthermal microwave effect. MD simulation revealed that electron dipole moments of OH groups ordered along the oscillated electric field decreased the entropy by absorbing microwave energy and simultaneously increased the entropy by dissipating energy to the solution as the thermal and nonthermal microwave effect. Ordered polar molecules interact to increase hydrogen bonds between OH groups as the nonthermal microwave effect, which explains the lower CSC-temperature of the OH protons. The nonthermal microwave effects contribute to the intrinsic acceleration of the organic reaction.

Hydration and Diffusion of H+, Li+, Na+, Cs+ Ions in Cation-Exchange Membranes Based on Polyethylene- and Sulfonated-Grafted Polystyrene Studied by NMR Technique and Ionic Conductivity Measurements.

The main particularities of sulfonate groups hydration, water molecule and alkaline metal cation translation mobility as well as ionic conductivity were revealed by NMR and impedance spectroscopy techniques. Cation-exchange membranes MSC based on cross-linked sulfonated polystyrene (PS) grafted on polyethylene with ion-exchange capacity of 2.5 mg-eq/g were investigated. Alkaline metal cation hydration numbers (h) calculated from temperature dependences of 1H chemical shift of water molecule for membranes equilibrated with water vapor at RH = 95% are 5, 6, and 4 for Li+, Na+, and Cs+ ions, respectively. These values are close to h for equimolar aqueous salt solutions. Water molecules and counter ions Li+, Na+, and Cs+ diffusion coefficients were measured by pulsed field gradient NMR on the 1H, 7Li, 23Na, and 133Cs nuclei. For membranes as well as for aqueous chloride solutions, cation diffusion coefficients increased in the following sequence: Li+ < Na+ < Cs+. Cation and water molecule diffusion activation energies in temperature range from 20 °C to 80 °C were close to each other (about 20 kJ/mol). The cation conductivity of MSC membranes is in the same sequence, Li+ < Na+ < Cs+ << H+. The conductivity values calculated from the NMR diffusion coefficients with the use of the Nernst–Einstein equation are essentially higher than experimentally determined coefficients. The reason for this discrepancy is the heterogeneity of membrane pore and channel system. Ionic conductivity is limited by cation transfer in narrow channels, whereas the diffusion coefficient characterizes ion mobility in wide pores first of all.

Temperature dependence of 1H NMR chemical shifts and diffusivity of confined ethylammonium nitrate ionic liquid

Some ionic liquids (ILs) change their dynamic properties when placed in a confinement between polar surfaces (Filippov et al., Phys. Chem. Chem. Phys. 2018, 20, 6316). The diffusivities of ions and NMR relaxation times in these ILs also reversibly change under a strong static magnetic field. The mechanisms of these phenomena are not clear, but it has been suggested that they involve modified hydrogen-bonding networks formed in these ILs in the presence of polar surfaces. To obtain a better understanding of these effects, we performed temperature-dependent measurements of chemical shifts and diffusion coefficients for ethylammonium nitrate (EAN) IL in the bulk phase (IB) and confined in layers with a thickness of ~4 μm between quartz plates unexposed (I phase) and exposed (IMF phase) to a static magnetic field of 9.4 T. It was shown that the NMR chemical shift of NH3 protons of EAN in the I phase is strongly shifted upfield, ~0.0145 ppm/K, which is due to weakening of the hydrogen-bonding network of the confined EAN. Exposure to the magnetic field leads to restitution of the hydrogen-bonding (H-bonding network). The temperature dependences of diffusion coefficients follow the order D(I) > D(IB) > D(IMF) and can be described by a Vogel-Fulcher-Tammann approach with variation of the pre-exponential factor, which is determined by the strength of the H-bonding network. Confinement of EAN between plates (IB → I) is an endothermic process, while processes occurring in a magnetic field, I → IMF and IMF → I, are exothermic and endothermic, respectively.

Quantum mechanical reference spectrum simulation for precursors and degradation products of chemicals relevant to the Chemical Weapons Convention.

A selection of acidic, alkaline and neutral degradation products relevant to the Chemical Weapons Convention was studied in wide range of pH conditions to determine their spin systems as well as spectral parameters. The pH dependence of chemical shifts and J couplings was parameterized using Henderson-Hasselbalch-based functions using dichloromethane as additional shift reference in TSP-d4 referenced spectra. The resulting parameters allowed calculation of precise chemical shifts and J coupling constants in arbitrary pH conditions. The validity of the obtained spin system definitions and parameters as a source of quantum mechanically simulated reference data in chemical verification analysis is demonstrated.

Investigation of solvated calcium dication structure in pure water, methanol, and ethanol solutions by means of K and L2,3-edges X-ray absorption spectroscopy

The interaction of an ion with surrounding solvent molecules is revealed by X-ray absorption spectroscopy in transmission mode for the case of calcium dication dissolved in pure water, methanol and ethanol at the Ca K and L2,3-edges. The near-threshold K-shell excitation leads to two distinct solvent-dependent features, the Ca 1s to 3d and 4p electronic transitions, clearly indicating the ligand field decreases as water is changed to methanol and methanol to ethanol. On the other hand, the Ca L2,3-edge excitation leads to small changes in the absorption spectra and the peak assignment – mostly in terms of Ca 2p to eg and t2g transitions for each of the two spin orbit components – is strongly depending on the geometry of the first solvation shell around the dication and on the magnitude of the ligand-field multiplet effects. Our interpretation is supported by ab-initio quantum chemical calculations of the solvent dependence of chemical shifts for the Ca 3d and 4p unoccupied levels of the 1s and 2p excited states, performed on a small cluster including the nearest neighbor solvent molecules surrounding the Ca dication assuming a cubic (8-fold) solvent coordination.

Pulsed Third-Spin-Assisted Recoupling NMR for Obtaining Long-Range 13C-13C and 15N-13C Distance Restraints.

We present a class of pulsed third-spin-assisted recoupling (P-TSAR) magic-angle-spinning solid-state NMR techniques that achieve efficient polarization transfer over long distances to provide important restraints for structure determination. These experiments utilize second-order cross terms between strong 1H-13C and 1H-15N dipolar couplings to achieve 13C-13C and 15N-13C polarization transfer, similar to the principle of continuous-wave (CW) TSAR experiments. However, in contrast to the CW-TSAR experiments, these P-TSAR experiments require much less radiofrequency (rf) energy and allow a much simpler routine for optimizing the rf field strength. We call the technique PULSAR (pulsed proton-assisted recoupling) for homonuclear spin pairs. For heteronuclear spin pairs, we improve the recently introduced PERSPIRATIONCP (proton-enhanced rotor-echo short pulse irradiation cross-polarization) experiment by shifting the pulse positions and removing the z-filters, which significantly broaden the bandwidth and increase the efficiency of polarization transfer. We demonstrate the PULSAR and PERSPIRATIONCP techniques on the model protein GB1 and found cross peaks for distances as long as 10 and 8 Å for 13C-13C and 15N-13C spin pairs, respectively. We then apply these methods to the amyloid fibrils formed by the peptide hormone glucagon and show that long-range correlation peaks are readily observed to constrain intermolecular packing in this cross-β fibril. We provide an analytical model for the PULSAR and PERSPIRATIONCP experiments to explain the measured and simulated chemical shift dependence and pulse flip angle dependence of polarization transfer. These two techniques are useful for measuring long-range distance restraints to determine the three-dimensional structures of proteins and other biological macromolecules.

Conformational heterogeneity of Savinase from NMR, HDX-MS and X-ray diffraction analysis.

BackgroundSeveral examples have emerged of enzymes where slow conformational changes are of key importance for function and where low populated conformations in the resting enzyme resemble the conformations of intermediate states in the catalytic process. Previous work on the subtilisin protease, Savinase, from Bacillus lentus by NMR spectroscopy suggested that this enzyme undergoes slow conformational dynamics around the substrate binding site. However, the functional importance of such dynamics is unknown.MethodsHere we have probed the conformational heterogeneity in Savinase by following the temperature dependent chemical shift changes. In addition, we have measured changes in the local stability of the enzyme when the inhibitor phenylmethylsulfonyl fluoride is bound using hydrogen-deuterium exchange mass spectrometry (HDX-MS). Finally, we have used X-ray crystallography to compare electron densities collected at cryogenic and ambient temperatures and searched for possible low populated alternative conformations in the crystals.ResultsThe NMR temperature titration shows that Savinase is most flexible around the active site, but no distinct alternative states could be identified. The HDX shows that modification of Savinase with inhibitor has very little impact on the stability of hydrogen bonds and solvent accessibility of the backbone. The most pronounced structural heterogeneities detected in the diffraction data are limited to alternative side-chain rotamers and a short peptide segment that has an alternative main-chain conformation in the crystal at cryo conditions. Collectively, our data show that there is very little structural heterogeneity in the resting state of Savinase and hence that Savinase does not rely on conformational selection to drive the catalytic process.

Pressure dependence of side chain 1H and 15N-chemical shifts in the model peptides Ac-Gly-Gly-Xxx-Ala-NH2

For interpreting the pressure induced shifts of resonance lines of folded as well as unfolded proteins the availability of data from well-defined model systems is indispensable. Here, we report the pressure dependence of 1H and 15N chemical shifts of the side chain atoms in the protected tetrapeptides Ac-Gly-Gly-Xxx-Ala-NH2 (Xxx is one of the 20 canonical amino acids) measured at 800 MHz proton frequency. As observed earlier for other nuclei the chemical shifts of the side chain nuclei have a nonlinear dependence on pressure in the range from 0.1 to 200 MPa. The pressure response is described by a second degree polynomial with the pressure coefficients B1 and B2 that are dependent on the atom type and type of amino acid studied. A number of resonances could be assigned stereospecifically including the 1H and 15N resonances of the guanidine group of arginine. In addition, stereoselectively isotope labeled SAIL amino acids were used to support the stereochemical assignments. The random-coil pressure coefficients are also dependent on the neighbor in the sequence as an analysis of the data shows. For Hα and HN correction factors for different amino acids were derived. In addition, a simple correction of compression effects in thermodynamic analysis of structural transitions in proteins was derived on the basis of random-coil pressure coefficients.

Hydrogen-bonded complexes protocatechualdehyde - acetone in carbon tetrachloride: NMR-spectroscopy and molecular dynamics simulation

Abstract A comprehensive study of the hydrogen-bonded complexes of protocatechualdehyde – acetone formed in an inert solvent was carried out. Concentration dependences of the 1H NMR chemical shifts of protocatechualdehyde in carbon tetrachloride – acetone‑d6 mixtures were obtained. The most probable structures of hydrogen-bonded protocatechualdehyde-acetone complexes were determined by means of all-atomic molecular dynamics simulation. Using the GIAO/DFT (Gauge-Independent Atomic Orbital/Density Functional Theory) approach, the proton chemical shifts of the complexes were calculated within the framework of the (если она не единственная такая, то a) continuum model of the solvent, as well as with an explicit solvent. The parameters obtained from quantum-chemical calculations were used to interpret the experimentally observed dependences. It is shown that protocatechualdehyde forms predominantly hydrogen-bonded complexes with a composition of 1:1 in a binary solvent. The hydrogen atom of the hydroxyl group in the para-position takes part in the formation of the intermolecular hydrogen bond with acetone, while the hydrogen atom of the hydroxyl group in the meta-position is involved in the formation of an intramolecular hydrogen bond.

NMR “Crystallography” of Membrane Proteins Aligned in Native-Like Bilayers

Oriented-sample (OS) solid-state NMR can be used to determine structures of membrane proteins without crystallization and use of cryogenic temperatures. Here, the orientationally dependent dipolar couplings and chemical shifts provide direct input for structure calculations. However, so far only the 1H-15N dipolar couplings and 15N chemical shifts have been routinely assessed in 15N labeled protein samples. While highly efficient for determining the tilt angles of short α-helical peptides, these two measurements alone are likely insufficient for determining protein structures of arbitrary topology. Therefore, determination of tertiary structures of membrane proteins by OS NMR requires measurements of additional angular restraints. Here we have developed a new experimental triple-resonance NMR technique, which was applied to uniformly doubly (15N, 13C) labeled Pf1 coat protein reconstituted in magnetically aligned DMPC/DHPC bicelles. The previously inaccessible 1Hα-13Cα dipolar couplings, which represent powerful chiral angular restraints, have been measured for the resolved NMR resonances of the transmembrane domain of Pf1. Inclusion of the third restraint makes it possible to determine the torsion angles Φ and Ψ between the adjacent peptide planes without assuming α-helical structure a priori. Using our recently developed structure calculation algorithm [1], simultaneous fitting of three restraints per peptide plane has resulted in families of structures which have subsequently been filtered using various Rosetta scoring functions including the energetics within the membrane environment, hydrogen bonding, and potential steric clashes. Such structural validation protocol has led to a consensus α-helical transmembrane structure for Pf1 coat protein, albeit obtained in a de novo fashion, thus greatly enhancing utility of OS NMR as a “crystallography” in near-native membrane environments. [1]. Lapin, J. and Nevzorov A.A., Validation of protein backbone structures calculated from NMR angular restraints using Rosetta, J. Biomol. NMR, 2019; 73(5):229-244.

Increasing field strength versus advanced isotope labeling for NMR-based fluxomics.

Nuclear magnetic resonance (NMR)-based fluxomics seeks to measure the incorporation of isotope labels in selected metabolites to follow kinetically the synthesis of the latter. It can however equally be used to understand the biosynthetic origin of the same metabolites. We investigate here different NMR approaches to optimize such experiments in terms of resolution and time requirement. Using the isoleucine biosynthesis as an example, we explore the use of different field strengths ranging from 500 MHz to 1.1 GHz. Because of the different field dependence of chemical shift and heteronuclear J couplings, the spectra change at different field strengths. We equally explore the approach to silence the leucine/valine methyl signals through the use of a suitable deuterated precursor, thereby allowing selective observation of the Ile 13 C labeling pattern. Combining both approaches, we arrive at an efficient procedure for the NMR-based exploration of Ile biosynthesis.

Temperature dependence of NMR chemical shifts: Tracking and statistical analysis.

Isotropic chemical shifts measured by solution nuclear magnetic resonance (NMR) spectroscopy offer extensive insights into protein structure and dynamics. Temperature dependences add a valuable dimension; notably, the temperature dependences of amide proton chemical shifts are valuable probes of hydrogen bonding, temperature‐dependent loss of structure, and exchange between distinct protein conformations. Accordingly, their uses include structural analysis of both folded and disordered proteins, and determination of the effects of mutations, binding, or solution conditions on protein energetics. Fundamentally, these temperature dependences result from changes in the local magnetic environments of nuclei, but correlations with global thermodynamic parameters measured via calorimetric methods have been observed. Although the temperature dependences of amide proton and nitrogen chemical shifts are often well approximated by a linear model, deviations from linearity are also observed and may be interpreted as evidence of fast exchange between distinct conformational states. Here, we describe computational methods, accessible via the Shift‐T web server, including an automated tracking algorithm that propagates initial (single temperature) 1H—15N cross peak assignments to spectra collected over a range of temperatures. Amide proton and nitrogen temperature coefficients (slopes determined by fitting chemical shift vs. temperature data to a linear model) are subsequently calculated. Also included are methods for the detection of systematic, statistically significant deviation from linearity (curvature) in the temperature dependences of amide proton chemical shifts. The use and utility of these methods are illustrated by example, and the Shift‐T web server is freely available at http://meieringlab.uwaterloo.ca/shiftt.

Prototropic behavior of cyclohexane substituted urethane and urea compounds. Observation of H-bond mediated 4HJH1H3 coupling constants across urea fragments

Two 1-aryl-3-cyclohexylurea and two 1-aryl-3-cyclohexylurethane with and without alkyl tail in aryl fragments were synthesized and their variable-temperature 1H NMR spectra in chloroform, DMSO and DMF were recorded. The temperature dependences of chemical shift of NH protons for all compounds have been observed. The type of dependence has been explained by the aggregation of molecules and formation of ionic structures. The formation of intermolecular complexes leads to possible formation of symmetrical N1 … … H ……N3 intramolecular hydrogen bond in urea fragment. As a result,4hJH1H3 trans-hydrogen bond scalar coupling constants have been observed for the first time in low molecular weight compounds using a simple one-dimensional 1H NMR experiment. The formation of strong intramolecular H-bond leads to π-conjugation not only with spacer but with phenyl ring too. It support the Gilli conception that RAHB formation is a result of π-electron delocalization. The presence in urea fragment of such kind interaction leads to formation of ionic structure, which has been detected by NMR and UV spectroscopies. The formation of ionic structure can explain the catalytic activity of such compounds and the mechanism of transformation in organic and bioorganic reactions in which involved the urea compounds.

High pressure response of 1H NMR chemical shifts of purine nucleotides

The study of the pressure response by NMR spectroscopy provides information on the thermodynamics of conformational equilibria in proteins and nucleic acids. For obtaining a database for expected pressure effects on free nucleotides and nucleotides bound in macromolecular complexes, the pressure response of 1H chemical shifts and J-coupling constants of the purine 5′-ribonucleotides AMP, ADP, ATP, GMP, GDP, and GTP were studied in the absence and presence of Mg2+-ions. Experiments are supported by quantum-chemical calculations of populations and chemical shift differences in order to corroborate structural interpretations and to estimate missing data for AMP. The preference of the ribose S puckering obtained from the analysis of the experimental J-couplings is also confirmed by the calculations. In addition, the pressure response of the non-hydrolysable GTP analogues GppNHp, GppCH2p, and GTPγS was examined within a pressure range up to 200 MPa. As observed earlier for 31P NMR chemical shifts of these nucleotides the pressure dependence of chemical shifts is clearly non-linear in most cases. In di- and tri-phospho nucleosides, the resonances of the two protons bound to the ribose 5′ carbon are non-equivalent and can be observed separately. The gg-rotamer at C4′- C5′ bond is strongly preferred and the downfield shifted resonance can be assigned to the H5″ proton in the nucleotides. In contrast, in adenosine itself the frequencies of the two resonances are interchanged.

Modeling pH-Dependent NMR Chemical Shift Perturbations in Peptides

Modeling the pH dependence of protein and peptide chemical shifts outside the range of physiological values (6.5–7) is key to understanding structure-function relationships of these systems. These capabilities are largely not available in current chemical shift prediction software. In this study, we utilize a combination of molecular dynamics and quantum mechanics to investigate the through-space and through-bond contributions of protonation-dependent chemical shift perturbations (CSPs) in model tripeptides. By altering the protonation state of the titratable group in the tripeptides, we observe a notable difference in the conformational ensembles and attendantly compute significant CSPs for all nuclei near the site of protonation. We thus demonstrate the ability to recapitulate experimental pH-dependent CSPs with good agreement (R = 0.85, 0.99, and 0.98 for 13C, 15N, and 1H, respectively). Broadly, we provide the groundwork for incorporating pH effects into empirical and semiempirical chemical shift predictors.