Rotational diffusion tensor of nucleic acids from 13C NMR relaxation

Rotational diffusion properties have been derived for the DNA dodecamer d(CGCGAATTCGCG)(2) from (13)C R(1rho) and R(1) measurements on the C(1'), C(3'), and C(4') carbons in samples uniformly enriched in (13)C. The narrow range of C-H bond vector orientations relative to the DNA axis make the analysis particularly sensitive to small structural deviations. As a result, the R(1rho)/R(1) ratios are found to fit poorly to the crystal structures of this dodecamer, but well to a recent solution NMR structure, determined in liquid crystalline media, even though globally the structures are quite similar. A fit of the R(1rho)/R(1) ratios to the solution structure is optimal for an axially symmetric rotational diffusion model, with a diffusion anisotropy, D(\|)/D(perpendicular), of 2.1+/-0.4, and an overall rotational correlation time, (2D(\|)+4D(perpendicular))(-1), of 3.35 ns at 35 degrees C in D(2)O, in excellent agreement with values obtained from hydrodynamic modeling.

Direct observation of dipolar couplings between distant protons in weakly aligned nucleic acids

Under conditions where macromolecules are aligned very weakly with respect to an external magnetic field, Brownian diffusion no longer averages internuclear dipole-dipole interactions to zero. The resulting residual dipolar coupling, although typically 3 orders of magnitude weaker than in a fully aligned sample, can readily be measured by solution NMR methods. To date, application of this idea has focused primarily on pairs of nuclei separated by one or two covalent bonds, where the internuclear separation is known and the measured dipolar coupling provides direct information on the orientation of the internuclear vector. A method is described that allows observation of dipolar interactions over much larger distances. By decoupling nearest-neighbor interactions, it is readily possible to observe direct dipolar interactions between protons separated by up to 12 A. The approach is demonstrated for the DNA dodecamer d(CGCGAATTCGCG)2, where direct interactions are observed between protons up to three base pairs apart.

Overall structure and sugar dynamics of a DNA dodecamer from homo- and heteronuclear dipolar couplings and 31P chemical shift anisotropy

The solution structure of d(CGCGAATTCGCG)(2) has been determined on the basis of an exceptionally large set of residual dipolar couplings. In addition to the heteronuclear (13)C-(1)H and (15)N-(1)H and qualitative homonuclear (1)H-(1)H dipolar couplings, previously measured in bicelle medium, more than 300 quantitative (1)H-(1)H and 22 (31)P-(1)H dipolar restraints were obtained in liquid crystalline Pf1 medium, and 22 (31)P chemical shift anisotropy restraints. High quality DNA structures can be obtained solely on the basis of these new restraints, and these structures are in close agreement with those calculated previously on the basis of (13)C-(1)H and (15)N-(1)H dipolar couplings. In the newly calculated structures, (31)P-(1)H dipolar and (3)JsubH3(')Psub couplings and (31)P CSA data restrain the phosphodiester backbone torsion angles. The final structure represents a quite regular B-form helix with a modest bending of approximately 10 degrees, which is essentially independent of whether or not electrostatic terms are used in the calculation. Combined, the number of homo- and heteronuclear dipolar couplings significantly exceeds the number of degrees of freedom in the system. Results indicate that the dipolar coupling data cannot be fit by a single structure, but are compatible with the presence of rapid equilibria between C2(')-endo and C3(')-endo deoxyribose puckers (sugar switching). The C2(')-H2(')/H2(") dipolar couplings in B-form DNA are particularly sensitive to sugar pucker and yield the largest discrepancies when fit to a single structure. To resolve these discrepancies, we suggest a simplified dipolar coupling analysis that yields N/S equilibria for the ribose sugar puckers, which are in good agreement with previous analyses of NMR J(HH) couplings, with a population of the minor C3(')-endo form higher for pyrimidines than for purines.

Calcium dependence of the interaction between calmodulin and anthrax edema factor

Edema factor (EF), a toxin from Bacillus anthracis (anthrax), possesses adenylyl cyclase activity and requires the ubiquitous Ca2+-sensor calmodulin (CaM) for activity. CaM can exist in three major structural states: an apo state with no Ca2+ bound, a two Ca2+ state with its C-terminal domain Ca2+-loaded, and a four Ca2+ state in which the lower Ca2+ affinity N-terminal domain is also ligated. Here, the interaction of EF with the three Ca2+ states of CaM has been examined by NMR spectroscopy and changes in the Ca2+ affinity of CaM in the presence of EF have been determined by flow dialysis. Backbone chemical shift perturbations of CaM show that EF interacts weakly with the N-terminal domain of apoCaM. The C-terminal CaM domain only engages in the interaction upon Ca2+ ligation, rendering the overall interaction much tighter. In the presence of EF, the C-terminal domain binds Ca2+ with higher affinity, but loses binding cooperativity, whereas the N-terminal domain exhibits strongly reduced Ca2+ affinity. As judged by chemical shift differences, the N-terminal CaM domain remains bound to EF upon subsequent Ca2+ ligation. This Ca2+ dependence of the EF-CaM interaction differs from that observed for most other CaM targets, which normally interact only with the Ca2+-bound CaM domains and become active following the transition to the four Ca2+ state.

Insights into the mobility of methyl-bearing side chains in proteins from (3)J(CC) and (3)J(CN) couplings

Side-chain dynamics in proteins can be characterized by the NMR measurement of (13)C and (2)H relaxation rates. Evaluation of the corresponding spectral densities limits the slowest motions that can be studied quantitatively to the time scale on which the overall molecular tumbling takes place. A different measure for the degree of side-chain order about the C(alpha)-C(beta) bond (chi(1) angle) can be derived from (3)J(C)(')(-)(C)(gamma) and (3)J(N)(-)(C)(gamma) couplings. These couplings can be measured at high accuracy, in particular for Thr, Ile, and Val residues. In conjunction with the known backbone structures of ubiquitin and the third IgG-binding domain of protein G, and an extensive set of (13)C-(1)H side-chain dipolar coupling measurements in oriented media, these (3)J couplings were used to parametrize empirical Karplus relationships for (3)J(C)(')(-)(C)(gamma) and (3)J(N)(-)(C)(gamma). These Karplus curves agree well with results from DFT calculations, including an unusual phase shift, which causes the maximum (3)J(CC) and (3)J(CN) couplings to occur for dihedral angles slightly smaller than 180 degrees, particularly noticeable in Thr residues. The new Karplus curves permit determination of rotamer populations for the chi(1) torsion angles. Similar rotamer populations can be derived from side-chain dipolar couplings. Conversion of these rotamer populations into generalized order parameters, S(J)(2) and S(D)(2), provides a view of side-chain dynamics that is complementary to that obtained from (13)C and (2)H relaxation. On average, results agree well with literature values for (2)H-relaxation-derived S(rel)(2) values in ubiquitin and HIV protease, but also identify a fraction of residues for which S(J,D)(2) < S(rel)(2). This indicates that some of the rotameric averaging occurs on a time scale too slow to be observable in traditional relaxation measurements.

Thermodynamic rearrangement synthesis and NMR structures of C1, C3, and T isomers of C60H36

The structures of three C60H36 isomers, produced by high-temperature transfer hydrogenation of C(60) in a 9,10-dihydroanthracene melt, was accomplished by 2D (1)H-detected NMR experiments, recorded at 800 MHz. The unsymmetrical C(1) isomer is found to be the most abundant one (60-70%), followed by the C(3) isomer (25-30%) and the least abundant T isomer (2-5%). All three isomers are closely related in structure and have three vicinal hydrogens located on each of the 12 pentagons. Facile hydrogen migration on the fullerene surface during annealing at elevated temperatures is believed to be responsible for the preferential formation of these thermodynamically most stable C60H36 isomers. This hypothesis was further supported by thermal conversion of C60H36 isomers to a single C(3v) isomer of C60H18.

Quaternary structure of hemoglobin in solution

Many important proteins perform their physiological functions under allosteric control, whereby the binding of a ligand at a specific site influences the binding affinity at a different site. Allosteric regulation usually involves a switch in protein conformation upon ligand binding. The energies of the corresponding structures are comparable, and, therefore, the possibility that a structure determined by x-ray diffraction in the crystalline state is influenced by its intermolecular contacts, and thus differs from the solution structure, cannot be excluded. Here, we demonstrate that the quaternary structure of tetrameric human normal adult carbonmonoxy-hemoglobin can readily be determined in solution at near-physiological conditions of pH, ionic strength, and temperature by NMR measurement of (15)N-(1)H residual dipolar couplings in weakly oriented samples. The structure is found to be a dynamic intermediate between two previously solved crystal structures, known as the R and R2 states. Exchange broadening at the subunit interface points to a rapid equilibrium between different structures that presumably include the crystallographically observed states.

Weak alignment offers new NMR opportunities to study protein structure and dynamics

Protein solution nuclear magnetic resonance (NMR) can be conducted in a slightly anisotropic environment, where the orientational distribution of the proteins is no longer random. In such an environment, the large one-bond internuclear dipolar interactions no longer average to zero and report on the average orientation of the corresponding vectors relative to the magnetic field. The desired very weak ordering, on the order of 10(-3), can be induced conveniently by the use of aqueous nematic liquid crystalline suspensions or by anisotropically compressed hydrogels. The resulting residual dipolar interactions are scaled down by three orders of magnitude relative to their static values, but nevertheless can be measured at high accuracy. They are very precise reporters on the average orientation of bonds relative to the molecular alignment frame, and they can be used in a variety of ways to enrich our understanding of protein structure and function. Applications to date have focused primarily on validation of structures, determined by NMR, X-ray crystallography, or homology modeling, and on refinement of structures determined by conventional NMR approaches. Although de novo structure determination on the basis of dipolar couplings suffers from a severe multiple minimum problem, related to the degeneracy of dipolar coupling relative to inversion of the internuclear vector, a number of approaches can address this problem and potentially can accelerate the NMR structure determination process considerably. In favorable cases, where large numbers of dipolar couplings can be measured, inconsistency between measured values can report on internal motions.

Long-range magnetization transfer between uncoupled nuclei by dipole-dipole cross-correlated relaxation: a precise probe of beta-sheet geometry in proteins

Interference between dipolar interactions in covalently linked (13)C-(1)H and nonlinked (1)H-(1)H pairs can be used to generate antiphase magnetization between noncoupled spins. The buildup rate of such antiphase terms is highly sensitive to local geometry, in particular the interproton distance and the (13)C-(1)H-(1)H internuclear angle. These rates have been measured for opposing C(alpha)H(alpha) pairs in antiparallel beta-sheets in the third Igg-binding domain of protein G (GB3) and in HIV protease, complexed with the inhibitor DMP323. For GB3, good agreement with the 1.1-A crystal structure is found. However, this agreement rapidly deteriorates with decreasing resolution of the corresponding X-ray structure. For HIV protease, two separate crystal structures that differ by less than 0.2 A from one another exhibit lower agreement in their predicted cross-correlated relaxation rates relative to one another than is found between experimental rates and the average of the rates predicted for the two structures. These data indicate that quantitative measurement of these cross-correlated relaxation rates can provide highly accurate structural information in macromolecules.

Structure and orientation of a G protein fragment in the receptor bound state from residual dipolar couplings

Residual dipolar couplings for a ligand that is in fast exchange between a free state and a state where it is bound to a macroscopically ordered membrane protein carry precise information on the structure and orientation of the bound ligand. The couplings originate in the bound state but can be detected on the free ligand using standard high resolution NMR. This approach is used to study an analog of the C-terminal undecapeptide of the alpha-subunit of the heterotrimeric G protein transducin when bound to photo-activated rhodopsin. Rhodopsin is the major constituent of disk-shaped membrane vesicles from rod outer segments of bovine retinas, which align spontaneously in the NMR magnet. Photo-activation of rhodopsin triggers transient binding of the peptide, resulting in measurable dipolar contributions to 1J(NH) and 1J(CH) splittings. These dipolar couplings report on the time-averaged orientation of bond vectors in the bound peptide relative to the magnetic field, i.e. relative to the membrane normal. Approximate distance restraints of the bound conformation were derived from transferred NOEs, as measured from the difference of NOESY spectra recorded prior to and after photo-activation. The N-terminal eight residues of the bound undecapeptide adopt a near-ideal alpha-helical conformation. The helix is terminated by an alpha(L) type C-cap, with Gly9 at the C' position in the center of the reverse turn. The angle between the helix axis and the membrane normal is 40 degrees (+/-4) degrees. Peptide protons that make close contact with the receptor are identified by analysis of the NOESY cross-relaxation pattern and include the hydrophobic C terminus of the peptide.

Measurement of long-range 1H-1H dipolar couplings in weakly aligned proteins

Measurement of 1H-1H dipolar couplings in macromolecules, weakly oriented by a dilute liquid crystalline medium, is generally limited to the largest such interactions. By removing dipolar couplings to nearest neighbors, either by decoupling, deuteration, or both, more remote interactions become accessible. The approach is demonstrated for measurement of amide-amide interactions in the proteins calmodulin and ubiquitin and permits observation of direct dipolar couplings between protons up to 7 A apart. Quantitative evaluation of 1H-1H dipolar couplings measured in ubiquitin shows excellent agreement with its solution structure.

Evaluation of uncertainty in alignment tensors obtained from dipolar couplings

Residual dipolar couplings and their corresponding alignment tensors are useful for structural analysis of macromolecules. The error in an alignment tensor, derived from residual dipolar couplings on the basis of a known structure, is determined not only by the accuracy of the measured couplings but also by the uncertainty in the structure (structural noise). This dependence is evaluated quantitatively on the basis of simulated structures using Monte-Carlo type analyses. When large numbers of dipolar couplings are available, structural noise is found to result in a systematic underestimate of the magnitude of the alignment tensor. Particularly in cases where only few dipolar couplings are available, structural noise can cause significant errors in best-fitted alignment tensor values, making determination of the relative orientation of small fragments and evaluation of local backbone mobility from dipolar couplings difficult. An example for the protein ubiquitin demonstrates the inherent limitations in characterizing motions on the basis of local alignment tensor magnitudes.