Two-dimensional Nuclear Overhauser Enhancement Spectra Research Articles

Structure Revision of Plakotenin Based on Computational Investigation of Transition States and Spectroscopic Properties

We show that the previously [Tetrahedron Lett.1992, 33, 2579] proposed structure of natural plakotenin must be revised. Recently, the total synthesis of plakotenin was achieved via an intramolecular Diels-Alder reaction from a (E,E,Z,E)-tetraene as linear precursor. Using density functional theory, the computation of the four possible transition states for this reaction shows that the previously proposed structure could only have been formed via an energetically high-lying transition state, which is very unlikely. Instead, we suggest that the structure of plakotenin corresponds to the product formed via the lowest transition state. A comparison of experimental and theoretical optical rotation, circular dichroism, and two-dimensional nuclear Overhauser enhancement spectra conclusively proves that the structure of plakotenin is the one that is suggested by the transition state computations. Moreover, the simulation of the nuclear Overhauser enhancement spectra suggests that it is most likely that the misassignment of the (1)H chemical shifts of two methyl groups has led to the wrong structure prediction in the 1992 work. The previously proposed structure of iso-plakotenin remains unaffected by our structure revision, but the structures of homo- and nor-plakotenin must also be revised. The present work shows how the total synthesis of a natural product, together with the theoretical determination of the barrier heights of the reactions involved, can be of great help to assign its structure. It appears that intramolecular Diels-Alder reactions can be modeled accurately by today's first-principles methods of quantum chemistry.

Rapid Optimization of a Peptide Inhibitor of Malaria Parasite Invasion by Comprehensive N-Methyl Scanning

Apical membrane antigen 1 (AMA1) of the malaria parasite Plasmodium falciparum has been implicated in the invasion of host erythrocytes and is an important vaccine candidate. We have previously described a 20-residue peptide, R1, that binds to AMA1 and subsequently blocks parasite invasion. Because this peptide appears to target a site critical for AMA1 function, it represents an important lead compound for anti-malarial drug development. However, the effectiveness of this peptide inhibitor was limited to a subset of parasite isolates, indicating a requirement for broader strain specificity. Furthermore, a barrier to the utility of any peptide as a potential therapeutic is its susceptibility to rapid proteolytic degradation. In this study, we sought to improve the proteolytic stability and AMA1 binding properties of the R1 peptide by systematic methylation of backbone amides (N-methylation). The inclusion of a single N-methyl group in the R1 peptide backbone dramatically increased AMA1 affinity, bioactivity, and proteolytic stability without introducing global structural alterations. In addition, N-methylation of multiple R1 residues further improved these properties. Therefore, we have shown that modifications to a biologically active peptide can dramatically enhance activity. This approach could be applied to many lead peptides or peptide therapeutics to simultaneously optimize a number of parameters.

Separation of cross-relaxation and exchange in two-site spin systems with weak spin-spin couplings

A ratio relating the cross-relaxation rates, rate constants, and cross-peak and diagonal peak integrals is used for the analysis of exchanging two-site spin systems with weak scalar spin-spin couplings to extract the rate constants and cross-relaxation rates contributing to the same cross-peaks in two-dimensional nuclear Overhauser enhancement spectra (2-D NOESY). The method to separate contributions from cross-relaxation and chemical exchange into the cross-peaks in NOESY spectra was applied to investigate the intramolecular dynamics of trisulfide 4H, 8H-naphtho[1,8-ef]1,2,3-trithiocine in solutions.

Interaction of urea with an unfolded protein The DNA-binding domain of the 434-repressor

Experimental techniques are presented for the observation of the solvation of the unfolded form of a globular protein, the N-terminal 63-residue polypeptide from the 434 repressor, in 7 M aqueous urea solution by both water and urea. With the use of 15N-labelled urea it is demonstrated that the cross sections through two-dimensional nuclear Overhauser enhancement (NOE) spectra at the chemical shifts of H 2O and urea both contain direct NOEs with the protein, under conditions where exchange peaks are observed only in the water cross section. A preliminary analysis of the data showed that the residence times of urea molecules in solvation sites near the methyl groups of Val, Leu and Ile are significantly longer than those of water molecules in the same sites.

Refined structure in solution of the sea anemone neurotoxin ShI.

A high resolution structure in solution has been determined for the polypeptide neurotoxin I, from the sea anemone Stichodactyla helianthus, using NMR data, distance geometry calculations, and refinement by back-calculation of two-dimensional nuclear Overhauser enhancement (NOE) spectra. A set of 913 distance constraints derived from NOEs was used, together with a large set of lower distance bounds based on the absence of NOEs in the spectrum. Eight published structures for neurotoxin I were refined independently to give structures which agree better with the experimental data, as reflected in reduced R factors calculated over well resolved cross-peaks of the two-dimensional NOE spectra and a lower total volume of peaks in back-calculated spectra that are absent from experimental spectra. The refined structures are also more precisely defined, with mean pairwise root mean square differences over backbone heavy atoms of 0.62 A for well defined residues and 1.14 A for all residues, compared with previous values of 1.09 and 2.41 A, respectively. The consensus constraint set from the 8 refined structures was also used to generate 12 new structures, with corresponding root mean square differences of 0.76 and 1.26 A. In all 20 structures the loop linking the first and second strands of the beta-sheet is considerably better defined than before. A type I beta-turn encompassing the functionally important residues Asp6, Asp7, and Glu8 is evident in the refined structures.

Chiral ligands: unambiguous assignment of absolute configuration by NMR spectroscopy

Two-dimensional rotating-frame nuclear Overhauser enhancement (ROESY) spectra are used to determine the relative configurations at the chiral carbon atoms in the two diastereoisomers of [o-[graphic omittted]d(Ph2PCHMeCHMePPh2)]+PF6– derived from (+) Ph2PCHMeCHMePPh21 and [o-[graphic omittted]dCl]2 of known configurations, and hence to show that the (+) form of 1 has the R,Rconfiguration.

Two-dimensional 1H NMR studies of 32-base-pair synthetic immobile Holliday junctions: complete assignments of the labile protons and identification of the base-pairing scheme.

The labile protons of two 32-base-pair, four-arm models of immobile Holliday junctions have been studied by two-dimensional 1H nuclear magnetic resonance (NMR) spectroscopy. Overlap of resonances in the imino-imino region of two-dimensional nuclear Overhauser enhancement (NOE) spectra necessitates the use of a multi-pathway approach for obtaining sequence-specific assignments wherein all possible NOE connectivities to the labile protons are utilized, including those from the 2H of adenine, 5CH3 of thymine, and 5H of cytosine. Resonance assignments are obtained for all slowly exchanging imino and cytosine amino protons. Base-pairing up to and including the junction point is found in all four arms of both Holliday junctions. Several cross-arm NOE connectivities are identified and can be used to infer the geometry of the helical stacking domains. The two Holliday junctions studied, which differ only by the exchange of two base pairs at the branch point, appear to have opposite arm stacking geometries. These assignments form an important part of the critical background for detailed NMR analysis of Holliday junction structure and dynamics.

Solution structure of [d(GCGTATACGC)] 2

The solution structure of the alternating pyrimidine-purine DNA duplex [d(GCGTATACGC)] 2 has been determined using two-dimensional nuclear magnetic resonance techniques and distance geometry methods. Backbone distance constraints derived from experimental nuclear Overhauser enhancement and J-coupling torsion angle constraints were required to adequately define the conformation of the inter-residue backbone linkages and to avoid underwinding of the duplex. The distance geometry structures were further refined by back-calculation of the two-dimensional nuclear Overhauser enhancement spectra to correct spin-diffusion distance errors. Fifteen final structures for [d(GCGTATACGC)] 2 were generated from the refined experimental distance bounds. These structures all exhibit fully wound B-form geometry with small penalty values (< 1.5 Å) against the distance bounds and small pair-wise root-mean-square deviation values (typically 0.6 Å to 1.5 Å). The final structures exhibit positive base-pair inclination with respect to the helix axis, a marked alternation in rise and twist, and are shorter and wider than classical fiber B-form DNA. The purines were found to adopt a sugar pucker close to the C-2′- endo conformation while pyrimidine sugars exhibited significantly lower pseudorotation phase angles in the C-1′- exo to C-2′- endo range. The minor groove crossstrand steric clashes at pyrimidine-purine steps that would exist in pure B-DNA are attenuated by an increased rise at these steps (and an increased roll angle at TpA steps). Concomitantly the backbone torsion angles of the pyrimidine moieties have larger γ values, larger ε values, and smaller ζ values than the purines. The structures generated by distance geometry methods were also compared with those obtained from restrained molecular dynamics with empirical force-field potentials. The results indicate that the nuclear magnetic resonance/distance geometry approach alone is capable of elucidating most of the salient structural features of double-stranded helical nucleic acids in solution without resorting to empirical energy potentials and without using any structural assumptions from crystallographic data.

Solution structure of [d(ATGAGCGAATA)] 2: Adjacent G: A mismatches stabilized by cross-strand base-stacking and B II phosphate groups

The solution structure of a rather unusual B-form duplex [d(ATGAGCGAATA)] 2 has been determined using two-dimensional nuclear magnetic resonance (2D-NMR) and distance geometry methods. This sequence forms a stable ten base-pair B-form duplex with 3′ overhangs and two pairs of adjacent G:A mismatches paired via a sheared hydrogen-bonding scheme. All non-exchangeable protons, including the stereo-specific H-5′S/H-5′R of the 3G and 7G residues, were assigned by 2D-NMR. The phosphorus spectrum was assigned using heteronuclear correlation with H-3′ and H-4′ resonances. The complete assignments reveal several unusual nuclear Overhauser enhancements (NOEs) and unusual chemical shifts for the neighboring G:A mismatch pairs and their adjacent nucleotides. Inter-proton distances were derived from time-dependent NOEs and used to generate initial structures, which were further refined by iterative back-calculation of the two-dimensional nuclear Overhauser enhancement spectra; 22 final structures were calculated from the refined distance bounds. All these final structures exhibit fully wound helical structures with small penalty values against the refined distance bounds and small pair-wise root-mean-square deviation values (typically 0.5 Å to 0.9 Å). The two helical strands exchange base stacking at both of the two G:A mismatch sites, resulting in base stacking down each side rather than down each strand of the twisted duplex. Very large twist angles (77 °) were found at the G:A mismatch steps. All the final structures were found to have B II phosphate conformations at the adjacent G:A mismatch sites, consistent with observed downfield 31P chemical shifts and Monte-Carlo conformational search results. Our results support the hypothesis that 31P chemical shifts are related to backbone torsion angles. These B II phosphate conformations in the adjacent G:A mismatch step suggest that hydrogen bonding of the G:A pair G-NH 2 to a nearby phosphate oxygen atom is unlikely. The unusual structure of the duplex may be stabilized by strong interstrand base stacking as well as intrastrand stacking, as indicated by excellent base overlap within the mismatch stacks.

The solution structure of a monomeric insulin

The solution conformation of des-(B26-B30)-insulin (DPI) has been investigated by 1H-NMR spectroscopy. A set of 250 approximate interproton distance restraints, derived from two-dimensional nuclear Overhauser enhancement spectra, were used as the basis of a structure determination using distance geometry (DG) and distance-bound driven dynamics (DDD). Sixteen DG structures were optimized using energy minimization (EM) and submitted to short 5-ps restrained molecular dynamics (RMD) simulations. A further refinement of the DDD structure with the lowest distance errors was done by energy minimization, a prolonged RMD simulation in vacuo and a time-averaged RMD simulation. An average structure was obtained from a trajectory generated during 20-ps RMD. The final structure was compared with the des-(B26-B30)-insulin crystal structure refined by molecular dynamics and the 2-Zn crystal structure of porcine insulin. This comparison shows that the overall structure of des-(B26-B30)-insulin is retained in solution with respect to the crystal structures with a high flexibility at the N-terminal part of the A chain and at the N-terminal and C-terminal parts of the B chain. In the RMD run a high mobility of Gly A1, Asn A21 and of the side chain of Phe B25 is noticed. One of the conformations adopted by des-(B26-B30)-insulin in solution is similar to that of molecule 1 (Chinese nomenclature) in the crystal structure of porcine insulin.

Two-dimensional proton NMR studies of synthetic immobile Holliday junctions

Two 32 base pair, four-arm immobile Holliday junctions have been prepared and studied by two-dimensional 1H nuclear magnetic resonance (NMR) spectroscopy. Two-quantum spectroscopy provides scalar (through bond) correlations for the 1'H, 2'H, and 2''H resonances of the deoxyribose sugar rings and the nonlabile cytosine and thymine base protons. Assignments in the deoxyribose sugars are extended to the 3'H resonances principally from relayed connectivities in total correlation spectra. Severe overlap of resonances in the standard regions of two-dimensional nuclear Overhauser enhancement (NOE) spectra necessitated the use of a unique approach for obtaining sequence-specific assignments of duplex DNA, wherein all possible NOE connectivities in the spectra are analyzed. These studies of 64-residue structures represent a substantial step forward with respect to the size of oligonucleotide for which virtually complete assignments are obtained. The assignments form the critical background for the detailed analysis of Holliday junction structure and dynamics that is required to address key issues in understanding the role of Holliday junctions in genetic recombination and repair.

Solution structure of the kringle 4 domain from human plasminogen by 1H nuclear magnetic resonance spectroscopy and distance geometry

Kringle 4 is an autonomous structural and folding domain within the proenzyme plasminogen. Homologous domains are found throughout the blood clotting and fibrinolytic proteins. In this paper, we present the almost complete assignment of the 1H nuclear magnetic resonance (n.m.r.) spectrum of the kringle 4 domain of human plasminogen. A detailed structural analysis has been completed. The sequential pattern of nuclear Overhauser enhancements indicated little regular secondary structure but rather a series of turns and loops connecting β-strands. A small stretch of antiparallel β-sheet was identified between the residues 61 to 63 and 71 to 73 and the close proximity of other strands was determined from two-dimensional nuclear Overhauser enhancement spectra. Slowly exchanging amide (NH) resonances were found to be associated with residues of the β-sheet and neighbouring strands that support the hydrophobic core of the domain. A total of 526 interproton distance constraints and two hydrogen bonds were specified as input to the distance geometry program DISGEO. Tertiary structures were produced that were consistent with the n.m.r. data. The structures were compared with that of our earlier model based on n.m.r. studies and with that of prothrombin fragment 1 determined crystallographically.

Consequences of magnetization transfer on the determination of solution structures of proteins

The complete two-dimensional nuclear Overhauser enhancement (NOESY) spectra of several proteins were theoretically simulated to analyze the effect that spin diffusion has on the cross-peak intensity-interproton distance relationship as a function of rotational time τ and mixing time t m. The NOESY cross-peak intensities were calculated by solving the generalized Bloch equations from the crystal coordinates of basic pancreatic trypsin inhibitor, cytochrome b562, domain 1 of phage lysozyme, and human deoxy hemoglobin (MW 6500, 11,459, 18,625, and 64,450, respectively). It was found that correlations between cross-peak intensity and internuclear distance, useful for obtaining solution structure of proteins from nuclear magnetic resonance data, can be defined for cross peaks observed at t m = 100 ms, for rotational correlation times τ < 2.3 ns (corresponding to MW < 11,459 at 68°C), and at t m = 50 ms, for τ < 5.9 ns (corresponding to MW < 30,000 at 68°C). Above these limits, the correlations start to break down, and two protons at any given distance may give rise to NOE cross peaks varying greatly in intensity, depending on the local structure of the protein.

1H 2D NMR and distance geometry study of the folding of Ecballium elaterium trypsin inhibitor, a member of the squash inhibitors family.

The solution conformation of synthetic Ecballium elaterium trypsin inhibitor II, a 28-residue peptide with 3 disulfide bridges, has been studied by 1H 2D NMR measurements. Secondary structure elements were determined: a miniantiparallel beta-sheet Met 7-Cys 9 and Gly 25-Cys 27, a beta-hairpin 20-28 with beta-turn 22-25, and two tight turns Asp 12-Cys 15 and Leu 16-Cys 19. A set of interproton distance restraints deduced from two-dimensional nuclear Overhauser enhancement spectra and 13 phi backbone torsion angles restraints were used as the basis of three-dimensional structure computations including disulfide bridges arrangement by using distance geometry calculations. Computations for the 15 possible S-S linkage combinations lead to the proposal of the array 2-19, 9-21, 15-27 as the most probable structure for EETI II.

Conformational studies of d-(AAAAATTTTT)2 using constraints from nuclear overhauser effects and from quantitative analysis of the cross-peak fine structures in two-dimensional 1H nuclear magnetic resonance spectra.

The conformation at the dA-dT junction in d-(AAAAATTTTT)2 was investigated by using a variety of phase-sensitive two-dimensional nuclear magnetic resonance experiments at 500 MHz for detailed studies of the deoxyribose ring puckers. Conformational constraints were collected from two-dimensional nuclear Overhauser enhancement spectra recorded with short mixing times and from quantitative simulations of the cross-peaks in two-dimensional correlated spectra. Overall, the decamer duplex adopts a conformation of the B-DNA type, and for dA4 and dA5 the pseudorotation phase angle P is in the standard range 150-180 degrees. The deoxyribose puckers for the other nucleotides deviate significantly from the standard B-DNA structure. Spectrum simulations assuming either static deviations from standard B-DNA or a simple two-state dynamic equilibrium between the C2'-endo and C3'-endo forms of the deoxyribose were used to analyze the experimental data. It was thus found that the ring pucker for dT6 deviates from the regular C2'-endo form of B-DNA by a static distortion, with the pseudorotation phase angle P in the range 100-130 degrees, and a similar value of P is indicated for dT7. For the peripheral base pairs dynamic distortions of the C2'-endo form of the deoxyribose were found. In agreement with recent papers on related duplexes containing (dA)n tracts, we observed prominent nuclear Overhauser effects between adenine-2H and deoxyribose-1'H, which could be largely due to pronounced propeller twisting as observed in the crystal structures of (dA)n-containing compounds.

Refinement of the solution structure of the DNA dodecamer 5'd(CGCGPATTCGCG)2 containing a stable purine-thymine base pair: combined use of nuclear magnetic resonance and restrained molecular dynamics.

The solution structure of the self-complementary dodecamer 5'd(CGCGPATTCGCG)2, containing a purine-thymine base pair within the hexameric canonical recognition site GAATTC for the restriction endonuclease EcoRI, is investigated by nuclear magnetic resonance spectroscopy and restrained molecular dynamics. Nonexchangeable and exchangeable protons are assigned in a sequential manner. A set of 228 approximate interproton distance restraints are derived from two-dimensional nuclear Overhauser enhancement spectra recorded at short mixing times. These distances are used as the basis for refinement using restrained molecular dynamics in which the interproton distance restraints are incorporated into the total energy function of the system in the form of effective potentials. Eight calculations are carried out, four starting from classical A-DNA and four from classical B-DNA. In all cases convergence to very similar B-type structures is achieved with an average atomic root mean square (rms) difference between the eight converged structures of 0.7 +/- 0.2 A, compared to a value of 6.5 A for that between the two starting structures. It is shown that the introduction of the purine-thymine mismatch does not result in any significant distortion of the structure. The variations in the helical parameters display a clear sequence dependence. The variation in helix twist and propeller twist follows Calladine's rules and can be attributed to the relief of interstrand purine-purine clash at adjacent base pairs. Overall the structure is straight. Closer examination, however, reveals that the central 5 base pair steps describe a smooth bend directed toward the major groove with a radius of curvature of approximately 38 A, which is compensated by two smaller kinks in the direction of the minor groove at base pair steps 3 and 9. These features can be explained in terms of the observed variation in roll and slide.

Location of an alpha-helix in fragment 96-133 from bovine somatotropin by 1H NMR spectroscopy.

By use of two-dimensional NMR techniques most of the proton resonances (greater than 90%) are assigned for the tryptic digest fragment 96-133 of bovine somatotropin in 30% 2,2,2-trifluoroethanol-d3/70% H2O. Qualitative analysis of the nuclear Overhauser enhancement (NOE) data indicates that a region of alpha-helix spans residues 106-128, while the N- and C-terminal regions assume nonregular structures. Amide-exchange rates and comparison of two-dimensional NOE spectra indicate that the most stable piece of helix spans residues 120-125 and that this piece of helix is stable in water at 25 degrees C. Evidence is given to support the fact that intermolecular association of the helical segments stabilizes the helix.

Refinement of the solution structure of the ribonucleotide 5'r(GCAUGC)2: combined use of nuclear magnetic resonance and restrained molecular dynamics

The solution structure of the self-complementary hexamer 5'r(GCAUGC)2 is investigated by means of nuclear magnetic resonance spectroscopy and restrained molecular dynamics. The proton resonances are assigned in a sequential manner, and a set of 110 approximate interproton distance restraints are derived from the two-dimensional nuclear Overhauser enhancement spectra. These distances are used as the basis of a structure refinement by restrained molecular dynamics in which the experimental restraints are incorporated into the total energy function of the system in the form of effective potentials. Eight restrained molecular dynamics simulations are carried out, four starting from a structure with regular A-type geometry and four from one with regular B-type geometry. The atomic root mean square (rms) difference between the initial structures is 3.2 A. In the case of all eight simulations, convergence is achieved both globally and locally to a set of very similar A-type structures with an average atomic rms difference between them of 0.8 +/- 0.2 A. Further, the atomic rms differences between the restrained dynamics structures obtained by starting out from the same initial structures but with different random number seeds for the assignment of the initial velocities are the same as those between the restrained dynamics structures starting out from the two different initial structures. These results suggest that the restrained dynamics structures represent good approximations of the solution structure. The converged structures exhibit clear sequence-dependent variation in some of the helical parameters, in particular helix twist, roll, slide, and propellor twist.(ABSTRACT TRUNCATED AT 250 WORDS)

Sequential NMR assignments of labile protons in DNA using two-dimensional nuclear-Overhauser-enhancement spectroscopy with three jump-and-return pulse sequences.

Two-dimensional nuclear Overhauser enhancement (NOESY) spectra of labile protons were recorded in H2O solutions of a protein and of a DNA duplex, using a modification of the standard NOESY experiment with all three 90 degree pulses replaced by jump-and-return sequences. For the protein as well as the DNA fragment the strategically important spectral regions could be recorded with good sensitivity and free of artifacts. Using this procedure, sequence-specific assignments were obtained for the imino protons, C2H of adenine, and C4NH2 of cytosine in a 23-base-pair DNA duplex which includes the 17-base-pair OR3 repressor binding site of bacteriophage lambda. Based on comparison with previously published results on the isolated OR3 binding site, these data were used for a study of chain termination effects on the chemical shifts of imino proton resonances of DNA duplexes.

The polypeptide fold of the globular domain of histone H5 in solution. A study using nuclear magnetic resonance, distance geometry and restrained molecular dynamics.

The polypeptide fold of the 79-residue globular domain of chicken histone H5 (GH5) in solution has been determined by the combined use of distance geometry and restrained molecular dynamics calculations. The structure determination is based on 307 approximate interproton distance restraints derived from n.m.r. measurements. The structure is composed of a core made up of residues 3-18, 23-34, 37-60 and 71-79, and two loops comprising residues 19-22 and 61-70. The structure of the core is well defined with an average backbone atomic r.m.s. difference of 2.3 +/- 0.3 A between the final eight converged restrained dynamics structures and the mean structure obtained by averaging their coordinates best fitted to the core residues. The two loops are also well defined locally but their orientation with respect to the core could not be determined as no long range ([i-j[ greater than 5) proton-proton contacts could be observed between the loop and core residues in the two-dimensional nuclear Overhauser enhancement spectra. The structure of the core is dominated by three helices and has a similar fold to the C-terminal DNA binding domain of the cAMP receptor protein.