A new conformational database potential involving dihedral angle relationships in databases of high-resolution highly refined protein crystal structures is presented as a method for improving the quality of structures generated from NMR data. The rationale for this procedure is based on the observation that uncertainties in the description of the nonbonded contacts present a key limiting factor in the attainable accuracy of protein NMR structures and that the nonbonded interaction terms presently used have poor discriminatory power between highand low-probability local conformations. The idea behind the conformational database potential is to restrict sampling during simulated annealing refinement to conformations that are likely to be energetically possible by effectively limiting the choices of dihedral angles to those that are known to be physically realizable. In this manner, the variability in the structures produced by this method is primarily a function of the experimental restraints, rather than an artifact of a poor nonbonded interaction model. We tested this approach with the experimental NMR data (comprising an average of about 30 restraints per residue and consisting of interproton distances, torsion angles, 'JHNa coupling constants, and I3C chemical shifts) used previously to calculate the solution structure of reduced human thioredoxin (Qin J, Clore GM, Gronenborn AM, 1994, Structure 2503-522). Incorporation of the conformational database potential into the target function used for refinement (which also includes terms for the experimental restraints, covalent geometry, and nonbonded interactions in the form of either a repulsive, repulsive-attractive, or 6-12 Lennard-Jones potential) results in a significant improvement in various quantitative measures of quality (Ramachandran plot, side-chain torsion angles, overall packing). This is achieved without compromising the agreement with the experimental restraints and the deviations from idealized covalent geometry that remain within experimental error, and the agreement between calculated and observed ' H chemical shifts that provides an independent NMR parameter of accuracy. The method is equally applicable to crystallographic refinement, and should be particular useful during the early stages of either an NMR or crystallographic structure determination and in cases where relatively few experimental restraints can be derived from the measured data (due, for example, to broad lines in the NMR spectra or to poorly diffracting crystals).Keywords: conformational databases; NMR; protein structure determination; refinement; X-rayThe determination and refinement of a protein structure by NMR or crystallography depends on the minimization of a target function (Etot) comprising the experimental restraints (E,,,), and two groups of a priori restraints consisting of covalent geometry (Ecou) and nonbonded interactions (Enb) (Jack &

No abstract

No abstract

The solution conformations of the protein hirudin have been investigated by the combined use of distance geometry and restrained molecular dynamics calculations. The basis for the structure determination comprised 359 approximate interproton distance restraints and 10 4 backbone torsion angle restraints derived from n.m.r. measurements. It is shown that hirudin is composed of three domains: a central core made up of residues 3-30, 37-46 and 56-57; a protruding 'finger' (residues 31-36) consisting of the tip of an antiparaliel 3 sheet, and an exposed loop (residues 47-55). The structure of each individual domain is relatively well defined with average backbone atomic r.m.s. differences of <2 A between the final seven converged restrained dynamic structures and the mean structure obtained by averaging their coordinates. The orientation of the two minor domains relative to the central core, however, could not be determined as no long-range (I i-jI >5) interdomain proton-proton contacts could be observed in the two-dimensional nuclear Overhauser enhancement spectra. From the restrained molecular dynamics calculations it appears that the two minor domains exhibit large rigid-body motions relative to the central core.