Improved analysis of NMR chemical shift perturbations through an error estimation method

In solution NMR, chemical shift perturbation (CSP) experiments are widely employed to study intermolecular interactions. However, excluding the nonsignificant peak shift is difficult because little is known about errors in CSP. Here, to address this issue, we introduce a method for estimating errors in CSP based on the noise level. First, we developed a technique that involves line shape fitting to estimate errors in peak position via Monte Carlo simulations. Second, this technique was applied to estimate errors in CSP. In intermolecular interaction analysis of VAP-A with SNX2, error estimation of CSP enabled the evaluation of small but significant changes in peak position and yielded detailed insights that are unattainable with conventional CSP analysis. Third, this technique was successfully applied to estimate errors in residual dipolar couplings. In conclusion, our error estimation method improves CSP analysis by excluding the nonsignificant peak shift.

Lipid- and substrate-induced conformational and dynamic changes in a glycosyltransferase involved in E. coli LPS synthesis revealed by 19F and 31P NMR

WaaG is a glycosyltransferase (GT) involved in the synthesis of the bacterial cell wall, and in Escherichia coli it catalyzes the transfer of a glucose moiety from the donor substrate UDP-glucose onto the nascent lipopolysaccharide (LPS) molecule which when completed constitutes the major component of the bacterium's outermost defenses. Similar to other GTs of the GT-B fold, having two Rossman-like domains connected by a short linker, WaaG is believed to undergo complex inter-domain motions as part of its function to accommodate the nascent LPS and UDP-glucose in the catalytic site located in the cleft between the two domains. As the nascent LPS is bulky and membrane-bound, WaaG is a peripheral membrane protein, adding to the complexity of studying the enzyme in a biologically relevant environment. Using specific 5-fluoro-Trp labelling of native and inserted tryptophans and 19F NMR we herein studied the dynamic interactions of WaaG with lipids using bicelles, and with the donor substrate. Line-shape changes when bicelles are added to WaaG show that the dynamic behavior is altered when binding to the model membrane, while a chemical shift change indicates an altered environment around a tryptophan located in the C-terminal domain of WaaG upon interaction with UDP-glucose or UDP. A lipid-bound paramagnetic probe was used to confirm that the membrane interaction is mediated by a loop region located in the N-terminal domain. Furthermore, the hydrolysis of the donor substrate by WaaG was quantified by 31P NMR.

The protease associated (PA) domain in ScpA from Streptococcus pyogenes plays a role in substrate recruitment

Annually, over 18 million disease cases and half a million deaths worldwide are estimated to be caused by Group A Streptococcus. ScpA (or C5a peptidase) is a well characterised member of the cell enveleope protease family, which possess a S8 subtilisin-like catalytic domain and a shared multi-domain architecture. ScpA cleaves complement factors C5a and C3a, impairing the function of these critical anaphylatoxins and disrupts complement-mediated innate immunity. Although the high resolution structure of ScpA is known, the details of how it recognises its substrate are only just emerging. Previous studies have identified a distant exosite on the 2nd fibronectin domain that plays an important role in recruitment via an interaction with the substrate core. Here, using a combination of solution NMR spectroscopy, mutagenesis with functional assays and computational approaches we identify a second exosite within the protease-associated (PA) domain. We propose a model in which the PA domain assists optimal delivery of the substrate's C terminus to the active site for cleavage.

1H, 15N, and 13C chemical shift backbone resonance NMR assignment of the accumulation-associated protein (Aap) lectin domain from Staphylococcus epidermidis

Staphylococcus epidermidis is the leading causative agent for hospital-acquired infections, especially device-related infections, due to its ability to form biofilms. The accumulation-associated protein (Aap) of S. epidermidis is primarily responsible for biofilm formation and consists of two domains, A and B. It was found that the A domain is responsible for the attachment to the abiotic/biotic surface, whereas the B domain is responsible for the accumulation of bacteria during biofilm formation. One of the parts of the A domain is the Aap lectin, which is a carbohydrate-binding domain having 222 amino acids in its structure. Here we report the near complete backbone chemical shift assignments for the lectin domain, as well as its predicted secondary structure. This data will provide a platform for future NMR studies to explore the role of lectin in biofilm formation.

Solvent accessibility of a GPCR transmembrane domain probed by in-membrane chemical modification (IMCM)

G protein-coupled receptors (GPCRs) transmit signals from drugs across cell membranes, leading to associated physiological effects. To study the structural basis of the transmembrane signaling, in-membrane chemical modification (IMCM) has previously been introduced for ¹⁹ F-labeling of GPCRs expressed in Spodoptera frugiperda (Sf9) insect cells. Here, IMCM is used with the A_2A adenosine receptor (A_2A AR) expressed in Pichia pastoris; ¹⁹ F-NMR revealed nearly complete solvent protection of the A_2A AR transmembrane domain in the membrane and in 2,2-didecylpropane-1,3-bis-β-D-maltopyranoside (LMNG)/cholesteryl hemisuccinate (CHS) micelles, and extensive solvent accessibility for A_2A AR in n-dodecyl β-D-maltoside (DDM)/CHS micelles. No Cys residue dominated non-specific labeling with 2,2,2-trifluoroethanethiol. These observations yield an improved protocol for IMCM ¹⁹ F-labeling of GPCRs and new insights into variable solvent accessibility for function-related characterization of GPCRs.

Structures and conformational dynamics of DNA minidumbbells in pyrimidine-rich repeats associated with neurodegenerative diseases

Expansions of short tandem repeats (STRs) are associated with approximately 50 human neurodegenerative diseases. These pathogenic STRs are prone to form non-B DNA structure, which has been considered as one of the causative factors for repeat expansions. Minidumbbell (MDB) is a relatively new type of non-B DNA structure formed by pyrimidine-rich STRs. An MDB is composed of two tetraloops or pentaloops, exhibiting a highly compact conformation with extensive loop-loop interactions. The MDB structures have been found to form in CCTG tetranucleotide repeats associated with myotonic dystrophy type 2, ATTCT pentanucleotide repeats associated with spinocerebellar ataxia type 10, and the recently discovered ATTTT/ATTTC repeats associated with spinocerebellar ataxia type 37 and familial adult myoclonic epilepsy. In this review, we first introduce the structures and conformational dynamics of MDBs with a focus on the high-resolution structural information determined by nuclear magnetic resonance spectroscopy. Then we discuss the effects of sequence context, chemical environment, and nucleobase modification on the structure and thermostability of MDBs. Finally, we provide perspectives on further explorations of sequence criteria and biological functions of MDBs.

Backbone and side chain resonance assignment of the intrinsically disordered human DBNDD1 protein

The dysbindin domain-containing protein 1 (DBNDD1) is a conserved protein among higher eukaryotes whose structure and function are poorly investigated so far. Here, we present the backbone and side chain nuclear magnetic resonance assignments for the human DBNDD1 protein. Our chemical-shift based secondary structure analysis reveals the human DBNDD1 as an intrinsically disordered protein.

Adsorption mechanisms of inositol hexakisphosphate in the presence of phosphate at the amorphous aluminum oxyhydroxide-water interface

Myo-inositol hexakisphosphate (myo-IHP) is one of the most common soil organic phosphorus (P) species in soil. Its retention in soil is often competed by phosphate, making bioavailability of P species difficult. In this study, the adsorption mechanism of myo-IHP at the amorphous aluminum (oxyhydr)oxide (AAH)-water interface was investigated at pH 6.5 in the presence of phosphate using batch adsorption experiments and solution ³¹P NMR spectroscopy. The ratio of [myo-IHP]_i/[phosphate]_i (Ri) was kept 0.33-3 while ligand addition was varied. In the absence of phosphate, myo-IHP forms inner-sphere surface complexes in AAH via P1,3, P2, P4,6, and P5 functional group coordination. When two ligands were simultaneously added, fewer P functional groups of myo-IHP coordinated to AAH and the surface complexes were altered with the coordination of mainly P1,3 and P2 functional groups. When phosphate was pre-adsorbed, myo-IHP adsorption decreased by 8.0-44% compared to the respective simultaneous addition system. P2 or P5 functional group was predominantly coordinated to the AAH surfaces at Ri = 0.33. Myo-IHP pre-adsorption resulted in an increase in the final myo-IHP adsorption compared to that in the simultaneous addition system under the respective Ri values (0.33-3). In this system, P1,3, P2, P4,6, and P5 functional groups were coordinated to form inner-sphere surface complexes regardless of Ri. The study revealed that the functional group specific adsorption mechanism of myo-IHP at the AAH-water interface was affected by addition sequence and Ri of two ligands. The competitive adsorption between organic P and phosphate plays an important role in the fate of P in soils.

Backbone NMR assignment of the nucleotide binding domain of the Bacillus subtilis ABC multidrug transporter BmrA in the post-hydrolysis state

ATP binding cassette (ABC) proteins are present in all phyla of life and form one of the largest protein families. The Bacillus subtilis ABC transporter BmrA is a functional homodimer that can extrude many different harmful compounds out of the cell. Each BmrA monomer is composed of a transmembrane domain (TMD) and a nucleotide binding domain (NBD). While the TMDs of ABC transporters are sequentially diverse, the highly conserved NBDs harbor distinctive conserved motifs that enable nucleotide binding and hydrolysis, interdomain communication and that mark a protein as a member of the ABC superfamily. In the catalytic cycle of an ABC transporter, the NBDs function as the molecular motor that fuels substrate translocation across the membrane via the TMDs and are thus pivotal for the entire transport process. For a better understanding of the structural and dynamic consequences of nucleotide interactions within the NBD at atomic resolution, we determined the ¹H, ¹³C and ¹⁵N backbone chemical shift assignments of the 259 amino acid wildtype BmrA-NBD in its post-hydrolytic, ADP-bound state.

Backbone chemical shift spectral assignments of SARS coronavirus-2 non-structural protein nsp9

As part of an International consortium aiming at the characterization by NMR of the proteins of the SARS-CoV-2 virus, we have obtained the virtually complete assignment of the backbone atoms of the non-structural protein nsp9. This small (12 kDa) protein is encoded by ORF1a, binds to RNA and seems to be essential for viral RNA synthesis. The crystal structures of the SARS-CoV-2 protein and other homologues suggest that the protein is dimeric as also confirmed by analytical ultracentrifugation and dynamic light scattering. Our data constitute the prerequisite for further NMR-based characterization, and provide the starting point for the identification of small molecule lead compounds that could interfere with RNA binding and prevent viral replication.

1H, 13C, and 15N NMR chemical shift assignment of the complex formed by the first EPEC EspF repeat and N-WASP GTPase binding domain

LEE-encoded effector EspF (EspF) is an effector protein part of enteropathogenic Escherichia coli's (EPEC's) arsenal for intestinal infection. This intrinsically disordered protein contains three highly conserved repeats which together compose over half of the protein's complete amino acid sequence. EPEC uses EspF to hijack host proteins in order to promote infection. In the attack EspF is translocated, together with other effector proteins, to host cell via type III secretion system. Inside host EspF stimulates actin polymerization by interacting with Neural Wiskott-Aldrich syndrome protein (N-WASP), a regulator in actin polymerization machinery. It is presumed that EspF acts by disrupting the autoinhibitory state of N-WASP GTPase binding domain. In this NMR spectroscopy study, we report the ¹H, ¹³C, and ¹⁵N resonance assignments for the complex formed by the first 47-residue repeat of EspF and N-WASP GTPase binding domain. These near-complete resonance assignments provide the basis for further studies which aim to characterize structure, interactions, and dynamics between these two proteins in solution.

Backbone chemical shift assignments for the SARS-CoV-2 non-structural protein Nsp9: intermediate (ms - μs) dynamics in the C-terminal helix at the dimer interface

The Betacoronavirus SARS-CoV-2 non-structural protein Nsp9 is a 113-residue protein that is essential for viral replication, and consequently, a potential target for the development of therapeutics against COVID19 infections. To capture insights into the dynamics of the protein's backbone in solution and accelerate the identification and mapping of ligand-binding surfaces through chemical shift perturbation studies, the backbone ¹H, ¹³C, and ¹⁵N NMR chemical shifts for Nsp9 have been extensively assigned. These assignments were assisted by the preparation of an ~ 70% deuterated sample and residue-specific, ¹⁵N-labelled samples (V, L, M, F, and K). A major feature of the assignments was the "missing" amide resonances for N96-L106 in the ¹H-¹⁵N HSQC spectrum, a region that comprises almost the complete C-terminal α-helix that forms a major part of the homodimer interface in the crystal structure of SARS-CoV-2 Nsp9, suggesting this region either undergoes intermediate motion in the ms to μs timescale and/or is heterogenous. These "missing" amide resonances do not unambiguously appear in the ¹H-¹⁵N HSQC spectrum of SARS-CoV-2 Nsp9 collected at a concentration of 0.0007 mM. At this concentration, at the detection limit, native mass spectrometry indicates the protein is exclusively in the monomeric state, suggesting the intermediate motion in the C-terminal of Nsp9 may be due to intramolecular dynamics. Perhaps this intermediate ms to μs timescale dynamics is the physical basis for a previously suggested "fluidity" of the C-terminal helix that may be responsible for homophilic (Nsp9-Nsp9) and postulated heterophilic (Nsp9-Unknown) protein-protein interactions.

Exploring interactions between lipids and amyloid-forming proteins: A review on applying fluorescence and NMR techniques

A hallmark of Alzheimer's, Parkinson's, and other amyloid diseases is the assembly of amyloid proteins into amyloid aggregates or fibrils. In many cases, the formation and cytotoxicity of amyloid assemblies are associated with their interaction with cell membranes. Despite studied for many years, the characterization of the interaction is challenged for reasons on the multiple aggregation states of amyloid-forming proteins, transient and weak interactions in the complex system. Although several strategies such as computation biology, spectroscopy, and imaging methods have been performed, there is an urgent need to detail the molecular mechanism in different time scales and high resolutions. This review highlighted the recent applications of fluorescence, solution and solid-state NMR in exploring the interactions between amyloid protein and membranes attributing to their advantages of high sensitivity and atomic resolution.

Efficient production of a functional G protein-coupled receptor in E. coli for structural studies

G protein-coupled receptors (GPCRs) are transmembrane signal transducers which regulate many key physiological process. Since their discovery, their analysis has been limited by difficulties in obtaining sufficient amounts of the receptors in high-quality, functional form from heterologous expression hosts. Albeit highly attractive because of its simplicity and the ease of isotope labeling for NMR studies, heterologous expression of functional GPCRs in E. coli has proven particularly challenging due to the absence of the more evolved protein expression and folding machinery of higher eukaryotic hosts. Here we first give an overview on the previous strategies for GPCR E. coli expression and then describe the development of an optimized robust protocol for the E. coli expression and purification of two mutants of the turkey β₁-adrenergic receptor (β₁AR) uniformly or selectively labeled in ¹⁵N or ²H,¹⁵N. These mutants had been previously optimized for thermal stability using insect cell expression and used successfully in crystallographic and NMR studies. The same sequences were then used for E. coli expression. Optimization of E. coli expression was achieved by a quantitative analysis of losses of receptor material at each step of the solubilization and purification procedure. Final yields are 0.2-0.3 mg receptor per liter culture. Whereas both expressed mutants are well folded and competent for orthosteric ligand binding, the less stable YY-β₁AR mutant also comprises the two native tyrosines Y^5.58 and Y^7.53, which enable G protein binding. High-quality ¹H-¹⁵N TROSY spectra were obtained for E. coli-expressed YY-β₁AR in three different functional states (antagonist, agonist, and agonist + G protein-mimicking nanobody-bound), which are identical to spectra obtained of the same forms of the receptor expressed in insect cells. NdeI and AgeI restriction sites introduced into the expression plasmid allow for the easy replacement of the receptor gene by other GPCR genes of interest, and the provided quantitative workflow analysis may guide the respective adaptation of the purification protocol.

Expression, purification and characterization of TMCO1 for structural studies

Transmembrane and coiled-coil domains 1 (TMCO1) has a highly conserved amino acid sequence among species, indicating a critical role of TMCO1 in cell physiology. The deficiency of TMCO1 in humans is associated with cerebrofaciothoracic dysplasia (CFTD), glaucoma, osteogenesis and the occurrence of cancer. TMCO1 was recently identified as an endoplasmic reticulum (ER) Ca²⁺ load-activated Ca²⁺ (CLAC) release channel, which prevents ER Ca²⁺ overload and maintains calcium homeostasis in the ER. However, the structural basis of the molecular function of TMCO1 channel remains elusive. To determine the structure of TMCO1, we screened the expression of TMCO1 in Escherichia coli and insect cell expression systems. TMCO1 from Dictyostelium discoideum (DdTMCO1) was successfully expressed in Escherichia coli with a high yield. The pure recombinant protein was obtained by affinity chromatography and size exclusion chromatography. The solution NMR of DdTMCO1 in DPC micelles showed three α-helical transmembrane regions.

1H, 13C, 15N resonance assignments and secondary structure of yeast oligosaccharyltransferase subunit Ost4 and its functionally important mutant Ost4V23D

Asparagine-linked glycosylation is an essential and highly conserved protein modification reaction that occurs in the endoplasmic reticulum of cells during protein synthesis at the ribosome. In the central reaction, a pre-assembled high-mannose sugar is transferred from a lipid-linked donor substrate to the side-chain of an asparagine residue in an -N-X-T/S- sequence (where X is any residue except proline). This reaction is carried by a membrane-bound multi-subunit enzyme complex, oligosaccharyltransferase (OST). In humans, genetic defects in OST lead to a group of rare metabolic diseases collectively known as Congenital Disorders of Glycosylation. Certain mutations are lethal for all organisms. In yeast, the OST is composed of nine non-identical protein subunits. The functional enzyme complex contains eight subunits with either Ost3 or Ost6 at any given time. Ost4, an unusually small protein, plays a very important role in the stabilization of the OST complex. It bridges the catalytic subunit Stt3 with Ost3 (or Ost6) in the Stt3-Ost4-Ost3 (or Ost6) sub-complex. Mutation of any residue from M18-I24 in the trans-membrane helix of yeast Ost4 negatively impacts N-linked glycosylation and the growth of yeast. Indeed, mutation of valine23 to an aspartate impairs OST function in vivo resulting in a lethal phenotype in yeast. To understand the structural mechanism of Ost4 in the stabilization of the enzyme complex, we have initiated a detailed investigation of Ost4 and its functionally important mutant, Ost4V23D. Here, we report the backbone ¹H, ¹³C, and ¹⁵N resonance assignments for Ost4 and Ost4V23D in dodecylphosphocholine micelles.

Principal component analysis of data from NMR titration experiment of uniformly 15N labeled amyloid beta (1-42) peptide with osmolytes and phenolic compounds

A simple NMR method to analyze the data obtained by NMR titration experiment of amyloid formation inhibitors against uniformly 15N-labeled amyloid-β 1-42 peptide (Aβ(1-42)) was described. By using solution nuclear magnetic resonance (NMR) measurement, the simplest method for monitoring the effects of Aβ fibrilization inhibitors is the NMR chemical shift perturbation (CSP) experiment using 15N-labeled Aβ(1-42). However, the flexible and dynamic nature of Aβ(1-42) monomer may hamper the interpretation of CSP data. Here we introduced principal component analysis (PCA) for visualizing and analyzing NMR data of Aβ(1-42) in the presence of amyloid inhibitors including high concentration osmolytes. We measured 1H-15N 2D spectra of Aβ(1-42) at various temperatures as well as of Aβ(1-42) with several inhibitors, and subjected all the data to PCA (PCA-HSQC). The PCA diagram succeeded in differentiating the various amyloid inhibitors, including epigallocatechin gallate (EGCg), rosmarinic acid (RA) and curcumin (CUR) from high concentration osmolytes. We hypothesized that the CSPs reflected the conformational equilibrium of intrinsically disordered Aβ(1-42) induced by weak inhibitor binding rather than the specific molecular interactions.

Analyzing multi-step ligand binding reactions for oligomeric proteins by NMR: Theoretical and computational considerations

Solution NMR spectroscopy is widely used to investigate the thermodynamics and kinetics of the binding of ligands to their biological receptors, as it provides detailed, atomistic information, potentially leading to microscopic affinities for each binding event, and, to the development of allosteric pathways describing how the binding at one site affects distal sites in the molecule. Importantly, weak interactions that are often invisible to other biophysical methods can also be probed. Methodological advancements in NMR have enabled the investigation of high molecular weight, homo-oligomeric complexes that bind multiple ligand molecules, with increasing numbers of studies of the structural dynamics and binding properties of these systems. It therefore becomes of interest to consider how binding and kinetics parameters can be extracted from experiments on these more complicated molecules. Here we present the theoretical framework for analyzing binding reactions of homo-oligomeric complexes by NMR, taking into account all of the chemical species in solution and their corresponding NMR observables. A number of simulations are presented to illustrate the utility of the derived expressions.

Optimizing the α1B-adrenergic receptor for solution NMR studies

Sample preparation for NMR studies of G protein-coupled receptors faces special requirements: Proteins need to be stable for prolonged measurements at elevated temperatures, they should ideally be uniformly labeled with the stable isotopes ¹³C, ¹⁵N, and all carbon-bound protons should be replaced by deuterons. In addition, certain NMR experiments require protonated methyl groups in the presence of a perdeuterated background. All these requirements are most easily satisfied when using Escherichia coli as the expression host. Here we describe a workflow, starting from a temperature-stabilized mutant of the α_1B-adrenergic receptor, obtained using the CHESS methodology, into an even more stable species, in which flexible parts from termini were removed and the intracellular loop 3 (ICL3) was stabilized against proteolytic cleavage. The yield after purification corresponds to 1-2 mg/L of D₂O culture. The final purification step is ligand-affinity chromatography to ensure that only well-folded ligand-binding protein is isolated. Proper selection of detergent has a remarkable influence on the quality of NMR spectra. All optimization steps of sequence and detergent are monitored on a small scale by monitoring the melting temperature and long-term thermal stability to allow for screening of many conditions. The stabilized mutant of the α_1B-adrenergic receptor was additionally incorporated in nanodiscs, but displayed slightly inferior spectra compared to a sample in detergent micelles. Finally, both [¹⁵N,¹H]- as well as [¹³C,¹H]-HSQC spectra are shown highlighting the high quality of the final NMR sample. Importantly, the quality of [¹³C,¹H]-HSQC spectra indicates that the so prepared receptor could be used for studying side-chain dynamics.

Complete sequential assignment and secondary structure prediction of the cannulae forming protein CanA from the hyperthermophilic archaeon Pyrodictium abyssi

CanA from Pyrodictium abyssi forms a heat-resistant organic hollow-fiber network together with CanB and CanC. An N-terminally truncated construct of CanA (K1-CanA) gave NMR spectra of good quality that could be assigned by three-dimensional NMR methods on 15N and 13C-15N enriched protein. We assigned the chemical shifts of 96% of all backbone 1HN atoms, 98% of all backbone 15N atoms, 100% of all 13Cα atoms, 100% of all 1Hα atoms, 90% of all 13C' atoms, and 100% of the 13Cβ atoms. Two short helices and 10 β-strands are estimated from an analysis of the chemical shifts leading to a secondary structure content of K1-CanA of 6% helices, 44% β-pleated sheets, and 50% coils.

Recording In-Cell NMR-Spectra in Living Mammalian Cells

At the foundation of many cellular processes as well as a large number of diseases is the (mis)folding of important intrinsically disordered proteins (IDPs). Despite tremendous scientific efforts, the factors driving their structural changes within the cellular context remain poorly understood. In-cell NMR spectroscopy enables investigation of IDPs directly in the living eukaryotic cell enabling investigation of its intermolecular interactions and ensuing modifications at an unprecedented atomic resolution. In the following protocol, we describe how to prepare in-cell NMR samples of IDPs within eukaryotic cells and how to measure these in-cell NMR samples of an IDP in its natural environment, the living mammalian cell. Furthermore, we outline a procedure to assess the intracellular recombinant protein concentration of the studied IDP based on in-cell NMR methods. We use α-synuclein as a model protein, but the presented approach is highly modular and therefore should be easily adapted and altered to the desired needs for the studies of different IDPs.

The evolution of solution state NMR pulse sequences through the 'eyes' of triple-resonance spectroscopy

Careful pulse sequence design and optimization is critical to the success of a given NMR experiment. Over the past several decades the level of sophistication of NMR pulse sequences has increased tremendously, leading to large spectral sensitivity and resolution improvements, to data sets with far fewer artifacts, and to much more rapid acquisition times, opening up a wide range of applications. Here I briefly highlight how pulse sequence 'engineering' has evolved, focusing on liquid state NMR, and, in particular, on the HNCA-class of triple-resonance experiment. In many respects, the evolution of triple-resonance NMR mirrors the evolution of solution state NMR experiments in general, with 'tricks' that first appeared in triple-resonance pulse sequences or that were motivated by them now incorporated into a broad range of experiments.

Structure determination using solution NMR: Is it worth the effort?

It has been almost 40 years since solution NMR joined X-ray crystallography as a technique for determining high-resolution structures of proteins. Since then NMR derived structure has contributed in fundamental ways to our understanding of the function of biomolecules. With the already existing mature field of X-ray crystallography and the emergence of cryo-EM as techniques to tackle high-resolution structures of large protein complexes, the role of NMR in structure determination has been questioned. However, NMR has the unique ability to recapitulate the dynamic motion of proteins in their structures, while size limitations of the biomolecular systems that can be routinely studied still present challenges. The field has continually developed methodology and instrumentation since its introduction, pushing its frontiers and redefining its limits. Here we present a brief overview of NMR-based structure determination over the past 40 years. We outline the current state of the field and look ahead to the challenges that still need to be addressed to realize the future potential of NMR as a structural technique.

Integrated BioNMR - "getting by with a little help from my friends"

Single types of methodologies are insufficient to adequately describe complex biological structures. As a result, integrated approaches that combine complementary data are being developed. Here, I describe the benefits of integrating solution and magic angle spinning BioNMR approaches to characterize structure and dynamics of protein assemblies.

Solution NMR Spectroscopy for the Determination of Structures of Membrane Proteins in a Lipid Environment

NMR spectroscopy has harnessed the recent technical advances to emerge as a competitive, elegant, and eminently viable technique for determining the solution structures of membrane proteins at the level of atomic resolution. Once a good level of cell-based or cell-free expression and purification of a suitably sized membrane protein has been achieved, then NMR offers a combination of several versatile strategies, for example choice of appropriate deuterated or nondeuterated detergents, temperature, and ionic strength; isotope labeling with ²H, ¹³C, ¹⁵N, with or without protonation of Ile (δ1), Leu, and Val methyl protons; combinatorial labeling or unlabeling of specific amino acids; TROSY based-, nonuniform sampling (NUS) based-, and other NMR experiments; measurement of residual dipolar couplings using stretched polyacrylamide gels or DNA nanotubes; spin labeling and paramagnetic relaxation enhancements (PRE). Strategic combinations of these advancements together with availability of highly sensitive cryogenically cooled-probes equipped high-field NMR spectrometers (up to 1 GHz ¹H frequency) have allowed the perseverant investigator to successfully overcome several of the conventional pitfalls associated with the NMR technique and membrane proteins, viz., low sensitivity, poor sample stability, spectral crowding, and a limited number of NOEs and other constraints for structure calculations. This has resulted in an unprecedented growth in the number of successfully determined NMR structures of large and complex membrane proteins over the last two decades, and this technique now holds great promise for the structure determination of an ever larger body of membrane proteins.

Succinyl-β-cyclodextrin: Influence of the substitution degree on albendazole inclusion complexes probed by NMR

Succinyl-β-CD derivatives were obtained by green synthesis with degrees of substitution (DS) 1.3 and 2.9. The spray-drying technique was used to obtain albendazole (ABZ):succinyl-β-CD inclusion complexes. Phase solubility diagrams indicated that both succinyl-β-CD derivatives formed 1:1 molar ratio ABZ complexes, but the complex with DS 2.9 has a lower formation constant. The presence of stable inclusion complexes in aqueous solution was confirmed by NMR. For both complexes the aromatic moiety is encapsulated into the host cavity. In the solid-state, ¹³C and ¹⁵N NMR spectral differences between ABZ and ABZ included in spray-dried systems showed that strong structural changes occurred in the systems. At least two different ABZ amorphous species were identified based on DS. ABZ species were stable over more than six months based on spectral data. Finally, the influence of DS in the number and type of the inclusion complexes was elucidated.

Practical aspects of high-pressure NMR spectroscopy and its applications in protein biophysics and structural biology

Pressure and temperature are the two fundamental variables of thermodynamics. Temperature and chemical perturbation are central experimental tools for the exploration of macromolecular structure and dynamics. Though it has long been recognized that hydrostatic pressure offers a complementary and often unique view of macromolecular structure, stability and dynamics, it has not been employed nearly as much. For solution NMR applications the limited use of high-pressure is undoubtedly traced to difficulties of employing pressure in the context of modern multinuclear and multidimensional NMR. Limitations in pressure tolerant NMR sample cells have been overcome and enable detailed studies of macromolecular energy landscapes, hydration, dynamics and function. Here we review the practical considerations for studies of biological macromolecules at elevated pressure, with a particular emphasis on applications in protein biophysics and structural biology.

On the use of Pichia pastoris for isotopic labeling of human GPCRs for NMR studies

NMR studies of human integral membrane proteins provide unique opportunities to probe structure and dynamics at specific locations and on multiple timescales, often with significant implications for disease mechanism and drug development. Since membrane proteins such as G protein-coupled receptors (GPCRs) are highly dynamic and regulated by ligands or other perturbations, NMR methods are potentially well suited to answer basic functional questions (such as addressing the biophysical basis of ligand efficacy) as well as guiding applications (such as novel ligand design). However, such studies on eukaryotic membrane proteins have often been limited by the inability to incorporate optimal isotopic labels for NMR methods developed for large protein/lipid complexes, including methyl TROSY. We review the different expression systems for production of isotopically labeled membrane proteins and highlight the use of the yeast Pichia pastoris to achieve perdeuteration and ¹³C methyl probe incorporation within isoleucine sidechains. We further illustrate the use of this method for labeling of several biomedically significant GPCRs.

Residue-specific mobility changes in soluble oligomers of the prion protein define regions involved in aggregation

Prion (PrP) diseases are neurodegenerative diseases characterized by the formation of β-sheet rich, insoluble and protease resistant protein deposits (called PrP^Sc) that occur throughout the brain. Formation of synthetic or in vitro PrP^Sc can occur through on-pathway toxic oligomers. Similarly, toxic and infectious oligomers identified in cell and animal models of prion disease indicate that soluble oligomers are likely intermediates in the formation of insoluble PrP^Sc. Despite the critical role of prion oligomers in disease progression, little is known about their structure. In order, to obtain structural insight into prion oligomers, we generated oligomers by shaking-induced conversion of recombinant, monomeric prion protein PrP^c (spanning residues 90-231). We then obtained two-dimensional solution NMR spectra of the PrP^c monomer, a 40% converted oligomer, and a 94% converted oligomer. Heteronuclear single-quantum correlation (¹H-¹⁵N) studies revealed that, in comparison to monomeric PrP^c, the oligomer has intense amide peak signals in the N-terminal (residues 90-114) and C-terminal regions (residues 226-231). Furthermore, a core region with decreased mobility is revealed from residues ~127 to 225. Within this core oligomer region with decreased mobility, there is a pocket of increased amide peak signal corresponding to the middle of α-helix 2 and the loop between α-helices 2 and 3 in the PrP^c monomer structure. Using high-resolution solution-state NMR, this work reveals detailed and divergent residue-specific changes in soluble oligomeric models of PrP.

Multiple frequency saturation pulses reduce CEST acquisition time for quantifying conformational exchange in biomolecules

Exchange between conformational states is required for biomolecular catalysis, allostery, and folding. A variety of NMR experiments have been developed to quantify motional regimes ranging from nanoseconds to seconds. In this work, we describe an approach to speed up the acquisition of chemical exchange saturation transfer (CEST) experiments that are commonly used to probe millisecond to second conformational exchange in proteins and nucleic acids. The standard approach is to obtain CEST datasets through the acquisition of a series of 2D correlation spectra where each experiment utilizes a single saturation frequency to ¹H, ¹⁵N or ¹³C. These pseudo 3D datasets are time consuming to collect and are further lengthened by reduced signal to noise stemming from the long saturation pulse. In this article, we show how usage of a multiple frequency saturation pulse (i.e., MF-CEST) changes the nature of data collection from series to parallel, and thus decreases the total acquisition time by an integer factor corresponding to the number of frequencies in the pulse. We demonstrate the applicability of MF-CEST on a Src homology 2 (SH2) domain from phospholipase Cγ and the secondary active transport protein EmrE as model systems by collecting ¹³C methyl and ¹⁵N backbone datasets. MF-CEST can also be extended to additional sites within proteins and nucleic acids. The only notable drawback of MF-CEST as applied to backbone ¹⁵N experiments occurs when a large chemical shift difference between the major and minor populations is present (typically greater than ~ 8 ppm). In these cases, ambiguity may arise between the chemical shift of the minor population and the multiple frequency saturation pulse. Nevertheless, this drawback does not occur for methyl group MF-CEST experiments or in cases where somewhat smaller chemical shift differences occur are present.

Structural elements of stromal interaction molecule function

Stromal interaction molecule (STIM)-1 and -2 are multi-domain, single-pass transmembrane proteins involved in sensing changes in compartmentalized calcium (Ca²⁺) levels and transducing this cellular signal to Orai1 channel proteins. Our understanding of the molecular mechanisms underlying STIM signaling has been dramatically improved through available X-ray crystal and solution NMR structures. This high-resolution structural data has revealed that intricate intramolecular and intermolecular protein-protein interactions are involved in converting STIMs from the quiescent to activation-competent states. This review article summarizes the current high resolution structural data on specific EF-hand, sterile α motif and coiled-coil interactions which drive STIM function in the activation of Orai1 channels. Further, the work discusses the effects of post-translational modifications on the structure and function of STIMs. Future structural studies on larger STIM:Orai complexes will be critical to fully defining the molecular bases for STIM function and how post-translational modifications influence these mechanisms.

Structural model of the SARS coronavirus E channel in LMPG micelles

Coronaviruses (CoV) cause common colds in humans, but are also responsible for the recent Severe Acute, and Middle East, respiratory syndromes (SARS and MERS, respectively). A promising approach for prevention are live attenuated vaccines (LAVs), some of which target the envelope (E) protein, which is a small membrane protein that forms ion channels. Unfortunately, detailed structural information is still limited for SARS-CoV E, and non-existent for other CoV E proteins. Herein, we report a structural model of a SARS-CoV E construct in LMPG micelles with, for the first time, unequivocal intermolecular NOEs. The model corresponding to the detergent-embedded region is consistent with previously obtained orientational restraints obtained in lipid bilayers and in vivo escape mutants. The C-terminal domain is mostly α-helical, and extramembrane intermolecular NOEs suggest interactions that may affect the TM channel conformation.

Mapping Antibody Epitopes by Solution NMR Spectroscopy: Practical Considerations

Identifying an epitope, the region of the antigen in contact with an antibody, is useful in both basic and pharmaceutical research, as well as in vaccine design. Solution NMR spectroscopy is particularly well suited to the residue level characterization of intermolecular interfaces, including antibody-antigen interactions, and thus to epitope identification. Here, we describe the use of NMR for residue level characterization of protein epitopes, focusing on experimental protocols and practical considerations, highlighting advantages and drawbacks of the approach.

Modern Technologies of Solution Nuclear Magnetic Resonance Spectroscopy for Three-dimensional Structure Determination of Proteins Open Avenues for Life Scientists

Nuclear magnetic resonance (NMR) spectroscopy is a powerful technique for structural studies of chemical compounds and biomolecules such as DNA and proteins. Since the NMR signal sensitively reflects the chemical environment and the dynamics of a nuclear spin, NMR experiments provide a wealth of structural and dynamic information about the molecule of interest at atomic resolution. In general, structural biology studies using NMR spectroscopy still require a reasonable understanding of the theory behind the technique and experience on how to recorded NMR data. Owing to the remarkable progress in the past decade, we can easily access suitable and popular analytical resources for NMR structure determination of proteins with high accuracy. Here, we describe the practical aspects, workflow and key points of modern NMR techniques used for solution structure determination of proteins. This review should aid NMR specialists aiming to develop new methods that accelerate the structure determination process, and open avenues for non-specialist and life scientists interested in using NMR spectroscopy to solve protein structures.

NMR structure and conformational dynamics of AtPDFL2.1, a defensin-like peptide from Arabidopsis thaliana

Plant defensins constitute the innate immune response against pathogens such as fungi and bacteria. Typical plant defensins are small, basic peptides that possess a characteristic three-dimensional fold stabilized by three or four disulfide bridges. In addition to known defensin genes, the Arabidopsis genome comprises >300 defensin-like genes coding for small cysteine-rich peptides. One of such genes encodes for AtPDFL2.1, a putative antifungal peptide of 55 amino acids, with six cysteine residues in its primary sequence. To understand the functional role of AtPDFL2.1, we carried out antifungal activity assays and determined its high-resolution three-dimensional structure using multidimensional solution NMR spectroscopy. We found that AtPDFL2.1 displays a strong inhibitory effect against Fusarium graminearum (IC₅₀≈4μM). This peptide folds in the canonical cysteine-stabilized αβ (CSαβ) motif, consisting of one α-helix and one triple-stranded antiparallel β-sheet stabilized by three disulfide bridges and a hydrophobic cluster of residues within its core where the α-helix packs tightly against the β-sheets. Nuclear spin relaxation measurements show that the structure of AtPDFL2.1 is essentially rigid, with the L3 loop located between β-strands 2 and 3 being more flexible and displaying conformational exchange. Interestingly, the dynamic features of loop L3 are conserved among defensins and are probably correlated to the antifungal and receptor binding activities.

Accessing long lived (1)H states via (2)H couplings

In this paper we demonstrate long-lived states involving a pair of chemically equivalent protons, with lifetimes ∼30 times T1 up to a total lifetime of ∼117s at high field (8.45T). This is demonstrated on trans-ethylene-d2 in solution, where magnetic inequivalence gives access to the long-lived states. It is shown that the remaining J-coupling between the two quadrupolar deuterium spins, JQQ, splits the conditions for optimally generating proton singlet states. Detailed simulations of the spin evolution are performed, shedding light on the coherent evolution during singlet-triplet conversion as well as on the incoherent evolution that causes relaxation. Subsequently, the simulations are compared with experimental results validating the theoretical insights. Possible applications include storage of hyperpolarization in the proton long-lived state. Of particular interest may be utilization of parahydrogen induced polarization to directly induce the examined long-lived states.

Nuclear magnetic resonance evidence for the dimer formation of beta amyloid peptide 1-42 in 1,1,1,3,3,3-hexafluoro-2-propanol

Alzheimer's disease involves accumulation of senile plaques in which filamentous aggregates of amyloid beta (Aβ) peptides are deposited. Recent studies demonstrate that oligomerization pathways of Aβ peptides may be complicated. To understand the mechanisms of Aβ(1-42) oligomer formation in more detail, we have established a method to produce (15)N-labeled Aβ(1-42) suited for nuclear magnetic resonance (NMR) studies. For physicochemical studies, the starting protein material should be solely monomeric and all Aβ aggregates must be removed. Here, we succeeded in fractionating a "precipitation-resistant" fraction of Aβ(1-42) from an "aggregation-prone" fraction by high-performance liquid chromatography (HPLC), even from bacterially overexpressed Aβ(1-42). However, both Aβ(1-42) fractions after 1,1,1,3,3,3-hexafluoro-2-propanol (HFIP) treatment formed amyloid fibrils. This indicates that the "aggregation seed" was not completely monomerized during HFIP treatment. In addition, Aβ(1-42) dissolved in HFIP was found to display a monomer-dimer equilibrium, as shown by two-dimensional (1)H-(15)N NMR. We demonstrated that the initial concentration of Aβ during the HFIP pretreatment altered the kinetic profiles of Aβ fibril formation in a thioflavin T fluorescence assay. The findings described here should ensure reproducible results when studying the Aβ(1-42) peptide.

Characterization of RNA-Protein Interactions: Lessons from Two RNA-Binding Proteins, SRSF1 and SRSF2

SR proteins are a class of RNA-binding proteins whose RNA-binding ability is required for both constitutive and alternative splicing. While members of the SR protein family were once thought to have redundant functions, in-depth biochemical analysis of their RNA-binding abilities has revealed distinct binding profiles for each SR protein, that often lead to either synergistic or antagonistic functions. SR protein family members SRSF1 and SRSF2 are two of the most highly studied RNA-binding proteins. Here we examine the various methods used to differentiate SRSF1 and SRSF2 RNA-binding ability. We discuss the benefits and type of information that can be determined using each method.

Effects of nucleotide binding to LmrA: A combined MAS-NMR and solution NMR study

ABC transporters are fascinating examples of fine-tuned molecular machines that use the energy from ATP hydrolysis to translocate a multitude of substrates across biological membranes. While structural details have emerged on many members of this large protein superfamily, a number of functional details are still under debate. High resolution structures yield valuable insights into protein function, but it is the combination of structural, functional and dynamic insights that facilitates a complete understanding of the workings of their complex molecular mechanisms. NMR is a technique well-suited to investigate proteins in atomic resolution while taking their dynamic properties into account. It thus nicely complements other structural techniques, such as X-ray crystallography, that have contributed high-resolution data to the architectural understanding of ABC transporters. Here, we describe the heterologous expression of LmrA, an ABC exporter from Lactococcus lactis, in Escherichia coli. This allows for more flexible isotope labeling for nuclear magnetic resonance (NMR) studies and the easy study of LmrA's multidrug resistance phenotype. We use a combination of solid-state magic angle spinning (MAS) on the reconstituted transporter and solution NMR on its isolated nucleotide binding domain to investigate consequences of nucleotide binding to LmrA. We find that nucleotide binding affects the protein globally, but that NMR is also able to pinpoint local dynamic effects to specific residues, such as the Walker A motif's conserved lysine residue.

Chemical shift assignments of zinc finger domain of methionine aminopeptidase 1 (MetAP1) from Homo sapiens

Methionine aminopeptidase Type I (MetAP1) cleaves the initiator methionine from about 70 % of all newly synthesized proteins in almost every living cell. Human MetAP1 is a two domain protein with a zinc finger on the N-terminus and a catalytic domain on the C-terminus. Here, we report the chemical shift assignments of the amino terminal zinc binding domain (ZBD) (1-83 residues) of the human MetAP1 derived by using advanced NMR spectroscopic methods. We were able to assign the chemical shifts of ZBD of MetAP1 nearly complete, which reveal two helical fragments involving residues P44-L49 (α1) and Q59-K82 (α2). The protein structure unfolds upon complex formation with the addition of 2 M excess EDTA, indicated by the appearance of amide resonances in the random coil chemical shift region of (15)NHSQC spectrum.

Fast conformational exchange between the sulfur-free and persulfide-bound rhodanese domain of E. coli YgaP

Rhodanese domains are abundant structural modules that catalyze the transfer of a sulfur atom from thiolsulfates to cyanide via formation of a covalent persulfide intermediate that is bound to an essential conserved cysteine residue. In this study, the three-dimensional structure of the rhodanese domain of YgaP from Escherichia coli was determined using solution NMR. A typical rhodanese domain fold was observed, as expected from the high homology with the catalytic domain of other sulfur transferases. The initial sulfur-transfer step and formation of the rhodanese persulfide intermediate were monitored by addition of sodium thiosulfate using two-dimensional (1)H-(15)N correlation spectroscopy. Discrete sharp signals were observed upon substrate addition, indicting fast exchange between sulfur-free and persulfide-intermediate forms. Residues exhibiting pronounced chemical shift changes were mapped to the structure, and included both substrate binding and surrounding residues.

Intracellular segment between transmembrane helices S0 and S1 of BK channel α subunit contains two amphipathic helices connected by a flexible loop

The BK channel, a tetrameric potassium channel with very high conductance, has a central role in numerous physiological functions. The BK channel can be activated by intracellular Ca(2+) and Mg(2+), as well as by membrane depolarization. Unlike other tetrameric potassium channels, the BK channel has seven transmembrane helices (S0-S6) including an extra helix S0. The intracellular segment between S0 and S1 (BK-IS1) is essential to BK channel functions and Asp99 in BK-IS1 is reported to be responsible for Mg(2+) coordination. In this study, BK-IS1 (44-113) was over-expressed using a bacterial system and purified in the presence of detergent micelles for multidimensional heteronuclear nuclear magnetic resonance (NMR) structural studies. Backbone resonance assignment and secondary structure analysis showed that BK-IS1 contains two amphipathic helices connected by a 36-residue loop. Amide (1)H-(15)N heteronuclear NOE analysis indicated that the loop is very flexible, while the two amphipathic helices are possibly stabilized through interaction with the membrane. A solution NMR-based titration assay of BK-IS1 was performed with various concentrations of Mg(2+). Two residues (Thr45 and Leu46) with chemical shift changes were observed but no, or very minor, chemical shift difference was observed for Asp99, indicating a possible site for binding divalent ions or other modulation partners.