The three-dimensional NMR solution structure of alpha-cobratoxin at pH 7.5 and conformational differences with the NMR solution structure at pH 3.2

The 3D solution structure of alpha-cobratoxin, a neurotoxin purified from the Naja naja siamensis snake venom, has been determined by Nuclear Magnetic Resonance spectroscopy, in conjunction with distance geometry and restrained molecular dynamics, at pH 7.5. A total of 490 distance restraints were obtained from NOE intensities and 25 phi dihedral angle restraints deduced from J-coupling data. The generated structures are well defined with root mean square deviations from a geometrical mean structure of 0.107 +/- 0.036 nm for the backbone atoms and 0.128 +/- 0.073 nm for the side-chain atoms (considering residues 1 to 66 minus 26 to 35). A comparison between the generated structures at pH 7.5 and the mean NMR solution structure at pH 3.2 revealed that the 3D structure of alpha-cobratoxin is more compact at neutral pH. This major difference is mainly due to the pH-dependent conformational variations of three residues His18, Thr44 and Thr59.     

An NMR study on the DNA-binding SPKK motif and a model for its interaction with DNA

The solution structure of one and two repeats of the 'SPKK' DNA-binding motif is reported on the basis of NMR measurements. In dimethylsulphoxide (DMSO) the major population (approximately 90%) of peptides, SPRKSPRK(S2) and GSPKKSPRK(S2b), adopts a conformation, which has two trans prolines. The two 'SP(R/K)K' units in these peptides are equivalent and each adopts a turn structure exchanging with an extended structure. This is suggested by an NOE connectivity of the beta-turn type, between the backbone amide protons of residues (i+2) and (i+3) and NOE connectivities of the Asx(sigma)-turn type, between protons of the ith Ser and the backbone amide proton on residue (i+2). This suggests that each SP(R/K)K unit has a structural intermediate between (or a combination of) a beta-turn and an Asx(sigma)-turn. In 90-10% DMSO/H2O at 4 degrees C the two units of S2 are connected more tightly by folding into a short 3(10) helix, broken at the second proline. For another peptide, Thr-Pro-Arg-Lys(T1), the major population (75%) in 100% DMSO comprises a beta-turn in rapid exchange with an extended structure. We did not observe an NOE connectivity of the Asx(sigma) type with the T1 peptide. A possible structure of the SPKK motif in the complex with DNA is discussed.     

Comparative conformational studies on cyclic hexapeptides corresponding to message sequence His-Phe-Arg-Trp of alpha-melanotropin by NMR

Solution conformation of cyclo(Gly1-His2-Phe3-Arg4-Trp5-Gly6) and its D-Phe analog corresponding to the message sequence [Gly-alpha-MSH5-10] of alpha-MSH has been studied by 1D and 2D proton magnetic resonance spectroscopy in dimethyl sulfoxide (DMSO)-d6 solution and in a DMSO-d6/H2O cryoprotective mixture. The NMR data for both the analogs in solution at 300 K cannot be interpreted based on a single ordered conformation, as evidenced by the broadening of only -NH resonances as well as the temperature coefficients of the amide protons. An analysis of the nuclear Overhauser effect (NOE) cross-peaks in conjunction with temperature coefficient data indicates an equilibrium of multiple conformers with a substantial population of particular conformational states at least in the D-analog. The molecular dynamics simulations without and with NOE constraints also reveal numerous low-energy conformers with two gamma-turns, a gamma-turn and a beta-turn, two beta-turns, etc. for both the analogs. The observed NMR spectra can be rationalized by a dynamic equilibrium of conformers characterized by a gamma-bend at Gly6, two gamma-bends at Phe3 and Gly6 and a conformer with a single beta-turn and a gamma-bend for the L-Phe analog. On the other hand, a conformation with two fused beta-turns around the two tetrads His2-D-Phe3-Arg4-Trp5 and Trp5-Gly6-Gly1-His2 dominates the equilibrium mixture for the D-Phe analog. For the D-Phe analog, the experimentally observed average conformation is corroborated by molecular dynamics simulations as well as by studies in cryoprotective solvent.     

Conformational interconversions in peptide beta-turns: analysis of turns in proteins and computational estimates of barriers

The two most important beta-turn features in peptides and proteins are the type I and type II turns, which differ mainly in the orientation of the central peptide unit. Facile conformational interconversion is possible, in principle, by a flip of the central peptide unit. Homologous crystal structures afford an opportunity to structurally characterize both possible conformational states, thus allowing identification of sites that are potentially stereochemically mobile. A representative data set of 250 high-resolution (相似文献   

NMR solution structure of the receptor binding domain of Pseudomonas aeruginosa pilin strain P1. Identification of a beta-turn

The solution structure of the peptide antigen from the receptor binding domain of Pseudomonas aeruginosa strain P1 has been determined using two-dimensional 1H NMR techniques. Ensembles of solution conformations for the trans form of this 23-residue disulfide bridged peptide have been generated using a simulated annealing procedure in conjunction with distance and torsion angle restraints derived from NMR data. Comparison of the NMR-derived solution structures of the P1 peptide with those previously determined for the 17-residue PAK, PAO and KB7 strain peptides [McInnes, C., et al. (1993) Biochemistry 32, 13432-13440; Campbell, A.P., et al. (1995) Biochemistry 34, 16255-16268] reveals the common structural motif of a beta-turn, which may be the necessary structural requirement for recognition of a common cell surface receptor and a common cross-reactive antibody to which all four strains bind. The importance of this conserved beta-turn in the PAK, PAO, KB7 and P1 peptides is discussed with regard to the design of a synthetic peptide vaccine effective against multiple strains of Pseudomonas aeruginosa infections.     

Structural redesign and stabilization of the overlapping tandem beta-turns of RNA polymerase II

Peptides representing single repeat units of the carboxy-terminal domain (CTD) of RNA polymerase II (Tyr-Ser-Pro-Thr-Ser-Pro-Ser-Tyr-NH2, 1) contain overlapping Ser-Pro-Xaa-Xaa beta-turn forming sites which permit their overall structure to closely resemble members of the quinoxaline class of antitumor DNA bisintercalators. We have modified this native sequence at the i+2 positions of each beta-turn unit by substituting Gly or D-Ala in an attempt to preorganize this structure in aqueous solution. CD and NMR spectroscopic investigations confirmed the presence of type II beta-turns within each of the substituted peptides in contrast to the native sequence which contains a relatively low population of turn structure. In addition, an examination of singly substituted peptides suggests that an increase in the population of beta-turn structure within the amino-terminal Ser-Pro-Xaa-Xaa site also increased the formation of beta-turn structure in the carboxy-terminal (unmodified) Ser-Pro-Xaa-Xaa site; in comparison, substitution in the carboxy-terminal site did not influence structure in the remaining portion of the peptide. Overall, these results suggest that the structures formed could provide unique, preorganized linkers for the construction of novel DNA-interactive bisintercalators.     

Assigning the NMR spectra of aromatic amino acids in proteins: analysis of two Ets pointed domains

The measurement of interproton nuclear Overhauser enhancements (NOEs) and dihedral angle restraints of aromatic amino acids is a critical step towards determining the structure of a protein. The complete assignment of the resonances from aromatic rings and the subsequent resolution and identification of their associated NOEs, however, can be a difficult task. Shown here is a strategy for assigning the 1H, 13C, and 15N signals from the aromatic side chains of histidine, tryptophan, tyrosine, and phenylalanine using a suite of homo- and hetero-nuclear scalar and NOE correlation experiments, as well as selective deuterium isotope labelling. In addition, a comparison of NOE information obtained from homonuclear NOE spectroscopy (NOESY) and 13C-edited NOESY-heteronuclear single quantum correlation experiments indicates that high-resolution homonuclear two-dimensional NOESY spectra of selectively deuterated proteins are invaluable for obtaining distance restraints to the aromatic residues.     

Investigation of the solution structure of chymotrypsin inhibitor 2 using molecular dynamics: comparison to x-ray crystallographic and NMR data

The native solution structure and dynamics of chymotrypsin inhibitor 2 (CI2) have been studied using a long (5.3 ns) molecular dynamics (MD) simulation without any imposed restraints. The majority of the experimentally observed spin-spin coupling constants, short- and long-range nuclear Overhauser effect (NOE) cross peaks and the amide hydrogen exchange behavior were reproduced by the MD simulation. This good correspondence suggests that the major structural features of the protein during the simulation are representative of the true protein structure in solution. Two water molecules formed hydrogen bond bridges between beta2 and beta3, in agreement with X-ray crystallographic data and a recent reassessment of the solution structure using time-averaged NMR restraints during MD refinement. The active-site loop of the protein displayed the greatest structural changes and the highest mobility. When this loop region was excluded, the average Calpha r.m.s. deviation of the simulated solution structures from the crystal structure was approximately 1.5 Angstrom from 0.5 to 5.3 ns. There is structural heterogeneity in particular regions of the NMR-derived solution structures, which could be a result of imprecision or true internal motion. A study of the distribution of mobility through the protein allows us to distinguish between these two alternatives. In particular, deviations in the active-site loop appear to be a result of heightened mobility, which is also supported by good correspondence between calculated and experimental S2 N-H order parameters. On the other hand, other ill-defined regions of the NMR-derived structures are well defined in the simulation and are probably the result of a lack of structural restraints (i.e. NOEs), as opposed to reflecting the true mobility.     

Solution structure of Lqh-8/6, a toxin-like peptide from a scorpion venom--structural heterogeneity induced by proline cis/trans isomerization

Lqh-8/6 is a minor fraction isolated from the venom of the scorpion Leiurus quinquestriatus hebraeus. Here we describe the purification, amino acid sequencing and solution structure determination by NMR and molecular modeling of this peptide. Lqh-8/6 is a small polypeptide (38 residues) which contains 8 half-cystines and is highly similar to another venom component, chlorotoxin. Standard homonuclear methods were used to sequentially assign the proton NMR spectra and to collect spatial restraints for structure determination. Two populations, identified early in the assignment step, are in slow interconversion on the NMR timescale. The two conformers were shown to originate from a cis/trans peptidyl-prolyl isomerization. Using a distance geometry program and simulated annealing protocol under the NMR restraints we obtained 10 final structures for the major conformation (trans isomer). None of the structures showed NOE violations larger than 0.05 nm, and the rmsd value relative to the mean structure (considering the main chain atoms in well-defined secondary structure) is 0.07 nm. The three-dimensional structure contains a short alpha-helix strapped on a small antiparallel beta-strand and an N-terminal extended fragment. The sequence/structure and structure/function relationships of the new scorpion toxin-like peptide are discussed in the context of the present structure determination. This toxin shows a stable, highly populated cis conformer of a peptidyl-prolyl peptide bond.     

"Ensemble" iterative relaxation matrix approach: a new NMR refinement protocol applied to the solution structure of crambin

The structure in solution of crambin, a small protein of 46 residues, has been determined from 2D NMR data using an iterative relaxation matrix approach (IRMA) together with distance geometry, distance bound driven dynamics, molecular dynamics, and energy minimization. A new protocol based on an "ensemble" approach is proposed and compared to the more standard initial rate analysis approach and a "single structure" relaxation matrix approach. The effects of fast local motions are included and R-factor calculations are performed on NOE build-ups to describe the quality of agreement between theory and experiment. A new method for stereospecific assignment of prochiral groups, based on a comparison of theoretical and experimental NOE intensities, has been applied. The solution structure of crambin could be determined with a precision (rmsd from the average structure) of 0.7 A on backbone atoms and 1.1 A on all heavy atoms and is largely similar to the crystal structure with a small difference observed in the position of the side chain of Tyr-29 which is determined in solution by both J-coupling and NOE data. Regions of higher structural variability (suggesting higher mobility) are found in the solution structure, in particular for the loop between the two helices (Gly-20 to Pro-22).     

Prediction of beta-turns

A residue-coupled model is proposed to predict the beta-turns in proteins. The rates of correct prediction for the 455 beta-turn tetrapeptides and 3807 non-beta-turn tetrapeptides in the training database are 94.7 and 81.3%, respectively. The rates of correct prediction for the 110 beta-turn tetrapeptides and 30,229 non-beta-turn tetrapeptides in the testing database are 80.0 and 80.2%, respectively. Compared with the rates of correct prediction based on the residue-independent model reported previously, the quality of prediction is significantly improved by the new model, implying that the residue-coupled effect along a polypeptide chain is important for the formation of reversal turns, such as beta-turns, during the process of protein folding.     

Insights into the dynamic nature of DNA duplex structure via analysis of nuclear Overhauser effect intensities

Sequence-dependent structures of DNA duplexes in solution can be reliably determined using NMR data if care is taken to determine restraint bounds accurately. This entails use of complete relaxation matrix methods to analyze nuclear Overhauser effect (NOE) spectroscopic cross-peak intensities, yielding accurate distance restraints. NMR studies of various DNA duplexes have suggested that there may be some limited internal motions. First, it is typically not possible to reconcile all vicinal proton coupling constants in deoxyribose rings with a single conformer. In addition, with the increased accuracy of interproton distance measurements afforded by the complete relaxation matrix algorithm MARDIGRAS, we find inconsistencies in certain distances which can most readily be ascribed to limited conformational flexibility, since conformational averaging is nonlinear. As base-sugar interproton distances depend on both sugar pucker and glycosidic torsion angle chi, motion involving these structural variables should be reflected by experimental data. Possible motional models have been considered to account for all of the data for three DNA duplexes. Analysis of intraresidue base-sugar interproton NOE bounds patterns suggests a motional model with individual sugars in equilibrium between S (2'-endo) and N (3'-endo) conformations, with S being preferred. As sugar repuckering is correlated with changes in glycosidic torsion angle chi, different sugar conformers imply different values for chi, but this is insufficient to account for all data. A two-state jump between anti and syn glycosidic conformers was considered, but it was incapable of accounting for all data. However, a model with restricted diffusion (rocking) about the glycosidic bond in addition to sugar repuckering was capable of accommodating all experimental data. This motional model is in qualitative agreement with experimental 13C relaxation-derived order parameter values in a DNA duplex.     

Localized solution structure refinement of an F45W variant of ubiquitin using stochastic boundary molecular dynamics and NMR distance restraints

The local structure within an 8-A radius around residue 45 of a recombinant F45W variant of human ubiquitin has been determined using 67 interproton distance restraints measured by two-dimensional proton NMR. Proton chemical shift evidence indicates that structural perturbations due to the F45W mutation are minimal and limited to the immediate vicinity of the site of mutation. Simulated annealing implemented with stochastic boundary molecular dynamics was applied to refine the structure of Trp 45 and 10 neighboring residues. The stochastic boundary method allowed the entire protein to be reassembled from the refined coordinates and the outlying unrefined coordinates with little distortion at the boundary. Refinement began with four low-energy indole ring orientations of F45W-substituted wild-type (WT) ubiquitin crystal coordinates. Distance restraints were derived from mostly long-range NOE cross peaks with 51 restraints involving the Trp 45 indole ring. Tandem refinements of 64 structures were done using either (1) upper and lower bounds derived from qualitative inspection of NOE crosspeak intensities or (2) quantitative analysis of cross-peak heights using the program MARDIGRAS. Though similar to those based on qualitative restraint, structures obtained using quantitative NOE analysis were superior in terms of precision and accuracy as measured by back-calculated sixth-root R factors. The six-membered portion of the indole ring is nearly coincident with the phenyl ring of the WT and the indole NH is exposed to solvent. Accommodation of the larger ring is accompanied by small perturbations in the backbone and a 120 degrees rotation of the chi 2 dihedral angle of Leu 50.     

NMR assignment and conformational analysis of the antigenic capsular polysaccharide from Streptococcus pneumoniae type 9N in aqueous solution

Complete 1H and 13C NMR assignments, determined by one- and two-dimensional homo- and hetero-nuclear experiments, are reported for the antigenic capsular polysaccharide (CPS) from Streptococcus pneumoniae serotype 9N (S9 in American nomenclature). Distance constraints derived from 1D NOE difference experiments were combined with energy minimisation (simulated annealing) and molecular dynamics (MD) calculations to determine the most favoured conformation of S9 in aqueous solution at 70 degrees C. NOE values simulated for several static conformational models using the NOEMOL program did not correlate well with experimental values, whereas time averaged interproton distances calculated from 500 ps of restrained MD (using the Tropp formalism to account for rapid internal mobility) were in close agreement with experimentally derived distance estimates.     

NMR and molecular dynamics studies of tachykinins: conformation of the C-terminal pentapeptide of substance P(SP7-11)

Substance P belongs to the tachykinin family of neuropeptides which exhibit diverse pharmacological activity. The conformation of Phe1-Phe2-Gly3-Leu4-Met5-NH2 the C-terminal pentapeptide of substance P (SP7-11) has been studied by NMR and molecular dynamics (MD) methods. NMR studies were carried out both in DMSO-d6 and 95% H2O. Based on the observed chemical shifts, 3JNH alpha coupling constants, temperature coefficients of chemical shifts of NH resonances and the pattern of inter- and intraresidue NOE's, a predominantly extended backbone conformation has been deduced for the peptide in both DMSO and H2O. MD calculations carried out in vacuo indicate that the global minimum energy conformation of the molecule is folded with an intramolecular hydrogen bond between the protonated N-terminal and the C-terminal CONH2 group. The simulation shows that beta-turns are energetically unfavourable, while alpha-helices are seen to be unstable for the peptide. gamma-Bends at either Gly3 or Leu4 are the most preferred ones. Simulations carried out in DMSO as well as in water show a preference for a nearly extended conformation.     

Conformational analyses of somatostatin-related cyclic hexapeptides containing peptoid residues

We report the conformational analysis by 1H NMR in DMSO and computer simulations involving distance geometry and molecular dynamics simulations of a series of peptoid analogues of the cyclic hexapeptide c[Phe11-Pro6-Phe7-d-Trp8-Lys9-Thr10] (1). The proline residue in compound 1 is replaced with the peptoid residues N-benzylglycine (Nphe) (compound 2), N-(S)-alpha-methylbenzylglycine [(S)-beta-MeNphe] (compound 3), and N-(R)-alpha-methylbenzylglycine [(R)-beta-MeNphe] (compound 4). The peptoid analogues 2 and 4 exhibit potent binding activities to the hsst2 receptor, while the binding affinities to the hsst5 and to the hsst3 receptors are reduced compared to that of the parent compound 1. Compound 3 shows reduced binding activities to the hsst2, hsst3, and hsst5 receptors compared to compound 1. The results of in vivo assays indicate that these compounds inhibit the growth hormone release but do not affect the insulin release. These peptoid-containing analogues show two sets of NMR signals corresponding to cis and trans conformations of the peptide bond between Phe11 and Nxaa6. We demonstrate that the backbone conformation and the orientation of the relevant side chains of compound 1 are maintained in the cis isomers of the peptoid analogues which adopt a type VI beta-turn centered around residues 11 and 6 and a type II' beta-turn with d-Trp in the i+1 position. The enhanced selectivity of the peptoid-containing analogues compared to compound 1 and the results of the conformational analysis suggest that the presence of a conformationally constrained hydrophobic group in position 6 in complementary topology to the Phe11 side chain enhances selective binding to the hsst2 receptor.     

Design of peptides using alpha,beta-dehydro-residues: synthesis, crystal structure and molecular conformation of Boc-L-Val-delta Phe-delta Phe-L-Val-OCH3

The peptide Boc-L-Val-delta Phe-delta Phe-L-Val-OCH3 was synthesized by the azlactone method in solution phase, and its crystal and molecular structures were determined by x-ray diffraction method. Single crystals were grown by slow evaporation from a methanol/water solution at 6 degrees C. The crystals belong to an orthorhombic space group P212121 with a = 10.478 (6) A, b = 13.953 (I), c = 24.347 (2) and Z = 4. The structure was determined by direct methods and refined by least squares procedure to an R value of 0.052. The structure consists of a peptide and a water molecule. The peptide adopts two overlapping beta-turn conformations of Types II and I' with torsion angles: phi 1 = -54.8 (6) psi 1 = 130.5 (4), phi 2 = 65.8 (5), psi 2 = 12.8 (6), phi 3 = 79.4 (5), psi 3 = 3.9 (7) degrees. The conformation is stabilized by intramolecular hydrogen bonds involving Boc CO and NH of delta Phe3 and CO of Val1 and NH of Val4. The molecules are tightly packed in the unit cell. The crystal structure is stabilized by hydrogen bonds involving NH of delta Phe2 and CO of a symmetry related (x-1/2, 1/2-y, -z) delta Phe2. The solvent-water molecule forms two hydrogen bonds with peptide molecule involving NH of Val1 as an acceptor and another with CO of a symmetry related (1-x, y-1/2, 1/2 -z) delta Phe3 as a donor. These studies indicate that a tetrapeptide with two consecutive delta Phe residues sequenced with valines on both ends adopts two overlapping beta-turns of Types II and I'.     

Refined structure of villin 14T and a detailed comparison with other actin-severing domains

Villin 14T is the amino terminal actin monomer binding domain from the actin-severing and bundling protein villin. Its structure has been determined in solution using heteronuclear multidimensional nuclear magnetic resonance (NMR) spectroscopy (Markus MA, Nakayama T, Matsudaira P, Wagner G. 1994. Solution structure of villin 14T, a domain conserved among actin-severing proteins. Protein Science 3:70-81). An additional nuclear Overhauser effect (NOE) spectroscopy data set, acquired using improved gradient techniques, and further detailed analysis of existing data sets, produced an additional 601 NOE restraints for structure calculations. The overall fold does not change significantly with the additional NOE restraints but the definition of the structure is improved, as judged by smaller deviations among an ensemble of calculated structures that adequately satisfy the NMR restraints. Some of the side chains, especially those in the hydrophobic core of the domain, are much more defined. This improvement in the detail of the solution structure of villin 14T makes it interesting to compare the structure with the crystal structure of gelsolin segment 1, which shares 58% sequence identity with villin 14T, in an effort to gain insight into villin 14T's weaker affinity for actin monomers. Villin 14T has smaller side chains at several positions that make hydrophobic contacts with actin in the context of gelsolin segment 1. The structure is also compared with the structure of the related actin-severing domain, severin domain 2.     

Structure of hen lysozyme in solution

The structure of the 129-residue protein hen lysozyme has been determined in solution by two-dimensional 1H nuclear magnetic resonance methods. 1158 NOE distance restraints, and 68 phi and 24 chi 1 dihedral angle restraints were employed in conjunction with distance geometry and simulated annealing procedures. The overall C alpha root-mean-square deviation from the average for 16 calculated structures is 1.8(+/- 0.2) A, but excluding 14 residues in exposed disordered regions, this value reduces to 1.3(+/- 0.2) A. Regions of secondary structure, and the four alpha-helices in particular, are well defined (C alpha root-mean-square deviation 0.8(+/- 0.3) A for helices). The main-chain fold is closely similar to structures of the protein in the crystalline state. Furthermore, many of the internal side-chains are found in well-defined conformational states in the solution structures, and these correspond well with the conformational states found in the crystal. The general high level of definition of mainchain and many internal side-chains in the solution structures is reinforced by the results of an analysis of coupling constants and ring current shifts. Many side-chains on the surface, however, are highly disordered amongst the set of solution structures. In certain cases this disorder has been shown to be dynamic in origin by the examination of 3J alpha beta coupling constants.     

A calculation strategy for the structure determination of symmetric dimers by 1H NMR

The structure determination of symmetric dimers by NMR is impeded by the ambiguity of inter- and intramonomer NOE crosspeaks. In this paper, a calculation strategy is presented that allows the calculation of dimer structures without resolving the ambiguity by additional experiments (like asymmetric labeling). The strategy employs a molecular dynamics-based simulated annealing approach to minimize a target function. The experimental part of the target function contains distance restraints that correctly describe the ambiguity of the NOE peaks, and a novel term that restrains the symmetry of the dimer without requiring the knowledge of the symmetry axis. The use of the method is illustrated by three examples, using experimentally obtained data and model data derived from a known structure. For the purpose of testing the method, it is assumed that every NOE crosspeak is ambiguous in all three cases. It is shown that the method is useful both in situations where the structure of a homologous protein is known and in ab initio structure determination. The method can be extended to higher order symmetric multimers.