Peter Güntert

Protein NMR Structure Determination with Automated NOE Assignment Using the New Software CANDID and the Torsion Angle Dynamics Algorithm DYANA

Combined automated NOE assignment and structure determination module (CANDID) is a new software for efficient NMR structure determination of proteins by automated assignment of the NOESY spectra. CANDID uses an iterative approach with multiple cycles of NOE cross-peak assignment and protein structure calculation using the fast DYANA torsion angle dynamics algorithm, so that the result from each CANDID cycle consists of exhaustive, possibly ambiguous NOE cross-peak assignments in all available spectra and a three-dimensional protein structure represented by a bundle of conformers. The input for the first CANDID cycle consists of the amino acid sequence, the chemical shift list from the sequence-specific resonance assignment, and listings of the cross-peak positions and volumes in one or several two, three or four-dimensional NOESY spectra. The input for the second and subsequent CANDID cycles contains the three-dimensional protein structure from the previous cycle, in addition to the complete input used for the first cycle. CANDID includes two new elements that make it robust with respect to the presence of artifacts in the input data, i.e. network-anchoring and constraint-combination, which have a key role in de novo protein structure determinations for the successful generation of the correct polypeptide fold by the first CANDID cycle. Network-anchoring makes use of the fact that any network of correct NOE cross-peak assignments forms a self-consistent set; the initial, chemical shift-based assignments for each individual NOE cross-peak are therefore weighted by the extent to which they can be embedded into the network formed by all other NOE cross-peak assignments. Constraint-combination reduces the deleterious impact of artifact NOE upper distance constraints in the input for a protein structure calculation by combining the assignments for two or several peaks into a single upper limit distance constraint, which lowers the probability that the presence of an artifact peak will influence the outcome of the structure calculation. CANDID test calculations were performed with NMR data sets of four proteins for which high-quality structures had previously been solved by interactive protocols, and they yielded comparable results to these reference structure determinations with regard to both the residual constraint violations, and the precision and accuracy of the atomic coordinates. The CANDID approach has further been validated by de novo NMR structure determinations of four additional proteins. The experience gained in these calculations shows that once nearly complete sequence-specific resonance assignments are available, the automated CANDID approach results in greatly enhanced efficiency of the NOESY spectral analysis. The fact that the correct fold is obtained in cycle 1 of a de novo structure calculation is the single most important advance achieved with CANDID, when compared with previously proposed automated NOESY assignment methods that do not use network-anchoring and constraint-combination.

The program XEASY for computer-supported NMR spectral analysis of biological macromolecules.

SummaryA new program package, XEASY, was written for interactive computer support of the analysis of NMR spectra for three-dimensional structure determination of biological macromolecules. XEASY was developed for work with 2D, 3D and 4D NMR data sets. It includes all the functions performed by the precursor program EASY, which was designed for the analysis of 2D NMR spectra, i.e., peak picking and support of sequence-specific resonance assignments, cross-peak assignments, cross-peak integration and rate constant determination for dynamic processes. Since the program utilizes the X-window system and the Motif widget set, it is portable on a wide range of UNIX workstations. The design objective was to provide maximal computer support for the analysis of spectra, while providing the user with complete control over the final resonance assignments. Technically important features of XEASY are the use and flexible visual display of ‘strips’, i.e., two-dimensional spectral regions that contain the relevant parts of 3D or 4D NMR spectra, automated sorting routines to narrow down the selection of strips that need to be interactively considered in a particular assignment step, a protocol of resonance assignments that can be used for reliable bookkeeping, independent of the assignment strategy used, and capabilities for proper treatment of spectral folding and efficient transfer of resonance assignments between spectra of different types and different dimensionality, including projected, reduced-dimensionality triple-resonance experiments.

Automated NMR structure calculation with CYANA.

This chapter gives an introduction to automated nuclear magnetic resonance (NMR) structure calculation with the program CYANA. Given a sufficiently complete list of assigned chemical shifts and one or several lists of cross-peak positions and columns from two-, three-, or four-dimensional nuclear Overhauser effect spectroscopy (NOESY) spectra, the assignment of the NOESY cross-peaks and the three-dimensional structure of the protein in solution can be calculated automatically with CYANA.

Efficient computation of three-dimensional protein structures in solution from nuclear magnetic resonance data using the program DIANA and the supporting programs CALIBA, HABAS and GLOMSA

A novel procedure for efficient computation of three-dimensional protein structures from nuclear magnetic resonance (n.m.r.) data in solution is described, which is based on using the program DIANA in combination with the supporting programs CALIBA, HABAS and GLOMSA. The first part of this paper describes the new programs DIANA. CALIBA and GLOMSA. DIANA is a new, fully vectorized implementation of the variable target function algorithm for the computation of protein structures from n.m.r. data. Its main advantages, when compared to previously available programs using the variable target function algorithm, are a significant reduction of the computation time, and a novel treatment of experimental distance constraints involving diastereotopic groups of hydrogen atoms that were not individually assigned. CALIBA converts the measured nuclear Overhauser effects into upper distance limits and thus prepares the input for the previously described program HABAS and for DIANA. GLOMSA is used for obtaining individual assignments for pairs of diastereotopic substituents by comparison of the experimental constraints with preliminary results of the structure calculations. With its general outlay, the presently used combination of the four programs is particularly user-friendly. In the second part of the paper, initial results are presented on the influence of the novel DIANA treatment of diastereotopic protons on the quality of the structures obtained, and a systematic study of the central processing unit times needed for the same protein structure calculation on a range of different, commonly available computers is described.

Optimal isotope labelling for NMR protein structure determinations.

Nuclear-magnetic-resonance spectroscopy can determine the three-dimensional structure of proteins in solution. However, its potential has been limited by the difficulty of interpreting NMR spectra in the presence of broadened and overlapping resonance lines and low signal-to-noise ratios. Here we present stereo-array isotope labelling (SAIL), a technique that can overcome many of these problems by applying a complete stereospecific and regiospecific pattern of stable isotopes that is optimal with regard to the quality and information content of the resulting NMR spectra. SAIL uses exclusively chemically and enzymatically synthesized amino acids for cell-free protein expression. We demonstrate for the 17-kDa protein calmodulin and the 41-kDa maltodextrin-binding protein that SAIL offers sharpened lines, spectral simplification without loss of information, and the ability to rapidly collect the structural restraints required to solve a high-quality solution structure for proteins twice as large as commonly solved by NMR. It thus makes a large class of proteins newly accessible to detailed solution structure determination.

Improved efficiency of protein structure calculations from NMR data using the program DIANA with redundant dihedral angle constraints.

SummaryA new strategy for NMR structure calculations of proteins with the variable target function method (Braun, W. and Go, N. (1985)J. Mol. Biol.,186, 611) is described, which makes use of redundant dihedral angle constraints (REDAC) derived from preliminary calculations of the complete structure. The REDAC approach reduces the computation time for obtaining a group of acceptable conformers with the program DIANA 5-100-fold, depending on the complexity of the protein structure, and retains good sampling of conformation space.

RECOORD: A Recalculated Coordinate Database of 500 Proteins from the PDB Using Restraints from the BioMagResBank

State‐of‐the‐art methods based on CNS and CYANA were used to recalculate the nuclear magnetic resonance (NMR) solution structures of 500+ proteins for which coordinates and NMR restraints are available from the Protein Data Bank. Curated restraints were obtained from the BioMagResBank FRED database. Although the original NMR structures were determined by various methods, they all were recalculated by CNS and CYANA and refined subsequently by restrained molecular dynamics (CNS) in a hydrated environment. We present an extensive analysis of the results, in terms of various quality indicators generated by PROCHECK and WHAT_CHECK. On average, the quality indicators for packing and Ramachandran appearance moved one standard deviation closer to the mean of the reference database. The structural quality of the recalculated structures is discussed in relation to various parameters, including number of restraints per residue, NOE completeness and positional root mean square deviation (RMSD). Correlations between pairs of these quality indicators were generally low; for example, there is a weak correlation between the number of restraints per residue and the Ramachandran appearance according to WHAT_CHECK (r = 0.31). The set of recalculated coordinates constitutes a unified database of protein structures in which potential user‐ and software‐dependent biases have been kept as small as possible. The database can be used by the structural biology community for further development of calculation protocols, validation tools, structure‐based statistical approaches and modeling. The RECOORD database of recalculated structures is publicly available from http://www.ebi.ac.uk/msd/recoord. Proteins 2005.

Protein structure determination in living cells by in-cell NMR spectroscopy

Investigating proteins ‘at work’ in a living environment at atomic resolution is a major goal of molecular biology, which has not been achieved even though methods for the three-dimensional (3D) structure determination of purified proteins in single crystals or in solution are widely used. Recent developments in NMR hardware and methodology have enabled the measurement of high-resolution heteronuclear multi-dimensional NMR spectra of macromolecules in living cells (in-cell NMR). Various intracellular events such as conformational changes, dynamics and binding events have been investigated by this method. However, the low sensitivity and the short lifetime of the samples have so far prevented the acquisition of sufficient structural information to determine protein structures by in-cell NMR. Here we show the first, to our knowledge, 3D protein structure calculated exclusively on the basis of information obtained in living cells. The structure of the putative heavy-metal binding protein TTHA1718 from Thermus thermophilus HB8 overexpressed in Escherichia coli cells was solved by in-cell NMR. Rapid measurement of the 3D NMR spectra by nonlinear sampling of the indirectly acquired dimensions was used to overcome problems caused by the instability and low sensitivity of living E. coli samples. Almost all of the expected backbone NMR resonances and most of the side-chain NMR resonances were observed and assigned, enabling high quality (0.96 ångström backbone root mean squared deviation) structures to be calculated that are very similar to the in vitro structure of TTHA1718 determined independently. The in-cell NMR approach can thus provide accurate high-resolution structures of proteins in living environments.

Processing of multi dimensional NMR data with the new software PROSA

SummaryThe new program PROSA is an efficient implementation of the common data-processing steps for multi-dimensional NMR spectra. PROSA performs linear prediction, digital filtering, Fourier transformation, automatic phase correction, and baseline correction. High efficiency is achieved by avoiding disk storage of intermediate data and by the absence of any graphics display, which enables calculation in the batch mode and facilitates porting PROSA on a variety of different computer systems; including supercomputers. Furthermore, all time-consuming routines are completely vectorized. The elimination of a graphics display was made possible by the use of a new, reliable automatic phase-correction routine. CPU times for complete processing of a typical heteronuclear three-dimensional NMR data set of a protein vary between less than 1 min on a NEC SX3 supercomputer and 40 min on a Sun-4 computer system.

NMR structure reveals intramolecular regulation mechanism for pheromone binding and release

Odorants are transmitted by small hydrophobic molecules that cross the aqueous sensillar lymph surrounding the dendrites of the olfactory neurons to stimulate the olfactory receptors. In insects, the transport of pheromones, which are a special class of odorants, is mediated by pheromone-binding proteins (PBPs), which occur at high concentrations in the sensillar lymph. The PBP from the silk moth Bombyx mori (BmPBP) undergoes a pH-dependent conformational transition between the forms BmPBPA present at pH 4.5 and BmPBPB present at pH 6.5. Here, we describe the NMR structure of BmPBPA, which consists of a tightly packed arrangement of seven α-helices linked by well defined peptide segments and knitted together by three disulfide bridges. A scaffold of four α-helices that forms the ligand binding site in the crystal structure of a BmPBP–pheromone complex is preserved in BmPBPA. The C-terminal dodecapeptide segment, which is in an extended conformation and located on the protein surface in the pheromone complex, forms a regular helix, α7, which is located in the pheromone-binding site in the core of the unliganded BmPBPA. Because investigations by others indicate that the pH value near the membrane surface is reduced with respect to the bulk sensillar lymph, the pH-dependent conformational transition of BmPBP suggests a novel physiological mechanism of intramolecular regulation of protein function, with the formation of α7 triggering the release of the pheromone from BmPBP to the membrane-standing receptor.