Apoflavodoxin (un)folding followed at the residue level by NMR

The denaturant-induced (un)folding of apoflavodoxin from Azotobacter vinelandii has been followed at the residue level by NMR spectroscopy. NH groups of 21 residues of the protein could be followed in a series of 1H-15N heteronuclear single-quantum coherence spectra recorded at increasing concentrations of guanidinium hydrochloride despite the formation of protein aggregate. These NH groups are distributed throughout the whole apoflavodoxin structure. The midpoints of unfolding determined by NMR coincide with the one obtained by fluorescence emission spectroscopy. Both techniques give rise to unfolding curves with transition zones at significantly lower denaturant concentrations than the one obtained by circular dichroism spectroscopy. The NMR (un)folding data support a mechanism for apoflavodoxin folding in which a relatively stable intermediate is involved. Native apoflavodoxin is shown to cooperatively unfold to a molten globule-like state with extremely broadened NMR resonances. This initial unfolding step is slow on the NMR chemical shift timescale. The subsequent unfolding of the molten globule is faster on the NMR chemical shift timescale and the limited appearance of 1H-15N HSQC cross peaks of unfolded apoflavodoxin in the denaturant range studied indicates that it is noncooperative.     

An NMR investigation of solution aggregation reactions preceding the misassembly of acid-denatured cold shock protein A into fibrils.

At pH 2.0, acid-denatured CspA undergoes a slow self-assembly process, which results in the formation of insoluble fibrils. 1H-15N HSQC, 3D HSQC-NOESY, and 15N T2 NMR experiments have been used to characterize the soluble components of this reaction. The kinetics of self-assembly show a lag phase followed by an exponential increase in polymerization. A single set of 1H-15N HSQC cross-peaks, corresponding to acid-denatured monomers, is observed during the entire course of the reaction. Under lag phase conditions, 15N resonances of residues that constitute the beta-strands of native CspA are selectively broadened with increasing protein concentration. The dependence of 15N T2 values on spin echo period duration demonstrates that line broadening is due to fast NMR exchange between acid-denatured monomers and soluble aggregates. Exchange contributions to T2 relaxation correlate with the squares of the chemical shift differences between native and acid-denatured CspA, and point to a stabilization of native-like structure upon aggregation. Time-dependent changes in 15N T2 relaxation accompanying the exponential phase of polymerization suggest that the first three beta-strands may be predominantly responsible for association interfaces that promote aggregate growth. CspA serves as a useful model system for exploring the conformational determinants of denatured protein misassembly.     

Altered active site flexibility and a structural metal-binding site in eukaryotic dUTPase: kinetic characterization, folding, and crystallographic studies of the homotrimeric Drosophila enzyme

dUTPase is responsible for preventive DNA repair via exclusion of uracil. Developmental regulation of the Drosophila enzyme is suggested to be involved in thymine-less apoptosis. Here we show that in addition to conserved dUTPase sequence motifs, the gene of Drosophila enzyme codes for a unique Ala-Pro-rich segment. Kinetic and structural analyses of the recombinant protein and a truncation mutant show that the Ala-Pro segment is flexible and has no regulatory role in vitro. The homotrimer enzyme unfolds reversibly as a trimeric entity with a melting temperature of 54 degrees C, 23 degrees C lower than Escherichia coli dUTPase. In contrast to the bacterial enzyme, Mg(2+) binding modulates conformation of fly dUTPase, as identified by spectroscopy and by increment in melting temperature. A single well folded, but inactive, homotrimeric core domain is generated through three distinct steps of limited trypsinolysis. In fly, but not in bacterial dUTPase, binding of the product dUMP induces protection against proteolysis at the tryptic site reflecting formation of the catalytically competent closed conformer. Crystallographic analysis argues for the presence of a stable monomer of Drosophila dUTPase in crystal phase. The significant differences between prototypes of eukaryotic and prokaryotic dUTPases with respect to conformational flexibility of the active site, substrate specificity, metal ion binding, and oligomerization in the crystal phase are consistent with alteration of the catalytic mechanism and hydropathy of subunit interfaces.     

Polypeptide backbone resonance assignments and secondary structure of Bacillus subtilis enzyme IIIglc determined by two-dimensional and three-dimensional heteronuclear NMR spectroscopy

The enzyme IIIglc-like domain of Bacillus subtilis IIglc (IIIglc, 162 residues, 17.4 kDa) has been cloned and overexpressed in Escherichia coli. Sequence-specific assignment of the backbone 1H and 15N resonances has been carried out with a combination of homonuclear and heteronuclear two-dimensional and heteronuclear three-dimensional (3D) NMR spectroscopy. Amide proton solvent exchange rate constants have been determined from a series of 1H-15N heteronuclear single-quantum coherence (HSQC) spectra acquired following dissolution of the protein in D2O. Major structural features of IIIglc have been inferred from the pattern of short-, medium- and long-range NOEs in 3D heteronuclear 1H nuclear Overhauser effect 1H-15N multiple-quantum coherence (3D NOESY-HMQC) spectra, together with the exchange rate constants. IIIglc contains three antiparallel beta-sheets comprised of eight, three, and two beta-strands. In addition, five beta-bulges were identified. No evidence of regular helical structure was found. The N-terminal 15 residues of the protein appear disordered, which is consistent with their being part of the Q-linker that connects the C-terminal enzyme IIIglc-like domain to the membrane-bound IIglc domain. Significantly, two histidine residues, His 68 and His 83, which are important for phosphotransferase function, are found from NOE measurements to be in close proximity at the ends of adjacent strands in the major beta-sheet.     

Structural characterization of an intrinsically unfolded mini-HBX protein from hepatitis B virus

The hepatitis B virus x protein (HBX) is expressed in HBVinfected liver cells and can interact with a wide range of cellular proteins. In order to understand such promiscuous behavior of HBX we expressed a truncated mini-HBX protein (named Tr-HBX) (residues 18-142) with 5 Cys → Ser mutations and characterized its structural features using circular dichroism (CD) spectropolarimetry, NMR spectroscopy as well as bioinformatics tools for predicting disorder in intrinsically unstructured proteins (IUPs). The secondary structural content of Tr-HBX from CD data suggests that Tr-HBX is only partially folded. The protein disorder prediction by IUPred reveals that the unstructured region encompasses its N-terminal ～30 residues of Tr-HBX. A two-dimensional (1)H-(15)N HSQC NMR spectrum exhibits fewer number of resonances than expected, suggesting that Tr-HBX is a hybrid type IUP where its folded C-terminal half coexists with a disordered N-terminal region. Many IUPs are known to be capable of having promiscuous interactions with a multitude of target proteins. Therefore the intrinsically disordered nature of Tr-HBX revealed in this study provides a partial structural basis for the promiscuous structure-function behavior of HBX.     

Comparison of cell-based and cell-free protocols for producing target proteins from the Arabidopsis thaliana genome for structural studies

We describe a comparative study of protein production from 96 Arabidopsis thaliana open reading frames (ORFs) by cell-based and cell-free protocols. Each target was carried through four pipeline protocols used by the Center for Eukaryotic Structural Genomics (CESG), one for the production of unlabeled protein to be used in crystallization trials and three for the production of 15N-labeled proteins to be analyzed by 1H-15N NMR correlation spectroscopy. Two of the protocols involved Escherichia coli cell-based and two involved wheat germ cell-free technology. The progress of each target through each of the protocols was followed with all failures and successes noted. Failures were of the following types: ORF not cloned, protein not expressed, low protein yield, no cleavage of fusion protein, insoluble protein, protein not purified, NMR sample too dilute. Those targets that reached the goal of analysis by 1H-15N NMR correlation spectroscopy were scored as HSQC+ (protein folded and suitable for NMR structural analysis), HSQC+/- (protein partially disordered or not in a single stable conformational state), HSQC- (protein unfolded, misfolded, or aggregated and thus unsuitable for NMR structural analysis). Targets were also scored as X- for failing to crystallize and X+ for successful crystallization. The results constitute a rich database for understanding differences between targets and protocols. In general, the wheat germ cell-free platform offers the advantage of greater genome coverage for NMR-based structural proteomics whereas the E. coli platform when successful yields more protein, as currently needed for crystallization trials for X-ray structure determination.     

Biophysical characterization of structural properties and folding of interleukin-1 receptor antagonist

Structural properties and folding of interleukin-1 receptor antagonist (IL-1ra), a therapeutically important cytokine with a symmetric beta-trefoil topology, are characterized using optical spectroscopy, high-resolution NMR, and size-exclusion chromatography. Spectral contributions of two tryptophan residues, Trp17 and Trp120, present in the wild-type protein, have been determined from mutational analysis. Trp17 dominates the emission spectrum of IL-1ra, while Trp120 is quenched presumably by the nearby cysteine residues in both folded and unfolded states. The same Trp17 gives rise to two characteristic negative peaks in the aromatic CD. Urea denaturation of the wild-type protein is probed by measuring intrinsic and extrinsic (binding of 1-anilinonaphthalene-8-sulfonic acid) fluorescence, near- and far-UV CD, and 1D and 2D ((1)H-(15)N heteronuclear single quantum coherence (HSQC)) NMR. Overall, the data suggest an essentially two-state equilibrium denaturation mechanism with small, but detectable structural changes within the pretransition region. The majority of the (1)H-(15)N HSQC cross-peaks of the folded state show only a limited chemical shift change as a function of the denaturant concentration. However, the amide cross-peak of Leu31 demonstrates a significant urea dependence that can be fitted to a two-state binding model with a dissociation constant of 0.95+/-0.04 M. This interaction has at least a five times higher affinity than reported values for nonspecific urea binding to denatured proteins and peptides, suggesting that the structural context around Leu31 stabilizes the protein-urea interaction. A possible role of denaturant binding in inducing the pretransition changes in IL-1ra is discussed. Urea unfolding of wild-type IL-1ra is sufficiently slow to enable HPLC separation of folded and unfolded states. Quantitative size-exclusion chromatography has provided a hydrodynamic view of the kinetic denaturation process. Thermodynamic stability and unfolding kinetics of IL-1ra resemble those of structurally and evolutionary close IL-1beta, suggesting similarity of their free energy landscapes.     

Mechanistic studies of dUTPases

The essential enzyme dUTPase is responsible for preventive DNA repair via exclusion of uracil. Lack or inhibition of the enzyme induces thymine-less cell death in cells performing active DNA synthesis, serving therefore as an important chemotherapeutic target. In the present work, employing differential circular dichroism spectroscopy, we show that D. mel. dUTPase, a recently described eukaryotic model, has a similar affinity of binding towards alpha,beta-imino-dUTP as compared to the prokaryotic E. coli enzyme. However, in contrast to the prokaryotic dUTPase, the nucleotide exerts significant protection against tryptic digestion at a specific tryptic site 20 A far from the active site in the fly enzyme. This result indicates that binding of the nucleotide in the active site induces an allosteric conformational change within the central threefold channel of the homotrimer exclusively in the eukaryotic enzyme. Nucleotide binding induced allosterism in the D. mel. dUTPase, but not in the E. coli enzyme, might be associated with the altered hydropathy of subunit interfaces in these two proteins.     

Backbone 1H and 15N resonance assignments of the N-terminal SH3 domain of drk in folded and unfolded states using enhanced-sensitivity pulsed field gradient NMR techniques

Summary The backbone ¹H and ¹⁵N resonances of the N-terminal SH3 domain of the Drosophila signaling adapter protein, drk, have been assigned. This domain is in slow exchange on the NMR timescale between folded and predominantly unfolded states. Data were collected on both states simultaneously, on samples of the SH3 in near physiological buffer exhibiting an approximately 1:1 ratio of the two states. NMR methods which exploit the chemical shift dispersion of the ¹⁵N resonances of unfolded states and pulsed field gradient water suppression approaches for avoiding saturation and dephasing of amide protons which rapidly exchange with solvent were utilized for the assignment.Abbreviations 2D, 3D two-, three-dimensional - drkN SH3 N-terminal SH3 domain of Drosophila drk - HSQC heteronuclear single-quantum spectroscopy - NOE nuclear Overhauser enhancement - SH3 Src homology domain 3 - TOCSY total correlation spectroscopy     

Characterization of the structure and dynamics of a near-native equilibrium intermediate in the unfolding pathway of an all beta-barrel protein

The structure and dynamics of equilibrium intermediate in the unfolding pathway of the human acidic fibroblast growth factor (hFGF-1) are investigated using a variety of biophysical techniques including multidimensional NMR spectroscopy. Guanidinium hydrochloride (GdnHCl)-induced unfolding of hFGF-1 proceeds with the accumulation of a stable intermediate state. The transition from the intermediate state to the unfolded state(s) is cooperative without the accumulation of additional intermediate(s). The intermediate state induced maximally in 0.96 m GdnHCl is found to be obligatory in the folding/unfolding pathway of hFGF-1. Most of the native tertiary structure interactions are preserved in the intermediate state. (1)H-(15)N chemical shift perturbation data suggest that the residues in the C-terminal segment including those located in the beta-strands IX, X, and XI undergo the most discernible structural change(s) in the intermediate state in 0.96 m GdnHCl. hFGF-1 in the intermediate state (0.96 m GdnHCl) does not bind to its ligand, sucrose octasulfate. Limited proteolytic digestion experiments and hydrogen-deuterium exchange monitored by (15)N heteronuclear single quantum coherence (HSQC) spectra show that the conformational flexibility of the protein in the intermediate state is significantly higher than in the native conformation. (15)N spin relaxation experiments show that many residues located in beta-strands IX, X, and XI exhibit conformational motions in the micro- to millisecond time scale. Analysis of (15)N relaxation data in conjunction with the amide proton exchange kinetics suggests that residues in the beta-strands II, VIII, and XII possibly constitute the stability core of the protein in the near-native intermediate state.     

A preliminary CD and NMR study of the Escherichia coli DNA polymerase III theta subunit.

The theta subunit of DNA polymerase III, the main replicative polymerase of Escherichia coli, has been examined by circular dichroism and by NMR spectroscopy. The polymerase core consists of three subunits: alpha, epsilon, and theta, with alpha possessing the polymerase activity, epsilon functioning as a proofreading exonuclease, and theta, a small subunit of 8.9 kD, of undetermined function. The theta subunit has been expressed in E. coli, and a CD analysis of theta indicates the presence of a significant amount of secondary structure: approximately 52% alpha helix, 9% beta sheet, 21% turns, and 18% random coil. However, at higher concentrations, theta yields a poorly-resolved 1D proton NMR spectrum in which both the amide protons and the methyl protons show poor chemical shift dispersion. Subsequent 1H-15N HSQC analysis of uniformly-15N-labeled theta supports the conclusion that approximately half of the protein is reasonably well-structured. Another quarter of the protein, probably including some of the N-terminal region, is highly mobile, exhibiting a chemical shift pattern indicative of random coil structure. The remaining amide resonances exhibit significant broadening, indicative of intermolecular and/or intramolecular exchange processes. Improved chemical shift dispersion and greater uniformity of resonance intensities in the 1H-15N HSQC spectra resulted when [U-15N]-theta was examined in the presence of epsilon186--the N-terminal domain of the epsilon-subunit. Further work is currently in progress to define the solution structure of theta and the theta-epsilon186 complex.     

Urea-dependent unfolding of murine adenosine deaminase: sequential destabilization as measured by 19F NMR

Murine adenosine deaminase (mADA) is a 40 kDa (beta/alpha)(8)-barrel protein consisting of eight central beta-strands and eight peripheral alpha-helices containing four tryptophan residues. In this study, we investigated the urea-dependent behavior of the protein labeled with 6-fluorotryptophan (6-(19)F-Trp). The (19)F NMR spectrum of 6-(19)F-Trp-labeled mADA reveals four distinct resonances in the native state and three partly overlapped resonances in the unfolded state. The resonances were assigned unambiguously by site-directed mutagenesis. Equilibrium unfolding of 6-(19)F-Trp-labeled mADA was monitored using (19)F NMR based on these assignments. The changes in intensity of folded and unfolded resonances as a function of urea concentration show transition midpoints consistent with data observed by far-UV CD and fluorescence spectroscopy, indicating that conformational changes in mADA during urea unfolding can be followed by (19)F NMR. Chemical shifts of the (19)F resonances exhibited different changes between 1.0 and 6.0 M urea, indicating that local structures around 6-(19)F-Trp residues change differently. The urea-induced changes in local structure around four 6-(19)F-Trp residues of mADA were analyzed on the basis of the tertiary structure and chemical shifts of folded resonances. The results reveal that different local regions in mADA have different urea-dependent behavior, and that local regions of mADA change sequentially from native to intermediate topologies on the unfolding pathway.     

Investigation of ligand binding and protein dynamics in Bacillus subtilis chorismate mutase by transverse relaxation optimized spectroscopy-nuclear magnetic resonance

The structural and dynamical consequences of ligand binding to a monofunctional chorismate mutase from Bacillus subtilis have been investigated by solution NMR spectroscopy. TROSY methods were employed to assign 98% of the backbone (1)H(N), (1)H(alpha), (15)N, (13)C', and (13)C(alpha) resonances as well as 86% of the side chain (13)C resonances of the 44 kDa trimeric enzyme at 20 degrees C. This information was used to map chemical shift perturbations and changes in intramolecular mobility caused by binding of prephenate or a transition state analogue to the X-ray structure. Model-free interpretation of backbone dynamics for the free enzyme and its complexes based on (15)N relaxation data measured at 600 and 900 MHz showed significant structural consolidation of the protein in the presence of a bound ligand. In agreement with earlier structural and biochemical studies, substantial ordering of 10 otherwise highly flexible residues at the C-terminus is particularly notable. The observed changes suggest direct contact between this protein segment and the bound ligand, providing support for the proposal that the C-terminus can serve as a lid for the active site, limiting diffusion into and out of the pocket and possibly imposing conformational control over substrate once bound. Other regions of the protein that experience substantial ligand-induced changes also border the active site or lie along the subunit interfaces, indicating that the enzyme adapts dynamically to ligands by a sort of induced fit mechanism. It is believed that the mutase-catalyzed chorismate-to-prephenate rearrangement is partially encounter controlled, and backbone motions on the millisecond time scale, as seen here, may contribute to the reaction barrier.     

Compact molten globule-like state of hUBF HMG Box1 at extremely low pH

Using far and near-UV CD, ANS fluorescence and 2D NMR spectroscopy, an acid-induced partly folded state (A state) at extremely low pH for hUBF HMG Box1 was identified and characterized. As compared to the native state (N), the A state has similar secondary structure, less compact pack with larger amounts of exposed hydrophobic surface, and narrower chemical shift dispersion in (1)H-(15)N HSQC spectrum, which implies that it is a molten globule (MG)-like species. On the other hand, substantial tertiary contacts and cooperative thermal denaturing transition indicate that the A state is closer-relative to the classic MG-to the native folded state. In addition, when the solution pH is adjusted to neutrality, the protein in the A state refolds to the native state easily. All these data suggest that the A state of hUBF HMG Box1 could represent a potential folding intermediate on protein folding pathway.     

Assignment of 1H NMR resonances of histidine and other aromatic residues in met-, cyano-, oxy-, and (carbon monoxy)myoglobins

The resolved 1H NMR resonances of the aromatic region in the 270-MHz NMR spectrum of sperm whale, horse, and pig metmyoglobin (metMb) have been assigned, including the observable H-2 and H-4 histidine resonances, the tryptophan H-2 resonances, and upfield-shifted resonances from one tyrosine residue. The use of different Mb species, carboxymethylation, and matching of pK values allows the assignment of the H-4 resonances, which agree in only three cases out of seven with scalar-correlated two-dimensional NMR spectroscopy assignments by others. The conversion to hydroxymyoglobin at high pH involves rearrangements throughout the molecule and is observed by many assigned residues. In sperm whale ferric cyanomyoglobin, nine H-2 and eight H-4 histidine resonances have been assigned, including the His-97 H-2 resonance and tyrosine resonances from residues 103 and 146. The hyperfine-shifted resonances from heme and near-heme protons observe a shift with a pK = 5.3 +/- 0.3 (probably due to deprotonation of His-97, pK = 5.6) and another shift at pK = 10.8 +/- 0.3. The spectrum of high-spin ferrous sperm whale deoxymyoglobin is very similar to that of metMb, which allows the assignment of seven surface histidine H-2 and H-4 resonances and also resonances from the two tryptophan residues and one tyrosine. In diamagnetic sperm whale (carbon monoxy)myoglobin (COMb), 10 His H-2 and 11 His H-4 resonances are observed, and 8 H-2 and 9 H-4 resonances are assigned, including His-64 H-4, the distal histidine. This important resonance is not observed in sperm whale oxymyoglobin, which in general shows very similar titration curves to COMb. Histidine-36 shows unusual titration behavior in the paramagnetic derivatives but normal behavior in the diamagnetic derivatives, which is discussed in the accompanying paper [Bradbury, J. H., & Carver, J. A. (1984) Biochemistry (following paper in this issue)].     

Mechanistic Studies of dUTPases

The essential enzyme dUTPase is responsible for preventive DNA repair via exclusion of uracil. Lack or inhibition of the enzyme induces thymine‐less cell death in cells performing active DNA synthesis, serving therefore as an important chemotherapeutic target. In the present work, employing differential circular dichroism spectroscopy, we show that D. mel. dUTPase, a recently described eukaryotic model, has a similar affinity of binding towards α,β‐imino‐dUTP as compared to the prokaryotic E. coli enzyme. However, in contrast to the prokaryotic dUTPase, the nucleotide exerts significant protection against tryptic digestion at a specific tryptic site 20 Å far from the active site in the fly enzyme. This result indicates that binding of the nucleotide in the active site induces an allosteric conformational change within the central threefold channel of the homotrimer exclusively in the eukaryotic enzyme. Nucleotide binding induced allosterism in the D. mel. dUTPase, but not in the E. coli enzyme, might be associated with the altered hydropathy of subunit interfaces in these two proteins.     

Backbone nuclear magnetic resonance assignment of human deoxyuridine 5′-triphosphate nucleotidohydrolase (dUTPase)

Nuclear-associated deoxyuridine 5′-triphosphate nucleotidohydrolase (dUTPase) is an enzyme that hydrolyses deoxyuridine 5′-triphosphate (dUTP) to the monophosphate, thereby controlling the dUTP levels of the organism, which is essential for survival. Further, dUTPase is up-regulated in many cancers. Thus, dUTPase is a highly interesting potential drug target. We report, for the first time, the near complete nuclear magnetic resonance (NMR) spectroscopy ¹⁵N/¹³C/¹H backbone assignment of the 3 × 164 amino acids homo-trimer human dUTPase. Previously, only a handful backbone resonances belonging to the flexible C-terminus has been published for any protein in the dUTPase family.     

Structural dynamics of the membrane translocation domain of colicin E9 and its interaction with TolB

In order for the 61 kDa colicin E9 protein toxin to enter the cytoplasm of susceptible cells and kill them by hydrolysing their DNA, the colicin must interact with the outer membrane BtuB receptor and Tol translocation pathway of target cells. The translocation function is located in the N-terminal domain of the colicin molecule. (1)H, (1)H-(1)H-(15)N and (1)H-(13)C-(15)N NMR studies of intact colicin E9, its DNase domain, minimal receptor-binding domain and two N-terminal constructs containing the translocation domain showed that the region of the translocation domain that governs the interaction of colicin E9 with TolB is largely unstructured and highly flexible. Of the expected 80 backbone NH resonances of the first 83 residues of intact colicin E9, 61 were identified, with 43 of them being assigned specifically. The absence of secondary structure for these was shown through chemical shift analyses and the lack of long-range NOEs in (1)H-(1)H-(15)N NOESY spectra (tau(m)=200 ms). The enhanced flexibility of the region of the translocation domain containing the TolB box compared to the overall tumbling rate of the protein was identified from the relatively large values of backbone and tryptophan indole (15)N spin-spin relaxation times, and from the negative (1)H-(15)N NOEs of the backbone NH resonances. Variable flexibility of the N-terminal region was revealed by the (15)N T(1)/T(2) ratios, which showed that the C-terminal end of the TolB box and the region immediately following it was motionally constrained compared to other parts of the N terminus. This, together with the observation of inter-residue NOEs involving Ile54, indicated that there was some structural ordering, resulting most probably from the interactions of side-chains. Conformational heterogeneity of parts of the translocation domain was evident from a multiplicity of signals for some of the residues. Im9 binding to colicin E9 had no effect on the chemical shifts or other NMR characteristics of the region of colicin E9 containing the TolB recognition sequence, though the interaction of TolB with intact colicin E9 bound to Im9 did affect resonances from this region. The flexibility of the translocation domain of colicin E9 may be connected with its need to recognise protein partners that assist it in crossing the outer membrane and in the translocation event itself.     

Automated evaluation of chemical shift perturbation spectra: New approaches to quantitative analysis of receptor-ligand interaction NMR spectra

This paper presents new methods designed for quantitative analysis of chemical shift perturbation NMR spectra. The methods automatically trace the displacements of cross peaks between a perturbed test spectrum and the reference spectrum (or among a series of titration spectra), and measure the changes of chemical shifts, heights, and widths of the altered peaks. The methods are primary aimed at the (1)H-(15)N HSQC spectra of relatively small proteins (<15 kDa) assuming fast exchange between free and ligand-bound states on the chemical shift time scale, or for comparing spectra of free and fully bound states in the slow exchange situation. Using the (1)H-(15)N HSQC spectra from a titration experiment of the 74-residue Pex13p SH3 domain with a Pex14p peptide ligand (14 residues, K (d)= approximately 40 microM), we demonstrate the scope and limits of our automatic peak tracing (APET) algorithm for efficient scoring of high-throughput SAR by NMR type HSQC spectra, and progressive peak tracing (PROPET) algorithm for detailed analysis of ligand titration spectra. Simulated spectra with low signal-to-noise ratios (S/N ranged from 20 to 1) were used to demonstrate the reliability and reproducibility of the results when dealing with poor quality spectra. These algorithms have been implemented in a new software module, FELIX-Autoscreen, for streamlined processing, analysis and visualization of SAR by NMR and other high-throughput receptor/ligand interaction experiments.