Biosynthetic site-specific <Superscript>13</Superscript> C labeling of the light-harvesting 2 protein complex: A model for solid state NMR structure determination of transmembrane proteins

Partly biosynthetic site-directed isotopically ¹³C enriched photosynthetic light-harvesting 2(LH2) complexes have been prepared from Rhodopseudomonas acidophila strain 10050 by using chemically labeled [1,2,3,4–¹³C], [1,4–¹³C] and [2,3–¹³C] succinic acid as a precursor in the growth medium. Two-dimensional proton driven spin diffusion (PDSD) solid state NMR correlation spectroscopy has been used to trace each individual ¹³C isotope from the labeled succinic acid precursor to its destination into the protein and into the embedded major light-absorbing bacteriochlorophyll cofactors. For both the residues of the protein and for the cofactors distinct labeling patterns have been deduced, for protein complexes prepared from [1,4–¹³C]-succinic acid or [2,3–¹³C]-succinic labeled media. All residues, except isoleucine and leucine, have been labeled almost homogeneously by the succinic acid precursor. Carbonyl carbons in the protein backbone were labeled by [1,4–¹³C]-succinic acid, while the C and C carbons of the residues were labeled by [2,3 ¹³C]-succinic acid. Leucine and isoleucine residues were labeled using a uniformly labeled amino acid mixture in the medium. The pattern labeling yields an increase of the resolution and less spectral crowding. The partial labeling technique in combination with conventional solid state NMR methods at ultra high magnetic fields provides an attractive route to resolve chemical shifts for -helical transmembrane protein structures.     

Assigning large proteins in the solid state: a MAS NMR resonance assignment strategy using selectively and extensively <Superscript>13</Superscript>C-labelled proteins

In recent years, solid-state magic-angle spinning nuclear magnetic resonance spectroscopy (MAS NMR) has been growing into an important technique to study the structure of membrane proteins, amyloid fibrils and other protein preparations which do not form crystals or are insoluble. Currently, a key bottleneck is the assignment process due to the absence of the resolving power of proton chemical shifts. Particularly for large proteins (approximately >150 residues) it is difficult to obtain a full set of resonance assignments. In order to address this problem, we present an assignment method based upon samples prepared using [1,3-¹³C]- and [2-¹³C]-glycerol as the sole carbon source in the bacterial growth medium (so-called selectively and extensively labelled protein). Such samples give rise to higher quality spectra than uniformly [¹³C]-labelled protein samples, and have previously been used to obtain long-range restraints for use in structure calculations. Our method exploits the characteristic cross-peak patterns observed for the different amino acid types in ¹³C-¹³C correlation and 3D NCACX and NCOCX spectra. An in-depth analysis of the patterns and how they can be used to aid assignment is presented, using spectra of the chicken α-spectrin SH3 domain (62 residues), αB-crystallin (175 residues) and outer membrane protein G (OmpG, 281 residues) as examples. Using this procedure, over 90% of the Cα, Cβ, C′ and N resonances in the core domain of αB-crystallin and around 73% in the flanking domains could be assigned (excluding 24 residues at the extreme termini of the protein). Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Backbone dynamics of bacteriorhodopsin as studied by <Superscript>13</Superscript>C solid-state NMR spectroscopy

The surface dynamics of bacteriorhodopsin was examined by measurements of site-specific ¹³C–¹H dipolar couplings in [3-¹³C]Ala-labeled bacteriorhodopsin. Motions of slow or intermediate frequency (correlation time <50 µs) scale down ¹³C–¹H dipolar couplings according to the motional amplitude. The two-dimensional dipolar and chemical shift (DIPSHIFT) correlation technique was utilized to obtain the dipolar coupling strength for each resolved peak in the ¹³C MAS solid-state NMR spectrum, providing the molecular order parameter of the respective site. In addition to the rotation of the Ala methyl group, which scales the dipolar coupling to 1/3 of the rigid limit value, fluctuations of the C–C vector result in additional motional averaging. Typical order parameters measured for mobile sites in bacteriorhodopsin are between 0.25 and 0.29. These can be assigned to Ala103 of the C–D loop and Ala235 at the C-terminal -helix protruded from the membrane surface, and Ala196 of the F–G loop, as well as to Ala228 and Ala233 of the C-terminal -helix and Ala51 from the transmembrane -helix. Such order parameters departing significantly from the value of 0.33 for rotating methyl groups are obviously direct evidence for the presence of fluctuation motions of the Ala C–C vectors of intact preparations of fully hydrated, wild-type bacteriorhodopsin at ambient temperature. The order parameter for Ala160 from the expectantly more flexible E–F loop, however, is unavailable under highest-field NMR conditions, probably because increased chemical shift anisotropy together with intrinsic fluctuation motions result in an unresolved ¹³C NMR signal.     

Optimization of amino acid type-specific <Superscript>13</Superscript>C and <Superscript>15</Superscript>N labeling for the backbone assignment of membrane proteins by solution- and solid-state NMR with the UPLABEL algorithm

We present a computational method for finding optimal labeling patterns for the backbone assignment of membrane proteins and other large proteins that cannot be assigned by conventional strategies. Following the approach of Kainosho and Tsuji (Biochemistry 21:6273–6279 (1982)), types of amino acids are labeled with ¹³C or/and ¹⁵N such that cross peaks between ¹³CO(i – 1) and ¹⁵NH(i) result only for pairs of sequentially adjacent amino acids of which the first is labeled with ¹³C and the second with ¹⁵N. In this way, unambiguous sequence-specific assignments can be obtained for unique pairs of amino acids that occur exactly once in the sequence of the protein. To be practical, it is crucial to limit the number of differently labeled protein samples that have to be prepared while obtaining an optimal extent of labeled unique amino acid pairs. Our computer algorithm UPLABEL for optimal unique pair labeling, implemented in the program CYANA and in a standalone program, and also available through a web portal, uses combinatorial optimization to find for a given amino acid sequence labeling patterns that maximize the number of unique pair assignments with a minimal number of differently labeled protein samples. Various auxiliary conditions, including labeled amino acid availability and price, previously known partial assignments, and sequence regions of particular interest can be taken into account when determining optimal amino acid type-specific labeling patterns. The method is illustrated for the assignment of the human G-protein coupled receptor bradykinin B2 (B₂R) and applied as a starting point for the backbone assignment of the membrane protein proteorhodopsin.     

<Superscript>31</Superscript>P NMR correlation maps of <Superscript>18</Superscript>O/<Superscript>16</Superscript>O chemical shift isotopic effects for phosphometabolite labeling studies

Intramolecular correlations among the ¹⁸O-labels of metabolic oligophosphates, mapped by J-decoupled ³¹P NMR 2D chemical shift correlation spectroscopy, impart stringent constraints to the ¹⁸O-isotope distributions over the whole oligophosphate moiety. The multiple deduced correlations of isotopic labels enable determination of site-specific fractional isotope enrichments and unravel the isotopologue statistics. This approach ensures accurate determination of ¹⁸O-labeling rates of phosphometabolites, critical in biochemical energy conversion and metabolic flux transmission. The biological usefulness of the J-decoupled ³¹P NMR 2D chemical shift correlation maps was validated on adenosine tri-phosphate fractionally ¹⁸O labeled in perfused mammalian hearts.     

Overexpression of a homogeneous oligosaccharide with <Superscript>13</Superscript>C labeling by genetically engineered yeast strain

This report describes a novel method for overexpression of ¹³C-labeled oligosaccharides using genetically engineered Saccharomyces cerevisiae cells, in which a homogeneous high-mannose-type oligosaccharide accumulates because of deletions of genes encoding three enzymes involved in the processing pathway of asparagine-linked oligosaccharides in the Golgi complex. Using uniformly ¹³C-labeled glucose as the sole carbon source in the culture medium of these engineered yeast cells, high yields of the isotopically labeled Man₈GlcNAc₂ oligosaccharide could be successfully harvested from glycoprotein extracts of the cells. Furthermore, ¹³C labeling at selected positions of the sugar residues in the oligosaccharide could be achieved using a site-specific ¹³C-enriched glucose as the metabolic precursor, facilitating NMR spectral assignments. The ¹³C-labeling method presented provides the technical basis for NMR analyses of structures, dynamics, and interactions of larger, branched oligosaccharides.     

<Superscript>1</Superscript>H, <Superscript>13</Superscript>C and <Superscript>15</Superscript>N NMR assignments of a calcium-binding protein from <Emphasis Type="Italic">Entamoeba histolytica</Emphasis>

We report almost complete sequence specific ¹H, ¹³C and ¹⁵N NMR assignments of a 150-residue long calmodulin-like calcium-binding protein from Entamoeba histolytica (EhCaBP6), as a prelude to its structural and functional characterization.     

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.     

<Superscript>1</Superscript>H, <Superscript>13</Superscript>C and <Superscript>15</Superscript>N NMR assignments of an unusual Ca<Superscript>2+</Superscript>-binding protein from <Emphasis Type="Italic">Entamoeba histolytica</Emphasis> in its apo form

We report almost complete sequence specific ¹H, ¹³C and ¹⁵N NMR assignments of an unusual Ca²⁺-binding protein from Entamoeba histolytica (EhCaBP6) in its apo form as a prelude to its structural and functional characterization.     

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.     

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.     

TEDOR with adiabatic inversion pulses: Resonance assignments of <Superscript>13</Superscript>C/<Superscript>15</Superscript>N labelled RNAs

We have examined via numerical simulations the performance characteristics of different 15N RF pulse schemes employed in the transferred echo double resonance (TEDOR) experimental protocol for generating 13C-15N dipolar chemical shift correlation spectra of isotopically labelled biological systems at moderate MAS frequencies (omega(r) approximately 10 kHz). With an 15N field strength of approximately 30-35 kHz that is typically available in 5 mm triple resonance MAS NMR probes, it is shown that a robust TEDOR sequence with significant tolerance to experimental imperfections sa as H1 inhomogeneity and resonance offsets can be effectively implemented using adiabatic heteronuclear dipolar recoupling pulse schemes. TEDOR-based 15N-13C and 15N-13C-13C chemical shift correlation experiments were carried out for obtaining 13C and 15N resonance assignments of an RNA composed of 97 (CUG) repeats which has been implicated in the neuromuscular disease myotonic dystrophy.     

<Superscript>1</Superscript>H, <Superscript>13</Superscript>C and <Superscript>15</Superscript>N backbone resonance assignment of the lamin C-terminal region specific to prelamin A

Lamins are the main components of the nucleoskeleton. They form a protein meshwork that underlies the inner nuclear membrane. Mutations in the LMNA gene coding for A-type lamins (lamins A and C) cause a large panel of human diseases, referred to as laminopathies. These diseases include muscular dystrophies, lipodystrophies and premature aging diseases. Lamin A exhibits a C-terminal region that is different from lamin C and is post-translationally modified. It is produced as prelamin A and it is then farnesylated, cleaved, carboxymethylated and cleaved again in order to become mature lamin A. In patients with the severe Hutchinson–Gilford progeria syndrome, a specific single point mutation in LMNA leads to an aberrant splicing of the LMNA gene preventing the post-translational processing of prelamin A. This leads to the accumulation of a permanently farnesylated lamin A mutant lacking 50 amino acids named progerin. We here report the NMR ¹H, ¹⁵N, ¹³CO, ¹³Cα and ¹³Cβ chemical shift assignment of the C-terminal region that is specific to prelamin A, from amino acid 567 to amino acid 664. We also report the NMR ¹H, ¹⁵N, ¹³CO, ¹³Cα and ¹³Cβ chemical shift assignment of the C-terminal region of the progerin variant, from amino acid 567 to amino acid 614. Analysis of these chemical shift data confirms that both prelamin A and progerin C-terminal domains are largely disordered and identifies a common partially populated α-helix from amino acid 576 to amino acid 585. This helix is well conserved from fishes to mammals.     

High resolution <Superscript>13</Superscript>C-detected solid-state NMR spectroscopy of a deuterated protein

High resolution ¹³C-detected solid-state NMR spectra of the deuterated beta-1 immunoglobulin binding domain of the protein G (GB1) have been collected to show that all ¹⁵N, ¹³C′, ¹³Cα and ¹³Cβ sites are resolved in ¹³C–¹³C and ¹⁵N–¹³C spectra, with significant improvement in T ₂ relaxation times and resolution at high magnetic field (750 MHz). The comparison of echo T ₂ values between deuterated and protonated GB1 at various spinning rates and under different decoupling schemes indicates that ¹³Cα T ₂′ times increase by almost a factor of two upon deuteration at all spinning rates and under moderate decoupling strength, and thus the deuteration enables application of scalar-based correlation experiments that are challenging from the standpoint of transverse relaxation, with moderate proton decoupling. Additionally, deuteration in large proteins is a useful strategy to selectively detect polar residues that are often important for protein function and protein–protein interactions.     

Structure determination of proteins in <Superscript>2</Superscript>H<Subscript>2</Subscript>O solution aided by a deuterium-decoupled 3D HCA(N)CO experiment

We developed an NMR pulse sequence, 3D HCA(N)CO, to correlate the chemical shifts of protein backbone ¹Hα and ¹³Cα to those of ¹³C′ in the preceding residue. By applying ²H decoupling, the experiment was accomplished with high sensitivity comparable to that of HCA(CO)N. When combined with HCACO, HCAN and HCA(CO)N, the HCA(N)CO sequence allows the sequential assignment using backbone ¹³C′ and amide ¹⁵N chemical shifts without resort to backbone amide protons. This assignment strategy was demonstrated for ¹³C/¹⁵N-labeled GB1 dissolved in ²H₂O. The quality of the GB1 structure determined in ²H₂O was similar to that determined in H₂O in spite of significantly smaller number of NOE correlations. Thus this strategy enables the determination of protein structures in ²H₂O or H₂O at high pH values.     

<Superscript>1</Superscript>H, <Superscript>13</Superscript>C, <Superscript>15</Superscript>N resonance assignment of recombinant <Emphasis Type="Italic">Euplotes raikovi</Emphasis> protein Er-23

Er-23 is a small, 51 amino acid, disulfide-rich pheromone protein used for cell signaling by Euplotes raikovi. Ten of the 51 amino acids are cysteine, allowing up to five disulfide bonds. Previous NMR work with Er-23 utilized homologously expressed protein, prohibiting isotopic labeling, and consequently the chemical shift assignments were incomplete. We have expressed uniformly ¹⁵N and ¹³C-labeled Er-23 in an E. coli expression system. Here we report the full backbone and side chain resonance assignments for recombinant Er-23.     

Resonance assignment of 13C/15N labeled solid proteins by two- and three-dimensional magic-angle-spinning NMR

The comprehensive structure determination of isotopically labeled proteins by solid-state NMR requires sequence-specific assignment of ¹³C and ¹⁵ N spectra. We describe several 2D and 3D MAS correlation techniques for resonance assignment and apply them, at 7.0 Tesla, to ¹³C and ¹⁵N labeled ubiquitin to examine the extent of resonance assignments in the solid state. Both interresidue and intraresidue assignments of the ¹³C and ¹⁵N resonances are addressed. The interresidue assignment was carried out by an N(CO)CA technique, which yields N_i-C_i–1 connectivities in protein backbones via two steps of dipolar-mediated coherence transfer. The intraresidue connectivities were obtained from a new 3D NCACB technique, which utilizes the well resolved C chemical shift to distinguish the different amino acids. Additional amino acid type assignment was provided by a ¹³C spin diffusion experiment, which exhibits ¹³C spin pairs as off-diagonal intensities in the 2D spectrum. To better resolve carbons with similar chemical shifts, we also performed a dipolar-mediated INADEQUATE experiment. By cross-referencing these spectra and exploiting the selective and extensive ¹³ C labeling approach, we assigned 25% of the amino acids in ubiquitin sequence-specifically and 47% of the residues to the amino acid types. The sensitivity and resolution of these experiments are evaluated, especially in the context of the selective and extensive ¹³C labeling approach.     

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.     

An intraresidual i(HCA)CO(CA)NH experiment for the assignment of main-chain resonances in <Superscript>15</Superscript>N, <Superscript>13</Superscript>C labeled proteins

An improved pulse sequence, intraresidual i(HCA)CO(CA)NH, is described for establishing solely ¹³C′(i), ¹⁵N(i), ¹H^N(i) connectivities in uniformly ¹⁵N/¹³C-labeled proteins. In comparison to the “out-and-back” style intra-HN(CA)CO experiment, the new pulse sequence offers at least two-fold higher experimental resolution in the ¹³C′ dimension and on average 1.6 times higher sensitivity especially for residues in α-helices. Performance of the new experiment was tested on a small globular protein ubiquitin and an intrinsically unfolded 110-residue cancer/testis antigen CT16/PAGE5. Use of intraresidual i(HCA)CO(CA)NH experiment in combination with the established HNCO experiment was crucial for the assignment of highly disordered CT16.     

Broadband <Superscript>15</Superscript>N–<Superscript>13</Superscript>C dipolar recoupling via symmetry-based RF pulse schemes at high MAS frequencies

An approach for generating efficient RN_n^{n_S, n_k} {\rm{RN}}_{n}^{\nu_{\rm{S}}, {\nu_{\rm{k}}}} symmetry-based dual channel RF pulse schemes for γ-encoded broadband ¹⁵N–¹³C dipolar recoupling at high magic angle spinning frequencies is presented. The method involves the numerical optimisation of the RF phase-modulation profile of the basic “R” element so as to obtain heteronuclear double quantum dipolar recoupling sequences with satisfactory magnetisation transfer characteristics. The basic “R” element was implemented as a sandwich of a small number of short pulses of equal duration with each pulse characterised by a RF phase and amplitude values. The performance characteristics of the sequences were evaluated via numerical simulations and ¹⁵N–¹³C chemical shift correlation experiments. Employing such ¹³C–¹⁵N double-quantum recoupling sequences and the multiple receiver capabilities available in the current generation of NMR spectrometers, the possibility to simultaneously acquire 3D NCC and CNH chemical shift correlation spectra is also demonstrated.