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.     

Chi(1) rotamer populations and angles of mobile surface side chains are accurately predicted by a torsion angle database potential of mean force

The equilibrium angles and distributions of chi(1) rotamers for mobile surface side chains of the small, 63-residue, B1 domain of protein L have been calculated from the static crystal structure by rigid body/torsion angle simulated annealing using a torsion angle database potential of mean force and compared to those deduced by Monte Carlo analysis of side chain residual dipolar couplings measured in solution. Good agreement between theory and experiment is observed, indicating that for side chains undergoing rotamer averaging that is fast on the chemical shift time scale, the equilibrium angles and distribution of chi(1) rotamers are largely determined by the backbone phi/psi torsion angles.     

Side chain orientation from methyl 1H-1H residual dipolar couplings measured in highly deuterated proteins

High-level deuteration is a prerequisite for the study of high molecular weight systems using liquid-state NMR. Here, we present new experiments for the measurement of proton-proton dipolar couplings in CH(2)D methyl groups of (13)C labeled, highly deuterated (70-80%) proteins. (1)H-(1)H residual dipolar couplings (RDCs) have been measured in two alignment media for 57 out of 70 possible methyl containing residues in the 167-residue flavodoxin-like domain of the E. coli sulfite reductase. These data yield information on the orientation of the methyl symmetry axis with respect to the molecular alignment frame. The alignment tensor characteristics were obtained very accurately from a set of backbone RDCs measured on the same protein sample. To demonstrate that accurate structural information is obtained from these data, the measured methyl RDCs for Valine residues are analyzed in terms of chi(1) torsion angles and stereospecific assignment of the prochiral methyl groups. On the basis of the previously determined backbone solution structure of this protein, the methyl RDC data proved sufficient to determine the chi(1) torsion angles in seven out of nine valines, assuming a single-rotamer model. Methyl RDCs are complementary to other NMR data, for example, methyl-methyl NOE, to determine side chain conformation in high molecular weight systems.     

Determination of methyl 13C-15N dipolar couplings in peptides and proteins by three-dimensional and four-dimensional magic-angle spinning solid-state NMR spectroscopy

We describe three- and four-dimensional semiconstant-time transferred echo double resonance (SCT-TEDOR) magic-angle spinning solid-state nuclear magnetic resonance (NMR) experiments for the simultaneous measurement of multiple long-range (15)N-(13)C(methyl) dipolar couplings in uniformly (13)C, (15)N-enriched peptides and proteins with high resolution and sensitivity. The methods take advantage of (13)C spin topologies characteristic of the side-chain methyl groups in amino acids alanine, isoleucine, leucine, methionine, threonine, and valine to encode up to three distinct frequencies ((15)N-(13)C(methyl) dipolar coupling, (15)N chemical shift, and (13)C(methyl) chemical shift) within a single SCT evolution period of initial duration approximately 1(1)J(CC) (where (1)J(CC) approximately 35 Hz, is the one-bond (13)C(methyl)-(13)C J-coupling) while concurrently suppressing the modulation of NMR coherences due to (13)C-(13)C and (15)N-(13)C J-couplings and transverse relaxation. The SCT-TEDOR schemes offer several important advantages over previous methods of this type. First, significant (approximately twofold to threefold) gains in experimental sensitivity can be realized for weak (15)N-(13)C(methyl) dipolar couplings (corresponding to structurally interesting, approximately 3.5 A or longer, distances) and typical (13)C(methyl) transverse relaxation rates. Second, the entire SCT evolution period can be used for (13)C(methyl) and/or (15)N frequency encoding, leading to increased spectral resolution with minimal additional coherence decay. Third, the experiments are inherently "methyl selective," which results in simplified NMR spectra and obviates the use of frequency-selective pulses or other spectral filtering techniques. Finally, the (15)N-(13)C cross-peak buildup trajectories are purely dipolar in nature (i.e., not influenced by J-couplings or relaxation), which enables the straightforward extraction of (15)N-(13)C(methyl) distances using an analytical model. The SCT-TEDOR experiments are demonstrated on a uniformly (13)C, (15)N-labeled peptide, N-acetyl-valine, and a 56 amino acid protein, B1 immunoglobulin-binding domain of protein G (GB1), where the measured (15)N-(13)C(methyl) dipolar couplings provide site-specific information about side-chain dihedral angles and the packing of protein molecules in the crystal lattice.     

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.     

Solid-state 19F NMR spectroscopy reveals that Trp41 participates in the gating mechanism of the M2 proton channel of influenza A virus

The integral membrane protein M2 of influenza A virus assembles as a tetrameric bundle to form a proton-conducting channel that is activated by low pH. The side chain of His37 in the transmembrane alpha-helix is known to play an important role in the pH activation of the proton channel. It has also been suggested that Trp41, which is located in an adjacent turn of the helix, forms part of the gating mechanism. Here, a synthetic 25-residue peptide containing the M2 transmembrane domain was labeled with 6F-Trp41 and studied in lipid membranes by solid-state 19F NMR. We monitored the pH-dependent differences in the 19F dipolar couplings and motionally narrowed chemical shift anisotropies of this 6F-Trp41 residue, and we discuss the pH activation mechanism of the H+ channel. At pH 8.0, the structural parameters implicate an inactivated state, while at pH 5.3 the tryptophan conformation represents the activated state. With the aid of COSMOS force field simulations, we have obtained new side-chain torsion angles for Trp41 in the inactivated state (chi1 = -100 degrees +/- 10 degrees , chi2 = +110 degrees +/- 10 degrees ), and we predict a most probable activated state with chi1 = -50 degrees +/- 10 degrees and chi2 = +115 degrees +/- 10 degrees . We have also validated the torsion angles of His37 in the inactivated state as chi1 = -175 degrees +/- 10 degrees and chi2 = -170 degrees +/- 10 degrees .     

Fluorescence of cis-1-amino-2-(3-indolyl)cyclohexane-1-carboxylic acid: a single tryptophan chi(1) rotamer model

A constrained derivative, cis-1-amino-2-(3-indolyl)cyclohexane-1-carboxylic acid, cis-W3, was designed to test the rotamer model of tryptophan photophysics. The conformational constraint enforces a single chi(1) conformation, analogous to the chi(1) = 60 degrees rotamer of tryptophan. The side-chain torsion angles in the X-ray structure of cis-W3 were chi(1) = 58.5 degrees and chi(2) = -88.7 degrees. Molecular mechanics calculations suggested two chi(2) rotamers for cis-W3 in solution, -100 degrees and 80 degrees, analogous to the chi(2) = +/-90 degrees rotamers of tryptophan. The fluorescence decay of the cis-W3 zwitterion was biexponential with lifetimes of 3.1 and 0.3 ns at 25 degrees C. The relative amplitudes of the lifetime components match the chi(2) rotamer populations predicted by molecular mechanics. The longer lifetime represents the major chi(2) = -100 degrees rotamer. The shorter lifetime represents the minor chi(2) = 80 degrees rotamer having the ammonium group closer to C4 of the indole ring (labeled C5 in the cis-W3 X-ray structure). Intramolecular excited-state proton transfer occurs at indole C4 in the tryptophan zwitterion (Saito, I.; Sugiyama, H.; Yamamoto, A.; Muramatsu, S.; Matsuura,T. J. Am. Chem. Soc. 1984, 106, 4286-4287). Photochemical isotope exchange experiments showed that H-D exchange occurs exclusively at C5 in the cis-W3 zwitterion, consistent with the presence of the chi(2) = 80 degrees rotamer in solution. The rates of two nonradiative processes, excited-state proton and electron transfer, were measured for individual chi(2) rotamers. The excited-state proton-transfer rate was determined from H-D exchange and fluorescence lifetime data. The excited-state electron-transfer rate was determined from the temperature dependence of the fluorescence lifetime. The major quenching process in the -100 degrees rotamer is electron transfer from the excited indole to carboxylate. Electron transfer also occurs in the 80 degrees rotamer, but the major quenching process is intramolecular proton transfer. Both quenching processes are suppressed by deprotonation of the amino group. The results for cis-W3 provide compelling evidence that the complex fluorescence decay of the tryptophan zwitterion originates in ground-state heterogeneity with the different lifetimes primarily reflecting different intramolecular excited-state proton- and electron-transfer rates in various rotamers.     

Statistically based reduced representation of amino acid side chains

Preferred conformations of amino acid side chains have been well established through statistically obtained rotamer libraries. Typically, these provide bond torsion angles allowing a side chain to be traced atom by atom. In cases where it is desirable to reduce the complexity of a protein representation or prediction, fixing all side-chain atoms may prove unwieldy. Therefore, we introduce a general parametrization to allow positions of representative atoms (in the present study, these are terminal atoms) to be predicted directly given backbone atom coordinates. Using a large, culled data set of amino acid residues from high-resolution protein crystal structures, anywhere from 1 to 7 preferred conformations were observed for each terminal atom of the non-glycine residues. Side-chain length from the backbone C(alpha) is one of the parameters determined for each conformation, which should itself be useful. Prediction of terminal atoms was then carried out for a second, nonredundant set of protein structures to validate the data set. Using four simple probabilistic approaches, the Monte Carlo style prediction of terminal atom locations given only backbone coordinates produced an average root mean-square deviation (RMSD) of approximately 3 A from the experimentally determined terminal atom positions. With prediction using conditional probabilities based on the side-chain chi(1) rotamer, this average RMSD was improved to 1.74 A. The observed terminal atom conformations therefore provide reasonable and potentially highly accurate representations of side-chain conformation, offering a viable alternative to existing all-atom rotamers for any case where reduction in protein model complexity, or in the amount of data to be handled, is desired. One application of this representation with strong potential is the prediction of charge density in proteins. This would likely be especially valuable on protein surfaces, where side chains are much less likely to be fixed in single rotamers. Prediction of ensembles of structures provides a method to determine the probability density of charge and atom location; such a prediction is demonstrated graphically.     

Evaluation of backbone proton positions and dynamics in a small protein by liquid crystal NMR spectroscopy

NMR measurements of a large set of protein backbone one-bond dipolar couplings have been carried out to refine the structure of the third IgG-binding domain of Protein G (GB3), previously solved by X-ray crystallography at a resolution of 1.1 A. Besides the commonly used bicelle, poly(ethylene glycol), and filamentous phage liquid crystalline media, dipolar couplings were also measured when the protein was aligned inside either positively or negatively charged stretched acrylamide gels. Refinement of the GB3 crystal structure against the (13)C(alpha)-(13)C' and (13)C'-(15)N dipolar couplings improves the agreement between experimental and predicted (15)N-(1)H(N) as well as (13)C(alpha)-(1)H(alpha) dipolar couplings. Evaluation of the peptide bond N-H orientations shows a weak anticorrelation between the deviation of the peptide bond torsion angle omega from 180 degrees and the angle between the N-H vector and the C'-N-C(alpha) plane. The slope of this correlation is -1, indicating that, on average, pyramidalization of the peptide N contributes to small deviations from peptide bond planarity ( = 179.3 +/- 3.1 degrees ) to the same degree as true twisting around the C'-N bond. Although hydrogens are commonly built onto crystal structures assuming the N-H vector orientation falls on the line bisecting the C'-N-C(alpha) angle, a better approximation adjusts the C(alpha)-C'-N-H torsion angle to -2 degrees. The (15)N-(1)H(N) dipolar data do not contradict the commonly accepted motional model where angular fluctuations of the N-H bond orthogonal to the peptide plane are larger than in-plane motions, but the amplitude of angular fluctuations orthogonal the C(alpha)(i-1)-N(i)-C(alpha)(i) plane exceeds that of in-plane motions by at most 10-15 degrees. Dipolar coupling analysis indicates that for most of the GB3 backbone, the amide order parameters, S, are highly homogeneous and vary by less than +/-7%. Evaluation of the H(alpha) proton positions indicates that the average C(alpha)-H(alpha) vector orientation deviates by less than 1 degrees from the direction that makes ideal tetrahedral angles with the C(alpha)-C(beta) and C(alpha)-N vectors.     

Self-consistent Karplus parametrization of 3J couplings depending on the polypeptide side-chain torsion chi1.

Recently proposed self-consistent 3J coupling analysis (Schmidt, J. M.; Blümel, M.; L?hr, F.; Rüterjans, H. J. Biomol. NMR 1999, 14, 1-12) has been carried out to calibrate Karplus parameters constituting the empirical dependence of 3J coupling constants on the chi1 dihedral angle in amino acid side chains. The procedure involves simultaneous least-squares optimization of six sets of three Karplus coefficients related to all six 3J coupling types accessible in 15N,13C-labeled proteins. A simple concept of fundamental and incremental component couplings is proposed to account for substituent effects, eventually yielding amino acid topology-specific Karplus parameters. The method is exemplified with recombinant Desulfovibrio vulgaris flavodoxin (147 amino acids, 16 kDa) with reference to a total of 749 experimental 3JHalpha,Hbeta, 3JN',Hbeta, 3JC',Hbeta, 3JHalpha,Cgamma, 3JN',Cgamma, and 3JC',Cgamma coupling constants. Unlike other parametrizations, the present method does not make reference to X-ray coordinates, so that the Karplus coefficients obtained are not influenced by differences between solution and crystal states. Cross validation using X-ray torsion angles demonstrates the improvement relative to previous parametrizations. The Karplus coefficients derived are applicable to other proteins, too. Parameter refinement also yields a series of chi1 torsion angles, providing valuable constraints for protein structure determination, as well as optional parameters of local angular mobility in the contexts of Gaussian random fluctuation or a three-site jump model. The procedure permits automatic stereospecific assignments of Hbeta and Cgamma chemical shifts. The majority of the flavodoxin side-chain conformations agrees with high-resolution X-ray structures of the protein. Marked deviations between NMR and X-ray datasets are attributed to different rotameric states due to crystal-packing effects and to conformational equilibria between multiple chi1 rotamers.     

Side-chain conformation of the M2 transmembrane peptide proton channel of influenza a virus from 19F solid-state NMR

The M2 transmembrane peptide (M2TMP) of the influenza A virus forms a tetrameric helical bundle that acts as a proton-selective channel important in the viral life cycle. The side-chain conformation of the peptide is largely unknown and is important for elucidating the proton-conducting mechanism and the channel stability. Using a 19F spin diffusion NMR technique called CODEX, we have measured the oligomeric states and interhelical side chain-side chain 19F-19F distances at several residues using singly fluorinated M2TMP bound to DMPC bilayers. 19F CODEX data at a key residue of the proton channel, Trp41, confirm the tetrameric state of the peptide and yield a nearest-neighbor interhelical distance of approximately 11 A under both neutral and acidic pH. Since the helix orientation is precisely known from previous 15N NMR experiments and the backbone channel diameter has a narrow allowed range, this 19F distance constrains the Trp41 side-chain conformation to t90 (chi1 approximately 180 degrees , chi2 approximately 90 degrees ). This Trp41 rotamer, combined with a previously measured 15N-13C distance between His37 and Trp411, suggests that the His37 rotamer is t-160. The implication of the proposed (His37, Trp41) rotamers to the gating mechanism of the M2 proton channel is discussed. Binding of the antiviral drug amantadine to the peptide does not affect the F-F distance at Trp41. Interhelical 19F-19F distances are also measured at residues 27 and 38, each mutated to 4-19F-Phe. For V27F-M2TMP, the 19F-19F distances suggest a mixture of dimers and tetramers, whereas the L38F-M2TMP data indicate two tetramers of different sizes, suggesting side chain conformational heterogeneity at this lipid-facing residue. This work shows that 19F spin diffusion NMR is a valuable tool for determining long-range intermolecular distances that shed light on the mechanism of action and conformational heterogeneity of membrane protein oligomers.     

Measurement of (13)C(alpha)-(13)C(beta) dipolar couplings in (15)N,(13)C,(2)H-labeled proteins: application to domain orientation in maltose binding protein.

TROSY-based HN(CO)CA 2D and 3D pulse schemes are presented for measurement of (13)C(alpha)-(13)C(beta) dipolar couplings in high molecular weight (15)N,(13)C,(2)H-labeled proteins. In one approach, (13)C(alpha)-(13)C(beta) dipolar couplings are obtained directly from the time modulation of cross-peak intensities in a set of 2D (15)N-(1)HN correlated spectra recorded in both the presence and absence of aligning media. In a second approach 3D data sets are recorded with (13)C(alpha)-(13)C(beta) couplings encoded in a frequency dimension. The utility of the experiments is demonstrated with an application to an (15)N,(13)C,(2)H-labeled sample of the ligand free form of maltose binding protein. A comparison of experimental dipolar couplings with those predicted from the X-ray structure of the apo form of this two-domain protein establishes that the relative orientation of the domains in solution and in the crystal state are very similar. This is in contrast to the situation for maltose binding protein in complex with beta-cyclodextrin where the solution structure can be generated from the crystal state via a 11 degrees domain closure.     

Docking of protein-protein complexes on the basis of highly ambiguous intermolecular distance restraints derived from 1H/15N chemical shift mapping and backbone 15N-1H residual dipolar couplings using conjoined rigid body/torsion angle dynamics

A simple and reliable method for docking protein-protein complexes using (1)H(N)/(15)N chemical shift mapping and backbone (15)N-(1)H residual dipolar couplings is presented and illustrated with three complexes (EIN-HPr, IIA(Glc)-HPr, and IIA(Mtl)-HPr) of known structure. The (1)H(N)/(15)N chemical shift mapping data are transformed into a set of highly ambiguous, intermolecular distance restraints (comprising between 400 and 3000 individual distances) with translational and some degree of orientational information content, while the dipolar couplings provide information on relative protein-protein orientation. The optimization protocol employs conjoined rigid body/torsion angle dynamics in simulated annealing calculations. The target function also comprises three nonbonded interactions terms: a van der Waals repulsion term to prevent atomic overlap, a radius of gyration term (E(rgyr)) to avoid expansion at the protein-protein interface, and a torsion angle database potential of mean force to bias interfacial side chain conformations toward physically allowed rotamers. For the EIN-HPr and IIA(Glc)-HPr complexes, all structures satisfying the experimental restraints (i.e., both the ambiguous intermolecular distance restraints and the dipolar couplings) converge to a single cluster with mean backbone coordinate accuracies of 0.7-1.5 A. For the IIA(Mtl)-HPr complex, twofold degeneracy remains, and the structures cluster into two distinct solutions differing by a 180 degrees rotation about the z axis of the alignment tensor. The correct and incorrect solutions which have mean backbone coordinate accuracies of approximately 0.5 and approximately 10.5 A, respectively, can readily be distinguished using a variety of criteria: (a) examination of the overall (1)H(N)/(15)N chemical shift perturbation map (because the incorrect cluster predicts the presence of residues at the interface that experience only minimal chemical shift perturbations; this information is readily incorporated into the calculations in the form of ambiguous intermolecular repulsion restraints); (b) back-calculation of dipolar couplings on the basis of molecular shape; or (c) the E(rgyr) distribution which, because of its global nature, directly reflects the interfacial packing quality. This methodology should be particularly useful for high throughput, NMR-based, structural proteomics.     

Determination of sugar structures in solution from residual dipolar coupling constants: methodology and application to methyl beta-D-xylopyranoside

We have developed methodology for the determination of solution structures of small molecules from residual dipolar coupling constants measured in dilute liquid crystals. The power of the new technique is demonstrated by the determination of the structure of methyl beta-d-xylopyranoside (I) in solution. An oriented sample of I was prepared using a mixture of C(12)E(5) and hexanol in D(2)O. Thirty residual dipolar coupling constants, ranging from -6.44 to 4.99 Hz, were measured using intensity-based J-modulated NMR techniques. These include 15 D(HH), 4 (1)D(CH), and 11 (n)D(CH) coupling constants. The accuracy of the dipolar coupling constants is estimated to be < +/- 0.02 Hz. New constant-time HMBC NMR experiments were developed for the measurement of (n)D(CH) coupling constants, the use of which was crucial for the successful structure determination of I, as they allowed us to increase the number of fitted parameters. The structure of I was refined using a model in which the directly bonded interatom distances were fixed at their ab initio values, while 16 geometrical and 5 order parameters were optimized. These included 2 CCC and 6 CCH angles, and 2 CCCC and 6 CCCH dihedral angles. Vibrationally averaged dipolar coupling constants were used during the refinement. The refined solution structure of I is very similar to that obtained by ab initio calculations, with 11 bond and dihedral angles differing by 0.8 degrees or less and the remaining 5 parameters differing by up to 3.3 degrees . Comparison with the neutron diffraction structure showed larger differences attributable to crystal packing effects. Reducing the degree of order by using dilute liquid crystalline media in combination with precise measurement of small residual dipolar coupling constants, as shown here, is a way of overcoming the limitation of strongly orienting liquid crystals associated with the complexity of (1)H NMR spectra for molecules with more than 12 protons.     

An investigation of solute-liquid crystal interaction in hydrogen bonded complexes of palmitic acid and p-octyloxybenzoic acid

The orientational order of a liquid crystalline phase which has a specific solute-liquid crystal interaction was investigated using nuclear magnetic resonance. Three isotopically substituted species of palmitic acid (palmitic acid-d₃₁, 1-¹³C-2.2-H₂-palmitic acid-d₂₉ and 2,2,3,3-H₄-palmitic acid-d₂₇) were dissolved in the liquid crystal p-octyloxybenzoic acid (p-OOBA) and the proton, deuteron and carbon 13 NMR spectra recorded as a function of temperature. ¹H-¹³H dipolar couplings were observed using a spin echo pulse sequence which removes heteronuclear dipolar couplings to the chain deuterons. In the case of the carbon 13 labelled compound, ¹H-¹³C dipolar couplings could be observed by applying an additional refocusing pulse to the ¹³C spins. The dipolar and quadrupolar couplings were used to calculate the complete orientational order matrix of the alpha methylene segment of palmitic acid in p-OOBA. The liquid crystal was shown to largely determine the orientational order of the head group and this was attributed to intermolecular hydrogen bonding. The dipolar and quadrupolar couplings for the rest of the chain were interpreted in terms of a mean field equilibrium statistical model, based on the Samulski Inertial Frame Model. Hydrogen bonding was shown to be of greater importance in the orientational ordering of the solutes in the liquid crystal than are electrostatic interactions in the ordering of the amphiphile in the potassium palmitate/water system.     

Revealing the configuration and crystal packing of organic compounds by solid-state NMR spectroscopy: methoxycarbonylurea, a case study

The molecular configuration and intermolecular arrangement of polycrystalline methoxycarbonylurea (MCU) has been studied by a combination of chemical editing, rotational echo double resonance (REDOR) spectroscopy and ab initio calculations. From the multispin IS(n) REDOR experiments several dipolar couplings were determined and converted into distance constraints. Intra- and intermolecular dipolar couplings were distinguished by isotope dilution. The configuration of the MCU molecule can be determined from three torsion angles Psi1, Psi2, and Psi3. Ab initio calculations showed that these angles are either 0 degrees or 180 degrees (Z or E). From the REDOR experiments, the E configuration was found for Psi1 and Psi2 and the Z configuration for Psi3. Thus the configuration of MCU in the solid state was determined to be EEZ. Distance constraints for the intermolecular arrangement of MCU were obtained by performing REDOR experiments on 13C15N2 MCU with different degrees of isotope dilution and on a cocrystallized 1:1 mixture of 13C(urea) MCU and 15N(amide) MCU. By combining these distance constraints with molecular modeling, three different possible packing motifs for MCU molecules were found. The molecules in these motifs are arranged as linear chains with methoxy groups at the borders of the chains. All the intermolecular hydrogen bond donors and acceptors in the interior of the chain are saturated.     

Determining valine side-chain rotamer conformations in proteins from methyl 13C chemical shifts: application to the 360 kDa half-proteasome

A method is presented for determining Val side-chain χ(1) rotamer distributions in proteins based exclusively on measured (13)C(γ1) and (13)C(γ2) chemical shifts. The approach selects an ensemble of 20 χ(1) values, calculates average methyl (13)C(γ1,γ2) chemical shifts via theoretical quantum chemical calculations and maximizes the agreement with the experimentally measured shifts using a genetic algorithm. The methodology is validated with an application involving six proteins for which (13)C(γ) chemical shifts and three-bond methyl-backbone scalar couplings are available. The utility of the methodology is demonstrated with an application to the 360 kDa 'half-proteasome' where the χ(1) rotameric distributions of Val residues are calculated on the basis of chemical shifts. For the most part the χ(1) profiles so obtained compare very well with those generated from the high-resolution (2.3 ?) X-ray structure of the proteasome. Both NMR and X-ray distributions are cross-validated by comparing calculated (1)H-(13)C methyl residual dipolar couplings with measured values, and the level of agreement is at least as good for the NMR derived χ(1) values. Notably, as the resolution of the X-ray data improves (rotamer distributions from 3.4 and 2.3 ? X-ray structures are compared with the NMR data), the agreement with the NMR gets significantly better. This emphasizes the importance of NMR approaches for the study of high molecular weight complexes that can be recalcitrant to high resolution X-ray analysis.     

(13)C-(1)H and (13)C-(13)C NMR J-couplings in (13)C-labeled N-acetyl-neuraminic acid: correlations with molecular structure

N-acetyl-neuraminic acid (Neu5Ac, 2) was prepared enzymatically containing single sites of (13)C-enrichment at C1, C2, and C3. Aqueous solutions of the three (13)C isotopomers were studied by (1)H and (13)C NMR spectroscopy at p(2)H 2 and pH 8 to obtain J(CH) and J(CC) values involving the labeled carbons. Experimental studies were complemented by DFT calculations of the same set of J-couplings in protonated and ionized structural mimics of 2 to determine how well theoretical predictions match the experimental findings in saccharides bearing ionizable functionality. Results show that: (a) (2)J(C2,H3ax/eq) values in 2 depend on anomeric configuration, thus complementing (3)J(C1,H3ax/eq) behavior, (b) J(CH) and J(CC) values involving C2 depend on anomeric configuration, the C1-C2 bond torsion, and solution pH, and (c) long-range (4)J(C2,H7) is sensitive to glycerol side-chain conformation. Intraring J(HH) and most (2)J(CH), (3)J(CH), (2)J(CC), and (3)J(CC) involving C1-C3 of 2 appear largely unaffected by the ionization state of the carboxyl group. In vacuo and solvated DFT calculations of geminal and vicinal J(CH) and J(CC) values are similar and reproduce the experimental data well, but better agreement with experiment was observed for (1)J(C1,C2) in the solvated calculations. The present work provides new information for future treatments of trans-glycoside couplings involving Neu5Ac residues by (a) providing new standard values of intraring J(CC) for coupling pathways that mimic those for trans-glycoside J(CC), (b) identifying potential effects of solution pH on trans-glycoside couplings inferred through the behavior of related intraring couplings, and (c) providing specific guidelines for more reliable DFT predictions of J(CH) and J(CC) values in ionizable saccharides.     

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.     

Measurement of Large Dipolar Couplings of a Liquid Crystal with Terminal Phenyl Rings and Estimation of the Order Parameters

NMR spectroscopy is a powerful means of studying liquid‐crystalline systems at atomic resolutions. Of the many parameters that can provide information on the dynamics and order of the systems, ¹H–¹³C dipolar couplings are an important means of obtaining such information. Depending on the details of the molecular structure and the magnitude of the order parameters, the dipolar couplings can vary over a wide range of values. Thus the method employed to estimate the dipolar couplings should be capable of estimating both large and small dipolar couplings at the same time. For this purpose, we consider here a two‐dimensional NMR experiment that works similar to the insensitive nuclei enhanced by polarization transfer (INEPT) experiment in solution. With the incorporation of a modification proposed earlier for experiments with low radio frequency power, the scheme is observed to enable a wide range of dipolar couplings to be estimated at the same time. We utilized this approach to obtain dipolar couplings in a liquid crystal with phenyl rings attached to either end of the molecule, and estimated its local order parameters.