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.     

The effect of noncollinearity of 15N-1H dipolar and 15N CSA tensors and rotational anisotropy on 15N relaxation, CSA/dipolar cross correlation, and TROSY

Current approaches to 15N relaxation in proteins assume that the 15N-1H dipolar and 15N CSA tensors are collinear. We show theoretically that, when there is significant anisotropy of molecular rotation, different orientations of the two tensors, experimentally observed in proteins, nucleic acids, and small peptides, will result in differences in site- specific correlation functions and spectral densities. The standard treatments of the rates of longitudinal and transverse relaxation of amide 15N nuclei, of the 15N CSA/15N-1H dipolar cross correlation, and of the TROSY experiment are extended to account for the effect of noncollinearity of the 15N-1H dipolar and 15N CSA (chemical shift anisotropy) tensors. This effect, proportional to the degree of anisotropy of the overall motion, (D/D–1), is sensitive to the relative orientation of the two tensors and to the orientation of the peptide plane with respect to the diffusion coordinate frame. The effect is negligible at small degrees of anisotropy, but is predicted to become significant for D/D1.5, and at high magnetic fields. The effect of noncollinearity of 15N CSA and 15N-1H dipolar interaction is sensitive to both gross (hydrodynamic) properties and atomic-level details of protein structure. Incorporation of this effect into relaxation data analysis is likely to improve both precision and accuracy of the derived characteristics of protein dynamics, especially at high magnetic fields and for molecules with a high degree of anisotropy of the overall motion. The effect will also make TROSY efficiency dependent on local orientation in moderately anisotropic systems.     

Temperature dependence of fast carbonyl backbone dynamics in chicken villin headpiece subdomain

Temperature-dependence of protein dynamics can provide information on details of the free energy landscape by probing the characteristics of the potential responsible for the fluctuations. We have investigated the temperature-dependence of picosecond to nanosecond backbone dynamics at carbonyl carbon sites in chicken villin headpiece subdomain protein using a combination of three NMR relaxation rates: ¹³C′ longitudinal rate, and two cross-correlated rates involving dipolar and chemical shift anisotropy (CSA) relaxation mechanisms, ¹³C′/¹³C′-¹³C^α CSA/dipolar and ¹³C′/¹³C′–¹⁵N CSA/dipolar. Order parameters have been extracted using the Lipari-Szabo model-free approach assuming a separation of the time scales of internal and molecular motions in the 2–16°C temperature range. There is a gradual deviation from this assumption from lower to higher temperatures, such that above 16°C the separation of the time scales is inconsistent with the experimental data and, thus, the Lipari-Szabo formalism can not be applied. While there are variations among the residues, on the average the order parameters indicate a markedly steeper temperature dependence at backbone carbonyl carbons compared to that probed at amide nitrogens in an earlier study. This strongly advocates for probing sites other than amide nitrogen for accurate characterization of the potential and other thermodynamics characteristics of protein backbone.     

Influence of N-H...O and O-H...O hydrogen bonds on the (17)O, (15)N and (13)C chemical shielding tensors in crystalline acetaminophen: a density functional theory study

A computational investigation was carried out to characterize the (17)O, (15)N and (13)C chemical shielding tensors in crystalline acetaminophen. We found that N-H...O and O-H...O hydrogen bonds around the acetaminophen molecule in the crystal lattice have different influences on the calculated (17)O, (15)N and (13)C chemical shielding eigenvalues and their orientations in the molecular frame of axes. The calculations were performed with the B3LYP method and 6-311++G(d, p) and 6-311+G(d) standard basis sets using the Gaussian 98 suite of programs. Calculated chemical shielding tensors were used to evaluate the (17)O, (15)N, and (13)C NMR chemical shift tensors in crystalline acetaminophen, which are in reasonable agreement with available experimental data. The difference between the calculated NMR parameters of the monomer and molecular clusters shows how much hydrogen-bonding interactions affect the chemical shielding tensors of each nucleus. The computed (17)O chemical shielding tensor on O(1), which is involved in two intermolecular hydrogen bonds, shows remarkable sensitivity toward the choice of the cluster model, whereas the (17)O chemical shielding tensor on O(2) involved in one N-H...O hydrogen bond, shows smaller improvement toward the hydrogen-bonding interactions. Also, a reasonably good agreement between the experimentally obtained solid-state (15)N and (13)C NMR chemical shifts and B3LYP/6-311++G(d, p) calculations is achievable only in molecular cluster model where a complete hydrogen-bonding network is considered. Moreover, at the B3LYP/6-311++G(d, p) level of theory, the calculated (17)O, (15)N and (13)C chemical shielding tensor orientations are able to reproduce the experimental values to a reasonably good degree of accuracy.     

Influence of the O-phosphorylation of serine, threonine and tyrosine in proteins on the amidic 15N chemical shielding anisotropy tensors

Density functional theory was employed to study the influence of O-phosphorylation of serine, threonine, and tyrosine on the amidic ¹⁵N chemical shielding anisotropy (CSA) tensor in the context of the complex chemical environments of protein structures. Our results indicate that the amidic ¹⁵N CSA tensor has sensitive responses to the introduction of the phosphate group and the phosphorylation-promoted rearrangement of solvent molecules and hydrogen bonding networks in the vicinity of the phosphorylated site. Yet, the calculated ¹⁵N CSA tensors in phosphorylated model peptides were in range of values experimentally observed for non-phosphorylated proteins. The extent of the phosphorylation induced changes suggests that the amidic ¹⁵N CSA tensor in phosphorylated proteins could be reasonably well approximated with averaged CSA tensor values experimentally determined for non-phosphorylated amino acids in practical NMR applications, where chemical surrounding of the phosphorylated site is not known a priori in majority of cases. Our calculations provide estimates of relative errors to be associated with the averaged CSA tensor values in interpretations of NMR data from phosphorylated proteins.     

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.     

Nuclear magnetic resonance of the filamentous bacteriophage fd.

The filamentous bacteriophage fd and its major coat protein are being studied by nuclear magnetic resonance (NMR) spectroscopy. 31P NMR shows that the chemical shielding tensor of the DNA phosphates of fd in solution is only slightly reduced in magnitude by motional averaging, indicating that DNA-protein interactions substantially immobilize the DNA packaged in the virus. There is no evidence of chemical interactions between the DNA backbone and the coat protein, since experiments on solid virus show the 31P resonances to have the same principle elements of its chemical shielding tensor as DNA. 1H and 13C NMR spectra of fd virus in solution indicate that the coat proteins are held rigidly in the structure except for some aliphatic side chains that undergo relatively rapid rotations. The presence of limited mobility in the viral coat proteins is substantiated by finding large quadrupole splittings in 2H NMR of deuterium labeled virions. The structure of the coat protein in a lipid environment differs significantly from that found for the assembled virus. Data from 1H and 13C NMR chemical shifts, amide proton exchange rates, and 13C relaxation measurements show that the coat protein in sodium dodecyl sulfate micelles has a native folded structure that varies from that of a typical globular protein or the coat protein in the virus by having a partially flexible backbone and some rapidly rotating aromatic rings.     

Quantification of protein backbone hydrogen-deuterium exchange rates by solid state NMR spectroscopy

We present the quantification of backbone amide hydrogen-deuterium exchange rates (HDX) for immobilized proteins. The experiments make use of the deuterium isotope effect on the amide nitrogen chemical shift, as well as on proton dilution by deuteration. We find that backbone amides in the microcrystalline α-spectrin SH3 domain exchange rather slowly with the solvent (with exchange rates negligible within the individual ¹⁵N–T ₁ timescales). We observed chemical exchange for 6 residues with HDX exchange rates in the range from 0.2 to 5 s⁻¹. Backbone amide ¹⁵N longitudinal relaxation times that we determined previously are not significantly affected for most residues, yielding no systematic artifacts upon quantification of backbone dynamics (Chevelkov et al. 2008b). Significant exchange was observed for the backbone amides of R21, S36 and K60, as well as for the sidechain amides of N38, N35 and for W41ε. These residues could not be fit in our previous motional analysis, demonstrating that amide proton chemical exchange needs to be considered in the analysis of protein dynamics in the solid-state, in case D₂O is employed as a solvent for sample preparation. Due to the intrinsically long ¹⁵N relaxation times in the solid-state, the approach proposed here can expand the range of accessible HDX rates in the intermediate regime that is not accessible so far with exchange quench and MEXICO type experiments.     

Cross correlations between 13C-1H dipolar interactions and 15N chemical shift anisotropy in nucleic acids

Two sets of cross-correlated relaxation rates involving chemical shift anisotropy and dipolar interactions have been measured in an RNA kissing complex. In one case, both the CSA and dipolar interaction tensors are located on the same nucleotide base and are rigidly fixed with respect to each other. In the other case, the CSA tensor is located on the nucleotide base whereas the dipolar interaction is located on the adjoining ribose unit. Analysis of the measured rates in terms of isotropic or anisotropic rotational diffusion has been carried out for both cases. A marked difference between the two models is observed for the cross-correlation rates involving rigidly fixed spin interactions. The influence of internal motions about the glycosidic linkage between the nucleotide base and the ribose unit on cross-correlated relaxation rates has been estimated by applying a model of restricted rotational diffusion. Local motions seem to have a more pronounced effect on cross-correlated relaxation rates when the two spin interactions are not rigidly fixed with respect to each other.     

Comparison of fast backbone dynamics at amide nitrogen and carbonyl sites in dematin headpiece C-terminal domain and its S74E mutant

We perform a detailed comparison of fast backbone dynamics probed at amide nitrogen versus carbonyl carbon sites for dematin headpiece C-terminal domain (DHP) and its S74E mutant (DHPS74E). Carbonyl dynamics is probed via auto-correlated longitudinal rates and transverse C′/C′-C^α CSA/dipolar and C′/C′–N CSA/dipolar cross-correlated rates, while ¹⁵N data are taken from a previous study. Resulting values of effective order parameters and internal correlation times support the conclusion that C′ relaxation reports on a different subset of fast motions compared to those probed at N–H bond vectors in the same peptide planes. ¹³C′ order parameters are on the average 0.08 lower than ¹⁵N order parameters with the exception of the flexible loop region in DHP. The reduction of mobility in the loop region upon the S74E mutation can be seen from the ¹⁵N order parameters but not from the ¹³C order parameters. Internal correlation times at ¹³C′ sites are on the average an order of magnitude longer than those at ¹⁵N sites for the well-structured C-terminal subdomains, while the more flexible N-terminal subdomains have more comparable average internal correlation times.     

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.     

Theoretical investigation of hydrogen bonding effects on oxygen, nitrogen, and hydrogen chemical shielding and electric field gradient tensors of chitosan/HI salt

A density functional theory study has been carried out to calculate the (17)O, (15)N, (13)C, and (1)H chemical shielding as well as (17)O, (14)N, and (2)H electric field gradient tensors of chitosan/HI type I salt. These calculations were performed using the B3LYP functional and 6-311++G (d,p) and 6-31++G (d,p) basis sets. Calculated EFG and chemical shielding tensors were used to evaluate the (17)O, (14)N, and (2)H nuclear quadruple resonance, NQR, and (17)O, (15)N, (13)C, and (1)H nuclear magnetic resonance, NMR, parameters in the cluster model, which are in good agreement with the available experimental data. The difference in the isotropic shielding (sigma(iso)) and quadrupole coupling constant (C(Q)) between monomer and target molecule in the cluster was analyzed in detail. It was shown that both EFG and CS tensors are sensitive to hydrogen-bonding interactions, and calculating both tensors is an advantage. A different influence of various hydrogen bond types, N-Hcdots, three dots, centeredI, O-Hcdots, three dots, centeredI, and N-Hcdots, three dots, centeredO was observed on the calculated CS and EFG tensors. On the basis of this study, nitrogen and O-6 are the most important nuclei to confirm crystalline structure of chitosan/HI. These nuclei have large change in their CS and EFG tensors because of forming intermolecular hydrogen bonds. Moreover, the quantum chemical calculations indicated that the intermolecular hydrogen-bonding interactions play an essential role in determining the relative orientation of CS and EFG tensors of O-6 and nitrogen atoms in the molecular frame axes.     

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.     

Measurement of ¹⁵N relaxation rates in perdeuterated proteins by TROSY-based methods

While extracting dynamics parameters from backbone (15)N relaxation measurements in proteins has become routine over the past two decades, it is increasingly recognized that accurate quantitative analysis can remain limited by the potential presence of systematic errors associated with the measurement of (15)N R(1) and R(2) or R(1ρ) relaxation rates as well as heteronuclear (15)N-{(1)H} NOE values. We show that systematic errors in such measurements can be far larger than the statistical error derived from either the observed signal-to-noise ratio, or from the reproducibility of the measurement. Unless special precautions are taken, the problem of systematic errors is shown to be particularly acute in perdeuterated systems, and even more so when TROSY instead of HSQC elements are used to read out the (15)N magnetization through the NMR-sensitive (1)H nucleus. A discussion of the most common sources of systematic errors is presented, as well as TROSY-based pulse schemes that appear free of systematic errors to the level of <1 %. Application to the small perdeuterated protein GB3, which yields exceptionally high S/N and therefore is an ideal test molecule for detection of systematic errors, yields relaxation rates that show considerably less residue by residue variation than previous measurements. Measured R(2)'/R(1)' ratios fit an axially symmetric diffusion tensor with a Pearson's correlation coefficient of 0.97, comparable to fits obtained for backbone amide RDCs to the Saupe matrix.     

Simultaneous measurement of 1H�C15N and Methyl 1Hm�C13Cm residual dipolar couplings in large proteins

A two-dimensional TROSY-based SIM-(13)C(m)-(1)H(m)/(1)H-(15)N NMR experiment for simultaneous measurements of methyl (1) D (CH) and backbone amide (1) D (NH) residual dipolar couplings (RDC) in {U-[(15)N,(2)H]; Ileδ1-[(13)CH(3)]; Leu,Val-[(13)CH(3)/(12)CD(3)]}-labeled samples of large proteins is described. Significant variation in the alignment tensor of the 82-kDa enzyme Malate synthase G is observed as a function of only slight changes in experimental conditions. The SIM-(13)C(m)-(1)H(m)/(1)H-(15)N data sets provide convenient means of establishing the alignment tensor characteristics via the measurement of (1) D (NH) RDCs in the same protein sample.     

Probing Slow Backbone Dynamics in Proteins Using TROSY-Based Experiments to Detect Cross-Correlated Time-Modulation of Isotropic Chemical Shifts

The difference in the relaxation rates of zero-quantum (ZQ) and double-quantum (DQ) coherences is the result of three principal mechanisms. These include the cross-correlation between the chemical shift anisotropies of the two participating nuclei, dipolar interactions with remote protons as well as interference effects due to the time-modulation of their isotropic chemical shifts as a consequence of slow micros-ms dynamics. The last effect when present, dominates the others resulting in large differences between the relaxation rates of ZQ and DQ coherences. We present here four sets of TROSY-based (Salzmann et al., 1998) experiments that measure this effect for several pairs of backbone nuclei including (15)N, (13)C(alpha) and (13)C'. These experiments allow the detection of the presence of slow dynamic processes in the protein backbone including correlated motion over two and three bonds. Further, we define a new parameter chi which represents the extent of correlated motion on the slow (micros-ms) timescale. This methodology has been applied to (15)N,(13)C,REDPRO-(2)H-labeled (Shekhtman et al., 2002) human ubiquitin. The ubiquitin backbone is seen to exhibit extensive dynamics on the slow timescale. This is most pronounced in several residues in the N-terminal region of the alpha-helix and in the loop connecting the strands beta(4) and beta(5). These residues which include Glu24, Asn25, Glu51 and Asp 52 form a continuous surface. As an additional benefit, the measured rates confirm the dependence of the (13)C(alpha) chemical shift tensor on local secondary structure of the protein backbone.     

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.     

Backbone dynamics measurements on leukemia inhibitory factor, a rigid four-helical bundle cytokine

The backbone dynamics of the four-helical bundle cytokine leukemia inhibitory factor (LIF) have been investigated using 15N NMR relaxation and amide proton exchange measurements on a murine-human chimera, MH35-LIF. For rapid backbone motions (on a time scale of 10 ps to 100 ns), as probed by 15N relaxation measurements, the dynamics parameters were calculated using the model-free formalism incorporating the model selection approach. The principal components of the inertia tensor of MH35-LIF, as calculated from its NMR structure, were 1:0.98:0.38. The global rotational motion of the molecule was, therefore, assumed to be axially symmetric in the analysis of its relaxation data. This yielded a diffusion anisotropy D(parallel)/D(perpendicular) of 1.31 and an effective correlation time (4D(perpendicular) + 2D(parallel))(-1) of 8.9 ns. The average values of the order parameters (S2) for the four helices, the long interhelical loops, and the N-terminus were 0.91, 0.84, and 0.65, respectively, indicating that LIF is fairly rigid in solution, except at the N-terminus. The S2 values for the long interhelical loops of MH35-LIF were higher than those of their counterparts in short-chain members of the four-helical bundle cytokine family. Residues involved in LIF receptor binding showed no consistent pattern of backbone mobilities, with S2 values ranging from 0.71 to 0.95, but residues contributing to receptor binding site III had relatively lower S2 values, implying higher amplitude motions than for the backbone of sites I and II. In the relatively slow motion regime, backbone amide exchange measurements showed that a number of amides from the helical bundle exchanged extremely slowly, persisting for several months in 2H2O at 37 degrees C. Evidence for local unfolding was considered, and correlations among various structure-related parameters and the backbone amide exchange rates were examined. Both sets of data concur in showing that LIF is one of the most rigid four-helical bundle cytokines.     

Density functional calculations of chemical shielding of backbone <Superscript>15</Superscript>N in helical residues of protein G

We performed density functional calculations of backbone ¹⁵N chemical shielding tensors in selected helical residues of protein G. Here we describe a computationally efficient methodology to include most of the important effects in the calculation of chemical shieldings of backbone ¹⁵N. We analyzed the role of long-range intra-protein electrostatic interactions by comparing models with different complexity in vacuum and in charge field. Our results show that the dipole moment of the α-helix can cause significant deshielding of ¹⁵N; therefore, it needs to be considered when calculating ¹⁵N chemical shielding. We found that it is important to include interactions with the side chains that are close in space when the charged form for ionizable side chains is adopted in the calculation. We also illustrate how the ionization state of these side chains can affect the chemical shielding tensor elements. Chemical shielding calculations using a 8-residue fragment model in vacuum and adopting the charged form of ionizable side chains yield a generally good agreement with experimental data.     

Protein global fold determination using site-directed spin and isotope labeling

We describe a simple experimental approach for the rapid determination of protein global folds. This strategy utilizes site-directed spin labeling (SDSL) in combination with isotope enrichment to determine long-range distance restraints between amide protons and the unpaired electron of a nitroxide spin label using the paramagnetic effect on relaxation rates. The precision and accuracy of calculating a protein global fold from only paramagnetic effects have been demonstrated on barnase, a well-characterized protein. Two monocysteine derivatives of barnase, (H102C) and (H102A/Q15C), were 15N enriched, and the paramagnetic nitroxide spin label, MTSSL, attached to the single Cys residue of each. Measurement of amide 1H longitudinal relaxation times, in both the oxidized and reduced states, allowed the determination of the paramagnetic contribution to the relaxation processes. Correlation times were obtained from the frequency dependence of these relaxation processes at 800, 600, and 500 MHz. Distances in the range of 8 to 35 A were calculated from the magnitude of the paramagnetic contribution to the relaxation processes and individual amide 1H correlation times. Distance restraints from the nitroxide spin to amide protons were used as restraints in structure calculations. Using nitroxide to amide 1H distances as long-range restraints and known secondary structure restraints, barnase global folds were calculated having backbone RMSDs <3 A from the crystal structure. This approach makes it possible to rapidly obtain the overall topology of a protein using a limited number of paramagnetic distance restraints.