XLSY: Extra-Large NMR Spectroscopy.

NMR studies of intrinsically disordered proteins and other complex biomolecular systems require spectra with the highest resolution and dimensionality. An efficient approach, extra-large NMR spectroscopy, is presented for experimental data collection, reconstruction, and handling of very large NMR spectra by a combination of the radial and non-uniform sampling, a new processing algorithm, and rigorous statistical validation. We demonstrate the first high-quality reconstruction of a full seven-dimensional HNCOCACONH and two five-dimensional HACACONH and HN(CA)CONH experiments for a representative intrinsically disordered protein α-synuclein. XLSY will significantly enhance the NMR toolbox in challenging biomolecular studies.

Fast multi-dimensional NMR acquisition and processing using the sparse FFT.

Increasing the dimensionality of NMR experiments strongly enhances the spectral resolution and provides invaluable direct information about atomic interactions. However, the price tag is high: long measurement times and heavy requirements on the computation power and data storage. We introduce sparse fast Fourier transform as a new method of NMR signal collection and processing, which is capable of reconstructing high quality spectra of large size and dimensionality with short measurement times, faster computations than the fast Fourier transform, and minimal storage for processing and handling of sparse spectra. The new algorithm is described and demonstrated for a 4D BEST-HNCOCA spectrum.

Combined NMR and molecular dynamics conformational filter identifies unambiguously dynamic ensembles of Dengue protease NS2B/NS3pro.

The dengue protease NS2B/NS3pro has been reported to adopt either an 'open' or a 'closed' conformation. We have developed a conformational filter that combines NMR with MD simulations to identify conformational ensembles that dominate in solution. Experimental values derived from relaxation parameters for the backbone and methyl side chains were compared with the corresponding back-calculated relaxation parameters of different conformational ensembles obtained from free MD simulations. Our results demonstrate a high prevalence for the 'closed' conformational ensemble while the 'open' conformation is absent, indicating that the latter conformation is most probably due to crystal contacts. Conversely, conformational ensembles in which the positioning of the co-factor NS2B results in a 'partially' open conformation, previously described in both MD simulations and X-ray studies, were identified by our conformational filter. Altogether, we believe that our approach allows for unambiguous identification of true conformational ensembles, an essential step for reliable drug discovery.

1H, 13C and 15N resonance assignment of backbone and IVL-methyl side chain of the S135A mutant NS3pro/NS2B protein of Dengue II virus reveals unique secondary structure features in solution.

The serotype II Dengue (DENV 2) virus is the most prevalent of all four known serotypes. Herein, we present nearly complete 1H, 15N, and 13C backbone and 1H, 13C isoleucine, valine, and leucine methyl resonance assignment of the apo S135A catalytically inactive variant of the DENV 2 protease enzyme folded as a tandem formed between the serine protease domain NS3pro and the cofactor NS2B, as well as the secondary structure prediction of this complex based on the assigned chemical shifts using the TALOS-N software. Our results provide a solid ground for future elucidation of the structure and dynamic of the apo NS3pro/NS2B complex, key for adequate development of inhibitors, and a thorough molecular understanding of their function(s).

Assignment of IVL-Methyl side chain of the ligand-free monomeric human MALT1 paracaspase-IgL3 domain in solution.

Mucosa-associated lymphoid tissue protein 1 (MALT1) plays a key role in adaptive immune responses by modulating specific intracellular signalling pathways that control the development and proliferation of both T and B cells. Dysfunction of these pathways is coupled to the progress of highly aggressive lymphoma as well as to potential development of an array of different immune disorders. In contrast to other signalling mediators, MALT1 is not only activated through the formation of the CBM complex together with the proteins CARMA1 and Bcl10, but also by acting as a protease that cleaves multiple substrates to promote lymphocyte proliferation and survival via the NF-κB signalling pathway. Herein, we present the partial 1H, 13C Ile/Val/Leu-Methyl resonance assignment of the monomeric apo form of the paracaspase-IgL3 domain of human MALT1. Our results provide a solid ground for future elucidation of both the three-dimensional structure and the dynamics of MALT1, key for adequate development of inhibitors, and a thorough molecular understanding of its function(s).

Resolution enhancement in NMR spectra by deconvolution with compressed sensing reconstruction.

NMR spectroscopy is one of the basic tools for molecular structure elucidation. Unfortunately, the resolution of the spectra is often limited by inter-nuclear couplings. The existing workarounds often alleviate the problem by trading it for another deficiency, such as spectral artefacts or difficult sample preparation and, thus, are rarely used. We suggest an approach using the coupling deconvolution in the framework of compressed sensing (CS) spectra processing that leads to a major increase in resolution, sensitivity, and overall quality of NUS reconstruction. A new mathematical description of the decoupling by deconvolution explains the effects of thermal noise and reveals a relation with the underlying assumption of the CS. The gain in resolution and sensitivity for challenging molecular systems is demonstrated for the key HNCA experiment used for protein backbone assignment applied to two large proteins: intrinsically disordered 441-residue Tau and a 509-residue globular bacteriophytochrome fragment. The approach will be valuable in a multitude of chemistry applications, where NMR experiments are compromised by the homonuclear scalar coupling.

Structural characterization and biological function of bivalent binding of CD2AP to intrinsically disordered domain of chikungunya virus nsP3 protein.

Alphavirus nsP3 proteins contain long, intrinsically disordered, hypervariable domains, HVD, which serve as hubs for interaction with many cellular proteins. Here, we have deciphered the mechanism and function of HVD interaction with host factors in alphavirus replication. Using NMR spectroscopy, we show that CHIKV HVD contains two SH3 domain-binding sites. Using an innovative chemical shift perturbation signature approach, we demonstrate that CD2AP interaction with HVD is mediated by its SH3-A and SH3-C domains, and this leaves the SH3-B domain available for interaction with other cellular factor(s). This cooperative interaction with two SH3 domains increases binding affinity to CD2AP and possibly induces long-range allosteric effects in HVD. Our data demonstrate that BIN1, CD2AP and SH3KBP1 play redundant roles in initiation of CHIKV replication. Point mutations in both CHIKV HVD binding sites abolish its interaction with all three proteins, CD2AP, BIN1 and SH3KBP1. This results in strong inhibition of viral replication initiation.

Hyperdimensional NMR spectroscopy with nonlinear sampling.

An approach is described for joint interleaved recording, real-time processing, and analysis of NMR data sets. The method employs multidimensional decomposition to find common information in a set of conventional triple-resonance spectra recorded in the nonlinear sampling mode, and builds a model of hyperdimensional (HD) spectrum. While preserving sensitivity per unit of measurement time and allowing for maximal spectral resolution, the approach reduces data collection time on average by 2 orders of magnitude compared to the conventional method. The 7-10 dimensional HD spectrum, which is represented as a set of deconvoluted 1D vectors, is easy to handle and amenable for automated analysis. The method is exemplified by automated assignment for two protein systems of low and high spectral complexity: ubiquitin (globular, 8 kDa) and zetacyt (naturally disordered, 13 kDa). The collection and backbone assignment of the data sets are achieved in real time after approximately 1 and 10 h, respectively. The approach removes the most critical time bottlenecks in data acquisition and analysis. Thus, it can significantly increase the value of NMR spectroscopy in structural biology, for example, in high-throughput structural genomics applications.

Propagation of dynamic changes in barnase upon binding of barstar: an NMR and computational study.

NMR spectroscopy and computer simulations were used to examine changes in chemical shifts and in dynamics of the ribonuclease barnase that result upon binding to its natural inhibitor barstar. Although the spatial structures of free and bound barnase are very similar, binding results in changes of the dynamics of both fast side-chains, as revealed by (2)H relaxation measurements, and NMR chemical shifts in an extended beta-sheet that is located far from the binding interface. Both side-chain dynamics and chemical shifts are sensitive to variations in the ensemble populations of the inter-converting molecular states, which can escape direct structural observation. Molecular dynamics simulations of free barnase and barnase in complex with barstar, as well as a normal mode analysis of barnase using a Gaussian network model, reveal relatively rigid domains that are separated by the extended beta-sheet mentioned above. The observed changes in NMR parameters upon ligation can thus be rationalized in terms of changes in inter-domain dynamics and in populations of exchanging states, without measurable structural changes. This provides an alternative model for the propagation of a molecular response to ligand binding across a protein that is based exclusively on changes in dynamics.

Binding of intrinsically disordered proteins is not necessarily accompanied by a structural transition to a folded form.

There are a large number of protein domains and even entire proteins, lacking ordered structure under physiological conditions. Intriguingly, a highly flexible, random coil-like conformation is the native and functional state for many proteins known to be involved in cell signaling. An example is a key component of immune signaling, the cytoplasmic region of the T cell receptor zeta subunit. This domain exhibits specific dimerization that is distinct from non-specific aggregation behavior seen in many systems. In this work, we use diffusion and chemical shift mapping NMR data to show that the protein does not undergo a transition between disordered and ordered states upon dimerization. This finding opposes the generally accepted view on the behavior of intrinsically disordered proteins, provides evidence for the existence of specific dimerization interactions for intrinsically disordered protein species and opens a new line of research in this new and quickly developing field.

An oxygen-sensitive toxin-antitoxin system.

The Hha and TomB proteins from Escherichia coli form an oxygen-dependent toxin-antitoxin (TA) system. Here we show that YmoB, the Yersinia orthologue of TomB, and its single cysteine variant [C117S]YmoB can replace TomB as antitoxins in E. coli. In contrast to other TA systems, [C117S]YmoB transiently interacts with Hha (rather than forming a stable complex) and enhances the spontaneous oxidation of the Hha conserved cysteine residue to a -SOxH-containing species (sulfenic, sulfinic or sulfonic acid), which destabilizes the toxin. The nuclear magnetic resonance structure of [C117S]YmoB and the homology model of TomB show that the two proteins form a four-helix bundle with a conserved buried cysteine connected to the exterior by a channel with a diameter comparable to that of an oxygen molecule. The Hha interaction site is located on the opposite side of the helix bundle.

Gated electron transfers and electron pathways in azurin: a NMR dynamic study at multiple fields and temperatures.

Dynamic properties of electron transfer pathways in a small blue copper cupredoxin are explored using an extensive 15N NMR relaxation study of reduced Pseudomonas aeruginosa azurin at four magnetic fields (500-900 MHz) and at two temperatures chosen well below the melting point of the protein. Following a careful model-free analysis, several protein regions with different dynamic regimes are identified. Nanosecond time-scale mobility characterizes various residues of the hydrophobic surface patch believed to mark the natural entry point for electrons, notably the surface-exposed copper-ligand His117. These findings are consistent with a gated electron transfer process according to the "dynamic docking" model. Residues 47-49 along intramolecular pathways of electrons show rigidity that is remarkably conserved when increasing the temperature. Three different conformational exchange processes were observed in the millisecond range, one near the only disulfide bridge in the molecule and two near the copper ion. The latter two processes are consistent with previous data such as crystal structures at various pH values and NMR relaxation dispersion experiments; they may indicate an additional gated electron transfer mechanism at slower time-scales.

Accuracy and robustness of three-way decomposition applied to NMR data.

Three-way decomposition is a very versatile analysis tool with applications in a variety of protein NMR fields. It has been used to extract structural data from 3D NOESYs, to determine relaxation rates in large proteins, to identify ligand binding in screening for lead compounds, and to complement non-uniformly recorded (sparse) spectra. All applications so far concerned experimental data sets; it thus remains to address questions of accuracy and robustness of the method using simulated data where the correct answer is known. Systematic tests are presented for relaxation and NOESY data sets. Mixtures of real and synthetic data are used to allow control of various parameters and comparisons with correct reference data, while working with input that is as realistic as possible. The influence of the following parameters is evaluated: signal-to-noise, overlap of signals and the use of a regularization procedure within the algorithm. The main criteria used for the evaluation are accuracy and precision. It is shown that deterioration of accuracy is indicated by internal checks such as decrease of precision. Both with relaxation data and when interpreting NOESY spectra, three-way decomposition exhibits a robust behavior in situations with severe signal overlap and/or poor signal-to-noise, e.g., by avoiding false positives in the NOE shapes of NOESY decompositions. As a complement to this study, three-way decomposition is compared to other methods that achieve the same type of results.

NMR detection of multiple transitions to low-populated states in azurin.

Transitions to conformational states with very low populations were detected for the reduced blue copper protein azurin from Pseudomonas aeruginosa by applying constant relaxation time CPMG measurements to the backbone (15)N nuclei at three magnetic fields (11.7, 14.1, and 18.8 T) and three temperatures (25.7, 35.4, and 44.8 degrees C). Two exchange processes with different rate constants could be discriminated despite populations of the excited states below 1% and spatial neighborhood of the two processes. The group of (15)N nuclei involved in the faster process exhibits at 44.8 degrees C a forward rate constant of 11.7+/-2.4 s(-1) and a population of the exited state of 0.39+/-0.07%. They surround the aromatic ring of histidine 35 whose protonation state is coupled to the flipping of a neighboring peptide plane. For the slower process, the forward rate constant and population of the exited state at 44.8 degrees C are 4.1+/-0.1 s(-1) and 0.45+/-0.02%, respectively. The residues involved cluster nearby the copper ion, which is separated from the protonation site of histidine 35 by about 8 A, indicating conformational rearrangements involving the copper coordinating loops. The dependence of the equilibrium constant on the temperature is consistent with an enthalpy-dominated transition around the copper, but an entropy-controlled transition near histidine 35. The detection by nuclear magnetic resonance of millisecond to second conformational transitions near the copper ion suggests a low energy-cost rearrangement of the copper-binding site that may be necessary for efficient electron transfer.

Automated analysis of large sets of heteronuclear correlation spectra in NMR-based drug discovery.

Drug discovery procedures based on NMR typically require the analysis of thousands of NMR spectra. For example, in "SAR by NMR", two-dimensional NMR spectra are recorded for a target protein mixed with ligand candidates from a comprehensive library of small molecules and are compared to the corresponding spectrum for the protein alone. We present an automated procedure for the comparative analysis of large sets of heteronuclear single quantum coherence spectra, which is based on three-way decomposition and implemented as the software package MUNIN. In a single step, spectra with differences in the peak positions (indicating ligand binding) and the affected peaks are identified. By omission of peak picking, ad hoc scoring of the quality of doubtful peaks is avoided. The procedure has been tested on the bacterial ribonuclease barnase, with a protein concentration of only 50 microM, using several small molecules including the substrate analogue 3'-GMP. Sets of 51 spectra were processed simultaneously, and it is concluded that spectra with binding ligands can be unambiguously identified from much larger sets of spectra.

Accurate relaxation parameters for large proteins.

The practical applicability, performance, and robustness of three-way decomposition (TWD) for the extraction of relaxation parameters are demonstrated for a large protein with 370 residues, the maltose binding protein. An ordinary set of seven relaxation-modulated (15)N HSQC spectra, recorded at another site, is systematically analyzed. For all 341 assigned backbone amide groups, including 21 pairs and one group of three overlapped peaks, T1 decay values were determined. On isolated peaks, TWD extracts T1 values with systematically lower error bounds compared to conventional tools, although for these simple cases the improvements remain limited. However, in the presence of spectral artifacts, the decrease in errors can become significant, demonstrating the higher robustness of TWD. For about half of the peaks in overlapped regions, the decomposition allowed separation of the signals, yielding significantly different T1 values between overlapping signals. For the rest, similarity of the decay times for the two or three overlapping signals could be confirmed within usually low error bounds. The use of TWD thus leads to a significant increase in the number of accessible relaxation probes in large proteins. With a newly implemented graphical user interface, the application of TWD requires merely a peak list, and thus no additional effort compared to conventional approaches is needed.

Effects of NMR spectral resolution on protein structure calculation.

Adequate digital resolution and signal sensitivity are two critical factors for protein structure determinations by solution NMR spectroscopy. The prime objective for obtaining high digital resolution is to resolve peak overlap, especially in NOESY spectra with thousands of signals where the signal analysis needs to be performed on a large scale. Achieving maximum digital resolution is usually limited by the practically available measurement time. We developed a method utilizing non-uniform sampling for balancing digital resolution and signal sensitivity, and performed a large-scale analysis of the effect of the digital resolution on the accuracy of the resulting protein structures. Structure calculations were performed as a function of digital resolution for about 400 proteins with molecular sizes ranging between 5 and 33 kDa. The structural accuracy was assessed by atomic coordinate RMSD values from the reference structures of the proteins. In addition, we monitored also the number of assigned NOESY cross peaks, the average signal sensitivity, and the chemical shift spectral overlap. We show that high resolution is equally important for proteins of every molecular size. The chemical shift spectral overlap depends strongly on the corresponding spectral digital resolution. Thus, knowing the extent of overlap can be a predictor of the resulting structural accuracy. Our results show that for every molecular size a minimal digital resolution, corresponding to the natural linewidth, needs to be achieved for obtaining the highest accuracy possible for the given protein size using state-of-the-art automated NOESY assignment and structure calculation methods.

A comparison of convex and non-convex compressed sensing applied to multidimensional NMR.

The resolution of multidimensional NMR spectra can be severely limited when regular sampling based on the Nyquist-Shannon theorem is used. The theorem binds the sampling rate with a bandwidth of a sampled signal and thus implicitly creates a dependence between the line width and the time of experiment, often making the latter one very long. Recently, Candès et al. (2006) [25] formulated a non-linear sampling theorem that determines the required number of sampling points to be dependent mostly on the number of peaks in a spectrum and only slightly on the number of spectral points. The result was pivotal for rapid development and broad use of signal processing method called compressed sensing. In our previous work, we have introduced compressed sensing to multidimensional NMR and have shown examples of reconstruction of two-dimensional spectra. In the present paper we discuss in detail the accuracy and robustness of two compressed sensing algorithms: convex (iterative soft thresholding) and non-convex (iteratively re-weighted least squares with local ℓ(0)-norm) in application to two- and three-dimensional datasets. We show that the latter method is in many terms more effective, which is in line with recent works on the theory of compressed sensing. We also present the comparison of both approaches with multidimensional decomposition which is one of the established methods for processing of non-linearly sampled data.