Evaluation of site-directed spin labeling for characterizing protein-ligand complexes using simulated restraints

Simulation studies have been performed to evaluate the utility of site-directed spin labeling for determining the structures of protein-ligand complexes, given a known protein structure. Two protein-ligand complexes were used as model systems for these studies: a 1.9-A-resolution x-ray structure of a dihydrofolate reductase mutant complexed with methotrexate, and a 1.5-A-resolution x-ray structure of the V-Src tyrosine kinase SH2 domain complexed with a five-residue phosphopeptide. Nitroxide spin labels were modeled at five dihydrofolate reductase residue positions and at four SH2 domain residue positions. For both systems, after energy minimization, conformational ensembles of the spin-labeled residues were generated by simulated annealing while holding the remainder of the protein-ligand complex fixed. Effective distances, simulating those that could be obtained from (1)H-NMR relaxation measurements, were calculated between ligand protons and the spin labels. These were converted to restraints with several different levels of precision. Restrained simulated annealing calculations were then performed with the aim of reproducing target ligand-binding modes. The effects of incorporating a few supplementary short-range (< or =5.0 A) distance restraints were also examined. For the dihydrofolate reductase-methotrexate complex, the ligand-binding mode was reproduced reasonably well using relatively tight spin-label restraints, but methotrexate was poorly localized using loose spin-label restraints. Short-range and spin-label restraints proved to be complementary. For the SH2 domain-phosphopeptide complex without the short-range restraints, the peptide did not localize to the correct depth in the binding groove; nevertheless, the orientation and internal conformation of the peptide was reproduced moderately well. Use of the spin-label restraints in conjunction with the short-range restraints resulted in relatively well defined structural ensembles. These results indicate that restraints derived from site-directed spin labeling can contribute significantly to defining the orientations and conformations of bound ligands. Accurate ligand localization appears to require either a few supplementary short-range distance restraints, or relatively tight spin-label restraints, with at least one spin label positioned so that some of the restraints draw the ligand into the binding pocket in the latter case.

Backbone and side chain dynamics of uncomplexed human adipocyte and muscle fatty acid-binding proteins

Adipocyte lipid-binding protein (A-LBP) and muscle fatty acid-binding protein (M-FABP) are members of a family of small ( approximately 15 kDa) cytosolic proteins that are involved in the metabolism of fatty acids and other lipid-soluble molecules. Although highly homologous (65%) and structurally very similar, A-LBP and M-FABP display distinct ligand binding characteristics. Since ligand binding may be influenced by intrinsic protein dynamical properties, we have characterized the backbone and side chain dynamics of uncomplexed (apo) human A-LBP and M-FABP. Backbone dynamics were characterized by measurements of 15N T1 and T2 values and ¿1H¿-15N NOEs. These data were analyzed using model-free spectral density functions and reduced spectral density mapping. The dynamics of methyl-containing side chains were charaterized by measurements of 2H T1 and T1rho relaxation times of 13C1H22H groups. The 2H relaxation data were analyzed using the model-free approach. For A-LBP, 15N relaxation data were obtained for 111 residues and 2H relaxation data were obtained for 42 methyl groups. For M-FABP, 15N relaxation data were obtained for 111 residues and 2H relaxation data were obtained for 53 methyl groups. The intrinsic flexibilities of these two proteins are compared, with particular emphasis placed on binding pocket residues. There are a number of distinct dynamical differences among corresponding residues between the two proteins. In particular, many residues display greater backbone picosecond to nanosecond and/or microsecond to millisecond time scale mobility in A-LBP relative to M-FABP, including F57, K58, and most residues in alpha-helix 2 (residues 28-35). Variations in the dynamics of this region may play a role in ligand selectivity. The side chains lining the fatty acid binding pocket display a wide range of motional restriction in both proteins. Side chains showing distinct dynamical differences between the two proteins include those of residues 20, 29, and 51. This information provides a necessary benchmark for determining dynamical changes induced by ligand binding and may ultimately lead to an enhanced understanding of ligand affinity and selectivity among fatty acid-binding proteins.

Solution structure of human CTLA-4 and delineation of a CD80/CD86 binding site conserved in CD28

The structure of human CTLA-4 reveals that residues Met 99, Tyr 100 and Tyr 104 of the M99YPPPY104 motif are adjacent to a patch of charged surface residues on the A'GFCC' face of the protein. Mutation of these residues, which are conserved in the CTLA-4/CD28 family, significantly reduces binding to CD80 and/or CD86, implicating this patch as a ligand binding site.

Characterization of NADP+ binding to perdeuterated MurB: backbone atom NMR assignments and chemical-shift changes

Backbone-atom resonances have been assigned for both the substrate-free and the NADP+-complexed forms of UDP-N-acetylenolpyruvylglucosamine reductase (MurB), a monomeric, 347-residue (38.5 kDa) flavoenzyme essential for bacterial cell-wall biosynthesis. NMR studies were performed using perdeuterated, uniformly 13C/15N-labeled samples of MurB. In the case of substrate-free MurB, one or more backbone atoms have been assigned for 334 residues (96%). The assigned backbone atoms include 309 1HN and 15N atoms (94%), 315 13CO atoms (91%), 331 13C(alpha) atoms (95%), and 297 13C(beta) atoms (93%). For NADP+-complexed MurB, one or more backbone atoms have been assigned for 313 residues (90%); these include 283 1HN and 15N atoms (86%), 305 13CO atoms (88%), 310 13C(alpha) atoms (89%), and 269 13C(beta) atoms (84%). The strategies used for obtaining resonance assignments are described in detail. Information on the secondary structure in solution for both the substrate-free and NADP+-complexed forms of the enzyme has been derived both from 13C(alpha) and 13C(beta) chemical-shift deviations from random-coil values and from 1HN-1HN NOEs. These data are compared to X-ray crystallographic structures of substrate-free MurB and MurB complexed with the UDP-N-acetylglucosamine enolpyruvate (UNAGEP) substrate. NADP+ binding induces significant chemical-shift changes in residues both within the known UNAGEP and FAD binding pockets and within regions known to undergo conformational changes upon UNAGEP binding. The NMR data indicate that NADP+ and UNAGEP utilize the same binding pocket and, furthermore, that the binding of NADP+ induces structural changes in MurB. Finally, many of the residues within the UNAGEP/NADP+ binding pocket were difficult to assign due to dynamic processes which weaken and/or broaden the respective resonances. Overall, our results are consistent with MurB having a flexible active site.

High-resolution solution structure of siamycin II: novel amphipathic character of a 21-residue peptide that inhibits HIV fusion

The 21-amino acid peptides siamycin II (BMY-29303) and siamycin I (BMY-29304), derived from Streptomyces strains AA3891 and AA6532, respectively, have been found to inhibit HIV-1 fusion and viral replication in cell culture. The primary sequence of siamycin II is CLGIGSCNDFAGCGYAIVCFW. Siamycin I differs by only one amino acid; it has a valine residue at position 4. In both peptides, disulfide bonds link Cys1 with Cys13 and Cys7 with Cys19, and the side chain of Asp9 forms an amide bond with the N-terminus. Siamycin II, when dissolved in a 50:50 mixture of DMSO and H2O, yields NOESY spectra with exceptional numbers of cross peaks for a peptide of this size. We have used 335 NOE distance constraints and 13 dihedral angle constraints to generate an ensemble of 30 siamycin II structures; these have average backbone atom and all heavy atom rmsd values to the mean coordinates of 0.24 and 0.52 A, respectively. The peptide displays an unusual wedge-shaped structure, with one face being predominantly hydrophobic and the other being predominantly hydrophilic. Chemical shift and NOE data show that the siamycin I structure is essentially identical to siamycin II. These peptides may act by preventing oligomerization of the HIV transmembrane glycoprotein gp41, or by interfering with interactions between gp41 and the envelope glycoprotein gp120, the cell membrane or membrane-bound proteins [Frèchet, D. et al. (1994) Biochemistry, 33, 42-50]. The amphipathic nature of siamycin II and siamycin I suggests that a polar (or apolar) site on the target protein may be masked by the apolar (or polar) face of the peptide upon peptide/protein complexation.

Refined solution structure of human profilin I

Profilin is a ubiquitous eukaryotic protein that binds to both cytosolic actin and the phospholipid phosphatidylinositol-4,5-bisphosphate. These dual competitive binding capabilities of profilin suggest that profilin serves as a link between the phosphatidyl inositol cycle and actin polymerization, and thus profilin may be an essential component in the signaling pathway leading to cytoskeletal rearrangement. The refined three-dimensional solution structure of human profilin I has been determined using multidimensional heteronuclear NMR spectroscopy. Twenty structures were selected to represent the solution conformational ensemble. This ensemble of structures has root-mean-square distance deviations from the mean structure of 0.58 A for the backbone atoms and 0.98 A for all non-hydrogen atoms. Comparison of the solution structure of human profilin to the crystal structure of bovine profilin reveals that, although profilin adopts essentially identical conformations in both states, the solution structure is more compact than the crystal structure. Interestingly, the regions that show the most structural diversity are located at or near the actin-binding site of profilin. We suggest that structural differences are reflective of dynamical properties of profilin that facilitate favorable interactions with actin. The global folding pattern of human profilin also closely resembles that of Acanthamoeba profilin I, reflective of the 22% sequence identity and approximately 45% sequence similarity between these two proteins.

Sequential 1H, 13C, and 15N NMR assignments and solution conformation of apokedarcidin

Kedarcidin is a recently discovered antitumor antibiotic chromoprotein. The solution conformation of the kedarcidin apoprotein (114 residues) has been characterized by heteronuclear multidimensional NMR spectroscopy. Sequence-specific backbone atom resonance assignments were obtained for a uniformly 13C/15N-enriched sample of apokedarcidin via a semiautomated analysis of 3D HNCACB, 3D CBCA-(CO)NH, 4D HNCAHA, 4D HN(CO)CAHA, 3D HBHA(CO)NH, and 3D HNHA(Gly) spectra. Side-chain assignments were subsequently obtained by analysis of (primarily) 3D HCCH-TOCSY and HCCH-COSY spectra. A qualitative analysis of the secondary structure is presented on the basis of 3J alpha NH coupling constants, deviations of 13C alpha and 13C beta chemical shifts from random coil values, and NOEs observed in 3D 15N- and 13C-edited NOESY-HSQC spectra. This analysis revealed a four-stranded antiparallel beta-sheet, a three-stranded antiparallel beta-sheet, and two two-standed antiparallel beta-sheets. The assignments of cross-peaks in the 3D NOESY spectra were assisted by reference to a preliminary model of apokedarcidin built using the program CONGEN starting from the X-ray structure of the homologous protein aponeocarzinostatin. An ensemble of 15 apokedarcidin solution structures has been generated by variable target function minimization (DIANA program) and refined by simulated annealing (X-PLOR program). The average backbone atom root-mean-square difference between the individual structures and the mean coordinates is 0.68 +/- 0.08 A. The overall fold of apokedarcidin is well-defined; it is composed of an immunoglobulin-like seven-stranded antiparallel beta-barrel and a subdomain containing two antiparallel beta-ribbons. Highly similar tertiary structures have been previously reported for the related proteins neocarzinostatin, macromomycin, and actinoxanthin. Important structural features are revealed, including the dimensions of the chromophore-binding pocket and the locations of side chains that are likely to be involved in chromophore stabilization.

Solution structure of an isolated antibody VL domain

The solution structure of the isolated VL domain of the anti-digoxin antibody 26-10 has been determined using data derived from heteronuclear multi-dimensional nuclear magnetic resonance (n.m.r.) experiments. Analytical ultracentrifugation and n.m.r. data demonstrate that the VL domain is only weakly associating (Kd = 2.5 (+/- 0.7) mM) and that it experiences a rapid monomer/dimer equilibrium under the n.m.r. experimental conditions. Therefore, the results reported here represent the first structure determination of an antibody VL domain in the absence of fixed quaternary interactions. The structure determination is based on 930 proton-proton distance constraints, 113 dihedral angle constraints, and 46 hydrogen bond constraints. Eighty initial structures were calculated with the variable target function program DIANA; of these, 31 were accepted on the basis of satisfaction of constraints (no distance constraint violations > 0.5 A; target function < 3.0 A2). Accepted DIANA structures were refined by restrained energy minimization using the X-PLOR program. The 15 best energy-minimized DIANA structures were chosen as a representative ensemble of solution conformations. The average root-mean-square differences (r.m.s.d.) between the individual structures of this ensemble and the mean coordinates is 0.85 (+/- 0.10) A for all backbone atoms and 1.29 (+/- 0.10) A for all heavy atoms. For beta-strands A, B, C, D, E and F, the average backbone atom r.m.s.d. to the mean structure is 0.46 (+/- 0.06) A. A higher-resolution ensemble, with all backbone atom and all heavy atom r.m.s.d.s. to the mean coordinates of 0.54 (+/- 0.08) A and 0.98 (+/- 0.12) A, respectively, was obtained by X-PLOR simulated annealing refinement of the 15 energy-minimized DIANA structures. A detailed analysis of the original ensemble of 15 energy-minimized DIANA structures is presented, as this ensemble retains a broader, and possibly more realistic, sampling of conformation space. The backbone atom and all heavy atom r.m.s.d.s between the mean energy-minimized DIANA structure and the X-ray derived coordinates of the VL domain within the Fab/digoxin complex are 1.05 A and 1.56 A, respectively. Subtle differences between the solution and X-ray structures occur primarily in CDR2, CDR3, beta-strands A, F and G, and localized regions of hydrophobic packing. Overall, these results demonstrate that the 26-10 VL domain conformation is determined primarily by intradomain interactions, and that quaternary VL-VH association induces relatively minor conformational adjustments.

Relaxation study of the backbone dynamics of human profilin by two-dimensional 1H-15N NMR

The dynamic properties of 111 backbone HN sites in uncomplexed human profilin, a protein of 139 residues, have been characterized by two-dimensional inverse-detected 1H-15N NMR spectroscopy. Heteronuclear (1H)-15N nuclear Overhauser effects and 15N longitudinal and transverse relaxation rates have been analyzed in terms of model-free spectral density functions and exchange contributions to transverse relaxation rates. Relatively high mobilities on the nanosecond time-scale are observed for Asp26 and Ser27, which form part of a loop connecting beta-strands A and B, and for Thr92 through Ala95, which are in a loop connecting beta-strands E and F. Significant exchange contributions, indicative of motions on the microsecond to millisecond time-scale, have been obtained for 30 residues. These include Leu77, Asp80 and Gly81 of a loop between beta-strands D and E, Ser84 and Met85 of beta-strand E, Gly121 of a loop connecting beta-strand G and the C-terminal helix, and Gln138, which is next to the C-terminal residue Tyr139. Some of the regions showing high flexibility in profilin are known to be involved in poly-L-proline binding.

Characterization of the three-dimensional solution structure of human profilin: 1H, 13C, and 15N NMR assignments and global folding pattern

Human profilin is a 15-kDa protein that plays a major role in the signaling pathway leading to cytoskeletal rearrangement. Essentially complete assignment of the 1H, 13C, and 15N resonances of human profilin have been made by analysis of multidimensional, double- and triple-resonance nuclear magnetic resonance (NMR) experiments. The deviation of the 13C alpha and 13C beta chemical shifts from their respective random coil values were analyzed and correlate well with the secondary structure determined from the NMR data. Twenty structures of human profilin were refined in the program X-PLOR using a total of 1186 experimentally derived conformational restraints. The structures converged to a root mean squared distance deviation of 1.5 A for the backbone atoms. The resultant conformational ensemble indicates that human profilin is an alpha/beta protein comprised of a seven-stranded, antiparallel beta-sheet and three helices. The secondary structure elements for human profilin are quite similar to those found in Acanthamoeba profilin I [Archer, S. J., Vinson, V. K., Pollard, T. D., & Torchia, D. A. (1993), Biochemistry 32, 6680-6687], suggesting that the three-dimensional structure of Acanthamoeba profilin I should be analogous to that determined here for human profilin. The structure determination of human profilin has facilitated the sequence alignment of lower eukaryotic and human profilins and provides a framework upon which the various functionalities of profilin can be explored. At least one element of the actin-binding region of human profilin is an alpha-helix. Two mechanisms by which phosphatidylinositol 4,5-bisphosphate can interfere with actin-binding by human profilin are proposed.

Micelle-bound conformational preferences of a peptide derived from a murine major histocompatibility complex class I molecule

Models of the micelle-bound conformation of a 17-residue major histocompatibility complex-derived peptide, [Ala85]Dk(69-85), have been determined by NMR spectroscopy and simulated annealing calculations. This peptide is a truncated, substituted version of Dk(61-85), which is a fragment of the murine major histocompatibility complex class I molecule H-2Dk. Dk(61-85) has been shown to adopt an ordered conformation required for augmentation of insulin-stimulated glucose uptake (Stagsted, J., Baase, W. A., Goldstein, A., and Olsson, L. (1991) J. Biol. Chem. 266, 12844-12847). [Ala85]Dk(69-85) retains full biological activity. Thirty-eight converged NMR structures of [Ala85]Dk(69-85) bound to dodecyl phosphocholine micelles have been generated. The NMR-derived models display a propensity for a type-I beta-bend involving residues 73-76 and an amphipathic helical region involving residues 77-84. CD spectra yield a helical content (8% at 20-25 degrees C) consistent with transient, partial helix formation. The relative orientation of the beta-bend region with respect to the helical region is not well defined by the NMR data. This may reflect true heterogeneity of the micelle-bound conformation. The NMR structures were compared with a model of [Ala85]Dk(69-85) derived from the x-ray coordinates of the human major histocompatibility complex class I allele HLA-Aw68 (Garrett, T. P. J., Saper, M. A., Bjorkmann, P. J., Strominger, T. L., and Wiley, D. C. (1989) Nature 342, 692-696). Structural features that are important for the bioactivity of [Ala85]Dk(69-85) are discussed with reference to reported structure-activity relationships (Stagsted, J., Mapelli, C., Myers, C., Matthews, B. W., Anfinsen, C. B., Goldstein, A., and Olsson, L. (1993) Proc. Natl. Acad. Sci. U.S.A., in press). A general description of the structural properties of the putative receptor site(s) that are likely to be required for binding [Ala85]Dk(69-85) is given.

Production and characterization of an antibody Fv fragment 15N-labeled in the VL domain only

Procedures for separating and recombining the light-chain variable (VL) and heavy-chain variable (VH) domains of antibody Fv fragments have been applied to produce a recombinant Fv of the anti-digoxin antibody 26-10 that is 15N-labeled exclusively in the VL domain. Comparison of a two-dimensional 1H-15N heteronuclear single-quantum correlation (HSQC) spectrum of the reconstituted Fv with a HSQC spectrum of a fully 15N-labeled Fv sample reveals that all 1H-15N correlations of the VL domain align precisely in both spectra. Assignments for 105 of the 106 backbone HN groups of the VL domain within the reconstituted Fv have been obtained by analysis of three-dimensional nuclear Overhauser effect spectroscopy-HSCQ and total correlation spectroscopy-HSQC spectra with reference to assignments previously reported for the isolated VL domain. Chemical shift differences between the isolated VL domain and the VL domain within the Fv are moderately correlated with proximity to the surfaces of the VH domain and bound hapten (ouabain) as defined by X-ray crystallography and molecular modeling. These results demonstrate that nuclear magnetic resonance studies of reconstituted antibody Fv fragments, in conjunction with investigations of isolated antibody domains, can yield extensive resonance assignments for the Fv. This will facilitate detailed studies of antigen-antibody and domain-domain interactions.

Characterization of the backbone dynamics of an anti-digoxin antibody VL domain by inverse detected 1H-15N NMR: comparisons with X-ray data for the Fab

The dynamic behavior of the polypeptide backbone of a recombinant antidigoxin antibody VL domain has been characterized by measurements of 15NT1 and T2 relaxation times, 1H-15N NOE values, and 1H-2H exchange rates. These data were acquired with 2D inverse detected heteronuclear 1H-15N NMR methods. The relaxation data are interpreted in terms of model free spectral density functions and exchange contributions to transverse relaxation rates R2 (= 1/T2). All characterized residues display low-amplitude picosecond time-scale librational motions. Fifteen residues undergo conformational changes on the nanosecond timescale, and 24 residues have significant R2 exchange contributions, which reflect motions on the microsecond to millisecond time-scale. For several residues, microsecond to millisecond motions of nearby aromatic rings are postulated to account for some or all of their observed R2 exchange contributions. The measured 1H-2H exchange rates are correlated with hydrogen bonding patterns and distances from the solvent accessible surface. The degree of local flexibility indicated by the NMR measurements is compared to crystallographic B-factors derived from X-ray analyses of the native Fab and the Fab/digoxin complex. In general, both the NMR and X-ray data indicate enhanced flexibility in the turns, hypervariable loops, and portions of beta-strands A, B, and G. However, on a residue-specific level, correlations among the various NMR data, and between the NMR and X-ray data, are often absent. This is attributed to the different dynamic processes and environments that influence the various observables. The combined data indicate that certain regions of the VL domain, including the three hypervariable loops, undergo dynamic changes upon VL:VH association and/or complexation with digoxin. Overall, the 26-10 VL domain exhibits relatively low flexibility on the ps-ns timescale. The possible functional consequences of this result are considered.

Aliphatic 1H and 13C resonance assignments for the 26-10 antibody VL domain derived from heteronuclear multidimensional NMR spectroscopy

Extensive 1H and 13C assignments have been obtained for the aliphatic resonances of a uniformly 13C- and 15N-labeled recombinant VL domain from the anti-digoxin antibody 26-10. Four-dimensional triple resonance NMR data acquired with the HNCAHA and HN(CO)CAHA pulse sequences [Kay et al. (1992) J. Magn. Reson., 98, 443-450] afforded assignments for the backbone HN, N, H alpha and C alpha resonances. These data confirm and extend HN, N and H alpha assignments derived previously from three-dimensional 1H-15N NMR studies of uniformly 15N-labeled VL domain [Constantine et al. (1992), Biochemistry, 31, 5033-5043]. The identified H alpha and C alpha resonances provided a starting point for assigning the side-chain aliphatic 1H and 13C resonances using three-dimensional HCCH-COSY and HCCH-TOCSY experiments [Clore et al. (1990), Biochemistry, 29, 8172-8184]. The C alpha and C beta chemical shifts are correlated with the VL domain secondary structure. The extensive set of side-chain assignments obtained will allow a detailed comparison to be made between the solution structure of the isolated VL domain and the X-ray structure of the VL domain within the 26-10 Fab.

Sequential 1H and 15N NMR assignments and secondary structure of a recombinant anti-digoxin antibody VL domain

A uniformly 15N-labeled recombinant light-chain variable (VL) domain from the anti-digoxin antibody 26-10 has been investigated by heteronuclear two-dimensional (2D) and three-dimensional (3D) NMR spectroscopy. Complementary homonuclear 2D NMR studies of the unlabeled VL domain were also performed. Sequence-specific assignments for 97% of the main-chain and 70% of the side-chain proton resonances have been obtained. Patterns of nuclear Overhauser effects observed in 2D NOESY, 3D NOESY-HSQC, and 3D NOESY-TOCSY-HSQC spectra afford a detailed characterization of the VL domain secondary structure in solution. The observed secondary structure--a nine-stranded antiparallel beta-barrel--corresponds to that observed crystallographically for VL domains involved in quaternary associations. The locations of slowly exchanging amide protons have been discerned from a 2D TOCSY spectrum recorded after dissolving the protein in 2H2O. Strands B, C, E, and F are found to be particularly stable. The possible consequences of these results for domain-domain interactions are discussed.

1H-n.m.r. studies of the fibronectin 13 kDa collagen-binding fragment. Evidence for autonomous conserved type I and type II domain folds

A 1H-n.m.r. study of a 117-residue (13 kDa) gelatin-binding fragment of human fibronectin, which contains the sixth (from the N-terminus) type I domain and the first type II domain, was undertaken. The resolution of the 1H-n.m.r. spectrum indicates that the domains are independent and mobile relative to each other. Analysis of two-dimensional 1H-n.m.r. experiments recorded at 500 MHz afforded spin-system identifications for all aromatic and a number of aliphatic residues. Utilizing the fact that phenylalanine residues occur only in the type II portion of this fragment, many spin systems were localized to either the type I or the type II module via analysis of two-dimensional nuclear-Overhauser-effect (NOESY) experiments. This allowed unambiguous assignment of the two tryptophan residues, as they occur singly in each domain. Patterns of NOESY connectivities are found to be consistent with known type I and type II domain structures; this affords a number of tentative sequence-specific assignments. For both domains, evidence of conserved hydrophobic cores and secondary-structure elements is obtained. In addition, 1H-n.m.r.-monitored thermal-melting studies demonstrate conclusively that the domains are independently folded and that the type I domain has high thermal stability relative to the type II domain. This is consistent with the results of calorimetric studies, and also confirms the localization of spin systems determined from the NOESY data.

Refined solution structure and ligand-binding properties of PDC-109 domain b. A collagen-binding type II domain

We have determined, via 1H-n.m.r., the solution conformation of the collagen-binding b-domain of the bovine seminal fluid protein PDC-109 (PDC-109/b). The structure determination is based on 341 interproton distance estimates and 42 dihedral angle estimates: a set of 24 initial structures were computed; 12 using the variable target function program DIANA, and 12 using the metric matrix program DISGEO. These structures were optimized by restrained energy minimization and dynamic simulated annealing using the CHARMM and X-PLOR programs. The average pairwise root-mean-square difference (r.m.s.d) between the optimized DIANA (DISGEO) structures is 0.71 A (0.82 A) for the backbone atoms, and 1.73 A (2.03 A) for all atoms. Both sets of structures exhibit the same global fold, secondary structure and placement of most non-polar side-chains. Two central antiparallel beta-sheets, which lie roughly perpendicular to each other, and two irregular loops support a large, partially exposed, hydrophobic surface that defines a putative binding site. A test of a hybrid relaxation matrix-based distance refinement protocol (MIDGE program) was performed using a normalized 250 millisecond NOESY spectrum. The resulting distances were input to the molecular mechanics/dynamics procedures mentioned above in order to optimize the DIANA structures. Our results indicate that relaxation matrix refinement of distances is most useful when used conservatively for identifying underestimated distance constraints. 1H-n.m.r. monitored ligand titration experiments revealed definite, albeit weak, binding interactions for phenethylamine and leucine analogs (Ka less than or equal to 25 M-1). Residues perturbed by ligand binding include Tyr7, Trp26, Tyr33, Asp34 and Trp39. These results suggest that PDC-109/b may recognize specific leucine and/or isoleucine-containing sequences within collagen.

Sequence-specific 1H NMR assignments and structural characterization of bovine seminal fluid protein PDC-109 domain b

Sequence-specific resonance assignments for the isolated second or b domain of the bovine seminal fluid protein PDC-109 have been obtained from analysis of two-dimensional 1H NMR experiments recorded at 500 MHz. These assignments include the identification of all aromatic and most aliphatic amino acid resonances. Stereospecific assignment of resonances stemming from the Val2 CH3 gamma,gamma' groups and from seven CH beta,beta' geminal pairs has been accomplished by analysis of 3J alpha beta coupling constants in conjunction with patterns of cross-peak intensities observed in two-dimensional nuclear Overhauser effect (NOESY) spectra. Analysis of NOESY and 3J alpha NH data reveals a small antiparallel beta-sheet involving stretches containing residues 25-28 and 39-42, a cis-proline residue (Pro4), antiparallel strands consisting of residues 1-3, 5-7, and 10-13, and an aromatic cluster composed of Tyr7, Trp26, and Tyr33. The results of distance geometry and restrained molecular dynamics calculations indicate that the global fold of the PDC-109 b domain, a type II module related to those found in fibronectin, is somewhat different from that predicted by modeling the structure on the basis of homology between type II and kringle units. A shallow depression in the molecular surface which presents a solvent-exposed hydrophobic area--a potential ligand-binding site-is identified in the NMR-based models.

The solution conformations of ferrichrome and deferriferrichrome determined by 1H-NMR spectroscopy and computational modeling

We have applied computational procedures that utilize nmr data to model the solution conformation of ferrichrome, a rigid microbial iron transport cyclohexapeptide of known x-ray crystallographic structure [D. van der Helm et al. (1980) J. Am. Chem. Soc. 102, 4224-4231]. The Al3+ and Ga3+ diamagnetic analogues, alumichrome and gallichrome, dissolved in d6-dimethylsulfoxide (d6-DMSO), were investigated via one- and two-dimensional 1H-nmr spectroscopy at 300, 600, and 620 MHz. Interproton distance constraints derived from proton Overhauser experiments were input to a distance geometry algorithm [T. F. Havel and K. Wüthrich (1984) Bull. Math. Biol. 46, 673-691] in order to generate a family of ferrichrome structures consistent with the experimental data. These models were subsequently optimized through restrained molecular dynamics/energy minimization [B. R. Brooks et al. (1983) J. Comp. Chem. 4, 187-217]. The resulting structures were characterized in terms of relative energies and conformational properties. Computations based on integration of the generalized Bloch equations for the complete molecule, which include the 14N-1H dipolar interaction, demonstrate that the x-ray coordinates reproduce the experimental nuclear Overhauser effect time courses very well, and indicate that there are no significant differences between the crystalline and solution conformations of ferrichrome. A similar study of the metal free peptide, deferriferrichrome, suggests that at least two conformers are present in d6-DMSO at 23 degrees C. Both are different from the ferrichrome structure and explain, through conformational averaging, the observed amide NH and CH alpha multiplet splittings. The occurrence of interconverting peptide backbone conformations yields an increased number of sequential NH-CH alpha and NH-NH Overhauser connectivities, which reflects the mean value of r-6 dependence of the dipolar interaction. Our results support the idea that, in the case of structurally rigid peptides, moderately accurate distance constraints define a conformational subspace encompassing the "true" structure, and that energy considerations reduce the size of this subspace. For flexible peptides, however, the straight-forward approach can be misleading since the nmr parameters are averaged over substantially different conformational states.