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.     

Direct structure refinement of high molecular weight proteins against residual dipolar couplings and carbonyl chemical shift changes upon alignment: an application to maltose binding protein

The global fold of maltose binding protein in complex with -cyclodextrin has been determined using a CNS-based torsion angle molecular dynamics protocol involving direct refinement against dipolar couplings and carbonyl chemical shift changes that occur upon alignment. The shift changes have been included as structural restraints using a new module, CANI, that has been incorporated into CNS. Force constants and timesteps have been determined that are particularly effective in structure refinement applications involving high molecular weight proteins with small to moderate numbers of NOE restraints. Solution structures of the N- and C-domains of MBP calculated with this new protocol are within 2 Å of the X-ray conformation.     

Overall structure and sugar dynamics of a DNA dodecamer from homo- and heteronuclear dipolar couplings and 31P chemical shift anisotropy

The solution structure of d(CGCGAATTCGCG)₂ has been determined on the basis of an exceptionally large set of residual dipolar couplings. In addition to the heteronuclear ¹³C-¹H and ¹⁵N-¹H and qualitative homonuclear ¹H-¹H dipolar couplings, previously measured in bicelle medium, more than 300 quantitative ¹H-¹H and 22 ³¹P-¹H dipolar restraints were obtained in liquid crystalline Pf1 medium, and 22 ³¹P chemical shift anisotropy restraints. High quality DNA structures can be obtained solely on the basis of these new restraints, and these structures are in close agreement with those calculated previously on the basis of ¹³C-¹H and ¹⁵N-¹H dipolar couplings. In the newly calculated structures, ³¹P-¹H dipolar and ³Jsub H3 P sub couplings and ³¹P CSA data restrain the phosphodiester backbone torsion angles. The final structure represents a quite regular B-form helix with a modest bending of 10°, which is essentially independent of whether or not electrostatic terms are used in the calculation. Combined, the number of homo- and heteronuclear dipolar couplings significantly exceeds the number of degrees of freedom in the system. Results indicate that the dipolar coupling data cannot be fit by a single structure, but are compatible with the presence of rapid equilibria between C2-endo and C3-endo deoxyribose puckers (sugar switching). The C2-H2/H2 dipolar couplings in B-form DNA are particularly sensitive to sugar pucker and yield the largest discrepancies when fit to a single structure. To resolve these discrepancies, we suggest a simplified dipolar coupling analysis that yields N/S equilibria for the ribose sugar puckers, which are in good agreement with previous analyses of NMR J_HH couplings, with a population of the minor C3-endo form higher for pyrimidines than for purines.     

Network of long-range concerted chemical shift displacements upon ligand binding to human angiogenin

Molecular recognition models of both induced fit and conformational selection rely on coupled networks of flexible residues and/or structural rearrangements to promote protein function. While the atomic details of these motional events still remain elusive, members of the pancreatic ribonuclease superfamily were previously shown to depend on subtle conformational heterogeneity for optimal catalytic function. Human angiogenin, a structural homologue of bovine pancreatic RNase A, induces blood vessel formation and relies on a weak yet functionally mandatory ribonucleolytic activity to promote neovascularization. Here, we use the NMR chemical shift projection analysis (CHESPA) to clarify the mechanism of ligand binding in human angiogenin, further providing information on long-range intramolecular residue networks potentially involved in the function of this enzyme. We identify two main clusters of residue networks displaying correlated linear chemical shift trajectories upon binding of substrate fragments to the purine- and pyrimidine-specific subsites of the catalytic cleft. A large correlated residue network clusters in the region corresponding to the V₁ domain, a site generally associated with the angiogenic response and structural stability of the enzyme. Another correlated network (residues 40–42) negatively affects the catalytic activity but also increases the angiogenic activity. ¹⁵N-CPMG relaxation dispersion experiments could not reveal the existence of millisecond timescale conformational exchange in this enzyme, a lack of flexibility supported by the very low-binding affinities and catalytic activity of angiogenin. Altogether, the current report potentially highlights the existence of long-range dynamic reorganization of the structure upon distinct subsite binding events in human angiogenin.     

De novo protein structure generation from incomplete chemical shift assignments

NMR chemical shifts provide important local structural information for proteins. Consistent structure generation from NMR chemical shift data has recently become feasible for proteins with sizes of up to 130 residues, and such structures are of a quality comparable to those obtained with the standard NMR protocol. This study investigates the influence of the completeness of chemical shift assignments on structures generated from chemical shifts. The Chemical-Shift-Rosetta (CS-Rosetta) protocol was used for de novo protein structure generation with various degrees of completeness of the chemical shift assignment, simulated by omission of entries in the experimental chemical shift data previously used for the initial demonstration of the CS-Rosetta approach. In addition, a new CS-Rosetta protocol is described that improves robustness of the method for proteins with missing or erroneous NMR chemical shift input data. This strategy, which uses traditional Rosetta for pre-filtering of the fragment selection process, is demonstrated for two paramagnetic proteins and also for two proteins with solid-state NMR chemical shift assignments. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

A simple method to adjust inconsistently referenced <Superscript>13</Superscript>C and <Superscript>15</Superscript>N chemical shift assignments of proteins

Inconsistent ¹³C and ¹⁵N chemical shift referencing is a continuing problem associated with protein chemical shift assignments deposited in BioMagResBank (BMRB). Here we describe a simple and robust approach that can quantitatively determine the ¹³C and ¹⁵N referencing offsets solely from chemical shift assignment data and independently of 3D coordinate data. This novel structure-independent approach permitted the assessment and determination of ¹³C and ¹⁵N reference offsets for all protein entries deposited in the BMRB. Tests on 452 proteins with known 3D structures show that this structure-independent approach yields ¹³C and ¹⁵N referencing offsets that exhibit excellent agreement with those calculated on the basis of 3D structures. Furthermore, this protocol appears to improve the accuracy of chemical shift-derived secondary structural identification, and has been formally incorporated into a computer program called PSSI (http//www.pronmr.com).Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1007/s10858-004-7441-3     

An evaluation of chemical shift index-based secondary structure determination in proteins: Influence of random coil chemical shifts

Random coil chemical shifts are commonly used to detect protein secondary structural elements in chemical shift index (CSI) calculations. Though this technique is widely used and seems reliable for folded proteins, the choice of reference random coil chemical shift values can significantly alter the outcome of secondary structure estimation. In order to evaluate these effects, we present a comparison of secondary structure content calculated using CSI, based on five different reference random coil chemical shift value sets, to that derived from three-dimensional structures.Our results show that none of the reference random coil data sets chosen for evaluation fully reproduces the actual secondary structures. Among the reference values generally available to date, most tend to be good estimators only of helices. Based on our evaluation, we recommend the experimental values measured by Schwarzinger et al.(2000), and statistical values obtained by Lukin et al. (1997), as good estimators of both helical and sheet content.     

Effects of Protein–pheromone Complexation on Correlated Chemical Shift Modulations

Major urinary protein (MUP) is a pheromone-carrying protein of the lipocalin family. Previous studies by isothermal titration calorimetry (ITC) show that the affinity of MUP for the pheromone 2-methoxy-3-isobutylpyrazine (IBMP) is mainly driven by enthalpy, with a small unfavourable entropic contribution. Entropic terms can be attributed in part to changes in internal motions of the protein upon binding. Slow internal motions can lead to correlated or anti-correlated modulations of the isotropic chemical shifts of carbonyl C′ and amide N nuclei. Correlated chemical shift modulations (CSM/CSM) in MUP have been determined by measuring differences of the transverse relaxation rates of zero- and double-quantum coherences ZQC{C′N} and DQC{C′N}, and by accounting for the effects of correlated fluctuations of dipole–dipole couplings (DD/DD) and chemical shift anisotropies (CSA/CSA). The latter can be predicted from tensor parameters of C′ and N nuclei that have been determined in earlier work. The effects of complexation on slow time-scale protein dynamics can be determined by comparing the temperature dependence of the relaxation rates of APO-MUP (i.e., without ligand) and HOLO-MUP (i.e., with IBMP as a ligand). Electronic supplementary material Electronic supplementary material is available for this article at and accessible for authorised users.     

13C solid-state NMR chemical shift anisotropy analysis of the anomeric carbon in carbohydrates

(13)C NMR solid-state structural analysis of the anomeric center in carbohydrates was performed on six monosaccharides: glucose (Glc), mannose (Man), galactose (Gal), galactosamine hydrochloride (GalN), glucosamine hydrochloride (GlcN), and N-acetyl-glucosamine (GlcNAc). In the 1D (13)C cross-polarization/magic-angle spinning (CP/MAS) spectrum, the anomeric center C-1 of these carbohydrates revealed two well resolved resonances shifted by 3-5ppm, which were readily assigned to the anomeric alpha and beta forms. From this experiment, we also extracted the (13)C chemical shift anisotropy (CSA) tensor elements of the two forms from their spinning sideband intensities, respectively. It was found out that the chemical shift tensor for the alpha anomer was more axially symmetrical than that of the beta form. A strong linear correlation was obtained when the ratio of the axial asymmetry of the (13)C chemical shift tensors of the two anomeric forms was plotted in a semilogarithmic plot against the relative population of the two anomers. Finally, we applied REDOR spectroscopy to discern whether or not there were any differences in the sugar ring conformation between the anomers. Identical two-bond distances of 2.57A (2.48A) were deduced for both the alpha and beta forms in GlcNAc (GlcN), suggesting that the two anomers have essentially identical sugar ring scaffolds in these sugars. In light of these REDOR distance measurements and the strong correlation observed between the ratio of the axial asymmetry parameters of the (13)C chemical shift tensors and the relative population between the two anomeric forms, we concluded that the anomeric effect arises principally from interaction of the electron charge clouds between the C-1-O-5 and the C-1-O-1 bonds in these monosaccharides.     

Estimation of protein secondary structure content directly from NMR spectra using an improved empirical correlation with averaged chemical shift

We have recently shown that the averaged chemical shift (ACS) of a nucleus in the protein backbone correlates well empirically to its secondary structure content (SSC). This allows the estimation of SSC directly from the NMR spectrum without the time intensive process of chemical shift assignment. Here, we present an empirical correlation that accounts both for contributions to the relevant protein and chemical shift databases made subsequent to the original analysis, and for missing or inconsistently referenced resonances. Our results affirm that this method provides a significant tool for initial structural prediction from NMR data prior to complete chemical shift assignment.     

A device for the measurement of residual chemical shift anisotropy and residual dipolar coupling in soluble and membrane-associated proteins

Residual dipolar coupling (RDC) and residual chemical shift anisotropy (RCSA) report on orientational properties of a dipolar bond vector and a chemical shift anisotropy principal axis system, respectively. They can be highly complementary in the analysis of backbone structure and dynamics in proteins as RCSAs generally include a report on vectors out of a peptide plane while RDCs usually report on in-plane vectors. Both RDC and RCSA average to zero in isotropic solutions and require partial orientation in a magnetic field to become observable. While the alignment and measurement of RDC has become routine, that of RCSA is less common. This is partly due to difficulties in providing a suitable isotopic reference spectrum for the measurement of the small chemical shift offsets coming from RCSA. Here we introduce a device (modified NMR tube) specifically designed for accurate measurement of reference and aligned spectra for RCSA measurements, but with a capacity for RDC measurements as well. Applications to both soluble and membrane anchored proteins are illustrated.     

Prediction algorithm for amino acid types with their secondary structure in proteins (PLATON) using chemical shifts

The algorithm PLATON is able to assign sets of chemical shifts derived from a single residue to amino acid types with its secondary structure (amino acid species). A subsequent ranking procedure using optionally two different penalty functions yields predictions for possible amino acid species for the given set of chemical shifts. This was demonstrated in the case of the -spectrin SH3 domain and applied to 9 further protein data sets taken from the BioMagRes database. A database consisting of reference chemical shift patterns (reference CSPs) was generated from assigned chemical shifts of proteins with known 3D-structure. This reference CSP database is used in our approach for extracting distributions of amino acid types with their most likely secondary structure elements (namely -helix, -sheet, and coil) for single amino acids by comparison with query CSPs. Results obtained for the 10 investigated proteins indicates that the percentage of correct amino acid species in the first three positions in the ranking list, ranges from 71.4% to 93.2% for the more favorable penalty function. Where only the top result of the ranking list for these 10 proteins is considered, 36.5% to 83.1% of the amino acid species are correctly predicted. The main advantage of our approach, over other methods that rely on average chemical shift values is the ability to increase database content by incorporating newly derived CSPs, and therefore to improve PLATON's performance over time.     

31P NMR relaxation measurements of the phosphate backbone of a double stranded hexadeoxynucleotide in solution: determination of the chemical shift anisotropy

The five phosphates of the deoxynucleotide d(CpGpTpApCpG)₂ have been assigned by two-dimensional heteronuclear NMR spectroscopy. The chemical shift anisotropy and correlation time of each phosphate group has been determined from measurements of the spin-lattice, spin-spin relaxation rate constants and the ³¹P-{¹H} nuclear Overhauser enhancement (NOE) at three magnetic field strengths (4.7 T, 9.4 T, and 11.75 T) and two temperatures (288 K and 298 K). As expected, the relaxation data require two mechanisms to account for the observed rate constants, i.e. dipole-dipole and chemical shift anisotropy. At 9.4 T and 11.75 T, the latter mechanism dominates the relaxation, leading to insignificant NOE intensities. The correlation time, chemical shift anisotropy and effective P-H distance were obtained from least-squares fitting to the data. Comparison of the fitted value for the correlation time with that obtained from ¹H measurements shows that the molecule behaves essentially as rigid rotor on the nanosecond timescale. Large amplitude motions observed in long segments of DNA are due to bending motions that do not contribute significantly to relaxation in short oligonucleotides.Abbreviations CSA chemical shift anisotropy - NOE nuclear Overhauser enhancement Offprint requests to: A. N. Lane     

Improved simulation of NOESY spectra by RELAX-JT2 including effects of J-coupling, transverse relaxation and chemical shift anisotropy

RELAX-JT2 is an extension of RELAX, a program for the simulation of 1H 2D NOESY spectra and (15)N or (13)C edited 3D NOESY-HSQC spectra of biological macromolecules. In addition to the already existing NOE-simulation it allows the proper simulation of line shapes by the integrated calculation of T(2) times and multiplet structures caused by J-couplings. Additionally the effects of relaxation mediated by chemical shift anisotropy are taken into account. The new routines have been implemented in the program AUREMOL, which aims at the automated NMR structure determination of proteins in solution. For a manual or automatic assignment of experimental spectra that is based on the comparison with the corresponding simulated spectra, the additional line shape information now available is a valuable aid. The new features have been successfully tested with the histidine-containing phosphocarrier protein HPr from Staphylococcus carnosus.     

Secondary chemical shifts in immobilized peptides and proteins: a qualitative basis for structure refinement under magic angle spinning

Resonance assignments recently obtained on immobilized polypeptides and a membrane protein aggregate under Magic Angle Spinning are compared to random coil values in the liquid state. The resulting chemical shift differences (secondary chemical shifts) are evaluated in light of the backbone torsion angle previously reported using X-ray crystallography. In all cases, a remarkable correlation is found suggesting that the concept of secondary chemical shifts, well established in the liquid state, can be of similar importance in the context of multiple-labelled polypeptides studied under MAS conditions.     

Assessing potential bias in the determination of rotational correlation times of proteins by NMR relaxation

The various factors that influence the reliable and efficient determination of the correlation time describing molecular reorientation of proteins by NMR relaxation methods are examined. Nuclear Overhauser effects, spin-lattice, and spin-spin relaxation parameters of 15N NMR relaxation in ubiquitin have been determined at 17.6, 14.1, 11.7 and 9.4 Tesla. This unusually broad set of relaxation parameters has allowed the examination of the influence of chemical shift anisotropy, the functional form of the model-free spectral density, and the reliability of determined spin- spin relaxation parameters on the characterization of global tumbling of the protein. Treating the 15N chemical shift anisotropy (CSA) as an adjustable parameter, a consensus value of –170 ± 15ppm for the breadth of the chemical shift tensor and a global isotropic correlation time of 4.1ns are found when using the model-free spectral density to fit T1 and NOE data from all fields. The inclusion of T2 relaxation parameters in the determination of the global correlation time results in its increase to 4.6ns. This apparent inconsistency may explain a large portion of the discrepancy often found between NMR- and fluorescence-derived m values for proteins. The near identity of observed T2 and T1 values suggests that contributions from slow motions are not the origin of the apparent inconsistency with obtained T1 and NOE data. Various considerations suggest that the origin of this apparent discrepancy may reside in a contribution to the spectral density at zero frequency that is not represented by the simple model-free formalism in addition to the usual experimental difficulties associated with the measurement of these relaxation parameters. Finally, an axially symmetric diffusion tensor for ubiquitin is obtained using exclusively T1 and NOE data. A recommendation is reached on the types and combinations of relaxation data that can be used to reliably determine m values. It is also noted that the reliable determination of m values from 15N T1 and NOE relaxation parameters will become increasingly difficult as m increases.     

Evidence for extensive anisotropic local motions in a small enzyme using a new method to determine NMR cross-correlated relaxation rates in the absence of resolved scalar coupling

Transverse ¹³CO-¹HN (dipole-dipole)/¹³CO (CSA) cross-correlated relaxation rates were measured for the ¹³CO resonances of the protein ribonuclease Binase from Bacillus intermedius (12.3 kDa). This was carried out with a novel E.COSY-type triple-resonance experiment, which allows the measurement of cross-correlated transverse relaxation rate from multiplet effects in the absence of resolved scalar coupling. The ¹³CO-¹HN (dipole-dipole)/¹³CO (CSA) cross-correlated relaxation rates were determined with an average precision of ±5% and cover a range of values between –1.5 and +0.6 Hz. The average (–0.44 Hz) is to be compared with the computed value of –0.83 Hz for this interaction. Mechanisms that potentially can cause the average to be smaller than the theoretical value and the unexpected large spread in observed values are discussed. It is suggested that large contributions to the variations are due to large amplitude local anisotropic motions.     

Orientational constraints as three-dimensional structural constraints from chemical shift anisotropy: the polypeptide backbone of gramicidin A in a lipid bilayer.

Chemical shifts observed from samples that are uniformly aligned with respect to the magnetic field can be used as very high-resolution structural constraints. This constraint takes the form of an orientational constraint rather than the more familiar distance constraint. The accuracy of these constraints is dependent upon the quality of the tensor characterization. Both tensor element magnitudes and tensor orientations with respect to the molecular frame need to be considered. Here these constraints have been used to evaluate models for the channel conformation of gramicidin A. Of the three models used, the one experimentally derived model of gramicidin in sodium dodecyl sulfate micelles fits the data least well.