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.

We present a structural model for amyloid fibrils formed by the 40-residue ␤-amyloid peptide associated with Alzheimer's disease (A␤ 1-40), based on a set of experimental constraints from solid state NMR spectroscopy. The model additionally incorporates the cross-␤ structural motif established by x-ray fiber diffraction and satisfies constraints on A␤ 1-40 fibril dimensions and mass-per-length determined from electron microscopy. Approximately the first 10 residues of A␤ 1-40 are structurally disordered in the fibrils. Residues 12-24 and 30 -40 adopt ␤-strand conformations and form parallel ␤-sheets through intermolecular hydrogen bonding. Residues 25-29 contain a bend of the peptide backbone that brings the two ␤-sheets in contact through sidechain-sidechain interactions. A single cross-␤ unit is then a double-layered ␤-sheet structure with a hydrophobic core and one hydrophobic face. The only charged sidechains in the core are those of D23 and K28, which form salt bridges. Fibrils with minimum mass-per-length and diameter consist of two cross-␤ units with their hydrophobic faces juxtaposed.A myloid fibrils are filamentous structures, with typical diameters of Ϸ10 nm and lengths up to several micrometers, formed by numerous peptides and proteins with disparate sequences and molecular weights. Biomedical interest in amyloid fibrils arises from their occurrence in amyloid diseases (1), including Alzheimer's disease, type 2 diabetes, Huntington's disease, and prion diseases. Current interest in the molecular structures of amyloid fibrils additionally arises from fundamental questions regarding the molecular mechanism of amyloid formation and the nature of the intermolecular interactions that stabilize these structures for an extremely diverse class of polypeptides.No high-resolution molecular structure of an amyloid fibril has yet been determined experimentally because amyloid fibrils are noncrystalline solid materials and are therefore incompatible with x-ray crystallography and liquid state NMR. X-ray fiber diffraction shows that amyloid fibrils contain cross-␤ structural motifs, i.e., extended ␤-sheets in which the ␤-strand segments run approximately perpendicular to, and the intermolecular hydrogen bonds run approximately parallel to, the long axis of the fibril (2, 3). Other molecular-level structural features of amyloid fibrils are not well established.In the case of fibrils formed by the full-length ␤-amyloid peptide associated with Alzheimer's disease (A␤), which ranges from 39 to 43 residues in length in vivo (4, 5), several molecular models have been proposed (6-10). These models exhibit many qualitative and quantitative differences, reflecting the paucity of experimental constraints. All of these models are inconsistent with recent measurements of 13 C-13 C nuclear magnetic dipole-dipole couplings (i.e., intermolecular distances) by solid state NMR (11-13), which imply an in-register parallel alignment of peptide chains within the cross-␤ motif in A␤ 1-40 and A␤ 1-42 fibrils (A␤ mϪn denotes residues m t...

NMR chemical shifts in proteins depend strongly on local structure. The program TALOS establishes an empirical relation between 13 C, 15 N and 1 H chemical shifts and backbone torsion angles φ and ψ (G. 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 percent. Addition of a twolayer 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 φ and ψ angles, equals ±13°. 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.

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.

Protein NMR chemical shifts are highly sensitive to local structure. A robust protocol is described that exploits this relation for de novo protein structure generation, using as input experimental parameters the 13 C ␣ , 13 C ␤ , 13 C , 15 N, 1 H ␣ and 1 H N NMR chemical shifts. These shifts are generally available at the early stage of the traditional NMR structure determination process, before the collection and analysis of structural restraints. The chemical shift based structure determination protocol uses an empirically optimized procedure to select protein fragments from the Protein Data Bank, in conjunction with the standard ROSETTA Monte Carlo assembly and relaxation methods. Evaluation of 16 proteins, varying in size from 56 to 129 residues, yielded full-atom models that have 0.7-1.8 Å root mean square deviations for the backbone atoms relative to the experimentally determined x-ray or NMR structures. The strategy also has been successfully applied in a blind manner to nine protein targets with molecular masses up to 15.4 kDa, whose conventional NMR structure determination was conducted in parallel by the Northeast Structural Genomics Consortium. This protocol potentially provides a new direction for high-throughput NMR structure determination. molecular fragment replacement ͉ protein structure prediction ͉ ROSETTA ͉ structural genomics

No abstract

NMR measurements of a large set of protein backbone one-bond dipolar couplings have been carried out to refine the structure of the third IgG-binding domain of Protein G (GB3), previously solved by X-ray crystallography at a resolution of 1.1 A. Besides the commonly used bicelle, poly(ethylene glycol), and filamentous phage liquid crystalline media, dipolar couplings were also measured when the protein was aligned inside either positively or negatively charged stretched acrylamide gels. Refinement of the GB3 crystal structure against the (13)C(alpha)-(13)C' and (13)C'-(15)N dipolar couplings improves the agreement between experimental and predicted (15)N-(1)H(N) as well as (13)C(alpha)-(1)H(alpha) dipolar couplings. Evaluation of the peptide bond N-H orientations shows a weak anticorrelation between the deviation of the peptide bond torsion angle omega from 180 degrees and the angle between the N-H vector and the C'-N-C(alpha) plane. The slope of this correlation is -1, indicating that, on average, pyramidalization of the peptide N contributes to small deviations from peptide bond planarity ( = 179.3 +/- 3.1 degrees ) to the same degree as true twisting around the C'-N bond. Although hydrogens are commonly built onto crystal structures assuming the N-H vector orientation falls on the line bisecting the C'-N-C(alpha) angle, a better approximation adjusts the C(alpha)-C'-N-H torsion angle to -2 degrees. The (15)N-(1)H(N) dipolar data do not contradict the commonly accepted motional model where angular fluctuations of the N-H bond orthogonal to the peptide plane are larger than in-plane motions, but the amplitude of angular fluctuations orthogonal the C(alpha)(i-1)-N(i)-C(alpha)(i) plane exceeds that of in-plane motions by at most 10-15 degrees. Dipolar coupling analysis indicates that for most of the GB3 backbone, the amide order parameters, S, are highly homogeneous and vary by less than +/-7%. Evaluation of the H(alpha) proton positions indicates that the average C(alpha)-H(alpha) vector orientation deviates by less than 1 degrees from the direction that makes ideal tetrahedral angles with the C(alpha)-C(beta) and C(alpha)-N vectors.

No abstract