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.     

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.     

<Superscript>15</Superscript>N and <Superscript>13</Superscript>C- SOFAST-HMQC editing enhances 3D-NOESY sensitivity in highly deuterated,selectively [<Superscript>1</Superscript>H,<Superscript>13</Superscript>C]-labeled proteins

The ongoing NMR method development effort strives for high quality multidimensional data with reduced collection time. Here, we apply ‘SOFAST-HMQC’ to frequency editing in 3D NOESY experiments and demonstrate the sensitivity benefits using highly deuterated and ¹⁵N, methyl labeled samples in H₂O. The experiments benefit from a combination of selective T ₁ relaxation (or L-optimized effect), from Ernst angle optimization and, in certain types of experiments, from using the mixing time for both NOE buildup and magnetization recovery. This effect enhances sensitivity by up to 2.4× at fast pulsing versus reference HMQC sequences of same overall length and water suppression characteristics. Representative experiments designed to address interesting protein NMR challenges are detailed. Editing capabilities are exploited with heteronuclear ¹⁵N,¹³C-edited, or with diagonal-free ¹³C aromatic/methyl-resolved 3D-SOFAST-HMQC–NOESY–HMQC. The latter experiment is used here to elucidate the methyl-aromatic NOE network in the hydrophobic core of the 19 kDa FliT-FliJ flagellar protein complex. Incorporation of fast pulsing to reference experiments such as 3D-NOESY–HMQC boosts digital resolution, simplifies the process of NOE assignment and helps to automate protein structure determination.     

J-UNIO protocol used for NMR structure determination of the 206-residue protein NP_346487.1 from <Emphasis Type="Italic">Streptococcus pneumoniae TIGR4</Emphasis>

The NMR structure of the 206-residue protein NP_346487.1 was determined with the J-UNIO protocol, which includes extensive automation of the structure determination. With input from three APSY-NMR experiments, UNIO-MATCH automatically yielded 77 % of the backbone assignments, which were interactively validated and extended to 97 %. With an input of the near-complete backbone assignments and three 3D heteronuclear-resolved [¹H,¹H]-NOESY spectra, automated side chain assignment with UNIO-ATNOS/ASCAN resulted in 77 % of the expected assignments, which was extended interactively to about 90 %. Automated NOE assignment and structure calculation with UNIO-ATNOS/CANDID in combination with CYANA was used for the structure determination of this two-domain protein. The individual domains in the NMR structure coincide closely with the crystal structure, and the NMR studies further imply that the two domains undergo restricted hinge motions relative to each other in solution. NP_346487.1 is so far the largest polypeptide chain to which the J-UNIO structure determination protocol has successfully been applied.     

<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.     

Efficient condensed-phase production of perdeuterated soluble and membrane proteins

Protein perdeuteration approaches have tremendous value in protein NMR studies, but are limited by the high cost of perdeuterated media. Here, we demonstrate that E. coli cultures expressing proteins using either the condensed single protein production method (cSPP), or conventional pET expression plasmids, can be condensed prior to protein expression, thereby providing high-quality ²H, ¹³C, ¹⁵N-enriched protein samples at 2.5–10% the cost of traditional methods. As an example of the value of such inexpensively-produced perdeuterated proteins, we produced ²H, ¹³C, ¹⁵N-enriched E. coli cold shock protein A (CspA) and EnvZb in 40× condensed phase media, and obtained NMR spectra suitable for 3D structure determination. The cSPP system was also used to produce ²H, ¹³C, ¹⁵N-enriched E. coli plasma membrane protein YaiZ and outer membrane protein X (OmpX) in condensed phase. NMR spectra can be obtained for these membrane proteins produced in the cSPP system following simple detergent extraction, without extensive purification or reconstitution. This allows a membrane protein’s structural and functional properties to be characterized prior to reconstitution, or as a probe of the effects of subsequent purification steps on the structural integrity of membrane proteins. We also provide a standardized protocol for production of perdeuterated proteins using the cSPP system. The 10–40 fold reduction in costs of fermentation media provided by using a condensed culture system opens the door to many new applications for perdeuterated proteins in spectroscopic and crystallographic studies.     

NMR characterization of HtpG,the <Emphasis Type="Italic">E. coli</Emphasis> Hsp90, using sparse labeling with <Superscript>13</Superscript>C-methyl alanine

A strategy for acquiring structural information from sparsely isotopically labeled large proteins is illustrated with an application to the E. coli heat-shock protein, HtpG (high temperature protein G), a 145 kDa dimer. It uses ¹³C-alanine methyl labeling in a perdeuterated background to take advantage of the sensitivity and resolution of Methyl-TROSY spectra, as well as the backbone-centered structural information from ¹H–¹³C residual dipolar couplings (RDCs) of alanine methyl groups. In all, 40 of the 47 expected crosspeaks were resolved and 36 gave RDC data. Assignments of crosspeaks were partially achieved by transferring assignments from those made on individual domains using triple resonance methods. However, these were incomplete and in many cases the transfer was ambiguous. A genetic algorithm search for consistency between predictions based on domain structures and measurements for chemical shifts and RDCs allowed 60% of the 40 resolved crosspeaks to be assigned with confidence. Chemical shift changes of these crosspeaks on adding an ATP analog to the apo-protein are shown to be consistent with structural changes expected on comparing previous crystal structures for apo- and complex- structures. RDCs collected on the assigned alanine methyl peaks are used to generate a new solution model for the apo-protein structure.     

Direct amide <Superscript>15</Superscript>N to <Superscript>13</Superscript>C transfers for solid-state assignment experiments in deuterated proteins

The assignment of protein backbone and side-chain NMR chemical shifts is the first step towards the characterization of protein structure. The recent introduction of proton detection in combination with fast MAS has opened up novel opportunities for assignment experiments. However, typical 3D sequential-assignment experiments using proton detection under fast MAS lead to signal intensities much smaller than the theoretically expected ones due to the low transfer efficiency of some of the steps. Here, we present a selective 3D experiment for deuterated and (amide) proton back-exchanged proteins where polarization is directly transferred from backbone nitrogen to selected backbone or sidechain carbons. The proposed pulse sequence uses only ¹H–¹⁵N cross-polarization (CP) transfers, which are, for deuterated proteins, about 30% more efficient than ¹H–¹³C CP transfers, and employs a dipolar version of the INEPT experiment for N–C transfer. By avoiding H^N–C (H^N stands for amide protons) and C–C CP transfers, we could achieve higher selectivity and increased signal intensities compared to other pulse sequences containing long-range CP transfers. The REDOR transfer is designed with an additional selective π pulse, which enables the selective transfer of the polarization to the desired ¹³C spins.     

<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>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.     

<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.     

<Superscript>1</Superscript>H, <Superscript>13</Superscript>C and <Superscript>15</Superscript>N chemical shift assignments of Saccharomyces cerevisiae type 1 thioredoxin in the oxidized state by solution NMR spectroscopy

Thioredoxins (Trx) are ubiquitous proteins that regulate several biochemical processes inside the cell. Trx is an important player, displaying oxidoreductase activity and helping to keep and regulate the oxidative state of the cellular environment. Trx also participates in the regulation of many cellular functions, such as DNA synthesis, protection against oxidative stress, cell cycle and signal transduction. The oxidized Trx is the target for another set of proteins, such as thioredoxin reductase (TrR), which used the reductive potential of NADPH. The oxidized state of Trx also plays important role in regulation of redox state in the cells. In this regard, the oxidized form of Trx is a putative conformer that contributes to the cellular redox environment. Here we report the chemical shift assignments (¹H, ¹³C and ¹⁵N) in solution at 15 °C. We also showed the secondary structure analysis of the oxidized form of yeast thioredoxin (yTrx1) as basis for future NMR studies of protein–target interactions and dynamics. The assignment was done at low concentration (200 µM) because it is important to keep intact the water cavity.     

<Superscript>1</Superscript>H, <Superscript>13</Superscript>C and <Superscript>15</Superscript>N resonance assignments and secondary structure analysis of translation initiation factor 1 from <Emphasis Type="Italic">Pseudomonas aeruginosa</Emphasis>

Pseudomonas aeruginosa is a Gram-negative opportunistic pathogen and a primary cause of infection in humans. P. aeruginosa can acquire resistance against multiple groups of antimicrobial agents, including β-lactams, aminoglycosides and fluoroquinolones, and multidrug resistance is increasing in this organism which makes treatment of the infections difficult and expensive. This has led to the unmet need for discovery of new compounds distinctly different from present antimicrobials. Protein synthesis is an essential metabolic process and a validated target for the development of new antibiotics. Translation initiation factor 1 from P. aeruginosa (Pa-IF1) is the smallest of the three initiation factors that acts to establish the 30S initiation complex to initiate translation during protein biosynthesis, and its structure is unknown. Here we report the ¹H, ¹³C and ¹⁵N chemical shift assignments of Pa-IF1 as the basis for NMR structure determination and interaction studies. Secondary structure analyses deduced from the NMR chemical shift data have identified five β-strands with an unusually extended β-strand at the C-terminal end of the protein and one short α-helix arranged in the sequential order β1–β2–β3–α1–β4–β5. This is further supported by ¹⁵N–{¹H} hetero NOEs. These secondary structure elements suggest the Pa-IF1 adopts the typical β-barrel structure and is composed of an oligomer-binding motif.     

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.     

Empirical correlation between protein backbone <Superscript>15</Superscript>N and <Superscript>13</Superscript>C secondary chemical shifts and its application to nitrogen chemical shift re-referencing

The linear analysis of chemical shifts (LACS) has provided a robust method for identifying and correcting ¹³C chemical shift referencing problems in data from protein NMR spectroscopy. Unlike other approaches, LACS does not require prior knowledge of the three-dimensional structure or inference of the secondary structure of the protein. It also does not require extensive assignment of the NMR data. We report here a way of extending the LACS approach to ¹⁵N NMR data from proteins, so as to enable the detection and correction of inconsistencies in chemical shift referencing for this nucleus. The approach is based on our finding that the secondary ¹⁵N chemical shift of the backbone nitrogen atom of residue i is strongly correlated with the secondary chemical shift difference (experimental minus random coil) between the alpha and beta carbons of residue i − 1. Thus once alpha and beta ¹³C chemical shifts are available (their difference is referencing error-free), the ¹⁵N referencing can be validated, and an appropriate offset correction can be derived. This approach can be implemented prior to a structure determination and can be used to analyze potential referencing problems in database data not associated with three-dimensional structure. Application of the LACS algorithm to the current BMRB protein chemical shift database, revealed that nearly 35% of the BMRB entries have δ ¹⁵N values mis-referenced by over 0.7 ppm and over 25% of them have δ ¹H^N values mis-referenced by over 0.12 ppm. One implication of the findings reported here is that a backbone ¹⁵N chemical shift provides a better indicator of the conformation of the preceding residue than of the residue itself. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Measurement of protein backbone <Superscript>13</Superscript>CO and <Superscript>15</Superscript>N relaxation dispersion at high resolution

The ³¹P NMR pressure response of guanine nucleotides bound to proteins has been studied in the past for characterizing the pressure perturbation of conformational equilibria. The pressure response of the ³¹P NMR chemical shifts of the phosphate groups of GMP, GDP, and GTP as well as the commonly used GTP analogs GppNHp, GppCH₂p and GTPγS was measured in the absence and presence of Mg²⁺-ions within a pressure range up to 200 MPa. The pressure dependence of chemical shifts is clearly non-linear. For all nucleotides a negative first order pressure coefficient B ₁ was determined indicating an upfield shift of the resonances with pressure. With exception of the α-phosphate group of Mg²⁺·GMP and Mg²⁺·GppNHp the second order pressure coefficients are positive. To describe the data of Mg²⁺·GppCH₂p and GTPγS a Taylor expansion of 3rd order is required. For distinguishing pH effects from pressure effects a complete pH titration set is presented for GMP, as well as GDP and GTP in absence and presence of Mg²⁺ ions using indirect referencing to DSS under identical experimental conditions. By a comparison between high pressure ³¹P NMR data on free Mg²⁺-GDP and Mg²⁺-GDP in complex with the proto-oncogene Ras we demonstrate that pressure induced changes in chemical shift are clearly different between both forms.     

NMR <Superscript>1</Superscript>H,<Superscript>13</Superscript>C, <Superscript>15</Superscript>N backbone and <Superscript>13</Superscript>C side chain resonance assignment of the G12C mutant of human K-Ras bound to GDP

K-Ras is a key driver of oncogenesis, accounting for approximately 80% of Ras-driven human cancers. The small GTPase cycles between an inactive, GDP-bound and an active, GTP-bound state, regulated by guanine nucleotide exchange factors and GTPase activating proteins, respectively. Activated K-Ras regulates cell proliferation, differentiation and survival by signaling through several effector pathways, including Raf-MAPK. Oncogenic mutations that impair the GTPase activity of K-Ras result in a hyperactivated state, leading to uncontrolled cellular proliferation and tumorogenesis. A cysteine mutation at glycine 12 is commonly found in K-Ras associated cancers, and has become a recent focus for therapeutic intervention. We report here ¹H^{N, 15}N, and ¹³C resonance assignments for the 19.3 kDa (aa 1–169) human K-Ras protein harboring an oncogenic G12C mutation in the GDP-bound form (K-RAS^G12C-GDP), using heteronuclear, multidimensional NMR spectroscopy. Backbone ¹H–¹⁵N correlations have been assigned for all non-proline residues, except for the first methionine residue.     

<Superscript>1</Superscript>H, <Superscript>13</Superscript>C, <Superscript>15</Superscript>N backbone assignment of the human heat-labile enterotoxin B-pentamer and chemical shift mapping of neolactotetraose binding

The major virulence factor of enterotoxigenic Escherichia coli is the heat-labile enterotoxin (LT), an AB₅ toxin closely related to the cholera toxin. LT consists of six subunits, the catalytically active A-subunit and five B-subunits arranged as a pentameric ring (LTB), which enable the toxin to bind to the epithelial cells in the intestinal lumen. LTB has two recognized binding sites; the primary binding site is responsible for anchoring the toxin to its main receptor, the GM₁-ganglioside, while the secondary binding site recognizes blood group antigens. Herein, we report the ¹H, ¹³C, ¹⁵N main chain assignment of LTB from human isolates (hLTB; 103 a.a. per subunit, with a total molecular mass of 58.5 kDa). The secondary structure was predicted based on ¹³C′, ¹³C^{α, 13}C^{β, 1}H^N and ¹⁵N chemical shifts and compared to a published crystal structure of LTB. Neolactotetraose (NEO) was titrated to hLTB and chemical shift perturbations were measured. The chemical shift perturbations were mapped onto the crystal structure, confirming that NEO binds to the primary binding site of hLTB and competes with GM₁-binding. Our new data further lend support to the hypothesis that binding at the primary binding site is transmitted to the secondary binding site of the toxin, where it may influence the binding to blood group antigens.     

Near-complete <Superscript>1</Superscript>H, <Superscript>13</Superscript>C, <Superscript>15</Superscript>N resonance assignments of dimethylsulfoxide-denatured TGFBIp FAS1-4 A546T

The transforming growth factor beta induced protein (TGFBIp) is a major protein component of the human cornea. Mutations occurring in TGFBIp may cause corneal dystrophies, which ultimately lead to loss of vision. The majority of the disease-causing mutations are located in the C-terminal domain of TGFBIp, referred as the fourth fascilin-1 (FAS1-4) domain. In the present study the FAS1-4 Ala546Thr, a mutation that causes lattice corneal dystrophy, was investigated in dimethylsulfoxide using liquid-state NMR spectroscopy, to enable H/D exchange strategies for identification of the core formed in mature fibrils. Isotope-labeled fibrillated FAS1-4 A546T was dissolved in a ternary mixture 95/4/1 v/v/v% dimethylsulfoxide/water/trifluoroacetic acid, to obtain and assign a reference 2D ¹H–¹⁵N HSQC spectrum for the H/D exchange analysis. Here, we report the near-complete assignments of backbone and aliphatic side chain ¹H, ¹³C and ¹⁵N resonances for unfolded FAS1-4 A546T at 25 °C.