Determination of the residue-specific 15N CSA tensor principal components using multiple alignment media

The individual components of the backbone ¹⁵N CSA tensor, σ₁₁, σ₂₂, σ₃₃, and the orientation of σ₁₁ relative to the NH bond described by the angle β have been determined for uniformly labeled ¹⁵N, ¹³C ubiquitin from partial alignment in phospholipid bicelles, Pf1 phage, and poly(ethylene glycol) by measuring the residue-specific residual dipolar couplings and chemical shift deviations. No strong correlation between any of the CSA tensor components is observed with any single structural feature. However, the experimentally determined tensor components agree with the previously determined average CSA principal components [Cornilescu and Bax (2000) J. Am. Chem. Soc. 122, 10143–10154]. Significant deviations from the averages coincide with residues in β-strand or extended regions, while α-helical residue tensor components cluster close to the average values.Electronic Supplementary Material Supplementary material is available to authorised users in the online version of this article at .     

Refinement of protein structure against non-redundant carbonyl <Superscript>13</Superscript>C NMR relaxation

Carbonyl ¹³C′ relaxation is dominated by the contribution from the ¹³C′ chemical shift anisotropy (CSA). The relaxation rates provide useful and non-redundant structural information in addition to dynamic parameters. It is straightforward to acquire, and offers complimentary structural information to the ¹⁵N relaxation data. Furthermore, the non-axial nature of the ¹³C′ CSA tensor results in a T₁/T₂ value that depends on an additional angular variable even when the diffusion tensor of the protein molecule is axially symmetric. This dependence on an extra degree of freedom provides new geometrical information that is not available from the NH dipolar relaxation. A protocol that incorporates such structural restraints into NMR structure calculation was developed within the program Xplor-NIH. Its application was illustrated with the yeast Fis1 NMR structure. Refinement against the ¹³C′ T₁/T₂ improved the overall quality of the structure, as evaluated by cross-validation against the residual dipolar coupling as well as the ¹⁵N relaxation data. In addition, possible variations of the CSA tensor were addressed. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Multidimensional oriented solid-state NMR experiments enable the sequential assignment of uniformly <Superscript>15</Superscript>N labeled integral membrane proteins in magnetically aligned lipid bilayers

Oriented solid-state NMR is the most direct methodology to obtain the orientation of membrane proteins with respect to the lipid bilayer. The method consists of measuring ¹H-¹⁵N dipolar couplings (DC) and ¹⁵N anisotropic chemical shifts (CSA) for membrane proteins that are uniformly aligned with respect to the membrane bilayer. A significant advantage of this approach is that tilt and azimuthal (rotational) angles of the protein domains can be directly derived from analytical expression of DC and CSA values, or, alternatively, obtained by refining protein structures using these values as harmonic restraints in simulated annealing calculations. The Achilles’ heel of this approach is the lack of suitable experiments for sequential assignment of the amide resonances. In this Article, we present a new pulse sequence that integrates proton driven spin diffusion (PDSD) with sensitivity-enhanced PISEMA in a 3D experiment ([¹H,¹⁵N]-SE-PISEMA-PDSD). The incorporation of 2D ¹⁵N/¹⁵N spin diffusion experiments into this new 3D experiment leads to the complete and unambiguous assignment of the ¹⁵N resonances. The feasibility of this approach is demonstrated for the membrane protein sarcolipin reconstituted in magnetically aligned lipid bicelles. Taken with low electric field probe technology, this approach will propel the determination of sequential assignment as well as structure and topology of larger integral membrane proteins in aligned lipid bilayers.     

Measurement of 15N relaxation in deuterated amide groups in proteins using direct nitrogen detection

15N chemical shielding tensors contain useful structural information, and their knowledge is essential for accurate analysis of protein backbone dynamics. The anisotropic component (CSA) of 15N chemical shielding can be obtained from 15N relaxation measurements in solution. However, the predominant contribution to nitrogen relaxation from 15N-(1)H dipolar coupling in amide groups limits the sensitivity of these measurements to the actual CSA values. Here we present nitrogen-detected NMR experiments for measuring 15N relaxation in deuterated amide groups in proteins, where the dipolar contribution to 15N relaxation is significantly reduced by the deuteration. Under these conditions nitrogen spin relaxation becomes a sensitive probe for variations in 15N chemical shielding tensors. Using the nitrogen direct-detection experiments we measured the rates of longitudinal and transverse 15N relaxation for backbone amides in protein G in D(2)O at 11.7 T. The measured relaxation rates are validated by comparing the overall rotational diffusion tensor obtained from these data with that from the conventional 15N relaxation measurements in H(2)O. This analysis revealed a 17-24 degree angle between the NH-bond and the unique axis of the 15N chemical shielding tensor.     

Vibrational averaging of chemical shift anisotropies in model peptides

The effects of chemical shift anisotropy (CSA) are evident in line-shapes or side-band analysis in solid-state NMR, in the observed line positions in partially oriented samples, and in relaxation effects in liquid-state studies. In all of these cases, the effective shielding tensor is influenced by fast vibrational averaging in addition to larger-amplitude internal motions and to overall libration or rotation. Here we compute the contributions of vibrational averaging (including zero-point motions) to the CSA relaxation strengths for the nitrogen and carbonyl carbon in two simple peptide models, and for snapshots taken from a path-integral simulation of a small protein. Because the ¹⁵N shielding tensor is determined by all the atoms of the peptide group, it is less influenced by vibrational motion than (for example) the N–H dipolar interaction, which is more sensitive to the motion of the light hydrogen atom. Computed order parameters for CSA averaging are hence much closer to unity than are N–H dipolar order parameters. This leads to a reduction by about 9% in the magnitude of the amide nitrogen CSA that is needed to fit liquid-state relaxation data. Similar considerations apply to the carbonyl carbon shielding tensor, but in this case the differences between dipolar and CSA averaging are smaller. These considerations will be important for making comparisons between CSA tensors extracted from various NMR experiments, and for comparisons to quantum chemical calculations carried out on static conformers.     

Analysis of the amide <Superscript>15</Superscript>N chemical shift tensor of the C<Subscript>α</Subscript> tetrasubstituted constituent of membrane-active peptaibols,the α-aminoisobutyric acid residue,compared to those of di- and tri-substituted proteinogenic amino acid residues

In protein NMR spectroscopy the chemical shift provides important information for the assignment of residues and a first structural evaluation of dihedral angles. Furthermore, angular restraints are obtained from oriented samples by solution and solid-state NMR spectroscopic approaches. Whereas the anisotropy of chemical shifts, quadrupolar couplings and dipolar interactions have been used to determine the structure, dynamics and topology of oriented membrane polypeptides using solid-state NMR spectroscopy similar concepts have been introduced to solution NMR through the measurements of residual dipolar couplings. The analysis of ¹⁵N chemical shift spectra depends on the accuracy of the chemical shift tensors. When investigating alamethicin and other peptaibols, i.e. polypeptides rich in α-aminoisobutyric acid (Aib), the ¹⁵N chemical shift tensor of this C_α-tetrasubstituted amino acid exhibits pronounced differences when compared to glycine, alanine and other proteinogenic residues. Here we present an experimental investigation on the ¹⁵N amide Aib tensor of N-acetyl-Aib-OH and for the Aib residues within peptaibols. Furthermore, a statistical analysis of the tensors published for di- (glycine) and tri-substituted residues has been performed, where for the first time the published data sets are compiled using a common reference. The size of the isotropic chemical shift and main tensor elements follows the order di- < tri- < tetra-substituted amino acids. A ¹⁵N chemical shift-¹H-¹⁵N dipolar coupling correlation NMR spectrum of alamethicin is used to evaluate the consequences of variations in the main tensor elements for the structural analysis of this membrane peptide.     

Quantitative estimation of magnitude and orientation of the CSA tensor from field dependence of longitudinal NMR relaxation rates

A method is presented that makes it possible to estimate both the orientation and the magnitude of the chemical shift anisotropy (CSA) tensor in molecules with a pair of spin 1/2 nuclei, typically ¹³C-¹H or ¹⁵ N-¹H. The method relies on the fact that the longitudinal cross-correlation rate as well as a linear combination of the autorelaxation rates of longitudinal heterospin magnetization, longitudinal two-spin order and longitudinal proton magnetization are proportional to the spectral density at the Larmor frequency of the heterospin. Therefore the ratio between the cross-correlation rate and the above linear combination is independent of the dynamics. From the field dependence of the ratio both the magnitude and the orientation of the CSA tensor can be estimated. The method is applicable to molecules in all motional regimes and is not limited to molecules in extreme narrowing or slow tumbling, nor is it sensitive to chemical exchange broadening. It is tested on the 22 amino acid residue peptide motilin, selectively ¹³ C labeled in the ortho positions in the ring of the single tyrosine residue. In the approximation of an axially symmetric ¹³C CSA tensor, the symmetry axis of the CSA tensor makes an angle of 23° ± 1° to the ¹³ C-¹H bond vector, and has a magnitude of 156 ± 5 ppm. This is in close agreement with solid-state NMR data on tyrosine powder [Frydman et al. (1992) Isr. J. Chem., 32, 161–164].     

Density functional calculations of 15N chemical shifts in solvated dipeptides

We performed density functional calculations to examine the effects of solvation, hydrogen bonding, backbone conformation, and the side chain on 15N chemical shielding in proteins. We used N-methylacetamide (NMA) and N-formyl-alanyl-X (with X being one of the 19 naturally occurring amino acids excluding proline) as model systems. In addition, calculations were performed for selected fragments from protein GB3. The conducting polarizable continuum model was employed to include the effect of solvent in the density functional calculations. Our calculations for NMA show that the augmentation of the polarizable continuum model with the explicit water molecules in the first solvation shell has a significant influence on isotropic 15N chemical shift but not as much on the chemical shift anisotropy. The difference in the isotropic chemical shift between the standard beta-sheet and alpha-helical conformations ranges from 0.8 to 6.2 ppm depending on the residue type, with the mean of 2.7 ppm. This is in good agreement with the experimental chemical shifts averaged over a database of 36 proteins containing >6100 amino acid residues. The orientation of the 15N chemical shielding tensor as well as its anisotropy and asymmetry are also in the range of values experimentally observed for peptides and proteins.     

Density functional calculations of backbone 15N shielding tensors in beta-sheet and turn residues of protein G

We performed density functional calculations of backbone ¹⁵N shielding tensors in the regions of beta-sheet and turns of protein G. The calculations were carried out for all twenty-four beta-sheet residues and eight beta-turn residues in the protein GB3 and the results were compared with the available experimental data from solid-state and solution NMR measurements. Together with the alpha-helix data, our calculations cover 39 out of the 55 residues (or 71%) in GB3. The applicability of several computational models developed previously (Cai et al. in J Biomol NMR 45:245–253, 2009) to compute ¹⁵N shielding tensors of alpha-helical residues is assessed. We show that the proposed quantum chemical computational model is capable of predicting isotropic ¹⁵N chemical shifts for an entire protein that are in good correlation with experimental data. However, the individual components of the predicted ¹⁵N shielding tensor agree with experiment less well: the computed values show much larger spread than the experimental data, and there is a profound difference in the behavior of the tensor components for alpha-helix/turns and beta-sheet residues. We discuss possible reasons for this.     

Estimates of methyl 13C and 1H CSA values (Δσ) in proteins from cross-correlated spin relaxation

Simple pulse schemes are presented for the measurement of methyl ¹³C and ¹H CSA values from ¹H–¹³C dipole/¹³C CSA and ¹H–¹³C dipole/¹H CSA cross-correlated relaxation. The methodology is applied to protein L and malate synthase G. Average ¹³C CSA values are considerably smaller for Ile than Leu/Val (17 vs 25 ppm) and are in good agreement with previous solid state NMR studies of powders of amino acids and dipeptides and in reasonable agreement with quantum-chemical DFT calculations of methyl carbon CSA values in peptide fragments. Small averaged ¹H CSA values on the order of 1 ppm are measured, consistent with a solid state NMR determination of the methyl group ¹H CSA in dimethylmalonic acid.     

Calculation of chemical shift anisotropy in proteins

Individual peptide groups in proteins must exhibit some variation in the chemical shift anisotropy (CSA) of their constituent atoms, but not much is known about the extent or origins of this dispersion. Direct spectroscopic measurement of CSA remains technically challenging, and theoretical methods can help to overcome these limitations by estimating shielding tensors for arbitrary structures. Here we use an automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) approach to compute ¹⁵N, ¹³C′ and ¹H chemical shift tensors for human ubiquitin and the GB1 and GB3 fragments of staphylococcal protein G. The average and range of variation of the anisotropies is in good agreement with experimental estimates from solid-state NMR, and the variation among residues is somewhat smaller than that estimated from solution-state measurements. Hydrogen-bond effects account for much of the variation, both between helix and sheet regions, and within elements of secondary structure, but other effects (including variations in torsion angles) may play a role as well.     

Chemical shifts in denatured proteins: Resonance assignments for denatured ubiquitin and comparisons with other denatured proteins

Chemical shift assignment is reported for the protein ubiquitin denatured in 8M urea at pH 2. The variations in ¹⁵N chemical shifts of three different proteins (ubiquitin, disulfide reduced, carboxymethylated lysozyme, all-Ala--lactalbumin), all without disulfides and denatured in 8M urea at pH 2 are compared to `random coil shifts' of small model peptides (Braun et al., 1994) and to the averaged native chemical shifts taken from the BMRB database. Both parameterizations show a remarkable agreement with the averaged measured ¹⁵N chemical shifts in the three denatured proteins. Detailed analysis of these experimental ¹⁵N chemical shifts provides an estimate of the influence of nearest neighbors and conformational preferences on the chemical shift and provides a direct means to identify non-random structural preferences in denatured proteins.     

Side-chain dynamics of the SAP SH2 domain correlate with a binding hot spot and a region with conformational plasticity

X-linked lymphoproliferative disease is caused by mutations in the protein SAP, which consists almost entirely of a single SH2 domain. SAP interacts with the Tyr281 site of the T<-->B cell signaling protein SLAM via its SH2 domain. Interestingly, binding is not dependent on phosphorylation but does involve interactions with residues N-terminal to the Tyr. We have used 15N and 2H NMR relaxation experiments to investigate the motional properties of the SAP SH2 domain backbone amides and side-chain methyl groups in the free protein and complexes with phosphorylated and non-phosphorylated peptides derived from the Tyr281 site of SLAM. The most mobile methyl groups are in side-chains with large RMSD values between the three crystal structures of SAP, suggesting that fast time-scale dynamics in side-chains is associated with conformational plasticity. The backbone amides of two residues which interact with the C-terminal part of the peptides experience fast time-scale motions in the free SH2 domain that are quenched upon binding of either the phosphorylated or non-phosphorylated peptide. Of most importance, the mobility of methyl groups in and around the binding site for residues in the N-terminus of the peptide is significantly restricted in the complexes, underscoring the dominance of this interaction with SAP and demonstrating a correlation between changes in rapid side-chain motion upon binding with local binding energy.     

NMR Resonance Assignment Methodology: Characterizing Large Sparsely Labeled Glycoproteins

Characterization of proteins using NMR methods begins with assignment of resonances to specific residues. This is usually accomplished using sequential connectivities between nuclear pairs in proteins uniformly labeled with NMR active isotopes. This becomes impractical for larger proteins, and especially for proteins that are best expressed in mammalian cells, including glycoproteins. Here an alternate protocol for the assignment of NMR resonances of sparsely labeled proteins, namely, the ones labeled with a single amino acid type, or a limited subset of types, isotopically enriched with ¹⁵N or ¹³C, is described. The protocol is based on comparison of data collected using extensions of simple two-dimensional NMR experiments (correlated chemical shifts, nuclear Overhauser effects, residual dipolar couplings) to predictions from molecular dynamics trajectories that begin with known protein structures. Optimal pairing of predicted and experimental values is facilitated by a software package that employs a genetic algorithm, ASSIGN_SLP_MD. The approach is applied to the 36-kDa luminal domain of the sialyltransferase, rST6Gal1, in which all phenylalanines are labeled with ¹⁵N, and the results are validated by elimination of resonances via single-point mutations of selected phenylalanines to tyrosines. Assignment allows the use of previously published paramagnetic relaxation enhancements to evaluate placement of a substrate analog in the active site of this protein. The protocol will open the way to structural characterization of the many glycosylated and other proteins that are best expressed in mammalian cells.     

Unraveling a phosphorylation event in a folded protein by NMR spectroscopy: phosphorylation of the Pin1 WW domain by PKA

The Pin1 protein plays a critical role in the functional regulation of the hyperphosphorylated neuronal Tau protein in Alzheimer’s disease and is by itself regulated by phosphorylation. We have used Nuclear Magnetic Resonance (NMR) spectroscopy to both identify the PKA phosphorylation site in the Pin1 WW domain and investigate the functional consequences of this phosphorylation. Detection and identification of phosphorylation on serine/threonine residues in a globular protein, while mostly occurring in solvent-exposed flexible loops, does not lead to chemical shift changes as obvious as in disordered proteins and hence does not necessarily shift the resonances outside the spectrum of the folded protein. Other complications were encountered to characterize the extent of the phosphorylation, as part of the ¹H,¹⁵N amide resonances around the phosphorylation site are specifically broadened in the unphosphorylated state. Despite these obstacles, NMR spectroscopy was an efficient tool to confirm phosphorylation on S16 of the WW domain and to quantify the level of phosphorylation. Based on this analytical characterization, we show that WW phosphorylation on S16 abolishes its binding capacity to a phosphorylated Tau peptide. A reduced conformational heterogeneity and flexibility of the phospho-binding loop upon S16 phosphorylation could account for part of the decreased affinity for its phosphorylated partner. Additionally, a structural model of the phospho-WW obtained by molecular dynamics simulation and energy minimization suggests that the phosphate moiety of phospho-S16 could compete with the phospho-substrate.     

Rapid automated determination of chemical shift anisotropy values in the carbonyl and carboxyl groups of fd-y21m bacteriophage using solid state NMR

Determination of chemical shift anisotropy (CSA) in immobilized proteins and protein assemblies is one of several tools to determine protein dynamics on the timescales of microseconds and faster. The large CSA values of C=O groups in the rigid limit makes them in particular attractive for measurements of large amplitude motions, or their absence. In this study, we implement a 3D R-symmetry-based sequence that recouples the second spatial component of the ¹³C CSA with the corresponding isotropic ¹³C′–¹³C cross-peaks in order to probe backbone and sidechain dynamics in an intact fd-y21m filamentous phage viral capsid. The assignment of the isotropic cross-peaks and the analysis were conducted automatically using a new software named ‘Raven’. The software can be utilized to auto-assign any 2D ¹³C–¹³C or ¹⁵N–¹³C spectrum given a previously-determined assignment table and generates simultaneously all intensity curves acquired in the third dimension. Here, all CSA spectra were automatically generated, and subsequently matched against a simulated set of CSA curves to yield their values. For the multi-copy, 50-residue-long protein capsid of fd-y21m, the backbone of the helical region is rigid, with reduced CSA values of ~?12.5 kHz (~?83 ppm). The N-terminus shows motionally-averaged CSA lineshapes and the carboxylic sidechain groups of four residues indicate large amplitude motions for D4, D5, D12 and E20. The current results further strengthen our previous studies of ¹⁵N CSA values and are in agreement with qualitative analysis of ¹³C–¹³C dipolar build-up curves, which were automatically obtained using our software. Our automated analysis technique is general and can be applied to study protein structure and dynamics, with data resulting from experiments that probe different variables such as relaxation rates and scaled anisotropic interactions.

Graphical abstract

Open image in new window

     

A novel experiment for the quantitative measurement of CSA(1HN)/CSA(15N) cross-correlated relaxation in 15N-labeled proteins*

An experiment is presented which allows for the quantitative measurement of the relaxation interference between the ¹H_N CSA and ¹⁵N CSA interactions in ¹⁵N labeled proteins. A constant-time buildup scheme is used to measure the differential relaxation rate, , between double-quantum (DQ) and zero-quantum (ZQ) ¹H_N-¹⁵N coherences. The CSA/CSA experiment was recorded at three different B_o field strengths. The CSA(¹H_N)/CSA(¹⁵N) cross-correlation rate was obtained from the linear fit of the measured rate, , versus B_o ² for 77 residues of the EH2 domain from mouse Eps15.     

15N relaxation study of the cold shock protein CspB at various solvent viscosities

For a detailed NMR study of the dynamics of the cold shock protein CspB from Bacillus subtilis, we determined ¹⁵N transverse and longitudinal relaxation rates and heteronuclear nuclear Overhauser effects at different solvent viscosities. Up to a relative viscosity of 2, which is equivalent to 27% ethylene glycol (EG), the overall correlation time follows the linear Stokes-Einstein equation. At a relative viscosity of 6 (70% EG) the correlation time deviates from linearity by 30%, indicating that CspB tumbles at a higher rate as expected from the solvent viscosity probably due to a preferential binding of water molecules at the protein surface. The corresponding hydrodynamic radii, determined by NMR diffusion experiments, show no variation with viscosity. The amplitudes of intramolecular motions on a sub-nanosecond time scale revealed by an extended Lipari–Szabo analysis were mainly independent of the solvent viscosity. The lower limit of the NMR `observation window' for the internal correlation time shifts above 0.5 ns at 70% EG, which is directly reflected in the experimentally derived internal correlation times. Chemical exchange contributions to the transverse relaxation rates derived from the Lipari-Szabo approach coincide with the experimentally determined values from the transverse ¹H-¹⁵N dipolar/¹⁵N chemical shift anisotropy relaxation interference. These contributions originate from fast protein folding reactions on a millisecond timescale, which get retarded at increased solvent viscosities.     

Pressure effect on the dynamics of an isolated α-helix studied by 15N-1H NMR relaxation

Dynamics and structure of (1–36)bacteriorhodopsin solubilized in chloroform/methanol mixture (1:1) were investigated by ¹H-¹⁵N NMR spectroscopy under a hydrostatic pressure of 2000 bar. It was shown that the peptide retains its spatial structure at high pressure. ¹⁵N transverse and longitudinal relaxation times, ¹⁵N{¹H} nuclear Overhauser effects, chemical shifts and the translation diffusion rate of the peptide at 2000 bar were compared with the respective data at ambient pressure [Orekhov et al. (1999) J. Biomol. NMR, 14, 345–356]. The model free analysis of the relaxation data for the helical 9–31 fragment revealed that the high pressure decreases the overall rotation and translation diffusion, as well as apparent order parameters of fast picosecond internal motions (S² _f) but has no effect on internal nanosecond motions (S² _s and _s) of the peptide. The decrease of translation and overall rotation diffusion was attributed to the increase in solvent viscosity and the decrease of apparent order parameters S² _f to a compression of hydrogen bonds. It is suggested that this compression causes an elongation of H-N bonds and a decrease of absolute values of chemical shift anisotropy (CSA). In particular, the observed decrease of S² _f at 2000 bar can be explained by 0.001 nm increase of N-H bond lengths and 10 ppm decrease of ¹⁵N CSA values.     

Mammalian production of an isotopically enriched outer domain of the HIV-1 gp120 glycoprotein for NMR spectroscopy

NMR spectroscopic characterization of the structure or the dynamics of proteins generally requires the production of samples isotopically enriched in ¹⁵N, ¹³C, or ²H. The bacterial expression systems currently in use to obtain isotopic enrichment, however, cannot produce a number of eukaryotic proteins, especially those that require post-translational modifications such as N-linked glycosylation for proper folding or activity. Here, we report the use of an adenovirus vector-based mammalian expression system to produce isotopically enriched ¹⁵N or ¹⁵N/¹³C samples of an outer domain variant of the HIV-1 gp120 envelope glycoprotein with 15 sites of N-linked glycosylation. Yields for the ¹⁵N- and ¹⁵N/¹³C-labeled gp120s after affinity chromatography were 45 and 44 mg/l, respectively, with an average of over 80% isotope incorporation. Recognition of the labeled gp120 by cognate antibodies that recognize complex epitopes showed affinities comparable to the unlabeled protein. NMR spectra, including ¹H-¹⁵N and ¹H-¹³C HSQCs, ¹⁵N-edited NOESY-HSQC, and 3D HNCO, were of high quality, with signal-to-noise consistent with an efficient level of isotope incorporation, and with chemical shift dispersion indicative of a well-folded protein. The exceptional protein yields, good isotope incorporation, and ability to obtain well-folded post-translationally modified proteins make this mammalian system attractive for the production of isotopically enriched eukaryotic proteins for NMR spectroscopy.