William J. Wedemeyer

Recent improvements in prediction of protein structure by global optimization of a potential energy function.

Recent improvements of a hierarchical ab initio or de novo approach for predicting both α and β structures of proteins are described. The united-residue energy function used in this procedure includes multibody interactions from a cumulant expansion of the free energy of polypeptide chains, with their relative weights determined by Z-score optimization. The critical initial stage of the hierarchical procedure involves a search of conformational space by the conformational space annealing (CSA) method, followed by optimization of an all-atom model. The procedure was assessed in a recent blind test of protein structure prediction (CASP4). The resulting lowest-energy structures of the target proteins (ranging in size from 70 to 244 residues) agreed with the experimental structures in many respects. The entire experimental structure of a cyclic α-helical protein of 70 residues was predicted to within 4.3 Å α-carbon (Cα) rms deviation (rmsd) whereas, for other α-helical proteins, fragments of roughly 60 residues were predicted to within 6.0 Å Cα rmsd. Whereas β structures can now be predicted with the new procedure, the success rate for α/β- and β-proteins is lower than that for α-proteins at present. For the β portions of α/β structures, the Cα rmsds are less than 6.0 Å for contiguous fragments of 30–40 residues; for one target, three fragments (of length 10, 23, and 28 residues, respectively) formed a compact part of the tertiary structure with a Cα rmsd less than 6.0 Å. Overall, these results constitute an important step toward the ab initio prediction of protein structure solely from the amino acid sequence.

Exact analytical loop closure in proteins using polynomial equations

Loop closure in proteins has been studied actively for over 25 years. Using spherical geometry and polynomial equations, several loop‐closure problems in proteins are solved exactly by reducing them to the determination of the real roots of a polynomial. Loops of seven, eight, and nine atoms are treated explicitly, including the tripeptide and disulfide‐bonded loop‐closure problems. The number of valid loop closures can be evaluated by the method of Sturm chains, which counts the number of real roots of a polynomial. Longer loops can be treated by three methods: by sampling enough dihedral angles to reduce the problem to a soluble loop‐closure problem; by applying the loop‐closure algorithm hierarchically; or by decimating the chain into independently moving rigid elements that can be reconnected using loop‐closure algorithms. Applications of the methods to docking, homology modeling and NMR problems are discussed. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 819–844, 1999

Structural determinants of oxidative folding in proteins

A method for determining the kinetic fate of structured disulfide species (i.e., whether they are preferentially oxidized or reshuffle back to an unstructured disulfide species) is introduced. The method relies on the sensitivity of unstructured disulfide species to low concentrations of reducing agents. Because a structured des species that preferentially reshuffles generally first rearranges to an unstructured species, a small concentration of reduced DTT (e.g., 260 μM) suffices to distinguish on-pathway intermediates from dead-end species. We apply this method to the oxidative folding of bovine pancreatic ribonuclease A (RNase A) and show that des[40–95] and des[65–72] are productive intermediates, whereas des[26–84] and des[58–110] are metastable dead-end species that preferentially reshuffle. The key factor in determining the kinetic fate of these des species is the relative accessibility of both their thiol groups and disulfide bonds. Productive intermediates tend to be disulfide-secure, meaning that their structural fluctuations preferentially expose their thiol groups, while keeping their disulfide bonds buried. By contrast, dead-end species tend to be disulfide-insecure, in that their structural fluctuations expose their disulfide bonds in concert with their thiol groups. This distinction leads to four generic types of oxidative folding pathways. We combine these results with those of earlier studies to suggest a general three-stage model of oxidative folding of RNase A and other single-domain proteins with multiple disulfide bonds.

A role for intermolecular disulfide bonds in prion diseases

The key event in prion diseases seems to be the conversion of the prion protein PrP from its normal cellular isoform (PrPC) to an aberrant “scrapie” isoform (PrPSc). Earlier studies have detected no covalent modification in the scrapie isoform and have concluded that the PrPC → PrPSc conversion is a purely conformational transition involving no chemical reactions. However, a reexamination of the available biochemical data suggests that the PrPC → PrPSc conversion also involves a covalent reaction of the (sole) intramolecular disulfide bond of PrPC. Specifically, the data are consistent with the hypothesis that infectious prions are composed of PrPSc polymers linked by intermolecular disulfide bonds. Thus, the PrPC → PrPSc conversion may involve not only a conformational transition but also a thiol/disulfide exchange reaction between the terminal thiolate of such a PrPSc polymer and the disulfide bond of a PrPC monomer. This hypothesis seems to account for several unusual features of prion diseases.

Evolution of physics-based methodology for exploring the conformational energy landscape of proteins

The evolution of our physics‐based computational methods for determining protein conformation without the introduction of secondary‐structure predictions, homology modeling, threading, or fragment coupling is described. Initial use of a hard‐sphere potential captured much of the structural properties of polypeptide chains, and subsequent more refined force fields, together with efficient methods of global optimization provide indications that progress is being made toward an understanding of the interresidue interactions that underlie protein folding.

Acceleration of oxidative folding of bovine pancreatic ribonuclease A by anion-induced stabilization and formation of structured native-like intermediates

Phosphate anions accelerate the oxidative folding of reduced bovine pancreatic ribonuclease A with dithiothreitol at several temperatures and ionic strengths. The addition of 400 mM phosphate at pH 8.1 increased the regeneration rate of native protein 2.5‐fold at 15°C, 3.5‐fold at 25°C, and 20‐fold at 37°C, compared to the rate in the absence of phosphate. In addition, the effects of other ions on the oxidative folding of RNase A were examined. Fluoride was found to accelerate the formation of native protein under the same oxidizing conditions. In contrast, cations of high charge density or ions with low charge density appear to have an opposite effect on the folding of RNase A. The catalysis of oxidative folding results largely from an anion‐dependent stabilization and formation of tertiary structure in productive disulfide intermediates (des‐species). Phosphate and fluoride also accelerate the initial equilibration of unstructured disulfide ensembles, presumably due to non‐specific electrostatic and hydrogen bonding effects on the protein and solvent.