Two folded conformers of ubiquitin revealed by high-pressure NMR.

High-pressure 15N/1H two-dimensional NMR spectroscopy has been utilized to study conformational fluctuation of a 76-residue protein ubiquitin at pH 4.5 at 20 degrees C. The on-line variable pressure cell technique is used in conjunction with a high-field NMR spectrometer operating at 750 MHz for 1H in the pressure range between 30 and 3500 bar. Large, continuous and reversible pressure-induced 1H and 15N chemical shifts were observed for 68 backbone amide groups, including the 7.52 ppm 15N shift of Val70 at 3500 bar, indicating a large-scale conformational change of ubiquitin with pressure. On the basis of the analysis of sigmoid-shaped pressure shifts, we conclude that ubiquitin exists as an equilibrium mixture of two major folded conformers mutually converting at a rate exceeding approximately 10(4) s(-1) at 20 degrees C at 2000 bar. The second conformer exists at a population of approximately 15% (DeltaG(0) = 4.2 kJ/mol) and is characterized with a significantly smaller partial molar volume (DeltaV(0) = -24 mL/mol) than that of the well-known basic native conformer. The analysis of 1H and 15N pressure shifts of individual amide groups indicates that the second conformer has a loosened core structure with weakened hydrogen bonds in the five-stranded beta-sheet. Furthermore, hydrogen bonds of residues 67-72 belonging to beta5 are substantially weakened or partially broken, giving increased freedom of motion for the C-terminal segment. The latter is confirmed by the significant decrease in 15N[1H] nuclear Overhauser effect for residues beyond 70 at high pressure. Since the C-terminal carboxyl group constitutes the reactive site for producing a multi-ubiquitin structure, the finding of the second folded conformer with a substantially altered conformation and mobility in the C-terminal region will shed new light on the reaction mechanism of ubiquitin.     

The geometry of alpha-sheet: Implications for its possible function as amyloid precursor in proteins

alpha-sheet has been proposed as the main constituent of the prefibrillar intermediate during amyloid formation. Here the helical parameters of the alpha-sheet strand are calculated from average main-chain dihedral angles reported from molecular dynamics simulations. It is an almost linear polypeptide that forms a right-handed helix of about 100 A diameter, with 100 residues and a rise of 30 A per turn. The strands are curved but untwisted, which implies that neighboring strands need not coil to make interstrand hydrogen bonds. This suggests that compared to beta-sheets in native folded proteins, alpha-sheets can be larger and stack more easily to create extensive 3D blocks. It is shown that alpha-sheet is related to a category of structures termed "mirror" structures. Mirror structures have repetitive pairs of main-chain dihedral angles at residues i and i+1 that satisfy the condition phi(i) (+1) = -psi(i), psi(i) (+1) = -phi(i). They are uniquely identified by the two orientations of their peptide planes, specified by phi(i) and psi(i). Their side chains point alternately in opposite directions. Interestingly, their conformations are insensitive to phi(i) and psi(i) in that the pseudo dihedral angle formed by four consecutive C(alpha) atoms is always close to 180 degrees . There are two types: "beta-mirror" and "alpha-mirror" structure; beta-mirror structures relate to beta-sheet by small peptide plane rotations, of less than 90 degrees , while alpha-mirror structures are close to alpha-sheet and relate to beta-sheet by approximately 180 degrees peptide plane flips. Most mirror structures, and in particular the alpha-mirror, form wide helices with diameters 50-70 A. Their gentle curvature, and therefore that of the alpha-sheet, arises from the orientation of successive peptide units causing the difference in the bond angles at the C and N atoms of the peptide unit to gradually change the direction of the chain.     

Dihydrofolate reductase: sequential resonance assignments using 2D and 3D NMR and secondary structure determination in solution

Three-dimensional (3D) heteronuclear NMR techniques have been used to make sequential 1H and 15N resonance assignments for most of the residues of Lactobacillus casei dihydrofolate reductase (DHFR), a monomeric protein of molecular mass 18,300 Da. A uniformly 15N-labeled sample of the protein was prepared and its complex with methotrexate (MTX) studied by 3D 15N/1H nuclear Overhauser-heteronuclear multiple quantum coherence (NOESY-HMQC), Hartmann-Hahn-heteronuclear multiple quantum coherence (HOHAHA-HMQC), and HMQC-NOESY-HMQC experiments. These experiments overcame most of the spectral overlap problems caused by chemical shift degeneracies in 2D spectra and allowed the 1H-1H through-space and through-bond connectivities to be identified unambiguously, leading to the resonance assignments. The novel HMQC-NOESY-HMQC experiment allows NOE cross peaks to be detected between NH protons even when their 1H chemical shifts are degenerate as long as the amide 15N chemical shifts are nondegenerate. The 3D experiments, in combination with conventional 2D NOESY, COSY, and HOHAHA experiments on unlabelled and selectively deuterated DHFR, provide backbone assignments for 146 of the 162 residues and side-chain assignments for 104 residues of the protein. Data from the NOE-based experiments and identification of the slowly exchanging amide protons provide detailed information about the secondary structure of the binary complex of the protein with methotrexate. Sequential NHi-NHi+1 NOEs define four regions with helical structure. Two of these regions, residues 44-49 and 79-89, correspond to within one amino acid to helices C and E in the crystal structure of the DHFR.methotrexate.NADPH complex [Bolin et al. (1982) J. Biol. Chem. 257, 13650-13662], while the NMR-determined helix formed by residues 26-35 is about one turn shorter at the N-terminus than helix B in the crystal structure, which spans residues 23-34. Similarly, the NMR-determined helical region comprising residues 102-110 is somewhat offset from the crystal structure's helix F, which encompasses residues 97-107. Regions of beta-sheet structure were characterized in the binary complex by strong alpha CHi-NHi+1 NOEs and by slowly exchanging amide protons. In addition, several long-range NOEs were identified linking together these stretches to form a beta-sheet. These elements align perfectly with corresponding elements in the crystal structure of the DHFR.methotrexate.NADPH complex, which contains an eight-stranded beta-sheet, indicating that the main body of the beta-sheet is preserved in the binary complex in solution.     

Structures of N-termini of helices in proteins.

We have surveyed 393 N-termini of alpha-helices and 156 N-termini of 3(10)-helices in 85 high resolution, non-homologous protein crystal structures for N-cap side-chain rotamer preferences, hydrogen bonding patterns, and solvent accessibilities. We find very strong rotamer preferences that are unique to N-cap sites. The following rules are generally observed for N-capping in alpha-helices: Thr and Ser N-cap side chains adopt the gauche - rotamer, hydrogen bond to the N3 NH and have psi restricted to 164 +/- 8 degrees. Asp and Asn N-cap side chains either adopt the gauche - rotamer and hydrogen bond to the N3 NH with psi = 172 +/- 10 degrees, or adopt the trans rotamer and hydrogen bond to both the N2 and N3 NH groups with psi = 1-7 +/- 19 degrees. With all other N-caps, the side chain is found in the gauche + rotamer so that the side chain does not interact unfavorably with the N-terminus by blocking solvation and psi is unrestricted. An i, i + 3 hydrogen bond from N3 NH to the N-cap backbone C = O in more likely to form at the N-terminus when an unfavorable N-cap is present. In the 3(10)-helix Asn and Asp remain favorable N-caps as they can hydrogen bond to the N2 NH while in the trans rotamer; in contrast, Ser and Thr are disfavored as their preferred hydrogen bonding partner (N3 NH) is inaccessible. This suggests that Ser is the optimum choice of N-cap when alpha-helix formation is to be encouraged while 3(10)-helix formation discouraged. The strong energetic and structural preferences found for N-caps, which differ greatly from positions within helix interiors, suggest that N-caps should be treated explicitly in any consideration of helical structure in peptides or proteins.     

Dynamic aspects of extracellular loop region as a proton release pathway of bacteriorhodopsin studied by relaxation time measurements by solid state NMR

Local dynamics of interhelical loops in bacteriorhodopsin (bR), the extracellular BC, DE and FG, and cytoplasmic AB and CD loops, and helix B were determined on the basis of a variety of relaxation parameters for the resolved 13C and 15N signals of [1-13C]Tyr-, [15N]Pro- and [1-13C]Val-, [15N]Pro-labeled bR. Rotational echo double resonance (REDOR) filter experiments were used to assign [1-13C]Val-, [15N]Pro signals to the specific residues in bR. The previous assignments of [1-13C]Val-labeled peaks, 172.9 or 171.1 ppm, to Val69 were revised: the assignment of peak, 172.1 ppm, to Val69 was made in view of the additional information of conformation-dependent 15N chemical shifts of Pro bonded to Val in the presence of 13C-15N correlation, although no assignment of peak is feasible for 13C nuclei not bonded to Pro. 13C or 15N spin-lattice relaxation times (T1), spin-spin relaxation times under the condition of CP-MAS (T2), and cross relaxation times (TCH and TNH) for 13C and 15N nuclei and carbon or nitrogen-resolved, 1H spin-lattice relaxation times in the rotating flame (1H T1 rho) for the assigned signals were measured in [1-13C]Val-, [15N]Pro-bR. It turned out that V69-P70 in the BC loop in the extracellular side has a rigid beta-sheet in spite of longer loop and possesses large amplitude motions as revealed from 13C and 15N conformation-dependent chemical shifts and T1, T2, 1H T1 rho and cross relaxation times. In addition, breakage of the beta-sheet structure in the BC loop was seen in bacterio-opsin (bO) in the absence of retinal.     

15N and 1H NMR study of histidine containing protein (HPr) from Staphylococcus carnosus at high pressure

The pressure-induced changes in 15N enriched HPr from Staphylococcus carnosus were investigated by two-dimensional (2D) heteronuclear NMR spectroscopy at pressures ranging from atmospheric pressure up to 200 MPa. The NMR experiments allowed the simultaneous observation of the backbone and side-chain amide protons and nitrogens. Most of the resonances shift downfield with increasing pressure indicating generalized pressure-induced conformational changes. The average pressure-induced shifts for amide protons and nitrogens are 0.285 ppm GPa(-1) at 278 K and 2.20 ppm GPa(-1), respectively. At 298 K the corresponding values are 0.275 and 2.41 ppm GPa(-1). Proton and nitrogen pressure coefficients show a significant but rather small correlation (0.31) if determined for all amide resonances. When restricting the analysis to amide groups in the beta-pleated sheet, the correlation between these coefficients is with 0.59 significantly higher. As already described for other proteins, the amide proton pressure coefficients are strongly correlated to the corresponding hydrogen bond distances, and thus are indicators for the pressure-induced changes of the hydrogen bond lengths. The nitrogen shift changes appear to sense other physical phenomena such as changes of the local backbone conformation as well. Interpretation of the pressure-induced shifts in terms of structural changes in the HPr protein suggests the following picture: the four-stranded beta-pleated sheet of HPr protein is the least compressible part of the structure showing only small pressure effects. The two long helices a and c show intermediary effects that could be explained by a higher compressibility and a concomitant bending of the helices. The largest pressure coefficients are found in the active center region around His15 and in the regulatory helix b which includes the phosphorylation site Ser46 for the HPr kinase. This suggests that this part of the structure occurs in a number of different structural states whose equilibrium populations are shifted by pressure. In contrast to the surrounding residues of the active center loop that show large pressure effects, Ile14 has a very small proton and nitrogen pressure coefficient. It could represent some kind of anchoring point of the active center loop that holds it in the right place in space, whereas other parts of the loop adapt themselves to changing external conditions.     

Structural and dynamic characterization of an unfolded state of poplar apo-plastocyanin formed under nondenaturing conditions

Plastocyanin is a predominantly beta-sheet protein containing a type I copper center. The conformational ensemble of a denatured state of apo-plastocyanin formed in solution under conditions of low salt and neutral pH has been investigated by multidimensional heteronuclear NMR spectroscopy. Chemical shift assignments were obtained by using three-dimensional triple-resonance NMR experiments to trace through-bond heteronuclear connectivities along the backbone and side chains. The (3)J(HN,Halpha) coupling constants, (15)N-edited proton-proton nuclear Overhauser effects (NOEs), and (15)N relaxation parameters were also measured for the purpose of structural and dynamic characterization. Most of the residues corresponding to beta-strands in the folded protein exhibit small upfield shifts of the (13)C(alpha) and (13)CO resonances relative to random coil values, suggesting a slight preference for backbone dihedral angles in the beta region of (phi,psi) space. This is further supported by the presence of strong sequential d(alphaN)(i, i + 1) NOEs throughout the sequence. The few d(NN)(i, i + 1) proton NOEs that are observed are mostly in regions that form loops in the native plastocyanin structure. No medium or long-range NOEs were observed. A short sequence, between residues 59 and 63, was found to populate a nonnative helical conformation in the unfolded state, as indicated by the shift of the (13)C(alpha), (13)CO, and (1)H(alpha) resonances relative to random coil values and by the decreased values of the (3)J(HN,Halpha) coupling constants. The (15)N relaxation parameters indicate restriction of motions on a nanosecond timescale in this region. Intriguingly, this helical conformation is present in a sequence that is close to but not in the same location as the single short helix in the native folded protein. The results are consistent with earlier NMR studies of peptide fragments of plastocyanin and confirm that the regions of the sequence that form beta-strands in the native protein spontaneously populate the beta-region of (phi,psi) space under folding conditions, even in the absence of stabilizing tertiary interactions. We conclude that the state of apo-plastocyanin present under nondenaturing conditions is a noncompact unfolded state with some evidence of nativelike and nonnative local structuring that may be initiation sites for folding of the protein.     

Response of native and denatured hen lysozyme to high pressure studied by (15)N/(1)H NMR spectroscopy.

High-pressure (15)N/(1)H NMR techniques were used to characterize the conformational fluctuations of hen lysozyme, in its native state and when denatured in 8 M urea, over the pressure range 30--2000 bar. Most (1)H and (15)N signals of native lysozyme show reversible shifts to low field with increasing pressure, the average pressure shifts being 0.069 +/- 0.101 p.p.m. ((1)H) and 0.51 +/- 0.36 p.p.m. ((15)N). The shifts indicate that the hydrogen bonds formed to carbonyl groups or water molecules by the backbone amides are, on average, shortened by approximately 0.02 A as a result of pressure. In native lysozyme, six residues in the beta domain or at the alpha/beta domain interface have anomalously large nonlinear (15)N and (1)H chemical-shift changes. All these residues lie close to water-containing cavities, suggesting that there are conformational changes involving these cavities, or the water molecules within them, at high pressure. The pressure-induced (1)H and (15)N shifts for lysozyme denatured in 8 M urea are much more uniform than those for native lysozyme, with average backbone amide shifts of 0.081 +/- 0.029 p.p.m. ((1)H) and 0.57 +/- 0.14 p.p.m. ((15)N). The results show that overall there are no significant variations in the local conformational properties of denatured lysozyme with pressure, although larger shifts in the vicinity of a persistent hydrophobic cluster indicate that interactions in this part of the sequence may rearrange. NMR diffusion measurements demonstrate that the effective hydrodynamic radius of denatured lysozyme, and hence the global properties of the denatured ensemble, do not change detectably at high pressure.     

The solution structure of bovine pancreatic trypsin inhibitor at high pressure

The solution structure of bovine pancreatic trypsin inhibitor (BPTI) at a pressure of 2 kbar is presented. The structure was calculated as a change from an energy-minimized low-pressure structure, using (1)H chemical shifts as restraints. The structure has changed by 0.24 A RMS, and has almost unchanged volume. The largest changes as a result of pressure are in the loop 10-16, which contains the active site of BPTI, and residues 38-42, which are adjacent to buried water molecules. Hydrogen bonds are compressed by 0.029 +/- 0.117 A, with the longer hydrogen bonds, including those to internal buried water molecules, being compressed more. The hydrophobic core is also compressed, largely from reduction of packing defects. The parts of the structure that have the greatest change are close to buried water molecules, thus highlighting the importance of water molecules as the nucleation sites for volume fluctuation of proteins in native conditions.     

High pressure NMR reveals a variety of fluctuating conformers in beta-lactoglobulin

High pressure 1H/15N two-dimensional NMR spectroscopy has been used to study conformational fluctuation in bovine beta-lactoglobulin at pH 2.0 and 36 degrees C. Pressure dependencies of 1H and 15N chemical shifts and cross-peak intensities were analyzed at more than 80 independent atom sites between 30 and 2000 bar. Unusually large and non-linear chemical shift pressure dependencies are found for residues centering in the hydrophobic core region, suggesting the existence of low-lying excited native states (N') of the protein. Measurement of 1H/15N cross-peak intensities at individual amide sites as a function of pressure suggests that unfolding events occur independently in two sides of the beta-barrel, i.e. the hydrophobic core side (betaF-H) (producing I2) and the non-core side (betaB-E) (producing I1). At 1 bar the stability is higher for the core region (DeltaG0 = 6.5(+/-2.0) kcal/mol) than for the non-core region (4.6(+/-1.3) kcal/mol), but at high pressure the stability is reversed due to a larger DeltaV value of unfolding for the core region (90.0(+/-35.2) ml/mol) than that for the non-core region (57.4(+/-14.4) ml/mol), possibly due to an uneven distribution of cavities. The DeltaG0 profile along the amino acid sequence obtained from the pressure experiment is found to coincide well with that estimated from hydrogen exchange experiments. Altogether, the high pressure NMR experiment has revealed a variety of fluctuating conformers of beta-lactoglobulin, notably N, N', I1, I2 and the totally unfolded conformer U. Fluctuation of N to I1 and I2 conformers with open barrel structures could be a common design of lipocalin family proteins which bind various hydrophobic compounds in its barrel structure.     

External reflection FTIR of peptide monolayer films in situ at the air/water interface: experimental design, spectra-structure correlations, and effects of hydrogen-deuterium exchange.

A Fourier transform infrared spectrometer has been interfaced with a surface balance and a new external reflection infrared sampling accessory, which permits the acquisition of spectra from protein monolayers in situ at the air/water interface. The accessory, a sample shuttle that permits the collection of spectra in alternating fashion from sample and background troughs, reduces interference from water vapor rotation-vibration bands in the amide I and amide II regions of protein spectra (1520-1690 cm-1) by nearly an order of magnitude. Residual interference from water vapor absorbance ranges from 50 to 200 microabsorbance units. The performance of the device is demonstrated through spectra of synthetic peptides designed to adopt alpha-helical, antiparallel beta-sheet, mixed beta-sheet/beta-turn, and unordered conformations at the air/water interface. The extent of exchange on the surface can be monitored from the relative intensities of the amide II and amide I modes. Hydrogen-deuterium exchange may lower the amide I frequency by as much as 11-12 cm-1 for helical secondary structures. This shifts the vibrational mode into a region normally associated with unordered structures and leads to uncertainties in the application of algorithms commonly used for determination of secondary structure from amide I contours of proteins in D2O solution.     

Folding of peptide fragments comprising the complete sequence of proteins. Models for initiation of protein folding. II. Plastocyanin.

In an attempt to understand the earliest events in the protein folding pathway, the complete sequence of French bean plastocyanin has been synthesized as a series of short peptide fragments, and the conformational preferences of each peptide examined in aqueous solution using proton n.m.r. methods. Plastocyanin consists largely of beta-sheet, with reverse turns and loops between the strands of the sheet, and one short helix. The n.m.r. experiments indicate that most of the peptides derived from the plastocyanin sequence have remarkably little propensity to adopt folded conformations in aqueous solution, in marked contrast to the peptides derived from the helical protein, myohemerythrin (accompanying paper). For most plastocyanin peptides, the backbone dihedral angles are predominantly in the beta-region of conformational space. Some of the peptides show weak NOE connectivities between adjacent amide protons, indicative of small local populations of backbone conformations in the a region of (phi,psi) space. A conformational preference for a reverse turn is seen in the sequence Ala65-Pro-Gly-Glu68, where a turn structure is found in the folded protein. Significantly, the peptide sequences that populate the alpha-region of (phi,psi) space are mostly derived from turn and loop regions in the protein. The addition of trifluoroethanol does not drive the peptides into helical conformations. In one region of the sequence, the n.m.r. spectra provide evidence of the formation of a hydrophobic cluster involving aromatic and aliphatic side-chains. These results have significance for understanding the initiation of protein folding. From these studies of the fragments of plastocyanin (this paper) and myohemerythrin (accompanying paper), it appears that there is a pre-partitioning of the conformational space sampled by the polypeptide backbone that is related to the secondary structure in the final folded state.     

The nature of the turn in omega loops of proteins

An analysis of Omega loops in a nonredundant set of protein structures from the Protein Data Bank has been carried out to determine the nature of the "turn elements" present. Because Omega loops essentially reverse their direction in three-dimensional space, this analysis was made with respect to four turn elements identified as (1) Gly; (2) Pro; (3) a residue with alpha-helical phi,psi angles, termed a helical residue; and (4) a cis peptide. A set of 1079 Omega loops from a set of 680 proteins were used for the analysis. Apart from other criteria that define Omega loops, the selection of an Omega loop from a cluster of loops is based on an exposure index. In this study, analyses have been made with two sets of data: (1) Omega loops arising from a minimum exposure index indicative of a less exposed loop (xmin set) and (2) Omega loops with a maximum exposure index indicative of a relatively exposed loop (xmax set). Overall residue preferences and positional preferences have been examined. Positions of the turn elements for Omega loops of varying length have also been studied. Specific positional preferences are observed for particular turn elements with regard to the length of Omega loops. Analysis in terms of the turn elements can provide guidelines for modeling of loops in proteins. Apart from Pro, which has the natural tendency to form cis peptide bonds, a higher occurrence of non-Pro cis peptide bonds is observed. Torsion angles in Omega loops also indicate the occurrence of a large number of residues with helical phi,psi angles, necessary for the turn in the loop structures.     

Thermodynamic and structural characterization of Asn and Ala residues in the disallowed II' region of the Ramachandran plot

Residue Asn47 at position L1 of a type II' beta-turn of the alpha-spectrin SH3 domain is located in a disallowed region of the Ramachandran plot (phi = 56 +/- 12, psi = -118 +/- 17). Therefore, it is expected that replacement of Asn47 by Gly should result in a considerable stabilization of the protein. Thermodynamic analysis of the N47G and N47A mutants shows that the change in free energy is small (approximately 0.7 kcal/mol; approximately 3 kJ/mol) and comparable to that found when mutating a Gly to Ala in a alpha-helix or beta-sheet. X-ray structural analysis of these mutants shows that the conformation of the beta-turn does not change upon mutation and, therefore, that there is no relaxation of the structure, nor is there any gain or loss of interactions that could explain the small energy change. Our results indicate that the energetic definition of II' region of the Ramachandran plot (phi = 60 +/- 30, psi = -115 +/- 15) should be revised for at least Ala and Asn in structure validation and protein design.     

Protective effect of lipidic surfaces against pressure-induced conformational changes of poly(L-lysine)

Poly(L-lysine) bound to phosphatidylglycerol or phosphatidic acid bilayers was submitted to hydrostatic pressure in a diamond anvil cell to investigate whether the lipidic surfaces can protect the polypeptide against pressure-induced conformational transformations. The amide I region of the infrared spectrum of dimyristoylphosphatidic acid bound polylysine shows that most of the polypeptide retains its beta-sheet structure up to 19 kbar, while it is known to convert entirely to alpha-helix at approximately 2 kbar in the absence of the lipid [Carrier, D., Mantsch, H.H., & Wong, P.T.T. (1989) Biopolymers (in press)]. The simultaneous binding of the polypeptidic molecules to two opposing bilayers appears to be required in order to preserve the beta-sheet structure at pressures over approximately 9 kbar: a small proportion of the polypeptide, most likely the molecules at the surface of the aggregated bilayers, was found to convert to unordered and eventually to alpha-helical conformations in the pressure range 9-19 kbar. The decrease from 1612 to 1606 cm-1 of the frequency of the major beta-sheet component of the infrared amide I band as the pressure is raised to 6 kbar indicates a strengthening of the interchain hydrogen bonds. The high-pressure infrared spectra of polylysine bound to dimyristoyl- and dipalmitoylphosphatidylglycerol show that the polypeptide remains alpha-helical up to approximately 12 kbar, though the changes in the bandshape indicate an increase in hydrogen bond strength. The formation of a small amount of beta-sheet was observed during decompression and is attributed to the effect of dehydration on the polypeptidic molecules located at the surface of the aggregates.(ABSTRACT TRUNCATED AT 250 WORDS)     

15N-labeled Escherichia coli tRNAfMet, tRNAGlu, tRNATyr, and tRNAPhe. Double resonance and two-dimensional NMR of N1-labeled pseudouridine

The N1 imino units in Escherichia coli tRNAfMet, tRNAGlu, tRNAPhe, and tRNATyr were studied by 1H-15N NMR using three different techniques to suppress signals of protons not attached to 15N. Two of the procedures, Fourier internuclear difference spectroscopy and two-dimensional forbidden echo spectroscopy permitted 1H and 15N chemical shifts to be measured simultaneously at 1H sensitivity. The tRNAs were labeled by fermentation of the uracil auxotroph S phi 187 on a minimal medium containing [1-15N]uracil. 1H and 15N resonances were detected for all of the N1 psi imino units except psi 13 at the end of the dihydrouridine stem in tRNAGlu. Chemical shifts for imino units in the tRNAs were compared with "intrinsic" values in model systems. The comparisons show that the A X psi pairs at the base of the anticodon stem in E. coli tRNAPhe and tRNATyr have psi in an anti conformation. The N1 protons of psi in other locations, including psi 32 in the anticodon loop of tRNAPhe, form internal hydrogen bonds to bridging water molecules or 2'-hydroxyl groups in nearby ribose units. These interactions permit psi to stabilize the tertiary structure of a tRNA beyond what is provided by the U it replaces.     

Conformational investigation of alpha,beta-dehydropeptides. XV: N-acetyl-alpha,beta-dehydroamino acid N 'N '-dimethylamides: conformational properties from infrared and theoretical studies.

The FTIR spectra were analysed in the region of the nu(s)(N-H), AI(C=O) and nu(s)(Calpha=Cbeta) bands for a series of Ac-DeltaXaa-NMe2, where DeltaXaa = DeltaAla, (Z)-DeltaAbu, (Z)-DeltaLeu, (Z)-DeltaPhe and DeltaVal, to determine a predominant solution conformation of these alpha,beta-dehydropeptide-related molecules. Measurements were taken in CCl4, DCM and MeCN solutions. In the same way, spectra of saturated analogues Ac-Xaa-NMe2, where Xaa = Ala, Abu, Leu, Phe and Val, were investigated. To help interpret the spectroscopic results, conformational maps were calculated by the B3LYP/6-31+G** method. Also, the relative energies of all conformers of the dehydro compounds in vacuo as well as in the studied solvents in addition to the theoretical IR frequencies of these conformers were calculated. For comparison, molecules of two saturated analogues, Ac-L-Ala-NMe2 and Ac-L-Phe-NMe2, were calculated in a similar way. Both unsaturated and saturated compounds, which have an aliphatic side chain, occur in CCl4 and DCM mainly as a mixture of extended conformers with the C5 H-bond and open conformers. As solvent polarity increases, participation of the open conformers also increases, and in MeCN, the model amides are almost exclusively in the open form, except Ac-DeltaAla-NMe2, which shows a small amount of the H-bonded conformer. Ac-DeltaAla-NMe2 and Ac-DeltaAbu-NMe2 have stronger C5 hydrogen bonds than those of their saturated counterparts. As the calculations indicate, the open conformation of the unsaturated amides is conformer H/F with phi, psi -44 +/- 5 degrees, 127 +/- 4 degrees. This is the second lowest in energy conformer in vacuo and in CCl4 and the lowest one in more polar solvents. The open conformation of Ac-L-Ala-NMe2 constitutes conformer C with phi, psi -101.5 degrees, 112.7 degrees. For Ac-DeltaAla-NMe2 and Ac-DeltaAbu-NMe2, FTIR also reveals the presence of a third conformer. Calculations indicate that is the semiextended conformer D with the N1-H1...N2 hydrogen bond/contact. In all solvents, Ac-L-Phe-NMe2 and Ac-(Z)-DeltaPhe-NMe2 show only the extended E and the open H/F, respectively. In both there is an amide/pi(Ph) interaction.     

Backbone dynamics of calcium-loaded calbindin D9k studied by two-dimensional proton-detected 15N NMR spectroscopy.

Backbone dynamics of calcium-loaded calbindin D9k have been investigated by two-dimensional proton-detected heteronuclear nuclear magnetic resonance spectroscopy, using a uniformly 15N enriched protein sample. Spin-lattice relaxation rate constants, spin-spin relaxation rate constants, and steady-state [1H]-15N nuclear Overhauser effects were determined for 71 of the 72 backbone amide 15N nuclei. The relaxation parameters were analyzed using a model-free formalism that incorporates the overall rotational correlation time of the molecule, and a generalized order parameter (S2) and an effective internal correlation time for each amide group. Calbindin D9k contains two helix-loop-helix motifs joined by a linker loop at one end of the protein and a beta-type interaction between the two calcium-binding loops at the other end. The amplitude of motions for the calcium-binding loops and the helices are similar, as judged from the average S2 values of 0.83 +/- 0.05 and 0.85 +/- 0.04, respectively. The linker region joining the two calcium-binding subdomains of the molecule has a significantly higher flexibility, as indicated by a substantially lower average S2 value of 0.59 +/- 0.23. For residues in the linker loop and at the C-terminus, the order parameter is further decomposed into separate order parameters for motional processes on two distinct time scales. The effective correlation times are significantly longer for helices I and IV than for helices II and III or for the calcium-binding loops. Residue by residue comparisons reveal correlations of the order parameters with both the crystallographic B-factors and amide proton exchange rates, despite vast differences in the time scales to which these properties are sensitive. The order parameters are also utilized to distinguish regions of the NMR-derived three-dimensional structure of calbindin D9k that are poorly defined due to inherently high flexibility, from poorly defined regions with average flexibility but a low density of structural constraints.     

Predicting 15N chemical shifts in proteins using the preceding residue-specific individual shielding surfaces from φ, ψi−1, and χ1torsion angles

Empirical shielding surfaces are most commonly used to predict chemical shifts in proteins from known backbone torsion angles, phi and psi. However, the prediction of (15)N chemical shifts using this technique is significantly poorer, compared to that for the other nuclei such as (1)H(alpha), (13)C(alpha), and (13)C(beta). In this study, we investigated the effects from the preceding residue and the side-chain geometry, chi(1), on (15)N chemical shifts by statistical methods. For an amino acid sequence XY, the (15)N chemical shift of Y is expressed as a function of the amino acid types of X and Y, as well as the backbone torsion angles, phi and psi(i-1). Accordingly, 380 empirical 'Preceding Residue Specific Individual (PRSI)' (15)N chemical shift shielding surfaces, representing all the combinations of X and Y (except for Y=Pro), were built and used to predict (15)N chemical shift from phi and psi(i-1). We further investigated the chi(1) effects, which were found to account for differences in (15)N chemical shifts by approximately 5 ppm for amino acids Val, Ile, Thr, Phe, His, Tyr, and Trp. Taking the chi(1) effects into account, the chi(1)-calibrated PRSI shielding surfaces (XPRSI) were built and used to predict (15)N chemical shifts for these amino acids. We demonstrated that (15)N chemical shift predictions are significantly improved by incorporating the preceding residue and chi(1) effects. The present PRSI and XPRSI shielding surfaces were extensively compared with three recently published programs, SHIFTX (Neal et al., 2003), SHIFTS (Xu and Case, 2001 and 2002), and PROSHIFT (Meiler, 2003) on a set of ten randomly selected proteins. A set of Java programs using XPRSI shielding surfaces to predict (15)N chemical shifts in proteins were developed and are freely available for academic users at http://www.pronmr.com or by sending email to one of the authors Yunjun Wang (yunjunwang@yahoo.com).