Comprehensive analysis of NMR data using advanced line shape fitting

NMR spectroscopy is uniquely suited for atomic resolution studies of biomolecules such as proteins, nucleic acids and metabolites, since detailed information on structure and dynamics are encoded in positions and line shapes of peaks in NMR spectra. Unfortunately, accurate determination of these parameters is often complicated and time consuming, in part due to the need for different software at the various analysis steps and for validating the results. Here, we present an integrated, cross-platform and open-source software that is significantly more versatile than the typical line shape fitting application. The software is a completely redesigned version of PINT (https://pint-nmr.github.io/PINT/). It features a graphical user interface and includes functionality for peak picking, editing of peak lists and line shape fitting. In addition, the obtained peak intensities can be used directly to extract, for instance, relaxation rates, heteronuclear NOE values and exchange parameters. In contrast to most available software the entire process from spectral visualization to preparation of publication-ready figures is done solely using PINT and often within minutes, thereby, increasing productivity for users of all experience levels. Unique to the software are also the outstanding tools for evaluating the quality of the fitting results and extensive, but easy-to-use, customization of the fitting protocol and graphical output. In this communication, we describe the features of the new version of PINT and benchmark its performance.     

INFOS: spectrum fitting software for NMR analysis

Software for fitting of NMR spectra in MATLAB is presented. Spectra are fitted in the frequency domain, using Fourier transformed lineshapes, which are derived using the experimental acquisition and processing parameters. This yields more accurate fits compared to common fitting methods that use Lorentzian or Gaussian functions. Furthermore, a very time-efficient algorithm for calculating and fitting spectra has been developed. The software also performs initial peak picking, followed by subsequent fitting and refinement of the peak list, by iteratively adding and removing peaks to improve the overall fit. Estimation of error on fitting parameters is performed using a Monte-Carlo approach. Many fitting options allow the software to be flexible enough for a wide array of applications, while still being straightforward to set up with minimal user input.     

Simulation of NOESY spectra of DNA segments: A new scaling procedure for iterative comparison of calculated and experimental NOE intensities

Summary A new algorithm for simulation of two-dimensional NOESY spectra of DNA segments has been developed. For any given structure, NOE intensities are calculated using the relaxation matrix approach and a new realistic procedure is suggested for 1:1 comparison of calculated and experimental intensities. The procedure involves a novel method for scaling of calculated NOE intensities to represent volumes of digitised cross peaks in NOESY spectra. A data base of fine structures of all the relevant cross peaks with Lorentzian line shapes and in-phase components, is generated in a digitised manner by two-dimensional Fourier transformation of simulated time domain data, assuming a total intensity of 1.0 for each of the cross peaks. With this procedure, it is shown that the integrated volumes of these digitised cross peaks above any given threshold scale exactly as the total intensity of the respective peaks. This procedure eliminates the repetitive generation of digitised cross peaks by two-dimensional Fourier transformation during the iterative process of structure alteration and NOE intensity calculation and thus enhances the speed of DNA structure optimization. Illustrative fits of experimental and calculated spectra obtained using the new procedure are shown.[/p]     

Biochemical analysis of human PIF1 helicase and functions of its N-terminal domain

The evolutionary conserved PIF1 DNA helicase family appears to have largely nonoverlapping cellular functions. To better understand the functions of human PIF1, we investigated biochemical properties of this protein. Analysis of single-stranded (ss) DNA-dependent ATPase activity revealed nonstructural ssDNA to greatly stimulate ATPase activity due to a high affinity for PIF1, even though PIF1 preferentially unwinds forked substrates. This suggests that PIF1 needs a ssDNA region for loading and a forked structure for translocation entrance into a double strand region. Deletion analysis demonstrated novel functions of a unique N-terminal portion, named the PIF1 N-terminal (PINT) domain. When the PINT domain was truncated, apparent affinity for ssDNA and unwinding activity were much reduced, even though the maximum velocity of ATPase activity and K_m value for ATP were not affected. We suggest that the PINT domain contributes to enhancing the interaction with ssDNA through intrinsic binding activity. In addition, we found DNA strand-annealing activity, also residing in the PINT domain. Notably, the unwinding and annealing activities were inhibited by replication protein A. These results suggest that the functions of PIF1 might be restricted with particular situations and DNA structures.     

Line-shape analysis of NMR difference spectra of an anti-spin-label antibody

Specifically deuteriated Fab fragments of the anti-spin-label antibody AN02 were prepared. NMR difference spectra were obtained, in which the spectrum of Fab with some fraction of the binding sites occupied with spin-label hapten was subtracted from the spectrum of Fab with no spin-label. The peak heights were analyzed as a function of the fractional occupation of the binding site, using a computer program that calculates a best fit to the observed spectra. This method treats all of the peaks in the spectra simultaneously. Analyzing all peaks at once allows for the interdependencies in the spectra arising from overlap of positive and negative signals from different peaks. The fitting program calculates line widths for the peaks arising from protons in the binding site region. Almost all of the line widths calculated for the spectrum of the Fab complex with diamagnetic hapten dinitrophenyldiglycine were found to be narrower than the line widths of the corresponding resonances in the spectrum of Fab with an empty binding site. The distances of the binding site region protons from the unpaired electron of the hapten were also obtained from this calculation. Two tyrosine protons were found to be close (less than A) to this electron. These line-width and distance results are discussed with respect to the structure and dynamics of the antibody binding site.     

Precision enhancement of MALDI-TOF MS using high resolution peak detection and label-free alignment

We have developed an automated procedure for aligning peaks in multiple TOF spectra that eliminates common timing errors and small variations in spectrometer output. Our method incorporates high-resolution peak detection, re-binning, and robust linear data fitting in the time domain. This procedure aligns label-free (uncalibrated) peaks to minimize the variation in each peak's location from one spectrum to the next, while maintaining a high number of degrees of freedom. We apply our method to replicate pooled-serum spectra from multiple laboratories and increase peak precision (t/sigma(t)) to values limited only by small random errors (with sigma(t) less than one time count in 89 out of 91 instances, 13 peaks in seven datasets). The resulting high precision allowed for an order of magnitude improvement in peak m/z reproducibility. We show that the CV for m/z is 0.01% (100 ppm) for 12 out of the 13 peaks that were observed in all datasets between 2995 and 9297 Da.     

Solution structure of the coiled-coil trimerization domain from lung surfactant protein D

Surfactant protein D (SP-D) is one of four known protein components of the pulmonary surfactant lining the lung alveoli. It is involved in immune and allergic responses. SP-D occurs as a tetramer of trimers. Trimerization is thought to be initiated by a coiled coil domain. We have determined the solution structure of a 64-residue peptide encompassing the coiled coil domain of human SP-D. As predicted, the domain forms a triple-helical parallel coiled coil. As with all symmetric oligomers, the structure calculation was complicated by the symmetry degeneracy in the NMR spectra. We used the symmetry-ADR (ambiguous distance restraint) structure calculation method to solve the structure. The results demonstrate that the leucine zipper region of SP-D is an autonomously folded domain. The structure is very similar to the independently determined X-ray crystal structure, differing mainly at a single residue, Tyr248. This residue is completely symmetric in the solution structure, and markedly asymmetric in the crystalline phase. This difference may be functionally important, as it affects the orientation of the antigenic surface presented by SP-D.     

Fluorescence-quenching-resolved spectroscopy of proteins

A new procedure is described for using fluorescence-quenching data of tryptophan residues in proteins to resolve their fluorescence emission spectra. In this concept the Stern-Volmer quenching plot is determined at each particular emission wavelength and iterative non-linear least-squares fitting procedure allowed to resolve the steady-state emission spectra into components. The resolved components, attributed to each of tryptophan residue, can be characterized by different accessibility to the quencher. The ability to resolve fluorescence emission spectra can be improved by using different kinds of efficient quenchers, which can selectively quench the emission of exposed or both exposed and buried fluorophores. The method was used to decompose emission fluorescence spectra in two-tryptophan-containing proteins; horse liver dehydrogenase, sperm whale apomyoglobin and metalloprotease from Staphylococcus aureus. The resolved spectra of alcohol dehydrogenase and metalloprotease are in excellent agreement with those previously obtained by single-photon counting or phase methods. The method presented here is technically simple and does not require expensive instrumentation.     

Characterization of the photoactive GAF domain of the CikA homolog (SyCikA, Slr1969) of the cyanobacterium Synechocystis sp. PCC 6803

Phytochromes and bacteriophytochromes in plants and some species of bacteria, respectively, are photoreceptors that bind linear tetrapyrroles and can respond to red and far-red light signals in a reversible manner. A related but distinct photoreceptor candidate, CikA (denoted ScCikA), has been reported to reset the circadian clock in the cyanobacterium Synechococcus elongatus PCC 7942 after a dark pulse. However, recent studies have indicated that ScCikA does not function as a photoreceptor but as a redox sensor. Moreover, the Cys residue that covalently ligates the chromophore in phytochromes is not conserved in the ScCikA protein. On the other hand, the CikA homolog in Synechocystis sp. PCC 6803 (Slr1969, denoted SyCikA) retains this conserved Cys residue. In our present study, we have isolated the putative chromophore-binding GAF domain of SyCikA from Synechocystis and phycocyanobilin-producing Escherichia coli. Absorption spectra of both preparations showed two peaks in the UV and violet regions. Irradiation of these proteins with violet light yielded a broad peak in a yellow region at the expense of the peaks in the UV and violet regions. Interestingly, successive irradiation with yellow light did not revert these absorption spectra but a partial dark reversion to the original form was detected. These results suggest that SyCikA may function as a violet light sensor in Synechocystis.     

Global analysis of three-state protein unfolding data

A new method for analyzing three-state protein unfolding equilibria is described that overcomes the difficulties created by direct effects of denaturants on circular dichroism (CD) and fluorescence spectra of the intermediate state. The procedure begins with a singular value analysis of the data matrix to determine the number of contributing species and perturbations. This result is used to choose a fitting model and remove all spectra from the fitting equation. Because the fitting model is a product of a matrix function which is nonlinear in the thermodynamic parameters and a matrix that is linear in the parameters that specify component spectra, the problem is solved with a variable projection algorithm. Advantages of this procedure are perturbation spectra do not have to be estimated before fitting, arbitrary assumptions about magnitudes of parameters that describe the intermediate state are not required, and multiple experiments involving different spectroscopic techniques can be simultaneously analyzed. Two tests of this method were performed: First, simulated three-state data were analyzed, and the original and recovered thermodynamic parameters agreed within one standard error, whereas recovered and original component spectra agreed within 0.5%. Second, guanidine-induced unfolding titrations of the human retinoid-X-receptor ligand-binding domain were analyzed according to a three-state model. The standard unfolding free energy changes in the absence of guanidine and the guanidine concentrations at zero free-energy change for both transitions were determined from a joint analysis of fluorescence and CD spectra. Realistic spectra of the three protein states were also obtained.     

A Novel Preprocessing Method Using Hilbert Huang Transform for MALDI-TOF and SELDI-TOF Mass Spectrometry Data

Motivation

Mass spectrometry is a high throughput, fast, and accurate method of protein analysis. Using the peaks detected in spectra, we can compare a normal group with a disease group. However, the spectrum is complicated by scale shifting and is also full of noise. Such shifting makes the spectra non-stationary and need to align before comparison. Consequently, the preprocessing of the mass data plays an important role during the analysis process. Noises in mass spectrometry data come in lots of different aspects and frequencies. A powerful data preprocessing method is needed for removing large amount of noises in mass spectrometry data.

Results

Hilbert-Huang Transformation is a non-stationary transformation used in signal processing. We provide a novel algorithm for preprocessing that can deal with MALDI and SELDI spectra. We use the Hilbert-Huang Transformation to decompose the spectrum and filter-out the very high frequencies and very low frequencies signal. We think the noise in mass spectrometry comes from many sources and some of the noises can be removed by analysis of signal frequence domain. Since the protein in the spectrum is expected to be a unique peak, its frequence domain should be in the middle part of frequence domain and will not be removed. The results show that HHT, when used for preprocessing, is generally better than other preprocessing methods. The approach not only is able to detect peaks successfully, but HHT has the advantage of denoising spectra efficiently, especially when the data is complex. The drawback of HHT is that this approach takes much longer for the processing than the wavlet and traditional methods. However, the processing time is still manageable and is worth the wait to obtain high quality data.     

A series of lysyldipeptide derivatives for racemization studies in peptide synthesis.

The series of diastereomeric peptide derivatives N-benzoyl-D,L-X-N epsilon-benzyloxycarbonyl-L-lysine methyl ester where X = alanyl, valyl, leucyl, phenylalanyl and isoleucyl are submitted as model systems for studying racemization in peptide synthesis. The diastereomers can be analyzed by quantitation of the separated ester methyl proton peaks of their nuclear magnetic resonance spectra. Data on the tendency to racemize of the different residues are presented. In polar solvents, valyl and isoleucyl residues racemize more readily than the other residues.     

Secondary structural analysis of proteins based on 13C chemical shift assignments in unresolved solid-state NMR spectra enhanced by fragmented structure database

Magic-angle-spinning solid-state ¹³C NMR spectroscopy is useful for structural analysis of non-crystalline proteins. However, the signal assignments and structural analysis are often hampered by the signal overlaps primarily due to minor structural heterogeneities, especially for uniformly-¹³C,¹⁵N labeled samples. To overcome this problem, we present a method for assigning ¹³C chemical shifts and secondary structures from unresolved two-dimensional ¹³C–¹³C MAS NMR spectra by spectral fitting, named reconstruction of spectra using protein local structures (RESPLS). The spectral fitting was conducted using databases of protein fragmented structures related to ¹³C^α, ¹³C^β, and ¹³C′ chemical shifts and cross-peak intensities. The experimental ¹³C–¹³C inter- and intra-residue correlation spectra of uniformly isotope-labeled ubiquitin in the lyophilized state had a few broad peaks. The fitting analysis for these spectra provided sequence-specific C^α, C^β, and C′ chemical shifts with an accuracy of about 1.5 ppm, which enabled the assignment of the secondary structures with an accuracy of 79 %. The structural heterogeneity of the lyophilized ubiquitin is revealed from the results. Test of RESPLS analysis for simulated spectra of five different types of proteins indicated that the method allowed the secondary structure determination with accuracy of about 80 % for the 50–200 residue proteins. These results demonstrate that the RESPLS approach expands the applicability of the NMR to non-crystalline proteins exhibiting unresolved ¹³C NMR spectra, such as lyophilized proteins, amyloids, membrane proteins and proteins in living cells.     

Two-dimensional 1H nuclear magnetic resonance studies on the gene V-encoded single-stranded DNA-binding protein of the filamentous bacteriophage IKe. II. Characterization of the DNA-binding wing with the aid of spin-labelled oligonucleotides

The DNA-binding domain of the single-stranded DNA-binding protein IKe GVP was studied by means of 1H nuclear magnetic resonance, through use of oligonucleotides of two and three adenyl residues in length, that were spin-labelled at their 3' and/or 5' termini. These spin-labelled ligands were found to cause line broadening of specific protein resonances when bound to the protein, although they were present in small quantities, i.e. of the order of 0.04 molar equivalent and less. The line broadening of protein resonances was made manifest by means of difference one and two-dimensional spectroscopy. Difference one-dimensional experiments revealed line broadening of the same protein resonances upon binding of either 3' or 5' spin-labelled oligonucleotides. Evidence in favour of the existence of a fixed 5' to 3' orientation in the binding of oligonucleotides to the protein surface was therefore not obtained from the spin-labelled oligonucleotide binding studies. Residue-specific assignments of broadened resonances could not, or could only sparsely, be derived from the difference one-dimensional spectra, because of the tremendous overlap in the aliphatic region of the spectrum. In contrast, such assignments were easily obtained from the difference two-dimensional spectra, which were recorded by means of both total correlated spectroscopy and nuclear Overhauser effect spectroscopy. Difference signals were detected for 15 spin systems; ten out of these were assigned to the residues I29, Y27, S20, G18, R16, T28, K22, Q21, V19 and S17 in the amino acid sequence of IKe GVP; the other five spin systems could be assigned to a phenylalanyl residue, an arginyl or lysyl residue, an aspartic acid or asparagyl residue, a glycyl residue and a glutamic acid or glutamyl residue. From the evaluation of the relative difference signals, it was concluded that the direct surroundings of the spin-label group of the labelled oligonucleotide in the bound state is composed of the first five residues in the former group of residues and the five residues in the latter group.(ABSTRACT TRUNCATED AT 400 WORDS)     

Quantification of two-dimensional NOE spectra via a combined linear and nonlinear least-squares fit

Summary Determining the volumes of peaks in 2D NMR spectra can be prohibitively difficult in cases of overlapping, broad lines. Deconvolution and parameter estimation can be attempted on either the time-domain or the frequency-domain data. We present a method of estimating spectral parameters from frequency-domain data, using a combination of Lorentzian and Gaussian lineshapes for reference lines. This approach combines a previously published method of projecting the data on a linear space spanned by reference lines with a nonlinear least-squares fitting algorithm. Comparison of this method with other published methods of frequency-domain deconvolution shows that it is both more precise and more accurate when estimating 2D volumes.     

Side chain and backbone dynamics of phospholamban in phospholipid bilayers utilizing 2H and 15N solid-state NMR spectroscopy

2H and 15N solid-state NMR spectroscopic techniques were used to investigate both the side chain and backbone dynamics of wild-type phospholamban (WT-PLB) and its phosphorylated form (P-PLB) incorporated into 1-palmitoyl-2-oleoyl-sn-glycerophosphocholine (POPC) phospholipid bilayers. 2H NMR spectra of site-specific CD3-labeled WT-PLB (at Leu51, Ala24, and Ala15) in POPC bilayers were similar under frozen conditions (-25 degrees C). However, significant differences in the line shapes of the 2H NMR spectra were observed in the liquid crystalline phase at and above 0 degrees C. The 2H NMR spectra indicate that Leu51, located toward the lower end of the transmembrane (TM) helix, shows restricted side chain motion, implying that it is embedded inside the POPC lipid bilayer. Additionally, the line shape of the 2H NMR spectrum of CD3-Ala24 reveals more side chain dynamics, indicating that this residue (located in the upper end of the TM helix) has additional backbone and internal side chain motions. 2H NMR spectra of both WT-PLB and P-PLB with CD3-Ala15 exhibit strong isotropic spectral line shapes. The dynamic isotropic nature of the 2H peak can be attributed to side chain and backbone motions to residues located in an aqueous environment outside the membrane. Also, the spectra of 15N-labeled amide WT-PLB at Leu51 and Leu42 residues showed only a single powder pattern component indicating that these two 15N-labeled residues located in the TM helix are motionally restricted at 25 degrees C. Conversely, 15N-labeled amide WT-PLB at Ala11 located in the cytoplasmic domain showed both powder and isotropic components at 25 degrees C. Upon phosphorylation, the mobile component contribution increases at Ala11. The 2H and 15N NMR data indicate significant backbone motion for the cytoplasmic domain of WT-PLB when compared to the transmembrane section.     

Processing of ND NMR spectra sampled in polar coordinates: a simple Fourier transform instead of a reconstruction

In order to reduce the acquisition time of multidimensional NMR spectra of biological macromolecules, projected spectra (or in other words, spectra sampled in polar coordinates) can be used. Their standard processing involves a regular FFT of the projections followed by a reconstruction, i.e. a non-linear process. In this communication, we show that a 2D discrete Fourier transform can be implemented in polar coordinates to obtain directly a frequency domain spectrum. Aliasing due to local violations of the Nyquist sampling theorem gives rise to base line ridges but the peak line-shapes are not distorted as in most reconstruction methods. The sampling scheme is not linear and the data points in the time domain should thus be weighted accordingly in the polar FT; however, artifacts can be reduced by additional data weighting of the undersampled regions. This processing does not require any parameter tuning and is straightforward to use. The algorithm written for polar sampling can be adapted to any sampling scheme and will permit to investigate better compromises in terms of experimental time and lack of artifacts.     

Solution NMR of a 463-residue phosphohexomutase: domain 4 mobility, substates, and phosphoryl transfer defect

Phosphomannomutase/phosphoglucomutase contributes to the infectivity of Pseudomonas aeruginosa, retains and reorients its intermediate by 180°, and rotates domain 4 to close the deep catalytic cleft. Nuclear magnetic resonance (NMR) spectra of the backbone of wild-type and S108C-inactivated enzymes were assigned to at least 90%. (13)C secondary chemical shifts report excellent agreement of solution and crystallographic structure over the 14 α-helices, C-capping motifs, and 20 of the 22 β-strands. Major and minor NMR peaks implicate substates affecting 28% of assigned residues. These can be attributed to the phosphorylation state and possibly to conformational interconversions. The S108C substitution of the phosphoryl donor and acceptor slowed transformation of the glucose 1-phosphate substrate by impairing k(cat). Addition of the glucose 1,6-bisphosphate intermediate accelerated this reaction by 2-3 orders of magnitude, somewhat bypassing the defect and apparently relieving substrate inhibition. The S108C mutation perturbs the NMR spectra and electron density map around the catalytic cleft while preserving the secondary structure in solution. Diminished peak heights and faster (15)N relaxation suggest line broadening and millisecond fluctuations within four loops that can contact phosphosugars. (15)N NMR relaxation and peak heights suggest that domain 4 reorients slightly faster in solution than domains 1-3, and with a different principal axis of diffusion. This adds to the crystallographic evidence of domain 4 rotations in the enzyme, which were previously suggested to couple to reorientation of the intermediate, substrate binding, and product release.     

Site-specific NMR monitoring of cis-trans isomerization in the folding of the proline-rich collagen triple helix

Understanding the folding of the proline-rich collagen triple helix requires consideration of the effects of proline cis-trans isomerization and may shed light on the misfolding of collagen in connective tissue diseases. Folding was monitored in real time by heteronuclear 2D NMR spectroscopy for the (15)N labeled positions in the triple-helical peptide T1-892 [GPAGPAGPVGPAGARGPAGPOGPOGPOGPOGV]. In the equilibrium unfolded monomer form, each labeled residue showed multiple peaks with interconversion rates consistent with cis-trans isomerization of Gly-Pro and Pro-Hyp bonds. Real-time NMR studies on the folding of T1-892 showed slow decay of monomer peaks and a concomitant increase in trimer peaks. Gly25 in the C-terminal rich (Gly-Pro-Hyp)(4) domain folds first, consistent with its being a nucleation domain. Analysis of the kinetics indicates that the folding of Gly25 is biphasic and the slower step represents cis-trans isomerization of imino acids. This illustrates that nucleation is limited by cis-trans isomerization. Monitoring Gly6, Gly10, Ala12, and Gly13 monomer and trimer peaks captures the C- to N-terminal propagation of the triple helix, which is also limited by Gly-Pro cis-trans isomerization events. The zipper-like nature of the propagation process is confirmed by the slower rate of folding of Ala6 compared to Gly13, reflecting the larger number of isomerization events encountered by the more N-terminal Ala6. The cis-trans isomerization events at multiple proline residues is a complex statistical process which can be visualized by these NMR studies.     

Isolation and Identification of 2,3-Dimethyl-4-hydroxy-2- nonenoic Acid Lactone from Shubi

The heating of 7S globulin caused changes in the intensities, but hardly affected the positions of the peaks and troughs of the second derivative absorption spectra at wavelengths below 270 nm. On the other hand, above 271 nm, changes were reflected both in the intensities and in the positions of peaks and troughs. The difference-second derivative absorption spectra indicated that 60 and 70 percent, respectively, of phenylalanine and tyrosine residues buried in the native 7S globulin remained as the buried form even after heating.

A spectrofluorimetry and fluorescence-quenching study suggested that one residue of tryptophan in the 7S globulin tended to be transferred to the more hydrophobic interior on heating.