Characterizing weak protein–protein complexes by NMR residual dipolar couplings

Protein–protein interactions occur with a wide range of affinities from tight complexes characterized by femtomolar dissociation constants to weak, and more transient, complexes of millimolar affinity. Many of the weak and transiently formed protein–protein complexes have escaped characterization due to the difficulties in obtaining experimental parameters that report on the complexes alone without contributions from the unbound, free proteins. Here, we review recent developments for characterizing the structures of weak protein–protein complexes using nuclear magnetic resonance spectroscopy with special emphasis on the utility of residual dipolar couplings.     

Protein backbone motions viewed by intraresidue and sequential HN–Hα residual dipolar couplings

Triple resonance E.COSY-based techniques were used to measure intra-residue and sequential H(N)-H(alpha) residual dipolar couplings (RDCs) for the third IgG-binding domain of protein G (GB3), aligned in Pf1 medium. Measurements closely correlate with values predicted on the basis of an NMR structure, previously determined on the basis of a large number of one-bond backbone RDCs measured in five alignment media. However, in particular the sequential H(N)-H(alpha) RDCs are smaller than predicted for a static structure, suggesting a degree of motion for these internuclear vectors that exceeds that of the backbone amide N-H vectors. Of all experimentally determined GB3 structures available, the best correlation between experimental (1)H-(1)H couplings is observed for a GB3 ensemble, previously derived to generate a realistic picture of the conformational space sampled by GB3 (Clore and Schwieters, J Mol Biol 355:879-886, 2006). However, for both NMR and X-ray-derived structures the (1)H-(1)H couplings are found to be systematically smaller than expected on the basis of alignment tensors derived from (15)N-(1)H amide RDCs, assuming librationally corrected N-H bond lengths of 1.041 A.     

Measurement of 1H�C15N and 1H�C13C residual dipolar couplings in nucleic acids from TROSY intensities

Analogous to the recently introduced ARTSY method for measurement of one-bond (1)H-(15)N residual dipolar couplings (RDCs) in large perdeuterated proteins, we introduce methods for measurement of base (13)C-(1)H and (15)N-(1)H RDCs in protonated nucleic acids. Measurements are based on quantitative analysis of intensities in (1)H-(15)N and (13)C-(1)H TROSY-HSQC spectra, and are illustrated for a 71-nucleotide adenine riboswitch. Results compare favorably with those of conventional frequency-based measurements in terms of completeness and convenience of use. The ARTSY method derives the size of the coupling from the ratio of intensities observed in two TROSY-HSQC spectra recorded with different dephasing delays, thereby minimizing potential resonance overlap problems. Precision of the RDC measurements is limited by the signal-to-noise ratio, S/N, achievable in the 2D TROSY-HSQC reference spectrum, and is approximately given by 30/(S/N) Hz for (15)N-(1)H and 65/(S/N) Hz for (13)C-(1)H. The signal-to-noise ratio of both (1)H-(15)N and (1)H-(13)C spectra greatly benefits when water magnetization during the experiments is not perturbed, such that rapid magnetization transfer from bulk water to the nucleic acid, mediated by rapid amino and hydroxyl hydrogen exchange coupled with (1)H-(1)H NOE transfer, allows for fast repetition of the experiment. RDCs in the mutated helix 1 of the riboswitch are compatible with nucleotide-specifically modeled, idealized A-form geometry and a static orientation relative to the helix 2/3 pair, which differs by ca 6° relative to the X-ray structure of the native riboswitch.     

Measurement of residual dipolar couplings from 1Hα to 13Cα and 15N using a simple HNCA-based experiment

Novel NMR pulse schemes for simultaneous measurement of ¹ D _CHand ² D _NHresidual dipolar couplings in proteins is presented. We show that ² D _NHcoupling can be very useful for protein structure determination. The ² D _NHcoupling can be measured from ¹⁵N dimension with good accuracy on a slowly relaxing TROSY resonance, utilizing HNCA-TROSY-based experiments, which concomitantly supply large ¹ D _CHcoupling. The dynamic range of ² D _NHcoupling is comparable to ¹ D _NC coupling, but instead, it also serves non-redundant information on the course of protein backbone, thanks to rotational degree of freedom with respect to peptide bond. The HNCA-TROSY-based experiments are optimal for measuring residual dipolar couplings at high magnetic fields owing to absence of rapid transverse relaxation of carbonyl carbon. The reliability of the proposed approach was tested on ¹⁵N/¹³C human ubiquitin. A very good correlation with ubiquitin solution as well as crystal structure, for both ¹ D _CHand ² D _NHcouplings, was obtained.     

Improved detection of long-range residual dipolar couplings in weakly aligned samples by Lee–Goldburg decoupling of homonuclear dipolar truncation

Homonuclear ¹H residual dipolar couplings (RDCs) truncate the evolution of transverse ¹H magnetization of weakly aligned molecules in high-resolution NMR experiments. This leads to losses in sensitivity or resolution in experiments that require extended ¹H evolution times. Lee–Goldburg decoupling schemes have been shown to remove the effects of homonuclear dipolar couplings, while preserving chemical shift evolution in a number of solid-state NMR applications. Here, it is shown that the Lee-Goldburg sequence can be effectively incorporated into INEPT- or HMQC-type transfer schemes in liquid state weak alignment experiments in order to increase the efficiency of the magnetization transfer. The method is applied to the sensitive detection of ¹H^N–¹³C long-range RDCs in a three-dimensional HCN experiment. As compared to a conventional HCN experiment, an average sensitivity increase by a factor of 2.4 is obtained for a sample of weakly aligned protein G. This makes it possible to detect 170 long-range ¹H^N–¹³C RDCs for distances up to 4.9 Å     

Improvement of hydrogen bond geometry in protein NMR structures by residual dipolar couplings – an assessment of the interrelation of NMR restraints

We have examined how the hydrogen bond geometry in three different proteins is affected when structural restraints based on measurements of residual dipolar couplings are included in the structure calculations. The study shows, that including restraints based solely on (1)H(N)-(15)N residual dipolar couplings has pronounced impact on the backbone rmsd and Ramachandran plot but does not improve the hydrogen bond geometry. In the case of chymotrypsin inhibitor 2 the addition of (13)CO-(13)C(alpha) and (15)N-(13)CO one bond dipolar couplings as restraints in the structure calculations improved the hydrogen bond geometry to a quality comparable to that obtained in the 1.8 A resolution X-ray structure of this protein. A systematic restraint study was performed, in which four types of restraints, residual dipolar couplings, hydrogen bonds, TALOS angles and NOEs, were allowed in two states. This study revealed the importance of using several types of residual dipolar couplings to get good hydrogen bond geometry. The study also showed that using a small set of NOEs derived only from the amide protons, together with a full set of residual dipolar couplings resulted in structures of very high quality. When reducing the NOE set, it is mainly the side-chain to side-chain NOEs that are removed. Despite of this the effect on the side-chain packing is very small when a reduced NOE set is used, which implies that the over all fold of a protein structure is mainly determined by correct folding of the backbone.     

Simultaneous measurement of 1H�C15N and Methyl 1Hm�C13Cm residual dipolar couplings in large proteins

A two-dimensional TROSY-based SIM-(13)C(m)-(1)H(m)/(1)H-(15)N NMR experiment for simultaneous measurements of methyl (1) D (CH) and backbone amide (1) D (NH) residual dipolar couplings (RDC) in {U-[(15)N,(2)H]; Ileδ1-[(13)CH(3)]; Leu,Val-[(13)CH(3)/(12)CD(3)]}-labeled samples of large proteins is described. Significant variation in the alignment tensor of the 82-kDa enzyme Malate synthase G is observed as a function of only slight changes in experimental conditions. The SIM-(13)C(m)-(1)H(m)/(1)H-(15)N data sets provide convenient means of establishing the alignment tensor characteristics via the measurement of (1) D (NH) RDCs in the same protein sample.     

Three-dimensional structure of α-conotoxin EI determined by1H NMR spectroscopy

Summary α-conotoxin EI is an 18-residue peptide (RDOCCYHPTCNMSNPQIC; 4–10, 5–18) isolated from the venom ofConus ermineus, the only fish-hunting cone snail of the Atlantic Ocean. This peptide targets specifically the nicotinic acetylcholine receptor (nAChR) found in mammalian skeletal muscle and the electric organTorpedo, showing a novel selectivity profile when compared to other α-conotoxins. The 3D structure of EI has been determined by 2D-NMR methods in combination with dynamical simulated annealing protocols. A total of 133 NOE-derived distances were used to produce 13 structures with minimum energy that complied with the NOE restraints. The structure of EI is characterized by a helical loop between THr⁹ and Met¹² that is stabilized by the Cys⁴-Cys¹⁰ disulfide bond and turns involving Cys⁴-Cys⁵ and Asn¹⁴-Pro¹⁵. Other regions of the peptide appear to be flexible. The overall fold of EI is similar to that of other α4/7-conotoxins (PnIA/B, MII, EpI). However, unlike these other α4/7-conotoxins, EI targets the muscular type nAChR. The differences in selectivity can be attributed to differences in the surface charge distribution among these α4/7-conotoxins. The implications for binding of EI to the muscular nAChR are discussed with respect to the current NMR structure of EI. Supplementary material available:¹H resonance assignments of α-conotoxin EI.     

Three-dimensional structure of α-conotoxin EI determined by 1H NMR spectroscopy H resonance assignments of α-conotoxin EI

α-Conotoxin EI is an 18-residue peptide (RDOCCYHPTCNMSNPQIC; 4–10, 5–18) isolated from the venom of Conus ermineus, the only fish-hunting cone snail of the Atlantic Ocean. This peptide targets specifically the nicotinic acetylcholine receptor (nAChR) found in mammalian skeletal muscle and the electric organ Torpedo, showing a novel selectivity profile when compared to other α-conotoxins. The 3D structure of EI has been determined by 2D-NMR methods in combination with dynamical simulated annealing protocols. A total of 133 NOE-derived distances were used to produce 13 structures with minimum energy that complied with the NOE restraints. The structure of EI is characterized by a helical loop between Thr⁹ and Met¹² that is stabilized by the Cys⁴-Cys¹⁰ disulfide bond and turns involving Cys⁴-Cys⁵ and Asn¹⁴-Pro¹⁵. Other regions of the peptide appear to be flexible. The overall fold of EI is similar to that of other α4/7-conotoxins (PnIA/B, MII, EpI). However, unlike these other α4/7-conotoxins, EI targets the muscular type nAChR. The differences in selectivity can be attributed to differences in the surface charge distribution among these α4/7-conotoxins. The implications for binding of EI to the muscular nAChR are discussed with respect to the current NMR structure of EI.     

Facile measurement of <Superscript>1</Superscript>H–<Superscript>15</Superscript>N residual dipolar couplings in larger perdeuterated proteins

We present a simple method, ARTSY, for extracting ¹J_NH couplings and ¹H–¹⁵N RDCs from an interleaved set of two-dimensional ¹H–¹⁵N TROSY-HSQC spectra, based on the principle of quantitative J correlation. The primary advantage of the ARTSY method over other methods is the ability to measure couplings without scaling peak positions or altering the narrow line widths characteristic of TROSY spectra. Accuracy of the method is demonstrated for the model system GB3. Application to the catalytic core domain of HIV integrase, a 36 kDa homodimer with unfavorable spectral characteristics, demonstrates its practical utility. Precision of the RDC measurement is limited by the signal-to-noise ratio, S/N, achievable in the 2D TROSY-HSQC spectrum, and is approximately given by 30/(S/N) Hz.     

Measurement of imino <Superscript>1</Superscript>H–<Superscript>1</Superscript>H residual dipolar couplings in RNA

Imino ¹H–¹⁵N residual dipolar couplings (RDCs) provide additional structural information that complements standard ¹H–¹H NOEs leading to improvements in both the local and global structure of RNAs. Here, we report measurement of imino ¹H–¹H RDCs for the Iron Responsive Element (IRE) RNA and native E. coli tRNA^Val using a BEST-Jcomp-HMQC2 experiment. ¹H–¹H RDCs are observed between the imino protons in G–U wobble base pairs and between imino protons on neighboring base pairs in both RNAs. These imino ¹H–¹H RDCs complement standard ¹H–¹⁵N RDCs because the ¹H–¹H vectors generally point along the helical axis, roughly perpendicular to ¹H–¹⁵N RDCs. The use of longitudinal relaxation enhancement increased the signal-to-noise of the spectra by ~3.5-fold over the standard experiment. The ability to measure imino ¹H–¹H RDCs offers a new restraint, which can be used in NMR domain orientation and structural studies of RNAs. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Three-dimensional structure of α-conotoxin EI determined by 1H NMR spectroscopy H resonance assignments of α-conotoxin EI

-Conotoxin EI is an 18-residue peptide (RDOCCYHPTCNMSNPQIC; 4–10, 5–18) isolated from the venom of Conus ermineus, the only fish-hunting cone snail of the Atlantic Ocean. This peptide targets specifically the nicotinic acetylcholine receptor (nAChR) found in mammalian skeletal muscle and the electric organ Torpedo, showing a novel selectivity profile when compared to other -conotoxins. The 3D structure of EI has been determined by 2D-NMR methods in combination with dynamical simulated annealing protocols. A total of 133 NOE-derived distances were used to produce 13 structures with minimum energy that complied with the NOE restraints. The structure of EI is characterized by a helical loop between Thr⁹ and Met¹² that is stabilized by the Cys⁴-Cys¹⁰ disulfide bond and turns involving Cys⁴-Cys⁵ and Asn¹⁴-Pro¹⁵. Other regions of the peptide appear to be flexible. The overall fold of EI is similar to that of other 4/7-conotoxins (PnIA/B, MII, EpI). However, unlike these other 4/7-conotoxins, EI targets the muscular type nAChR. The differences in selectivity can be attributed to differences in the surface charge distribution among these 4/7-conotoxins. The implications for binding of EI to the muscular nAChR are discussed with respect to the current NMR structure of EI.     

Improved accuracy of <Superscript>15</Superscript>N–<Superscript>1</Superscript>H scalar and residual dipolar couplings from gradient-enhanced IPAP-HSQC experiments on protonated proteins

The presence of dipole-dipole cross-correlated relaxation as well as unresolved E.COSY effects adversely impacts the accuracy of ¹ J _NH splittings measured from gradient-enhanced IPAP-HSQC spectra. For isotropic samples, the size of the systematic errors caused by these effects depends on the values of ² J _NHα , ³ J _NHβ and ³ J _HNHα . Insertion of band-selective ¹H decoupling pulses in the IPAP-HSQC experiment eliminates these systematic errors and for the protein GB3 yields ¹ J _NH splittings that agree to within a root-mean-square difference of 0.04 Hz with values measured for perdeuterated GB3. Accuracy of the method is also highlighted by a good fit to the GB3 structure of the ¹H-¹⁵N RDCs extracted from the minute differences in ¹J_NH splitting measured at 500 and 750 MHz ¹H frequencies, resulting from magnetic susceptibility anisotropy. A nearly complete set of ² J _NHα couplings was measured in GB3 in order to evaluate whether the impact of cross-correlated relaxation is dominated by the ¹⁵N–¹H ^α or ¹⁵N–¹H ^β dipolar interaction. As expected, we find that ² J _NHα  ≤ 2 Hz, with values in the α-helix (0.86 ± 0.52 Hz) slightly larger than in β-sheet (0.66 ± 0.26 Hz). Results indicate that under isotropic conditions, N–H^N/N–H ^β cross-correlated relaxation often dominates. Unresolved E.COSY effects under isotropic conditions involve ³ J _HNHα and J _NHα , but when weakly aligned any aliphatic proton proximate to both N and H^N can contribute. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Pressure effect on the dynamics of an isolated α-helix studied by 15N-1H NMR relaxation

Dynamics and structure of (1–36)bacteriorhodopsin solubilized in chloroform/methanol mixture (1:1) were investigated by ¹H-¹⁵N NMR spectroscopy under a hydrostatic pressure of 2000 bar. It was shown that the peptide retains its spatial structure at high pressure. ¹⁵N transverse and longitudinal relaxation times, ¹⁵N{¹H} nuclear Overhauser effects, chemical shifts and the translation diffusion rate of the peptide at 2000 bar were compared with the respective data at ambient pressure [Orekhov et al. (1999) J. Biomol. NMR, 14, 345–356]. The model free analysis of the relaxation data for the helical 9–31 fragment revealed that the high pressure decreases the overall rotation and translation diffusion, as well as apparent order parameters of fast picosecond internal motions (S² _f) but has no effect on internal nanosecond motions (S² _s and _s) of the peptide. The decrease of translation and overall rotation diffusion was attributed to the increase in solvent viscosity and the decrease of apparent order parameters S² _f to a compression of hydrogen bonds. It is suggested that this compression causes an elongation of H-N bonds and a decrease of absolute values of chemical shift anisotropy (CSA). In particular, the observed decrease of S² _f at 2000 bar can be explained by 0.001 nm increase of N-H bond lengths and 10 ppm decrease of ¹⁵N CSA values.     

Probing the equilibrium unfolding of ketosteroid isomerase through xenon-perturbed 1H–15N multidimensional NMR spectroscopy

We used xenon-perturbed ¹H–¹⁵N multidimensional NMR to investigate the structural changes in the urea-induced equilibrium unfolding of the dimeric ketosteroid isomerase (KSI) from Pseudomonas putida biotype B. Three limited regions located on the β3-, β5- and β6-strands of dimeric interface were significantly perturbed by urea in the early stage of KSI unfolding, which could lead to dissociation of the dimer into structured monomers at higher denaturant concentration as the interactions in these regions are weakened. The results indicate that the use of xenon as an indirect probe for multidimensional NMR can be a useful method for the equilibrium unfolding study of protein at residue level.     

1H, 13C and 15N NMR assignments and solution secondary structure of rat Apo-S100β

Summary The ¹H, ¹³C and ¹⁵N NMR assignments of the backbone and side-chain resonances of rat S100 were made at pH 6.5 and 37°C using heteronuclear multidimensional NMR spectroscopy. Analysis of the NOE correlations, together with amide exchange rate and ¹H, ¹³C and ¹³C chemical shift data, provided extensive secondary structural information. Thus, the secondary structure of S100 was determined to comprise four helices (Leu³-Ser¹⁸, helix I; Lys²⁹-Leu⁴⁰, helix II; Gln⁵⁰-Glu⁶², helix III; and Phe⁷⁰-Ala⁸³, helix IV), four loops (Gly¹⁹-His²⁵, loop I; Ser⁴¹-Glu⁴⁹, loop II; Asp⁶³-Gly⁶⁶, loop III; and Cys⁸⁴-Glu⁹¹, loop IV) and two -strands (Lys²⁶-Lys²⁸, -strand I and Glu⁶⁷-Asp⁶⁹, -strand II). The -strands were found to align in an antiparallel manner to form a very small -sheet. This secondary structure is consistent with predictions that S100 contains two helix-loop-helix Ca²⁺-binding motifs known as EF-hands. The alignment of the -sheet, which brings the two EF-hand domains of S100 into close proximity, is similar to that of several other Ca²⁺ ion-binding proteins.     

Measurement of 1H3′-31P dipolar couplings in a DNA oligonucleotide by constant-time NOESY difference spectroscopy

The ratios of cross peak intensities in a selective constant-time NOESY experiment, recorded with and without ³¹P decoupling, yield values for the sum of the H3-P scalar and dipolar couplings. The selective refocusing of H3 resonances in this experiment results in excellent resolution and sensitivity, even in the liquid crystalline phase where the ¹H spectrum is broadened by unresolved homonuclear dipolar couplings. The vicinal H3-P scalar and dipolar couplings in the DNA oligomer d(CGCGAATTCGCG)₂ were measured in both isotropic solution, and in a liquid crystalline phase. Isotropic values are in good agreement with values reported previously. Dipolar couplings are in excellent agreement with the NMR structure for this dodecamer, and to a somewhat lesser extent with the X-ray structures.     

Accurate measurement of <Superscript>15</Superscript>N–<Superscript>13</Superscript>C residual dipolar couplings in nucleic acids

New 3D HCN quantitative J (QJ) pulse schemes are presented for the precise and accurate measurement of one-bond ¹⁵N_1/9–¹³C₁, ¹⁵N_1/9–¹³C_6/8, and ¹⁵N_1/9–¹³C_2/4 residual dipolar couplings (RDCs) in weakly aligned nucleic acids. The methods employ ¹H–¹³C multiple quantum (MQ) coherence or TROSY-type pulse sequences for optimal resolution and sensitivity. RDCs are obtained from the intensity ratio of H₁–C₁–N_1/9 (MQ-HCN-QJ) or H_6/8–C_6/8–N_1/9 (TROSY-HCN-QJ) correlations in two interleaved 3D NMR spectra, with dephasing intervals of zero (reference spectrum) and 1/(2J_NC) (attenuated spectrum). The different types of ¹⁵N–¹³C couplings can be obtained by using either the 3D MQ-HCN-QJ or TROSY-HCN-QJ pulse scheme, with the appropriate setting of the duration of the constant-time ¹⁵N evolution period and the offset of two frequency-selective ¹³C pulses. The methods are demonstrated for a uniformly ¹³C, ¹⁵N-enriched 24-nucleotide stem-loop RNA sequence, helix-35, aligned in the magnetic field using phage Pf1. For measurements of RDCs systematic errors are found to be negligible, and experiments performed on a 1.5 mM helix-35 sample result in an estimated precision of ca. 0.07 Hz for ¹D_NC, indicating the utility of the measured RDCs in structure validation and refinement. Indeed, for a complete set of ¹⁵N_1/9–¹³C₁, ¹⁵N_1/9–¹³C_6/8, and ¹⁵N_1/9–¹³C_2/4 dipolar couplings obtained for the stem nucleotides, the measured RDCs are in excellent agreement with those predicted for an NMR structure of helix-35, refined using independently measured observables, including ¹³C–¹H, ¹³C–¹³C and ¹H–¹H dipolar couplings.Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1007/s10858-005-0646-2.     

Measurement and application of 1H-19F dipolar couplings in the structure determination of 2′-fluorolabeled RNA

Residual dipolar couplings can provide powerful restraints for determination and refinement of the solution structure of macromolecules. The application of these couplings in nucleic acid structure elucidation can have an especially dramatic impact, since they provide long-range restraints, typically absent in NOE and J-coupling measurements. Here we describe sensitive X-filtered-E.COSY-type methods designed to measure both the sign and magnitude of long-range ¹H-¹⁹F dipolar couplings in selectively fluorine labeled RNA oligonucleotides oriented in solution by a liquid crystalline medium. The techniques for measuring ¹H-¹⁹F dipolar couplings are demonstrated on a 21-mer RNA hairpin, which has been specifically labeled with fluorine at the 2-hydroxyl position of three ribose sugars. Experimentally measured ¹H-¹⁹F dipolar couplings for the 2-deoxy-2-fluoro-sugars located in the helical region of the RNA hairpin were found to be in excellent agreement with values predicted using canonical A-form helical geometry, demonstrating that these couplings can provide accurate restraints for the refinement of RNA structures determined by NMR.     

[15N,1H]/[13C,1H]-TROSY for simultaneous detection of backbone 15N–1H, aromatic 13C–1H and side-chain 15N–1H2 correlations in large proteins

This paper describes a [¹⁵N,¹H]/[¹³C,¹H]-TROSY experiment for the simultaneous acquisition of the heteronuclear chemical shift correlations of backbone amide ¹⁵N–¹H groups, side chain ¹⁵N–¹H₂ groups and aromatic ¹³C–¹H groups in otherwise highly deuterated proteins. The ¹⁵N–¹H and ¹³C–¹H correlations are extracted from two subspectra of the same data set, thus preventing possible spectral overlap of aromatic and amide protons in the ¹H dimension. The side-chain ¹⁵N–¹H₂ groups, which are suppressed in conventional [¹⁵N,¹H)-TROSY, are observed with high sensitivity in the ¹⁵N–¹H subspectrum. [¹⁵N,¹H]/[¹³C,¹H]-TROSY was used as the heteronuclear correlation block in a 3D [¹H,¹H]-NOESY-[¹⁵N,¹H]/[¹³C,¹H]-TROSY experiment with the membrane protein OmpA reconstituted in detergent micelles of molecular weight 80000 Da, which enabled the detection of numerous NOEs between backbone amide protons and both aromatic protons and side chain ¹⁵N–¹H₂ groups.