
==== Front
Phys Chem Chem Phys
Phys Chem Chem Phys
CP
PPCPFQ
Physical Chemistry Chemical Physics
1463-9076
1463-9084
The Royal Society of Chemistry

39230359
d4cp02484b
10.1039/d4cp02484b
Chemistry
Conformational dependence of chemical shifts in the proline rich region of TAU protein†
† Electronic supplementary information (ESI) available: Pdf with results and figures of all evaluated methods. See DOI: https://doi.org/10.1039/d4cp02484b

https://orcid.org/0000-0003-0466-5432
Stöckelmaier Johannes a
https://orcid.org/0000-0002-4232-2556
Oostenbrink Chris a
a Institute of Molecular Modeling and Simulation (MMS), University of Natural Resources and Life Sciences Vienna Austria chris.oostenbrink@boku.ac.at

30 8 2024
18 9 2024
30 8 2024
26 36 2385623870
21 6 2024
15 8 2024
This journal is © the Owner Societies
2024
The Royal Society of Chemistry
https://creativecommons.org/licenses/by/3.0/ This article is licensed under a Creative Commons Attribution 3.0 Unported Licence. You can use material from this article in other publications without requesting further permissions from the RSC, provided that the correct acknowledgement is given.
Nuclear magnetic resonance (NMR) is an important method for structure elucidation of proteins, as it is an easily accessible and well understood method. To characterize intrinsically disordered proteins (IDPs) using computational models it is often necessary to analyze and integrate calculated observables with measurements derived from solution NMR experiments. In this case study, we investigate whether and which chemical shifts of the proline-rich region of Tau protein (residues 210–240) offer information about the conformational state to distinguish two different microscopic conformers. Using multiple computational methods, the chemical shifts of these two conformationally distinct structures are calculated. The different methods are compared regarding their ability to compute chemical shifts that are sensitive to conformational change. The analysis of the data shows significant differences between the available methods and gives suggestions for an improved pathway for ensemble reweighting. Nevertheless, the variation in the chemical shifts which are predicted for configurations that are commonly considered to belong to the same conformation is such that this obscures a comparison between distinct conformations. Conformational sensitivity is found for up to ∼26% of calculated chemical shifts. It is found to be unrelated to the atom element and has a minor relationship with the change in the corresponding ϕ dihedral angle.

Chemical shifts of a stretched and a globular conformer of one TAU-protein fragment are calculated. The conformational sensitivity of each chemical shift is evaluated.

Österreichische Nationalstiftung für Forschung, Technologie und Entwicklung 10.13039/100010132 Unassigned Austrian Science Fund 10.13039/501100002428 I-4588 Christian Doppler Forschungsgesellschaft 10.13039/501100006012 MIB pubstatusPaginated Article
==== Body
pmc1 Introduction

For decades, the well-known structure–function paradigm had its firm place in the understanding of biochemistry.1,2 The discovery of much more flexible, disordered, proteins led to a rethinking of established theories in the early 2000s.3 Fully disordered proteins are named intrinsically disordered proteins (IDPs), while partly disordered proteins contain intrinsically disordered regions (IDRs).4

During the last two decades, IDPs were the subject of significant scientific interest as they can be physiologically active and are predicted to play a role in understanding diseases such as Alzheimer's and Parkinson's disease.5–7 It is estimated that more than one third of proteins in eukaryotic organisms feature intrinsically disordered regions.8

At room temperature, IDPs can access many different conformational states within less than one microsecond. Therefore, they need to be described by a set of structures, called the conformational ensemble.9 To calculate the conformational ensemble of an IDP, molecular dynamics simulation can be performed.10,11 Due to the flexible nature of IDPs, they usually feature flatter potential energy surface regions that may span multiple conformations,12 which makes computer simulations very sensitive to inaccuracies of the force field. To overcome these limitations, conformational ensembles can be optimized with reweighting algorithms by combining experimental and simulated data.13

A chemical shift is a measurement of the resonance frequency change of a nucleus in reference to a standard in an NMR experiment. This corresponds to a change in the magnetic shielding tensor of the atomic core in reference to a standard.14 In organic chemistry, it is well understood that the local geometry of a molecule has a strong influence on the chemical shift. In biochemistry, secondary chemical shifts are applied to differentiate between α-helix or β-strand regions in structured proteins.15–17 While empirical predictors are designed to handle such cases, it is a topic of discussion to which extent they can capture conformational dynamics.18–20 As measured chemical shifts do not represent single conformations, but only the ensemble average,21,22 it is observed that chemical shifts calculated from ensembles show increased agreement with experiments.23,24

The quality of a reweighted conformational ensemble is crucially dependent on the quality of the input data. To allow reweighting tools to yield the best results, the chosen experimental and simulated datasets should fulfill the following characteristics:

(1) The measured signal must be sensitive to the overall conformation of the sample.

(2) The expected error of measurement must be less than the expected conformational sensitivity.

(3) The measured data must not have assignment errors.

(4) The measured properties should be easily accessible both computationally and experimentally.

To study IDPs, it has become a common practice to use data from residual dipolar coupling, NOE, SAXS and FRET experiments.13 While chemical shifts are a property relatively easy to measure and understand, it is disputed if they are an appropriate observable for use with protein ensemble reweighting. Here, we aim to determine if predicted chemical shifts for individual configurations differentiate between distinct conformations. We analysed two distinct conformations of an intrinsically disordered TAU-protein fragment, by selecting five highly similar configurations for each of the conformations. We compare the variation in the predicted chemical shifts within a single conformation to the variation between conformations to determine the conformational sensitivity of the predictions. We subsequently ask the question whether conformationally sensitive chemical shifts can be predicted from the molecular properties, or if a specific algorithm is more suitable to obtain conformationally sensitive chemical shifts. While our work is not aimed to validate the predicted chemical shifts against experimental data, we do perform a comparison against the experimental data to determine if systematic errors in the predictions occur.

2 Materials and methods

2.1 Molecular dynamics simulation

A 31-amino-acid long fragment of TAU-protein (aa220–aa240) as described in Lasorsa et al.25,26 was chosen as a test case. The fragment was simulated with acetylated and N-methylated N- and C termini, respectively, at 310 K for 500 ns using the AMBER99SB-ILDN force field27 with PME calculation of nonbonded interactions using the OpenMM molecular dynamics engine.28 We have used the capped ends to reflect the fact that this is just a fragment from a larger protein. The fragment was solvated using 1.5 nm padding in explicit OPC water.29 The LangevinMiddle integrator using a 2 fs time-step was used together with SHAKE constraints on the bond lengths. A MonteCarlo barostat was used to keep the pressure at 1 atm, creating an NpT-ensemble.

Analysis of the radius of gyration of the backbone showed reversible collapse and extension during the entire time span (Fig. 1). From the trajectory, both the frame with the maximal radius of gyration (RGYR) as well as the frame with minimal RGYR were selected, representing the most stretched and the most globular conformation. We selected two of the most differing conformations from the entire simulation to increase the probability of observing large differences in structurally sensitive properties. For each of the two structures, the four most similar conformations were also selected, summing up to two groups of five conformations (Fig. 1). In reference to the corresponding central structure, the maximum RMSD of the globular ensemble is ∼1.4 Å, while it is ∼2.0 Å for the stretched ensemble.

Fig. 1 (a) A set of five structures representing the most stretched conformation (blue) and a set of five representing the most globular conformation (red). (b) The radius of gyration of the polypeptide shows reversible fluctuations between ∼9.7 Å and ∼27.5 Å.

2.2 Geometry optimization

As it is known that chemical shifts are very sensitive to changes in the local chemical environment, the selected frames obtained from molecular dynamics were geometry optimized using MOPAC v22.1.0 utilizing the MOZYME protocol and the PM7 semiempirical method.30 For each frame, the 31 amino acid long polypeptide plus the two residues of the N- and C-termini as well as the first solvation shell of water were selected. The geometry of all peptide atoms plus termini was optimized, while the position of the water molecules in the solvation shell was frozen to prevent the macroscopic conformation of the peptide from changing.

2.3 Chemical shift prediction

Many different methods to predict chemical shifts have been developed in the past decades. Some of those have been selected to be investigated in this work and are listed in Table 1. These methods can be separated into two groups – empirical models using statistical approaches and DFT-based calculations using first principles methods. The class of DFT-based methods applying the gauge-independent atomic orbital (GIAO) theory can be further split using pure quantum (QM) and combined quantum/molecular-mechanics (QM/MM) approaches. Further subdivision can be made into subgroups depending on the solvation method.

Overview of the tested methods to predict chemical shifts from static protein structures. The column 'Solvation' describes how solvation effects are included into the calculation of the chemical shifts. Calculations that do not take solvation effects into account (solvation 'None') are called vacuum calculations in this article

Method	Type	Solvation	Annotations/settings	Ref.	
DFT-based (NWChem)	QM	None	2 functional/basis-set combinations.	31 and 32	
DFT-based (NWChem)	QM	Implicit	12 functional/basis-set combinations.	31 and 32	
DFT-based (NWChem)	QM	Explicit	3 functional/basis-set combinations.	31 and 32	
DFT-based (NWChem)	QM/MM	None	2 functional/basis-set combinations.	31 and 32	
DFT-based (NWChem)	QM/MM	Implicit	12 functional/basis-set combinations.	31 and 32	
SHIFTX2	Empirical	Empirical		33	
SPARTA+	Empirical	Empirical	-first 2 -last 32	34	
UCBShift-X	Empirical	Empirical	–shiftx_only –pH 6.5	35	
PPM	Empirical	Empirical	-model ann	36	

After the geometry optimization, all 2 × 5 conformations of the peptide were saved as PDB files and then converted into the corresponding input format of each method using a MDAnalysis software37,38 based custom toolchain.

For use with the DFT-based QM algorithm, the polypeptide had to be preprocessed as the system would otherwise be too large. The polypeptide was split into 33 smaller fragments. Each resulting fragment contains a central amino-acid as well as all other amino-acids and end-groups which feature at least one atom within 4 Å of the central amino-acid. At the cutoff, where the backbone of the protein was cut, the new fragments were missing atoms to replace the broken bonds. To fix the fragments, they were saturated with ACE/NME end-groups.

Implicit solvent simulations applied the continuum solvation model COSMO (York–Karplus formulation39) as implemented in NWChem. The dielectric constant of the medium was set to 78.4 to represent water. Explicit solvent simulations were set-up to feature one solvation shell around the central amino-acid; itself embedded in implicit solvation. Fig. 2 shows a visualization of the setup for one of the DFT-based QM calculations. The micro-solvation contained all water molecules within 4 Å of the central fragment and was taken from the MD-simulation with preserved geometry and orientation of the water molecules.

Fig. 2 The QM approach splits the entire polypeptide into smaller fragments. One amino-acid is selected for property calculation which then constitutes the center of the created fragment (shown in color). The environment is modeled by all amino-acids with at least one atom closer than 4 Å to the central residue and, in the case of a calculation with explicit solvent, the first solvation shell of water molecules. The visualization shows the setup for calculation of the 22nd amino acid of the polypeptide in its globular form with explicit solvation.

In the case of the DFT-based QM/MM approach, the fragmentation was handled by the NWChem software. Similar to the pure QM approach, all amino-acids close to the central amino-acid were included into the QM region of the QM/MM calculation. The QM-region includes the point charges of the MM-region according to the modified AMBER95 force field as implemented in NWChem. Broken bonds between the QM and the MM region were automatically repaired using hydrogen link atoms. In both vacuum and implicit solvent calculations, the interaction zone between the QM-region and the MM-region as well as the MM-cutoff was set to half the box-size. Fig. 3 shows a visualization of the setup for one of the QM/MM calculations.

Fig. 3 The QM/MM approach splits the molecule into two regions. The QM-region is composed of the central residue (shown in color) and all residues with at least one atom closer than 4 Å to the central residue (shown in brown). The rest of the protein is modeled as the MM-region, shown as gray spheres. Shielding properties are calculated for the atoms of the central residue. The visualization shows the setup for calculation of the 22nd amino acid of the polypeptide in its globular form in vacuum.

To evaluate the influence of functional/basis-set combinations, three basis-sets (6-31G*,40–42 cc-pdvz,43 and pcSseg-144) and four functionals (B3LYP,45 Becke97-2,46 Becke97-D47 and wb97x-d348) were tested. Empirical GD3 dispersion49 has been applied where available. The resulting 12 combinations were used with implicit solvent both by the pure QM and with the QM/MM approach. Taking the observations from these calculations into account, a smaller sub-set of three combinations was chosen to test the influence of vacuum and explicit solvent on the results.

To obtain a full set of chemical shifts for the entire polypeptide, the chemical shifts of all 33 central residues were combined. While the empirical methods yield chemical shifts directly, the DFT-based methods yield isotropic nuclear magnetic shieldings. To obtain the chemical shifts, the magnetic shieldings were referenced against standards as recommended in the study by Pavlíková Pecechtlová et al.50 The chemical shifts of 1H and 13C were referenced against tetramethylsilane (TMS) (eqn (1)) while the 15N chemical shifts were referenced against methylamine (eqn (2)) using a secondary standard referencing scheme, where δcalcX is the chemical shift of atom X and σcalcX is the isotropic nuclear magnetic shielding of atom X. The values of (244.6 ppm,51) and (249.5 ppm,52) were taken from the literature. All DFT-based methods were referenced against standards using implicit solvent (TMS) and vacuum simulations (CH3NH2).1 δcalcX = σcalcTMS − σcalcX

2

Empirical predictors, which have been trained to predict chemical shifts directly from a three-dimensional structure, are also available. As these models were trained on structured proteins it remains to be tested if they are able to calculate chemical shifts that show conformational sensitivity toward IDPs.

The four empirical predictors work in a straight forward way and take PDBs of the entire polypeptide as input. The water molecules and ions of the system were removed before calculation of the chemical shifts, which are calculated directly without need for a reference calculation. As these predictors are only trained for specific tasks, they are not able to reproduce the chemical shifts of all atoms. Most chemical shifts are calculated with SHIFTX2, which includes the backbone and most of the sidechains (401 atoms in total). The PPM software calculates the chemical shifts of the backbone, Cβ atoms and most of the sidechain hydrogen atoms (317 atoms in total). Both SPARTA+ and UCBShiftX yield chemical shifts only for the backbone and Cβ atoms (170 atoms).

2.4 Analysis

To analyze the sensitivity of chemical shifts with regards to conformational change, the five stretched and the five globular conformations of the polypeptide were considered to be equivalent in both cases. Thus, the simulation is regarded as five independent measurements of both the stretched and globular conformation.

For each evaluated atom, five chemical-shift values for both the stretched and globular conformation were calculated, respectively. Using the two times five samples of the observable, two Gaussian shaped probability distributions were obtained (Fig. 4a). It was assumed that the standard deviations of the two distributions were similar, thus allowing the calculation of a pooled standard deviation. Each of the two distributions has one expectation value and the difference between them, in multiples of the pooled standard deviation, is regarded as the sensitivity of the chemical shift to changes in the overall conformation. The chemical shift is regarded as conformationally sensitive if the overlap of the two probability distributions is less than 10%, which is equal to a difference in expectation value of 3.29σ. This process was repeated with all tested atoms (Fig. 4b).

Fig. 4 The central methodology of this work can be described using the three plots presented above. All calculations shown in these figures were calculated using the UCBShiftX method. (a) The chemical shift of the Cα atom of residue 217 was calculated for each of the 10 conformations. The result is displayed as blue and red dots, depending on whether the conformation was stretched or globular. Dots of the same color are part of the same probability distribution (PDF), shown as a dotted curve. The pooled standard distribution of the two PDFs is noted on the lower x-axis. The value zero is set to be in the middle of the two peaks. The distance between the two expectation values, in multiples of pooled σ, is used as a metric of the conformational sensitivity and quality criteria of the prediction methods. An overlap of 10% is equal to a distance of 3.29σ. The chemical shift yielded from Cα, residue 217, has a conformational sensitivity of 3.0σ and thus is just not conformationally sensitive. (b) Average chemical shifts of both the stretched (blue) and globular (red) population of all 31 Cα atoms can be seen. The vertical lines display a confidence interval of ±2σ. (c) The conformational sensitivities of all Cα atoms are displayed as a box plot. The majority of Cα atoms show a conformational sensitivity of below 3.29σ. Some shifts are more sensitive, with the highest sensitivity at just over 8σ.

The conformational sensitivities were categorized – regarding their atom of origin – into seven groups: The backbone atoms Cα, Ccarbonyl, Hamide and Namide as well as atoms from side-chains Cβ, Cother, Hother and Nother. For each of the groups, the conformational sensitivity (difference in chemical shift expectation value counted in pooled σ) is displayed as a box plot in Fig. 4c. The cutoff of 3.29σ is presented as a dashed horizontal red line. Chemical shifts which originated from the capped end groups as well as from oxygen atoms were excluded from further analysis as it is uncommon for them to be measured experimentally and are thus rarely used for ensemble reweighting.

In addition to the chemical sensitivity, the agreement with experimental data was evaluated to assess if any systematic biases are observed. If chemical shifts predicted for any conformation are systematically under- or over-estimated with respect to the experiment, a reweighing of the ensemble becomes difficult. The experimental chemical shifts were obtained from the study by Lasorsa et al.25 in which they were measured at 600 MHz, 5 °C and pH 7.3. Backbone assignment was performed using triple-resonance solution state NMR experiments. Sidechains were assigned using additional 3D NOESY and TOCSY experiments. The mean relative error (eqn (3)) was chosen as a metric of compliance with the experiment.3

where, is the mean relative error of the tested method, Oexpi is the experimental measurement of the observable i, and Osimi,mean is the simulated mean of observable i.

To rule out a significant share of non-random coil secondary structure in the tested polypeptide, DSSP analysis53 using MDAnalysis as well as the secondary chemical shifts were analyzed. The secondary chemical shift is defined as the difference between the random coil chemical shift and the measured chemical shift of the same atom (Δδ = δobs − δrc).54 To follow the more modern convention used in key publications,16,17,55 the sign of the equation has been inverted compared to in the study by Dalgarno et al. Due to the contrasting behavior of Cα and Cβ secondary chemical shifts with respect to stable secondary structures,17 they can be subtracted56 to create the secondary structure identifier ΔΔδαβ = ΔδCα − ΔδCβ. A positive ΔΔδαβ indicates an α-structure, while a negative ΔΔδαβ indicates a β-structure.17,26 To calculate the random coil chemical shifts the POTENCI software57 was used. The settings were chosen as pH = 7.3, T = 293 K and ionic strength = 0.13 mol L−1 to represent the experimental conditions as closely as possible.

3 Results

3.1 Conformational sensitivity

The overall conformational sensitivity of all approaches is summarized in Fig. 5. The sensitivity in terms of mean chemical shift difference between the stretched and globular conformations (in measures of Δσ) is shown in Fig. 5a while Fig. 5b shows the percentage of chemical shifts that is regarded as conformationally sensitive. While both metrics confirm that the overall sensitivity of chemical shifts with respect to conformation is rather limited, the empirical predictors UCBShiftX, SPARTA+ and PPM were found to be the most sensitive. The SHIFTX2 predictor performed equally to the DFT-based methods calculated in vacuum. DFT-based calculations, both QM and QM/MM, with implicit solvation differentiated very little depending on the choice of functional and basis-set. The introduction of micro-solvation as an explicit solvent model removes most of the remaining conformational sensitivity.

Fig. 5 Measurements of the mean conformational sensitivity (a) as well as the share of chemical shifts above the threshold (b) show that three empirical methods (PPM, SPARTA+ and UCBShiftX) are more sensitive to conformational change compared to DFT-based methods. It has to be noted, that there is no experimental reference and the true value is unknown. It is expected, that a higher value represents a more sensitive model.

To reweight a conformational ensemble into agreement with experimental data, it is required that the observable of interest is sensitive to the overall molecular conformation. While reweighting algorithms may be able to disregard some insensitive data that does not give information about the conformation, we emphasize that a major share of chemical shifts in these calculations seems insensitive. It may, hence, be necessary to find a sub-set of chemical shifts that contains most of the information in regard to the conformational state. We investigated whether the four empirical methods and the DFT-based predictions with a selected functional/basis-set combination can find agreement regarding the conformational sensitivity of chemical shifts. Fig. 6 shows whether the methods regard individual chemical shifts as conformationally sensitive or not. A green mark represents a sensitive chemical shift, while a red one is expected not to be sensitive. Out of the total 497 chemical shifts, there are 42 cases with at least four of nine predictions agreeing on the chemical shift to be conformational sensitive. Another two sensitive cases are chemical shifts of oxygen atoms, which are not considered for further analysis.

Fig. 6 The matrix shows whether a chemical shift is predicted to be conformationally sensitive by each of the nine compared methods. A green mark (○) represents a conformational sensitive chemical shift while a red mark identifies a non-sensitive shift. If a chemical shift has not been calculated, it is marked as white. Chemical shifts with at least four out of the nine methods agreeing on conformational sensitivity are shown. A complete overview of all 497 shifts can be found in the ESI,† Chapter 1.2.

3.1.1 Relationship of atom category and conformational sensitivity

Fig. 7 shows the conformational sensitivity of the seven atom-categories. Visual assessment of the plots generated from DFT-based calculations with implicit solvent and from vacuum simulations, as well as half of the empirical predictors, seem to show minor increased median sensitivity on the amide-nitrogen and the accompanying proton compared to the other atom categories. Another common theme of the DFT-based calculation is increased maximum sensitivities of some side-chain atoms. The choice of functional and basis leads to only few noteworthy differences in conformational sensitivity (comparisons are given in the ESI,† Chapter 1.4).

Fig. 7 Each column shows the conformational sensitivity of an atom category. On the x-axis, seven atom categories can be viewed, to check whether some atoms are more prone to sensitive chemical shifts than others. The y-axis shows the conformational sensitivity as described in Fig. 4a. The orange line of each box represents the median value while the lower and upper edges represent the end of the first and third quartiles, respectively. Data points outside of that range are displayed as dots with the largest being annotated for each atom category. The left panel shows a representative result from a DFT-based calculation while the right panel shows results obtained with the empirical method PPM.

The empirical models PPM, SPARTA+ and SHIFTX2 agree on the slightly increased sensitivity of the amide-hydrogen (Fig. 7 (right) and Fig. S4.3 and S4.4, ESI†). In contrast, UCBShiftX evaluates Cα and the amide-nitrogen atoms to contain the most conformational information (Fig. S4.2, ESI†). Most empirical chemical shift predictors show less data as they are not designed to calculate the chemical shifts of all side chain atoms. Therefore, the effect of the overall conformation on the chemical shifts of sidechains cannot be observed.

3.1.2 Relationship of chemical shift and backbone torsion

It can be hypothesized that the change in the chemical shift is related to changes in the ψ and ϕ angles of the protein backbone. To check this hypothesis, the conformational sensitivities of the residues were plotted against the change in backbone torsion Δψ and Δϕ in Fig. 8.

Fig. 8 The left and middle columns show Ramachandran plots for both the stretched and globular conformation, respectively. The right column shows the change in ψ and ϕ angles when switching from the globular to the stretched conformation. The chemical sensitivity is represented in the color of the dots. Conformational sensitivities in this figure were obtained using the DFT-based QM/MM method in vacuum with the wb97x-d3/cc-pvdz theory. The Ramachandran background was plotted using data from ref. 58 and 59.

As the conformational sensitivity is a per-atom metric, it has to be transformed into a per-residue metric according to eqn (4). To calculate the average sensitivity of the residue, the sensitivities of all backbone atoms and Cβ were averaged.4

where <svg xmlns="http://www.w3.org/2000/svg" version="1.0" width="13.666667pt" height="16.000000pt" viewBox="0 0 13.666667 16.000000" preserveAspectRatio="xMidYMid meet"><metadata> Created by potrace 1.16, written by Peter Selinger 2001-2019 </metadata><g transform="translate(1.000000,15.000000) scale(0.014583,-0.014583)" fill="currentColor" stroke="none"><path d="M160 920 l0 -40 -40 0 -40 0 0 -40 0 -40 -40 0 -40 0 0 -80 0 -80 40 0 40 0 0 -80 0 -80 40 0 40 0 0 -40 0 -40 80 0 80 0 0 -40 0 -40 80 0 80 0 0 -120 0 -120 -120 0 -120 0 0 40 0 40 -40 0 -40 0 0 40 0 40 -40 0 -40 0 0 -120 0 -120 280 0 280 0 0 40 0 40 40 0 40 0 0 160 0 160 -40 0 -40 0 0 80 0 80 -80 0 -80 0 0 40 0 40 -80 0 -80 0 0 40 0 40 -40 0 -40 0 0 40 0 40 160 0 160 0 0 -40 0 -40 80 0 80 0 0 80 0 80 -40 0 -40 0 0 40 0 40 -240 0 -240 0 0 -40z m0 -160 l0 -40 40 0 40 0 0 -80 0 -80 120 0 120 0 0 -40 0 -40 40 0 40 0 0 -40 0 -40 40 0 40 0 0 -120 0 -120 -40 0 -40 0 0 80 0 80 -40 0 -40 0 0 40 0 40 -80 0 -80 0 0 40 0 40 -80 0 -80 0 0 80 0 80 -40 0 -40 0 0 80 0 80 40 0 40 0 0 -40z"/></g></svg> res is the conformational sensitivity of the residue res, i is the conformational sensitivity of atom i, and N is the number of single sensitivities averaged into res.

If, instead of mean sensitivities, the influence of Δψ and Δϕ on specific atoms is of interest, the residue sensitivity can be set equal to the sensitivity of that specific atom (res = i). Plots for Ccarb, Cα, Cβ, Hamide and Namide atoms as well as the remaining methods can be found in the ESI,† Chapter 1.5.

In the case of the DFT-based QM/MM method with wb97x-d3/cc-pvdz theory, residue 237 shows the highest sensitivity and has also a significant change in ψ angle. Two more residues show sensitivity and have a significant change in ϕ angle. On the other hand, there are plenty of residues that show changes in the angles but no sensitivity.

3.1.3 Random forest analysis

In the previous sections the visual analysis of the influence of dihedral angles and atom category on the conformational sensitivity was assessed to be very minor with unclear statistical significance. To evaluate methodically if there is any feature that is related with a chemical shift being conformational sensitive or not, a permutation feature selection has been performed. Table 2 shows both geometrical and biochemical features that can be expected to influence the chemical sensitivity. As label, the conformational sensitivity was chosen.

Overview of the features tested to evaluate their influence on the conformational sensitivity of chemical shifts. The atom of interest is the atom for which the conformational sensitivity is evaluated

Feature	Annotations	
Atom category	Categorical value according to the groups in Fig. 7.	
Atom name	Categorical value according to the name of the atom.	
Atom element	Categorical value according to the element of the atom.	
Residue name	Categorical value according to PDB residue name.	
Residue number	Number of the associated residue in the peptide sequence.	
Δ-distance atom (*)	Change in distance from the atom of interest to the next closest atom of type (*) associated with a residue at least three amino acids away from the selected atom. Atom (*) can either be oxygen (O), nitrogen (N) or the center of geometry of the residue (RES).	
Δψ and Δϕ angle	Change in backbone dihedral angle.	
ΔSASA	Change in the solvent accessible surface areas per atom between stretched and globular conformation. Calculated with FreeSASA.60	
Is sidechain?	Categorical value if the atom is part of the backbone or sidechain.	
Per-atom alignability	The five equal samples of both the stretched and globular conformation are aligned. The per-atom alignability measures the mean self distance of the atoms to their own copies.	
Negative control (continuous)	Uniform random number (float) between 0.0 and 20.0.	
Negative control (categorical)	Categorical random number (int) between 0 and 6.	

To train the random forest regressor,61 categorical features like atom-category and residue-name had to be converted into representative, numerical dummies using one-hot or ordinal encoding. While one-hot-encoding is expected to yield higher quality results, it is difficult to reverse the encoding to obtain the importance of whole feature categories. To evaluate whole features, for example the importance of the amino acid type, ordinal encoding was applied. To evaluate the importance of single elements within a feature, to answer for example if proline or valine amino acids are related with conformational sensitivity, one-hot encoding was used.

The random forest was trained on the features, using a 80/20 split between training and testing data, 100 trees and squared error as a measurement of the split quality. To prevent overfitting on the training data, min_samples_leaf was set to 10, min_samples_split to 15 and max_depth to 10. To measure the quality of the regression, artificial control features were added in addition to the geometrical and biophysical features. The two negative controls were always unrelated to the conformational sensitivity and were added to confirm the validity of the method. For each regression, the R2-score was calculated. After building the model, single features were randomly shuffled to calculate a permutation feature importance. If a feature has influence on the model, shuffling it will reduce the R2-score significantly. The features that yield the biggest loss if shuffled, are regarded as the most important. The process was repeated 50 times with different seeding of the random number generator to create a set of models. Single trainings, which yielded a negative coefficient of determination (R2) on the testing data, were removed from the set. The remaining feature importance values and model scores of the set were then averaged. Overall, only weak relationships were found with testing R2-scores up to 0.12 in the case of DFT-based calculations (Fig. S6.23, ESI†).

Feature importances obtained from DFT-based calculations with different functional/basis-set combinations show comparable results. Fig. 9 shows feature importances obtained with the wb97x-d3/cc-pvdz theory, demonstrating that the atom-category and element is not of importance whether a chemical shift is conformational sensitive or not. The most important feature that has an influence is the Δϕ angle, the per-atom alignability of the five replicas per conformation and the change in distance to an oxygen atom. The finding of increased importance of the Δϕ dihedral angle was also confirmed by feature importance calculations using the empirical methods (Fig. S6.1–S6.3, ESI†).

Fig. 9 Permutation importance of geometrical and biophysical features in regard to the conformational sensitivity of chemical shifts show that the Δϕ-angle is of importance to predict whether a chemical shift is sensitive to conformational change or not. An error bar represents ± one standard deviation.

3.2 Agreement between simulation and experimental results

Not only is the sensitivity in regard to conformational change relevant, but so is the absolute agreement with experimentally obtained values. Due to the liquid nature of the NMR sample and the measurement time of the NMR methods used, it must be assumed that the experimental chemical shifts are averages of all accessible conformations. Therefore, the experimental chemical shifts have to be compared with simulated observables calculated from a reasonably complete conformational ensemble. A set of two conformational extreme cases, as discussed here, does not constitute a complete conformational ensemble so it is not expected that the average matches the experiment exactly. In order to properly validate the predicted chemical shifts against the experimental values, much more extensive simulations would be required, and the chemical shifts would need to be computed to all relevant conformations. Still, systematic over- or under-prediction of chemical shifts can be assessed using an incomplete ensemble and should be avoided.

Fig. 10 shows a comparison between experiment and simulation for both a DFT-based QM/MM calculation (wb97x-d3/cc-pdvz) in vacuum and an empirical (PPM) prediction method. To compare the accuracy of the methods, the mean relative error of the simulation compared to the experiment was calculated using eqn (3) and is shown in Fig. 11.

Fig. 10 A comparison of the experimental chemical shifts (x-axis) and simulated chemical shifts (y-axis) allows an overview of the prediction. The diagonal gray line represents a perfect agreement between experiment and simulation, the blue dots indicate chemical shifts of the stretched conformation and the red dots indicate chemical shifts of the globular one. (a) The results obtained with DFT-calculations using the QM/MM method in vacuum using the wb97x-d3/cc-pdvz theory. Aliphatic carbons and hydrogen show a good quality of prediction, whereas the chemical shifts of amid protons is too low. (b) Chemical shifts calculated with the PPM software show very good agreement with the experiment both for the stretched as well as for the globular case.

Fig. 11 Overview of the mean relative error of the different approaches to predict chemical shifts. DFT-based methods achieved the best accuracies using the wb97x-d3 functional and cc-pdvz basis-set.

Even though the empirical methods show the strongest conformational sensitivity, the accuracy is also remarkably good. While the conformational sensitivity showed little dependence on the choice of functional and basis set, the influence on accuracy is slightly greater. Using DFT-based methods, the best accuracies were achieved using the wb97x-d3 functional and cc-pdvz basis-set independent of the fragmentation and solvation method.

3.3 Secondary chemical shifts

For each non-terminal residue, a secondary structure identifier ΔΔδαβ was calculated using secondary chemical shift data. In Fig. 12, the calculated secondary structure identifier values for the experimental dataset (gray bars) and both the globular (red dots) and stretched (blue dots) conformations can be seen. The experimentally obtained data show no indication that either alpha- or beta-structured conformers make up a significant share of the ensemble, a finding supported by the DSSP analysis of the molecular dynamics trajectory (Fig. S3.1, ESI†). Applying the same secondary chemical shifts analysis to both the stretched and globular conformers with data obtained from the UCBShiftX method, stronger derivations from the random coil can be observed. The globular conformer shows mostly slightly higher secondary structure identifier values but the majority of data points for both conformers remain in the region attributed to the random coil. It should be mentioned that secondary chemical shift based analysis is very sensitive to systematic offsets of chemical shift prediction and measurement, as it is a comparison with tabulated random coil chemical shifts. Thus, it may be easy to wrongly declare parts of the peptide to be either alpha- or beta-structured. Therefore, only simulated data from the UCBShiftX method is discussed here, as the smallest mean relative error to the experiment was observed with this method.

Fig. 12 The calculation of a secondary structure identifier ΔΔδαβ gives indications about the occurrence of secondary structures. To differentiate between random coil and non-random coil secondary structures, literature advises a secondary chemical shift of |Δδ| > 0.7 which is more or less continuous for more than four residues.55 As ΔΔδαβ is the difference between two secondary chemical shifts, which are expected to be additive in behavior, a threshold of 1.0 has been set arbitrarily and marked as a horizontal line. Gray bars represent data obtained from the experiment, blue from predictions for the stretched conformer and red from the globular conformer.

4 Discussion

4.1 Conformational sensitivity

The results of the calculations show that the conformational sensitivity is mainly limited by the precision of the estimator. Even small differences between the five equally treated overall conformations (Fig. 1) yield different chemical shifts with the DFT-based methods. This results in wide probability distribution peaks and thus weak conformational sensitivities measured as multiples of pooled standard deviation. The mean probability density function width of all methods can be seen in the ESI,† Chapter 1.1 which shows a significant difference between DFT-based and empirical methods. The use of explicit solvent leads to an even stronger decrease in conformational sensitivity. With explicit solvent, all five replicas of the same overall conformation feature micro-solvation with water molecules at different positions and orientations, taken from the MD-trajectory. Wide probability distributions can be observed especially with side chain atoms and protons, thus weakening the distinction between the expectation values of the stretched and globular group.

Regarding DFT-based prediction methods, the tested QM/MM approach which models protein parts in greater distance as MM-region show slightly better conformational sensitivity compared to the QM approach but the difference is small. A slightly bigger improvement was made when replacing the implicit solvent with a vacuum simulation for both the QM and QM/MM cases. It has to be noted that the vacuum simulations showed weaker convergence behavior, with some functional/basis-set combinations (Becke97-2 and Becke97-D functionals with the pcSseg-1 basis-set) unable to yield converged shieldings tensors for some atoms.

Empirical predictors yield results that are much closer to the experimental average for both the stretched and globular conformation and the influence of conformation on the absolute value of the chemical shift is much smaller than with the DFT-based calculations. Nevertheless, the results are more capable of differentiating between the two conformations.

4.2 Influence of features

Visual interpretation of the influence of atom categories and backbone torsion on the conformational sensitivity could only partially be reproduced with random forest permutation feature selection. It was not possible to find a relationship between atom-category and conformational sensitivity but the chemical shift prediction models could find agreement that a change in an amino acid's ϕ angle has influence on the conformational sensitivity of chemical shifts. It has to be mentioned that all four empirical models have been parameterized using the ϕ angle or a derived property as input feature. Nevertheless, the influence of the Δϕ angle was also witnessed using ab initio methods in vacuum and with implicit solvent.

The R2-scores of the regression models were shown to be weak, and only slightly better than those of a constant model. Many evaluated biophysical features were shown to be unrelated to the predicted conformational sensitivity. Still, the results of the feature importance are consistent over the methods. It has to be noted, that the random forest had to be trained on a very small data-set while the feature space was wide. The weak regression scores also explain the difficulties with visual interpretation of the results and confirms the necessity to perform a similar study with a bigger set of proteins to conclude whether the findings can be generalized.

4.3 Agreement between experiment and simulation

When comparing simulated and experimental chemical shifts to assess the quality of prediction, it has to be taken into account that experimental chemical shifts are an averaged property of a molecular ensemble. Thus, simulated chemical shifts must be calculated using a reasonably complete conformational ensemble. As this is not the case when using only two conformational extreme cases, it is not unexpected that the predicted chemical shifts of single conformations can show deviations from the experimental means.

In the ideal case, chemical shifts calculated from single conformations can be found to be clustered around the experimental mean. In practice, most DFT-based calculations show either a systematic under- or over-prediction of values (example: Fig. S8.14, S8.22 and S8.34, ESI†). The error can not only be explained by faulty single-point referencing, as chemical shifts of aliphatic carbons are often predicted to be slightly too high while carbonyl carbon chemical shifts may be too low using the same referencing. As publications from Rablen et al. and the Tantillo group show,62,63 it may be necessary to reference chemical shifts calculated using the DFT-based method not just by one single reference point (intercept) but also to calculate a scaling factor. The choice of a fitting basis-set is key to minimizing these systematic offsets, so that predicted values are evenly clustered around the experimental means (Fig. 10a).

With DFT-based approaches using vacuum simulations or implicit solvent, the absolute chemical shift values of the Hamide atoms were constantly underpredicted, even when other shifts were systematically overpredicted. The addition of discrete water molecules as a micro solvent helped to improve these absolute values and the error due to underestimation could be considerably reduced.

Compared to the DFT-based methods, empirical predictors do an excellent job reproducing experimental chemical shifts. There are no obvious systematical errors and both chemical shifts of the stretched and globular conformation match the experimental average very well.

5 Conclusions

This study explored different methods to calculate chemical shifts of proteins and their sensitivity in regard to protein conformation. Judging by the results of the evaluated test system, the majority of the chemical shifts are expected not to be sensitive to changes in overall conformation. We find that many chemical shifts that are predicted for very similar configurations, which would generally be considered as being in the same conformation, actually differ in a similar amount as chemical shifts predicted for two really distinct conformations. It should be noted that the TAU-protein fragment evaluated in this case study is unlikely to feature a significant secondary structure, neither in the experimental data nor observed in the simulation or in the selected conformers. While it is very possible that the choice of this fragment makes conformational differentiation using chemical shifts more difficult, it has to be expected that such regions showing no major secondary structure propensities occur often in the context of ensemble reweighting.

There is likely no relationship between atom type, atom name and element with regard to the conformational sensitivity of the chemical shift. Up to ∼26% of the calculated chemical shifts (UCBShiftX software) show conformational sensitivity in the case of the tested peptide but it remains difficult to predict why exactly those are sensitive while others are not. Compared to established chemical shift-based structure elucidation methods targeting conformations with stable secondary structures, particular attention must be paid to only select data that offers information about the conformation with reasonable probability when working with IDPs.

When comparing empirical methods with DFT-based ones, the most obvious difference is the efficiency and compute-time needed to fulfill the task. Empirical methods remain orders of magnitude faster than DFT-based calculations. While most of the empirical chemical shift predictors were designed for and trained by globular proteins, they are still capable of including most conformational sensitivity in the predicted chemical shifts in this case study.

Taking efficiency, time spent and unmatched accuracy into account, empirical predictors will remain the method of choice for most researchers to calculate chemical shifts even if they can only be applied to a subset of atoms.

Author contributions

JS conceptualized and managed the design and methodology of the project, implemented the necessary scripts, ran the calculations and wrote the manuscript. CO acted as the supervisor, organized the funding of the project, and edited and reviewed the manuscript. All authors have read and agreed to the published version of the manuscript.

Data availability

Data supporting the findings of this study are available at: https://zenodo.org/doi/10.5281/zenodo.11086149. The archive contains the 2 × 5 evaluated structures of the tested TAU-fragment, scripts to set-up DFT-based and empirical calculations and data analysis for all tested methods.

Conflicts of interest

The authors declare that there are no conflicts of interest.

Supplementary Material

CP-026-D4CP02484B-s001

The authors thank Krishnendu Bera and Jozef Hritz from CEITEC MU, Masaryk University, Brno for sharing experimental NMR data. The authors thank Peter Poliak from BOKU (University of Natural Resources and Life Sciences, Vienna) for sharing general advice. Financial support by the Austrian Science Fund (FWF; grant number I-4588) and by the Austrian Federal Ministry for Digital and Economic Affairs, the National Foundation for Research, Technology and Development and the Christian Doppler Research Association is gratefully acknowledged. Preparation and analysis of the study has been performed using various Python scripts. Besides other already mentioned software, the open source packages BioTite,64 RDKit,65 SciPy,66 NumPy,67 Pandas,68,69 Pillow,70 Matplotlib71 and Seaborn72 have been used in this study.
==== Refs
Fischer E. Einfluss der Configuration auf die Wirkung der Enzyme Ber. Dtsch. Chem. Ges. 1894 27 2985 2993 10.1002/cber.18940270364
Anfinsen C. B. Principles that Govern the Folding of Protein Chains Science 1973 181 223 230 10.1126/science.181.4096.223 4124164
Uversky V. N. Kulkarni P. Intrinsically disordered proteins: Chronology of a discovery Biophys. Chem. 2021 279 106694 10.1016/j.bpc.2021.106694 34607199
Kulkarni P. Leite V. B. P. Roy S. Bhattacharyya S. Mohanty A. Achuthan S. Singh D. Appadurai R. Rangarajan G. Weninger K. Orban J. Srivastava A. Jolly M. K. Onuchic J. N. Uversky V. N. Salgia R. Intrinsically disordered proteins: Ensembles at the limits of Anfinsen's dogma Biophys. Rev. 2022 3 011306 10.1063/5.0080512 38505224
Goedert M. Klug A. Crowther R. A. Tau protein, the paired helical filament and Alzheimer's disease J. Alzheimer's Dis. 2006 9 195 207
Lei P. Ayton S. Finkelstein D. I. Adlard P. A. Masters C. L. Bush A. I. Tau protein: Relevance to Parkinson's disease Int. J. Biochem. Cell Biol. 2010 42 1775 1778 10.1016/j.biocel.2010.07.016 20678581
Zhang X. Gao F. Wang D. Li C. Fu Y. He W. Zhang J. Tau Pathology in Parkinson’s Disease Front. Neurol. 2018 9 809 10.3389/fneur.2018.00809 30333786
Ward J. Sodhi J. McGuffin L. Buxton B. Jones D. Prediction and Functional Analysis of Native Disorder in Proteins from the Three Kingdoms of Life J. Mol. Biol. 2004 337 635 645 10.1016/j.jmb.2004.02.002 15019783
Fisher C. K. Stultz C. M. Constructing ensembles for intrinsically disordered proteins Curr. Opin. Struct. Biol. 2011 21 426 431 10.1016/j.sbi.2011.04.001 21530234
Lindorff-Larsen K. Best R. B. DePristo M. A. Dobson C. M. Vendruscolo M. Simultaneous determination of protein structure and dynamics Nature 2005 433 128 132 10.1038/nature03199 15650731
Chan-Yao-Chong M. Durand D. Ha-Duong T. Molecular Dynamics Simulations Combined with Nuclear Magnetic Resonance and/or Small-Angle X-ray Scattering Data for Characterizing Intrinsically Disordered Protein Conformational Ensembles J. Chem. Inf. Model. 2019 59 1743 1758 10.1021/acs.jcim.8b00928 30840442
Jensen M. R. Zweckstetter M. Huang J.-R. Blackledge M. Exploring Free-Energy Landscapes of Intrinsically Disordered Proteins at Atomic Resolution Using NMR Spectroscopy Chem. Rev. 2014 114 6632 6660 10.1021/cr400688u 24725176
Gama Lima Costa R. Fushman D. Reweighting methods for elucidation of conformation ensembles of proteins Curr. Opin. Struct. Biol. 2022 77 102470 10.1016/j.sbi.2022.102470 36183447
Facelli J. C. Chemical shift tensors: Theory and application to molecular structural problems Prog. Nucl. Magn. Reson. Spectrosc. 2011 58 176 201 10.1016/j.pnmrs.2010.10.003 21397119
Markley J. L. Meadows D. H. Jardetzky O. Nuclear magnetic resonance studies of helix-coil transitions in polyamino acids J. Mol. Biol. 1967 27 25 40 10.1016/0022-2836(67)90349-X 6033611
Wishart D. Sykes B. Richards F. Relationship between nuclear magnetic resonance chemical shift and protein secondary structure J. Mol. Biol. 1991 222 311 333 10.1016/0022-2836(91)90214-Q 1960729
Spera S. Bax A. Empirical correlation between protein backbone conformation and C.alpha. and C.beta. 13C nuclear magnetic resonance chemical shifts J. Am. Chem. Soc. 1991 113 5490 5492 10.1021/ja00014a071
Robustelli P. Stafford K. A. Palmer A. G. I. I. I. Interpreting Protein Structural Dynamics from NMR Chemical Shifts J. Am. Chem. Soc. 2012 134 6365 6374 10.1021/ja300265w 22381384
Sumowski C. V. Hanni M. Schweizer S. Ochsenfeld C. Sensitivity of ab Initio vs Empirical Methods in Computing Structural Effects on NMR Chemical Shifts for the Example of Peptides J. Chem. Theory Comput. 2014 10 122 133 10.1021/ct400713t 26579896
, PMID: 26579896

Christensen A. S. Linnet T. E. Borg M. Boomsma W. Lindorff-Larsen K. Hamelryck T. Jensen J. H. Protein Structure Validation and Refinement Using Amide Proton Chemical Shifts Derived from Quantum Mechanics PLoS One 2014 8 1 10
van Gunsteren W. F. Allison J. R. Daura X. Dolenc J. Hansen N. Mark A. E. Oostenbrink C. Rusu V. H. Smith L. J. Deriving Structural Information from Experimentally Measured Data on Biomolecules Angew. Chem., Int. Ed. 2016 55 15990 16010 10.1002/anie.201601828 27862777
Grutsch S. Brüschweiler S. Tollinger M. NMR Methods to Study Dynamic Allostery PLoS Comput. Biol. 2016 12 1 20 10.1371/journal.pcbi.1004620 26964042
Li D.-W. Brüschweiler R. Certification of Molecular Dynamics Trajectories with NMR Chemical Shifts J. Phys. Chem. Lett. 2010 1 246 248 10.1021/jz9001345
Markwick P. R. L. Cervantes C. F. Abel B. L. Komives E. A. Blackledge M. McCammon J. A. Enhanced Conformational Space Sampling Improves the Prediction of Chemical Shifts in Proteins J. Am. Chem. Soc. 2010 132 1220 1221 10.1021/ja9093692 20063881
Lasorsa A. Malki I. Cantrelle F.-X. Merzougui H. Boll E. Lambert J.-C. Landrieu I. Structural Basis of Tau Interaction With BIN1 and Regulation by Tau Phosphorylation Front. Mol. Neurosci. 2018 11 421 10.3389/fnmol.2018.00421 30487734
Lasorsa A. Bera K. Malki I. Dupré E. Cantrelle F.-X. Merzougui H. Sinnaeve D. Hanoulle X. Hritz J. Landrieu I. Conformation and Affinity Modulations by Multiple Phosphorylation Occurring in the BIN1 SH3 Domain Binding Site of the Tau Protein Proline-Rich Region Biochemistry 2023 62 1631 1642 10.1021/acs.biochem.2c00717 37167199
Lindorff-Larsen K. Piana S. Palmo K. Maragakis P. Klepeis J. L. Dror R. O. Shaw D. E. Improved side-chain torsion potentials for the Amber ff99SB protein force field Proteins 2010 78 1950 1958 10.1002/prot.22711 20408171
Eastman P. Swails J. Chodera J. D. McGibbon R. T. Zhao Y. Beauchamp K. A. Wang L.-P. Simmonett A. C. Harrigan M. P. Stern C. D. Wiewiora R. P. Brooks B. R. Pande V. S. OpenMM 7: Rapid development of high performance algorithms for molecular dynamics PLoS Comput. Biol. 2017 13 e1005659 10.1371/journal.pcbi.1005659 28746339
Izadi S. Anandakrishnan R. Onufriev A. V. Building Water Models: A Different Approach J. Phys. Chem. Lett. 2014 5 3863 3871 10.1021/jz501780a 25400877
Moussa J. E. and Stewart J. J. P. , MOPAC, version v22.1.0, 2023
Aprà E. Bylaska E. J. de Jong W. A. Govind N. Kowalski K. Straatsma T. P. Valiev M. van Dam H. J. J. Alexeev Y. Anchell J. Anisimov V. Aquino F. W. Atta-Fynn R. Autschbach J. Bauman N. P. Becca J. C. Bernholdt D. E. Bhaskaran- Nair K. Bogatko S. Borowski P. Boschen J. Brabec J. Bruner A. Cauët E. Chen Y. Chuev G. N. Cramer C. J. Daily J. Deegan M. J. O. Dunning, Jr T. H. Dupuis M. Dyall K. G. Fann G. I. Fischer S. A. Fonari A. Früchtl H. Gagliardi L. Garza J. Gawande N. Ghosh S. Glaesemann K. Götz A. W. Hammond J. Helms V. Hermes E. D. Hirao K. Hirata S. Jacquelin M. Jensen L. Johnson B. G. Jónsson H. Kendall R. A. Klemm M. Kobayashi R. Konkov V. Krishnamoorthy S. Krishnan M. Lin Z. Lins R. D. Littlefield R. J. Logsdail A. J. Lopata K. Ma W. Marenich A. V. Martin del Campo J. Mejia-Rodriguez D. Moore J. E. Mullin J. M. Nakajima T. Nascimento D. R. Nichols J. A. Nichols P. J. Nieplocha J. Otero-de-la-Roza A. Palmer B. Panyala A. Pirojsirikul T. Peng B. Peverati R. Pittner J. Pollack L. Richard R. M. Sadayappan P. Schatz G. C. Shelton W. A. Silverstein D. W. Smith D. M. A. Soares T. A. Song D. Swart M. Taylor H. L. Thomas G. S. Tipparaju V. Truhlar D. G. Tsemekhman K. Van Voorhis T. Vázquez-Mayagoitia Á. Verma P. Villa O. Vishnu A. Vogiatzis K. D. Wang D. Weare J. H. Williamson M. J. Windus T. L. Woliski K. Wong A. T. Wu Q. Yang C. Yu Q. Zacharias M. Zhang Z. Zhao Y. Harrison R. J. NWChem: Past, present, and future J. Chem. Phys. 2020 152 184102 10.1063/5.0004997 32414274
Wolinski K. Hinton J. F. Pulay P. Efficient implementation of the gauge-independent atomic orbital method for NMR chemical shift calculations J. Am. Chem. Soc. 1990 112 8251 8260 10.1021/ja00179a005
Han B. Liu Y. Ginzinger S. W. Wishart D. S. SHIFTX2: significantly improved protein chemical shift prediction J. Biomol. NMR 2011 50 43 57 10.1007/s10858-011-9478-4 21448735
Shen Y. Bax A. SPARTA +: a modest improvement in empirical NMR chemical shift prediction by means of an artificial neural network J. Biomol. NMR 2010 48 13 22 10.1007/s10858-010-9433-9 20628786
Li J. Bennett K. C. Liu Y. Martin M. V. Head-Gordon T. Accurate prediction of chemical shifts for aqueous protein structure on Real World data Chem. Sci. 2020 11 3180 3191 10.1039/C9SC06561J 34122823
Li D.-W. Brüschweiler R. PPM: a side-chain and backbone chemical shift predictor for the assessment of protein conformational ensembles J. Biomol. NMR 2012 54 257 265 10.1007/s10858-012-9668-8 22972619
Michaud-Agrawal N. Denning E. J. Woolf T. B. Beckstein O. MDAnalysis: A toolkit for the analysis of molecular dynamics simulations J. Comput. Chem. 2011 32 2319 2327 10.1002/jcc.21787 21500218
Gowers R. J. , Linke M. , Barnoud J. , Reddy T. J. E. , Melo M. N. , Seyler S. L. , Domaski J. , Dotson D. L. , Buchoux S. , Kenney I. M. and Beckstein O. , in Proceedings of the 15th Python in Science Conference, ed. S. Benthall and S. Rostrup, 2016, pp. 98–105
York D. M. Karplus M. A Smooth Solvation Potential Based on the Conductor-Like Screening Model J. Phys. Chem. A 1999 103 11060 11079 10.1021/jp992097l
Ditchfield R. Hehre W. J. Pople J. A. SelfConsistent MolecularOrbital Methods. IX. An Extended GaussianType Basis for MolecularOrbital Studies of Organic Molecules J. Chem. Phys. 1971 54 724 728 10.1063/1.1674902
Hehre W. J. Ditchfield R. Pople J. A. SelfConsistent Molecular Orbital Methods. XII. Further Extensions of GaussianType Basis Sets for Use in Molecular Orbital Studies of Organic Molecules J. Chem. Phys. 1972 56 2257 2261 10.1063/1.1677527
Hariharan P. C. Pople J. A. The influence of polarization functions on molecular orbital hydrogenation energies Theor. Chim. Acta 1973 28 213 222 10.1007/BF00533485
Dunning J. Thom H. Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen J. Chem. Phys. 1989 90 1007 1023 10.1063/1.456153
Jensen F. Segmented Contracted Basis Sets Optimized for Nuclear Magnetic Shielding J. Chem. Theory Comput. 2015 11 132 138 10.1021/ct5009526 26574211
Becke A. D. Densityfunctional thermochemistry. III. The role of exact exchange J. Chem. Phys. 1993 98 5648 5652 10.1063/1.464913
Wilson P. J. Bradley T. J. Tozer D. J. Hybrid exchange-correlation functional determined from thermochemical data and ab initio potentials J. Chem. Phys. 2001 115 9233 9242 10.1063/1.1412605
Grimme S. Semiempirical GGA-type density functional constructed with a long-range dispersion correction J. Comput. Chem. 2006 27 1787 1799 10.1002/jcc.20495 16955487
Lin Y.-S. Li G.-D. Mao S.-P. Chai J.-D. Long-Range Corrected Hybrid Density Functionals with Improved Dispersion Corrections J. Chem. Theory Comput. 2013 9 263 272 10.1021/ct300715s 26589028
Grimme S. Antony J. Ehrlich S. Krieg H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu J. Chem. Phys. 2010 132 154104 10.1063/1.3382344 20423165
Pavlíková Pecechtlová J. Mládek A. Zapletal V. Hritz J. Quantum Chemical Calculations of NMR Chemical Shifts in Phosphorylated Intrinsically Disordered Proteins J. Chem. Theory Comput. 2019 15 5642 5658 10.1021/acs.jctc.8b00257 31487161
Jameson C. J. Jameson A. K. Oppusunggu D. Wille S. Burrell P. M. Mason J. 15N nuclear magnetic shielding scale from gas phase studies J. Chem. Phys. 1981 74 81 88 10.1063/1.440797
Cramer C. , Essentials of Computational Chemistry: Theories and Models, Wiley, 2005, p. 347
Kabsch W. Sander C. Dictionary of protein secondary structure: Pattern recognition of hydrogen-bonded and geometrical features Biopolymers 1983 22 2577 2637 10.1002/bip.360221211 6667333
Dalgarno D. C. Levine B. A. Williams R. J. Structural information from NMR secondary chemical shifts of peptide alpha C-H protons in proteins Biosci. Rep. 1983 3 443 452 10.1007/BF01121955 6882888
Wishart D. S. Sykes B. D. The 13C chemical-shift index: a simple method for the identification of protein secondary structure using 13C chemical-shift data J. Biomol. NMR 1994 4 171 180 10.1007/BF00175245 8019132
Luca S. Filippov D. V. van Boom J. H. Oschkinat H. de Groot H. J. M. Baldus M. Secondary chemical shifts in immobilized peptides and proteins: A qualitative basis for structure refinement under Magic Angle Spinning J. Biomol. NMR 2001 20 325 331 10.1023/A:1011278317489 11563556
Nielsen J. T. Mulder F. A. A. POTENCI: prediction of temperature, neighbor and pH-corrected chemical shifts for intrinsically disordered proteins J. Biomol. NMR 2018 70 141 165 10.1007/s10858-018-0166-5 29399725
Richardson Laboratory, Duke University, Top8000 rotamer data (reference_data), Github, 2023
Williams C. J. Headd J. J. Moriarty N. W. Prisant M. G. Videau L. L. Deis L. N. Verma V. Keedy D. A. Hintze B. J. Chen V. B. Jain S. Lewis S. M. Arendall III W. B. Snoeyink J. Adams P. D. Lovell S. C. Richardson J. S. Richardson D. C. MolProbity: More and better reference data for improved all-atom structure validation Protein Sci. 2018 27 293 315 10.1002/pro.3330 29067766
Mitternacht S. FreeSASA: An open source C library for solvent accessible surface area calculations F1000Research 2016 5 189 26973785
Pedregosa F. Varoquaux G. Gramfort A. Michel V. Thirion B. Grisel O. Blondel M. Prettenhofer P. Weiss R. Dubourg V. Vanderplas J. Passos A. Cournapeau D. Brucher M. Perrot M. Duchesnay E. Scikit-learn: Machine Learning in Python J. Mach. Learn. Res. 2011 12 2825 2830
Rablen P. R. Pearlman S. A. Finkbiner J. A Comparison of Density Functional Methods for the Estimation of Proton Chemical Shifts with Chemical Accuracy J. Phys. Chem. A 1999 103 7357 7363 10.1021/jp9916889
Lodewyk M. W. Siebert M. R. Tantillo D. J. Computational Prediction of 1H and 13C Chemical Shifts: A Useful Tool for Natural Product, Mechanistic, and Synthetic Organic Chemistry Chem. Rev. 2012 112 1839 1862 10.1021/cr200106v 22091891
Kunzmann P. Hamacher K. Biotite: a unifying open source computational biology framework in Python BMC Bioinf. 2018 19 346 10.1186/s12859-018-2367-z 30285630
Landrum G. , Tosco P. , Kelley B. , Ric , Cosgrove D. , Sriniker G. , Vianello R. , Schneider N. , Kawashima E. , Jones D. N. G. , Dalke A. , Cole B. , Swain M. , Turk S. , Savelyev A. , Vaucher A. , Wójcikowski M. , Take I. , Probst D. , Ujihara K. , Scalfani V. F. , Godin G. , Pahl A. , Berenger F. , Varjo J. L. and Walker R. , jasondbiggs and strets123, rdkit/rdkit: 2023_03_1 (Q1 2023) Release, version Release_2023_03_1, 2023
Virtanen P. Gommers R. Oliphant T. E. Haberland M. Reddy T. Cournapeau D. Burovski E. Peterson P. Weckesser W. Bright J. van der Walt S. J. Brett M. Wilson J. Millman K. J. Mayorov N. Nelson A. R. J. Jones E. Kern R. Larson E. Carey Polat C. J. Feng Y. Moore E. W. VanderPlas J. Laxalde D. Perktold J. Cimrman R. Henriksen I. Quintero E. A. Harris C. R. Archibald A. M. Ribeiro A. H. Pedregosa F. van Mulbregt P. Vijaykumar A. Bardelli A. P. Rothberg A. Hilboll A. Kloeckner A. Scopatz A. Lee A. Rokem A. Woods C. N. Fulton C. Masson C. Häggström C. Fitzgerald C. Nicholson D. A. Hagen D. R. Pasechnik D. V. Olivetti E. Martin E. Wieser E. Silva F. Lenders F. Wilhelm F. Young G. Price G. A. Ingold G.-L. Allen G. E. Lee G. R. Audren H. Probst I. Dietrich J. P. Silterra J. Webber J. T. Slavi J. Nothman J. Buchner J. Kulick J. Schönberger J. L. de Miranda Cardoso J. V. Reimer J. Harrington J. Rodríguez J. L. C. Nunez-Iglesias J. Kuczynski J. Tritz K. Thoma M. Newville M. Kümmerer M. Bolingbroke M. Tartre M. Pak M. Smith N. J. Nowaczyk N. Shebanov N. Pavlyk O. Brodtkorb P. A. Lee P. McGibbon R. T. Feldbauer R. Lewis S. Tygier S. Sievert S. Vigna S. Peterson S. More S. Pudlik T. Oshima T. Pingel T. J. Robitaille T. P. Spura T. Jones T. R. Cera T. Leslie T. Zito T. Krauss T. Upadhyay U. Halchenko Y. O. Vázquez-Baeza Y. Contributors S. l SciPy 1.0: fundamental algorithms for scientific computing in Python Nat. Methods 2020 17 261 272 10.1038/s41592-019-0686-2 32015543
Harris C. R. Millman K. J. van der Walt S. J. Gommers R. Virtanen P. Cournapeau D. Wieser E. Taylor J. Berg S. Smith N. J. Kern R. Picus M. Hoyer S. van Kerkwijk M. H. Brett M. Haldane A. del Río J. F. Wiebe M. Peterson P. Gérard-Marchant P. Sheppard K. Reddy T. Weckesser W. Abbasi H. Gohlke C. Oliphant T. E. Array programming with NumPy Nature 2020 585 357 362 10.1038/s41586-020-2649-2 32939066
McKinney W. , in Proceedings of the 9th Python in Science Conference, ed. S. van der Walt and J. Millman, 2010, pp. 56–61
T. pandas development team, pandas-dev/pandas: Pandas, version v2.0.0, 2023
Murray A. , van Kemenade H. , Wiredfool , (Alex) J. A. C. , Karpinsky A. , Baranovi O. , Gohlke C. , Dufresne J. , Dwesl , Schmidt D. , Kopachev K. , Houghton A. , Mani S. , Landey S. , Ware J. , Vashek , Piolie , Douglas J. , Stanislau T. , Caro D. , Martinez U. , Kossouho S. , Lahd R. , Lee A. , Brown E. W. , Tonnhofer O. and Bonfill M. , python-pillow/Pillow: 9.5.0, version 9.5.0, 2023
Team T. M. D. , Matplotlib: Visualization with Python, version v3.8.1, 2023
Waskom M. L. seaborn: statistical data visualization J. Open Source Software 2021 6 3021 10.21105/joss.03021
