APSY-NMR for protein backbone assignment in high-throughput structural biology

A standard set of three APSY-NMR experiments has been used in daily practice to obtain polypeptide backbone NMR assignments in globular proteins with sizes up to about 150 residues, which had been identified as targets for structure determination by the Joint Center for Structural Genomics (JCSG) under the auspices of the Protein Structure Initiative (PSI). In a representative sample of 30 proteins, initial fully automated data analysis with the software UNIO-MATCH-2014 yielded complete or partial assignments for over 90 % of the residues. For most proteins the APSY data acquisition was completed in less than 30 h. The results of the automated procedure provided a basis for efficient interactive validation and extension to near-completion of the assignments by reference to the same 3D heteronuclear-resolved [¹H,¹H]-NOESY spectra that were subsequently used for the collection of conformational constraints. High-quality structures were obtained for all 30 proteins, using the J-UNIO protocol, which includes extensive automation of NMR structure determination.     

Reliability of exclusively NOESY-based automated resonance assignment and structure determination of proteins

Protein structure determination by NMR can in principle be speeded up both by reducing the measurement time on the NMR spectrometer and by a more efficient analysis of the spectra. Here we study the reliability of protein structure determination based on a single type of spectra, namely nuclear Overhauser effect spectroscopy (NOESY), using a fully automated procedure for the sequence-specific resonance assignment with the recently introduced FLYA algorithm, followed by combined automated NOE distance restraint assignment and structure calculation with CYANA. This NOESY-FLYA method was applied to eight proteins with 63–160 residues for which resonance assignments and solution structures had previously been determined by the Northeast Structural Genomics Consortium (NESG), and unrefined and refined NOESY data sets have been made available for the Critical Assessment of Automated Structure Determination of Proteins by NMR project. Using only peak lists from three-dimensional ¹³C- or ¹⁵N-resolved NOESY spectra as input, the FLYA algorithm yielded for the eight proteins 91–98 % correct backbone and side-chain assignments if manually refined peak lists are used, and 64–96 % correct assignments based on raw peak lists. Subsequent structure calculations with CYANA then produced structures with root-mean-square deviation (RMSD) values to the manually determined reference structures of 0.8–2.0 Å if refined peak lists are used. With raw peak lists, calculations for 4 proteins converged resulting in RMSDs to the reference structure of 0.8–2.8 Å, whereas no convergence was obtained for the four other proteins (two of which did already not converge with the correct manual resonance assignments given as input). These results show that, given high-quality experimental NOESY peak lists, the chemical shift assignments can be uncovered, without any recourse to traditional through-bond type assignment experiments, to an extent that is sufficient for calculating accurate three-dimensional structures.     

Automated solid-state NMR resonance assignment of protein microcrystals and amyloids

Solid-state NMR is an emerging structure determination technique for crystalline and non-crystalline protein assemblies, e.g., amyloids. Resonance assignment constitutes the first and often very time-consuming step to a structure. We present ssFLYA, a generally applicable algorithm for automatic assignment of protein solid-state NMR spectra. Application to microcrystals of ubiquitin and the Ure2 prion C-terminal domain, as well as amyloids of HET-s(218–289) and α-synuclein yielded 88–97 % correctness for the backbone and side-chain assignments that are classified as self-consistent by the algorithm, and 77–90 % correctness if also assignments classified as tentative by the algorithm are included.     

<Superscript>1</Superscript>H, <Superscript>13</Superscript>C and <Superscript>15</Superscript>N backbone NMR chemical shift assignments of the C-terminal P4 domain of Ahnak

Ahnak is a ~?700 kDa polypeptide that was originally identified as a tumour-related nuclear phosphoprotein, but later recognized to play a variety of diverse physiological roles related to cell architecture and migration. A critical function of Ahnak is modulation of Ca²⁺ signaling in cardiomyocytes by interacting with the β subunit of the L-type Ca²⁺ channel (Ca_V1.2). Previous studies have identified the C-terminal region of Ahnak, designated as P3 and P4 domains, as a key mediator of its functional activity. We report here the nearly complete ¹H, ¹³C and ¹⁵N backbone NMR chemical shift assignments of the 11 kDa C-terminal P4 domain of Ahnak. This study lays the foundations for future investigations of functional dynamics, structure determination and interaction site mapping of the Ca_V1.2-Ahnak complex.     

Fully automated high-quality NMR structure determination of small <Superscript>2</Superscript>H-enriched proteins

Determination of high-quality small protein structures by nuclear magnetic resonance (NMR) methods generally requires acquisition and analysis of an extensive set of structural constraints. The process generally demands extensive backbone and sidechain resonance assignments, and weeks or even months of data collection and interpretation. Here we demonstrate rapid and high-quality protein NMR structure generation using CS-Rosetta with a perdeuterated protein sample made at a significantly reduced cost using new bacterial culture condensation methods. Our strategy provides the basis for a high-throughput approach for routine, rapid, high-quality structure determination of small proteins. As an example, we demonstrate the determination of a high-quality 3D structure of a small 8 kDa protein, E. coli cold shock protein A (CspA), using <4 days of data collection and fully automated data analysis methods together with CS-Rosetta. The resulting CspA structure is highly converged and in excellent agreement with the published crystal structure, with a backbone RMSD value of 0.5 Å, an all atom RMSD value of 1.2 Å to the crystal structure for well-defined regions, and RMSD value of 1.1 Å to crystal structure for core, non-solvent exposed sidechain atoms. Cross validation of the structure with ¹⁵N- and ¹³C-edited NOESY data obtained with a perdeuterated ¹⁵N, ¹³C-enriched ¹³CH₃ methyl protonated CspA sample confirms that essentially all of these independently-interpreted NOE-based constraints are already satisfied in each of the 10 CS-Rosetta structures. By these criteria, the CS-Rosetta structure generated by fully automated analysis of data for a perdeuterated sample provides an accurate structure of CspA. This represents a general approach for rapid, automated structure determination of small proteins by NMR.     

Automatic <Superscript>13</Superscript>C chemical shift reference correction for unassigned protein NMR spectra

Poor chemical shift referencing, especially for ¹³C in protein Nuclear Magnetic Resonance (NMR) experiments, fundamentally limits and even prevents effective study of biomacromolecules via NMR, including protein structure determination and analysis of protein dynamics. To solve this problem, we constructed a Bayesian probabilistic framework that circumvents the limitations of previous reference correction methods that required protein resonance assignment and/or three-dimensional protein structure. Our algorithm named Bayesian Model Optimized Reference Correction (BaMORC) can detect and correct ¹³C chemical shift referencing errors before the protein resonance assignment step of analysis and without three-dimensional structure. By combining the BaMORC methodology with a new intra-peaklist grouping algorithm, we created a combined method called Unassigned BaMORC that utilizes only unassigned experimental peak lists and the amino acid sequence. Unassigned BaMORC kept all experimental three-dimensional HN(CO)CACB-type peak lists tested within ±?0.4 ppm of the correct ¹³C reference value. On a much larger unassigned chemical shift test set, the base method kept ¹³C chemical shift referencing errors to within ±?0.45 ppm at a 90% confidence interval. With chemical shift assignments, Assigned BaMORC can detect and correct ¹³C chemical shift referencing errors to within ±?0.22 at a 90% confidence interval. Therefore, Unassigned BaMORC can correct ¹³C chemical shift referencing errors when it will have the most impact, right before protein resonance assignment and other downstream analyses are started. After assignment, chemical shift reference correction can be further refined with Assigned BaMORC. These new methods will allow non-NMR experts to detect and correct ¹³C referencing error at critical early data analysis steps, lowering the bar of NMR expertise required for effective protein NMR analysis.     

Automated combined assignment of NOESY spectra and three-dimensional protein structure determination

A procedure for automated protein structure determination is presented that is based on an iterative procedure during which the NOESY peak list assignment and the structure calculation are performed simultaneously. The input consists of a list of NOESY peak positions and a list of chemical shifts as obtained from sequence-specific resonance assignment. For the present applications of this approach the previously introduced NOAH routine was implemented in the distance geometry program DIANA. As an illustration, experimental 2D and 3D NOESY cross-peak lists of six proteins have been analyzed, for which complete sequence-specific ¹H assignments are available for the polypeptide backbone and the amino acid side chains. The automated method assigned 70–90% of all NOESY cross peaks, which is on average 10% less than with the interactive approach, and only between 0.8% and 2.4% of the automatically assigned peaks had a different assignment than in the corresponding manually assigned peak lists. The structures obtained with NOAH/DIANA are in close agreement with those from manually assigned peak lists, and with both approaches the residual constraint violations correspond to high-quality NMR structure determinations. Systematic comparisons of the bundles of conformers that represent corresponding automatically and interactively determined structures document the absence of significant bias in either approach, indicating that an important step has been made towards automation of structure determination from NMR spectra.     

<Superscript>1</Superscript>H, <Superscript>13</Superscript>C,and <Superscript>15</Superscript>N resonance assignments of FAS1-IV domain of human periostin,a component of extracellular matrix proteins

Periostin, an extracellular matrix protein, is secreted by fibroblasts and is overexpressed in various types of cancers. The four internal repeat fasciclin 1 (FAS1) domains of human periostin play crucial roles in promoting tumor metastasis and progression via interaction with cell surface integrins. Among four FAS1 domains of human periostin, the fourth FAS1 domain (FAS1-IV) was prepared for NMR study, since only FAS1-IV was highly soluble, and showed a well-dispersed 2D ¹H-¹⁵N HSQC spectrum. Here, we report nearly complete backbone and side chain resonance assignments and a secondary structural analysis of the FAS1-IV domain as first steps toward the structure determination of FAS1-IV of human periostin.     

Automated amino acid side-chain NMR assignment of proteins using <Superscript>13</Superscript>C- and <Superscript>15</Superscript>N-resolved 3D [<Superscript>1</Superscript>H,<Superscript>1</Superscript>H]-NOESY

ASCAN is a new algorithm for automatic sequence-specific NMR assignment of amino acid side-chains in proteins, which uses as input the primary structure of the protein, chemical shift lists of (1)H(N), (15)N, (13)C(alpha), (13)C(beta) and possibly (1)H(alpha) from the previous polypeptide backbone assignment, and one or several 3D (13)C- or (15)N-resolved [(1)H,(1)H]-NOESY spectra. ASCAN has also been laid out for the use of TOCSY-type data sets as supplementary input. The program assigns new resonances based on comparison of the NMR signals expected from the chemical structure with the experimentally observed NOESY peak patterns. The core parts of the algorithm are a procedure for generating expected peak positions, which is based on variable combinations of assigned and unassigned resonances that arise for the different amino acid types during the assignment procedure, and a corresponding set of acceptance criteria for assignments based on the NMR experiments used. Expected patterns of NOESY cross peaks involving unassigned resonances are generated using the list of previously assigned resonances, and tentative chemical shift values for the unassigned signals taken from the BMRB statistics for globular proteins. Use of this approach with the 101-amino acid residue protein FimD(25-125) resulted in 84% of the hydrogen atoms and their covalently bound heavy atoms being assigned with a correctness rate of 90%. Use of these side-chain assignments as input for automated NOE assignment and structure calculation with the ATNOS/CANDID/DYANA program suite yielded structure bundles of comparable quality, in terms of precision and accuracy of the atomic coordinates, as those of a reference structure determined with interactive assignment procedures. A rationale for the high quality of the ASCAN-based structure determination results from an analysis of the distribution of the assigned side chains, which revealed near-complete assignments in the core of the protein, with most of the incompletely assigned residues located at or near the protein surface.     

The J-UNIO protocol for automated protein structure determination by NMR in solution

The J-UNIO (JCSG protocol using the software UNIO) procedure for automated protein structure determination by NMR in solution is introduced. In the present implementation, J-UNIO makes use of APSY-NMR spectroscopy, 3D heteronuclear-resolved [(1)H,(1)H]-NOESY experiments, and the software UNIO. Applications with proteins from the JCSG target list with sizes up to 150 residues showed that the procedure is highly robust and efficient. In all instances the correct polypeptide fold was obtained in the first round of automated data analysis and structure calculation. After interactive validation of the data obtained from the automated routine, the quality of the final structures was comparable to results from interactive structure determination. Special advantages are that the NMR data have been recorded with 6-10 days of instrument time per protein, that there is only a single step of chemical shift adjustments to relate the backbone signals in the APSY-NMR spectra with the corresponding backbone signals in the NOESY spectra, and that the NOE-based amino acid side chain chemical shift assignments are automatically focused on those residues that are heavily weighted in the structure calculation. The individual working steps of J-UNIO are illustrated with the structure determination of the protein YP_926445.1 from Shewanella amazonensis, and the results obtained with 17 JCSG targets are critically evaluated.     

<Superscript>1</Superscript>H, <Superscript>15</Superscript>N, <Superscript>13</Superscript>C backbone resonance assignment of the C-terminal domain of enzyme I from <Emphasis Type="Italic">Thermoanaerobacter tengcongensis</Emphasis>

Phosphoenolpyruvate binding to the C-terminal domain (EIC) of enzyme I of the bacterial phosphotransferase system (PTS) initiates a phosphorylation cascade that results in sugar translocation across the cell membrane and controls a large number of essential pathways in bacterial metabolism. EIC undergoes an expanded to compact conformational equilibrium that is regulated by ligand binding and determines the phosphorylation state of the overall PTS. Here, we report the backbone ¹H, ¹⁵N and ¹³C chemical shift assignments of the 70 kDa EIC dimer from the thermophilic bacterium Thermoanaerobacter tengcongensis. Assignments were obtained at 70 °C by heteronuclear multidimensional NMR spectroscopy. In total, 90% of all backbone resonances were assigned, with 264 out of a possible 299 residues assigned in the ¹H–¹⁵N TROSY spectrum. The secondary structure predicted from the assigned backbone resonance using the program TALOS+ is in good agreement with the X-ray crystal structure of T. tengcongensis EIC. The reported assignments will allow detailed structural and thermodynamic investigations on the coupling between ligand binding and conformational dynamics in EIC.     

Backbone <Superscript>1</Superscript>H, <Superscript>13</Superscript>C and <Superscript>15</Superscript>N chemical shift assignment of full-length human uracil DNA glycosylase UNG2

Human uracil N-glycosylase isoform 2—UNG2 consists of an N-terminal intrinsically disordered regulatory domain (UNG2 residues 1–92, 9.3 kDa) and a C-terminal structured catalytic domain (UNG2 residues 93–313, 25.1 kDa). Here, we report the backbone ¹H, ¹³C, and ¹⁵N chemical shift assignment as well as secondary structure analysis of the N-and C-terminal domains of UNG2 representing the full-length UNG2 protein.     

Fast Assignment of <Superscript>15</Superscript>N-HSQC Peaks using High-Resolution 3D HNcocaNH Experiments with Non-Uniform Sampling

We describe an efficient NMR triple resonance approach for fast assignment of backbone amide resonance peaks in the ¹⁵N-HSQC spectrum. The exceptionally high resolutions achieved in the 3D HncocaNH and hNcocaNH experiments together with non-uniform sampling facilitate error-free sequential connection of backbone amides. Data required for the complete backbone amide assignment of the 56-residue protein GB1 domain were obtained in 14 h. Data analysis was vastly streamlined using a ‘backbone NH walk’ method to determine sequential connectivities without the need for ¹³C chemical shifts comparison. Amino acid residues in the sequentially connected NH chains are classified into two groups by a simple variation of the NMR pulse sequence, and the resulting ‘ZeBra’ stripe patterns are useful for mapping these chains to the protein sequence. In addition to resolving ambiguous assignments derived from conventional backbone experiments, this approach can be employed to rapidly assign small proteins or flexible regions in larger proteins, and to transfer assignments to mutant proteins or proteins in different ligand-binding states.     

Protein side-chain resonance assignment and NOE assignment using RDC-defined backbones without TOCSY data

One bottleneck in NMR structure determination lies in the laborious and time-consuming process of side-chain resonance and NOE assignments. Compared to the well-studied backbone resonance assignment problem, automated side-chain resonance and NOE assignments are relatively less explored. Most NOE assignment algorithms require nearly complete side-chain resonance assignments from a series of through-bond experiments such as HCCH-TOCSY or HCCCONH. Unfortunately, these TOCSY experiments perform poorly on large proteins. To overcome this deficiency, we present a novel algorithm, called Nasca (NOE Assignment and Side-Chain Assignment), to automate both side-chain resonance and NOE assignments and to perform high-resolution protein structure determination in the absence of any explicit through-bond experiment to facilitate side-chain resonance assignment, such as HCCH-TOCSY. After casting the assignment problem into a Markov Random Field (MRF), Nasca extends and applies combinatorial protein design algorithms to compute optimal assignments that best interpret the NMR data. The MRF captures the contact map information of the protein derived from NOESY spectra, exploits the backbone structural information determined by RDCs, and considers all possible side-chain rotamers. The complexity of the combinatorial search is reduced by using a dead-end elimination (DEE) algorithm, which prunes side-chain resonance assignments that are provably not part of the optimal solution. Then an A* search algorithm is employed to find a set of optimal side-chain resonance assignments that best fit the NMR data. These side-chain resonance assignments are then used to resolve the NOE assignment ambiguity and compute high-resolution protein structures. Tests on five proteins show that Nasca assigns resonances for more than 90% of side-chain protons, and achieves about 80% correct assignments. The final structures computed using the NOE distance restraints assigned by Nasca have backbone RMSD 0.8–1.5 Å from the reference structures determined by traditional NMR approaches.     

Residual backbone and side-chain <Superscript>13</Superscript>C and <Superscript>15</Superscript>N resonance assignments of the intrinsic transmembrane light-harvesting 2 protein complex by solid-state Magic Angle Spinning NMR spectroscopy

This study reports the sequence specific chemical shifts assignments for 76 residues of the 94 residues containing monomeric unit of the photosynthetic light-harvesting 2 transmembrane protein complex from Rhodopseudomonas acidophila strain 10050, using Magic Angle Spinning (MAS) NMR in combination with extensive and selective biosynthetic isotope labeling methods. The sequence specific chemical shifts assignment is an essential step for structure determination by MAS NMR. Assignments have been performed on the basis of 2-dimensional proton-driven spin diffusion ¹³C–¹³C correlation experiments with mixing times of 20 and 500 ms and band selective ¹³C–¹⁵N correlation spectroscopy on a series of site-specific biosynthetically labeled samples. The decreased line width and the reduced number of correlation signals of the selectively labeled samples with respect to the uniformly labeled samples enable to resolve the narrowly distributed correlation signals of the backbone carbons and nitrogens involved in the long -helical transmembrane segments. Inter-space correlations between nearby residues and between residues and the labeled BChl a cofactors, provided by the ¹³C–¹³C correlation experiments using a 500 ms spin diffusion period, are used to arrive at sequence specific chemical shift assignments for many residues in the protein complex. In this way it is demonstrated that MAS NMR methods combined with site-specific biosynthetic isotope labeling can be used for sequence specific assignment of the NMR response of transmembrane proteins.     

<Superscript>1</Superscript>H, <Superscript>15</Superscript>N, <Superscript>13</Superscript>C resonance assignments for <Emphasis Type="Italic">Saccharomyces cerevisiae</Emphasis> Rad23 UBL domain

Rad23 functions in nucleotide excision repair and proteasome-mediated protein degradation. It has four distinct structural domains that are connected by flexible linker regions, including an N-terminal ubiquitin-like (UBL) domain that binds proteasomes. We report in this NMR study the ¹H, ¹⁵N and ¹³C resonance assignments for the backbone and side chain atoms of the Rad23 UBL domain (Rad23^UBL) with BioMagResBank accession number 25825. We find that a Rad23 proline amino acid (P20) located in a loop undergoes isomerization. The secondary structural elements predicted from the NMR data fit well to that of the Rad23^UBL when complexed with E4 ubiquitin ligase Ufd2, as reported in a crystallographic structure. These complete assignments can be used to study the protein dynamics of the Rad23^UBL and its interaction of with other ubiquitin receptors or proteasome subunits.     

<Superscript>1</Superscript>H, <Superscript>15</Superscript>N and <Superscript>13</Superscript>C resonance assignments of the C-terminal domain of <Emphasis Type="Italic">Vibrio cholerae</Emphasis> TolA protein

Vibrio cholerae is the bacterial causative agent of the human disease cholera. Non-pathogenic bacterium can be converted to pathogenic following infection by a filamentous phage, CTXΦ, that carries the cholera toxin encoding genes. A crucial step during phage infection requires a direct interaction between the CTXΦ minor coat protein (pIII^CTX) and the C-terminal domain of V. cholerae TolA protein (TolAIIIvc). In order to get a better understanding of TolA function during the infection process, we have initiated a study of the V. cholerae TolAIII domain by 2D and 3D heteronuclear NMR. With the exception of the His-tag (H123–H128), 97 % of backbone ¹H, ¹⁵N and ¹³C resonances were assigned and the side chain assignments for 92 % of the protein were obtained (BMRB deposit with accession number 25689).     

<Superscript>1</Superscript>H, <Superscript>15</Superscript>N,and <Superscript>13</Superscript>C resonance assignments of <Emphasis Type="Italic">Staphylococcus aureus</Emphasis> extracellular adherence protein domain 4

The pathogenic bacterium Staphylococcus aureus has evolved to actively evade many aspects of the human innate immune system by expressing a series of secreted inhibitory proteins. Among these, the extracellular adherence protein (Eap) has been shown to inhibit the classical and lectin pathways of the complement system. By binding to complement component C4b, Eap is able to inhibit formation of the CP/LP C3 pro-convertase. Secreted full-length, mature Eap consists of four ~98 residue domains, all of which adopt a similar beta-grasp fold, and are connected through a short linker region. Through multiple biochemical approaches, it has been determined that the third and fourth domains of Eap are responsible for C4b binding. Here we report the backbone and side-chain resonance assignments of the 11.3 kDa fourth domain of Eap. The assignment data has been deposited in the BMRB database under the accession number 26726.