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.     

The trans labilization of cis-[PtCl2(13CH3NH2)2] by glutathione can be monitored at physiological pH by [1H,13C] HSQC NMR

In order to monitor the trans labilization of cisplatin at physiological pH we have prepared the complex cis-[PtCl₂(¹³CH₃NH₂)₂] and studied its interactions with excess glutathione in aqueous solution at neutral pH by two-dimensional [1H,13C] heteronuclear single-quantum correlation (HSQC) NMR spectroscopy. [1H,13C] HSQC spectroscopy is a good method for following the release of ¹³CH₃NH₂ but is not so good for characterizing the Pt species in solution. In the reaction of cisplatin with glutathione, Pt–S bonds are formed and Pt–NH₃ bonds are broken. The best technique for following the formation of Pt–S bonds of cisplatin is by UV spectroscopy. [1H,13C] HSQC spectroscopy is the best method for following the breaking of the Pt–N bonds. [1H,15N] HSQC spectroscopy is the best method for characterizing the different species in solution. However, the intensity of the peaks in the ¹⁵NH₃–Pt–S region, in [1H,15N] HSQC, reflects a balance between the formation of Pt–S bonds, which increases the signal intensity, and the trans labilization, which decreases the signal intensity. [1H,15N] HSQC spectroscopy and [1H,13C] HSQC spectroscopy are complementary techniques that should be used in conjunction in order to obtain the most accurate information on the interaction of platinum complexes with sulfur-containing ligands.     

Temperature-induced reversible conformational change in the first 100 residues of alpha-synuclein

Natively disordered proteins are a growing class of anomalies to the structure-function paradigm. The natively disordered protein alpha-synuclein is the primary component of Lewy bodies, the cellular hallmark of Parkinson's disease. We noticed a dramatic difference in dilute solution 1H-15N Heteronuclear Single Quantum Coherence (HSQC) spectra of wild-type alpha-synuclein and two disease-related mutants (A30P and A53T), with spectra collected at 35 degrees C showing fewer cross-peaks than spectra acquired at 10 degrees C. Here, we show the change to be the result of a reversible conformational exchange linked to an increase in hydrodynamic radius and secondary structure as the temperature is raised. Combined with analytical ultracentrifugation data showing alpha-synuclein to be monomeric at both temperatures, we conclude that the poor quality of the 1H-15N HSQC spectra obtained at 35 degrees C is due to conformational fluctuations that occur on the proton chemical shift time scale. Using a truncated variant of alpha-synuclein, we show the conformational exchange occurs in the first 100 amino acids of the protein. Our data illustrate a key difference between globular and natively disordered proteins. The properties of globular proteins change little with solution conditions until they denature cooperatively, but the properties of natively disordered proteins can vary dramatically with solution conditions.     

Assignment of cytosine N3 resonances in nucleic acids via intrabase three-bond coupling to amino protons

Coherences were observed between ¹⁵N3 of cytosine and its trans amino proton (H42) using a modified gradient-based heteronuclear single quantum coherence (HSQC) pulse sequence optimized for three-bond proton-nitrogen couplings. The method is demonstrated with a 22-nucleotide RNA fragment of the P5abc region of a group I intron uniformly labeled with ¹⁵N. Use of intraresidue ¹⁵ N3-amino proton couplings to assign cytosine ¹⁵ N3 signals complements the recently proposed J_NN HNN COSY [Dingley, A.J. and Grzesiek, S. (1998) J. Am. Chem. Soc., 120, 8293–8297] method of identifying hydrogen-bonded base pairs in RNA.     

Refocusing revisited: An optimized,gradient-enhanced refocused HSQC and its applications in 2D and 3D NMR and in deuterium exchange experiments

Summary 2D ¹⁵N-¹H correlation spectra are ideal for measuring backbone amide populations to determine amide exchange protection factors in studies of protein folding or other structural features. Most protein NMR spectroscopists use HSQC, which has been shown to be generally superior to HMQC in both resolution and sensitivity. The refocused HSQC experiment is intrinsically less sensitive than the regular HSQC, due to T₂ relaxation during the refocusing delays. However, we show here that, when high ¹⁵N resolution is needed, an optimized refocused HSQC sequence that utilizes a semi-constant time evolution period and pulsed field gradients has better signal-to-noise ratio and resolution, and integrates more accurately, than a similar HSQC. The differences are demonstrated on a 20 kDa protein. The technique can also be applied to 3D NOESY experiments to eliminate strong NH₂ geminal peaks and their truncation artefacts at a modest cost in sensitivity.     

Correlation of nucleotide base and sugar protons in a 15N-labeled HIV-1 RNA oligonucleotide by 1H-15N HSQC experiments

Summary The advent of methods for preparing ¹⁵N- and ¹³C-labeled RNA oligonucleotides holds promise for extending the size of RNA molecules that can be studies by NMR spectroscopy. A practical limitation is the expense of the ¹³C label. It may therefore sometimes be desirable to prepare a relatively inexpensive ¹⁵N-labeled sample only. Here we show that the two-bond ¹H-¹⁵N HSQC experiment can be used on ¹⁵N-labeled RNA to correlate the intranucleotide H1 and H8,H6,H5 resonances indirectly through the shared glycosidic nitrogen. The nonrefocused version of a standard HSQC experiment for 2D proton-detected ¹H-¹⁵N chemical-shift correlation is applied in order to minimize the sensitivity loss due to the relatively fast spin-spin relaxation of RNA oligonucleotides. The experiment is applied to the 30-nucleotide RNA RBE3 which contains the high-affinity binding site of the RRE (rev response element) for the Rev protein of HIV. The results indicate that this simple experiment allows a straightforward identification of the base proton resonances CH5, CH6, UH5, UH6, purine H8, and AH2 as well as the intranucleotide H1 and H8,H6,H5 connectivities. When combined with a NOESY experiment, complete sequential assignments can be obtained.     

Simultaneous CT-13C and VT-15N chemical shift labelling: Application to 3D NOESY-CH3NH and 3D 13C,15N HSQC-NOESY-CH3NH

Based on the HSQC scheme, we have designed a 2D heterocorrelated experiment which combines constant time (CT) ¹³C and variable time (VT) ¹⁵N chemical shift labelling. Although applicable to all carbons, this mode is particularly suitable for simultaneous recording of methyl-carbon and nitrogen chemical shifts at high digital resolution. The methyl carbon magnetisation is in the transverse plane during the whole CT period (1/J_CC=28.6 ms). The magnetisation originating from NH protons is initially stored in the 2H_zN_z state, then prior to the VT chemical shift labelling period is converted into 2H_zN_y coherence. The VT -¹⁵N mode eliminates the effect of ¹ J _N,CO and ^1,2 J _N,CA coupling constants without the need for band-selective carbon pulses. An optional editing procedure is incorporated which eliminates signals from CH₂ groups, thus removing any potential overlap with the CH₃ signals. The CT-¹³CH₃,VT-¹⁵N HSQC building block is used to construct two 3D experiments: 3D NOESY-CH₃NH and 3D ¹³C,¹⁵N HSQC-NOESY-CH₃NH. Combined use of these experiments yields proton and heteronuclear chemical shifts for moieties experiencing NOEs with CH₃ and NH protons. These NOE interactions are resolved as a consequence of the high digital resolution in the carbon and nitrogen chemical shifts of CH₃ and NH groups, respectively. The techniques are illustrated using a double labelled sample of the CH domain from calponin.     

Characterisation of a mobile protein-binding epitope in the translocation domain of colicin E9

The 61 kDa colicin E9 protein toxin enters the cytoplasm of susceptible cells by interacting with outer membrane and periplasmic helper proteins, and kills them by hydrolysing their DNA. The membrane translocation function is located in the N-terminal domain of the colicin, with a key signal sequence being a pentapeptide region that governs the interaction with the helper protein TolB (the TolB box). Previous NMR studies (Collins et al., 2002 J. Mol. Biol. 318, 787-804) have shown that the N-terminal 83 residues of colicin E9, which includes the TolB box, is largely unstructured and highly flexible. In order to further define the properties of this region we have studied a fusion protein containing residues 1-61 of colicin E9 connected to the N-terminus of the E9 DNase by an eight-residue linking sequence. 53 of the expected 58 backbone NH resonances for the first 61 residues and all of the expected 7 backbone NH resonances of the linking sequence were assigned with 3D (1)H-(13)C-(15)N NMR experiments, and the backbone dynamics of these regions investigated through measurement of (1)H-(15)N relaxation properties. Reduced spectral density mapping, extended Lipari-Szabo modelling, and fitting backbone R(2) relaxation rates to a polymer dynamics model identifies three clusters of interacting residues, each containing a tryptophan. Each of these clusters is perturbed by TolB binding to the intact colicin, showing that the significant region for TolB binding extends beyond the recognized five amino acids of the TolB box and demonstrating that the binding epitope for TolB involves a considerable degree of order within an otherwise disordered and flexible domain. Abbreviations : Im9, the immunity protein for colicin E9; E9 DNase, the endonuclease domain of colicin E9; HSQC, heteronuclear single quantum coherence; ppm, parts per million; DSS, 2,2-(dimethylsilyl)propanesulfonic acid; TSP, sodium 3-trimethylsilypropionate; T(1 - 61)-DNase fusion protein, residues 1-61 of colicin E9 connected to the N-terminus of the E9 DNase by an eight residue thrombin cleavage sequence.     

Robust refocusing of 13C magnetization in multidimensional NMR experiments by adiabatic fast passage pulses

We show that adiabatic fast passage (AFP) pulses are robust refocusing elements of transverse ¹³C magnetization in multidimensional NMR experiments. A pair of identical AFP pulses can refocus selected parts or a complete ¹³ C chemical shift range in ¹³C spectra. In the constant time ¹³C-¹H HSQC, replacement of attenuated rectangular pulses by selective AFP pulses results in a sensitivity enhancement of up to a factor of 1.8. In the 3D CBCA(CO)NH the signal-to-noise ratio is increased by a factor of up to 1.6.     

Mapping protein–protein interaction by 13C′-detected heteronuclear NMR spectroscopy

The copper-mediated protein–protein interaction between yeast Atx1 and Ccc2 has been examined by protonless heteronuclear NMR and compared with the already available ¹H–¹⁵N HSQC information. The observed chemical shift variations are analyzed with respect to the actual solution structure, available through intermolecular NOEs. The advantage of using the CON-IPAP spectrum with respect to the ¹H–¹⁵N HSQC resides in the increased number of signals observed, including those of prolines. CBCACO-IPAP experiments allow us to focus on the interaction region and on side-chain carbonyls, while a newly designed CEN-IPAP experiment on side-chains of lysines. An attempt is made to rationalize the chemical shift variations on the basis of the structural data involving the interface between the proteins and the nearby regions. It is here proposed that protonless ¹³C direct-detection NMR is a useful complement to ¹H based NMR spectroscopy for monitoring protein–protein and protein–ligand interactions. Electronic Supplementary Material Supplementary material is available to authorised users in the online version of this article at     

Time-saving methods for heteronuclear multidimensional NMR of (13C, 15N) doubly labeled proteins

Summary Heteronuclear 2D (¹³C, ¹H) and (¹⁵N, ¹H) correlation spectra of (¹³C, ¹⁵N) fully enriched proteins can be acquired simultaneously with virtually no sensitivity loss or increase in artefact levels. Three pulse sequences are described, for 2D time-shared or TS-HSQC, 2D TS-HMQC and 2D TS-HSMQC spectra, respectively. Independent spectral widths can be sampled for both heteronuclei. The sequences can be greatly improved by combining them with field-gradient methods. By applying the sequences to 3D and 4D NMR spectroscopy, considerable time savings can be obtained. The method is demonstrated for the 18 kDa HU protein.Abbreviations HMQC heteronuclear multiple-quantum coherence spectroscopy - HSQC heteronuclear single-quantum coherence spectroscopy - HSMQC heteronuclear single- and multiple-quantum coherence spectroscopy - NOESY nuclear Overhauser enhancement spectroscopy     

Insights into signal transduction revealed by the low resolution structure of the FixJ response regulator

P19(INK4d) is a tumor suppressing protein and belongs to a family of cyclin D-dependent kinase inhibitors of CDK4 and CDK6, which play a key role in human cell cycle control. P19 comprises ten alpha-helices arranged sequentially in five ankyrin repeats forming an elongated structure. This rather simple topology, combined with its physiological function, makes p19 an interesting model protein for folding studies. Urea-induced unfolding transitions monitored by far-UV CD and phenylalanine fluorescence coincide and suggest a two-state mechanism for equilibrium unfolding. Unfolding of p19 followed by 2D (1)H-(15)N HSQC spectra revealed a third species at moderate urea concentrations with a maximum population of about 30 % near 3.2 M urea. It shows poor chemical shift dispersion, but cross-peaks emerge for some residues that are distinct from the native or unfolded state. This equilibrium intermediate either arises only at high protein concentrations (as in the NMR experiment) or has similar optical properties to the unfolded state. Stopped-flow far-UV CD experiments at various urea concentrations revealed that alpha-helical structure is formed in three phases, of which only the fastest phase (10 s(-1)) depends upon the urea concentration. The kinetic of the slowest phase (0.017 s(-1)) can be resolved by 1D real-time NMR and accelerated by cyclophilin. It is limited in rate by prolyl isomerization, and native-like ordered structure cannot form prior to this isomerization. The two fast phases lead to 83 % native protein within the dead time of the NMR experiment. In contrast to p16(INK4a), which exhibits only a marginal stability and high unfolding rates, p19 shows the expected stability for a protein of this size with a clear kinetic barrier between the unfolded and folded state. Therefore, p19 might complement the function of less stable INK4 inhibitors in cell cycle control under unfavorable conditions.     

Complete 1H, 13C and 15N NMR assignments and secondary structure of the 269-residue serine protease PB92 from Bacillus alcalophilus

Summary The ¹H, ¹³C and ¹⁵N NMR resonances of serine protease PB92 have been assigned using 3D tripleresonance NMR techniques. With a molecular weight of 27 kDa (269 residues) this protein is one of the largest monomeric proteins assigned so far. The side-chain assignments were based mainly on 3D H(C)CH and 3D (H)CCH COSY and TOCSY experiments. The set of assignments encompasses all backbone carbonyl and CH_n carbons, all amide (NH and NH₂) nitrogens and 99.2% of the amide and CH_n protons. The secondary structure and general topology appear to be identical to those found in the crystal structure of serine protease PB92 [Van der Laan et al. (1992) Protein Eng., 5, 405–411], as judged by chemical shift deviations from random coil values, NH exchange data and analysis of NOEs between backbone NH groups.Abbreviations 2D/3D/4D two-/three-/four-dimensional - HSQC heteronuclear single-quantum coherence - HMQC heteronuclear multiple-quantum coherence - COSY correlation spectroscopy - TOCSY total correlation spectroscopy - NOE nuclear Overhauser enhancement (connectivity) - NOESY 2D NOE spectroscopy Experiment nomenclature (H(C)CH, etc.) follows the conventions used elsewhere [e.g. Ikura et al. (1990) Biochemistry, 29, 4659–4667].     

A1H-15N NM R study of human c-Ha-ras protein: Biosynthetic incorporation of15N-labeled amino acids

Summary A¹H-¹⁵N NMR study was performed on the GDP-bound form of a truncated human c-Ha-ras oncogene product (171 amino acid residues). Resonance cross peaks of the backbone amide¹H-¹⁵N nuclei of a uniformly¹⁵N-labeled protein were observed with heteronuclear-single-quantum coherence spectroscopy (HSQC). In order to resolve overlapping cross peaks, selective¹⁵N-labeling of one or two types of amino acid residues (Ala, Arg, Asx, Glx, Gly, His, lie, Leu, Lys, Met, Phe, Ser; Thr, Tyr and/or Val) was carried out using appropriateE. coli mutant strains. By this procedure, all the backbone¹H-¹⁵N cross peaks were classified into amino acid types.This work was supported in part by a Grant-in-Aid for Cancer Research from the Ministry of Education, Science and Culture, Japan     

Resonance assignment and structural analysis of acid denatured E. coli [U-15N]-glutaredoxin 3

Virtually complete sequence specific ¹H and ¹⁵N resonance assignments are presented for acid denatured reduced E. coli glutaredoxin 3. The sequential resonance assignments of the backbone rely on the combined use of 3D F₁-decoupled ROESY-¹⁵N-HSQC and 3D ¹⁵N-HSQC-(TOCSY-NOESY)-¹⁵N-HSQC using a single uniformly ¹⁵N labelled protein sample. The sidechain resonances were assigned from a 3D TOCSY-¹⁵N-HSQC and a homonouclear TOCSY spectrum. The presented assignment strategy works in the absence of chemical exchange peaks with signals from the native conformation and without ¹³C/¹⁵N double labelling. Chemical shifts, ³J(H, NH) coupling constants and NOEs indicate extensive conformational averaging of both backbone and side chains in agreement with a random coil conformation. The only secondary structure element persisting at pH 3.5 appears to be a short helical segment comprising residues 37 to 40.Abbreviations HSQC heteronuclear single quantum coherence - NMR nuclear magnetic resonance - NOE nuclear Overhauser effect - NOESY two-dimensional NOE spectroscopy - ROE nuclear Overhauser effect in the rotating frame - ROESY two-dimensional ROE spectroscopy - TOCSY total correlation spectroscopy - TPPI time proportional phase incrementation Correspondence to: G. Otting     

Dictyostelium discoideum as Expression Host: Isotopic Labeling of a Recombinant Glycoprotein for NMR Studies

The advantages of the organism Dictyostelium discoideum as an expression host for recombinant glycoproteins have been exploited for the production of an isotopically labeled cell surface protein for NMR structure studies. Growth medium containing [¹⁵N]NH₄Cl and [¹³C]glycerol was used to generate isotopically labeled Escherichia coli, which was subsequently introduced to D. discoideum cells in simple Mes buffer. A variety of growth conditions were screened to establish minimal amounts of nitrogen and carbon metabolites for a cost-effective protocol. Following single-step purification by anion-exchange chromatography, 8 mg of uniformly ¹³C,¹⁵N-labeled protein secreted by approximately 10¹⁰D. discoideum cells was isolated from 3.3 liters of supernatant. Mass spectrometry showed the recombinant protein of 16 kDa to have incorporated greater than 99.9% isotopic label. The two-dimensional ¹H-¹³C HSQC spectrum confirms ¹³C labeling of both glycan and amino acid residues of the glycoprotein. All heteronuclear NMR spectra showed a good dispersion of cross-peaks essential for high-quality structure determination.     

Synthesis and structural analysis of five novel oligosaccharides prepared by glucosyltransfer from beta-D-glucose 1-phosphate to isokestose and nystose using Thermoanaerobacter brockii kojibiose phosphorylase

Five novel oligosaccharides (tetra-, penta- and hexa-saccharides) were synthesized by glucosyltransfer from beta-D-glucose 1-phosphate to isokestose (O-beta-D-fructofuranosyl-(2-->1)-O-beta-D-fructofuranosyl-(2-->1)-alpha-D-glucopyranoside) or nystose (O-beta-D-fructofuranosyl-(2-->1)-O-beta-D-fructofuranosyl-(2-->1)-O-beta-D-fructofuranosyl-(2-->1)-alpha-D-glucopyranoside) using Thermoanaerobacter brockii kojibiose phosphorylase. The oligosaccharides were identified as 2(2-alpha-D-glucopyranosyl)(m)isokestose; [O-alpha-D-glucopyranosyl-(1-->2)](m)-O-[beta-D-fructofuranosyl-(2-->1)](2)-alpha-D-glucopyranoside: m=1, 2, and 3, and 2(2-alpha-D-glucopyranosyl)(n)nystose; [O-alpha-D-glucopyranosyl-(1-->2)](n)-O-[beta-D-fructofuranosyl-(2-->1)](3)-alpha-D-glucopyranoside: n=1 and 2 using gas liquid chromatography analysis of the methyl derivatives, and MALDI-TOF-MS and NMR measurements of the newly formed oligosaccharides. 1H, 13C NMR signals of each saccharide were assigned using 2D-NMR techniques, including COSY, HSQC, HSQC-TOCSY, HMBC, CH(2)-selected E-HSQC, and CH(2)-selected E-HSQC-TOCSY.     

Plant and soil natural abundance <Emphasis Type="Italic">δ</Emphasis><Superscript>15</Superscript>N: indicators of relative rates of nitrogen cycling in temperate forest ecosystems

Watersheds within the Catskill Mountains, New York, receive among the highest rates of nitrogen (N) deposition in the northeastern United States and are beginning to show signs of N saturation. Despite similar amounts of N deposition across watersheds within the Catskill Mountains, rates of soil N cycling and N retention vary significantly among stands of different tree species. We examined the potential use of δ ¹⁵N of plants and soils as an indicator of relative forest soil N cycling rates. We analyzed the δ ¹⁵N of foliage, litterfall, bole wood, surface litter layer, fine roots and organic soil from single-species stands of American beech (Fagus grandifolia), eastern hemlock (Tsuga canadensis), red oak (Quercus rubra), and sugar maple (Acer saccharum). Fine root and organic soil δ ¹⁵N values were highest within sugar maple stands, which correlated significantly with higher rates of net mineralization and nitrification. Results from this study suggest that fine root and organic soil δ ¹⁵N can be used as an indicator of relative rates of soil N cycling. Although not statistically significant, δ ¹⁵N was highest within foliage, wood and litterfall of beech stands, a tree species associated with intermediate levels of soil N cycling rates and forest N retention. Our results show that belowground δ ¹⁵N values are a better indicator of relative rates of soil N cycling than are aboveground δ ¹⁵N values.     

Influence of light intensity on the distribution of carbon and consequent effects on mineralization of soil nitrogen in a barley (Hordeum vulgare L.)-soil system

A pot experiment was conducted in a ¹⁴C-labelled atmosphere to study the influence of living plants on organic-N mineralization. The soil organic matter had been labelled, by means of a 200-days incubation, with ¹⁵N. The influence of the carbon input from the roots on the formation of microbial biomass was evaluated by using two different light intensities (I). Mineralization of ¹⁵N-labelled soil N was examined by following its fate in both the soil biomass and the plants. Less dry matter accumulated in shoots and roots at the lower light intensity. Furthermore, in all the plant-soil compartments examined, with the exception of rhizosphere respiration, the proportion of net assimilated ¹⁴C was lower in the low-I treatment than in the high-I treatment. The lower rates of ¹⁴C and ¹⁵N incorporation into the soil biomass were associated with less root-derived ¹⁴C. During the chamber period (¹⁴CO₂-atmosphere), mineralized amounts of ¹⁵N (measured as plant uptake of ¹⁵N) were small and represented about 6.8 to 7.8% of the initial amount of organic ¹⁵N in the soil. Amounts of unlabelled N found in the plants, as a percentage of total soil N, were 2.5 to 3.3%. The low availability of labelled N to microorganisms was the result of its stabilization during the 210 days of soil incubation. Differences in carbon supply resulted in different rates of N mineralization which is consistent with the hypothesis that roots induce N mineralization. N mineralization was higher in the high-I treatment. On the other hand, the rate of mineralization of unlabelled stable soil N was lower than labelled soil ¹⁵N which was stabilized. The amounts of ¹⁵N mineralized in planted soil during the chamber period (43 days) which were comparable with those mineralized in unplanted soil incubated for 210 days, also suggested that living plants increased the turnover rate of soil organic matter.