NMRPipe: a multidimensional spectral processing system based on UNIX pipes

The NMRPipe system is a UNIX software environment of processing, graphics, and analysis tools designed to meet current routine and research-oriented multidimensional processing requirements, and to anticipate and accommodate future demands and developments. The system is based on UNIX pipes, which allow programs running simultaneously to exchange streams of data under user control. In an NMRPipe processing scheme, a stream of spectral data flows through a pipeline of processing programs, each of which performs one component of the overall scheme, such as Fourier transformation or linear prediction. Complete multidimensional processing schemes are constructed as simple UNIX shell scripts. The processing modules themselves maintain and exploit accurate records of data sizes, detection modes, and calibration information in all dimensions, so that schemes can be constructed without the need to explicitly define or anticipate data sizes or storage details of real and imaginary channels during processing. The asynchronous pipeline scheme provides other substantial advantages, including high flexibility, favorable processing speeds, choice of both all-in-memory and disk-bound processing, easy adaptation to different data formats, simpler software development and maintenance, and the ability to distribute processing tasks on multi-CPU computers and computer networks.

Protein backbone angle restraints from searching a database for chemical shift and sequence homology

Chemical shifts of backbone atoms in proteins are exquisitely sensitive to local conformation, and homologous proteins show quite similar patterns of secondary chemical shifts. The inverse of this relation is used to search a database for triplets of adjacent residues with secondary chemical shifts and sequence similarity which provide the best match to the query triplet of interest. The database contains 13C alpha, 13C beta, 13C', 1H alpha and 15N chemical shifts for 20 proteins for which a high resolution X-ray structure is available. The computer program TALOS was developed to search this database for strings of residues with chemical shift and residue type homology. The relative importance of the weighting factors attached to the secondary chemical shifts of the five types of resonances relative to that of sequence similarity was optimized empirically. TALOS yields the 10 triplets which have the closest similarity in secondary chemical shift and amino acid sequence to those of the query sequence. If the central residues in these 10 triplets exhibit similar phi and psi backbone angles, their averages can reliably be used as angular restraints for the protein whose structure is being studied. Tests carried out for proteins of known structure indicate that the root-mean-square difference (rmsd) between the output of TALOS and the X-ray derived backbone angles is about 15 degrees. Approximately 3% of the predictions made by TALOS are found to be in error.

TALOS+: a hybrid method for predicting protein backbone torsion angles from NMR chemical shifts

NMR chemical shifts in proteins depend strongly on local structure. The program TALOS establishes an empirical relation between 13C, 15N and 1H chemical shifts and backbone torsion angles phi and psi (Cornilescu et al. J Biomol NMR 13 289-302, 1999). Extension of the original 20-protein database to 200 proteins increased the fraction of residues for which backbone angles could be predicted from 65 to 74%, while reducing the error rate from 3 to 2.5%. Addition of a two-layer neural network filter to the database fragment selection process forms the basis for a new program, TALOS+, which further enhances the prediction rate to 88.5%, without increasing the error rate. Excluding the 2.5% of residues for which TALOS+ makes predictions that strongly differ from those observed in the crystalline state, the accuracy of predicted phi and psi angles, equals +/-13 degrees . Large discrepancies between predictions and crystal structures are primarily limited to loop regions, and for the few cases where multiple X-ray structures are available such residues are often found in different states in the different structures. The TALOS+ output includes predictions for individual residues with missing chemical shifts, and the neural network component of the program also predicts secondary structure with good accuracy.

Backbone dynamics of proteins as studied by 15N inverse detected heteronuclear NMR spectroscopy: application to staphylococcal nuclease

This paper describes the use of novel two-dimensional nuclear magnetic resonance (NMR) pulse sequences to provide insight into protein dynamics. The sequences developed permit the measurement of the relaxation properties of individual nuclei in macromolecules, thereby providing a powerful experimental approach to the study of local protein mobility. For isotopically labeled macromolecules, the sequences enable measurements of heteronuclear nuclear Overhauser effects (NOE) and spin-lattice (T1) and spin-spin (T2) 15N or 13C relaxation times with a sensitivity similar to those of many homonuclear 1H experiments. Because T1 values and heteronuclear NOEs are sensitive to high-frequency motions (10(8)-10(12) s-1) while T2 values are also a function of much slower processes, it is possible to explore dynamic events occurring over a large time scale. We have applied these techniques to investigate the backbone dynamics of the protein staphylococcal nuclease (S. Nase) complexed with thymidine 3',5'-bisphosphate (pdTp) and Ca2+ and labeled uniformly with 15N. T1, T2, and NOE values were obtained for over 100 assigned backbone amide nitrogens in the protein. Values of the order parameter (S), characterizing the extent of rapid 1H-15N bond motions, have been determined. These results suggest that there is no correlation between these rapid small amplitude motions and secondary structure for S. Nase. In contrast, 15N line widths suggest a possible correlation between secondary structure and motions on the millisecond time scale. In particular, the loop region between residues 42 and 56 appears to be considerably more flexible on this slow time scale than the rest of the protein.

Direct measurement of distances and angles in biomolecules by NMR in a dilute liquid crystalline medium

In isotropic solution, internuclear dipolar couplings average to zero as a result of rotational diffusion. By dissolving macromolecules in a dilute aqueous nematic discotic liquid-crystalline medium containing widely spaced magnetically oriented particles, a tunable degree of solute alignment with the magnetic field can be created while retaining the high resolution and sensitivity of the regular isotropic nuclear magnetic resonance (NMR) spectrum. Dipolar couplings between 1H-1H, 1H-13C, 1H-15N, and 13C-13C pairs in such an oriented macromolecule no longer average to zero, and are readily measured. Distances and angles derived from dipolar couplings in human ubiquitin are in excellent agreement with its crystal structure. The approach promises to improve the accuracy of structures determined by NMR, and extend the size limit.

Solution structure of a calmodulin-target peptide complex by multidimensional NMR

The three-dimensional solution structure of the complex between calcium-bound calmodulin (Ca(2+)-CaM) and a 26-residue synthetic peptide comprising the CaM binding domain (residues 577 to 602) of skeletal muscle myosin light chain kinase, has been determined using multidimensional heteronuclear filtered and separated nuclear magnetic resonance spectroscopy. The two domains of CaM (residues 6 to 73 and 83 to 146) remain essentially unchanged upon complexation. The long central helix (residues 65 to 93), however, which connects the two domains in the crystal structure of Ca(2+)-CaM, is disrupted into two helices connected by a long flexible loop (residues 74 to 82), thereby enabling the two domains to clamp residues 3 to 21 of the bound peptide, which adopt a helical conformation. The overall structure of the complex is globular, approximating an ellipsoid of dimensions 47 by 32 by 30 angstroms. The helical peptide is located in a hydrophobic channel that passes through the center of the ellipsoid at an angle of approximately 45 degrees with its long axis. The complex is mainly stabilized by hydrophobic interactions which, from the CaM side, involve an unusually large number of methionines. Key residues of the peptide are Trp4 and Phe17, which serve to anchor the amino- and carboxyl-terminal halves of the peptide to the carboxyl- and amino-terminal domains of CaM, respectively. Sequence comparisons indicate that a number of peptides that bind CaM with high affinity share this common feature containing either aromatic residues or long-chain hydrophobic ones separated by a stretch of 12 residues, suggesting that they interact with CaM in a similar manner.

Protein backbone and sidechain torsion angles predicted from NMR chemical shifts using artificial neural networks

A new program, TALOS-N, is introduced for predicting protein backbone torsion angles from NMR chemical shifts. The program relies far more extensively on the use of trained artificial neural networks than its predecessor, TALOS+. Validation on an independent set of proteins indicates that backbone torsion angles can be predicted for a larger, ≥90 % fraction of the residues, with an error rate smaller than ca 3.5 %, using an acceptance criterion that is nearly two-fold tighter than that used previously, and a root mean square difference between predicted and crystallographically observed (ϕ, ψ) torsion angles of ca 12º. TALOS-N also reports sidechain χ(1) rotameric states for about 50 % of the residues, and a consistency with reference structures of 89 %. The program includes a neural network trained to identify secondary structure from residue sequence and chemical shifts.

Measurement of J and dipolar couplings from simplified two-dimensional NMR spectra

Simple procedures are described for recording complementary in-phase and antiphase J-coupled NMR spectra. The sum and difference of these spectra contain only the upfield and the downfield components of a doublet, making it possible to measure the J splitting directly from these combinations without an increase in resonance overlap relative to the decoupled spectrum. The approach is demonstrated for measurement of 1JNH splittings and 2JHNC splittings in oriented and isotropic ubiquitin. Dipolar couplings obtained from differences in the splittings measured in the oriented and isotropic phases are in excellent agreement with dipolar couplings obtained from direct measurement of the splitting or from a conventional E. COSY-type measurement.

Structure and dynamics of micelle-bound human alpha-synuclein

Misfolding of the protein alpha-synuclein (aS), which associates with presynaptic vesicles, has been implicated in the molecular chain of events leading to Parkinson's disease. Here, the structure and dynamics of micelle-bound aS are reported. Val3-Val37 and Lys45-Thr92 form curved alpha-helices, connected by a well ordered, extended linker in an unexpected anti-parallel arrangement, followed by another short extended region (Gly93-Lys97), overlapping the recently identified chaperone-mediated autophagy recognition motif and a highly mobile tail (Asp98-Ala140). Helix curvature is significantly less than predicted based on the native micelle shape, indicating a deformation of the micelle by aS. Structural and dynamic parameters show a reduced helical content for Ala30-Val37. A dynamic variation in interhelical distance on the microsecond timescale is complemented by enhanced sub-nanosecond timescale dynamics, particularly in the remarkably glycine-rich segments of the helices. These unusually rich dynamics may serve to mitigate the effect of aS binding on membrane fluidity. The well ordered conformation of the helix-helix connector indicates a defined interaction with lipidic surfaces, suggesting that, when bound to larger diameter synaptic vesicles, it can act as a switch between this structure and a previously proposed uninterrupted helix.

Backbone dynamics of calmodulin studied by 15N relaxation using inverse detected two-dimensional NMR spectroscopy: the central helix is flexible

The backbone dynamics of Ca(2+)-saturated recombinant Drosophila calmodulin has been studied by 15N longitudinal and transverse relaxation experiments, combined with 15N(1H) NOE measurements. Results indicate a high degree of mobility near the middle of the central helix of calmodulin, from residue K77 through S81, with order parameters (S2) in the 0.5-0.6 range. The anisotropy observed in the motion of the two globular calmodulin domains is much smaller than expected on the basis of hydrodynamic calculations for a rigid dumbbell type structure. This indicates that, for the purposes of 15N relaxation, the tumbling of the N-terminal (L4-K77) and C-terminal (E82-S147) lobes of calmodulin is effectively independent. A slightly shorter motional correlation time (tau c approximately 6.3 ns) is obtained for the C-terminal domain compared to the N-terminal domain (tau c approximately 7.1 ns), in agreement with the smaller size of the C-terminal domain. A high degree of mobility, with order parameters of approximately 0.5, is also observed in the loop that connects the first with the second EF-hand type calcium binding domain and in the loop connecting the third and fourth calcium binding domain.

Overcoming the overlap problem in the assignment of 1H NMR spectra of larger proteins by use of three-dimensional heteronuclear 1H-15N Hartmann-Hahn-multiple quantum coherence and nuclear Overhauser-multiple quantum coherence spectroscopy: application to interleukin 1 beta

The application of three-dimensional (3D) heteronuclear NMR spectroscopy to the sequential assignment of the 1H NMR spectra of larger proteins is presented, using uniformly labeled (approximately 95%) [15N]interleukin 1 beta, a protein of 153 residues and molecular mass of 17.4 kDa, as an example. The two-dimensional (2D) 600-MHz spectra of interleukin 1 beta are too complex for complete analysis, owing to extensive cross-peak overlap and chemical shift degeneracy. We show that the combined use of 3D 1H-15N Hartmann-Hahn-multiple quantum coherence (HOHAHA-HMQC) and nuclear Overhauser-multiple quantum coherence (NOESY-HMQC) spectroscopy, designed to provide the necessary through-bond and through-space correlations for sequential assignment, provides a practical general-purpose method for resolving ambiguities which severely limit the analysis of conventional 2D NMR spectra. The absence of overlapping cross-peaks in these 3D spectra allows the unambiguous identification of C alpha H(i)-NH(i+1) and NH(i)-NH(i+1) through-space nuclear Overhauser connectivities necessary for connecting a particular C alpha H(i)-NH(i) through-bond correlation with its associated through-space sequential cross-peak The problem of amide NH chemical shift degeneracy in the 1H NMR spectrum is therefore effectively removed, and the assignment procedure simply involves inspecting a series of 2D 1H-1H slices edited by the chemical shift of the directly bonded 15N atom. Connections between residues can be identified almost without any knowledge of the spin system types involved, though this type of information is clearly required for the eventual placement of the connected residues within the primary sequence.

A novel approach for sequential assignment of 1H, 13C, and 15N spectra of proteins: heteronuclear triple-resonance three-dimensional NMR spectroscopy. Application to calmodulin

A novel approach is described for obtaining sequential assignment of the backbone 1H, 13C, and 15N resonances of larger proteins. The approach is demonstrated for the protein calmodulin (16.7 kDa), uniformly (approximately 95%) labeled with 15N and 13C. Sequential assignment of the backbone residues by standard methods was not possible because of the very narrow chemical shift distribution range of both NH and C alpha H protons in this largely alpha-helical protein. We demonstrate that the combined use of four new types of heteronuclear 3D NMR spectra together with the previously described HOHAHA-HMQC 3D experiment [Marion, D., et al. (1989) Biochemistry 28, 6150-6156] can provide unambiguous sequential assignment of protein backbone resonances. Sequential connectivity is derived from one-bond J couplings and the procedure is therefore independent of the backbone conformation. All the new 3D NMR experiments use 1H detection and rely on multiple-step magnetization transfers via well-resolved one-bond J couplings, offering high sensitivity and requiring a total of only 9 days for the recording of all five 3D spectra. Because the combination of 3D spectra offers at least two and often three independent pathways for determining sequential connectivity, the new assignment procedure is easily automated. Complete assignments are reported for the proton, carbon, and nitrogen backbone resonances of calmodulin, complexed with calcium.

Solution structure of calcium-free calmodulin

The three-dimensional structure of calmodulin in the absence of Ca2+ has been determined by three- and four-dimensional heteronuclear NMR experiments, including ROE, isotope-filtering combined with reverse labelling, and measurement of more than 700 three-bond J-couplings. In analogy with the Ca(2+)-ligated state of this protein, it consists of two small globular domains separated by a flexible linker, with no stable, direct contacts between the two domains. In the absence of Ca2+, the four helices in each of the two globular domains form a highly twisted bundle, capped by a short anti-parallel beta-sheet. This arrangement is qualitatively similar to that observed in the crystal structure of the Ca(2+)-free N-terminal domain of troponin C.