Determination of a high precision structure of a novel protein, Linum usitatissimum trypsin inhibitor (LUTI), using computer-aided assignment of NOESY cross-peaks

The solution structure of a novel 69 residue proteinase inhibitor, Linum usitatissimum trypsin inhibitor (LUTI), was determined using a method based on computer aided assignment of nuclear Overhauser enhancement spectroscopy (NOESY) data. The approach applied uses the program NOAH/DYANA for automatic assignment of NOESY cross-peaks. Calculations were carried out using two unassigned NOESY peak lists and a set of determined dihedral angle restraints. In addition, hydrogen bonds involving amide protons were identified during calculations using geometrical criteria and values of HN temperature coefficients. Stereospecific assignment of beta-methylene protons was carried out using a standard procedure based on nuclear Overhauser enhancement intensities and 3J(alpha)(beta) coupling constants. Further stereospecific assignment of methylene protons and diastereotopic methyl groups were established upon structure-based method available in the program GLOMSA and chemical shift calculations. The applied algorithm allowed us to assign 1968 out of 2164 peaks (91%) derived from NOESY spectra recorded in H2O and 2H2O. The final experimental data input consisted of 1609 interproton distance restraints, 88 restraints for 44 hydrogen bonds, 63 torsion angle restraints and 32 stereospecifically assigned methylene proton pairs and methyl groups. The algorithm allowed the calculation of a high precision protein structure without the laborious manual assignment of NOESY cross-peaks. For the 20 best conformers selected out of 40 refined ones in the program CNS, the calculated average pairwise rmsd values for residues 3 to 69 were 0.38 A (backbone atoms) and 1.02 A (all heavy atoms). The three-dimensional LUTI structure consists of a mixed parallel and antiparallel beta-sheet, a single alpha-helix and shows the fold of the potato 1 family of proteinase inhibitors. Compared to known structures of the family, LUTI contains Arg and Trp residues at positions P6' and P8', respectively, instead of two Arg residues, involved in the proteinase binding loop stabilization. A consequence of the ArgTrp substitution at P8' is a slightly more compact conformation of the loop relative to the protein core.     

Study of bradykinin conformation in a dimethylsulfoxide solution by two-dimensional 1H-NMR spectroscopy

The role of charged groups of the nonapeptide bradykinin in stabilization of its spatial structure in dimethyl sulfoxide solution was investigated. The signal assignment in the 1H-NMR spectra was achieved by means of two dimensional correlated spectroscopy (COSY) and nuclear Overhauser enhancement spectroscopy (NOESY). The changes in the NH and C alpha H proton chemical shifts of the Arg1 and Arg9 residues, variations both in temperature coefficients of chemical shifts of NH-resonances and coupling constants, as well as the appearance of additional NOE cross-peaks in NOESY spectra for d alpha N and d beta N 1H-1H distances were revealed by comparing the NMR spectra for two states--with the protonated C-terminal carboxyl group and deprotonated one. The experimental results are in agreement with the assumption that the conformation of the peptide in (CD3)2SO is stabilized by electrostatic interaction between the oppositely charged N- and C-terminal groups. The conformation with deprotonated alpha-carboxyl group is characterized by two beta-turns in the sequences Pro2-Pro-Gly-Phe5 and Ser6-Pro-Phe-Arg9.     

Sequence-specific 1H NMR assignments and secondary structure of porcine motilin

The solution structure of the 22-residue peptide hormone motilin has been studied by circular dichroism and two-dimensional 1H nuclear magnetic resonance spectroscopy. Circular dichroism spectra indicate the presence of alpha-helical secondary structure in aqueous solution, and the secondary structure can be stabilized with hexafluoro-2-propanol. Sequence-specific assignments of the proton NMR spectrum of porcine motilin in 30% hexafluoro-2-propanol have been made by using two-dimensional NMR techniques. All backbone proton resonances (NH and alpha CH) and most of the side-chain resonances have been assigned by using double-quantum-filtered COSY, RELAYED-COSY, and NOESY experiments. Simulations of NOESY cross-peak intensities as a function of mixing time indicate that spin diffusion has a relatively small effect in peptides the size of motilin, thereby allowing the use of long mixing times to confidently make assignments and delineate secondary structure. Sequential alpha CH-NH and NH-NH NOESY connectivities were observed over a significant portion of the length of the peptide. A number of medium-range NOESY cross-peaks indicate that the peptide is folded into alpha-helix from Glu9 to Lys20, which agrees favorably with the 50% helical content determined from CD measurements. The intensities of selected NOESY cross-peaks relative to corresponding diagonal peaks were used to estimate a rotational correlation time of approximately 2.5 ns for the peptide, indicating that the peptide exists as a monomer in solution under the conditions used here.     

A 2D NMR study of the internal flexibility of the antifungal peptide stendomycin

A 2-D ¹H NMR study (NOESY, COSY, HOHAHA and ROESY experiments) of the antifungal peptide stendomycin is presented. The variation of the NOESY cross peak intensities is measured as a function of temperature in order to discriminate between constant and fluctuating interproton distances. It is shown that among 71 NOESY cross peaks, only 12 correspond to well defined interproton distances and their correlation time is determined. The other cross peaks cannot be translated accurately in terms of distances owing to internal molecular motions. ¹H}–¹³C nOe measurements confirm the internal mobility of the molecule. Finally a flexibility map of stendomycin can be established. Offprint requests to: J. P. Simorre     

Experimental and theoretical studies of the conformational perturbations induced by an abasic site

The three-dimensional structure of the natural undecamer duplex d(CGCACACACGC). d(GCGTGTGTGCG) has been determined by the combined use of NMR spectroscopy and restrained molecular dynamics (rMD) and also by molecular mechanics calculations using the JUMNA program without experimental distance constraints. Both procedures have also been used to model the abasic structure d(CGCACOCACGC).d(GCGTGTGTGCG), where 'O' indicates a modified abasic site: 3-hydroxy-2-(hydroxymethyl) tetrahydrofuran. For the natural duplex, 134 interproton distances have been obtained by complete relaxation matrix analysis of the NOESY cross-peaks intensities, using MARDIGRAS software. These distances along with 100 torsion angles for sugar ring and additional data derived from canonical A and B-DNA, have been used for structures refinement by restrained molecular dynamics. Comparison of the natural oligomer with the abasic structure obtained earlier by NMR/rMD (Y. Coppel, N. Berthet, C. Coulombeau, Ce. Coulombeau, J. Garcia and J. Lhomme, Biochemistry 36, 4817-4830, 1997) confirms that the creation of an abasic site, in this sequence context, leads to marked helix kinking. It is also shown that the JUMNA procedure is capable of reproducing the overall structural features of the natural and damaged DNA conformations without the use of experimental constraints.     

A joint 2D NMR and theoretical investigation of Ca2+ binding loops III and IV of calmodulin

A 2D NMR NOESY spectrum of integral CaM in water(148 residues) reveals a series of downfield-shifted crosspeaks stemming from the NH protons of the Ca2(+)-binding loops III and IV. Their attribution, with the help of already assigned proton resonances of isolated tryptic fragments, was complemented by means of energy-minimizations on the Ca2+ complexes of loops III and IV. From these calculations, a set of two alternative, related conformations was obtained for each loop. The first type of conformation provides a coordination pattern for Ca2+ that is similar to that found in loop EF of parvalbumin. The computed interproton distances in both loops are fully compatible with the inferences from the sets of NOESY cross-peaks. Evidence is also provided for interloop interactions.     

NMR structural refinement of an extrahelical adenosine tridecamer d(CGCAGAATTCGCG)2 via a hybrid relaxation matrix procedure

Until very recently interproton distances from NOESY experiments have been derived solely from the two-spin approximation method. Unfortunately, even at short mixing times, there is a significant error in many of these distances. A complete relaxation matrix approach employing a matrix eigenvalue/eigenvector solution to the Bloch equations avoids the approximation of the two-spin method. We have calculated the structure of an extrahelical adenosine tridecamer oligodeoxyribonucleotide duplex, d(CGCAGAATTCGCG)2, by an iterative refinement approach using a hybrid relaxation matrix method combined with restrained molecular dynamics calculations. Distances from the 2D NOESY spectra have been calculated from the relaxation rate matrix which has been evaluated from a hybrid NOESY volume matrix comprising elements from the experiment and those calculated from an initial structure. The hybrid matrix derived distances have then been used in a restrained molecular dynamics procedure to obtain a new structure that better approximates the NOESY spectra. The resulting partially refined structure is then used to calculate an improved theoretical NOESY volume matrix which is once again merged with the experimental matrix until refinement is complete. Although the crystal structure of the tridecamer clearly shows the extrahelical adenosine looped out way from the duplex, the NOESY distance restrained hybrid matrix/molecular dynamics structural refinement establishes that the extrahelical adenosine stacks into the duplex.     

A Joint 2D NMR and Theoretical Investigation of Ca2+ Binding Loops III and IV of Calmodulin

Abstract

A 2D NMR NOESY spectrum of integral CaM in water(148 residues) reveals a series of downfield-shifted crosspeaks stemming from the NH protons of the Ca²⁺ -binding loops III and IV. Their attribution, with the help of already assigned proton resonances of isolated tryptic fragments, was complemented by means of energy-minimizations on the Ca²⁺ complexes of loops III and IV. From these calculations, a set of two alternative, related conformations was obtained for each loop. The first type of conformation provides a coordination pattern for Ca²⁺ that is similar to that found in loop EF of parvalbumin. The computed interproton distances in both loops are fully compatible with the inferences from the sets of NOESY cross-peaks. Evidence is also provided for interloop interactions.     

Solution structure of the phosphocarrier protein HPr from Bacillus subtilis by two-dimensional NMR spectroscopy.

The solution structure of the phosphocarrier protein, HPr, from Bacillus subtilis has been determined by analysis of two-dimensional (2D) NMR spectra acquired for the unphosphorylated form of the protein. Inverse-detected 2D (1H-15N) heteronuclear multiple quantum correlation nuclear Overhauser effect (HMQC NOESY) and homonuclear Hartmann-Hahn (HOHAHA) spectra utilizing 15N assignments (reported here) as well as previously published 1H assignments were used to identify cross-peaks that are not resolved in 2D homonuclear 1H spectra. Distance constraints derived from NOESY cross-peaks, hydrogen-bonding patterns derived from 1H-2H exchange experiments, and dihedral angle constraints derived from analysis of coupling constants were used for structure calculations using the variable target function algorithm, DIANA. The calculated models were refined by dynamical simulated annealing using the program X-PLOR. The resulting family of structures has a mean backbone rmsd of 0.63 A (N, C alpha, C', O atoms), excluding the segments containing residues 45-59 and 84-88. The structure is comprised of a four-stranded antiparallel beta-sheet with two antiparallel alpha-helices on one side of the sheet. The active-site His 15 residue serves as the N-cap of alpha-helix A, with its N delta 1 atom pointed toward the solvent to accept the phosphoryl group during the phosphotransfer reaction with enzyme I. The existence of a hydrogen bond between the side-chain oxygen atom of Tyr 37 and the amide proton of Ala 56 is suggested, which may account for the observed stabilization of the region that includes the beta-turn comprised of residues 37-40. If the beta alpha beta beta alpha beta (alpha) folding topology of HPr is considered with the peptide chain polarity reversed, the protein fold is identical to that described for another group of beta alpha beta beta alpha beta proteins that include acylphosphatase and the RNA-binding domains of the U1 snRNP A and hnRNP C proteins.     

Intermolecular relaxation has little effect on intra-peptide exchange-transferred NOE intensities

Exchange-transferred nuclear Overhauser enhancement (etNOE) provides a useful method for determining the 3-dimensional structure of a ligand bound to a high-molecular-weight complex. Some concern about the accuracy of such structures has arisen because indirect relaxation can occur between the ligand and macromolecule. Such indirect relaxation, or spin diffusion, would lead to errors in interproton distances used as restraints in structure determination. We address this concern by assessing the extent of intermolecular spin diffusion in nineteen peptide-protein complexes of known structure. Transferred NOE intensities were simulated with the program CORONA (Calculated OR Observed NOESY Analysis) using the rate-matrix approach to include contributions from indirect relaxation between protein-ligand and intraligand proton pairs. Intermolecular spin diffusion contributions were determined by comparing intensities calculated with protonated protein to those calculated with fully deuterated protein. The differences were found to be insignificant overall, and to diminish at short mixing times and high mole ratios of peptide to protein. Spin diffusion between the peptide ligand and the protein contributes less to the etNOE intensities and alters fewer cross peaks than the well-studied intramolecular spin diffusion effects. Errors in intraligand interproton distances due to intermolecular relaxation effects were small on average and can be accounted for with the restraint functions commonly used in NMR structure determination methods. In addition, a rate-matrix approach to calculate distances from etNOESY intensities using a volume matrix comprising only intraligand intensities was found to give accurate values. Based on these results, we conclude that structures determined from etNOESY data are no less accurate due to spin diffusion than structures determined from conventional NOE intensities.     

The effect of selective deuteration on magnetization transfer in larger proteins

Summary The effects of selective deuteration on calculated NOESY intensities have been analyzed for the structure of theE. coli trp aporepressor, a 25 kDa protein. It is shown that selectively deuteratedtrp aporepressor proteins display larger calculated NOESY intensities than those for the same interproton distances in the natural abundance protein. The relatively larger magnetization transfer is demonstrated by a comparison of the NOE build-up curves for specific proton pairs, and for the calculated NOE intensities of short-range NOEs to backbone amide protons. This increase in intensity is especially pronounced for the NH¹–NH¹⁺¹ cross peaks in the -helical regions, and particularly for amide protons of two sequential deuterated residues. The effect is shown to be further intensified for longer mixing times. It is also shown that in all cases, each amide proton exhibits stronger NOEs to its own side chain, with an enhanced effect for deuterated derivatives. This theoretical analysis demonstrates that an evaluation of the relative NOE intensities for different selectively deuterated analogs may be an important tool in assigning NMR spectra of large proteins. These results also serve as a guide for the interpretation of NOEs in terms of distances for structure calculations based on data using selectively deuterated proteins.     

Nuclear Overhauser effect and cross-relaxation rate determinations of dihedral and transannular interproton distances in the decapeptide tyrocidine A.

The following interproton distances are reported for the decapeptide tyrocidine A in solution: (a) r(phi) distances between NH(i) and H alpha (i), (b) r(psi) distances between NH (i + 1) and H alpha (i), (c) r(phi psi) distances between NH(i + 1) and NH(i), (d) NH in equilibrium NH transannular distances, (e) H alpha in equilibrium H alpha transannular distances, (f) r x 1 distances between H alpha and H beta protons, (g) NH(i) in equilibrium H beta (i) distances, (h) NH (i + 1) in equilibrium H beta (i) distances, (i) carboxamide-backbone protons and carboxamide-side chain proton distances, (j) side chain proton-side chain proton distances. The procedures for distance calculations were: NOE ratios and calibration distances, sigma ratios and calibration distances, and correlation times and sigma parameters. The cross-relaxation parameters were obtained from the product, say, of NOE 1 leads to 2 and the monoselective relaxation rate of proton 2; the NOEs were measured by NOE difference spectroscopy. The data are consistent with a type I beta-turn/ type II' beta-turn/ approximately antiparallel beta-pleated sheet conformation of tyrocidine A in solution and the NOEs, cross-relaxation parameters, and interproton distances serve as distinguishing criteria for beta-turn and beta-pleated sheet conformations. It should be borne in mind that measurement of only r phi and r psi distances for a decapeptide only defines the ( phi, psi)-space in terms of 4(10) possible conformations; the distances b-j served to reduce the degeneracy in possible (phi, psi)-space to one tyrocidine A conformation. The latter conformation is consistent with that derived from scalar coupling constants, hydrogen bonding studies, and proton-chromophore distance measurement, and closely resembles the conformation of gramicidin S.     

Removal of zero-quantum interference in NOESY spectra of proteins by utilizing the natural inhomogeneity of the radiofrequency field

Summary A single scan method for the suppression of signals arising from zero-quantum coherences (ZQC) is analysed with respect to its application to NMR experiments on proteins. The ZQC are dephased during a spinlock period due to the natural RF inhomogeneity of a commercial probe. A quantitative analysis of a ZQC-compensated NOESY experiment is given. Although the build-up curve for the cross peaks in ZQC-compensated NOESY experiments differ from those in uncompensated experiments, interproton distances in medium-sized proteins can be evaluated with high accuracy. The proposed method is compared with other techniques for ZQC suppression.     

Conformational variations in d(TGACGTCA) and its reverse sequence d(ACTGCAGT): a joint circular dichroism and nuclear magnetic resonance study

Circular dichroism (CD) and nuclear magnetic resonance (NMR) techniques have been used to characterize the structural properties of the two self-complementary DNA octamers d(TGACGTCA) (I) and d(ACTGCAGT) (II). These display as distinctive features reverse sequences and central steps CpG and GpC, respectively. CD experiments lead to B-form DNA spectra characterized by larger magnitude signals in the case of octamer (I). NMR COSY spectra indicate that in the two octamers all the residues are predominantly south, S, (2'-endo) sugar conformation. NMR NOESY spectra show most of the glycosidic angles confined in the range predicted for B-form DNA although important heterogeneity is noticed along the chains, more pronounced in the case of octamer (I). Both the increase of north, N, (3'-endo) sugar conformation and P (pseudorotation phase angle) deviation from its standard B-form DNA value (162 degrees) express local sequence dependent structure distortions, remarkably visible in CpG step of octamer (I) and agreeing with NOESY cross-peaks intensities. Results interpreted according to Calladine's rules indicate higher cross-chain strains in octamer (I) than in octamer (II). All together, we find evidence to support for octamer (I) in solution local structures with A-DNA properties likely dictated by the central CpG step. Such structures could be involved in the DNA recognition by proteins and anticancerous drugs.     

Determination of cross relaxation rates by 3D NOE-NOE spectroscopy

Summary A method for quantitative determination of cross-relaxation rates of macromolecules in solution is developed. The method is based on the analysis of the intensities of cross peaks in 3D NOE-NOE spectra. The linear combination of the intensities of 3D peaks (spin-diffusion peaks, back-transfer peaks) results in an expression directly proportional to the cross-relaxation rate. The proposed approach allows to determine interproton distances in macromolecules more accurately.     

Relaxation matrix analysis of two-dimensional nuclear Overhauser effect spectra

The large number of interproton distances extracted from two-dimensional nuclear Overhauser effect spectra has enabled determination of biomolecular structures in solution. The accuracy of those distances is increased substantially and the number of distances increased significantly by analysis of the experimental peak intensities using a complete relaxation matrix approach. More distances and more accurate distances both lead to a higher resolution structure. A complete relaxation matrix analysis also enables simulation of peak intensities for any postulated structure; comparison of these intensities with experimental intensities can provide a guide for structure refinement as well as a measure of the quality of the structure derived.     

反向DPPC脂囊泡中短杆菌肽S的二维氢核磁共振研究

本文用300MHz核磁共振仪对反向DPPC脂囊泡中短杆菌肽S进行了研究.用二维相关谱(COSY)和二维NOE谱(NOESY)对反向脂囊泡中短杆菌肽S的共振峰进行了识别,短杆菌肽S所有的共振峰都被指定.并在此基础上,通过NOESY及一维自旋回波谱对短杆菌肽S在反向脂囊泡中的构象进行了分析.结果表明,短杆菌肽S在反向脂囊泡中不再为平面环状结构,而是有某种程度的折叠.     

Conformation of NADP+ bound to a type II dihydrofolate reductase

Type II dihydrofolate reductases (DHFRs) encoded by the R67 and R388 plasmids are sequence and structurally different from known chromosomal DHFRs. These plasmid-derived DHFRs are responsible for confering trimethoprim resistance to the host strain. A derivative of R388 DHFR, RBG200, has been cloned and its physical properties have been characterized. This enzyme has been shown to transfer the pro-R hydrogen of NADPH to its substrate, dihydrofolate, making it a member of the A-stereospecific class of dehydrogenases [Brito, R. M. M., Reddick, R., Bennett, G. N., Rudolph, F. B., & Rosevear, P. R. (1990) Biochemistry 29,9825]. Two distinct binary RBG200.NADP+ complexes were detected. Addition of NADP+ to RBG200 DHFR results in formation of an initial binary complex, conformation I, which slowly interconverts to a second more stable binary complex, conformation II. The binding of NADP+ to RBG200 DHFR in the second binary complex was found to be weak, KD = 1.9 +/- 0.4 mM. Transferred NOEs were used to determine the conformation of NADP+ bound to RBG200 DHFR. The initial slope of the NOE buildup curves, measured from the intensity of the cross-peaks as a function of the mixing time in NOESY spectra, allowed interproton distances on enzyme-bound NADP+ to be estimated. The experimentally measured distances were used to define upper and lower bound distance constraints between proton pairs in distance geometry calculations. All NADP+ structures consistent with the experimental distance bounds were found to have a syn conformation about the nicotinamide-ribose (X = 94 +/- 26 degrees) and an anti conformation about the adenine-ribose (X = -92 +/- 32 degrees) glycosidic bonds.(ABSTRACT TRUNCATED AT 250 WORDS)     

Quantification of DNA structure from NMR data: conformation of d-ACATCGATGT

Resonance assignments of nonexchangeable base and sugar protons have been obtained in double-helical d-ACATCGATGT by using two-dimensional correlated spectroscopy (COSY) and nuclear Overhauser enhancement spectroscopy (NOESY). The exchangeable imino protons have been assigned on the basis of their chemical shifts. The characteristic phase-sensitive multiplet patterns of the intrasugar cross-peaks in the omega 1-scaled COSY spectrum have been used to estimate several scalar coupling constants (J). The information on the J values combined with the intranucleotide COSY cross-peak intensities has been used to identify sugar puckers of individual nucleotide units. In most cases, the deoxyribofuranose rings are found to adopt a conformation close to O4'-endo. Spin diffusion has been monitored from the buildup of the normalized volumes of NOE cross-peaks in NOESY spectra as a function of mixing time. A set of 52 intranucleotide and internucleotide proton-proton distances have been estimated by using low mixing time NOESY spectra (tau m = 40 and 80 ms). The estimated intrasugar proton-proton distances rule out possibilities of existence of a fast equilibrium between C2'-endo and C3'-endo conformations. Intranucleotide proton-proton distances combined with the knowledge of sugar puckers have been used to fix the glycosidic bond torsion angle (chi). For this purpose, simulated distance contours depicting the dependence of intranucleotide proton-proton distances on pseudorotational phase angle (P) and glycosidic bond torsion angle (chi) have been used. Further, the proton homonuclear (J, delta) spectrum has been used to monitor the 31P-1H heteronuclear couplings, which are preserved in the omega 2 projection.(ABSTRACT TRUNCATED AT 250 WORDS)     

The structure of the double-stranded RNA pentamer 5'(CACAG) . 5'(CUGUG) determined by nuclear Overhauser enhancement measurements: interproton distance determination and structure refinement on the basis of X-ray coordinates

The structure in solution of the duplex RNA pentamer 5'(CACAG) . 5'(CUGUG), comprising the stem of the T psi C loop of yeast tRNAPhe, has been investigated by means of one- and two-dimensional nuclear Overhauser enhancement measurements. All non-exchangeable base and sugar proton resonances with the exception of the H5'/H5" sugar resonances are assigned in a sequential manner. From the relative intensities of the cross-peaks obtained in the pure-phase absorption two-dimensional nuclear Overhauser enhancement spectra at several mixing times, it is deduced that the RNA pentamer adopts an A-type conformation in solution. Cross-relaxation rates and interproton distances are determined from the time dependence of the nuclear Overhauser effects, principally by one-dimensional measurements. The structure of the RNA pentamer is then refined by restrained least-squares minimization on the basis of both distance and planarity restraints using fibre diffraction data as an initial model. The refined structure of the RNA pentamer is of the A type but exhibits local structural variations in glycosidic bond and backbone torsion angles as well as in propeller twist, base roll and base tilt angles.