Distribution of chemical shifts in 1H nuclear magnetic resonance spectra of proteins

Neural Networks in Chemosystematic Studies of Asteraceae: A Classification Based on a Dichotomic Approach

This paper describes the application of artificial neural nets as an alternative and efficient method for the classification of botanical taxa based on chemical data (chemosystematics). A total of 28,000 botanical occurrences of chemical compounds isolated from the Asteraceae family were chosen from the literature, and grouped by chemical class for each species. Four tests were carried out to differentiate and classify different botanical taxa. The qualifying capacity of the artificial neural nets was dichotomically tested at different hierarchical levels of the family, such as subfamilies and groups of Heliantheae subtribes. Furthermore, two specific subtribes of the Heliantheae and two genera of one of these subtribes were also tested. In general, the artificial neural net gave rise to good results, with multiple-correlation values R>0.90. Hence, it was possible to differentiate the dichotomic character of the botanical taxa studied.

Prediction of occurrences of diverse chemical classes in the Asteraceae through artificial neural networks

The training and the application of a neural network system for the prediction of occurrences of secondary metabolites belonging to diverse chemical classes in the Asteraceae is described. From a database containing about 604 genera and 28,000 occurrences of secondary metabolites in the plant family, information was collected encompassing nine chemical classes and their respective occurrences for training of a multi-layer net using the back-propagation algorithm. The net supplied as output the presence or absence of the chemical classes as well as the number of compounds isolated from each taxon. The results provided by the net from the presence or absence of a chemical class showed a 89% hit rate; by excluding triterpenes from the analysis, only 5% of the genera studied exhibited errors greater than 10%.

Conformational study of [Met5]enkephalin-Arg-Phe in the presence of phosphatidylserine vesicles

The interaction of [Met5]enkephalin-Arg.Phe with phosphatidylserine (PtdSer) was studied by circular dichroism (CD), two-dimensional nuclear magnetic resonance spectroscopy, hybrid distance geometry simulated annealing (DG-SA) and molecular dynamics (MD) calculations. The very low solubility of [Met5]enkephalin-Arg-Phe and the instability of the solution containing PtdSer vesicles at low pH values did not allow us to observe the amide proton resonances in the usual two-dimensional NMR work. NOESY cross-peaks of protons of side chains from two-dimensional NMR were converted into distances which were used as restraints for modelling with DG-SA and MD. Our results indicate that, in aqueous solutions at pH 7.68 [Met5]enkephalin-Arg-Phe exists in the absence of PtdSer as a random distribution of conformers, whereas in the presence of PtdSer it adopts conformations containing a common orientation of the bonds of C alpha 2, C alpha 3, C alpha 4, and C alpha 5, although different orientations of the peptide planes are consistent with the results. Two of the reported conformers from MD simulations are characterized by the presence of a 2<--4 gamma and inverse gamma turns centered on Gly3. A gradual decline of order was observed when moving from the central moiety of the peptide to both the N-terminus and C-terminus. Finally, the DG-SA and MD calculations resulted in a structure such that the orientation of the Phe4 and Met5 side chains favours hydrophobic interactions with the apolar portion of the PtdSer vesicle to form a hydrophobic cluster. These data support the hypothesis of a role of lipids to modify the conformation of [Met5]enkephalin-Arg-Phe to permit the interactions with the receptor site.

Solution conformation of the Pseudomonas syringae pv. syringae phytotoxic lipodepsipeptide syringopeptin 25-A. Two-dimensional NMR, distance geometry and molecular dynamics

Syringopeptin 25-A is a phytotoxic amphiphilic lipodepsipeptide containing 25 amino acid residues, produced by some isolates of the plant pathogenic bacterium Pseudomonas syringae pv. syringae. Previous papers have reported its covalent structure and some of its biological properties. Attention has now been directed to define its conformation in solution, a structural feature regarded as important for understanding its possible role in the bacterial colonization of host plants, and its toxic action on the plant cell. Here we report the stereochemistry of its amino acid components, the complete interpretation of the two-dimensional NMR spectra and NOE data, and finally the structure obtained by computer simulations applying distance geometry and molecular dynamics procedures. The conformation of syringopeptin 25-A in aqueous solution includes three different structural regions interrupted by rigid 2,3-dehydro-2-aminobutyric acid residues: a loop from residue 2 to 6, a helicoidal zone from 8 to 15, and the lactone ring from 18 to 25. The three-dimensional structure of the lactone moiety is very similar to that of two previously studied bioactive lipodepsinonapeptides. Preliminary circular dichroism evidence of conformational variations in solution of trifluoroethanol, which stimulates a membrane-like environment, are also reported.

The program XEASY for computer-supported NMR spectral analysis of biological macromolecules

A new program package, XEASY, was written for interactive computer support of the analysis of NMR spectra for three-dimensional structure determination of biological macromolecules. XEASY was developed for work with 2D, 3D and 4D NMR data sets. It includes all the functions performed by the precursor program EASY, which was designed for the analysis of 2D NMR spectra, i.e., peak picking and support of sequence-specific resonance assignments, cross-peak assignments, cross-peak integration and rate constant determination for dynamic processes. Since the program utilizes the X-window system and the Motif widget set, it is portable on a wide range of UNIX workstations. The design objective was to provide maximal computer support for the analysis of spectra, while providing the user with complete control over the final resonance assignments. Technically important features of XEASY are the use and flexible visual display of 'strips', i.e., two-dimensional spectral regions that contain the relevant parts of 3D or 4D NMR spectra, automated sorting routines to narrow down the selection of strips that need to be interactively considered in a particular assignment step, a protocol of resonance assignments that can be used for reliable bookkeeping, independent of the assignment strategy used, and capabilities for proper treatment of spectral folding and efficient transfer of resonance assignments between spectra of different types and different dimensionality, including projected, reduced-dimensionality triple-resonance experiments.

A two-dimensional NMR study of bovine Cu, Co superoxide dismutase. Further assignments in the region surrounding the active site

Bovine Cu, Co superoxide dismutase has been investigated by two-dimensional NMR with regard to the resonances arising from protons that surround the copper site. These protons have magnetic properties that are intermediate between those belonging to residues coordinated to the paramagnetic metal center and the diamagnetic residues. Nuclear Overhauser enhancement spectroscopy spectra have been recorded by using a new procedure and with different mixing times and different spectral widths, in order to observed dipolar connectivities between isotropically shifted and diamagnetic resonances. Scalar-correlated spectra were obtained with both correlation spectroscopy and total-correlation spectroscopy experiments. The original X-ray coordinates of the Cu, Co enzyme were used in order to obtain the appropriate interproton distances. The data allowed us to assign more than 20 new resonances to protons which cover a wide region around the copper ion, including two protons belonging to the catalytically important Arg141 residue. The results represent significant progress in the effort to elucidate the three-dimensional features of the region surrounding the active site of Cu, Zn superoxide dismutase in solution and a tool for a deeper investigation of the reaction mechanism of the enzyme.

1H, 13C and 15N random coil NMR chemical shifts of the common amino acids. I. Investigations of nearest-neighbor effects

In this study we report on the 1H, 13C and 15N NMR chemical shifts for the random coil state and nearest-neighbor sequence effects measured from the protected linear hexapeptide Gly-Gly-X-Y-Gly-Gly (where X and Y are any of the 20 common amino acids). We present data for a set of 40 peptides (of the possible 400) including Gly-Gly-X-Ala-Gly-Gly and Gly-Gly-X-Pro-Gly-Gly, measured under identical aqueous conditions. Because all spectra were collected under identical experimental conditions, the data from the Gly-Gly-X-Ala-Gly-Gly series provide a complete and internally consistent set of 1H, 13C and 15N random coil chemical shifts for all 20 common amino acids. In addition, studies were also conducted into nearest-neighbor effects on the random coil shift arising from a variety of X and Y positional substitutions. Comparisons between the chemical shift measurements obtained from Gly-Gly-X-Ala-Gly-Gly and Gly-Gly-X-Pro-Gly-Gly reveal significant systematic shift differences arising from the presence of proline in the peptide sequence. Similarly, measurements of the chemical shift changes occurring for both alanine and proline (i.e., the residues in the Y position) are found to depend strongly on the type of amino acid substituted into the X position. These data lend support to the hypothesis that sequence effects play a significant role in determining peptide and protein chemical shifts.

Phosphophoryn, an "acidic" biomineralization regulatory protein: conformational folding in the presence of Cd(II)

The divalent cation-induced protein folding properties of the template macromolecule, bovine dentine phosphophoryn (BDPP), have been examined by 1H/31P/13C/113Cd-nmr spectroscopy. Cd(II) was employed, exploiting the sensitivity of 113Cd-nmr to ligand-binding interactions and kinetics. Cation binding was studied over the stoichiometric range of 0-50: 1 Cd(II): protein (mole ratio), well below the range of Cd(II) concentration required to induce protein precipitation. The stepwise titration of divalent cation-depleted phosphophoryn at pH 7.2 in H2O/D2O with 113CdCl2 revealed that (PSer)n, (PSerAsp)n, and (Asp)n polyelectrolyte cation-binding domains undergo two major transitions in their secondary and tertiary structures: the first transition, occurring between 1:10 and 1:1 Cd(II): protein stoichiometry, and the second, between 10:1 and 50:1. By monitoring the amide NH intensities, 31P-nmr chemical shift, and 13C Asp-C, resonances, it was concluded that Cd(II) ions exhibit a binding-site preference for polyelectrolyte cation-binding domains, in the order (PSer)n > (PSerAsp)n > (Asp)n This preference correlates with the degree of negative charge density for each sequence motif. Accompanying the backbone conformational transitions at the polyelectrolyte regions were conformational transitions in the flanking hinge domains, indicating that the hinge domains participate in the folding of the phosphophoryn molecule as divalent cation binding occurs at the polyelectrolyte domains. We were unsuccessful in detecting phosphophoryn-bound Cd(II) species by 113Cd-nmr because of chemical exchange modulation. However, using a smaller 21-residue peptide mimetic of phosphophoryn, we have observed three stoichiometric-dependent 113Cd resonances that differ in terms of the oxoanion coordination number. Our observation of multiple Cd(II) species in the presence of the peptide supports our contention that Cd(II) has many chemically distinct coordination sites on phosphophoryn, each in multiple equilibria with H2O, Cl-, and side-chain oxoatoms.

1H-NMR assignments and secondary structure of human plasminogen kringle 1

The 1H-NMR spectrum of the kringle 1 domain of human plasminogen complexed with 6-aminohexanoic acid, an antifibrinolytic drug, has been assigned. Elements of secondary structure have been identified on the basis of sequential, medium and long-range dipolar interactions, back-bone amide spin-spin couplings (3JHN-H alpha) and 1H-2H exchange rates. The kringle contains scarcely any repetitive secondary structure: eight reverse turns and two short beta-sheets. These comprise 40% and 12% of the domain, respectively. No alpha-helix was found. An aromatic cluster formed by His31, Phe36, Trp62, Phe64, Tyr72 and Tyr74 is indicated by several inter-residue Overhauser connectivities. Contacts between the methyl groups of Leu46 and the side chains of Phe36, Trp62 and Trp25 are observed. A second hydrophobic cluster formed by Tyr9, Ile77 and Leu78 is also indicated. A comparison of secondary structure elements among plasminogen kringles 1 and 4 and tissue-type plasminogen activator kringle 2 suggests that there is variability in the position and number of reverse turns on going from one kringle to another; however, the beta-sheets are conserved among the homologs.

Assignment of the 13C and 13CO resonances for Rhodobacter capsulatus ferrocytochrome c2 using double-resonance and triple-resonance NMR spectroscopy

Rhodobacter capsulatus cytochrome c2 uniformly labelled with 13C/15N has been prepared. The 13C resonances of the reduced state, including those of the carbonyl and heme 13C, have been assigned using a combination of various two- and three-dimensional correlated NMR experiments. Assignment of the sidechain 13C resonances facilitated correction of a small number of previously misassigned sidechain 1H and led to the additional assignment of 32 1H. It was found that 13C alpha and 13CO secondary shifts were better indicators of secondary structure than 1H alpha and 13C beta secondary shifts. Moreover, it was demonstrated that, despite the significant ring current effects present in heme proteins, 13C alpha and 13CO secondary shifts can be employed to accurately identify secondary structure in heme proteins, independently of NOE experiments.

Phosphophoryn, a biomineralization template protein: pH-dependent protein folding experiments

The protein folding behavior of a polyelectrolyte protein, bovine dentine phosphophoryn (BDPP), in the pH range of 1.82-11.0 has been investigated. One- and two-dimensional nmr spectroscopy has been utilized to obtain proton spin assignments for amino acid residues in D2O and in H2O. One-dimensional 31P-nmr experiments verify the existence of three separate classes of O-phosphoserine (PSer) resonances in BDPP (alpha, beta, chi), representing three distinct PSer residue populations at pH 6.94. By means of pH titration and 1H-nmr, five populations of Asp residues can be identified. Three of these populations exhibit secondary inflection points on their pH titration curves that correspond to an observed pKa of 6.17-6.95. The presence or absence of secondary inflection points for Asp populations and the 31P-nmr chemical shift dispersion for the three PSer residue populations indicate that BDPP may be comprised of homologous (Asp-Asp)n. (PSer-PSer)n, and heterologous (PSer-Asp)n sequences arranged into polyelectrolyte cluster regions. The pH titration also revealed that certain populations of Ser, Gly, and Pro residues in BDPP exhibit pH-dependent resonance frequency shifts. The "apparent" pKa for the transition points of these frequency shifts corresponds to either the pK1a of Pser monophosphate ester and/or the pKa of Asp COOH group of BDPP polyelectrolyte regions. On the basis of these transition points, we can assign four types of Ser, Gly, or Pro-containing "intervening" regions in BDPP, based on their sensitivity to protonation and deprotonation events occurring at (Asp)n, (PSer)n, or (PSer-Asp)n anionic cluster regions that flank the intervening regions. Our 1H-nmr experiments also reveal that BDPP assumes a folded conformation at low pH. As the pH increases, this conformation undergoes several unfolding transitions as the BDPP molecule assumes more open conformations in response to increased electrostatic repulsion between polyelectrolyte anionic regions in the protein. These folding-unfolding transitions are mediated by the intervening regions, which act as "hinges" to allow the polyelectrolyte regions to fold relative to one another.

pH-induced structural changes in human serum apotransferrin. pKa values of histidine residues and N-terminal amino group determined by 1H-NMR spectroscopy

The binding of apotransferrin (80 kDa) to the transferrin receptor is known to be highly pH-dependent. We have investigated pH-induced structural changes in human serum apotransferrin over the pH* (meter reading in D2O solutions) range 2.5-11 using 1H-NMR spectroscopy. The pKa values of 14 (possibly 15) of the 19 His residues in the protein have been determined as well as that of the terminal amino group (Val1, 7.75). About eight His residues deprotonate when the pH* is raised from the endosomal value of about 5.5 to the blood plasma value (7.4). Four His residues have pKa < 6. Sharp discontinuities in the His titration curves were observed below pH 4.3 and at pH 3.5 molten globule states were detected.

Proton NMR studies of bovine serum albumin. Assignment of spin systems

A variety of one- and two-dimensional 1H-NMR methods have been applied to the study of defatted 66.5-kDa bovine serum albumin in solution. 1. The majority of the protons gave rise to broad unresolved resonances and spectral enhancement methods for one-dimensional spectra were investigated in detail. A combination of exponential and sine-bell functions was particularly effective. 2. The presence of contaminating glycoproteins in some commercial samples of bovine serum albumin was readily detectable from their N-acetyl resonances at about 2.1 ppm. 3. The release of bound Cys (from mixed disulphide at Cys34) was observed after addition of dithiothreitol. 4. Through the use of two-dimensional shift-correlated spectroscopy, assignments of some 80 spin systems to amino acid type were made. 5. The pKa of the N-terminal Asp was measured as 7.8 (0.1 M phosphate buffer, 310 K). 6. 1H NMR spectra of bovine, human, porcine and rat serum albumins have been compared. Using sequence comparisons, specific assignments have been made for the N-terminal residues of bovine (Asp-Thr-His), human (Asp-Ala-His), porcine (Asp-Thr-Tyr) and rat (Glu-Ala-His) albumins, and for Thr189, Tyr155 and His59/377 of bovine albumin. 7. These NMR data suggest that certain local regions of bovine serum albumin are highly mobile yet structured in solution, and demonstrate that the application of both one- and and two-dimensional NMR methods will allow more detailed investigations of structural transitions in serum albumins induced by, for example, pH, drug and metal binding.

Hydrophobic clustering in nonnative states of a protein: interpretation of chemical shifts in NMR spectra of denatured states of lysozyme

Chemical shifts of resonances of specific protons in the 1H NMR spectrum of thermally denatured hen lysozyme have been determined by exchange correlation with assigned native state resonances in 2D NOESY spectra obtained under conditions where the two states are interconverting. There are subtle but widespread deviations of the measured shifts from the values which would be anticipated for a random coil; in the case of side chain protons these are virtually all net upfield shifts and it is shown that this may be the averaged effect of interactions with aromatic rings in a partially collapsed denatured state. In a very few cases, notably that of two sequential tryptophan residues, it is possible to interpret these effects in terms of specific, local interresidue interactions. Generally, however, there is no correlation with either native state shift perturbations or with sequence proximity to aromatic groups. Diminution of most of the residual shift perturbations on reduction of the disulfide cross-links confirms that they are not simply effects of residues adjacent in the sequence. Similar effects of chemical denaturants, with the disulfides intact, demonstrate that the shift perturbations reflect an enhanced tendency to side chain clustering in the thermally denatured state. The temperature dependences of the shift perturbations suggest that this clustering is noncooperative and is driven by small, favorable enthalpy changes. While the extent of conformational averaging is clearly much greater than that observed for a homologous protein, alpha-lactalbumin, in its partially folded "molten globule" state, the results clearly show that thermally denatured lysozyme differs substantially from a random coil, principally in that it is partially hydrophobically collapsed.

Tertiary structure of two-electron reduced Megasphaera elsdenii flavodoxin and some implications, as determined by two-dimensional 1H-NMR and restrained molecular dynamics

The tertiary structure of the non-crystallizable two-electron-reduced Megasphaera elsdenii flavodoxin (15 kDa, 137 amino acid residues) has been determined using nuclear Overhauser enhancement restraints extracted from two-dimensional 1H-NMR spectra. A tertiary structure satisfying the experimental restraints very well (maximum NOE violation of 66 pm) was obtained with use of restrained molecular dynamics, using 509 distance restraints (including one non-NOE) on a starting structure modeled from the crystal structure of one-electron-reduced Clostridium MP flavodoxin. The protein consists of a central parallel beta-sheet surrounded on both sides by two alpha-helices. The flavin is positioned at the periphery of the molecule. The tertiary structure of the protein is highly defined with the exception of the flavin. The latter is expected to result from performing the restrained molecular dynamics simulation without water molecules and without proper charges on the flavin. The flavin, including the phosphate, the ribityl side chain and the isoalloxazine ring, is solvent accessible under the experimental conditions used and evidenced by a two-dimensional amide exchange experiment. This accessibility is expected to be important in the redox potential regulation of the semiquinone/hydroquinone couple of the protein. The amide exchange against deuterons and several typical line shapes in the two-dimensional NMR spectra are consistent with the structure generated. The structure is discussed in detail.

Secondary-structure dependent chemical shifts in proteins

Chemical shift data have been collected on eight proteins that have the same conformation in solution as in their crystal structures. Ring-current shifts have been calculated and subtracted from the exerimentally measured shifts, to leave shifts that depend only on local conformation. Overall, the shifts show an approximately normal distribution with no appreciable skewness, thus confirming that ring-current shifts have the overall effect of skewing the distribution to high field. In helices, NH and C(alpha)H have a high significant tendency to resonate to high field, whereas they resonate to low field in beta-sheets. Side-chain protons resonate slightly to high field in beta-sheets. Chemical shift distributions are narrowest for side-chain protons, and widest for amide protons. When only slowly exchanging amide protons are considered, the high field shift for amide protons in helices is more pronounced, but there is only a small difference in sheets. C(alpha)H signals at the N-terminal end of helices tend to resonate to higher field than those at the C-terminal end, whereas for NH signals it is the C-terminal end that resonates to higher field. There is no significant effect of position within the helix on side-chain signals, implying that the helix dipole has little effect on shifts within the helix.

A two-dimensional 1H NMR study on Megasphaera elsdenii flavodoxin in the reduced state. Sequential assignments

Assignments for the 137 amino acid residues of Megasphaera elsdenii flavodoxin in the reduced state have been made using the sequential resonance assignment procedure. Several hydroxyl and sulfhydryl protons were observed at 41 degrees C at pH 8.3. Spin systems were sequentially assigned using phase-sensitive two-dimensional-correlated spectroscopy and phase-sensitive nuclear Overhauser enhancement spectroscopy. Spectra of the protein in H2O and of protein preparations either completely or partly exchanged against 2H2O were obtained. Use of the fast electron shuttle between the paramagnetic semiquinone and the diamagnetic hydroquinone state greatly simplified the NMR spectra, making it possible to assign easily the 1H resonances of amino acid residues located in the immediate neighbourhood of the isoalloxazine ring. The majority of the nuclear Overhauser effect contracts between the flavin and the apoprotein correspond to the crystal structure of the flavin domain of Clostridium MP flavodoxin, but differences are also observed. The assignments provide the basis for the structure determination of M. elsdenii flavodoxin in the reduced state as well as for assigning the resonances of the oxidized flavodoxin.