Structures of protein-protein complexes are docked using only NMR restraints from residual dipolar coupling and chemical shift perturbations

NMR structures of protein-protein and protein-ligand complexes rely heavily on intermolecular NOEs. Recent work has shown that if no significant conformational changes occur upon complex formation residual dipolar coupling can replace most of the NOE restraints in protein-protein complexes, while restraints derived from chemical shift perturbations can largely replace intermolecular NOEs in protein-ligand structures. By combining restraints from chemical shift perturbations with orientation restraints derived from measurements of residual dipolar couplings, we show that the structure of the EIN-HPr complex can be calculated without NOE restraints. The final structure, built from the crystal structures of EIN and HPr in their uncomplexed form and docked only with NMR restraints, places HPr within 2.5 A of the position determined from the mean NMR structure of the complex.     

15N,(1)H Heteronuclear correlation NMR of gramicidin A in DMPC-d(67)

We have studied gramicidin A, an environmentally sensitive polymorphic pentadecapeptide, fully 15N-labelled and dispersed in a highly deuterated phospholipid bilayer system. By submitting the sample to fast magic angle spinning, we were able to reduce the polypeptide amide hydrogen linewidths to 160 Hz, and hence to partially resolve them. By correlating these resonances with the 40 Hz wide dipolar coupled 15N in a 2D-CROPSY (cross-polarization spectroscopy) experiment, it was possible to observe the 20 partially overlapping 1H-15N signal pairs from the polypeptide backbone and sidechains. Both chemical shift distributions closely match those of the same peptide in SDS micelles, but only poorly match those of conformationally different gramicidin A in trifluoroethanol, dimethylsulfoxide, or methanol/chloroform mixture. Our results are indicative of the N-to-N right-handed beta6.3-helix conformation of gramicidin A and offer sufficient resolution to encourage development of experiments to measure orientational or distance restraints using through-space dipolar couplings.     

Determination of natural abundance 15N-1H and 13C-1H dipolar couplings of molecules in a strongly orienting media using two-dimensional inverse experiments

NMR spectra of molecules oriented in liquid crystals provide homo- and heteronuclear dipolar couplings and thereby the geometry of the molecules. Several inequivalent dilute spins such as 13C and 15N coupled to protons form different coupled spin systems in their natural abundance and appear as satellites in the proton spectra. Identification of transitions belonging to each spin system is essential to determine heteronuclear dipolar couplings, which is a formidable task. In the present study, using 15N-1H and 13C-1H HSQC, and HMQC experiments we have selectively detected spectra of each rare spin coupled to protons. The 15N-1H and 13C-1H dipolar couplings have been determined in the natural abundance of 13C and 15N for the molecules pyrazine, pyrimidine and pyridazine oriented in a thermotropic liquid crystal.     

Determination of multiple torsion-angle constraints in U-(13)C,(15)N-labeled peptides: 3D (1)H-(15)N-(13)C-(1)H dipolar chemical shift NMR spectroscopy in rotating solids

We demonstrate constraint of peptide backbone and side-chain conformation with 3D (1)H-(15)N-(13)C-(1)H dipolar chemical shift, magic-angle spinning NMR experiments. In these experiments, polarization is transferred from (15)N[i] by ramped SPECIFIC cross polarization to the (13)C(alpha)[i], (13)C(beta)[i], and (13)C(alpha)[i - 1] resonances and evolves coherently under the correlated (1)H-(15)N and (1)H-(13)C dipolar couplings. The resulting set of frequency-labeled (15)N(1)H-(13)C(1)H dipolar spectra depend strongly upon the molecular torsion angles phi[i], chi1[i], and psi[i - 1]. To interpret the data with high precision, we considered the effects of weakly coupled protons and differential relaxation of proton coherences via an average Liouvillian theory formalism for multispin clusters and employed average Hamiltonian theory to describe the transfer of (15)N polarization to three coupled (13)C spins ((13)C(alpha)[i], (13)C(beta)[i], and (13)C(alpha)[i - 1]). Degeneracies in the conformational solution space were minimized by combining data from multiple (15)N(1)H-(13)C(1)H line shapes and analogous data from other 3D (1)H-(13)C(alpha)-(13)C(beta)-(1)H (chi1), (15)N-(13)C(alpha)-(13)C'-(15)N (psi), and (1)H-(15)N[i]-(15)N[i + 1]-(1)H (phi, psi) experiments. The method is demonstrated here with studies of the uniformly (13)C,(15)N-labeled solid tripeptide N-formyl-Met-Leu-Phe-OH, where the combined data constrains a total of eight torsion angles (three phi, three chi1, and two psi): phi(Met) = -146 degrees, psi(Met) = 159 degrees, chi1(Met) = -85 degrees, phi(Leu) = -90 degrees, psi(Leu) = -40 degrees, chi1(Leu) = -59 degrees, phi(Phe) = -166 degrees, and chi1(Phe) = 56 degrees. The high sensitivity and dynamic range of the 3D experiments and the data analysis methods provided here will permit immediate application to larger peptides and proteins when sufficient resolution is available in the (15)N-(13)C chemical shift correlation spectra.     

Measurement of multiple psi torsion angles in uniformly 13C,15N-labeled alpha-spectrin SH3 domain using 3D 15N-13C-13C-15N MAS dipolar-chemical shift correlation spectroscopy

We demonstrate the simultaneous measurement of several backbone torsion angles psi in the uniformly (13)C,(15)N-labeled alpha-Spectrin SH3 domain using two different 3D 15N-13C-13C-15N dipolar-chemical shift magic-angle spinning (MAS) NMR experiments. The first NCCN experiment utilizes double quantum (DQ) spectroscopy combined with the INADEQUATE type 13C-13C chemical shift correlation. The decay of the DQ coherences formed between 13C'(i) and 13C(alphai) spin pairs is determined by the "correlated" dipolar field due to 15N(i)-13C(alphai) and 13C'(i)-15N(i+1) dipolar couplings and is particularly sensitive to variations of the torsion angle in the regime |psi| > 140 degrees. However, the ability of this experiment to constrain multiple psi-torsion angles is limited by the resolution of the 13C(alpha)-(13)CO correlation spectrum. This problem is partially addressed in the second approach described here, which is an NCOCA NCCN experiment. In this case the resolution is enhanced by the superior spectral dispersion of the 15N resonances present in the 15N(i+1)-13C(alphai) part of the NCOCA chemical shift correlation spectrum. For the case of the 62-residue alpha-spectrin SH3 domain, we determined 13 psi angle constraints with the INADEQUATE NCCN experiment and 22 psi constraints were measured in the NCOCA NCCN experiment.     

1H, 13C and 15N NMR spectral characterization of twenty-seven 1,2-diaryl-(4E)-arylidene-2-imidazolin-5-ones

1H, 13C and 15N NMR chemical shifts and couplings (n)J(H,C) in DMSO-d(6) at 30 degrees C have been determined for 1,2-diaryl-(4E)-arylidene-2-imidazolin-5-one derivatives 1-27. Their chemical shift assignments are based on PFG DQF 1H,1H COSY, PFG 1H,13C HMQC as well as PFG 1H,13C and 1H,15N HMBC experiments. For compounds 1-10 including aryl fluorine substituent(s) also the couplings (n)J(F,C) (n = 1 - 4) are reported.     

NMR “Crystallography” for Uniformly (13C, 15N)-Labeled Oriented Membrane Proteins

In oriented-sample (OS) solid-state NMR of membrane proteins, the angular-dependent dipolar couplings and chemical shifts provide a direct input for structure calculations. However, so far only ¹H–¹⁵N dipolar couplings and ¹⁵N chemical shifts have been routinely assessed in oriented ¹⁵N-labeled samples. The main obstacle for extending this technique to membrane proteins of arbitrary topology has remained in the lack of additional experimental restraints. We have developed a new experimental triple-resonance NMR technique, which was applied to uniformly doubly (¹⁵N, ¹³C)-labeled Pf1 coat protein in magnetically aligned DMPC/DHPC bicelles. The previously inaccessible ¹H_α–¹³C_α dipolar couplings have been measured, which make it possible to determine the torsion angles between the peptide planes without assuming α-helical structure a priori. The fitting of three angular restraints per peptide plane and filtering by Rosetta scoring functions has yielded a consensus α-helical transmembrane structure for Pf1 protein.     

NMR “Crystallography” for Uniformly (13C, 15N)‐Labeled Oriented Membrane Proteins

In oriented‐sample (OS) solid‐state NMR of membrane proteins, the angular‐dependent dipolar couplings and chemical shifts provide a direct input for structure calculations. However, so far only ¹H–¹⁵N dipolar couplings and ¹⁵N chemical shifts have been routinely assessed in oriented ¹⁵N‐labeled samples. The main obstacle for extending this technique to membrane proteins of arbitrary topology has remained in the lack of additional experimental restraints. We have developed a new experimental triple‐resonance NMR technique, which was applied to uniformly doubly (¹⁵N, ¹³C)‐labeled Pf1 coat protein in magnetically aligned DMPC/DHPC bicelles. The previously inaccessible ¹H_α–¹³C_α dipolar couplings have been measured, which make it possible to determine the torsion angles between the peptide planes without assuming α‐helical structure a priori. The fitting of three angular restraints per peptide plane and filtering by Rosetta scoring functions has yielded a consensus α‐helical transmembrane structure for Pf1 protein.     

Experimental and theoretical investigation of the 13C and 15N chemical shift tensors in melanostatin-exploring the chemical shift tensor as a structural probe

The determination of backbone conformations in powdered peptides using 13C and 15N shift tensor information is explored. The 13C and 15N principal shift values in natural abundance 13C and 15N melanostatin (L-Pro-L-Leu-Gly amide) are measured using the FIREMAT technique. Furthermore, the orientation of the C-N bond in the 13C shift principal axis system for the backbone carbons is obtained from the presence of the 13C-14N dipolar coupling. The Ramachandran angles for the title compound are obtained from solid-state NMR data by comparing the experimentally determined shift tensor information to systematic theoretical shielding calculations on N-formyl-L-amino acid-amide models. The effects of geometry optimization and neglect of intermolecular interactions on the theoretical shielding values in the model compounds are investigated. The sets of NMR derived Ramachandran angles are assembled in a set of test structures that are compared to the available single-crystal X-ray structure. Shift tensor calculations on the test structures and the X-ray structure are used to further assess the importance of intermolecular interactions when the shift tensor is used as a structural probe in powdered peptides.     

Double quantum 1H MAS NMR studies of hydrogen-bonded protons and water dynamics in materials

Two-dimensional double quantum (DQ) 1H MAS NMR was used to investigate different proton environments in a series of alkali (Na, K, Rb, Cs) [Nb6O19]8- Lindqvist salts, with the water and hydrogen-bound intercluster protons being clearly resolved. Through the analysis of the DQ 1H NMR spinning sideband pattern, it is possible to extract both the mean and distribution of the motionally averaged intramolecular homonuclear 1H-1H dipolar coupling for the different water environments and the intercluster protons. Motional order parameters for the water environments were then calculated from the averaged dipolar couplings. The influence of additional intermolecular dipolar couplings due to multispin interactions were simulated and discussed.     

Protein Backbone 1H(N)-13Calpha and 15N-13Calpha residual dipolar and J couplings: new constraints for NMR structure determination

A simple, sensitivity-enhanced experiment was devised for accurate measurement of backbone 15N-13Calpha and 1HN-13Calpha couplings in proteins. The measured residual dipolar couplings 2DHCA, 1DNCA, 3DHCA, and 2DNCA for protein GB1 display very good agreement with the refined NMR structure (PDB code: 3GB1). A Karplus-type relationship between the one-bond 1JNCA couplings and the backbone dihedral psi angles holds, and on the basis of the two-bond 2JNCA couplings a secondary structure index can be established.     

13C/15N distance determination by CPMAS NMR in uniformly 13C labeled molecules

The REDOR and CPMAS techniques are applied for measuring 13C-15N dipolar coupling constants in glycine. It is shown that the selective CP or SPECIFIC CP technique removes the coherent evolution of the spin system under homonuclear 13C-13C J couplings. While the large coupling constant (approximately 900 Hz) is readily determined because of the presence of large oscillations in the CPMAS dynamics, their absence precludes the measurement of the small coupling constant (approximately 200 Hz). The experimental results and numerical simulations demonstrate that the determination of 13C-15N coupling constants of medium size (<1 kHz) by the CPMAS technique is mainly limited by the strength of the 1H decoupling field and the size of the 13C and 15N chemical shift anisotropies.     

应用相敏HMQC和HSQC谱准确测定异核碳-氢偶合常数

准确测定各种同核和异核偶合常数是核磁共振（NMR）方法研究的一个非常跃的领域。首先,各种三键偶合常数通过Karplus关系式^[1]反映了相应二面角的大小,因此,多键偶合常数的准确测定直接影响分子结构确定的精确性。其次,由于稀液晶溶剂体系NMR方法的发展^[2],准确测定各种异核键偶合常数也显得非常重要,特别是应用场相关偶合常数研究分子在磁场中的取向时,对异核－键偶合常数测定的准确性要求更加严格^[3]。异核－键偶合常数的最准确的测定方法是异核偶合调制的HSQC（Heteronuclear Single-Quantum Coherence)实验^[3],它通过测定一系列异核耦合调制的二维HSQC谱,对交叉峰的强度进行分析来精确确定相应的异核－键偶合常数。这一方法的缺点是比较费时。作者在异核多键偶合常数的准确测定方面也做了一些有意义的工作^[4-6 α]。在前文^[5]工作的基础上,本文提出了二维相敏HMQC（Heteronuclear Multiple-Quantum Coherence)和HSQC（Heteronuclear Single-Quantum Coherence)实验,用于准确测定异核－键偶合常数。     

Determination of calpha chemical shift tensor orientation in peptides by dipolar-modulated chemical shift recoupling NMR spectroscopy

We present a new method for determining the orientation of chemical shift tensors in polycrystalline solids with site resolution and demonstrate its application to the determination of the Calpha chemical shift tensor orientation in a model peptide with beta-sheet torsion angles. The tensor orientation is obtained under magic angle spinning by modulating a recoupled chemical shift anisotropy (CSA) pattern with various dipolar couplings. These dipolar-modulated chemical shift patterns constitute the indirect dimension of a 2D spectrum and are resolved according to the isotropic chemical shifts of different sites in the direct dimension. These dipolar-modulated CSA spectra are equivalent to the projection of a 2D static separated-local-field spectrum onto its chemical shift dimension, except that its dipolar dimension is multiplied with a modulation function. Both (13)C-(1)H and (13)C-(15)N dipolar couplings can modulate the CSA spectra of the Calpha site in an amino acid and yield the relative orientations of the chemical shift principal axes to the C-H and C-N bonds. We demonstrate the C-H and C-N modulated CSA experiments on methylmalonic acid and N-tBoc-glycine, respectively. The MAS results agree well with the results of the 2D separated-local-field spectra, thus confirming the validity of this MAS dipolar-modulation approach. Using this technique, we measured the Val Calpha tensor orientation in N-acetylvaline, which has beta-sheet torsion angles. The sigma(11) axis is oriented at 158 degrees (or 22 degrees) from the C-H bond, while the sigma(22) axis is tilted by 144 degrees (or 36 degrees) from the C-N bond. Both the orientations and the magnitude of this chemical shift tensor are in excellent agreement with quantum chemical calculations.     

Dipolar waves as NMR maps of protein structure

The anisotropy of nuclear spin interactions results in a unique mapping of structure to the resonance frequencies and split tings observed in NMR spectra, however, the determination of molecular structure from experimentally measured spectral parameters is complicated by angular ambiguities resulting from the symmetry properties of dipole-dipole and chemical shift interactions. This issue can be addressed through the periodicity inherent in secondary structure elements, which can be used as an index of topology. Distinctive wheel-like patterns are observed in two-dimensional 1H-15N heteronuclear dipolar/15N chemical shift PISEMA (polarization inversion spin-exchange at the magic angle) spectra of helical membrane proteins in highly aligned lipid bilayer samples. One-dimensional dipolar waves are an extension of two-dimensional PISA (polarity index slant angle) wheels to map protein structure in NMR spectra of both highly and weakly aligned samples. Dipolar waves describe the periodic wavelike variations of the magnitudes of the static heteronuclear dipolar couplings as a function of residue number in the absence of chemical shift effects. Weakly aligned samples of proteins display these same effects, primarily as residual dipolar couplings (RDCs), in solution NMR spectra. The corresponding properties of the RDCs in solution NMR spectra of weakly aligned helices represent a convergence of solid-state and solution NMR approaches to structure determination.     

15N chemical shifts and one bond 15N?1H coupling constants in simple amides

An indirect method is employed for determining the ¹⁵N parameters at the natural abundance level in a series of simple acyclic and cyclic amides. The one bond coupling constant, ¹J(¹⁵N¹H), and the ¹⁵N chemical shift are measured as a function of the carbonyl substituent group or the ring size and the nature of the solvent (CCl₄ or H₂O). These ¹⁵N parameters are related to the amide bond structure, the nitrogen configuration and possible intermolecular hydrogen bonding (amide-amide or amide-water).     

HADDOCK: a protein-protein docking approach based on biochemical or biophysical information

The structure determination of protein-protein complexes is a rather tedious and lengthy process, by both NMR and X-ray crystallography. Several methods based on docking to study protein complexes have also been well developed over the past few years. Most of these approaches are not driven by experimental data but are based on a combination of energetics and shape complementarity. Here, we present an approach called HADDOCK (High Ambiguity Driven protein-protein Docking) that makes use of biochemical and/or biophysical interaction data such as chemical shift perturbation data resulting from NMR titration experiments or mutagenesis data. This information is introduced as Ambiguous Interaction Restraints (AIRs) to drive the docking process. An AIR is defined as an ambiguous distance between all residues shown to be involved in the interaction. The accuracy of our approach is demonstrated with three molecular complexes. For two of these complexes, for which both the complex and the free protein structures have been solved, NMR titration data were available. Mutagenesis data were used in the last example. In all cases, the best structures generated by HADDOCK, that is, the structures with the lowest intermolecular energies, were the closest to the published structure of the respective complexes (within 2.0 A backbone RMSD).     

Structure calculation from unambiguous long-range amide and methyl 1H-1H distance restraints for a microcrystalline protein with MAS solid-state NMR spectroscopy

Magic-angle spinning (MAS) solid-state NMR becomes an increasingly important tool for the determination of structures of membrane proteins and amyloid fibrils. Extensive deuteration of the protein allows multidimensional experiments with exceptionally high sensitivity and resolution to be obtained. Here we present an experimental strategy to measure highly unambiguous spatial correlations for distances up to 13 ?. Two complementary three-dimensional experiments, or alternatively a four-dimensional experiment, yield highly unambiguous cross-peak assignments, which rely on four encoded chemical shift dimensions. Correlations to residual aliphatic protons are accessible via synchronous evolution of the (15)N and (13)C chemical shifts, which encode valuable amide-methyl distance restraints. On average, we obtain six restraints per residue. Importantly, 50% of all restraints correspond to long-range distances between residues i and j with |i - j| > 5, which are of particular importance in structure calculations. Using ARIA, we calculate a high-resolution structure for the microcrystalline 7.2 kDa α-spectrin SH3 domain with a backbone precision of ～1.1 ?.