NOE Restraints Research Articles

Solution structures of the R6 human insulin hexamer,.

The three-dimensional solution structure of the phenol-stabilized 36 kDa R6 insulin hexamer was determined by NMR spectroscopy and restrained molecular dynamics. The hexamer structures were derived using a stepwise procedure. Initially, 60 monomers were obtained by distance geometry from 665 NOE-derived distance restraints and three disulfide bridges. Subsequently, the hexamer structures were calculated by simulated annealing, using 30 hexamers constructed from the best 36 monomer structures as the starting models. The NMR data show that the aromatic ring of residue Phe(B25) can take two different orientations in the solution hexamer: one in which it points inward (molecule 1, about 90%) and one in which it points outward from the surface of the monomer (molecule 2, about 10%). Therefore, two hexamer structures were calculated: a symmetric hexamer consisting of six molecule 1 monomers and a nonsymmetric hexamer consisting of five molecule 1 monomers and one molecule 2 monomer. For each of the six monomers, the restraints used in the calculations of the hexamer structures include, in addition to the intramonomeric restraints, 25 NOEs between insulin and phenol, 23 NOEs and two hydrogen bonds across the dimer interface, nine NOEs across the trimer interface, and five intramonomeric or two intermonomeric NOEs, respectively, specifying the different orientations of the Phe(B25) ring. The coordination of the two Zn atoms was defined by eight distance restraints. Thus, a total of 4394 and 4391 distance restraints, respectively, were used in the two hexamer calculations. The NOE restraints were classified in an iterative process as intra- or intermonomeric on the basis of their consistency or inconsistency with the structure of the monomer. The assignment of the dimer- and trimer-specific NOEs was made using the crystal structure of the R6 hexamer as the starting model. For both solution hexamers, the average backbone rms deviation is 0.81 A, if the less well-defined N- and C-terminal residues are excluded. The corresponding rms deviations for all heavy atoms are 1.17 and 1.19 A for the nonsymmetric and symmetric hexamer, respectively. The overall solution structure of the R6 insulin hexamer is compact, rigid, and symmetric and resembles the corresponding crystal structure. However, the extension of the B-chain alpha-helix, which characterizes the R state, is shorter in the solution structure than in the crystal structure. Also, the study shows that the orientation of the Phe(B25) ring has no effect on the structure of the rest of the molecule, within the uncertainty of the structure determination. The importance of these findings for the current model for the insulin-receptor interaction is discussed.

NMR solution structure of the antitumor compound PT523 and NADPH in the ternary complex with human dihydrofolate reductase.

The antitumor compound PT523 [N(alpha)-(4-amino-4-deoxypteroyl)-N(delta)-hemiphthaloyl-L- ornithine] was found to have an inhibition constant (K(i)) of 0.35 +/- 0.10 pM against human dihydrofolate reductase (hDHFR), 15-fold lower than that of the classical antifolate drug methotrexate (MTX). The structure of PT523 bound to hDHFR and hDHFR-NADPH was investigated using multinuclear NMR techniques. NMR data indicate that the binary complex has two distinct conformations in solution which are in slow exchange and that the addition of NADPH stabilizes the ternary complex in a single bound state. Comparison of resonance assignments in the PT523 and MTX ternary complexes revealed that substantial protein chemical shift differences are limited to small regions of hDHFR tertiary structure. A restrained molecular dynamics and energy minimization protocol was performed for the hDHFR-PT523-NADPH complex, using 185 NOE restraints (33 intermolecular) to define the ligand-binding region. The positions of the pteridine and pABA rings of PT523 and the nicotinamide and ribose rings of NADPH are well defined in the solution structures (RMSD = 0.59 A) and are consistent with previously determined structures of DHFR complexes. The N(delta)-hemiphthaloyl-L-ornithine group of PT523 is less well defined, and the calculated model structures suggest the hemiphthaloyl ring may adopt more than one conformation in solution. Contacts between the hemiphthaloyl ring and hDHFR, which are not possible in the hDHFR-MTX-NADPH complex, may explain the greater inhibition potency of PT523.

NMR structure determination of a novel conotoxin, [Pro 7,13] alpha A-conotoxin PIVA.

A high-resolution solution conformation of a novel conotoxin, [Pro 7,13] alpha A-conotoxin PIVA, GCCGSYPNAACHPCSCKDROSYCGQ-NH2, has been determined by two-dimensional 1H NMR methods and distance geometry calculations. The total of 324 NOE-derived interproton distance restraints including 33 long-range NOE restraints as well as 11 phi and 7 chi 1 torsion angle restraints was used for computation of structures. Back-calculation from the experimental NOE spectrum has provided 49 new NOE restraints and yielded the final R-factors of Ra = 0.641 and Rb = 0.157. The final RMSD values are 0.90 and 1.16 A for the backbone and the heavy atoms, respectively. The C-terminal half of the molecule involving the residues 12-24 is extremely well-defined with a backbone RMSD value of 0.56 A, whereas the N-terminal 3-11 disulfide loop is relatively flexible, possessing a backbone RMSD value of 1.09 A. The [Pro 7,13] alpha A-conotoxin PIVA does not contain any significant secondary structure although the 21S-24G nearly completes one turn of a 3(10) helix. The overall protein fold is largely maintained by the three disulfide bridges of 2-16, 3-11, and 14-23. The presence of the three disulfide bridges imposes geometric constraints that force the molecule to form six continuous bends involving the following residues: 3C-5S, 7P-10A, 12H-14C, 15S-17K, 17K-19R, and 21S-25Q. The overall shape of the [Pro 7,13] alpha A-conotoxin PIVA can be described as an "iron". Residues 15S-19R form a loop that protrudes out of the "bottom plate" formed by the rest of the protein and constitute the handle of the iron. The N-terminal tip of the molecule is relatively immobile due to attractive electrostatic interactions between the gamma-hydroxyl group of 20 Hyp and the phenolic hydroxyl group of 22Y. The flexible 3-11 disulfide loop consists mostly of hydrophobic residues, while the best-defined 14-23 disulfide loop contains the highly charged hydrophilic 15S-19R "handle" domain exposed to the exterior of the protein. Binding to nicotinic acetylcholine receptor can be mediated through two different types of interactions: one involving the aromatic hydrophobic residues such as 6Y and 12H and the other involving the positively charged hydrophilic side chain of the 19R. The side chain of the 19R in the [Pro 7, 13] alpha A-conotoxin PIVA and that of the 9R of the alpha-conotoxin G1, and also the side chains of the 12H and 6Y in the former and those of 10H and 11Y in the latter can be aligned to point to the same direction when the corresponding backbone atoms are superimposed to an RMSD value of 2.5 A.

An approach to global fold determination using limited NMR data from larger proteins selectively protonated at specific residue types

A combination of calculation and experiment is used to demonstrate that the global fold of larger proteins can be rapidly determined using limited NMR data. The approach involves a combination of heteronuclear triple resonance NMR experiments with protonation of selected residue types in an otherwise completely deuterated protein. This method of labelling produces proteins with alpha-specific deuteration in the protonated residues, and the results suggest that this will improve the sensitivity of experiments involving correlation of side-chain ((1)H and (13)C) and backbone ((1)H and (15)N) amide resonances. It will allow the rapid assignment of backbone resonances with high sensitivity and the determination of a reasonable structural model of a protein based on limited NOE restraints, an application that is of increasing importance as data from the large number of genome sequencing projects accumulates. The method that we propose should also be of utility in extending the use of NMR spectroscopy to determine the structures of larger proteins.

Practical applications of time-averaged restrained molecular dynamics to ligand-receptor systems: FK506 bound to the Q50R,A95H,K98I triple mutant of FKBP-13

The ability of time-averaged restrained molecular dynamics (TARMD) to escape local low-energy conformations and explore conformational space is compared with conventional simulated-annealing methods. Practical suggestions are offered for performing TARMD calculations with ligand-receptor systems, and are illustrated for the complex of the immunosuppressant FK506 bound to Q50R,A95H,K98I triple mutant FKBP-13. The structure of (13)C-labeled FK506 bound to triple-mutant FKBP-13 was determined using a set of 87 NOE distance restraints derived from HSQC-NOESY experiments. TARMD was found to be superior to conventional simulated-annealing methods, and produced structures that were conformationally similar to FK506 bound to wild-type FKBP-12. The individual and combined effects of varying the NOE restraint force constant, using an explicit model for the protein binding pocket, and starting the calculations from different ligand conformations were explored in detail.

Folded conformations of the delta-selective opioid dermenkephalin with head-to-tail interactions. A simulated annealing study through NMR restraints.

Despite similar tripeptide N-termini, dermorphin (Tyr-D-Ala-Phe-Gly-Tyr-Pro-Ser-NH2) and dermenkephalin (Tyr-D-Met-Phe-His-Leu-Met-Asp-NH2), naturally occuring opioid peptides from frog skin, exhibit high affinity but contrasting selectivity for the mu- and delta-opioid receptors, respectively. Structure-activity relationship studies have shown that the N-terminal tripeptide, Tyr-D-Xaa-Phe (where Xaa is either Ala or Met), is necessary for binding with both the mu- and delta-receptors while the nature and/or the conformation of the C-terminus His-Leu-Met-Asp-NH2 of dermenkephalin are responsible for addressing the peptide to the delta-receptor. In order to examine the conformational characteristics that are related to the selectivity of dermenkephalin towards the delta-receptor, 50 NOE restraints (10 between non-adjacent residues), and 7 dihedral angles, derived from a two-dimensional 1H-NMR study of dermenkephalin in dimethyl sulfoxide, were used in simulated annealing and energy minimization procedures. Twenty-four resulting conformers (60% of the generated structures) with no severe distance restraint violation were pooled into seven groups and three related families. These 24 conformers show close proximity between the two methionine residues, S-shaped structures, mean planes of N-terminal and C-terminal moieties almost at right angles to each other, a C-terminus region above the plane of the N-terminal region and g- as preferential orientation in the side chain of Phe. Aside these similarities, families of conformers differ by the preferential orientation in the side chain of Tyr (t or g-) and proximity between Tyr and Asp, or Tyr and the C-terminus. In contrast to previous models, practically no beta-turn structures exist for dermenkephalin, most of the NH hydrogen bonds participating to gamma-turns. The possible relationship between the conformational characteristics of dermenkephalin and the delta-opioid receptor selectivity is discussed.

Solution structure of the E-domain of staphylococcal protein A.

The E-domain of staphylococcal protein A is one of five homologous IgG-binding domains designated E, D, A, B, and C that comprise the extracellular portion of protein A. The E-domain binds tightly to Fc fragments of IgG and binds certain Fv fragments with micromolar affinity. To explore further the structural features of Fc binding by protein A, and as a first step in developing a structural understanding of E-domain/Fv complex formation, we have determined the solution structure of the uncomplexed E-domain using 2D homonuclear and heteronuclear NMR spectroscopy. Complete 1H and 15N resonance assignments were obtained, and the structure was determined from 383 NOE-derived distance restrains, 34 phi and 19 chi 1 dihedral angle restraints, and 54 restraints for 27 H-bonds. 3JH alpha-H beta coupling constants and long-range NOEs involving Phe11 indicate the side chain exists in more than one conformation with differing chi 1 values. NOE restraints that were incompatible with chi 1 = -60 degrees were removed from one set of structure calculations, and those incompatible with chi 1 = 180 degrees were removed from a second set to allow Phe11 to explore both rotamer wells. Thus, two sets of 20 final structures, having no distance or dihedral angle restraint violations greater than 0.12 A or 1.6 degrees, respectively, represent the solution structure of the E-domain. Backbone atomic rms differences with respect to the mean coordinates for each set of 20 structures for residues 8-53 averaged 0.41 +/- 0.06 and 0.35 +/- 0.06 A. No significant differences in the overall structure result from the different orientations of Phe11. The solution structure of the E-domain consists of three alpha-helices that pack together to form a compact helical bundle. A detailed comparison between the E-domain ensembles and the previously determined structure for the B-domain in complex with Fc indicates that only the 180 degrees chi 1 rotamer of Phe11 is competent for binding; the -60 degrees chi 1 rotamer must reorient to 180 degrees to create a cavity that is filled by Ile253 from the CH2 domain of Fc in the Fc-bound complex.

Assignment of aliphatic side-chain 1HN/15N resonances in perdeuterated proteins

The perdeuteration of aliphatic sites in large proteins has been shown to greatly facilitate the process of sequential backbone and side-chain 13C assignments and has also been utilized in obtaining long-range NOE distance restraints for structure calculations. To obtain the maximum information from a 4D 15N/15N-separated NOESY, as many main-chain and side-chain 1HN/15N resonances as possible must be assigned. Traditionally, only backbone amide 1HN/15N resonances are assigned by correlation experiments, whereas slowly exchanging side-chain amide, amino, and guanidino protons are assigned by NOEs to side-chain aliphatic protons. In a perdeuterated protein, however, there is a minimal number of such protons. We have therefore developed several gradient-enhanced and sensitivity-enhanced pulse sequences, containing water-flipback pulses, to provide through-bond correlations of the aliphatic side-chain 1HN/15N resonances to side-chain 13C resonances with high sensitivity: NH2-filtered 2D 1H-15N HSQC(H2N-HSQC), 3D H2N(CO)C gamma/beta and 3D H2N(COC gamma/beta)C beta/alpha for glutamine and asparagine side-chain amide groups; 2D refocused H(N epsilon/zeta)C delta/epsilon and H(N epsilon/zeta C delta/epsilon)C gamma/delta for arginine side-chain amino groups and non-refocused versions for lysine side-chain amino groups; and 2D refocused H(N epsilon)C zeta and nonrefocused H(N epsilon, eta)C zeta for arginine side-chain guanidino groups. These pulse sequences have been applied to perdeuterated 13C-/15N-labeled human carbonic anhydrase II (2H-HCA II). Because more than 95% of all side-chain 13C resonances in 2H-HCA II have already been assigned with the C(CC)(CO)NH experiment, the assignment of the side-chain 1HN/15N resonances has been straightforward using the pulse sequences mentioned above. The importance of assigning these side-chain HN protons has been demonstrated by recent studies in which the calculation of protein global folds was simulated using only 1HN-1HN NOE restraints. In these studies, the inclusion of NOE restraints to side-chain HN protons significantly improved the quality of the global fold that could be determined for a perdeuterated protein [R.A. Venters et al. (1995) J. Am. Chem. Soc., 117, 9592-9593].

Do NOE distances contain enough information to assess the relative populations of multi-conformer structures?

The feasibility of determining the relative populations of multi-conformer structures from NOE-derived distances alone is assessed. Without cross-validation of the NOE restraints, any population ratio can be refined to a similar quality of the fit. Complete cross-validation provides a less biased measure of fit and allows the estimation of the correct population ratio when used in conjunction with very tight distance restraints. With the qualitative distance restraints most commonly used in NMR structure determination, cross-validation is unsuccessful in providing the correct answer. Other experimental sources are therefore needed to determine relative populations of multi-conformer structures.

Solution structure of the catalytic domain of human stromelysin complexed with a hydrophobic inhibitor.

Stromelysin, a representative matrix metalloproteinase and target of drug development efforts, plays a prominent role in the pathological proteolysis associated with arthritis and secondarily in that of cancer metastasis and invasion. To provide a structural template to aid the development of therapeutic inhibitors, we have determined a medium-resolution structure of a 20-kDa complex of human stromelysin's catalytic domain with a hydrophobic peptidic inhibitor using multinuclear, multidimensional NMR spectroscopy. This domain of this zinc hydrolase contains a mixed beta-sheet comprising one antiparallel strand and four parallel strands, three helices, and a methionine-containing turn near the catalytic center. The ensemble of 20 structures was calculated using, on average, 8 interresidue NOE restraints per residue for the 166-residue protein fragment complexed with a 4-residue substrate analogue. The mean RMS deviation (RMSD) to the average structure for backbone heavy atoms is 0.91 A and for all heavy atoms is 1.42 A. The structure has good stereochemical properties, including its backbone torsion angles. The beta-sheet and alpha-helices of the catalytic domains of human stromelysin (NMR model) and human fibroblast collagenase (X-ray crystallographic model of Lovejoy B et al., 1994b, Biochemistry 33:8207-8217) superimpose well, having a pairwise RMSD for backbone heavy atoms of 2.28 A when three loop segments are disregarded. The hydroxamate-substituted inhibitor binds across the hydrophobic active site of stromelysin in an extended conformation. The first hydrophobic side chain is deeply buried in the principal S'1 subsite, the second hydrophobic side chain is located on the opposite side of the inhibitor backbone in the hydrophobic S'2 surface subsite, and a third hydrophobic side chain (P'3) lies at the surface.

Refined solution structure of the oligomerization domain of the tumour suppressor p53

The NMR solution structure of the oligomerization domain of the tumour suppressor p53 (residues 319-360) has been refined. The structure comprises a dimer of dimers, oriented in an approximately orthogonal manner. The present structure determination is based on 4,472 experimental NMR restraints which represents a three and half fold increase over our previous work in the number of NOE restraints at the tetramerization interface. A comparison with the recently solved 1.7 A resolution X-ray structure shows that the structures are very similar and that the average angular root-mean-square difference in the interhelical angles is about 1 degree. The results of recent extensive mutagenesis data and the possible effects of mutations which have been identified in human cancers are discussed in the light of the present structure.

NMR observation of interactions in the combining site region of an antibody using a spin‐labeled peptide antigen and NOESY difference spectroscopy

A spin-labeled peptide antigen (TEMPOVEVPGSQHIDSQ) was used to measure NOESY difference spectra that show interactions in the binding site region of the Fab fragment of the anti-cholera toxin peptide antibody TE33. In addition to identification of peptide-Fab interactions and interactions within the bound peptide, these difference spectra show well-resolved cross peaks due to interactions within the large Fab fragment (50 kDa). These difference spectra indicate that the conformational changes in the Fab upon peptide binding are confined to the combining site region of the antibody. The NOESY difference spectra of selectively deuterated Fab molecules were used in combination with HOHAHA measurements to assign the interactions to amino acid type and to identify the interactions within the Fab as either inter- or intraresidue interactions. The assignment of interactions within the Fab to corresponding aromatic residues in the Fab sequence was facilitated by an earlier NMR-derived model calculated on the basis of NOE restraints on Fab-peptide and intra-bound-peptide distances. The new restraints on distances within the Fab, combined with the previously obtained restraints, were used to generate a refined NMR-derived model for the TE33-peptide complex.

Solution structure of FK506 bound to the R42K, H87V double mutant of FKBP-12.

The binding of the FK506/FKBP-12 complex to calcineurin (CN), its putative target for immunosuppression, involves recognition of solvent-exposed regions of the ligand as well as FKBP-12 residues near the active site. The R42K, H87V double mutation of FKBP-12 decreases the CN affinity of the complex by 550-fold [Aldape, R. A., Futer, O., DeCenzo, M. T., Jarrett, B. P., Murcko, M. A., & Livingston, D. J. (1992) J. Biol. Chem. 267, 16029-16032]. This work reports the solution structure of 13C-labeled FK506 bound to R42K, H87V FKBP-12. Assignments and NOE measurements at three mixing times were made from inverse-detected 1H-13C NMR experiments. Structures were calculated by several different methods, including distance geometry, restrained molecular dynamics, and molecular dynamics with time-averaged restraints. The NMR structures of the ligand are very well defined by the NOE restraints and differ slightly from the X-ray structure in regions that are involved in crystal packing. Comparison with the NMR structure of FK506 bound to wild-type FKBP-12 reveals that the R42K, H87V mutation causes the ligand backbone near C16 to move by 2.5 to 4.5 A, reorients 15-MeO by 90 degrees, and shifts 13-MeO by approximately 1.5 A. FK506 appears to undergo a concerted, mutationally induced shift in the binding pocket, with the greatest changes occurring in the effector region of the drug. The altered effector conformation of mutant-bound FK506 may perturb interactions between the drug and CN, thus accounting for the effect of the double mutation upon the CN inhibitory activity of the complex.

The solution conformation of hyaluronan: a combined NMR and molecular dynamics study.

Hyaluronan (HA) is a negatively charged glycosaminoglycan that exhibits a wide variety of biological effects mediated by binding to cell-surface and extracellular matrix proteins (hyaladherins). Short HA oligosaccharides have been shown to retain the specific interactions and biological effects of high molecular weight HA. Although it has a simple disaccharide repeating unit, the aqueous solution conformation of HA has been very difficult to determine because of strong coupling and overlapping resonances. In this study, we propose aqueous solution conformations for an octasaccharide of HA, derived from proton-proton NOE data and restrained molecular dynamics. To overcome spectral overlap and strong coupling, alternate methods for extracting distance restraints were employed. Restrained molecular dynamics calculations yielded one set of interglycosidic angle values for the beta (1,3) linkage (phi 13 = 46 degrees, psi 13 = 24 degrees). In contrast, two sets of values for the beta (1,4) linkage were consistent with the NOE restraints (phi 14 = 24 degrees, psi 14 = -53 degrees or phi 14 = 48 degrees, psi 14 = 8 degrees). The potential difference in flexibility for the two linkages is consistent with unrestrained as well as the restrained molecular dynamics trajectories described here. The conformational parameters obtained from restrained molecular dynamics are used to predict helical parameters of high molecular weight HA and will provide a basis for studies of HA binding to proteins.

Multiple conformations and proline cis-trans isomerization in salmon calcitonin: a combined nuclear magnetic resonance, distance geometry, and molecular mechanics study.

The relationship between multiple conformations and proline cis-trans isomerization in salmon calcitonin (sCT) has been investigated by 1H-NMR, distance geometry, and molecular mechanics. In different solvents, the isomerization of Pro23 induces resonance heterogeneity for amino acids adjacent to it. NOESY experiments have related such heterogeneity to two different isoforms of the hormone, which interconvert into each other in polar solvents as well as in SDS micelles. The cis-trans ratio was found to depend upon the structure of the hormone. In water, the random-coil sCT showed a 35% cis population, while 16% cis was observed in the folded SDS-bound hormone. Such a decrease contributes, per se, a Gibbs free energy stabilization of 3.10 kJ mol-1. Except for the 21-25 region, a common NOE pathway was found for both isomers. This is in agreement with the common secondary structure for both isoforms: a central helix and an extended C-terminal segment interacting with it. Calculations indicated that while the trans isomer is helical in the Thr6-Tyr22 region, a shorter helix (Thr6-Lys18) is present in the cis isoform. It is concluded that isomerization of Pro23 does not alter the three-dimensional structure of the hormone, although trans structures show a lower average NOE root mean square deviation and a higher relative stability. Both cis and trans structures satisfactorily reproduce the experimental data, although they do not fulfill the whole set of NOE restraints, pointing out the danger involved in structure calculation whenever the cis-trans equilibrium does not affect the global fold.

Three-dimensional structure in solution of neurotoxin III from the sea anemone Anemonia sulcata.

The three-dimensional structure in aqueous solution of the 27-residue polypeptide neurotoxin Anemonia sulcata toxin III (ATX III) has been determined from 1H NMR data. As ATX III self-associates in the millimolar concentration range, causing a marked concentration dependence for the chemical shifts of several residues [Norton, R. S., Cross, K., Braach-Maksvytis, V., & Wachter, E. (1993) Biochem. J. 293, 545-551], it was necessary to record NOESY spectra over a range of concentrations in order to eliminate any intermolecular interactions from the NOE restraint set. The pairings of the six half-cystine residues were also unknown and had to be determined (as 3-17, 4-11, and 6-22) from preliminary structure calculations performed using both upper bound distance restraints from NOESY data and a substantial number of lower bound restraints inferred from the absence of NOESY cross-peaks. Final structures were determined, using the program X-PLOR, from interproton distance restraints inferred from NOEs, backbone and side chain dihedral angle restraints from spin-spin coupling measurements, and a smaller number of lower bound restraints. Stereospecific assignments for 11 beta-methylene pairs were also included. The final set of 28 structures had an average pairwise RMS difference of 1.32 A over the backbone heavy atoms (N, C alpha, and C) and 2.18 A over all heavy atoms. For the well-defined region encompassing residues 3-22, the corresponding values were 0.62 and 1.28 A, respectively. ATX III adopts a compact structure containing four reverse turns (a distorted type I beta-turn at residues 6-9, a type I beta-turn at residues 8-11, and inverse gamma-turns at residues 12-14 and 15-17) and two other chain reversals, but no regular alpha-helix or beta-sheet. Several of the residues most affected by aggregation are located on the surface of the molecule, forming a hydrophobic patch which may constitute part of the sodium channel binding surface. Possible relationships between the structure of ATX III and those of other sea anemone toxins that interact with the same site on the voltage-gated sodium channel are considered.

Solution structure of the amino-terminal fragment of urokinase-type plasminogen activator.

The amino-terminal fragment (ATF) of urokinase-type plasminogen activator is a two domain protein which consists of a growth factor and a kringle domain. The 1H, 13C, and 15N chemical shifts of this protein have been assigned using heteronuclear two- and three-dimensional NMR experiments on selective and uniformly 15N- and 15N/13C-labeled protein isolated from mammalian cells that overexpress the protein. The chemical shift assignments were used to interpret the NOE data which resulted in a total of 1299 NOE restraints. The NOE restraints were used along with 27 phi angle restraints and 21 hydrogen-bonding restraints to produce 15 low energy structures. The individual domains in the structures are highly converged, but the two domains are structurally independent. The root mean square deviations (rmsd) between residues 11-46 in the growth factor domain and the mean atomic coordinates were 0.99 +/- 0.2 for backbone heavy atoms and 1.65 +/- 0.2 for all non-hydrogen atoms. For residues 55-130 in the kringle domain, the rmsd was 0.84 +/- 0.2 for backbone heavy atoms and 1.42 +/- 0.2 for all non-hydrogen atoms. The overall structures of the individual domains are very similar to the structures of homologous proteins. However, important structural differences between the growth factor and other homologous proteins were observed in the region which has been implicated in binding the urokinase receptor which may explain, in part, why other growth factors show no appreciable affinity for the urokinase receptor.

Peptide flexibility and calculations of an ensemble of molecules

Conformational averaging fast on the NMR time scale has been examined by computer simulations of multiple copies of the molecule and application of NOE and coupling constant restraints as an ensemble average. The calculation is illustrated for a model cyclic peptide, cyclo[-D-Pro-Ala 2 -Ala 3 -Ala 4 -Ala 5 -], for which conformational averaging is taking place. There is a well-defined type II' β-turn about the D-Pro-Ala 2 , while no single conformation can be ascribed to the other half of the molecule which fulfills the NMR observables

Time- and ensemble-averaged direct NOE restraints

NMR data are collected as time- and ensemble-averaged quantities. Yet, in commonly used methods for structure determination of biomolecules, structures are required to satisfy simultaneously a large number of constrainsts. Recently, however, methods have been developed that allow a better fit of the experimental data by the use of time- or ensemble-averaged restraints. Thus far, these methods have been applied to structure refinement using distance and J-coupling restraints. In this paper, time and ensemble averaging is extended to the direct refinement with experimental NOE data. The implementation of time- and ensemble-averaged NOE restraints in DINOSAUR is described and illustrated with experimental NMR data for crambin, a 46-residue protein. Structure refinement with both time- and ensemble-averaged NOE restraints results in lower R-factors, indicating a better fit of the experimental NOE data.

Conformational Variability of Chicken Cystatin: Comparison of Structures Determined by X-ray Diffraction and NMR Spectroscopy

The structural model derived from X-ray crystallography for unphosphorylated wild-type chicken cystatin is compared with two chicken cystatin structures derived from NMR spectroscopy: the phosphorylated wild-type and the genetically engineered variant AEF-S1M-M29I-M89L. The comparison shows the same overall fold, but also significant differences in structurally variable segments of the polypeptide chain. The largest such segment, comprising residues 71 to 89, is a region characteristic of the family 2 cystatin inhibitors which contains a disulphide bridge (71-81) and the phosphorylation site (Ser80) discussed in the accompanying article. In the crystal structure, the segment 71 to 76 is found as a flexible loop, 77 to 85 as an α-helical segment, and 86 to 89 is completely undefined. The solution NMR structures on the other hand are disordered in the initial segment 72 to 80, have an extended conformation at 81 to 83 in contact with the β-sheet, and clearly show a β-turn at residues 87 to 90. The segment comprising residues 53 to 57, with smaller variability, is of particular interest as the hairpin loop conserved throughout the cystatin superfamily which binds to the cysteine proteinase. In most of the solution NMR structures, this segment adopts a conformation more like that of stefin B, a family 1 cystatin inhibitor, as was observed in the crystal structure of its inhibitory complex with papain. The differences between the structures are rationalized by an examination of the crystal contacts generated by hypothetical crystal packing of the NMR structures. Additionally, the X-ray refinement shows evidence of conformational disorder in the crystal. Joint refinement with NOE restraints and reflection data does not produce a structure to satisfy the restraints of both methods.