Internal mobility of cyclic RGD hexapeptides studied by 13C NMR relaxation and the model-free approach

Summary The internal mobility of three isomeric cyclic RGD hexapeptides designed to contain two -turns in defined positions, cyclo(Arg-Gly-Asp-Gly-d-Pro-Pro) (I), cyclo(Arg-Gly-Asp-d-Pro-Gly-Pro) (II) and cyclo(Arg-Gly-Asp-d-Pro-Pro-Gly) (III), have been studied by ¹³C NMR longitudinal and transverse relaxation experiments and measurements of steady-state heteronuclear {¹H}-¹³C NOE enhancement with ¹³C at natural abundance. The data were interpreted according to the model-free formalism of Lipari and Szabo, which is usually applied to data from macromolecules or larger sized peptides with overall rotational correlation times exceeding 1 ns, to yield information about internal motions on the 10–100 ps time scale. The applicability of the model-free analysis with acceptable uncertainties to these small peptides, with overall rotational correlation times slightly below 0.3 ns, was demonstrated for this specific instance. Chemical exchange contributions to T₂ from slower motions were also identified in the process. According to the order parameters obtained for its backbone -carbon atoms, II has the most rigid backbone conformation on the 10–100 ps time scale, and I the most flexible. This result coincides with the results of earlier NMR-constrained conformational searches, which indicated greatest uncertainty in the structure of I and least in II.     

Internal motions of apo-neocarzinostatin as studied by 13C NMR methine relaxation at natural abundance

Summary Dynamics of the backbone and some side chains of apo-neocarzinostatin, a 10.7 kDa carrier protein, have been studied from ¹³C relaxation rates R1, R2 and steady-state ¹³C-{¹H} NOEs, measured at natural abundance. Relaxation data were obtained for 79 nonoverlapping C resonances and for 11 threonine C single resonances. Except for three C relaxation rates, all data were analysed from a simple two-parameter spectral density function using the model-free approach of Lipari and Szabo. The corresponding C–H fragments exhibit fast (_e < 40 ps) restricted libration motions (S²=0.73 to 0.95). Global examination of the microdynamical parameters S² and _e along the amino acid sequence gives no immediate correlation with structural elements. However, different trends for the three loops involved in the binding site are revealed. The -ribbon comprising residues 37 to 47 is spatially restricted, with relatively large _e values in its hairpin region. The other -ribbon (residues 72 to 87) and the large disordered loop ranging between residues 97–107 experience small-amplitude motions on a much faster (picosecond) time scale. The two N-terminal residues, Ala¹ and Ala², and the C-terminal residue Asn¹¹³, exhibit an additional slow motion on a subnanosecond time scale (400–500 ps). Similarly, the relaxation data for eight threonine side-chain C must be interpreted in terms of a three-parameter spectral density function. They exhibit slower motions, on the nanosecond time scale (500–3000 ps). Three threonine (Thr⁶⁵, Thr⁶⁸, Thr⁸¹) side chains do not display a slow component, but an exchange contribution to the observed transverse relaxation rate R2 could not be excluded at these sites. The microdynamical parameters (S², _e and R2_ex) or (S _infslow ^sup2 , S _inffast ^sup2 and _slow) were obtained from a straightforward solution of the equations describing the relaxation data. They were calculated assuming an overall isotropic rotational correlation time _e for the protein of 5.7 ns, determined using standard procedures from R2/R1 ratios. However, it is shown that the product (1–S²)× _e is nearly independent of _e for residues not exhibiting slow motions on the nanosecond time scale. In addition, this parameter very closely follows the heteronuclear NOEs, which therefore could be good indices for local fast motions on the picosecond time scale.     

13C NMR relaxation studies of RNA base and ribose nuclei reveal a complex pattern of motions in the RNA binding site for human U1A protein

The widespread importance of induced fit and order-disorder transition in RNA recognition by proteins and small molecules makes it imperative that RNA motional properties are characterized quantitatively. Until now, however, very few studies have been dedicated to the systematic characterization of RNA motion and to their changes upon protein or small-molecule binding. The U1A protein-RNA complexes provide some of the best-studied examples of the role of RNA motional changes upon protein binding. Here, we report (13)C NMR relaxation studies of base and ribose dynamics for the RNA internal loop target of human U1A protein located within the 3'-untranslated region (3'-UTR) of the mRNA coding for U1A itself. We also report the semi-quantitative analysis of both fast (nano- to picosecond) and intermediate (micro- to millisecond) motions for this paradigmatic RNA system. We measure (13)C T(1), T(1rho) and heteronuclear nuclear Overhauser effects (NOEs) for sugar and base nuclei, as well as the power dependence of T(1rho) at 500 MHz and 750 MHz, and analyze these results using the model-free formalism. The results provide a much clearer picture of the type of motions experienced by this RNA in the absence of the protein than was provided by the analysis of the structure based solely on NOEs and scalar couplings. They define a model where the RNA internal loop region "breathes" on a micro- to millisecond timescale with respect to the double-helical regions. Superimposed on this slower motion, the residues at the very tip of the loop undergo faster (nano- to picosecond) motions. We hypothesize that these motions allow the RNA to sample multiple conformations so that the protein can select a structure within the ensemble that optimizes intermolecular contacts.     

NMR studies of dynamics in RNA and DNA by 13C relaxation

RNA and DNA molecules experience motions on a wide range of time scales, ranging from rapid localized motions to much slower collective motions of entire helical domains. The many functions of RNA in biology very often require this molecule to change its conformation in response to biological signals in the form of small molecules, proteins or other nucleic acids, whereas local motions in DNA may facilitate protein recognition and allow enzymes acting on DNA to access functional groups on the bases that would otherwise be buried in Watson-Crick base pairs. Although these statements make a compelling case to study the sequence dependent dynamics in nucleic acids, there are few residue-specific studies of nucleic acid dynamics. Fortunately, NMR studies of dynamics of nucleic acids and nucleic acids-protein complexes are gaining increased attention. The aim of this review is to provide an update of the recent progress in studies of nucleic acid dynamics by NMR based on the application of solution relaxation techniques.     

Improved labeling strategy for 13C relaxation measurements of methyl groups in proteins

Selective incorporation of ¹³C into the methyl groupsof protein side chains is described as a means for simplifying themeasurement and interpretation of ¹³C relaxation parameters.High incorporation (>90%) is accomplished by using pyruvate(3-¹³C, 99%) as the sole carbon source in the growthmedia for protein overexpression in E. coli. This improved labeling schemeincreases the sensitivity of the relaxation experiments by approximatelyfivefold when compared to randomly fractionally ¹³C-labeledprotein, allowing high-quality measurements on relatively dilute (<1 mM)protein samples at a relatively low cost.     

Nucleotide sequence of phenylalanine transfer RNA from Schizosaccharomyces pombe: implications for transfer RNA recognition by yeast phenylalanyl-tRNA synthetase.

The nucleotide sequence of Schizosaccharomyces pombe tRNAPhe was determined to be pG-U-C-G-C-A-A-U-G**-G*-U-G-psi-A-G-D-D-G-G-G-A-G-C-A-psi-G*-A-C-A-G-A-Cm-U-Gm-A-A-Y-A-psi-m5C-U-G-U-U-G-m7G-U*-C-A-U-C-G-G-T-psi-C-G-A-U-C-C-C-G-G-U-U-U-G-U-G-A-C-A-C-C-AOH. This sequence differs from that of S. cerevisiae tRNAPhe in 27 nucleotides. Saccharomyces cerevisiae phenylalanyl-tRNA synthetase aminoacylates both the homologous tRNAPhe and S. pombe t-NAPhe; the reactions have similar Km and Vmax values. However, the nucleotide sequence in the D stem is different in the two tRNAs. This region was proposed by Roe, B., et al. [(1973) Biochemistry 12, 4146--4154] to be the major recognition site for yeast phenylalanyl-tRNA synthetase, but the present results cast doubt on the validity of this hypothesis.     

Backbone dynamics of a model membrane protein: 13C NMR spectroscopy of alanine methyl groups in detergent-solubilized M13 coat protein

The filamentous coliphage M13 possesses multiple copies of a 50-residue coat protein which is inserted into the inner membrane of Escherichia coli during infection. 13C nuclear magnetic resonance (NMR) spectroscopy has been used to probe the structure and dynamics of M13 coat protein solubilized in detergent micelles. A comparison of backbone dynamics within the hydrophobic core region and the hydrophilic terminal domains was obtained by biosynthetic incorporation of [3-13C]alanine. Alanine is distributed throughout the protein and accounts for 10 residues (i.e., 20% of the total). Similar 13C NMR spectra of the protein have been obtained in two anionic detergents, sodium deoxycholate and sodium dodecyl sulfate, although the structures and physical properties of these solubilizing agents are quite different. The N-terminal alanine residues, assigned by pH titration, and the penultimate residue, assigned by carboxypeptidase A digestion, give rise to analogous peaks in both detergent systems. The pKa of Ala-1 (approximately 8.8) and the relaxation parameters of individual carbon atoms (T1, T2, and the nuclear Overhauser enhancement) are also generally similar, suggesting a similarity in the overall protein structure. Relaxation data have been analyzed according to the model-free approach of Lipari and Szabo [Lipari, G., & Szabo, A. (1982) J. Am. Chem. Soc. 104, 4546-4559]. The overall correlation times were obtained by fitting the three experimental relaxation values for a given well-resolved single carbon atom to obtain a unique value for the generalized order parameter, S2, and the effective correlation time, tau e. The former parameter reflects the spatial restriction of motion, and the latter, the rate.(ABSTRACT TRUNCATED AT 250 WORDS)     

Orientation of double-helical segments in crystals of yeast phenylalanine transfer RNA

A search of the X-ray intensities of the P2₁ crystal form of yeast transfer RNA^Phe has revealed the orientation of the double-helical segments in the crystal. Because of the ambiguity imposed by the crystal symmetry on choosing the helices belonging to the same molecule, and because of the difficulty of determining lengths and positions of helices, a unique model cannot be deduced, but only a small number of types are possible. Among the possibilities are a “boot”-shaped molecule, which may be derived from an earlier model proposed by the author, by bending out the anticodon arm, and also an L-shaped molecule. The latter is, however, not oriented in the way proposed by Kim et al. (1973) for the case of the closely related orthorhombic crystal form.     

Internal motions of transfer RNA: a study of exchanging protons by magnetic resonance

Proton exchange is a probe of macromolecular structure and kinetics. Its value is enhanced when the exchanging protons can be identified by nmr. After dilution of tRNA-H2O samples in D2O, slowly exchanging imino protons are observed, with exchange times ranging from minutes to days. In many cases they originate from the dihydro-uracil region. Most slow exchangers are sensitive to buffer catalysis. Extrapolation to infinite buffer concentration yields the life-time of the closed form, in a two-state model of each base-pair. As predicted by the model, the lifetime obtained by extrapolation is independent of the buffer. Typical lifetimes are 14 minutes for CG11 of yeast tRNAPhe at 17 degrees C, or 5 minutes for U8-A14 of yeast tRNA(Asp) at 20 degrees C, without magnesium. For most slow exchangers, magnesium increases the lifetime of the closed form, but moderately, by factors never more than five. The exchange rates of other, fast-exchanging, imino protons, as determined by line-broadening, are found to depend on buffer concentration. Base-pair lifetimes are determined as above. For instance UA6 of yeast tRNA(Phe) has a lifetime of 14 ms at 17 degrees C. Base-pairs 4 and 6 have shorter lifetimes than the rest of the acceptor stem. Imidazole is a good catalyst for proton exchange of both the long-and the short-lived base-pairs, whereas phosphate is not. Tris is efficient except for cases where, possibly, access is impeded by its size; magnesium reduces the efficiency of catalysis by tris buffer. From the variation of exchange time vs buffer concentration, one determines the buffer concentration for which the exchange rate from the open state is equal to the closing rate. Remarquably, this concentration takes comparable values for most base-pairs, whether short-lived or long-lived. Buffer effects have also been observed in poly(rA).poly(rU), for which we derive a lifetime of 2.5 ms at 27 degrees C, and in other polynucleotides. Some of the exchange times identified in the literature as base-pair lifetimes may instead reflect incomplete catalysis.     

Location of a platinum binding site in the structure of yeast phenylalanine transfer RNA

Trans-dichlorodiammineplatinum (II) reacts with yeast phenylalanine transfer RNA to yield a major platinum binding site. The tightly bound platinum has been located on the oligonucleotide Gm-A-A-Y-A-ψp containing the anticodon by standard fingerprinting methods using ³²P-labelled tRNA^Phe. This site corresponds to a single major platinum site identified during an X-ray crystallographic analysis of yeast tRNA^Phe. The solution studies have given confidence to the assignment of part of the 3 Å electron density map to the anticodon region of the molecular structure of yeast tRNA^Phe.     

Rotational diffusion tensor of nucleic acids from 13C NMR relaxation

Rotational diffusion properties have been derived for the DNA dodecamer d(CGCGAATTCGCG)₂ from ¹³C R₁ and R₁ measurements on the C₁, C₃, and C₄ carbons in samples uniformly enriched in ¹³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₁/R₁ 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₁/R₁ ratios to the solution structure is optimal for an axially symmetric rotational diffusion model, with a diffusion anisotropy, D_||/D, of 2.1±0.4, and an overall rotational correlation time, (2D_||+4D)^–1, of 3.35 ns at 35 °C in D₂O, in excellent agreement with values obtained from hydrodynamic modeling.     

Quantifying the effects of long-range 13C-13C dipolar coupling on measured relaxation rates in RNA

Journal of Biomolecular NMR - Selective stable isotope labeling has transformed structural and dynamics analysis of RNA by NMR spectroscopy. These methods can remove 13C-13C dipolar couplings that...     

Assignment of the low field proton nuclear magnetic resonance spectrum of yeast phenylalanine transfer RNA to specific base pairs

High-resolution proton nuclear magnetic resonance spectra at 220 and 300 MHz have been used to investigate the base-pairing structure of fragments of yeast tRNA^Phe, of chemically modified tRNA^Phe and of intact tRNA^Phe. To a very good approximation the positions of the fragment spectra are additive within 0·2 part per million, indicating that factors responsible for certain structural features in the intact molecule are already present in the smaller fragments (half molecules, hairpins and $\text{3}\text{4}$ molecules). A simple first-order ring-current shift theory taken in conjunction with the cloverleaf model for tRNA^Phe (RajBhandary et al., 1967) has been used to predict the low-field (? 15 to ?11 part per million) nuclear magnetic resonance spectra and make assignments of the resolved resonances to ring NH protons of specific base pairs. The general agreement between the predicted and observed spectra to within 0·2 part per million confirms in detail the cloverleaf model for the secondary structure of tRNA^Phe in solution. It is also established that ring-current shifts are the principal factor responsible for the wide range of shifts observed in the low-field spectra. As a result it is evident that the resonances are very sensitive to small changes in the secondary structure and in some cases changes in the interbase distance as small as 0·2 Å could easily be detected. It is also clear from the analysis that certain of the resonances are sensitive to the tertiary structure of the molecule and specific examples are discussed. As with our previous study, we find no evidence for any strong Watson-Crick type base pairs beyond those predicted by the cloverleaf structure.     

Nuclear magnetic resonance signal assignments of purified [13C]methyl-enriched yeast phenylalanine transfer ribonucleic acid

Yeast tRNA Phe, enriched in carbon-13 specifically at the naturally occurring methyl groups, has been produced through biosynthesis, then purified, and analyzed. Transfer RNA Phe was purified from the [13C]methyl-enriched, unfractionated tRNA that had been extracted from a methionine auxotroph of Saccharomyces cerevisiae [Agris, P. F., Kovacs, S. A. H., Smith, C., Kopper, R. H., & Schmidt, P. G. (1983) Biochemistry 22, 1402-1408]. The yeast had been grown in minimal medium supplemented with [13C]methylmethionine. Transfer RNA Phe purity and the full extent of nucleoside modification were confirmed by high-performance liquid chromatography of constituent nucleosides with simultaneous UV spectral identification and quantitation. Mass spectometry of [13C]methyl-enriched nucleosides and NMR of the tRNA indicated an enrichment of at least 70 atom %. Twelve resolved and prominent carbon-13 NMR signals from the tRNA were seen between 10 and 60 ppm. These have been assigned to 13 of the 14 naturally occurring methyl groups. However, the partially resolved signals assigned to the two 5-methylcytidines could not be assigned to their specific nucleoside positions of either 40 or 49 in the molecule. In addition, the partially resolved signals of the two methyl esters of wybutosine could not be distinguished. The methyl group found not to be enriched with 13C is bound to the ring carbon in the hypermodified nucleoside wybutosine (Y). A 13th enriched signal downfield (120.9 ppm) has been assigned to one of the two carbons added to guanosine to form the third ring in the biosynthesis of Y. The 13C enrichment of this ring carbon demonstrates its origin from the methionine methyl group.(ABSTRACT TRUNCATED AT 250 WORDS)     

Internal motional amplitudes and correlated bond rotations in an alpha-helical peptide derived from 13C and 15N NMR relaxation

Peptide GFSKAELAKARAAKRGGY folds in an alpha-helical conformation that is stabilized by formation of a hydrophobic staple motif and an N-terminal capping box (Munoz V. Blanco FJ, Serrano L, 1995, Struct Biol 2:380-385). To investigate backbone and side-chain internal motions within the helix and hydrophobic staple, residues F2, A5, L7, A8, and A10 were selectively 13C- and 15N-enriched and NMR relaxation experiments were performed in water and in water/trifluoroethanol (TFE) solution at four Larmor frequencies (62.5, 125, 150, and 200 MHz for 13C). Relaxation data were analyzed using the model free approach and an anisotropic diffusion model. In water, angular variances of motional vectors range from 10 to 20 degrees and backbone phi,psi bond rotations for helix residues A5, L7, A8, and A10 are correlated indicating the presence of Calpha-H, Calpha-Cbeta, and N-H rocking-type motions along the helix dipole axis. L7 side-chain CbetaH2 and CgammaH motions are also correlated and as motionally restricted as backbone CalphaH, suggesting considerable steric hindrance with neighboring groups. In TFE which stabilizes the fold, internal motional amplitudes are attenuated and rotational correlations are increased. For the side chain of hydrophobic staple residue F2, wobbling-in-a-cone type motions dominate in water, whereas in TFE, the Cbeta-Cgamma bond and phenyl ring fluctuate more simply about the Calpha-Cbeta bond. These data support the Daragan-Mayo model of correlated bond rotations (Daragan VA, Mayo KH, 1996, J Phys Chem 100:8378-8388) and contribute to a general understanding of internal motions in peptides and proteins.     

Synthesis of methyl 4'-O-methyl-beta-D-cellobioside-13C12 from D-glucose-13C6. Part 2: solid-state NMR studies

Double Quantum (DQ) NMR, which utilizes the magnetic dipole interaction between the (13)C atoms, was used for the complete assignment of the (13)C NMR resonances to the corresponding carbon ring positions for the monoclinic and triclinic allomorphs of methyl 4'-O-methyl-beta-D-cellobioside-(13)C(12)(1-(13)C(12)), a cellodextrin model compound of cellulose (13)C-perlabeled at the cellobiose core. The through-space interactions were used to identify the direct chemical bonds between adjacent carbon atoms in the rings. More importantly, the (13)C NMR signals of the carbon sites C1' and C4 involved in the glycosidic bond were identified. This allowed for the complete (13)C chemical shift assignment, that when combined with the X-ray crystallography data provides a complete characterization.