Hyperdimensional NMR spectroscopy with nonlinear sampling

Non-uniform sampling to enhance the performance of compact NMR for characterizing new psychoactive substances

Efficient and robust analytical methods are needed to improve the identification and subsequent regulation of new psychoactive substances (NPS). NMR spectroscopy is a unique method able to determine the structure of small molecules such as NPS even in mixtures. However, high-field NMR analysis is associated with expensive purchase and maintenance costs. For more than a decade, compact NMR spectrometers have changed this paradigm. It was recently shown that a dedicated analytical workflow combining compact NMR and databases could identify the molecular structure of NPS, in spite of the lower spectral dispersion and sensitivity of compact spectrometers. This approach relies on 1H-13C HSQC to both recognize NPS and elucidate the structure of unknown substances. Still, its performance is limited by the need to compromise between resolution and experiment time. Here, we show that this strategy can be significantly improved by implementing non-uniform sampling (NUS) to improve spectral resolution in the 13C dimension of HSQC at no cost in terms of experiment time. Gains in the range of 3 to 4 in resolution are achieved for pure NPS and for a mixture. Finally, 2D HSQC with NUS was applied to improve the identification of NPS with the assistance of databases. The resulting method appears as a useful tool for the characterization of NPS in mixtures, which is essential for forensic laboratories.

Recruitment of trimeric eIF2 by phosphatase non-catalytic subunit PPP1R15B

Regulated protein phosphorylation controls most cellular processes. The protein phosphatase PP1 is the catalytic subunit of many holoenzymes that dephosphorylate serine/threonine residues. How these enzymes recruit their substrates is largely unknown. Here, we integrated diverse approaches to elucidate how the PP1 non-catalytic subunit PPP1R15B (R15B) captures its full trimeric eIF2 substrate. We found that the substrate-recruitment module of R15B is largely disordered with three short helical elements, H1, H2, and H3. H1 and H2 form a clamp that grasps the substrate in a region remote from the phosphorylated residue. A homozygous N423D variant, adjacent to H1, reducing substrate binding and dephosphorylation was discovered in a rare syndrome with microcephaly, developmental delay, and intellectual disability. These findings explain how R15B captures its 125 kDa substrate by binding the far end of the complex relative to the phosphosite to present it for dephosphorylation by PP1, a paradigm of broad relevance.

Experimental NOE, Chemical Shift, and Proline Isomerization Data Provide Detailed Insights into Amelotin Oligomerization

Amelotin is an intrinsically disordered protein (IDP) rich in Pro residues and is involved in hydroxyapatite mineralization. It rapidly oligomerizes under physiological conditions of pH and pressure but reverts to its monomeric IDP state at elevated pressure. We identified a 105-residue segment of the protein that becomes ordered upon oligomerization, and we used pressure-jump NMR spectroscopy to measure long-range NOE contacts that exist exclusively in the oligomeric NMR-invisible state. The kinetics of oligomerization and dissociation were probed at the residue-specific level, revealing that the oligomerization process is initiated in the C-terminal half of the segment. Using pressure-jump NMR, the degree of order in the oligomer at the sites of Pro residues was probed by monitoring changes in cis/trans equilibria relative to the IDP state after long-term equilibration under oligomerizing conditions. Whereas most Pro residues revert to trans in the oligomeric state, Pro-49 favors a cis configuration and three Pro residues retain an unchanged cis fraction, pointing to their local lack of order in the oligomeric state. NOE contacts and secondary 13C chemical shifts in the oligomeric state indicate the presence of an 11-residue α-helix, preceded by a small intramolecular antiparallel β-sheet, with slower formation of long-range intermolecular interactions to N-terminal residues. Although none of the models generated by AlphaFold2 for the amelotin monomer was consistent with experimental data, subunits of a hexamer generated by AlphaFold-Multimer satisfied intramolecular NOE and chemical shift data and may provide a starting point for developing atomic models for the oligomeric state.

Non-uniform sampling of similar NMR spectra and its application to studies of the interaction between alpha-synuclein and liposomes

The accelerated acquisition of multidimensional NMR spectra using sparse non-uniform sampling (NUS) has been widely adopted in recent years. The key concept in NUS is that a major part of the data is omitted during measurement, and then reconstructed using, for example, compressed sensing (CS) methods. CS requires spectra to be compressible, that is, they should contain relatively few "significant" points. The more compressible the spectrum, the fewer experimental NUS points needed in order for it to be accurately reconstructed. In this paper we show that the CS processing of similar spectra can be enhanced by reconstructing only the differences between them. Accurate reconstruction can be obtained at lower sampling levels as the difference is sparser than the spectrum itself. In many situations this method is superior to "conventional" compressed sensing. We exemplify the concept of "difference CS" with one such case-the study of alpha-synuclein binding to liposomes and its dependence on temperature. To obtain information on temperature-dependent transitions between different states, we need to acquire several dozen spectra at various temperatures, with and without the presence of liposomes. Our detailed investigation reveals that changes in the binding modes of the alpha-synuclein ensemble are not only temperature-dependent but also show non-linear behavior in their transitions. Our proposed CS processing approach dramatically reduces the number of NUS points required and thus significantly shortens the experimental time.

All Domains of SARS-CoV-2 nsp1 Determine Translational Shutoff and Cytotoxicity of the Protein

Replication of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) strongly affects cellular metabolism and results in rapid development of the cytopathic effect (CPE). The hallmarks of virus-induced modifications are inhibition of translation of cellular mRNAs and redirection of the cellular translational machinery to the synthesis of virus-specific proteins. The multifunctional nonstructural protein 1 (nsp1) of SARS-CoV-2 is a major virulence factor and a key contributor to the development of translational shutoff. In this study, we applied a wide range of virological and structural approaches to further analyze nsp1 functions. The expression of this protein alone was found to be sufficient to cause CPE. However, we selected several nsp1 mutants exhibiting noncytopathic phenotypes. The attenuating mutations were detected in three clusters, located in the C-terminal helices, in one of the loops of the structured domain and in the junction of the disordered and structured fragment of nsp1. NMR-based analysis of the wild type nsp1 and its mutants did not confirm the existence of a stable β5-strand that was proposed by the X-ray structure. In solution, this protein appears to be present in a dynamic conformation, which is required for its functions in CPE development and viral replication. The NMR data also suggest a dynamic interaction between the N-terminal and C-terminal domains. The identified nsp1 mutations make this protein noncytotoxic and incapable of inducing translational shutoff, but they do not result in deleterious effects on viral cytopathogenicity. IMPORTANCE The nsp1 of SARS-CoV-2 is a multifunctional protein that modifies the intracellular environment for the needs of viral replication. It is responsible for the development of translational shutoff, and its expression alone is sufficient to cause a cytopathic effect (CPE). In this study, we selected a wide range of nsp1 mutants exhibiting noncytopathic phenotypes. The attenuating mutations, clustered in three different fragments of nsp1, were extensively characterized via virological and structural methods. Our data strongly suggest interactions between the nsp1 domains, which are required for the protein's functions in CPE development. Most of the mutations made nsp1 noncytotoxic and incapable of inducing translational shutoff. Most of them did not affect the viability of the viruses, but they did decrease the rates of replication in cells competent in type I IFN induction and signaling. These mutations, and their combinations, in particular, can be used for the development of SARS-CoV-2 variants with attenuated phenotypes.

Juxtaposition of Bub1 and Cdc20 on phosphorylated Mad1 during catalytic mitotic checkpoint complex assembly

In response to improper kinetochore-microtubule attachments in mitosis, the spindle assembly checkpoint (SAC) assembles the mitotic checkpoint complex (MCC) to inhibit the anaphase-promoting complex/cyclosome, thereby delaying entry into anaphase. The MCC comprises Mad2:Cdc20:BubR1:Bub3. Its assembly is catalysed by unattached kinetochores on a Mad1:Mad2 platform. Mad1-bound closed-Mad2 (C-Mad2) recruits open-Mad2 (O-Mad2) through self-dimerization. This interaction, combined with Mps1 kinase-mediated phosphorylation of Bub1 and Mad1, accelerates MCC assembly, in a process that requires O-Mad2 to C-Mad2 conversion and concomitant binding of Cdc20. How Mad1 phosphorylation catalyses MCC assembly is poorly understood. Here, we characterized Mps1 phosphorylation of Mad1 and obtained structural insights into a phosphorylation-specific Mad1:Cdc20 interaction. This interaction, together with the Mps1-phosphorylation dependent association of Bub1 and Mad1, generates a tripartite assembly of Bub1 and Cdc20 onto the C-terminal domain of Mad1 (Mad1CTD). We additionally identify flexibility of Mad1:Mad2 that suggests how the Cdc20:Mad1CTD interaction brings the Mad2-interacting motif (MIM) of Cdc20 near O-Mad2. Thus, Mps1-dependent formation of the MCC-assembly scaffold functions to position and orient Cdc20 MIM near O-Mad2, thereby catalysing formation of C-Mad2:Cdc20.

Linear discriminant analysis reveals hidden patterns in NMR chemical shifts of intrinsically disordered proteins

NMR spectroscopy is key in the study of intrinsically disordered proteins (IDPs). Yet, even the first step in such an analysis-the assignment of observed resonances to particular nuclei-is often problematic due to low peak dispersion in the spectra of IDPs. We show that the assignment process can be aided by finding "hidden" chemical shift patterns specific to the amino acid residue types. We find such patterns in the training data from the Biological Magnetic Resonance Bank using linear discriminant analysis, and then use them to classify spin systems in an α-synuclein sample prepared by us. We describe two situations in which the procedure can greatly facilitate the analysis of NMR spectra. The first involves the mapping of spin systems chains onto the protein sequence, which is part of the assignment procedure-a prerequisite for any NMR-based protein analysis. In the second, the method supports assignment transfer between similar samples. We conducted experiments to demonstrate these cases, and both times the majority of spin systems could be unambiguously assigned to the correct residue types.

Bacterial divisome protein FtsA forms curved antiparallel double filaments when binding to FtsN

During bacterial cell division, filaments of tubulin-like FtsZ form the Z-ring, which is the cytoplasmic scaffold for divisome assembly. In Escherichia coli, the actin homologue FtsA anchors the Z-ring to the membrane and recruits divisome components, including bitopic FtsN. FtsN regulates the periplasmic peptidoglycan synthase FtsWI. To characterize how FtsA regulates FtsN, we applied electron microscopy to show that E. coli FtsA forms antiparallel double filaments on lipid monolayers when bound to the cytoplasmic tail of FtsN. Using X-ray crystallography, we demonstrate that Vibrio maritimus FtsA crystallizes as an equivalent double filament. We identified an FtsA-FtsN interaction site in the IA-IC interdomain cleft of FtsA using X-ray crystallography and confirmed that FtsA forms double filaments in vivo by site-specific cysteine cross-linking. FtsA-FtsN double filaments reconstituted in or on liposomes prefer negative Gaussian curvature, like those of MreB, the actin-like protein of the elongasome. We propose that curved antiparallel FtsA double filaments together with treadmilling FtsZ filaments organize septal peptidoglycan synthesis in the division plane.

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 ¹H, ¹³C Ile/Val/Leu-Methyl resonance assignment of the monomeric apo form of the paracaspase-IgL₃ 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).

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 ¹H, ¹⁵N, and ¹³C backbone and ¹H, ¹³C 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).

Dynamics in Fip1 regulate eukaryotic mRNA 3' end processing

In this study, Kumar et al. characterized the structure–function relationship of the essential poly(A) factor Fip1. Using in vitro reconstitution and structural studies, the authors report that Fip1 dynamics within the 3′ end processing machinery are required to coordinate cleavage and polyadenylation.

Cleavage and polyadenylation factor (CPF/CPSF) is a multiprotein complex essential for mRNA 3′ end processing in eukaryotes. It contains an endonuclease that cleaves pre-mRNAs, and a polymerase that adds a poly(A) tail onto the cleaved 3′ end. Several CPF subunits, including Fip1, contain intrinsically disordered regions (IDRs). IDRs within multiprotein complexes can be flexible, or can become ordered upon interaction with binding partners. Here, we show that yeast Fip1 anchors the poly(A) polymerase Pap1 onto CPF via an interaction with zinc finger 4 of another CPF subunit, Yth1. We also reconstitute a fully recombinant 850-kDa CPF. By incorporating selectively labeled Fip1 into recombinant CPF, we could study the dynamics of Fip1 within the megadalton complex using nuclear magnetic resonance (NMR) spectroscopy. This reveals that a Fip1 IDR that connects the Yth1- and Pap1-binding sites remains highly dynamic within CPF. Together, our data suggest that Fip1 dynamics within the 3′ end processing machinery are required to coordinate cleavage and polyadenylation.

Multiplexing experiments in NMR and multi-nuclear MRI

Multiplexing NMR experiments by direct detection of multiple free induction decays (FIDs) in a single experiment offers a dramatic increase in the spectral information content and often yields significant improvement in sensitivity per unit time. Experiments with multi-FID detection have been designed with both homonuclear and multinuclear acquisition, and the advent of multiple receivers on commercial spectrometers opens up new possibilities for recording spectra from different nuclear species in parallel. Here we provide an extensive overview of such techniques, designed for applications in liquid- and solid-state NMR as well as in hyperpolarized samples. A brief overview of multinuclear MRI is also provided, to stimulate cross fertilization of ideas between the two areas of research (NMR and MRI). It is shown how such techniques enable the design of experiments that allow structure elucidation of small molecules from a single measurement. Likewise, in biomolecular NMR experiments multi-FID detection allows complete resonance assignment in proteins. Probes with multiple RF microcoils routed to multiple NMR receivers provide an alternative way of increasing the throughput of modern NMR systems, effectively reducing the cost of NMR analysis and increasing the information content at the same time. Solid-state NMR experiments have also benefited immensely from both parallel and sequential multi-FID detection in a variety of multi-dimensional pulse schemes. We are confident that multi-FID detection will become an essential component of future NMR methodologies, effectively increasing the sensitivity and information content of NMR measurements.

Facilitated structure verification of the biopharmaceutical peptide exenatide by 2D heteronuclear NMR maps

Exenatide is a peptide based anti-diabetic prescription medication. Until now, the literature has lacked a comprehensive atom-specific molecular characterization for this complex large peptide by NMR spectroscopy that can be effortlessly and rapidly utilized for biopharmaceutical structural veracity. Peptide structure verification by NMR is challenging and cumbersome when reliant on traditional proton-based methodology (through-bond and through-space proton connectivity) alone due to increasing complexity, low signal dispersion, and overlap. These challenges are overcome by using 2D heteronuclear (¹H-¹³C and ¹H-¹⁵N) maps that not only allow unambiguous signal assignment, but also condense the structural verification information within simplified peptide amide and carbon fingerprint maps. Here we report such simplified amide and carbon fingerprint maps for exenatide; made possible by the first ever comprehensive heteronuclear (¹H,¹³C, and ¹⁵N) atom specific assignment of exenatide. These heteronuclear assignments were obtained without any isotopic enrichments i.e. at natural abundance, and hence are easily deployable as routine procedures. Furthermore, we compare the 2D heteronuclear maps of exenatide to a chemically identical peptide differing only in the isomerism of the Cα position of the first amino acid, [dHis1]-exenatide, to demonstrate the uniqueness of these maps. We show that despite deliberate changes in pH, temperature, and concentrations, the differences between the amide maps of exenatide and [dHis1]-exenatide are retained. The work presented here not only provides a facilitated structure verification of exenatide but also a framework for heteronuclear NMR data acquisition and signal assignment of large peptides, at natural abundance, in creating their respective unique 2D fingerprint maps.

Software for reconstruction of nonuniformly sampled NMR data

Nonuniform sampling (NUS) of multidimensional NMR experiments is a powerful tool to obtain high-resolution spectra with less instrument time. With NUS, only a subset of the points needed for conventional Fourier transformation is recorded, and sophisticated algorithms are needed to reconstruct the missing data points. During the last decade, several software packages implementing the reconstruction algorithms have emerged and been refined and now result in spectra of almost similar quality as spectra from conventionally recorded and processed data. However, from the number of literature references to the reconstruction methods, many more multidimensional NMR spectra could presumably be recorded with NUS. To help researchers considering to start using NUS for their NMR experiments, we here review 13 different reconstruction methods found in five software packages (CambridgeCS, hmsIST, MddNMR, NESTA-NMR, and SMILE). We have compared how the methods run with the provided example scripts for reconstructing a nonuniform sampled three-dimensional ¹⁵ N-NOESY-HSQC at sampling densities from 5% to 50%. Overall, the spectra are all of similar quality above 20% sampling density. Thus, without any particular knowledge about the details of the reconstruction algorithms, significant reduction in the experiment time can be achieved. Below 20% sampling density, the intensities of particular weak peaks start being affected. MddNMR's IST with virtual echo and the SMILE algorithms still reproduced the spectra with the highest accuracy of peak intensities.

Novel NMR Assignment Strategy Reveals Structural Heterogeneity in Solution of the nsP3 HVD Domain of Venezuelan Equine Encephalitis Virus

In recent years, intrinsically disordered proteins (IDPs) and disordered domains have attracted great attention. Many of them contain linear motifs that mediate interactions with other factors during formation of multicomponent protein complexes. NMR spectrometry is a valuable tool for characterizing this type of interactions on both amino acid (aa) and atomic levels. Alphaviruses encode a nonstructural protein nsP3, which drives viral replication complex assembly. nsP3 proteins contain over 200-aa-long hypervariable domains (HVDs), which exhibits no homology between different alphavirus species, are predicted to be intrinsically disordered and appear to be critical for alphavirus adaptation to different cells. Previously, we have shown that nsP3 HVD of chikungunya virus (CHIKV) is completely disordered with low tendency to form secondary structures in free form. In this new study, we used novel NMR approaches to assign the spectra for the nsP3 HVD of Venezuelan equine encephalitis virus (VEEV). The HVDs of CHIKV and VEEV have no homology but are both involved in replication complex assembly and function. We have found that VEEV nsP3 HVD is also mostly disordered but contains a short stable α-helix in its C-terminal fragment, which mediates interaction with the members of cellular Fragile X syndrome protein family. Our NMR data also suggest that VEEV HVD has several regions with tendency to form secondary structures.

Reducing the measurement time of exact NOEs by non-uniform sampling

We have previously reported on the measurement of exact NOEs (eNOEs), which yield a wealth of additional information in comparison to conventional NOEs. We have used these eNOEs in a variety of applications, including calculating high-resolution structures of proteins and RNA molecules. The collection of eNOEs is challenging, however, due to the need to measure a NOESY buildup series consisting of typically four NOESY spectra with varying mixing times in a single measurement session. While the 2D version can be completed in a few days, a fully sampled 3D-NOESY buildup series can take 10 days or more to acquire. This can be both expensive as well as problematic in the case of samples that are not stable over such a long period of time. One potential method to significantly decrease the required measurement time of eNOEs is to use non-uniform sampling (NUS) to decrease the number of points measured in the indirect dimensions. The effect of NUS on the extremely tight distance restraints extracted from eNOEs may be very pronounced. Therefore, we investigated the fidelity of eNOEs measured from three test cases at decreasing NUS densities: the 18.4 kDa protein human Pin1, the 4.1 kDa WW domain of Pin1 (both in 3D), and a 4.6 kDa 14mer RNA UUCG tetraloop (2D). Our results show that NUS imparted negligible error on the eNOE distances derived from good quality data down to 10% sampling for all three cases, but there is a noticeable decrease in the eNOE yield that is dependent upon the underlying sparsity, and thus complexity, of the sample. For Pin1, this transition occurred at roughly 40% while for the WW domain and the UUCG tetraloop it occurred at lower NUS densities of 20% and 10%, respectively. We rationalized these numbers through reconstruction simulations under various conditions. The extent of this loss depends upon the number of scans taken as well as the number of peaks to be reconstructed. Based on these findings, we have created guidelines for choosing an optimal NUS density depending on the number of peaks needed to be reconstructed in the densest region of a 2D or 3D NOESY spectrum.

Assigned NMR backbone resonances of the ligand-binding region domain of the pneumococcal serine-rich repeat protein (PsrP-BR) reveal a rigid monomer in solution

The pneumococcal serine rich repeat protein (PsrP) is displayed on the surface of Streptococcus pneumoniae with a suggested role in colonization in the human upper respiratory tract. Full-length PsrP is a 4000 residue-long multi-domain protein comprising a positively charged functional binding region (BR) domain for interaction with keratin and extracellular DNA during pneumococcal adhesion and biofilm formation, respectively. The previously determined crystal structure of the BR domain revealed a flat compressed barrel comprising two sides with an extended β-sheet on one side, and another β-sheet that is distorted by loops and β-turns on the other side. Crystallographic B-factors indicated a relatively high mobility of loop regions that were hypothesized to be important for binding. Furthermore, the crystal structure revealed an inter-molecular β-sheet formed between edge strands of two symmetry-related molecules, which could promote bacterial aggregation during biofilm formation. Here we report the near complete ¹⁵N/¹³C/¹H backbone resonance assignment of the BR domain of PsrP, revealing a secondary structure profile that is almost identical to the X-ray structure. Dynamic ¹⁵N-T₁, T₂ and NOE data suggest a monomeric and rigid structure of BR with disordered residues only at the N- and C-termini. The presented peak assignment will allow us to identify BR residues that are crucial for ligand binding.

Structural Basis for GTP versus ATP Selectivity in the NMP Kinase AK3

ATP and GTP are exceptionally important molecules in biology with multiple, and often discrete, functions. Therefore, enzymes that bind to either of them must develop robust mechanisms to selectively utilize one or the other. Here, this specific problem is addressed by molecular studies of the human NMP kinase AK3, which uses GTP to phosphorylate AMP. AK3 plays an important role in the citric acid cycle, where it is responsible for GTP/GDP recycling. By combining a structural biology approach with functional experiments, we present a comprehensive structural and mechanistic understanding of the enzyme. We discovered that AK3 functions by recruitment of GTP to the active site, while ATP is rejected and nonproductively bound to the AMP binding site. Consequently, ATP acts as an inhibitor with respect to GTP and AMP. The overall features with specific recognition of the correct substrate and nonproductive binding by the incorrect substrate bear a strong similarity to previous findings for the ATP specific NMP kinase adenylate kinase. Taken together, we are now able to provide the fundamental principles for GTP and ATP selectivity in the large NMP kinase family. As a side-result originating from nonlinearity of chemical shifts in GTP and ATP titrations, we find that protein surfaces offer a general and weak binding affinity for both GTP and ATP. These nonspecific interactions likely act to lower the available intracellular GTP and ATP concentrations and may have driven evolution of the Michaelis constants of NMP kinases accordingly.

Mapping protein-polymer conformations in bioconjugates with atomic precision

Rational design of protein-polymer bioconjugates is hindered by limited experimental data and mechanistic understanding on interactions between the two. In this communication, nuclear magnetic resonance (NMR) paramagnetic relaxation enhancement (PRE) reports on distances between paramagnetic spin labels and NMR active nuclei, informing on the conformation of conjugated polymers. 1H/15N-heteronuclear single quantum coherence (HSQC) NMR spectra were collected for ubiquitin (Ub) modified with block copolymers incorporating spin labels at different positions along their backbone. The resultant PRE data show that the conjugated polymers have conformations biased towards the nonpolar β-sheet face of Ub, rather than behaving as if in solution. The bioconjugates are stabilized against denaturation by guanidine-hydrochloride, as measured by circular dichroism (CD), and this stabilization is attributed to the interaction between the protein and conjugated polymer.

Escherichia coli O176 LPS structure and dynamics: A NMR spectroscopy and MD simulation study

A lipopolysaccharide (LPS) molecule is a key component of the bacterial outer membrane used to protect the bacterium and to interact with the environment. To gain insight into its function, the study of the LPS conformation and dynamics at the molecular and cellular levels is necessary, but these highly diverse and dynamic membrane-LPS systems are difficult to study. In this work, by using NMR spectroscopy and molecular dynamics (MD) simulations, we determined the conformational preferences of an E. coli O176 O-antigen polysaccharide at the atomic level. Moreover, we analyzed the use of non-uniform sampling (NUS) for the acquisition of high dynamic range spectra, like ¹H,¹H-NOESY NMR experiments. A comparison of the effective transglycosidic distances derived from conventional uniformly sampled and NUS ¹H,¹H-NOESY data showed high similarity under equal measuring time conditions. Furthermore, the experimentally derived internuclear distances of the O-antigen polysaccharide with ten repeating units (RUs) showed very good agreement to those calculated from the MD simulations of the same O-antigen polysaccharide in solution. Analysis of the LPS bilayer simulations with five and with ten RUs revealed that, although similar with respect to populated states in solution, the O-antigen in LPS bilayers had more extended chains as a result of spatial limitations due to close packing. Additional MD simulations of O-antigen polysaccharides from E. coli O6 (branched repeating unit) and O91 (negatively charged linear repeating unit) in solution and LPS bilayers were performed and compared to those of O176 (linear polymer). For all three O-antigens, the ensemble of structures present for the polysaccharides in solution were consistent with the results from their ¹H,¹H-NOESY experiments. In addition, the similarities between the O-antigen on its own and as a constituent of the full LPS in bilayer environment makes it possible to realistically describe the LPS conformation and dynamics from the MD simulations.

•

Uniform and non-uniform sampled NOESY NMR data yield similar internuclear distances.

•

O-antigen internuclear distances from NMR and MD show excellent agreement.

•

O-antigen ensemble structures from MD are consistent with NMR observations.

•

O-antigen structures are more extended in LPS bilayers than in solution.

•

MD simulations can describe realistic LPS conformation and dynamics.

Rapid Biophysical Characterization and NMR Spectroscopy Structural Analysis of Small Proteins from Bacteria and Archaea

Proteins encoded by small open reading frames (sORFs) have a widespread occurrence in diverse microorganisms and can be of high functional importance. However, due to annotation biases and their technically challenging direct detection, these small proteins have been overlooked for a long time and were only recently rediscovered. The currently rapidly growing number of such proteins requires efficient methods to investigate their structure-function relationship. Herein, a method is presented for fast determination of the conformational properties of small proteins. Their small size makes them perfectly amenable for solution-state NMR spectroscopy. NMR spectroscopy can provide detailed information about their conformational states (folded, partially folded, and unstructured). In the context of the priority program on small proteins funded by the German research foundation (SPP2002), 27 small proteins from 9 different bacterial and archaeal organisms have been investigated. It is found that most of these small proteins are unstructured or partially folded. Bioinformatics tools predict that some of these unstructured proteins can potentially fold upon complex formation. A protocol for fast NMR spectroscopy structure elucidation is described for the small proteins that adopt a persistently folded structure by implementation of new NMR technologies, including automated resonance assignment and nonuniform sampling in combination with targeted acquisition.

Enhancing Compression Level for More Efficient Compressed Sensing and Other Lessons from NMR Spectroscopy

Modern nuclear magnetic resonance spectroscopy (NMR) is based on two- and higher-dimensional experiments that allow the solving of molecular structures, i.e., determine the relative positions of single atoms very precisely. However, rich chemical information comes at the price of long data acquisition times (up to several days). This problem can be alleviated by compressed sensing (CS)—a method that revolutionized many fields of technology. It is known that CS performs the most efficiently when measured objects feature a high level of compressibility, which in the case of NMR signal means that its frequency domain representation (spectrum) has a low number of significant points. However, many NMR spectroscopists are not aware of the fact that various well-known signal acquisition procedures enhance compressibility and thus should be used prior to CS reconstruction. In this study, we discuss such procedures and show to what extent they are complementary to CS approaches. We believe that the survey will be useful not only for NMR spectroscopists but also to inspire the broader signal processing community.

Fast time-resolved NMR with non-uniform sampling

NMR spectroscopy is a versatile tool for studying time-dependent processes: chemical reactions, phase transitions or macromolecular structure changes. However, time-resolved NMR is usually based on the simplest among available techniques - one-dimensional spectra serving as "snapshots" of the studied process. One of the reasons is that multidimensional experiments are very time-expensive due to costly sampling of evolution time space. In this review we summarize efforts to alleviate the problem of limited applicability of multidimensional NMR in time-resolved studies. We focus on techniques based on sparse or non-uniform sampling (NUS), which lead to experimental time reduction by omitting a significant part of the data during measurement and reconstructing it mathematically, adopting certain assumptions about the spectrum. NUS spectra are faster to acquire than conventional ones and thus better suited to the role of "snapshots", but still suffer from non-stationarity of the signal i.e. amplitude and frequency variations within a dataset. We discuss in detail how these instabilities affect the spectra, and what are the optimal ways of sampling the non-stationary FID signal. Finally, we discuss related areas of NMR where serial experiments are exploited and how they can benefit from the same NUS-based approaches.

Latest trends in structure based drug design with protein targets

Structure based drug designing is an important endeavor in the field of structural bioinformatics. Previously the entire process was dependent on the wet-lab experiments to build libraries of ligand molecules. And the molecules used to be tested to determine their binding efficacies with protein target. However, the entire process is very lengthy and above all highly expensive. With the advent of supercomputers and increasing computational powers, the search process for finding suitable ligand molecules against target proteins have become more streamlined and cost-effective. Now the entire ligand search process is performed in-silico with the help of the techniques of virtual screening, molecular docking simulations and molecular dynamics studies. In the present chapter, a brief overview of the computational techniques involved in structure based drug designing is presented with a special emphasis on the thermodynamic principles behind the molecular interactions.

Investigation of backbone dynamics and local geometry of bio-molecules using calculated NMR chemical shifts and anisotropies

Prerequisite for chemical shift (CS) and CS tensor calculations are highly refined structures defining the molecular surroundings of the nuclei under study. Here, we present geometry optimizations with ¹³C and ¹⁵N CS constraints for large bio-molecules like peptides and proteins. The method discussed here provides both, refined structures and chemical shift tensors. Furthermore, since the experimental resonances of aligned systems are related to CS tensors, they strongly depend on the orientation and motion of molecules, their fragments, functional groups and moieties. For efficient CS calculations we apply a semi-empirical approach-the bond polarization theory (BPT). The BPT relies on linear bond polarization parameters and we present a new set of parameters based on ab initio second-order Møller-Plesset perturbation theory calculations. The new parametrization extends the applicability of the BPT approach to a wide range of organic molecules and bio-polymers. Here, the method has been applied to the protein ubiquitin and the membrane-active peptide gramicidin A (dimer) in oriented bilayers. The calculated ¹³C and ¹⁵N CS values of best-refined structures published until now gave a large scatter with respect to the experiment. It will be shown that BPT CS optimizations can reduce these errors to values near the experimental uncertainty. In combination with molecular dynamics with orientational constraints it is possible to study motional dynamics and BPT calculations can provide residual chemical shift anisotropies.

Modulation of Structural Heterogeneity Controls Phytochrome Photoswitching

Phytochromes sense red/far-red light and control many biological processes in plants, fungi, and bacteria. Although the crystal structures of dark- and light-adapted states have been determined, the molecular mechanisms underlying photoactivation remain elusive. Here, we demonstrate that the conserved tongue region of the PHY domain of a 57-kDa photosensory module of Deinococcus radiodurans phytochrome changes from a structurally heterogeneous dark state to an ordered, light-activated state. The results were obtained in solution by utilizing a laser-triggered activation approach detected on the atomic level with high-resolution protein NMR spectroscopy. The data suggest that photosignaling of phytochromes relies on careful modulation of structural heterogeneity of the PHY tongue.

Crystal structure of the putative cyclase IdmH from the indanomycin nonribosomal peptide synthase/polyketide synthase

Indanomycin is biosynthesized by a hybrid nonribosomal peptide synthase/polyketide synthase (NRPS/PKS) followed by a number of 'tailoring' steps to form the two ring systems that are present in the mature product. It had previously been hypothesized that the indane ring of indanomycin was formed by the action of IdmH using a Diels-Alder reaction. Here, the crystal structure of a selenomethionine-labelled truncated form of IdmH (IdmH-Δ99-107) was solved using single-wavelength anomalous dispersion (SAD) phasing. This truncated variant allows consistent and easy crystallization, but importantly the structure was used as a search model in molecular replacement, allowing the full-length IdmH structure to be determined to 2.7 Å resolution. IdmH is a homodimer, with the individual protomers consisting of an α+β barrel. Each protomer contains a deep hydrophobic pocket which is proposed to constitute the active site of the enzyme. To investigate the reaction catalysed by IdmH, 88% of the backbone NMR resonances were assigned, and using chemical shift perturbation of [¹⁵N]-labelled IdmH it was demonstrated that indanomycin binds in the active-site pocket. Finally, combined quantum mechanical/molecular mechanical (QM/MM) modelling of the IdmH reaction shows that the active site of the enzyme provides an appropriate environment to promote indane-ring formation, supporting the assignment of IdmH as the key Diels-Alderase catalysing the final step in the biosynthesis of indanomycin through a similar mechanism to other recently characterized Diels-Alderases involved in polyketide-tailoring reactions. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at https://proteopedia.org/w/Journal:IUCrJ:S2052252519012399.

Nuclear magnetic resonance (NMR)-based metabolomics for cancer research

Nuclear magnetic resonance (NMR) has emerged as an effective tool in various spheres of biomedical research, amongst which metabolomics is an important method for the study of various types of disease. Metabolomics has proved its stronghold in cancer research by the development of different NMR methods over time for the study of metabolites, thus identifying key players in the aetiology of cancer. A plethora of one-dimensional and two-dimensional NMR experiments (in solids, semi-solids and solution phases) are utilized to obtain metabolic profiles of biofluids, cell extracts and tissue biopsy samples, which can further be subjected to statistical analysis. Any alteration in the assigned metabolite peaks gives an indication of changes in metabolic pathways. These defined changes demonstrate the utility of NMR in the early diagnosis of cancer and provide further measures to combat malignancy and its progression. This review provides a snapshot of the trending NMR techniques and the statistical analysis involved in the metabolomics of diseases, with emphasis on advances in NMR methodology developed for cancer research.

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.

Experiments with direct detection of multiple FIDs

Pulse schemes with direct observation of multiple free induction decays (FIDs) offer a dramatic increase in the spectral information content of NMR experiments and often yield substantial improvement in measurement sensitivity per unit time. Availability of multiple receivers on the state-of-the-art commercial spectrometers allows spectra from different nuclear species to be recorded in parallel routinely. Experiments with multi-FID detection have been designed with both, homonuclear and multinuclear acquisition. We provide a brief overview of such techniques designed for applications in liquid- and solid- state NMR as well as in hyperpolarized samples. Here we show how these techniques have led to design of experiments that allow structure elucidation of small molecules and resonance assignment in proteins from a single measurement. Probes with multiple RF micro-coils routed to multiple NMR receivers provide an alternative way of increasing the throughput of modern NMR systems. Solid-state NMR experiments have also benefited immensely from both parallel and simultaneous FID acquisition in a variety of multi-dimensional pulse schemes. We believe that multi-FID detection will become an essential component of the future NMR methodologies effectively increasing the information content of NMR experiments and reducing the cost of NMR analysis.

Backbone 1H, 13C, and 15N resonance assignments of BoMan26A, a β-mannanase of the glycoside hydrolase family 26 from the human gut bacterium Bacteroides ovatus

Bacteroides ovatus is a member of the human gut microbiota. The importance of this microbial consortium involves the degradation of complex dietary glycans mainly conferred by glycoside hydrolases. In this study we focus on one such catabolic glycoside hydrolase from B. ovatus. The enzyme, termed BoMan26A, is a β-mannanase that takes part in the hydrolytic degradation of galactomannans. The crystal structure of BoMan26A has previously been determined to reveal a TIM-barrel like fold, but the relation between the protein structure and the mode of substrate processing has not yet been studied. Here we report residue-specific assignments for 95% of the 344 backbone amides of BoMan26A. The assignments form the basis for future studies of the relationship between substrate interactions and protein dynamics. In particular, the potential role of loops adjacent to glycan binding sites is of interest for such studies.

Monitoring Hydrogenation Reactions using Benchtop 2D NMR with Extraordinary Sensitivity and Spectral Resolution

Low-field benchtop nuclear magnetic resonance (BT-NMR) spectrometers with Halbach magnets are being increasingly used in science and industry as cost-efficient tools for the monitoring of chemical reactions, including hydrogenation. However, their use of low-field magnets limits both resolution and sensitivity. In this paper, we show that it is possible to alleviate these two problems through the combination of parahydrogen-induced polarization (PHIP) and fast correlation spectroscopy with time-resolved non-uniform sampling (TR-NUS). PHIP can enhance NMR signals so that substrates are easily detectable on BT-NMR spectrometers. The interleaved acquisition of one- and two-dimensional spectra with TR-NUS provides unique insight into the consecutive moments of hydrogenation reactions, with a spectral resolution unachievable in a standard approach. We illustrate the potential of the technique with two examples: the hydrogenation of ethylphenyl propiolate and the hydrogenation of a mixture of two substrates - ethylphenyl propiolate and ethyl 2-butynoate.

TReNDS-Software for reaction monitoring with time-resolved non-uniform sampling

NMR spectroscopy, used routinely for structure elucidation, has also become a widely applied tool for process and reaction monitoring. However, the most informative of NMR methods-correlation experiments-are often useless in this kind of applications. The traditional sampling of a multidimensional FID is usually time-consuming, and thus, the reaction-monitoring toolbox was practically limited to 1D experiments (with rare exceptions, e.g., single-scan or fast-sampling experiments). Recently, the technique of time-resolved non-uniform sampling (TR-NUS) has been proposed, which allows to use standard multidimensional pulse sequences preserving the temporal resolution close to that achievable in 1D experiments. However, the method existed only as a prototype and did not allow on-the-fly processing during the reaction. In this paper, we introduce TReNDS: free, user-friendly software kit for acquisition and processing of TR-NUS data. The program works on Bruker, Agilent, and Magritek spectrometers, allowing to carry out up to four experiments with interleaved TR-NUS. The performance of the program is demonstrated on the example of enzymatic hydrolysis of sucrose.

Measurement of protein backbone 13CO and 15N relaxation dispersion at high resolution

Peak overlap in crowded regions of two-dimensional spectra prevents characterization of dynamics for many sites of interest in globular and intrinsically disordered proteins. We present new three-dimensional pulse sequences for measurement of Carr-Purcell-Meiboom-Gill relaxation dispersions at backbone nitrogen and carbonyl positions. To alleviate increase in the measurement time associated with the additional spectral dimension, we use non-uniform sampling in combination with two distinct methods of spectrum reconstruction: compressed sensing and co-processing with multi-dimensional decomposition. The new methodology was validated using disordered protein CD79A from B-cell receptor and an SH3 domain from Abp1p in exchange between its free form and bound to a peptide from the protein Ark1p. We show that, while providing much better resolution, the 3D NUS experiments give the similar accuracy and precision of the dynamic parameters to ones obtained using traditional 2D experiments. Furthermore, we show that jackknife resampling of the spectra yields robust estimates of peak intensities errors, eliminating the need for recording duplicate data points.

Acceleration of protein backbone NMR assignment by combinatorial labeling: Application to a small molecule binding study

Selective labeling with stable isotopes has long been recognized as a valuable tool in protein NMR to alleviate signal overlap and sensitivity limitations. In this study, combinatorial 15 N-, 13 Cα -, and 13 C'-selective labeling has been used during the backbone assignment of human cyclophilin D to explore binding of an inhibitor molecule. Using a cell-free expression system, a scheme that involves 15 N, 1-13 C, 2-13 C, fully 15 N/13 C, and unlabeled amino acids was optimized to gain a maximum of assignment information from three samples. This scheme was combined with time-shared triple-resonance NMR experiments, which allows a fast and efficient backbone assignment by giving the unambiguous assignment of unique amino acid pairs in the protein, the identity of ambiguous pairs and information about all 19 non-proline amino acid types. It is therefore well suited for binding studies where de novo assignments of amide 1 H and 15 N resonances need to be obtained, even in cases where sensitivity is the limiting factor.

Accurate determination of rates from non-uniformly sampled relaxation data

The application of non-uniform sampling (NUS) to relaxation experiments traditionally used to characterize the fast internal motion of proteins is quantitatively examined. Experimentally acquired Poisson-gap sampled data reconstructed with iterative soft thresholding are compared to regular sequentially sampled (RSS) data. Using ubiquitin as a model system, it is shown that 25 % sampling is sufficient for the determination of quantitatively accurate relaxation rates. When the sampling density is fixed at 25 %, the accuracy of rates is shown to increase sharply with the total number of sampled points until eventually converging near the inherent reproducibility of the experiment. Perhaps contrary to some expectations, it is found that accurate peak height reconstruction is not required for the determination of accurate rates. Instead, inaccuracies in rates arise from inconsistencies in reconstruction across the relaxation series that primarily manifest as a non-linearity in the recovered peak height. This indicates that the performance of an NUS relaxation experiment cannot be predicted from comparison of peak heights using a single RSS reference spectrum. The generality of these findings was assessed using three alternative reconstruction algorithms, eight different relaxation measurements, and three additional proteins that exhibit varying degrees of spectral complexity. From these data, it is revealed that non-linearity in peak height reconstruction across the relaxation series is strongly correlated with errors in NUS-derived relaxation rates. Importantly, it is shown that this correlation can be exploited to reliably predict the performance of an NUS-relaxation experiment by using three or more RSS reference planes from the relaxation series. The RSS reference time points can also serve to provide estimates of the uncertainty of the sampled intensity, which for a typical relaxation times series incurs no penalty in total acquisition time.

Non-uniform sampling of NMR relaxation data

The use of non-uniform sampling of NMR spectra may give significant reductions in the data acquisition time. For quantitative experiments such as the measurement of spin relaxation rates, non-uniform sampling is however not widely used as inaccuracies in peak intensities may lead to errors in the extracted dynamic parameters. By systematic reducing the coverage of the Nyquist grid of (15)N Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion datasets for four different proteins and performing a full data analysis of the resulting non-uniform sampled datasets, we have compared the performance of the multi-dimensional decomposition and iterative re-weighted least-squares algorithms in reconstructing spectra with accurate peak intensities. As long as a single fully sampled spectrum is included in a series of otherwise non-uniform sampled two-dimensional spectra, multi-dimensional decomposition reconstructs the non-uniform sampled spectra with high accuracy. For two of the four analyzed datasets, a coverage of only 20% results in essentially the same results as the fully sampled data. As exemplified by other data, such a low coverage is in general not enough to produce reliable results. We find that a coverage level not compromising the final results can be estimated by recording a single full two-dimensional spectrum and reducing the spectrum quality in silico.

Backbone Assignment of the MALT1 Paracaspase by Solution NMR

Mucosa-associated lymphoid tissue lymphoma translocation protein 1 (MALT1) is a unique paracaspase protein whose protease activity mediates oncogenic NF-κB signalling in activated B cell-like diffuse large B cell lymphomas (ABC-DLBCLs). ABC-DLBCLs are aggressive lymphomas with high resistance to current chemotherapies. Low survival rate among patients emphasizes the urgent need for alternative treatment options. The characterization of the MALT1 will be an essential tool for developing new target-directed drugs against MALT1 dependent disorders. As the first step in the atomic-level NMR studies of the system, here we report, the (15)N/(13)C/(1)H backbone assignment of the apo form of the MALT1 paracaspase region together with the third immunoglobulin-like (Ig3) domain, 44 kDa, by high resolution NMR. In addition, the non-uniform sampling (NUS) based targeted acquisition procedure is evaluated as a mean of decreasing acquisition and analysis time for larger proteins.

Lorentzian sparsity based spectroscopic reconstruction for fast high-dimensional magnetic resonance spectroscopy

Two-dimensional magnetic resonance spectroscopy (2D MRS) is challenging, even with state-of-art compressive sensing methods, such as L1-sparsity method. In this work, using the prior that the 2D MRS can be regarded as a series of Lorentzian functions, we aim to develop a robust Lorentzian-sparsity based spectroscopy reconstruction method for high-dimensional MRS. The proposed method sparsifies 2D MRS in Lorentzian functions. Instead of thousands of pixel-wise variables, this Lorentzian-sparsity method significantly reduces the number of unknowns to several geometric variables, such as the center, magnitude and shape parameters for each Lorentzian function. The spectroscopy reconstruction is formulated as a nonlinear and nonconvex optimization problem, and the simulated annealing algorithm is developed to solve the problem. The proposed method was compared with inverse FFT method and L1-sparsity method, under various undersampling factors. While FFT and L1 results contained severe artifacts, the Lorentzian-sparsity results provided significantly improved spectroscopy. A new 2D MRS reconstruction method is proposed using the Lorentzian sparsity, with significantly improved MRS reconstruction quality, in comparison with standard inverse FFT method or state-of-art L1-sparsity method.

Applications of high dimensionality experiments to biomolecular NMR

High dimensionality NMR experiments facilitate resonance assignment and precise determination of spectral parameters such as coupling constants. Sparse non-uniform sampling enables acquisition of experiments of high dimensionality with high resolution in acceptable time. In this review we present and compare some significant applications of NMR experiments of dimensionality higher than three in the field of biomolecular studies in solution.

Non-uniform sampling: post-Fourier era of NMR data collection and processing

The invention of multidimensional techniques in the 1970s revolutionized NMR, making it the general tool of structural analysis of molecules and materials. In the most straightforward approach, the signal sampling in the indirect dimensions of a multidimensional experiment is performed in the same manner as in the direct dimension, i.e. with a grid of equally spaced points. This results in lengthy experiments with a resolution often far from optimum. To circumvent this problem, numerous sparse-sampling techniques have been developed in the last three decades, including two traditionally distinct approaches: the radial sampling and non-uniform sampling. This mini review discusses the sparse signal sampling and reconstruction techniques from the point of view of an underdetermined linear algebra problem that arises when a full, equally spaced set of sampled points is replaced with sparse sampling. Additional assumptions that are introduced to solve the problem, as well as the shape of the undersampled Fourier transform operator (visualized as so-called point spread function), are shown to be the main differences between various sparse-sampling methods.

Development of a method for reconstruction of crowded NMR spectra from undersampled time-domain data

NMR is a unique methodology for obtaining information about the conformational dynamics of proteins in heterogeneous biomolecular systems. In various NMR methods, such as transferred cross-saturation, relaxation dispersion, and paramagnetic relaxation enhancement experiments, fast determination of the signal intensity ratios in the NMR spectra with high accuracy is required for analyses of targets with low yields and stabilities. However, conventional methods for the reconstruction of spectra from undersampled time-domain data, such as linear prediction, spectroscopy with integration of frequency and time domain, and analysis of Fourier, and compressed sensing were not effective for the accurate determination of the signal intensity ratios of the crowded two-dimensional spectra of proteins. Here, we developed an NMR spectra reconstruction method, "conservation of experimental data in analysis of Fourier" (Co-ANAFOR), to reconstruct the crowded spectra from the undersampled time-domain data. The number of sampling points required for the transferred cross-saturation experiments between membrane proteins, photosystem I and cytochrome b 6 f, and their ligand, plastocyanin, with Co-ANAFOR was half of that needed for linear prediction, and the peak height reduction ratios of the spectra reconstructed from truncated time-domain data by Co-ANAFOR were more accurate than those reconstructed from non-uniformly sampled data by compressed sensing.

Reduced dimensionality (3,2)D NMR experiments and their automated analysis: implications to high-throughput structural studies on proteins

Protein NMR spectroscopy has expanded dramatically over the last decade into a powerful tool for the study of their structure, dynamics, and interactions. The primary requirement for all such investigations is sequence-specific resonance assignment. The demand now is to obtain this information as rapidly as possible and in all types of protein systems, stable/unstable, soluble/insoluble, small/big, structured/unstructured, and so on. In this context, we introduce here two reduced dimensionality experiments – (3,2)D-hNCOcanH and (3,2)D-hNcoCAnH – which enhance the previously described 2D NMR-based assignment methods quite significantly. Both the experiments can be recorded in just about 2-3 h each and hence would be of immense value for high-throughput structural proteomics and drug discovery research. The applicability of the method has been demonstrated using alpha-helical bovine apo calbindin-D9k P43M mutant (75 aa) protein. Automated assignment of this data using AUTOBA has been presented, which enhances the utility of these experiments. The backbone resonance assignments so derived are utilized to estimate secondary structures and the backbone fold using Web-based algorithms. Taken together, we believe that the method and the protocol proposed here can be used for routine high-throughput structural studies of proteins.

Modified OMP algorithm for exponentially decaying signals

A group of signal reconstruction methods, referred to as compressed sensing (CS), has recently found a variety of applications in numerous branches of science and technology. However, the condition of the applicability of standard CS algorithms (e.g., orthogonal matching pursuit, OMP), i.e., the existence of the strictly sparse representation of a signal, is rarely met. Thus, dedicated algorithms for solving particular problems have to be developed. In this paper, we introduce a modification of OMP motivated by nuclear magnetic resonance (NMR) application of CS. The algorithm is based on the fact that the NMR spectrum consists of Lorentzian peaks and matches a single Lorentzian peak in each of its iterations. Thus, we propose the name Lorentzian peak matching pursuit (LPMP). We also consider certain modification of the algorithm by introducing the allowed positions of the Lorentzian peaks' centers. Our results show that the LPMP algorithm outperforms other CS algorithms when applied to exponentially decaying signals.

Nonuniform sampling and non-Fourier signal processing methods in multidimensional NMR

Beginning with the introduction of Fourier Transform NMR by Ernst and Anderson in 1966, time domain measurement of the impulse response (the free induction decay, FID) consisted of sampling the signal at a series of discrete intervals. For compatibility with the discrete Fourier transform (DFT), the intervals are kept uniform, and the Nyquist theorem dictates the largest value of the interval sufficient to avoid aliasing. With the proposal by Jeener of parametric sampling along an indirect time dimension, extension to multidimensional experiments employed the same sampling techniques used in one dimension, similarly subject to the Nyquist condition and suitable for processing via the discrete Fourier transform. The challenges of obtaining high-resolution spectral estimates from short data records using the DFT were already well understood, however. Despite techniques such as linear prediction extrapolation, the achievable resolution in the indirect dimensions is limited by practical constraints on measuring time. The advent of non-Fourier methods of spectrum analysis capable of processing nonuniformly sampled data has led to an explosion in the development of novel sampling strategies that avoid the limits on resolution and measurement time imposed by uniform sampling. The first part of this review discusses the many approaches to data sampling in multidimensional NMR, the second part highlights commonly used methods for signal processing of such data, and the review concludes with a discussion of other approaches to speeding up data acquisition in NMR.

Protein interacting with C-kinase 1 (PICK1) binding promiscuity relies on unconventional PSD-95/discs-large/ZO-1 homology (PDZ) binding modes for nonclass II PDZ ligands

PDZ domain proteins control multiple cellular functions by governing assembly of protein complexes. It remains unknown why individual PDZ domains can bind the extreme C terminus of very diverse binding partners and maintain selectivity. By employing NMR spectroscopy, together with molecular modeling, mutational analysis, and fluorescent polarization binding experiments, we identify here three structural mechanisms explaining why the PDZ domain of PICK1 selectively binds >30 receptors, transporters, and kinases. Class II ligands, including the dopamine transporter, adopt a canonical binding mode with promiscuity obtained via differential packing in the binding groove. Class I ligands, such as protein kinase Cα, depend on residues upstream from the canonical binding sequence that are likely to interact with flexible loop residues of the PDZ domain. Finally, we obtain evidence that the unconventional ligand ASIC1a has a dual binding mode involving a canonical insertion and a noncanonical internal insertion with the two C-terminal residues forming interactions outside the groove. Together with an evolutionary analysis, the data show how unconventional binding modes might evolve for a protein recognition domain to expand the repertoire of functionally important interactions.

Tyrosine phosphorylation within the intrinsically disordered cytosolic domains of the B-cell receptor: an NMR-based structural analysis

Intrinsically disordered proteins are found extensively in cell signaling pathways where they often are targets of posttranslational modifications e.g. phosphorylation. Such modifications can sometimes induce or disrupt secondary structure elements present in the modified protein. CD79a and CD79b are membrane-spanning, signal-transducing components of the B-cell receptor. The cytosolic domains of these proteins are intrinsically disordered and each has an immunoreceptor tyrosine-based activation motif (ITAM). When an antigen binds to the receptor, conserved tyrosines located in the ITAMs are phosphorylated which initiate further downstream signaling. Here we use NMR spectroscopy to examine the secondary structure propensity of the cytosolic domains of CD79a and CD79b in vitro before and after phosphorylation. The phosphorylation patterns are identified through analysis of changes of backbone chemical shifts found for the affected tyrosines and neighboring residues. The number of the phosphorylated sites is confirmed by mass spectrometry. The secondary structure propensities are calculated using the method of intrinsic referencing, where the reference random coil chemical shifts are measured for the same protein under denaturing conditions. Our analysis revealed that CD79a and CD79b both have an overall propensity for α-helical structure that is greatest in the C-terminal region of the ITAM. Phosphorylation of CD79a caused a decrease in helical propensity in the C-terminal ITAM region. For CD79b, the opposite was observed and phosphorylation resulted in an increase of helical propensity in the C-terminal part.