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.     

Determination of protein global folds using backbone residual dipolar coupling and long-range NOE restraints

We report the determination of the global fold of human ubiquitin using protein backbone NMR residual dipolar coupling and long-range nuclear Overhauser effect (NOE) data as conformational restraints. Specifically, by use of a maximum of three backbone residual dipolar couplings per residue (N_i-H^N _i, N_i-C_i–1, H^N _i - C_i–1) in two tensor frames and only backbone H^N-H^N NOEs, a global fold of ubiquitin can be derived with a backbone root-mean-square deviation of 1.4 Å with respect to the crystal structure. This degree of accuracy is more than adequate for use in databases of structural motifs, and suggests a general approach for the determination of protein global folds using conformational restraints derived only from backbone atoms.     

Refinement of d(GCGAAGC) hairpin structure using one- and two-bond residual dipolar couplings

The structure of the ¹³C,¹⁵N-labeled d(GCGAAGC) hairpin, as determined by NMR spectroscopy and refined using molecular dynamics with NOE-derived distances, torsion angles, and residual dipolar couplings (RDCs), is presented. Although the studied molecule is of small size, it is demonstrated that the incorporation of diminutive RDCs can significantly improve local structure determination of regions undefined by the conventional restraints. Very good correlation between the experimental and back-calculated small one- and two-bond ¹H-¹³C, ¹H-¹⁵N, ¹³C-¹³C and ¹³C-¹⁵N coupling constants has been attained. The final structures clearly show typical features of the miniloop architecture. The structure is discussed in context of the extraordinary stability of the d(GCGAAGC) hairpin, which originates from a complex interplay between the aromatic base stacking and hydrogen bonding interactions.     

Determination of residual dipolar couplings in homonuclear MOCCA-SIAM experiments

In solutions with partial molecular alignment, anisotropic magnetic interactions such as the chemical shift anisotropy, the electric quadrupole interaction, and the magnetic dipole-dipole interaction are no longer averaged out to zero in contrast to isotropic solutions. The resulting residual anisotropic magnetic interactions are increasingly used in biological NMR studies for the determination of 3D structures of proteins and other biomolecules. In the present paper we propose a new approach allowing the measurement of residual H^N-H dipolar couplings of non-isotope enriched proteins based on the application of the MOCCA-SIAM experiment. This experiment allows the measurement of homonuclear coupling constants with an accuracy of ca. ±0.2 Hz and is therefore particularly well suited to determine residual dipolar couplings at relatively low degrees of molecular orientation. The agreement between experimentally determined residual H^N-H couplings and calculated values is demonstrated for BPTI.     

Argininamide binding arrests global motions in HIV-1 TAR RNA: comparison with Mg2+-induced conformational stabilization

The structure and dynamics of the stem-loop transactivation response element (TAR) RNA from the human immunodeficiency virus type-1 (HIV-1) bound to the ligand argininamide (ARG) has been characterized using a combination of a large number of residual dipolar couplings (RDCs) and trans-hydrogen bond NMR methodology. Binding of ARG to TAR changes the average inter-helical angle between the two stems from approximately 47 degrees in the free state to approximately 11 degrees in the bound state, and leads to the arrest of large amplitude (+/-46 degrees ) inter-helical motions observed previously in the free state. While the global structural dynamics of TAR-ARG is similar to that previously reported for TAR bound to Mg2+, there are substantial differences in the hydrogen bond alignment of bulge and neighboring residues. Based on a novel H5(C5)NN experiment for probing hydrogen-mediated 2hJ(N,N) scalar couplings as well as measured RDCs, the TAR-ARG complex is stabilized by a U38-A27.U23 base-triple involving an A27.U23 reverse Hoogsteen hydrogen bond alignment as well as by a A22-U40 Watson-Crick base-pair at the junction of stem I. These hydrogen bond alignments are not observed in either the free or Mg2+ bound forms of TAR. The combined conformational analysis of TAR under three states reveals that ligands and divalent ions can stabilize similar RNA global conformations through distinct interactions involving different hydrogen bond alignments in the RNA.     

NMR experiments for the measurement of proton-proton and carbon-carbon residual dipolar couplings in uniformly labelled oligosaccharides

A 2D-HSQC-carbon selective/proton selective-constant time COSY, 2D-HSQC-(sel C, sel H)-CT COSY experiment, which is applicable to uniformly ¹³C isotopically enriched samples (U-¹³C) of oligosaccharides or oligonucleotides is proposed for the measurement of proton–proton RDC in crowded regions of 2D-spectra. In addition, a heteronuclear constant time-COSY experiment, ¹³C-¹³C CT-COSY, is proposed for the measurement of one bond carbon–carbon RDC in these molecules. These two methods provide an extension, to U-¹³C molecules, of the original homonuclear constant time-COSY experiment proposed by Tian et al. (1999) for saccharides. The combination of a number of these RDC with NOE data may provide the method of choice to study oligosaccharide conformation in the free and receptor-bound state.     

Refinement of the structure of protein–RNA complexes by residual dipolar coupling analysis

The main limitation in NMR-determined structures of nucleic acids and their complexes with proteins derives from the elongated, non-globular nature of physiologically important DNA and RNA molecules. Since it is generally not possible to obtain long-range distance constraints between distinct regions of the structure, long-range properties such as bending or kinking at sites of protein recognition cannot be determined accurately nor precisely. Here we show that use of residual dipolar couplings in the refinement of the structure of a protein–RNA complex improves the definition of the long-range properties of the RNA. These features are often an important aspect of molecular recognition and biological function; therefore, their improved definition is of significant value in RNA structural biology.     

Molecular conformations of a disaccharide investigated using NMR spectroscopy

The molecular structure of -l-Rhap-(1→ 2)--l-Rhap-OMe has been investigated using conformation sensitive NMR parameters: cross-relaxation rates, scalar ³ J _CH couplings and residual dipolar couplings obtained in a dilute liquid crystalline phase. The order matrices of the two sugar residues are different, which indicates that the molecule cannot exist in a single conformation. The conformational distribution function, , related to the two glycosidic linkage torsion angles and was constructed using the APME method, valid in the low orientational order limit. The APME approach is based on the additive potential (AP) and maximum entropy (ME) models. The analyses of the trajectories generated in molecular dynamics and Langevin dynamics (LD) computer simulations gave support to the distribution functions constructed from the experimental NMR parameters. It is shown that at least two conformational regions are populated on the Ramachandran map and that these regions exhibit very different molecular order.Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Optimized set of two-dimensional experiments for fast sequential assignment,secondary structure determination,and backbone fold validation of 13C/15N-labelled proteins

NMR experiments are presented which allow backbone resonance assignment, secondary structure identification, and in favorable cases also molecular fold topology determination from a series of two-dimensional ¹H-¹⁵N HSQC-like spectra. The ¹H-¹⁵N correlation peaks are frequency shifted by an amount ± _X along the ¹⁵N dimension, where _X is the C, C, or H frequency of the same or the preceding residue. Because of the low dimensionality (2D) of the experiments, high-resolution spectra are obtained in a short overall experimental time. The whole series of seven experiments can be performed in typically less than one day. This approach significantly reduces experimental time when compared to the standard 3D-based methods. The here presented methodology is thus especially appealing in the context of high-throughput NMR studies of protein structure, dynamics or molecular interfaces.     

Measurement of small scalar and dipolar couplings in purine and pyrimidine bases*

A suite of spin-state-selective excitation (S³E) NMR experiments for the measurements of small one-bond (¹³C-¹³C, ¹⁵N-¹³C) and two-bond (¹H-¹³C, ¹H-¹⁵N) coupling constants in ¹³C,¹⁵N labeled purine and pyrimidine bases is presented. The incorporation of band-selective shaped pulses, elimination of the cross talk between and sub-spectra, and accuracy and precision of the proposed approach are discussed. Merits of using S³E rather than /-half-filter are demonstrated using results obtained on isotopically labeled DNA oligonucleotides.     

Structural characterization of a mannose-binding protein-trimannoside complex using residual dipolar couplings

The ligand-binding properties of a 53 kDa homomultimeric trimer from mannose-binding protein (MBP) have been investigated using residual dipolar couplings (RDCs) that are easily measured from NMR spectra of the ligand and isotopically labeled protein. Using a limited set of 1H-15N backbone amide NMR assignments for MBP and orientational information derived from the RDC measurements in aligned media, an order tensor for MBP has been determined that is consistent with symmetry-based predictions of an axially symmetric system. 13C-1H couplings for a bound trisaccharide ligand, methyl 3,6-di-O-(alpha-D-mannopyranosyl)-alpha-D-mannopyranoside (trimannoside) have been determined at natural abundance and used as orientational constraints. The bound ligand geometry and orientational constraints allowed docking of the trimannoside ligand in the binding site of MBP to produce a structural model for MBP-oligosaccharide interactions.     

NMR investigations of allosteric processes in a two-domain Thermus thermophilus Hsp70 molecular chaperone

Hsp70 chaperones are two-domain proteins that assist in intra-cellular protein (re) folding processes in all species. The protein folding activity of the substrate binding domain of the Hsp70s is regulated by nucleotide binding at the nucleotide-binding domain through an as yet undefined heterotropic allosteric mechanism. The available structures of the isolated domains of Hsp70s have given very limited indications of nucleotide-induced conformational changes that could modulate the affinity for substrate proteins. Here, we present a multi-dimensional NMR study of a prokaryotic Hsp70 homolog, Thermus thermophilus DnaK, using a 54kDa construct containing both nucleotide binding domain and most of the substrate binding domain. It is determined that the nucleotide binding domain and substrate binding domain are closely associated in all ligand states studied. Comparison of the assigned NMR spectra of the two-domain construct with those of the previously studied isolated nucleotide binding domain, allowed the identification of the nucleotide binding domain-substrate binding domain interface. A global three-dimensional structure was obtained for the two-domain construct on the basis of this information and of NMR residual dipolar couplings measurements. This is the first experimental elucidation of the relative positioning of the nucleotide binding domain and substrate binding domain for any Hsp70 chaperone. Comparisons of NMR data between various ligand states including nucleotide-free, ATP, ADP.Pi and ADP.Pi+ peptide bound, identified residues involved in the allosteric inter-domain communication. In particular, peptide binding to the substrate binding domain was found to cause conformational changes in the NBD extending to the nucleotide binding pocket. Detailed analysis suggests that the inter-domain interface becomes tighter in the (nucleotide binding domain ligation/substrate binding domain ligation) order ATP/apo, ADP.Pi/apo ADP.Pi/peptide.     

Accurate determination of leucine and valine side-chain conformations using U-[15N/13C/2H]/[1H-(methine/methyl)-Leu/Val] isotope labeling, NOE pattern recognition, and methine Cgamma-Hgamma/Cbeta-Hbeta residual dipolar couplings: application to the 34-kDa enzyme IIA(chitobiose)

An isotope labeling scheme is described in which specific protonation of methine and methyl protons of leucine and valine is obtained on a ¹⁵N/¹³C labeled background with uniform deuteration of all other non-exchangeable protons. The presence of a protonated methine group has little effect on the favorable relaxation properties of the methyl protons of Leu and Val. This labeling scheme permits the rotameric state of leucine side-chains to be readily determined by simple inspection of the pattern of Hγ(i)–H_N(i) and Hγ(i)–H_N(i+1) NOEs in a 3D ¹⁵N-separated NOE spectrum free of complications arising from spectral overlap and spin-diffusion. In addition, one-bond residual dipolar couplings for the methine ¹³C–¹H bond vectors of Leu and Val can be accurately determined from an intensity J-modulated constant-time HCCH-COSY experiment and used to accurately orient the side-chains of Leu and Val. Incorporation of these data into structure refinement improves the accuracy with which the conformations of Leu and Val side-chains can be established. This is important to ensure optimal packing both within the protein core and at intermolecular interfaces. The impact of the method on protein structure determination is illustrated by application to enzyme IIA^Chitobiose, a 34 kDa homotrimeric phosphotransferase protein.