Assignment of protein backbone resonances using connectivity, torsion angles and 13Cα chemical shifts

A program is presented which will return the most probable sequence location for a short connected set of residues in a protein given just (13)C(alpha) chemical shifts (delta((13)C(alpha))) and data restricting the phi and psi backbone angles. Data taken from both the BioMagResBank and the Protein Data Bank were used to create a probability density function (PDF) using a multivariate normal distribution in delta((13)C(alpha)), phi, and psi space for each amino acid residue. Extracting and combining probabilities for particular amino acid residues in a short proposed sequence yields a score indicative of the correctness of the proposed assignment. The program is illustrated using several proteins for which structure and (13)C(alpha) chemical shift data are available.     

Automated prediction of 15N, 13Cα, 13Cβ and 13C′ chemical shifts in proteins using a density functional database

A database of peptide chemical shifts, computed at the density functional level, has been used to develop an algorithm for prediction of ¹⁵N and ¹³C shifts in proteins from their structure; the method is incorporated into a program called SHIFTS (version 4.0). The database was built from the calculated chemical shift patterns of 1335 peptides whose backbone torsion angles are limited to areas of the Ramachandran map around helical and sheet configurations. For each tripeptide in these regions of regular secondary structure (which constitute about 40% of residues in globular proteins) SHIFTS also consults the database for information about sidechain torsion angle effects for the residue of interest and for the preceding residue, and estimates hydrogen bonding effects through an empirical formula that is also based on density functional calculations on peptides. The program optionally searches for alternate side-chain torsion angles that could significantly improve agreement between calculated and observed shifts. The application of the program on 20 proteins shows good consistency with experimental data, with correlation coefficients of 0.92, 0.98, 0.99 and 0.90 and r.m.s. deviations of 1.94, 0.97, 1.05, and 1.08 ppm for ¹⁵N, ¹³C, ¹³C and ¹³C, respectively. Reference shifts fit to protein data are in good agreement with `random-coil' values derived from experimental measurements on peptides. This prediction algorithm should be helpful in NMR assignment, crystal and solution structure comparison, and structure refinement.     

Differential deuterium isotope shifts and one-bond 1H−13C scalar couplings in the conformational analysis of protein glycine residues

Summary The one-bond deuterium isotope shift effect for glycine C resonances exhibits a conformational dependence comparable to that of the corresponding ¹J_HC scalar coupling in both magnitude (11 Hz at 14.1 T) and dihedral angle dependence. The similarity in the conformational dependence of the ¹J_HC and deuterium isotope shift values suggests a common physical basis. Given the known distribution of (,) main-chain dihedral angles for glycine residues, the deuterium isotope shifts and the ¹J_HC scalar couplings can determine conformations in the left-and right-handed helical-to-bridge regions of the (,) plane to an accuracy of approximately 13°. In the absence of stereochemical assignments, the differential deuterium isotope shifts and the ¹J_HC scalar couplings can be combined with limited independent structural information (e.g., the sign of ) to determine the chirality of the deuterium substitution.     

Measurement of one-bond 15N-13C′ dipolar couplings in medium sized proteins

A simple and accurate method is described for measurement of ¹ J _CN splittings in isotopically enriched proteins. The method is of the quantitative J correlation type, and the ¹ J _CN splitting is derived from the relative intensity in two 3D TROSY-HNCO spectra with ¹ J _CN dephasing intervals of 1/(2¹ J _CN) (reference intensity) and 1/¹ J _CN (residual intensity). If the two spectra are recorded under identical conditions and with the same number of scans, the random error in the ¹ J _CN value extracted in this manner is inversely related to the signal-to-noise (S/N) in the reference spectrum. A S/N of 30:1 in the reference spectrum yields random errors of less than 0.2 Hz in the extracted ¹ J _CN value. Dipolar couplings obtained from the difference in ¹ J _CN splitting in the isotropic and liquid crystalline phase for the C-terminal domain of calmodulin are in excellent agreement with its 1.68-Å crystal structure, but agree considerably less with the 2.2-Å structure.     

Pulse schemes for the measurement of3 JC′Cγ and3 JNCγ scalar couplings in 15N,13C uniformly labeled proteins

Pulse sequences are presented for the measurement of³J_CC and³J_NC scalar couplings for allC containing residues in¹⁵N,¹³C uniformly labeled proteins. The methodsdescribed are based on quantitative J correlation spectroscopy pioneered byBax and co-workers [Bax et al. (1994) Methods Enzymol., 239, 79–105].The combination of ³J_CC and³J_NC scalar coupling constants allows theassignment of discrete rotameric states about the ₁ torsion angle in cases where such states exist or, alternatively,facilitates the establishment of noncanonical ₁conformations or the presence of rotameric averaging. The methods areapplied to a 1.5 mM sample of staphylococcal nuclease.     

Characterization of the hydrogen bond network in guanosine quartets by internucleotide 3hJNC′ and 2hJNN scalar couplings

Scalar coupling correlations across hydrogen bonds with carbonyl groups as acceptors have been observed in a variety of proteins, but not in nucleic acids. Here we present a pulse scheme that allows such an observation and quantification of trans-hydrogen bond ^3hJ_NC correlations in nucleic acid base pairs, between the imino nitrogen ¹⁵N1 and the carbonyl ¹³C6 nuclei within the guanine quartets of the Oxy-1.5 DNA-quadruplex. Intra- and internucleotide N-H···O=C connectivities can be traced around each guanine quartet, allowing the hydrogen bonding partners to be unambiguously assigned. Absolute values of the ^3hJ_NC couplings are approximately 0.2 Hz as quantified by a selective long-range H(N)CO experiment and are thus on average smaller than the analogous ^3hJ_NC couplings observed in proteins. In addition, an improved version of the pseudo-heteronuclear H(N)N-COSY [Majumdar et al. (1999) J. Biomol. NMR, 14, 67–70] is presented which allows simultaneous detection of the ¹⁵N-donor and ¹⁵N-acceptor resonances connected by ^2hJ_NN couplings in hydrogen bonds involving amino groups. Using this experiment, values ranging between 6 and 8 Hz are determined for the ^2hJ_NN couplings between ¹⁵N2 and ¹⁵N7 nuclei in the guanine quartet. These values are not strongly influenced by the presence of a significant amount of chemical exchange broadening due to amino group rotations.     

Leucine side-chain rotamers in a glycophorin A transmembrane peptide as revealed by three-bond carbon—carbon couplings and 13C chemical shifts

Summary We have used a spin-echo difference NMR pulse sequence to measure three-bond J couplings between - and -carbons of the leucine residues in a micelle-associated helical peptide dimer that corresponds to residues 62–101 of the transmembrane erythrocyte protein glycophorin A. The observed ³J couplings correlate strongly with the ¹³C chemical shift of the -methyl groups, and within experimental error both the shift distribution of the methyl carbons and the variations in ³J can be accounted for by variations in side-chain rotamer populations. We infer that all leucine side chains in this peptide dimer are in fast exchange among X ² rotamers and sample two of the three possible rotameric states, even when the side chain forms part of the dimer interface. The observed correlation of chemical shift with couplings can be traced to a -gauche interaction of methyl and -carbons. This correlation may provide an alternate route to rotamer analysis in some protein systems.     

Deuterium isotope shifts for backbone 1H, 15N and 13C nuclei in intrinsically disordered protein α-synuclein

Intrinsically disordered proteins (IDPs) are abundant in nature and characterization of their potential structural propensities remains a widely pursued but challenging task. Analysis of NMR secondary chemical shifts plays an important role in such studies, but the output of such analyses depends on the accuracy of reference random coil chemical shifts. Although uniform perdeuteration of IDPs can dramatically increase spectral resolution, a feature particularly important for the poorly dispersed IDP spectra, the impact of deuterium isotope shifts on random coil values has not yet been fully characterized. Very precise ²H isotope shift measurements for ¹³C^??, ¹³C^??, ¹³C??, ¹⁵N, and ¹H^N have been obtained by using a mixed sample of protonated and uniformly perdeuterated ??-synuclein, a protein with chemical shifts exceptionally close to random coil values. Decomposition of these isotope shifts into one-bond, two-bond and three-bond effects as well as intra- and sequential residue contributions shows that such an analysis, which ignores conformational dependence, is meaningful but does not fully describe the total isotope shift to within the precision of the measurements. Random coil ²H isotope shifts provide an important starting point for analysis of such shifts in structural terms in folded proteins, where they are known to depend strongly on local geometry.     

Measurement of residual dipolar couplings from 1Hα to 13Cα and 15N using a simple HNCA-based experiment

Novel NMR pulse schemes for simultaneous measurement of ¹ D _CHand ² D _NHresidual dipolar couplings in proteins is presented. We show that ² D _NHcoupling can be very useful for protein structure determination. The ² D _NHcoupling can be measured from ¹⁵N dimension with good accuracy on a slowly relaxing TROSY resonance, utilizing HNCA-TROSY-based experiments, which concomitantly supply large ¹ D _CHcoupling. The dynamic range of ² D _NHcoupling is comparable to ¹ D _NC coupling, but instead, it also serves non-redundant information on the course of protein backbone, thanks to rotational degree of freedom with respect to peptide bond. The HNCA-TROSY-based experiments are optimal for measuring residual dipolar couplings at high magnetic fields owing to absence of rapid transverse relaxation of carbonyl carbon. The reliability of the proposed approach was tested on ¹⁵N/¹³C human ubiquitin. A very good correlation with ubiquitin solution as well as crystal structure, for both ¹ D _CHand ² D _NHcouplings, was obtained.     

TROSY-based HNCO pulse sequences for the measurement of 1HN-15N, 15N-13CO, 1HN-13CO, 13CO-13Cα and 1HN-13Cα dipolar couplings in 15N, 13C, 2H-labeled proteins

HNCO-based 3D pulse schemes are presented for measuring ¹HN-¹⁵N,¹⁵N-¹³CO, ¹HN-¹³CO,¹³CO-¹³C and ¹HN-¹³C dipolar couplings in ¹⁵N,¹³C,²-labeled proteins. The experiments are based on recently developed TROSY methodology for improving spectral resolution and sensitivity. Data sets recorded on a complex of Val, Leu, Ile (1 only) methyl protonated ¹⁵N,¹³C,²H-labeled maltose binding protein and -cyclodextrin as well as ¹⁵N,¹³C,²H-labeled human carbonic anhydrase II demonstrate that precise dipolar couplings can be obtained on proteins in the 30–40 kDa molecular weight range. These couplings will serve as powerful restraints for obtaining global folds of highly deuterated proteins.     

MQ-HCN-based pulse sequences for the measurement of 13C1′-1H1′, 13C1′-15N, 1H1′-15N, 13C1′-13C2′, 1H1′-13C2′, 13C6/8-1H6/8, 13C6/8-15N, 1H6/8-15N, 13C6-13C5, 1H6-13C5 dipolar couplings in 13C, 15N-labeled DNA (and RNA)

A suite of multiple quantum (MQ) HCN-based pulse sequences has been developed for the purpose of collecting dipolar coupling data in labeled nucleic acids. All the pulse sequences are based on the robust MQ-HCN experiment which has been utilized for assignment purposes in labeled nucleic acids for a number of years and provides much-needed resolution for the dipolar coupling measurements. We have attempted to collect multiple couplings centered on the 13C1' and 13C6/8 positions. Six pulse sequences are described, one each for measurement of one-bond 13C1'-1H1' and 13C6/8-1H6/8 couplings, one for measurement of one-bond 13C1'-15N and two-bond 1H1'-15N couplings, one for measurement of one-bond 13C6/8-15N and two-bond 1H6/8-15N couplings, one for measurement of one-bond 13C1'- 13C2' and two-bond 1H1'-13C2' couplings, and one for measurement of one-bond 13C6-13C5 and two-bond 1H6-13C5 couplings in the bases of C and T. These sequences are demonstrated for a labeled 18 bp DNA duplex in a 47 kDa ternary complex of DNA, CBFbeta, and the CBFalpha Runt domain, thus clearly demonstrating the robustness of the pulse sequences even for a very large complex.     

A 4D TROSY-based pulse scheme for correlating 1HNi,15Ni,13Cαi,13C′i−1 chemical shifts in high molecular weight, 15N,13C, 2H labeled proteins

A 4D TROSY-based triple resonance experiment, 4D-HNCO_i–1CA_i, is presented which correlates intra-residue ¹HN, ¹⁵N, ¹³ C chemical shifts with the carbonyl (¹³C) shift of the preceding residue. The experiment is best used in concert with recently described 4D TROSY-HNCOCA and -HNCACO experiments [Yang, D. and Kay, L.E. (1999) J. Am. Chem. Soc., 121, 2571–2575]. In cases where degeneracy of (¹HN,¹⁵N) spin pairs precludes assignment using the HNCOCA and HNCACO, the HNCO_i–1CA_i often allows resolution of the ambiguity by linking the ¹³C and ¹³C spins surrounding the (¹HN,¹⁵N) pair. The experiment is demonstrated on a sample of ¹⁵N, ¹³C, ² H labeled maltose binding protein in complex with -cyclodextrin that tumbles with a correlation time of 46 ns.     

1H, 13C and 15N resonance assignments of the bb′ domains of human protein disulfide isomerase

Protein disulfide isomerase (PDI) participates in protein folding and catalyses formation of disulfide bonds. The b′ domain of human PDI contributes to binding unfolded proteins; its structure is stabilized by the b domain. Here, we report NMR chemical shift assignments for the bb′ fragment.     

Heteronuclear relayed E.COSY applied to the determination of accurate 3J(HN,C′) and 3J(Hβ,C′) coupling constants in Desulfovibrio vulgaris flavodoxin

Summary A simple constant-time 3D heteronuclear NMR pulse sequence has been developed to quantitatively determine the heteronuclear three-bond couplings ³J(H^N,C) and ³J(H,C) in uniformly ¹³C-enriched proteins. The protocols for measuring accurate coupling constants are based on ¹H,¹³C-heteronuclear relayed E.COSY [Schmidt, J.M., Ernst, R.R., Aimoto, S. and Kainosho, M. (1995) J. Biomol. NMR, 6, 95–105] in combination with numerical least-squares spectrum evaluation. Accurate coupling constants are extracted from 2D spectrum projections using 2D multiplet simulation. Confidence intervals for the obtained three-bond coupling constants are calculated from F-statistics. The three-bond couplings are relevant to the determination of and X ₁ dihedral-angle conformations in the amino acid backbone and side chain. The methods are demonstrated on the recombinant ¹³C, ¹⁵N-doubly enriched 147-amino acid protein Desulfovibrio vulgaris flavodoxin with bound flavin mononucleotide in its oxidized form. In total, 109 ³J(H^N,C) and 100 ³J(H,C) coupling constants are obtained from a single spectrum.Abbreviations ANOVA analysis of variances - COSY correlated spectroscopy - E.COSY exclusive correlation spectroscopy - FMN flavin mononucleotide - HMQC heteronuclear multiple-quantum coherence - HSQC heteronuclear single-quantum coherence     

Backbone and Ile-δ1, Leu, Val Methyl 1H, 13C and 15N NMR chemical shift assignments for human interferon-stimulated gene 15 protein

Human interferon-stimulated gene 15 protein (ISG15), also called ubiquitin cross-reactive protein (UCRP), is the first identified ubiquitin-like protein containing two ubiquitin-like domains fused in tandem. The active form of ISG15 is conjugated to target proteins via the C-terminal glycine residue through an isopeptide bond in a manner similar to ubiquitin. The biological role of ISG15 is strongly associated with the modulation of cell immune function, and there is mounting evidence suggesting that many viral pathogens evade the host innate immune response by interfering with ISG15 conjugation to both host and viral proteins in a variety of ways. Here we report nearly complete backbone ¹H^N, ¹⁵N, ¹³C′, and ¹³C^α, as well as side chain ¹³C^β, methyl (Ile-δ1, Leu, Val), amide (Asn, Gln), and indole N–H (Trp) NMR resonance assignments for the 157-residue human ISG15 protein. These resonance assignments provide the basis for future structural and functional solution NMR studies of the biologically important human ISG15 protein.     

Intraresidue 1H-15N-13C′ and 1Hα-13Cα-13C′ dipole-CSA relaxation interference as a source of constraints for structural refinement of metal-binding sites in zinc-finger proteins

¹H(i)-¹⁵N(i)-¹³C(i) dipole-chemical shift anisotropy (CSA) relaxation interference was quantified for the ¹³C,¹⁵N labeled zinc-finger protein qCRP2(LIM2). The cross-correlation rates obtained for residues located in the metal coordination sites of qCRP2(LIM2) show a high degree of correlation with the peptide plane torsion angles and taken from the solution structure. ¹H(i)-¹⁵N(i)-¹³C(i) as well as ¹³C(i)-¹H(i)-¹³C(i) dipole-CSA cross-correlation rates were subsequently used to improve the geometry of the metal binding site. The optimized dihedral angles of the two zinc-binding sites in qCRP2(LIM2) are in better agreement with values obtained from crystal structures of other zinc-finger proteins and thus establish the utility of this approach to improve the metal-binding site geometry of zinc-finger proteins studied by NMR spectroscopy in solution.     

1H, 13C, and 15N backbone and side-chain chemical shift assignments for the 29 kDa human galectin-1 protein dimer

Galectin-1 is an important regulator of leukocyte function and tumor angiogenesis. Recently, this lectin has been identified as a molecular target for the potent angiogenesis inhibitor anginex. Here, we report ¹H, ¹³C, and ¹⁵N chemical shift assignments for human galectin-1 as determined by using heteronuclear triple resonance NMR spectroscopy.     

Accurate measurements of the effects of deuteration at backbone amide positions on the chemical shifts of 15N, 13Cα, 13Cβ, 13CO and 1Hα nuclei in proteins

An approach towards accurate NMR measurements of deuterium isotope effects on the chemical shifts of all backbone nuclei in proteins (¹⁵N, ¹³C_α, ¹³CO, ¹H_α) and ¹³C_β nuclei arising from ¹H-to-D substitutions at amide nitrogen positions is described. Isolation of molecular species with a defined protonation/deuteration pattern at successive backbone nitrogen positions in the polypeptide chain allows quantifying all deuterium isotope shifts of these nuclei from the first to the fourth order. Some of the deuterium isotope shifts measured in the proteins ubiquitin and GB1 can be interpreted in terms of backbone geometry via empirical relationships describing their dependence on (φ; ψ) backbone dihedral angles. Because of their relatively large variability and notable dependence on the protein secondary structure, the two- and three-bond ¹³C_α isotope shifts, ²ΔC_α(N_iD) and ³ΔC_α(N_i+1D), and three-bond ¹³C_β isotope shifts, ³ΔC_β(N_iD), are useful reporters of the local geometry of the protein backbone.     

Quantification of H/D Isotope Effects on Protein Hydrogen-bonds by h3JNC′ and 1JNC′ Couplings and Peptide Group 15N and 13C′ Chemical Shifts

The effect of hydrogen/deuterium exchange on protein hydrogen bond coupling constants (h3)J(NC') has been investigated in the small globular protein ubiquitin. The couplings across deuterated or protonated hydrogen bonds were measured by a long-range quantitative HA(CACO)NCO experiment. The analysis is combined with a determination of the H(N)/D(N) isotope effect on the amide group (1)J(NC') couplings and the (15)N and (13)C' chemical shifts. On average, H-bond deuteration exchange weakens (h3)J(NC') and strengthens (1)J(NC') couplings. A correlation is found between the size of the (15)N isotope shift, the (15)N chemical shift, and the (h3)J(NC') coupling constants. The data are consistent with a reduction of donor-acceptor overlap as expected from the classical Ubbelohde effect and the common understanding that H(N)/D(N) exchange leads to a shortening of the N-hydron bond length.