Search

A new computational method for simultaneously estimating all the proton-proton coupling constants in a molecule from COSY spectra [Yang, J.-X. and Havel, T.F. (1994) J. Biomol. NMR, 4, 807-826] is applied to experimental data from two polypeptides. The first of these is a cyclic hexapeptide denoted as VDA (-D-Ala1-Phe2-Trp3-Lys(Z)4-Val5-Phe6-), in deuterated DMSO, while the second is a 39-residue protein, called decorsin, in aqueous solution. The effect of different data processing strategies and different initial parameter values on the accuracy of the coupling constants was explored. In the case of VDA, most of the coupling constants did not depend strongly on the initial values chosen for the optimization or on how the data were processed. This, together with our previous experience using simulated data, implies strongly that these values are accurate estimates of the coupling constants. They also differ by an average of only 0.36 Hz from the values of the 14 coupling constants that could be measured independently by established methods. In the case of decorsin, many of the coupling constants exhibited a moderate dependence on their initial values and a strong dependence on how the data were processed. With the most successful data processing strategy, the amide-alpha coupling constants differed by an average of 1.11 Hz from the 21 values that could be measured by established methods, while two thirds of the three-bond coupling constants fell within 1.0 Hz of the ranges obtained by applying the Karplus relation to an independently computed ensemble of distance geometry structures. The averages of the coupling constants over multiple optimizations using random initial values were computed in order to obtain the best possible estimates of the coupling constants. Most clearly incorrect averages can be identified by large standard deviations in the coupling constants or the associated line widths and chemical shifts, and can be explained by strong coupling and/or overlap with the water signal, the diagonal peaks or other cross peaks.

A method is proposed for defining a probability distribution on an ensemble of protein conformations from a 2D NOE spectrum, while at the same time back-calculating the experimental spectrum from the ensemble. This enables one to assess the relative quality and significance of the conformations, and to test the consistency of the ensemble as a whole with the experimental spectrum. The method eliminates the need to integrate the cross-peak intensities and is surprisingly insensitive to random noise in the spectrum. In this communication, these advantages are demonstrated by applying the method to simulated data, for which the correct result is already known.

A computational method is described that takes an initial estimate of the chemical shifts, line widths and scalar coupling constants for the protons in a molecule, and refines this estimate so as to improve the least-squares fit between an experimental COSY spectrum and the spectrum simulated from these parameters in the weak-coupling approximation. In order to evaluate the potential of such refinements for estimating these parameters from COSY experiments, the method has been applied to a large number of sample problems which were themselves simulated from standard conformations of the amino acids, along with 25 near-native conformations of the protein bovine pancreatic trypsin inhibitor. The results of this evaluation show that: (i) if the chemical shifts are known to within ca. 0.01 ppm and no noise or artifacts are present in the data, the method is capable of recovering the correct coupling constants, starting from essentially arbitrary values, to within 0.1 Hz in almost all cases. (ii) Although the precision of these estimates of the coupling constants is degraded by the limited resolution, noise and artifacts present in most experimental spectra, the large majority of coupling constants can still be recovered to within 1.0 Hz; the local minimum problem is not made significantly worse by such defects in the data. (iii) The method assigns an 'effective' line width to all the resonances, and in the process can resolve overlapping cross peaks. (iv) The method is not capable of determining the chemical shifts a priori, due to the presence of numerous local minima in the least-squares residual as a function of these parameters.