Search

The ability to reconstruct non-uniformly sampled (NUS) NMR spectra has mostly been accepted. Still a concern is lingering regarding artifacts from sampling non-uniformly. As experienced, some sampling schedules yield better results than others. Finding a useful schedule is relatively trivial for a low dynamic range spectrum and a conservative sparsity, but not so when the dynamic range is large and/or when extreme sparsity is used. High dynamic range is typically found in NOESY and spectra of metabolites, where quantification of peak heights is desired at high fidelity. Extreme sparsity is desired when high throughput is a goal. In all cases, selecting a poor sampling schedule can create unnecessary artifacts. Effectively, it is important to select a sampling schedule that provides a signal-to-artifact apex ratio (SAAR) value in par or better than the signal-to-noise ratio (SNR) value. Notably, by signal-to-artifact apex ratio we consider reconstruction fidelity as the apex intensity likeness, i.e., as the true signal to the tallest artifact. We show that the quality of reconstruction depends on the particular sampling schedule. We evaluate the reconstruction quality in the frequency domain following a matched Lorentz-to-Gauss transform plus common apodization and Fourier Transform. As the Lorentz-to-Gauss transform improves resolution and reduces ridges we include this when defining the Signal-to-Artifact Apex Ratio (SAAR) metric. This metric measures the ratio of simulated reconstructed peak height to the tallest artifact of reconstruction in a spectrum without noise. Once a NUS schedule is found with an optimal SAAR it will be satisfactory for all spectra recorded with the same parameter set. Tables with good seed values are provided in the supplement., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 Elsevier Inc. All rights reserved.)

Although the concepts of nonuniform sampling (NUS​​​​​​​) and non-Fourier spectral reconstruction in multidimensional NMR began to emerge 4 decades ago , it is only relatively recently that NUS has become more commonplace. Advantages of NUS include the ability to tailor experiments to reduce data collection time and to improve spectral quality, whether through detection of closely spaced peaks (i.e., "resolution") or peaks of weak intensity (i.e., "sensitivity"). Wider adoption of these methods is the result of improvements in computational performance, a growing abundance and flexibility of software, support from NMR spectrometer vendors, and the increased data sampling demands imposed by higher magnetic fields. However, the identification of best practices still remains a significant and unmet challenge. Unlike the discrete Fourier transform, non-Fourier methods used to reconstruct spectra from NUS data are nonlinear, depend on the complexity and nature of the signals, and lack quantitative or formal theory describing their performance. Seemingly subtle algorithmic differences may lead to significant variabilities in spectral qualities and artifacts. A community-based critical assessment of NUS challenge problems has been initiated, called the "Nonuniform Sampling Contest" (NUScon), with the objective of determining best practices for processing and analyzing NUS experiments. We address this objective by constructing challenges from NMR experiments that we inject with synthetic signals, and we process these challenges using workflows submitted by the community. In the initial rounds of NUScon our aim is to establish objective criteria for evaluating the quality of spectral reconstructions. We present here a software package for performing the quantitative analyses, and we present the results from the first two rounds of NUScon. We discuss the challenges that remain and present a roadmap for continued community-driven development with the ultimate aim of providing best practices in this rapidly evolving field. The NUScon software package and all data from evaluating the challenge problems are hosted on the NMRbox platform., Competing Interests: Some authors are members of the editorial board of Magnetic Resonance. The peer-review process was guided by an independent editor, and the authors have also no other competing interests to declare., (Copyright: © 2021 Yulia Pustovalova et al.)

Nonuniform sampling was first proposed more than 40 years ago as an alternate method for sampling two-dimensional NMR data was initially pursued by only a small number of scientists. However, it has been gradually adopted after it was shown that major gains in measuring time and spectrum resolution can be obtained. Furthermore, migration of NMR software to the Unix environment facilitated development of new processing tools, and there is now a selection of programs available that yield high-quality reconstructions of NUS data. Moreover, it became obvious that recording high-resolution 3D and 4D protein NMR spectra at the resolution provided by modern high-field instruments was not possible with uniform sampling. It has become apparent that sparse, low dynamic-range NMR spectra, in particular, the triple resonance experiments are all best recorded with NUS. Optimal sampling schedules can yield benefits with respect to detecting weak peaks in high dynamic-range spectra and therefore a careful use of 2D HSQC-like spectra of mixed concentrations of small molecules is feasible. It is not yet clear whether crowded high dynamic-range spectra, such as NOESYs with many cross-peaks, benefit from NUS. On the other hand, NUS appears to be the best option for recording such high dynamic-range spectra if crowding can be reduced by shorten NOESY mixing times and/or by reducing overlap through isotopic labeling in high-resolution 3D and 4D experiments. Here we discuss principals and applications of the method, trying to provide a wide overview, but are biased by our own approaches. Thus, we will obviously miss important developments elsewhere and do not claim to be comprehensive., (© 2019 Elsevier Inc. All rights reserved.)

Non-Uniform Sampling has the potential to exploit the optimal resolution of high-field NMR instruments. This is not possible in 3D and 4D NMR experiments when using traditional uniform sampling due to the long overall measurement time. Nominally, uniformly sampled time domain data acquired to a maximum evolution time t max can be extended to high resolution via a virtual maximum evolution time t* max while extrapolating with linear prediction or iterative soft thresholding (IST). At the high resolution obtainable with extrapolation of US data, however, the accuracy of peak positions is compromised as observed when comparing inter- and intra-residue peaks in a 3D HNCA experiment. However, the accuracy of peak positions is largely improved by spreading the same number of acquired time domain data points non-uniformly over a larger evolution time to an optimal t max followed by extrapolation to a total t* max and processing the data with an appropriate reconstruction method, such as hmsIST. To explore the optimum value of experimentally measured t max to be reached non-uniformly with a given number of sampling points we have created test situations of time-equivalent experiments and evaluate sensitivity and accuracy of peak positions. Here we use signal-to-maximum-noise ratio as the decisive measure of sensitivity. We find that both sensitivity and resolution are optimal when PoissonGap sampling to a t max of about ½*T 2 * . Digital resolution is further enhanced by extrapolating the range of acquired time domain data to 2*T 2 * but without measuring experimental points beyond ½*T 2 * .

Isotopomer analysis using either 13 C NMR or LC/GC-MS has been an invaluable tool for studying metabolic activities in a variety of systems. Traditional challenges are, however, that 13 C-detected NMR is insensitive despite its high resolution, and that MS-based techniques cannot easily differentiate positional isotopomers. In addition, current 13 C NMR or LC/GC-MS has limitations in detecting metabolites in living cells. Here, we describe a non-uniform sampling-based 2D heteronuclear single quantum coherence (NUS HSQC) approach to measure metabolic isotopomers in both cell lysates and living cells. The method provides a high resolution that can resolve multiplet structures in the 13 C dimension while retaining the sensitivity of the 1 H-indirect detection. The approach was tested in L1210 mouse leukemia cells labeled with 13 C acetate by measuring NUS HSQC with 25% sampling density. The results gave a variety of metabolic information such as (1) higher usage of acetate in acetylation pathway than aspartate synthesis, (2) TCA cycle efficiency changes upon the inhibition of mitochondrial oxidative phosphorylation by pharmacological agents, and (3) position-dependent isotopomer patterns in fatty acids in living cells. In addition, we were able to detect fatty acids along with other hydrophilic molecules in one sample of live cells without extraction. Overall, the high sensitivity and resolution along with the application to live cells should make the NUS HSQC approach attractive in studying carbon flux information in metabolic studies.

Structural characterization of membrane proteins and other large proteins with NMR relies increasingly on perdeuteration combined with incorporation of specifically protonated amino acid moieties, such as methyl groups of isoleucines, valines, or leucines. The resulting proton dilution reduces dipolar broadening producing sharper resonance lines, ameliorates spectral crowding, and enables measuring of crucial distances between and to methyl groups. While incorporation of specific methyl labeling is now well established for bacterial expression using suitable precursors, corresponding methods are still lacking for cell-free expression, which is often the only choice for producing labeled eukaryotic membrane proteins in mg quantities. Here we show that we can express methyl-labeled human integral membrane proteins cost-effectively by cell-free expression based of crude hydrolyzed ILV-labeled OmpX inclusion bodies. These are obtained in Escherichia coli with very high quantity and represent an optimal intermediate to channel ILV precursors into the eukaryotic proteins.

We report acquisition of diagonal-compensated protein structural restraints from four-dimensional solid-state NMR spectra on extensively deuterated and (1)H back-exchanged proteins. To achieve this, we use homonuclear (1)H-(1)H correlations with diagonal suppression and nonuniform sampling (NUS). Suppression of the diagonal allows the accurate identification of cross-peaks which are otherwise obscured by the strong autocorrelation or whose intensity is biased due to partial overlap with the diagonal. The approach results in unambiguous spectral interpretation and relatively few but reliable restraints for structure calculation. In addition, the diagonal suppression produces a spectrum with low dynamic range for which ultrasparse NUS data sets can be readily reconstructed, allowing straightforward application of NUS with only 2% sampling density with the advantage of more heavily sampling time-domain regions of high signal intensity. The method is demonstrated here for two proteins, α-spectrin SH3 microcrystals and hydrophobin functional amyloids. For the case of SH3, suppression of the diagonal results in facilitated identification of unambiguous restraints and improvement of the quality of the calculated structural ensemble compared to nondiagonal-suppressed 4D spectra. For the only partly assigned hydrophobin rodlets, the structure is yet unknown. Applied to this protein of biological significance with large inhomogeneous broadening, the method allows identification of unambiguous crosspeaks that are otherwise obscured by the diagonal.

Multi-dimensional NMR spectra have traditionally been processed with the fast Fourier transformation (FFT). The availability of high field instruments, the complexity of spectra of large proteins, the narrow signal dispersion of some unstructured proteins, and the time needed to record the necessary increments in the indirect dimensions to exploit the resolution of the highfield instruments make this traditional approach unsatisfactory. New procedures need to be developed beyond uniform sampling of the indirect dimensions and reconstruction methods other than the straight FFT are necessary. Here we discuss approaches of non-uniform sampling (NUS) and suitable reconstruction methods. We expect that such methods will become standard for multi-dimensional NMR data acquisition with complex biological macromolecules and will dramatically enhance the power of modern biological NMR., (Copyright © 2013 Elsevier Inc. All rights reserved.)

It is well established that non-uniform sampling (NUS) allows acquisition of multi-dimensional NMR spectra at a resolution that cannot be obtained with traditional uniform acquisition through the indirect dimensions. However, the impact of NUS on the signal-to-noise ratio (SNR) and sensitivity are less well documented. SNR and sensitivity are essential aspects of NMR experiments as they define the quality and extent of data that can be obtained. This is particularly important for spectroscopy with low concentration samples of biological macromolecules. There are different ways of defining the SNR depending on how to measure the noise, and the distinction between SNR and sensitivity is often not clear. While there are defined procedures for measuring sensitivity with high concentration NMR standards, such as sucrose, there is no clear or generally accepted definition of sensitivity when comparing different acquisition and processing methods for spectra of biological macromolecules with many weak signals close to the level of noise. Here we propose tools for estimating the SNR and sensitivity of NUS spectra with respect to sampling schedule and reconstruction method. We compare uniformly acquired spectra with NUS spectra obtained in the same total measuring time. The time saving obtained when only 1/k of the Nyquist grid points are sampled is used to measure k-fold more scans per increment. We show that judiciously chosen NUS schedules together with suitable reconstruction methods can yield a significant increase of the SNR within the same total measurement time. Furthermore, we propose to define the sensitivity as the probability to detect weak peaks and show that time-equivalent NUS can considerably increase this detection sensitivity. The sensitivity gain increases with the number of NUS indirect dimensions. Thus, well-chosen NUS schedules and reconstruction methods can significantly increase the information content of multidimensional NMR spectra of challenging biological macromolecules.

Recognition of the proper start codon on mRNAs is essential for protein synthesis, which requires scanning and involves eukaryotic initiation factors (eIFs) eIF1, eIF1A, eIF2, and eIF5. The carboxyl terminal domain (CTD) of eIF5 stimulates 43S preinitiation complex (PIC) assembly; however, its precise role in scanning and start codon selection has remained unknown. Using nuclear magnetic resonance (NMR) spectroscopy, we identified the binding sites of eIF1 and eIF2β on eIF5-CTD and found that they partially overlapped. Mutating select eIF5 residues in the common interface specifically disrupts interaction with both factors. Genetic and biochemical evidence indicates that these eIF5-CTD mutations impair start codon recognition and impede eIF1 release from the PIC by abrogating eIF5-CTD binding to eIF2β. This study provides mechanistic insight into the role of eIF5-CTD's dynamic interplay with eIF1 and eIF2β in switching PICs from an open to a closed state at start codons., (Copyright © 2012 The Authors. Published by Elsevier Inc. All rights reserved.)

The fast Fourier transformation has been the gold standard for transforming data from time to frequency domain in many spectroscopic methods, including NMR. While reliable, it has as a drawback that it requires a grid of uniformly sampled data points. This needs very long measuring times for sampling in multidimensional experiments in all indirect dimensions uniformly and even does not allow reaching optimal evolution times that would match the resolution power of modern high-field instruments. Thus, many alternative sampling and transformation schemes have been proposed. Their common challenges are the suppression of the artifacts due to the non-uniformity of the sampling schedules, the preservation of the relative signal amplitudes, and the computing time needed for spectra reconstruction. Here we present a fast implementation of the Iterative Soft Thresholding approach (istHMS) that can reconstruct high-resolution non-uniformly sampled NMR data up to four dimensions within a few hours and make routine reconstruction of high-resolution NUS 3D and 4D spectra convenient. We include a graphical user interface for generating sampling schedules with the Poisson-Gap method and an estimation of optimal evolution times based on molecular properties. The performance of the approach is demonstrated with the reconstruction of non-uniformly sampled medium and high-resolution 3D and 4D protein spectra acquired with sampling densities as low as 0.8%. The method presented here facilitates acquisition, reconstruction and use of multidimensional NMR spectra at otherwise unreachable spectral resolution in indirect dimensions.

Modern high-field NMR instruments provide unprecedented resolution. To make use of the resolving power in multidimensional NMR experiment standard linear sampling through the indirect dimensions to the maximum optimal evolution times (~1.2T (2)) is not practical because it would require extremely long measurement times. Thus, alternative sampling methods have been proposed during the past 20 years. Originally, random nonlinear sampling with an exponentially decreasing sampling density was suggested, and data were transformed with a maximum entropy algorithm (Barna et al., J Magn Reson 73:69-77, 1987). Numerous other procedures have been proposed in the meantime. It has become obvious that the quality of spectra depends crucially on the sampling schedules and the algorithms of data reconstruction. Here we use the forward maximum entropy (FM) reconstruction method to evaluate several alternate sampling schedules. At the current stage, multidimensional NMR spectra that do not have a serious dynamic range problem, such as triple resonance experiments used for sequential assignments, are readily recorded and faithfully reconstructed using non-uniform sampling. Thus, these experiments can all be recorded non-uniformly to utilize the power of modern instruments. On the other hand, for spectra with a large dynamic range, such as 3D and 4D NOESYs, choosing optimal sampling schedules and the best reconstruction method is crucial if one wants to recover very weak peaks. Thus, this chapter is focused on selecting the best sampling schedules and processing methods for high-dynamic range spectra.

Previous experiments suggest a connection between the N-alpha-acetylation of proteins and sensitivity of cells to apoptotic signals. Here, we describe a biochemical assay to detect the acetylation status of proteins and demonstrate that protein N-alpha-acetylation is regulated by the availability of acetyl-CoA. Because the antiapoptotic protein Bcl-xL is known to influence mitochondrial metabolism, we reasoned that Bcl-xL may provide a link between protein N-alpha-acetylation and apoptosis. Indeed, Bcl-xL overexpression leads to a reduction in levels of acetyl-CoA and N-alpha-acetylated proteins in the cell. This effect is independent of Bax and Bak, the known binding partners of Bcl-xL. Increasing cellular levels of acetyl-CoA by addition of acetate or citrate restores protein N-alpha-acetylation in Bcl-xL-expressing cells and confers sensitivity to apoptotic stimuli. We propose that acetyl-CoA serves as a signaling molecule that couples apoptotic sensitivity to metabolism by regulating protein N-alpha-acetylation., (Copyright © 2011 Elsevier Inc. All rights reserved.)

The primary limitation of solution state NMR with larger, highly dynamic, or paramagnetic systems originates from signal losses due to fast transverse relaxation. This is related to the high gyromagnetic ratio gamma of protons, which are usually detected. Thus, it is attractive to consider detection of nuclei with lower gamma, such as (13)C, for extending the size limits of NMR. Here, we present an approach for complete assignment of C(alpha) and N resonances in fast relaxing proteins using a C(alpha) detected 3D CANCA experiment for perdeuterated proteins. The CANCA experiment correlates alpha carbons with the sequentially adjacent and succeeding nitrogen and alpha carbons. This enables elongation of the chain of assigned residues simply by navigating along both nitrogen and carbon dimensions using a "stairway" assignment procedure. The simultaneous use of both C(alpha) and N sequential connectivities makes the experiment more robust than conventional 3D experiments, which rely solely on a single (13)C indirect dimension for sequential information. The 3D CANCA experiment, which is very useful for mainchain assignments of higher molecular weight proteins at high magnetic field, also provides an attractive alterative for smaller proteins. Two versions of the experiment are described for samples that are (13)C labeled either uniformly or at alternate positions for removing one-bond (13)C-(13)C couplings. To achieve both high resolution and sensitivity, extensive nonuniform sampling was employed. Adding longitudinal relaxation enhancement agents can allow for shorter recycling delays, decreased measuring time, or enhanced sensitivity.

The Fourier transform has been the gold standard for transforming data from the time domain to the frequency domain in many spectroscopic methods, including NMR spectroscopy. While reliable, it has the drawback that it requires a grid of uniformely sampled data points, which is not efficient for decaying signals, and it also suffers from artifacts when dealing with nondecaying signals. Over several decades, many alternative sampling and transformation schemes have been proposed. Their common problem is that relative signal amplitudes are not well-preserved. Here we demonstrate the superior performance of a sine-weighted Poisson-gap distribution sparse-sampling scheme combined with forward maximum entropy (FM) reconstruction. While the relative signal amplitudes are well-preserved, we also find that the signal-to-noise ratio is enhanced up to 4-fold per unit of data acquisition time relative to traditional linear sampling.

Non-uniform sampling (NUS) enables recording of multidimensional NMR data at resolutions matching the resolving power of modern instruments without using excessive measuring time. However, in order to obtain satisfying results, efficient reconstruction methods are needed. Here we describe an optimized version of the Forward Maximum entropy (FM) reconstruction method, which can reconstruct up to three indirect dimensions. For complex datasets, such as NOESY spectra, the performance of the procedure is enhanced by a distillation procedure that reduces artifacts stemming from intense peaks.

Programmed Cell Death 4 (PDCD4) is a protein known to bind eukaryotic initiation factor 4A (eIF4A), inhibit translation initiation, and act as a tumor suppressor. PDCD4 contains two C-terminal MA3 domains, which are thought to be responsible for its inhibitory function. Here, we analyze the structures and inhibitory functions of these two PDCD4 MA3 domains by x-ray crystallography, NMR, and surface plasmon resonance. We show that both MA3 domains are structurally and functionally very similar and bind specifically to the eIF4A N-terminal domain (eIF4A-NTD) using similar binding interfaces. We found that the PDCD4 MA3 domains compete with the eIF4G MA3 domain and RNA for eIF4A binding. Our data provide evidence that PDCD4 inhibits translation initiation by displacing eIF4G and RNA from eIF4A. The PDCD4 MA3 domains act synergistically to form a tighter and more stable complex with eIF4A, which explains the need for two tandem MA3 domains.

To obtain a comprehensive assessment of metabolite levels from extracts of leukocytes, we have recorded ultrahigh-resolution 1H-13C HSQC NMR spectra of cell extracts, which exhibit spectral signatures of numerous small molecules. However, conventional acquisition of such spectra is time-consuming and hampers measurements on multiple samples, which would be needed for statistical analysis of metabolite concentrations. Here we show that the measurement time can be dramatically reduced without loss of spectral quality when using nonlinear sampling (NLS) and a new high-fidelity forward maximum-entropy (FM) reconstruction algorithm. This FM reconstruction conserves all measured time-domain data points and guesses the missing data points by an iterative process. This consists of discrete Fourier transformation of the sparse time-domain data set, computation of the spectral entropy, determination of a multidimensional entropy gradient, and calculation of new values for the missing time-domain data points with a conjugate gradient approach. Since this procedure does not alter measured data points, it reproduces signal intensities with high fidelity and does not suffer from a dynamic range problem. As an example we measured a natural abundance 1H-13C HSQC spectrum of metabolites from granulocyte cell extracts. We show that a high-resolution 1H-13C HSQC spectrum with 4k complex increments recorded linearly within 3.7 days can be reconstructed from one-seventh of the increments with nearly identical spectral appearance, indistinguishable signal intensities, and comparable or even lower root-mean-square (rms) and peak noise patterns measured in signal-free areas. Thus, this approach allows recording of ultrahigh resolution 1H-13C HSQC spectra in a fraction of the time needed for recording linearly sampled spectra.

The sterol regulatory element binding protein (SREBP) family of transcription activators are critical regulators of cholesterol and fatty acid homeostasis. We previously demonstrated that human SREBPs bind the CREB-binding protein (CBP)/p300 acetyltransferase KIX domain and recruit activator-recruited co-factor (ARC)/Mediator co-activator complexes through unknown mechanisms. Here we show that SREBPs use the evolutionarily conserved ARC105 (also called MED15) subunit to activate target genes. Structural analysis of the SREBP-binding domain in ARC105 by NMR revealed a three-helix bundle with marked similarity to the CBP/p300 KIX domain. In contrast to SREBPs, the CREB and c-Myb activators do not bind the ARC105 KIX domain, although they interact with the CBP KIX domain, revealing a surprising specificity among structurally related activator-binding domains. The Caenorhabditis elegans SREBP homologue SBP-1 promotes fatty acid homeostasis by regulating the expression of lipogenic enzymes. We found that, like SBP-1, the C. elegans ARC105 homologue MDT-15 is required for fatty acid homeostasis, and show that both SBP-1 and MDT-15 control transcription of genes governing desaturation of stearic acid to oleic acid. Notably, dietary addition of oleic acid significantly rescued various defects of nematodes targeted with RNA interference against sbp-1 and mdt-15, including impaired intestinal fat storage, infertility, decreased size and slow locomotion, suggesting that regulation of oleic acid levels represents a physiologically critical function of SBP-1 and MDT-15. Taken together, our findings demonstrate that ARC105 is a key effector of SREBP-dependent gene regulation and control of lipid homeostasis in metazoans.

The RAG1 and RAG2 proteins catalyze V(D)J recombination and are essential for generation of the diverse repertoire of antigen receptor genes and effective immune responses. RAG2 is composed of a "core" domain that is required for the recombination reaction and a C-terminal nonessential or "non-core" region. Recent evidence has emerged arguing that the non-core region plays a critical regulatory role in the recombination reaction, and mutations in this region have been identified in patients with immunodeficiencies. Here we present the first structural data for the RAG2 protein, using NMR spectroscopy to demonstrate that the C terminus of RAG2 contains a noncanonical PHD finger. All of the non-core mutations of RAG2 that are implicated in the development of immunodeficiencies are located within the PHD finger, at either zinc-coordinating residues or residues adjacent to an alpha-helix on the surface of the domain that participates in binding to the signaling molecules, phosphoinositides. Functional analysis of disease and phosphoinositide-binding mutations reveals novel intramolecular interactions within the non-core region and suggests that the PHD finger adopts two distinct states. We propose a model in which the equilibrium between these states modulates recombination activity. Together, these data identify the PHD finger as a novel and functionally important domain of RAG2.

We describe an efficient NMR triple resonance approach that correlates, at high resolution, protein side-chain and backbone resonances. It relies on the combination of two strategies: joint evolution of aliphatic side-chain proton/carbon coherences using a backbone N-H based HCcoNH reduced dimensionality (RD) experiment and non-uniform sampling (NUS) in two indirect dimensions. A typical data set containing such correlation information can be acquired in 2 days, at very high resolution unfeasible for conventional 4D HCcoNH-TOCSY experiments. The resonances of the aliphatic side-chain protons are unambiguously assigned to their attached carbons through the analysis of the 'sum' and 'difference' spectra. This approach circumvents the tedious process of manual resonance assignments using HCcH-TOCSY data, while providing additional resolving power of backbone N-H signals. A simple peak-list based algorithm has been implemented in the IBIS software for rapid automated backbone and side-chain assignments.

We have developed a tool for computer-assisted assignments of protein NMR spectra from triple resonance data. The program is designed to resemble established manual assignment procedures as closely as possible. IBIS exports its results in XEASY format. Thus, using IBIS the operator has continuous visual and accounting control over the progress of the assignment procedure. IBIS achieves complete assignments for those residues that exhibit sequential triple resonance connectivities within a few hours or days.

Bacteriophage T7 DNA primase recognizes 5'-GTC-3' in single-stranded DNA. The primase contains a single Cys4 zinc-binding motif that is essential for recognition. Biochemical and mutagenic analyses suggest that the Cys4 motif contacts cytosine of 5'-GTC-3' and may also contribute to thymine recognition. Residues His33 and Asp31 are critical for these interactions. Biochemical analysis also reveals that T7 primase selectively binds CTP in the absence of DNA. We propose that bound CTP selects the remaining base G, of 5'-GTC-3', by base pairing. Our deduced mechanism for recognition of ssDNA by Cys4 motifs bears little resemblance to the recognition of trinucleotides of double-stranded DNA by Cys2His2 zinc fingers.

Amino acids in the serine proteinase inhibitor eglin c important for its inhibitory specificity and activity have been investigated by site-directed mutagenesis. The specificity of eglin c could be changed from elastase to trypsin inhibition by the point mutation Leu45----Arg (L45R) in position P1 [nomenclature according to Schechter and Berger (1967) Biochem. Biophys. Res. Commun. 27, 157-162]. Model building studies based on the crystal structure of mutant L45R [Heinz et al. (1991) J. Mol. Biol. 217, 353-371] were used to rationalize this specificity change. Surprisingly, the double mutant L45R/D46S was found to be a substrate of trypsin and various other serine proteinases. Multidimensional NMR studies show that wild-type eglin c and the double mutant have virtually identical conformations. In the double mutant L45R/D46S, however, the N-H bond vector of the scissile peptide bond shows a much higher mobility, indicating that the internal rigidity of the binding loop is significantly weakened due to the loss or destabilization of the internal hydrogen bond of the P1' residue. Mutant T44P was constructed to examine the role of a proline in position P2, which is frequently found in serine proteinase inhibitors [Laskowski and Kato (1980) Annu. Rev. Biochem. 49, 593-626]. The mutant remains a potent elastase inhibitor but no longer inhibits subtilisin, which could be explained by model building. Both Arg51 and Arg53, located in the core of the molecule and participating in the hydrogen bonding network with residues in the binding loop to maintain rigidity around the scissile bond, were individually replaced with the shorter but equally charged amino acid lysine. Both mutants showed a decrease in their inhibitory potential. The crystal structure of mutant R53K revealed the loss of two hydrogen bonds between the core and the binding loop of the inhibitor, which are partially restored by a solvent molecule, leading to a decrease in inhibition of elastase by 2 orders of magnitude.

A high-precision solution structure of the elastase inhibitor eglin c was determined by NMR and distance geometry calculations. A large set of 947 nuclear Overhauser (NOE) distance constraints was identified, 417 of which were quantified from two-dimensional NOE spectra at short mixing times. In addition, a large number of homonuclear 1H-1H and heteronuclear 1H-15N vicinal coupling constants were used, and constraints on 42 chi 1 and 38 phi angles were obtained. Structure calculations were carried out using the distance geometry program DG-II. These calculations had a high convergence rate, in that 66 out of 75 calculations converged with maximum residual NOE violations ranging from 0.17 A to 0.47 A. The spread of the structures was characterized with average root mean square deviations () between the structures and a mean structure. To calculate the unbiased toward any single structure, a new procedure was used for structure alignment. A canonical structure was calculated from the mean distances, and all structures were aligned relative to that. Furthermore, an angular order parameter S was defined and used to characterize the spread of structures in torsion angle space. To obtain an accurate estimate of the precision of the structure, the number of calculations was increased until the and the angular order parameters stabilized. This was achieved after approximately 40 calculations. The structure consists of a well-defined core whose backbone deviates from the canonical structure ca. 0.4 A, a disordered N-terminal heptapeptide whose backbone deviates by 0.8-12 A, and a proteinase-binding loop whose backbone deviates up to 3.0 A. Analysis of the angular order parameters and inspection of the structures indicates that a hinge-bending motion of the binding loop may occur in solution. Secondary structures were analyzed by comparison of dihedral angle patterns. The high precision of the structure allows one to identify subtle differences with four crystal structures of eglin c determined in complexes with proteinases.

Structural studies of globular proteins by nmr can be enhanced by the use of isotope enrichment. We have been working with proteins enriched with 15N, and with both 15N and 13C. Due to the isotope enrichment we could assign several large proteins with up to 186 residues and could address structural questions. Furthermore, we can accurately measure heteronuclear and homonuclear vicinal coupling constants. This involves in part multidimensional multiple resonance experiments. This is important for characterization of minor conformational changes caused by mutations. We have also made use of isotope enrichment to study the internal mobility of proteins. We also have developed novel methods for measuring accurately 15N relaxation parameters, in particular transverse relaxation rates. This has led us toward a method for directly mapping spectral density functions of the rotational motions of N-H bond vectors in proteins. The protein systems that are discussed include the unlabeled proteins kistrin and cytochrome c551, and the labeled proteins eglin c, a flavodoxin, and human dihydrofolate reductase.

Sequence-specific nuclear magnetic resonance assignments were obtained for eglin c, a polypeptide inhibitor of the granulocytic proteinases elastase and cathepsin G and some other proteinases. The protein consists of a single polypeptide chain of 70 residues. All proton resonances were assigned except for some labile protons of arginine side chains. The patterns of nuclear Overhauser enhancements and coupling constants and the observation of slow hydrogen exchange were used to characterize the secondary structure of the protein. The results indicate that the solution structure of the free inhibitor is very similar to the crystal structure reported for the same protein in the complex with subtilisin Carlsberg. However, a part of the binding loop seems to have a significantly different conformation in the free protein.

Stereospecific assignments were obtained in the protein Eglin c for beta-methylene protons of 9 out of the 24 residues with AMX symmetry, for three residues with long side chains and for the gamma-methyl groups of 8 out of the 11 valines. This was achieved by analyzing the cross-peak multiplet structures in two-dimensional correlated spectra with spectral simulations and peak fitting, and by quantitative measurements of intraresidue nuclear Overhauser enhancements. In addition to the stereospecific assignments, information on the torsion angles phi and chi 1 was obtained. To estimate inferences of internal motions on this analysis, the effect of uniform averaging within a certain range of the torsion angle chi 1 on cross-peak multiplet structures and on relative intraresidue nuclear Overhauser enhancement is discussed.