Predicting15N amide chemical shifts in proteins. I. An additive model for the backbone contribution

Predicted and Experimental NMR Chemical Shifts at Variable Temperatures: The Effect of Protein Conformational Dynamics

NMR chemical shifts provide a sensitive probe of protein structure and dynamics but remain challenging to predict and interpret. We examine the effect of protein conformational distributions on 15N chemical shifts for dihydrofolate reductase (DHFR), comparing QM/MM predicted shifts with experimental shifts in solution as well as frozen distributions. Representative snapshots from MD trajectories exhibit variation in predicted 15N chemical shifts of up to 25 ppm. The average over the fluctuations is in significantly better agreement with room temperature solution experimental values than the prediction for any single optimal conformations. Meanwhile, solid-state NMR (SSNMR) measurements of frozen solutions at 105 K exhibit broad lines whose widths agree well with the widths of distributions of predicted shifts for samples from the trajectory. The backbone torsion angle ψi-1 varies over 60° on the picosecond time scale, compensated by φi. These fluctuations can explain much of the shift variation.

The influence of random-coil chemical shifts on the assessment of structural propensities in folded proteins and IDPs

In studying secondary structural propensities of proteins by nuclear magnetic resonance (NMR) spectroscopy, secondary chemical shifts (SCSs) serve as the primary atomic scale observables. For SCS calculation, the selection of an appropriate random coil chemical shift (RCCS) dataset is a crucial step, especially when investigating intrinsically disordered proteins (IDPs). The scientific literature is abundant in such datasets, however, the effect of choosing one over all the others in a concrete application has not yet been studied thoroughly and systematically. Hereby, we review the available RCCS prediction methods and to compare them, we conduct statistical inference by means of the nonparametric sum of ranking differences and comparison of ranks to random numbers (SRD-CRRN) method. We try to find the RCCS predictors best representing the general consensus regarding secondary structural propensities. The existence and the magnitude of resulting differences on secondary structure determination under varying sample conditions (temperature, pH) are demonstrated and discussed for globular proteins and especially IDPs.

Predicted and Experimental NMR Chemical Shifts at Variable Temperatures: The Effect of Protein Conformational Dynamics

NMR chemical shifts provide a sensitive probe of protein structure and dynamics. Prediction of shifts, and therefore interpretation of shifts, particularly for the frequently measured amidic 15 N sites, remains a tall challenge. We demonstrate that protein 15 N chemical shift prediction from QM/MM predictions can be improved if conformational variation is included via MD sampling, focusing on the antibiotic target, E. coli Dihydrofolate reductase (DHFR). Variations of up to 25 ppm in predicted 15 N chemical shifts are observed over the trajectory. For solution shifts the average of fluctuations on the low picosecond timescale results in a superior prediction to a single optimal conformation. For low temperature solid state measurements, the histogram of predicted shifts for locally minimized snapshots with specific solvent arrangements sampled from the trajectory explains the heterogeneous linewidths; in other words, the conformations and associated solvent are 'frozen out' at low temperatures and result in inhomogeneously broadened NMR peaks. We identified conformational degrees of freedom that contribute to chemical shift variation. Backbone torsion angles show high amplitude fluctuations during the trajectory on the low picosecond timescale. For a number of residues, including I60, ψ varies by up to 60º within a conformational basin during the MD simulations, despite the fact that I60 (and other sites studied) are in a secondary structure element and remain well folded during the trajectory. Fluctuations in ψ appear to be compensated by other degrees of freedom in the protein, including φ of the succeeding residue, resulting in "rocking" of the amide plane with changes in hydrogen bonding interactions. Good agreement for both room temperature and low temperature NMR spectra provides strong support for the specific approach to conformational averaging of computed chemical shifts.

A comparison of quantum chemical models for calculating NMR shielding parameters in peptides: mixed basis set and ONIOM methods combined with a complete basis set extrapolation

This article compares several quantum mechanical approaches to the computation of chemical shielding tensors in peptide fragments. First, we describe the effects of basis set quality up to the complete basis set (CBS) limit and level of theory (HF, MP2, and DFT) for four different atoms in trans N-methylacetamide. For both isotropic shielding and shielding anisotropy, the MP2 results in the CBS limit show the best agreement with experiment. The HF values show quite a different tendency to MP2, and even in the CBS limit they are far from experiment for not only the isotropic shielding of carbonyl carbon but also most shielding anisotropies. In most cases, the DFT values differ systematically from MP2, and small basis-set (double- or triple-zeta) results are often fortuitously in better agreement with the experiment than the CBS ones. Second, we compare the mixed basis set and ONIOM methods, combined with CBS extrapolation, for chemical shielding calculations at a DFT level using various model peptides. From the results, it is shown that the mixed basis set method provides better results than ONIOM, compared to CBS calculations using the nonpartitioned full systems. The information studied here will be useful in guiding the selection of proper quantum chemical models, which are in a tradeoff between accuracy and cost, for shielding studies of peptides and proteins.

How Does an Amide-N Chemical Shift Tensor Vary in Peptides?

This study addresses a void in the existing literature on the amide-(15)N chemical shift anisotropy (CSA) tensor of peptides: a systematic investigation of how the tensor varies in different peptides. Amide-(15)N CSA tensors for several dipeptides are obtained using quantum chemical calculations, as well as for a series of model Ala-X and X-Ala sequences in both α-helical and β-sheet conformations (where X is one of the naturally occurring amino acids). The calculated values show a significant variation in both isolated and extended peptide structures. Hydrogen bonding at both the carbonyl group and the N-H bond of the peptide plane is shown to affect the principal values of the tensor. Calculations on model peptides indicate that the amide-(15)N CSA tensor is dependent on atoms located within a distance of five bonds. Consequently, the tensor of a given peptide residue is unaffected by residues other than those adjacent to it, which implies that the amide-(15)N CSA tensor should be considered in the context of tripeptide sequences. This further suggests that the amide-(15)N CSA tensor of the second residue of a given tripeptide sequence may be extrapolated to the same sequence in any other polypeptide or protein, given the same backbone conformation and intermolecular environment. These conclusions will facilitate future NMR structural studies of proteins.

15N Chemical Shielding in Glycyl Tripeptides: Measurement by Solid-State NMR and Correlation with X-ray Structure

15N chemical shielding parameters are reported for central glycyl residues in crystallographically characterized tripeptides with alpha-helix, beta-strand, polyglycine II (3(1)-helix), and extended structures. Accurate values of the shielding components (2-5 ppm) are determined from MAS and stationary spectra of peptides containing [2-(13)C,(15)N]Gly. Two dipolar couplings, (1)H-(15)N and (13)C(alpha)-(15)N, are used to examine (15)N shielding tensor orientations in the molecular frame and the results indicate that the delta(11), delta(33) plane of the shielding tensor is not coincident with the peptide plane. The observed isotropic shifts, which vary over a range of 13 ppm, depend on hydrogen bonding (direct and indirect) and local conformation. Tensor spans, delta(span) = delta(11) - delta(33), and their deviations from axial symmetry, delta(dev) = delta(22) - delta(33), vary over a larger range and are grouped according to 2 degrees structure. Augmented by previously reported (13)C(alpha) shielding parameters, a prediction scheme for the 2 degrees structure of glycyl residues in proteins based on shielding parameters is proposed.

Probing multiple effects on 15N, 13C alpha, 13C beta, and 13C' chemical shifts in peptides using density functional theory

We have used density functional calculations on model peptides to study conformational effects on (15)N, (13)C alpha, (13)C beta, and (13)C' chemical shifts, associated with hydrogen bonding, backbone conformation, and side-chain orientation. The results show a significant dependence on the backbone torsion angles of the nearest three residues. Contributions to (15)N chemical shifts from hydrogen bonding (up to 8 ppm), backbone conformation (up to 13 ppm), side-chain orientation and neighborhood residue effects (up to 22 ppm) are significant, and a unified theory will be required to account for their behavior in proteins. In contrast to this, the dependence on sequence and hydrogen bonding is much less for (13)C alpha and (13)C beta chemical shifts (<0.5 ppm), and moderate for carbonyl carbon shifts (<2 ppm). The effects of side-chain orientation are mainly limited to the residue itself for both nitrogen and carbon, but the chi(1) effect is also significant for the nitrogen shift of the following residue and for the (13)C' shift of the preceding residue. The calculated results are used, in conjunction with an additive model of chemical shift contributions, to create an algorithm for prediction of (15)N and (13)C shifts in proteins from their structure; this includes a model to extrapolate results to regions of torsion angle space that have not been explicitly studied by density functional theory (DFT) calculations. Crystal structures of 20 proteins with measured shifts have been used to test the prediction scheme. Root mean square deviations between calculated and experimental shifts 2.71, 1.22, 1.31, and 1.28 ppm for N, C alpha, C beta, and C', respectively. This prediction algorithm should be helpful in NMR assignment, crystal and solution structure comparison, and structure refinement.