Sparsely-sampled,high-resolution 4-D omit spectra for detection and assignment of intermolecular NOEs of protein complexes

Unambiguous detection and assignment of intermolecular NOEs are essential for structure determination of protein complexes by NMR. Such information has traditionally been obtained with 3-D half-filtered experiments, where scalar coupling-based purging of intramolecular signals allows for selective detection of intermolecular NOEs. However, due to the large variation of ¹J_HC scalar couplings and limited chemical shift dispersion in the indirect proton dimension, it is difficult to obtain reliable and complete assignments of interfacial NOEs. Here, we demonstrate a strategy that combines selective labeling and high-resolution 4-D NOE spectroscopy with sparse sampling for reliable identification and assignment of intermolecular NOEs. Spectral subtraction of component-labeled complexes from a uniformly-labeled protein complex yields an “omit” spectrum containing positive intermolecular NOEs with little signal degeneracy. Such a strategy can be broadly applied to unbiased detection, assignment and presentation of intermolecular NOEs of protein complexes.     

Enhancement of bound-state residual dipolar couplings: conformational analysis of lactose bound to Galectin-3

Residual dipolar couplings (RDCs) have proven to be a valuable NMR tool that can provide long-range constraints for molecular structure determination. The constraints are orientational in nature and are, thus, highly complementary to conventional distance constraints from NOE data. This complementarity would seem to extend to the study of the geometry of ligands bound to proteins. However, unlike transferred NOEs, where collection, even with a large excess of free ligand, results in measurements dominated by bound contributions, RDCs of exchanging ligands can be dominated by free-state contributions. Here we present a strategy for enhancement of RDCs from bound states that is based on specifically enhancing the alignment of the protein to which a ligand will bind. The protein is modified by addition of a hydrophobic alkyl tail that anchors it to the bicelles that are a part of the ordering medium needed for RDC measurement. As an illustration, we have added a propyl chain to the C terminus of the carbohydrate recognition domain of the protein, Galectin-3, and report enhanced RDCs that prove consistent with known bound-ligand geometries for this protein.     

Structure determination of a Galectin-3-carbohydrate complex using paramagnetism-based NMR constraints

The determination of the location and conformation of a natural ligand bound to a protein receptor is often a first step in the rational design of molecules that can modulate receptor function. NMR observables, including NOEs, often provide the basis for these determinations. However, when ligands are carbohydrates, interactions mediated by extensive hydrogen-bonding networks often reduce or eliminate NOEs between ligand and protein protons. In these cases, it is useful to look to other distance- and orientation-dependent observables that can constrain the geometry of ligand-protein complexes. Here we illustrate the use of paramagnetism-based NMR constraints, including pseudo-contact shifts (PCS) and field-induced residual dipolar couplings (RDCs). When a paramagnetic center can be attached to the protein, field-induced RDCs and PCS reflect only bound-state properties of the ligand, even when averages over small fractions of bound states and large fractions of free states are observed. The effects can also be observed over a long range, making it possible to attach a paramagnetic center to a remote part of the protein. The system studied here is a Galectin-3-lactose complex. A lanthanide-binding peptide showing minimal flexibility with respect to the protein was integrated into the C terminus of an expression construct for the Galectin-3-carbohydrate-binding domain. Dysprosium ion, which has a large magnetic susceptibility anisotropy, was complexed to the peptide, making it possible to observe both PCSs and field-induced RDCs for the protein and the ligand. The structure determined from these constraints shows agreement with a crystal structure of a Galectin-3-N-acetyllactosamine complex.     

Generation of native-like protein structures from limited NMR data,modern force fields and advanced conformational sampling

Determining an accurate initial native-like protein fold is one of the most important and time-consuming steps of de novo NMR structure determination. Here we demonstrate that high-quality native-like models can be rapidly generated from initial structures obtained using limited NOE assignments, through replica exchange molecular dynamics refinement with a generalized Born implicit solvent (REX/GB). Conventional structure calculations using an initial sparse NOE set were unable to identify a unique topology for the zinc-bound C-terminal domain of E. coli chaperone Hsp33, due to a lack of unambiguous long range NOEs. An accurate overall topology was eventually obtained through laborious hand identification of long range NOEs. However we were able to obtain high-quality models with backbone RMSD values of about 2 Å with respect to the final structures, using REX/GB refinement with the original limited set of initial NOE restraints. These models could then be used to make further assignments of ambiguous NOEs and thereby speed up the structure determination process. The ability to calculate accurate starting structures from the limited unambiguous NOE set available at the beginning of a structure calculation offers the potential of a much more rapid and automated process for NMR structure determination. Jianhan Chen: Authors contributed equally to this work.Hyung-Sik Won: Authors contributed equally to this work.     

A complete algorithm to resolve ambiguity for intersubunit NOE assignment in structure determination of symmetric homo-oligomers

Assignment of nuclear Overhauser effect (NOE) data is a key bottleneck in structure determination by NMR. NOE assignment resolves the ambiguity as to which pair of protons generated the observed NOE peaks, and thus should be restrained in structure determination. In the case of intersubunit NOEs in symmetric homo-oligomers, the ambiguity includes both the identities of the protons within a subunit, and the identities of the subunits to which they belong. This paper develops an algorithm for simultaneous intersubunit NOE assignment and C(n) symmetric homo-oligomeric structure determinations, given the subunit structure. By using a configuration space framework, our algorithm guarantees completeness, in that it identifies structures representing, to within a user-defined similarity level, every structure consistent with the available data (ambiguous or not). However, while our approach is complete in considering all conformations and assignments, it avoids explicit enumeration of the exponential number of combinations of possible assignments. Our algorithm can draw two types of conclusions not possible under previous methods: (1) that different assignments for an NOE would lead to different structural classes, or (2) that it is not necessary to uniquely assign an NOE, since it would have little impact on structural precision. We demonstrate on two test proteins that our method reduces the average number of possible assignments per NOE by a factor of 2.6 for MinE and 4.2 for CCMP. It results in high structural precision, reducing the average variance in atomic positions by factors of 1.5 and 3.6, respectively.     

Calculation of symmetric multimer structures from NMR data using a priori knowledge of the monomer structure, co-monomer restraints, and interface mapping: The case of leucine zippers

Summary NMR studies of symmetric multimers are problematic due to the difficulty in distinguishing between intra-, inter-, and co-monomer (mixed) NOE signals. Previously, one of us described a general calculation strategy called dynamic assignment by which this difficulty can be overcome [Nilges, M. (1993) Proteins, 17, 297–309]. Here we describe extensions to the method for handling many co-monomer NOEs and for taking advantege of prior knowledge of monomer structures. The new protocol was developed for the particularly difficult case of leucine zipper (LZ) homodimers, for which the previous protocol proved inefficient. In addition to the problem of dimer symmetry, LZs have a particularly high proportion of co-monomer NOE signals and a high degree of repetition in sequence and structure, leading to significant spectral overlap. Furthermore, the leucine zipper is a rather extended (as opposed to globular) protein domain; accurately determining such a structure based only on the very short distances obtainable by NMR is clearly a challenge to the NMR structure determination method. We have previously shown that, for LZ homodimers, many of the backbone-backbone NOESY cross peaks can be unambiguously assigned as intra-monomer, enabling approximate monomer structures to be calculated. Using model and experimental data sets, we verified that the new protocol converges to the correct dimer structure. The results show that short-range NMR distance data can be sufficient to define accurately the extended LZ. The protocol has been used to derive a novel solution structure of the c-Jun LZ domain. Based on these calculations, we propose the protocol as a prototype for the general case of symmetric multimers where the monomer structure is known.Abbreviations 3D three-dimensional - GCN4-c crystal structure of the GCN4 LZ homodimer - GCN4-s solution structures of GCN4 - GSYM global symmetry - Jun-m model structure of the Jun LZ homodimer - Jun-s solution structure of Jun - LZ leucine zipper - MFP mean force potential - MDSA molecular dynamical simulated annealing - NCS noncrystallographic symmetry - NOE nuclear Overhauser enhancement - rmsd root-mean-square deviation - vdW van der Waals     

RDC-assisted modeling of symmetric protein homo-oligomers

Protein oligomerization serves an important function in biological processes, yet solving structures of protein oligomers has always been a challenge. For solution NMR, the challenge arises both from the increased size of these systems and, in the case of homo-oligomers, from ambiguities in assignment of intra- as opposed to intersubunit NOEs. In this study, we present a residual dipolar coupling (RDC)-assisted method for constructing models of homo-oligomers with purely rotational symmetry. Utilizing the fact that one of the principal axes of the tensor describing the alignment needed for RDC measurement is always parallel to the oligomer symmetry axis, it is possible to greatly restrict possible models for the oligomer. Here, it is shown that, if the monomer structure is known, all allowed dimer models can be constructed using a grid search algorithm and evaluated based on RDC simulations and the quality of the interface between the subunits. Using the Bacillus subtilis protein YkuJ as an example, it is shown that the evaluation criteria based on just two sets of NH RDCs are very selective and can unambiguously produce a model in good agreement with an existing X-ray structure of YkuJ.     

Relaxation data in NMR structure determination: model calculations for the lysozyme-Gd3+ complex

The effect of including paramagnetic relaxation data as additional restraints in the determination of protein tertiary structures from NMR data has been explored by a systematic series of model calculations. The system used for testing the method was the 2.0 A resolution tetragonal crystal structure of hen egg white lysozyme (129 amino acid residues) and structures were generated using a version of the hybrid "distance geometry-dynamic simulated annealing" procedure. A limited set of 769 NOEs was used as restraints in all the calculations; the strengths of these were categorized into three classes on the basis of distances observed in the crystal structure. The values of 50 phi angles were also restrained on the basis of amide-alpha coupling constants calculated from the X-ray structure. Five sets of 12 structures were determined using differing sets of paramagnetic relaxation data as restraints additional to those involving the NOE and coupling constant data. The paramagnetic relaxation data were modeled on the basis of the distances of defined protons from the crystallographic binding site of Gd3+ in lysozyme. Analysis of the results showed that the relaxation data significantly improved the correspondence between the set of generated structures and the crystal structure, and that the more well defined the relaxation data, the more significant the improvement in the quality of the structures. The results suggest that the inclusion of paramagnetic relaxation restraints could be of significant value for the experimental determination of protein structures from NMR data.     

Protein NMR structure determination with automated NOE-identification in the NOESY spectra using the new software ATNOS

Novel algorithms are presented for automated NOESY peak picking and NOE signal identification in homonuclear 2D and heteronuclear-resolved 3D [¹H,¹H]-NOESY spectra during de novoprotein structure determination by NMR, which have been implemented in the new software ATNOS (automated NOESY peak picking). The input for ATNOS consists of the amino acid sequence of the protein, chemical shift lists from the sequence-specific resonance assignment, and one or several 2D or 3D NOESY spectra. In the present implementation, ATNOS performs multiple cycles of NOE peak identification in concert with automated NOE assignment with the software CANDID and protein structure calculation with the program DYANA. In the second and subsequent cycles, the intermediate protein structures are used as an additional guide for the interpretation of the NOESY spectra. By incorporating the analysis of the raw NMR data into the process of automated de novoprotein NMR structure determination, ATNOS enables direct feedback between the protein structure, the NOE assignments and the experimental NOESY spectra. The main elements of the algorithms for NOESY spectral analysis are techniques for local baseline correction and evaluation of local noise level amplitudes, automated determination of spectrum-specific threshold parameters, the use of symmetry relations, and the inclusion of the chemical shift information and the intermediate protein structures in the process of distinguishing between NOE peaks and artifacts. The ATNOS procedure has been validated with experimental NMR data sets of three proteins, for which high-quality NMR structures had previously been obtained by interactive interpretation of the NOESY spectra. The ATNOS-based structures coincide closely with those obtained with interactive peak picking. Overall, we present the algorithms used in this paper as a further important step towards objective and efficient de novoprotein structure determination by NMR.     

A geometric arrangement algorithm for structure determination of symmetric protein homo-oligomers from NOEs and RDCs

Nuclear magnetic resonance (NMR) spectroscopy is a primary tool to perform structural studies of proteins in physiologically-relevant solution conditions. Restraints on distances between pairs of nuclei in the protein, derived from the nuclear Overhauser effect (NOE), provide information about the structure of the protein in its folded state. NMR studies of symmetric protein homo-oligomers present a unique challenge. Using X-filtered NOESY experiments, it is possible to determine whether an NOE restrains a pair of protons across different subunits or within a single subunit, but current experimental techniques are unable to determine in which subunits the restrained protons lie. Consequently, it is difficult to assign NOEs to particular pairs of subunits with certainty, thus hindering the structural analysis of the oligomeric state. Computational approaches are needed to address this subunit ambiguity, but traditional solutions often rely on stochastic search coupled with simulated annealing and simulations of simplified molecular dynamics, which have many tunable parameters that must be chosen carefully and can also fail to report structures consistent with the experimental restraints. In addition, these traditional approaches rarely provide guarantees on running time or solution quality. We reduce the structure determination of homo-oligomers with cyclic symmetry to computing geometric arrangements of unions of annuli in a plane. Our algorithm, disco, runs in expected O(n2) time, where n is the number of distance restraints, potentially assigned ambiguously. disco is guaranteed to report the exact set of oligomer structures consistent with the distance restraints and also with orientational restraints from residual dipolar couplings (RDCs). We demonstrate our method using two symmetric protein complexes: the trimeric E. coli diacylglycerol kinase (DAGK) and a dimeric mutant of the immunoglobulin-binding domain B1 of streptococcal protein G (GB1). In both cases, disco computes oligomer structures with high precision and also finds distance restraints that are either mutually inconsistent or inconsistent with the RDCs. The entire protocol DISCO has been completely automated in a software package that is freely available and open-source at www.cs.duke.edu/donaldlab/software.php.     

Fast methionine-based solution structure determination of calcium-calmodulin complexes

Here we present a novel NMR method for the structure determination of calcium-calmodulin (Ca²⁺-CaM)-peptide complexes from a limited set of experimental restraints. A comparison of solved CaM-peptide structures reveals invariability in CaM’s backbone conformation and a structural plasticity in CaM’s domain orientation enabled by a flexible linker. Knowing this, the collection and analysis of an extensive set of NOESY spectra is redundant. Although RDCs can define CaM domain orientation in the complex, they lack the translational information required to position the domains on the bound peptide and highlight the necessity of intermolecular NOEs. Here we employ a specific isotope labeling strategy in which the role of methionine in CaM-peptide interactions is exploited to collect these critical NOEs. By ¹H, ¹³C-labeling the methyl groups of deuterated methionine against a ²H, ¹²C background, we can acquire a ¹³C-edited NOESY characterized by simplified, easily analyzable spectra. Together with measured CaM backbone H_N-N RDCs and intrapeptide NOE-based distances, these intermolecular NOEs provide restraints for a low temperature torsion-angle dynamics and simulated annealing protocol used to calculate the complex structure. We have applied our method to a CaM complex previously solved through X-ray crystallography: Ca²⁺-CaM bound to the CaM kinase I peptide (PDB code: 1MXE). The resulting structure has a backbone RMSD of 1.6 Å to that previously published. We have also used this test complex to investigate the importance of homologous model selection on the calculated outcome. In addition to having application for fast complex structure determination, this method can be used to determine the structures of difficult complexes characterized by chemical shift overlap and broad signals for which the traditional method based on the use of fully ¹³C, ¹⁵N-labeled CaM fails.     

The EAGR box structure: a motif involved in mycoplasma motility

Mycoplasma genitalium is an emerging human pathogen with the smallest genome found among self‐replicating organisms. M. genitalium presents a complex cytoskeleton with a differentiated protrusion known as the terminal organelle. This polar structure plays a central role in functions essential for the virulence of the microorganism, such as motility and cell‐host adhesion. A well‐conserved Enriched in Aromatic and Glycine Residues motif, the EAGR box, is present in many of the proteins found in the terminal organelle. We determined the crystal structure of the globular domain from M. genitalium MG200 that contains an EAGR box. This structural information is the first at near atomic resolution for the components of the terminal organelle. The structure revealed a dimer stabilized by a compact hydrophobic core that extends throughout the dimer interface. Monomers present a new fold that contains an accurate intra‐subunit symmetry relating two conspicuous hairpins. Some features, such as the domain plasticity and the presence and organization of the intra‐ and inter‐subunit symmetry axes, support a role for the EAGR box in protein–protein interactions. Genetic, biochemical and microcinematography analyses of MG200 variants lacking the EAGR box containing domain confirm the relevant and specific association of this domain with cell motility.     

Dynamic features of homodimer interfaces calculated by normal‐mode analysis

Knowledge of the dynamic features of protein interfaces is necessary for a deeper understanding of protein–protein interactions. We performed normal‐mode analysis (NMA) of 517 nonredundant homodimers and their protomers to characterize dimer interfaces from a dynamic perspective. The motion vector calculated by NMA for each atom of a dimer was decomposed into internal and external motion vectors in individual component subunits, followed by the averaging of time‐averaged correlations between these vectors over atom pairs in the interface. This averaged correlation coefficient (ACC) was defined for various combinations of vectors and investigated in detail. ACCs decrease exponentially with an increasing interface area and r‐value, that is, interface area divided by the entire subunit surface area. As the r‐value reflects the nature of dimer formation, the result suggests that both the interface area and the nature of dimer formation are responsible for the dynamic properties of dimer interfaces. For interfaces with small or medium r‐values and without intersubunit entanglements, ACCs are found to increase on dimer formation when compared with those in the protomer state. In contrast, ACCs do not increase on dimer formation for interfaces with large r‐values and intersubunit entanglements such as in interwinding dimers. Furthermore, relationships between ACCs for intrasubunit atom pairs and for intersubunit atom pairs are found to significantly differ between interwinding and noninterwinding dimers for external motions. External motions are considered as an important factor for characterizing dimer interfaces.     

Active and inactive state structures of unliganded Lactobacillus casei allosteric L‐lactate dehydrogenase

Lactobacillus casei L ‐lactate dehydrogenase (LCLDH) is activated through the homotropic and heterotropic activation effects of pyruvate and fructose 1,6‐bisphosphate (FBP), respectively, and exhibits unusually high pH‐dependence in the allosteric effects of these ligands. The active (R) and inactive (T) state structures of unliganded LCLDH were determined at 2.5 and 2.6 Å resolution, respectively. In the catalytic site, the structural rearrangements are concerned mostly in switching of the orientation of Arg171 through the flexible intersubunit contact at the Q‐axis subunit interface. The distorted orientation of Arg171 in the T state is stabilized by a unique intra‐helix salt bridge between Arg171 and Glu178, which is in striking contrast to the multiple intersubunit salt bridges in Lactobacillus pentosus nonallosteric L ‐lactate dehydrogenase. In the backbone structure, major structural rearrangements of LCLDH are focused in two mobile regions of the catalytic domain. The two regions form an intersubunit linkage through contact at the P‐axis subunit interface involving Arg185, replacement of which with Gln severely decreases the homotropic and hetertropic activation effects on the enzyme. These two regions form another intersubunit linkage in the Q‐axis related dimer through the rigid NAD‐binding domain, and thus constitute a pivotal frame of the intersubunit linkage for the allosteric motion, which is coupled with the concerted structural change of the four subunits in a tetramer, and of the binding sites for pyruvate and FBP. The unique intersubunit salt bridges, which are observed only in the R state structure, are likely involved in the pH‐dependent allosteric equilibrium. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Sequence-specific 1H NMR assignments and secondary structure in solution of Escherichia coli trp repressor

Sequence-specific 1H NMR assignments are reported for the active L-tryptophan-bound form of Escherichia coli trp repressor. The repressor is a symmetric dimer of 107 residues per monomer; thus at 25 kDa, this is the largest protein for which such detailed sequence-specific assignments have been made. At this molecular mass the broad line widths of the NMR resonances preclude the use of assignment methods based on 1H-1H scalar coupling. Our assignment strategy centers on two-dimensional nuclear Overhauser spectroscopy (NOESY) of a series of selectively deuterated repressor analogues. A new methodology was developed for analysis of the spectra on the basis of the effects of selective deuteration on cross-peak intensities in the NOESY spectra. A total of 90% of the backbone amide protons have been assigned, and 70% of the alpha and side-chain proton resonances are assigned. The local secondary structure was calculated from sequential and medium-range backbone NOEs with the double-iterated Kalman filter method [Altman, R. B., & Jardetzky, O. (1989) Methods Enzymol. 177, 218-246]. The secondary structure agrees with that of the crystal structure [Schevitz, R., Otwinowski, Z., Joachimiak, A., Lawson, C. L., & Sigler, P. B. (1985) Nature 317, 782], except that the solution state is somewhat more disordered in the DNA binding region and in the N-terminal region of the first alpha-helix. Since the repressor is a symmetric dimer, long-range intersubunit NOEs were distinguished from intrasubunit interactions by formation of heterodimers between two appropriate selectively deuterated proteins and comparison of the resulting NOESY spectrum with that of each selectively deuterated homodimer. Thus, from spectra of three heterodimers, long-range NOEs between eight pairs of residues were identified as intersubunit NOEs, and two additional long-range intrasubunits NOEs were assigned.     

Determination of the NMR solution structure of a specific DNA complex of the Myb DNA-binding domain

Summary The solution structure of a specific DNA complex of the minimum DNA-binding domain of the mouse c-Myb protein was determined by distance geometry calculations using a set of 1732 nuclear Overhauser enhancement (NOE) distance restraints. In order to determine the complex structure independent of the initial guess, we have developed two different procedures for the docking calculation using simulated annealing in four-dimensional space (4D-SA). One is a multiple-step procedure, where the protein and the DNA were first constructed independently by 4D-SA using only the individual intramolecular NOE distance restraints. Here, the initial structure of the protein was a random coil and that of the DNA was a typical B-form duplex. Then, as the starting structure for the next docking procedure, the converged protein and DNA structures were placed in random molecular orientations, separated by 50 Å. The two molecules were docked by 4D-SA utilizing all the restraints, including the additional 66 intermolecular distance restraints. The second procedure comprised a single step, in which a random-coil protein and a typical B-form DNA duplex were first placed 70 Å from each other. Then, using all the intramolecular and intermolecular NOE distance restraints, the complex structure was constructed by 4D-SA. Both procedures yielded the converged complex structures with similar quality and structural divergence, but the multiple-step procedure has much better convergence power than the single-step procedure. A model study of the two procedures was performed to confirm the structural quality, depending upon the number of intermolecular distance restraints, using the X-ray structure of the engrailed homeodomain-DNA complex.Abbreviations rmsd root-mean-square deviation - NOE nuclear Overhauser enhancement - 4D-SA simulated annealing in four-dimensional space - Myb-R2R3 repeats 2 and 3 of the DNA-binding domain of the c-Myb protein - DNA 16 Myb-specific binding DNA duplex with 16 base pairs - IHDD-C residues 3 to 59 of the C-chain of the engrailed homeodomain-DNA complex - DNA11 DNA duplex with base pairs 9 to 19 of the engrailed homeodomain-DNA complex     

Determination of solution structures of paramagnetic proteins by NMR

Standard procedures for using nuclear Overhauser enhancements (NOE) between protons to generate structures for diamagnetic proteins in solution from NMR data may be supplemented by using dipolar shifts if the protein is paramagnetic. This is advantageous since the electron-nuclear dipolar coupling provides relatively long-range geometric information with respect to the paramagnetic centre which complements the short-range distance constraints from NOEs. Several different strategies have been developed to date, but none of these attempts to combine data from NOEs and dipolar shifts in the initial stages of structure calculation or to determine three dimensional protein structures together with their magnetic properties. This work shows that the magnetic and atomic structures are highly correlated and that it is important to have additional constraints both to provide starting parameters for the magnetic properties and to improve the definition of the best fit. Useful parameters can be obtained for haem proteins from Fermi contact shifts; this approach is compared with a new method based on the analysis of dipolar shifts in haem methyl groups with respect to data from horse and tuna ferricytochromes c. The methods developed for using data from NOEs and dipolar shifts have been incorporated in a new computer program, PARADYANA, which is demonstrated in application to a model data set for the sequence of the haem octapeptide known as microperoxidase-8. Received: 13 October 1997 / Accepted: 19 December 1997     

NMR Assignment of the mTOR Domain Responsible for Rapamycin Binding

A fully automated, NOE-based NMR structure determination of a uniformly ¹³C,¹⁵N-labeled protein was achieved in crude cell-extract, without purification of the overexpressed protein. Essentially complete sequence-specific assignments were obtained using triple resonance experiments, based on the high intensity of the resonances from the overexpressed protein relative to those of the background. For the collection of NOE distance constraints, efficient discrimination between NOE cross peaks from the target protein and background signals was achieved using the programs ATNOS and CANDID. In the iterative ATNOS/CANDID procedure, the identification of the desired protein NOEs is initially guided by the self-consistency of the protein NOE-network. Although the intensities of the signals in this network vary over a wide range, and are in many instances comparable to or smaller than those of the background, the first cycle of calculations resulted in the correct global polypeptide fold, and the structure was then refined in six subsequent cycles using the intermediate NMR structures for additional guidance. The experience gained with this work demonstrates that the ATNOS/CANDID procedure for automatic protein structure determination is highly robust and reliable in the presence of intense background signals, and might thus also represent a platform for future protein structure determinations in physiological fluids.     

1H NMR assignments of apo calcyclin and comparative structural analysis with calbindin D9k and S100 beta.

The homodimeric S100 protein calcyclin has been studied in the apo state by two-dimensional 1H NMR spectroscopy. Using a combination of scalar correlation and NOE experiments, sequence-specific 1H NMR assignments were obtained for all but one backbone and > 90% of the side-chain resonances. To our knowledge, the 2 x 90 residue (20 kDa) calcyclin dimer is the largest protein system for which such complete assignments have been made by purely homonuclear methods. Sequential and medium-range NOEs and slowly exchanging backbone amide protons identified directly the four helices and the short antiparallel beta-type interaction between the two binding loops that comprise each subunit of the dimer. Further analysis of NOEs enabled the unambiguous assignment of 556 intrasubunit distance constraints, 24 intrasubunit hydrogen bonding constraints, and 2 x 26 intersubunit distance constraints. The conformation of the monomer subunit was refined by distance geometry and restrained molecular dynamics calculations using the intrasubunit constraints only. Calculation of the dimer structure starting from this conformational ensemble has been reported elsewhere. The extent of structural homology among the apo calcyclin subunit, the monomer subunit of apo S100 beta, and monomeric apo calbindin D9k has been examined in detail by comparing 1H NMR chemical shifts and secondary structures. This analysis was extended to a comprehensive comparison of the three-dimensional structures of the calcyclin monomer subunit and calbindin D9k, which revealed greater similarity in the packing of their hydrophobic cores than was anticipated previously. Together, these results support the hypothesis that all members of the S100 family have similar core structures and similar modes of dimerization. Analysis of the amphiphilicity of Helix IV is used to explain why calbindin D9k is monomeric, but full-length S100 proteins form homodimers.