Local interactions and the optimization of protein folding

The role of local interactions in protein folding has recently been the subject of some controversy. Here we investigate an extension of Zwanzig's simple and general model of folding in which local and nonlocal interactions are represented by functions of single and multiple conformational degrees of freedom, respectively. The kinetics and thermodynamics of folding are studied for a series of energy functions in which the energy of the native structure is fixed, but the relative contributions of local and nonlocal interactions to this energy are varied over a broad range. For funnel shaped energy landscapes, we find that 1) the rate of folding increases, but the stability of the folded state decreases, as the contribution of local interactions to the energy of the native structure increases, and 2) the amount of native structure in the unfolded state and the transition state vary considerably with the local interaction strength. Simple exponential kinetics and a well-defined free energy barrier separating folded and unfolded states are observed when nonlocal interactions make an appreciable contribution to the energy of the native structure; in such cases a transition state theory type approximation yields reasonably accurate estimates of the folding rate. Bumps in the folding funnel near the native state, which could result from desolvation effects, side chain freezing, or the breaking of nonnative contacts, significantly alter the dependence of the folding rate on the local interaction strength: the rate of folding decreases when the local interaction strength is increased beyond a certain point. A survey of the distribution of strong contacts in the protein structure database suggests that evolutionary optimization has involved both kinetics and thermodynamics: strong contacts are enriched at both very short and very long sequence separations. Proteins 29:282–291, 1997. © 1997 Wiley-Liss, Inc.     

Perturbations of the denatured state ensemble: modeling their effects on protein stability and folding kinetics.

By considering the denatured state of a protein as an ensemble of conformations with varying numbers of sequence-specific interactions, the effects on stability, folding kinetics, and aggregation of perturbing these interactions can be predicted from changes in the molecular partition function. From general considerations, the following conclusions are drawn: (1) A perturbation that enhances a native interaction in denatured state conformations always increases the stability of the native state. (2) A perturbation that promotes a non-native interaction in the denatured state always decreases the stability of the native state. (3) A change in the denatured state ensemble can alter the kinetics of aggregation and folding. (4) The loss (or increase) in stability accompanying two mutations, each of which lowers (or raises) the free energy of the denatured state, will be less than the sum of the effects of the single mutations, except in cases where both mutations affect the same set of partially folded conformations. By modeling the denatured state as the ensemble of all non-native conformations of hydrophobic-polar (HP) chains configured on a square lattice, it can be shown that the stabilization obtained from enhancement of native interactions derives in large measure from the avoidance of non-native interactions in the D state. In addition, the kinetic effects of fixing single native contacts in the denatured state or imposing linear gradients in the HH contact probabilities are found, for some sequences, to significantly enhance the efficiency of folding by a simple hydrophobic zippering algorithm. Again, the dominant mechanism appears to be avoidance of non-native interactions. These results suggest stabilization of native interactions and imposition of gradients in the stability of local structure are two plausible mechanisms involving the denatured state that could play a role in the evolution of protein folding and stability.     

A generalized knowledge‐based discriminatory function for biomolecular interactions

Several novel and established knowledge‐based discriminatory function formulations and reference state derivations have been evaluated to identify parameter sets capable of distinguishing native and near‐native biomolecular interactions from incorrect ones. We developed the r·m·r function, a novel atomic level radial distribution function with mean reference state that averages over all pairwise atom types from a reduced atom type composition, using experimentally determined intermolecular complexes in the Cambridge Structural Database (CSD) and the Protein Data Bank (PDB) as the information sources. We demonstrate that r·m·r had the best discriminatory accuracy and power for protein‐small molecule and protein‐DNA interactions, regardless of whether the native complex was included or excluded, from the test set. The superior performance of the r·m·r discriminatory function compared with seventeen alternative functions evaluated on publicly available test sets for protein‐small molecule and protein‐DNA interactions indicated that the function was not over optimized through back testing on a single class of biomolecular interactions. The initial success of the reduced composition and superior performance with the CSD as the distribution set over the PDB implies that further improvements and generality of the function are possible by deriving probabilities from subsets of the CSD, using structures that consist of only the atom types to be considered for given biomolecular interactions. The method is available as a web server module at http://protinfo.compbio.washington.edu . Proteins 2009. © 2008 Wiley‐Liss, Inc.     

The presence of non‐native helical structure in the unfolding of a beta‐sheet protein MPT63

MPT63, a major secreted protein from Mycobacterium tuberculosis, has been shown to have immunogenic properties and has been implicated in virulence. MPT63 is a β‐sandwich protein containing 11 β strands and a very short stretch of 3₁₀ helix. The detailed experimental and computational study reported here investigates the equilibrium unfolding transition of MPT63. It is shown that in spite of being a complete β‐sheet protein, MPT63 has a strong propensity toward helix structures in its early intermediates. Far UV‐CD and FTIR spectra clearly suggest that the low‐pH intermediate of MTP63 has enhanced helical content, while fluorescence correlation spectroscopy suggests a significant contraction. Molecular dynamics simulation complements the experimental results indicating that the unfolded state of MPT63 traverses through intermediate forms with increased helical characteristics. It is found that this early intermediate contains exposed hydrophobic surface, and is aggregation prone. Although MPT63 is a complete β‐sheet protein in its native form, the present findings suggest that the secondary structure preferences of the local interactions in early folding pathway may not always follow the native conformation. Furthermore, the Gly25Ala mutant supports the proposed hypothesis by increasing the non‐native helical propensity of the protein structure.     

Extending the folding nucleus of ubiquitin with an independently folding beta-hairpin finger: hurdles to rapid folding arising from the stabilisation of local interactions

The N-terminal beta-hairpin sequence of ubiquitin has been implicated as a folding nucleation site. To extend and stabilise the ubiquitin folding nucleus, we have inserted an autonomously folding 14-residue peptide sequence beta4 which in isolation forms a highly populated beta-hairpin (>70%) stabilised by local interactions. NMR structural analysis of the ubiquitin mutant (Ubeta4) shows that the hairpin finger is fully structured and stabilises ubiquitin by approximately 8kJmol(-1). Protein engineering and kinetic (phi(F)-value) analysis of a series of Ubeta4 mutants shows that the hairpin extension of Ubeta4 is also significantly populated in the transition state (phi(F)-values >0.7) and has the effect of templating the formation of native contacts in the folding nucleus of ubiquitin. However, at low denaturant concentrations the chevron plot of Ubeta4 shows a small deviation from linearity (roll-over effect), indicative of the population of a compact collapsed state, which appears to arise from over-stabilisation of local interactions. Destabilising mutations within the native hairpin sequence and within the engineered hairpin extension, but not elsewhere, eliminate this non-linearity and restore apparent two-state behaviour. The pitfall to stabilising local interactions is to present hurdles to the rapid and efficient folding of small proteins down a smooth folding funnel by trapping partially folded or misfolded states that must unfold or rearrange before refolding.     

Dissecting the roles of individual interactions in protein stability: lessons from a circularized protein

A circular form of bovine pancreatic trypsin inhibitor (BPTI) has been prepared by introducing a peptide bond between the N- and C-termini, which are in close proximity in the native conformation. The pathway and energetics of the disulphide-coupled folding transition of the circular protein have been studied using methods applied previously to the unmodified protein. The cross-link between the termini was found not to significantly stabilize the native state in spite of the expected reduction in entropy of the unfolded protein. This unexpected result has led to a reexamination of the stabilization expected from a cross-link, considering effects on the native, as well as unfolded, states of the protein. The greatest stabilization is expected when the cross-linked groups are held rigidly in the native protein in the optimum orientation for forming the cross-link. Similar analyses, utilizing thermodynamic cycles, can be applied to other interactions that stabilize native proteins, including disulphide bonds, salt bridges, and hydrogen bonds and to modifications to the protein that remove them. In general, the contribution of an individual interaction to the stability of the native state depends on the extent to which the interaction is favored in the native conformation, which can vary greatly depending on the local environment of the interacting groups.     

Engineering stabilising beta-sheet interactions into a conformationally flexible region of the folding transition state of ubiquitin

Protein engineering studies suggest that the transition state for the folding of ubiquitin is highly polarised towards the N-terminal part of the sequence and involves a nucleus of residues within the beta-hairpin (residues 1-17) and main alpha-helix (residues 23-34). In contrast, the observation of small phi-values for residues in the C-terminal portion of the sequence (residues 35-76), coupled with a folding topology that results in a much higher contact order, suggests that fast folding of ubiquitin is dependent upon configurational flexibility in the C-terminal part of the polypeptide chain to ensure passage down a relatively smooth folding funnel to the native state. We show that the introduction of a small mini-hairpin motif as an extension of the native 43-50 hairpin stabilises local interactions in the C-terminal part of the sequence, resulting largely in a deceleration of the unfolding kinetics without perturbing the apparent two-state folding mechanism. However, a single-point Leu-->Phe substitution within the engineered hairpin sequence leads to the premature collapse of the denatured ensemble through the stabilisation of non-native interactions and the population of a compact intermediate. Non-linear effects in the kinetic data at low concentrations of denaturant suggest that the collapsed state, which is further stabilised in the presence of cosmotropic salts, may subsequently fold directly to the native state through a "triangular" reaction scheme involving internal rearrangement rather than unfolding and refolding.     

Improved Gō-like models demonstrate the robustness of protein folding mechanisms towards non-native interactions

The use of simple theoretical models has provided a considerable contribution to our present understanding of the means by which proteins adopt their native fold from the plethora of available unfolded states. A common assumption in building computationally tractable models has been the neglect of stabilizing non-native interactions in the class of models described as "Gō-like." The focus of this study is the characterization of the folding of a number of proteins via a Gō-like model, which aims to map a maximal amount of information reflecting the protein sequence onto a "minimalist" skeleton. This model is shown to contain sufficient information to reproduce the folding transition states of a number of proteins, including topologically analogous proteins that fold via different transition states. Remarkably, these models also demonstrate consistency with the general features of folding transition states thought to be stabilized by non-native interactions. This suggests that native interactions are the primary determinant of most protein folding transition states, and that non-native interactions lead only to local structural perturbations. A prediction is also included for an asymmetrical folding transition state of bacteriophage lambda protein W, which has yet to be subjected to experimental characterization.     

Thermodynamics and kinetics of non-native interactions in protein folding: a single point mutant significantly stabilizes the N-terminal domain of L9 by modulating non-native interactions in the denatured state

Comparatively little is known about the role of non-native interactions in protein folding and their role in both folding and stability is controversial. We demonstrate that non-native electrostatic interactions involving specific residues in the denatured state can have a significant effect upon protein stability and can persist in the transition state for folding. Mutation of a single surface exposed residue, Lys12 to Met, in the N-terminal domain of the ribosomal protein L9 (NTL9), significantly increased the stability of the protein and led to faster folding. Structural and energetic studies of the wild-type and K12M mutant show that the 1.9 kcal mol(-1) increase in stability is not due to native state effects, but rather is caused by modulation of specific non-native electrostatic interactions in the denatured state. pH dependent stability measurements confirm that the increased stability of the K12M is due to the elimination of favorable non-native interactions in the denatured state. Kinetic studies show that the non-native electrostatic interactions involving K12 persist in the transition state. The analysis demonstrates that canonical Phi-values can arise from the disruption of non-native interactions as well as from the development of native interactions.     

Local control of peptide conformation: Stabilization of cis proline peptide bonds by aromatic proline interactions

In the native state of proteins there is a marked tendency for an aromatic amino acid to precede a cis proline. There are also significant differences between the three aromatic amino acids with Tyr exhibiting a noticeably higher propensity than Phe or Trp to precede a cis proline residue. In order to study the role that local interactions play in these conformation preferences, a set of tetrapeptides of the general sequence acetyl-Gly-X-Pro-Gly-carboxamide (GXPG), where X = Tyr, Phe, Trp, Ala, or cyclohexyl alanine, were synthesized and studied by nmr. Analysis of the nmr data shows that none of the peptides adopt a specific backbone structure. Ring current shifts, the equilibrium constants, the Van't Hoff enthalpy, and the measured rate of cis-trans isomerization all indicate that the cis proline conformer is stabilized by favorable interactions between the aromatic ring and the proline residue. Analysis of the side chain conformation of the aromatic residue and analysis of the chemical shifts of the pyrrolidine ring protons shows that the aromatic side chain adopts a preferred conformation in the cis form. The distribution of rotamers and the effect of an aromatic residue on the cis-trans equilibrium indicate that the preferred conformation is populated to approximately 62% for the Phe containing peptide, 67% for the Tyr containing peptide, and between 75 and 80% for the Trp containing peptide. The interaction is unaffected by the addition of 8M urea. These local interactions favor an aromatic residue immediately preceding a cis proline, but they cannot explain the relative propensities for Phe-Pro, Tyr-Pro, and Trp-Pro cis peptide bonds observed in the native state of proteins. In the model peptides the percentage of the cis proline conformer is 21% GYPG while it is 17% for GFPG. This difference is considerably smaller than the almost three to one preponderance observed for cis Tyr-Pro peptide bonds vs cis Phe-Pro peptide bonds in the protein database. © 1998 John Wiley & Sons, Inc. Biopoly 45: 381–394, 1998     

Lens Crystallins and Their Microbial Homologs: Structure,Stability, and Function

Referee: Franz Schmid, Biochemicshes Laboratorium, Universitaet Bayeuth, D-95440 Bayeuth, Germany

abg-Crystallins are the major protein components in the vertebrate eye lens — a as a molecular chaperone and b and g as structural proteins. Surprisingly, the latter two share some structural characteristics with a number of microbial stress proteins. The common denominator is not only the Greek key topology of their polypeptide chains but also their high intrinsic stability, which, in certain microbial crystallin homologs, is further enhanced by high-affinity Ca²⁺-binding. Recent studies of natural and mutant vertebrate bg-crystallins as well as spherulin 3a from Physarum polycephalum and Protein S from Myxococcus xanthus allowed the correlation of structure and stability of crystallins to be elucidated in some detail. From the thermo-dynamic point of view, stability increments come from (1) local interactions involved in the close packing of the cooperative units, (2) the all-b secondary structure of the Greek-key motif, (3) intramolecular interactions between domains, (4) intermolecular domain interactions, including 3D domain swapping and (v) excluded volume effects due to “molecular crowding” at the high cellular protein concentrations. Apart from these contributions to the Gibbs free energy of stability, significant kinetic stabilization originates from the high activation energy barrier determining the rate of unfolding from the native to the unfolded state. From the functional point of view, the high stability is responsible for the long-term transparency of the eye lens, on the one hand, and the stress resistance of the microorganisms in their dormant state on the other. Local structural perturbations due to chemical modification, wrong protein interactions, or other irreversible processes may lead to protein aggregation. A leading cataract hypothesis is that only after a-crystallin, a member of the small heat-shock protein family, is titrated out does pathological opacity occur. Understanding the structural basis of protein stability in the healthy eye lens is the route to solve the enormous medical and economical problem of cataract.     

Structural basis of hierarchical multiple substates of a protein. I: Introduction

A computer experiment of protein dynamics is carried out, which consists of two steps: (1) A Monte Carlo simulation of thermal fluctuations in the native state of a globular protein, bovine pancreatic trypsin inhibitor; and (2) a simulation of the quick freezing of fluctuating conformations into energy minima by minimization of the energy of a number of conformations sampled in the Monte Carlo simulation. From the analysis of results of the computer experiment is obtained the following picture of protein dynamics: multiple energy minima exist in the native state, and they are distributed in clusters in the conformational space. The dynamics has a hierarchical structure which has at least two levels. In the first level, dynamics is restricted within one of the clusters of minima. In the second, transitions occur among the clusters. Local parts of a protein molecule, side chains and local main chain segments, can take multiple locally stable conformations in the native state. Many minima result from combinations of these multiple local conformations. The hierarchical structure in the dynamics comes from interactions among the local parts. Protein molecules have two types of flexibility, each associated with elastic and plastic deformations, respectively.     

Kinetic consequences of native state optimization of surface-exposed electrostatic interactions in the Fyn SH3 domain

Optimization of surface exposed charge-charge interactions in the native state has emerged as an effective means to enhance protein stability; but the effect of electrostatic interactions on the kinetics of protein folding is not well understood. To investigate the kinetic consequences of surface charge optimization, we characterized the folding kinetics of a Fyn SH3 domain variant containing five amino acid substitutions that was computationally designed to optimize surface charge-charge interactions. Our results demonstrate that this optimized Fyn SH3 domain is stabilized primarily through an eight-fold acceleration in the folding rate. Analyses of the constituent single amino acid substitutions indicate that the effects of optimization of charge-charge interactions on folding rate are additive. This is in contrast to the trend seen in folded state stability, and suggests that electrostatic interactions are less specific in the transition state compared to the folded state. Simulations of the transition state using a coarse-grained chain model show that native electrostatic contacts are weakly formed, thereby making the transition state conducive to nonspecific, or even nonnative, electrostatic interactions. Because folding from the unfolded state to the folding transition state for small proteins is accompanied by an increase in charge density, nonspecific electrostatic interactions, that is, generic charge density effects can have a significant contribution to the kinetics of protein folding. Thus, the interpretation of the effects of amino acid substitutions at surface charged positions may be complicated and consideration of only native-state interactions may fail to provide an adequate picture.     

The role of atomic level steric effects and attractive forces in protein folding

Protein folding into tertiary structures is controlled by an interplay of attractive contact interactions and steric effects. We investigate the balance between these contributions using structure‐based models using an all‐atom representation of the structure combined with a coarse‐grained contact potential. Tertiary contact interactions between atoms are collected into a single broad attractive well between the C^β atoms between each residue pair in a native contact. Through the width of these contact potentials we control their tolerance for deviations from the ideal structure and the spatial range of attractive interactions. In the compact native state dominant packing constraints limit the effects of a coarse‐grained contact potential. During folding, however, the broad attractive potentials allow an early collapse that starts before the native local structure is completely adopted. As a consequence the folding transition is broadened and the free energy barrier is decreased. Eventually two‐state folding behavior is lost completely for systems with very broad attractive potentials. The stabilization of native‐like residue interactions in non‐perfect geometries early in the folding process frequently leads to structural traps. Global mirror images are a notable example. These traps are penalized by the details of the repulsive interactions only after further collapse. Successful folding to the native state requires simultaneous guidance from both attractive and repulsive interactions. Proteins 2012. © 2011 Wiley Periodicals, Inc.     

How strong are side chain interactions in the folding intermediate?

Influence of 12 nonpolar amino acids residues from the hydrophobic core of apomyoglobin on stability of its native state and folding intermediate was studied. Six of the selected residues are from the A, G and H helices; these are conserved in structure of the globin family, although nonfunctional, that is, not involved in heme binding. The rest are nonconserved hydrophobic residues that belong to the B, C, D, and E helices. Each residue was substituted by alanine, and equilibrium pH‐induced transitions in apomyoglobin and its mutants were studied by circular dichroism and fluorescent spectroscopy. The obtained results allowed estimating changes in their free energy during formation of the intermediate state. It was first shown that the strength of side chain interactions in the apomyoglobin intermediate state amounts to 15–50% of that in its native state for conserved residues, and practically to 0% for nonconserved residues. These results allow a better understanding of interactions occurring in the intermediate state and shed light on involvement of certain residues in protein folding at different stages.     

Local interactions favor the native 8-residue disulfide loop in the oxidation of a fragment corresponding to the sequence Ser-50-Met-79 derived from bovine pancreatic ribonuclease A

A 30-residue peptide was obtained from ribonuclease A by chemical cleavage with cyanogen bromide, subsequent sulfitolysis with concomitant S-sulfonation, and finally enzymatic cleavage withStaphylococcus aureus protease. The peptide was converted to the free thiol form by reductive cleavage of the S-sulfo-protecting groups withd,l-dithiothreitol. This peptide consisted of residues 50–79 of the native sequence of ribonuclease A, with the exception that methionine-79 had been converted to homoserine. Included in this sequence are residues cysteine-65 and cysteine-72, which form a disulfide bond in the native enzyme, as well as cysteine-58. This molecule may form one of three possible intramolecular disulfide bonds upon thiol oxidation, viz. one loop of 15 and 2 of 8 residues each. These isomeric peptides were prepared by oxidation with cystamine, 2-aminoethanethiolation of residual thiols, and fractionation by reverse-phase high-performance liquid chromatography. Disulfide pairings were established by mapping the tryptic fragments and confirming their composition by amino acid analysis. After protracted incubation under oxidizing conditions at 25.0°C andp H 8.0, the 26-member ring incorporating the native disulfide bond between residues 65 and 72 is the dominant product. Assuming that equilibrium is established, we infer that local interactions in the sequence of ribonuclease A significantly stabilize the native 8-residue disulfide loop with respect to the non-native 8-residue loop (G°=–1.1±0.1 kcal mole^–1). The implications of this observation for the oxidative folding of the intact protein are discussed.     

Effects of ligand binding on the stability of aldo–keto reductases: Implications for stabilizer or destabilizer chaperones

Ligands such as enzyme inhibitors stabilize the native conformation of a protein upon binding to the native state, but some compounds destabilize the native conformation upon binding to the non‐native state. The former ligands are termed “stabilizer chaperones” and the latter ones “destabilizer chaperones.” Because the stabilization effects are essential for the medical chaperone (MC) hypothesis, here we have formulated a thermodynamic system consisting of a ligand and a protein in its native‐ and non‐native state. Using the differential scanning fluorimetry and the circular dichroism varying the urea concentration and temperature, we found that when the coenzyme NADP⁺ was absent, inhibitors such as isolithocholic acid stabilized the aldo–keto reductase AKR1A1 upon binding, which showed actually the three‐state folding, but destabilized AKR1B10. In contrast, in the presence of NADP⁺, they destabilized AKR1A1 and stabilized AKR1B10. To explain these phenomena, we decomposed the free energy of stabilization (ΔΔG) into its enthalpy (ΔΔH) and entropy (ΔΔS) components. Then we found that in a relatively unstable protein showing the three‐state folding, native conformation was stabilized by the negative ΔΔH in association with the negative ΔΔS, suggesting that the stabilizer chaperon decreases the conformational fluctuation of the target protein or increase its hydration. However, in other cases, ΔΔG was essentially determined by the delicate balance between ΔΔH and ΔΔS. The proposed thermodynamic formalism is applicable to the system including multiple ligands with allosteric interactions. These findings would promote the development of screening strategies for MCs to regulate the target conformations.     

Folding of a beta-sheet protein monitored by real-time NMR spectroscopy

At low ionic strength, apoplastocyanin forms an unfolded state under non-denaturing conditions. The refolding of this state is sufficiently slow to allow real-time NMR experiments to be performed. Folding of apoplastocyanin, initiated by the addition of salt and followed by real-time 2D 1H-15N heteronuclear single quantum coherence (HSQC) spectroscopy, is highly cooperative. A concomitant increase in the intensity of both sequential and long-range nuclear Overhauser effects (NOEs) between backbone amide protons in successive acquisitions of 1H-15N HSQC-NOESY-HSQC spectra provides the first direct observation of the development of structure-specific NOEs as a protein folds. Our results show that the local and long-range interactions in the native apoplastocyanin are formed simultaneously, consistent with highly cooperative formation of the native structure.     

Monte Carlo studies on equilibrium globular protein folding. I. Homopolymeric lattice models of beta-barrel proteins

Dynamic Monte Carlo studies have been performed on various diamond lattice models of β-proteins. Unlike previous work, no bias toward the native state is introduced; instead, the protein is allowed to freely hunt through all of phase space to find the equilibrium conformation. Thus, these systems may aid in the elucidation of the rules governing protein folding from a given primary sequence; in particular, the interplay of short- vs long-range interaction can be explored. Three distinct models (A? C) were examined. In model A, in addition to the preference for trans (t) over gauche states (g⁺ and g^?) (thereby perhaps favoring β-sheet formation), attractive interactions are allowed between all nonbonded, nearest neighbor pairs of segments. If the molecules possess a relatively large fraction of t states in the denatured form, on cooling spontaneous collapse to a well-defined β-barrel is observed. Unfortunately, in model A the denatured state exhibits too much secondary structure to correctly model the globular protein collapse transition. Thus in models B and C, the local stiffness is reduced. In model B, in the absence of long-range interactions, t and g states are equally weighted, and cooperativity is introduced by favoring formation of adjacent pairs of nonbonded (but not necessarily parallel) t states. While the denatured state of these systems behaves like a random coil, their native globular structure is poorly defined. Model C retains the cooperativity of model B but allows for a slight preference of t over g states in the short-range interactions. Here, the denatured state is indistinguishable from a random coil, and the globular state is a well-defined β-barrel. Over a range of chain lengths, the collapse is well represented by an all-or-none model. Hence, model C possesses the essential qualitative features observed in real globular proteins. These studies strongly suggest that the uniqueness of the globular conformation requires some residual secondary structure to be present in the denatured state.     

Reduced C(beta) statistical potentials can outperform all-atom potentials in decoy identification

We developed a series of statistical potentials to recognize the native protein from decoys, particularly when using only a reduced representation in which each side chain is treated as a single C(beta) atom. Beginning with a highly successful all-atom statistical potential, the Discrete Optimized Protein Energy function (DOPE), we considered the implications of including additional information in the all-atom statistical potential and subsequently reducing to the C(beta) representation. One of the potentials includes interaction energies conditional on backbone geometries. A second potential separates sequence local from sequence nonlocal interactions and introduces a novel reference state for the sequence local interactions. The resultant potentials perform better than the original DOPE statistical potential in decoy identification. Moreover, even upon passing to a reduced C(beta) representation, these statistical potentials outscore the original (all-atom) DOPE potential in identifying native states for sets of decoys. Interestingly, the backbone-dependent statistical potential is shown to retain nearly all of the information content of the all-atom representation in the C(beta) representation. In addition, these new statistical potentials are combined with existing potentials to model hydrogen bonding, torsion energies, and solvation energies to produce even better performing potentials. The ability of the C(beta) statistical potentials to accurately represent protein interactions bodes well for computational efficiency in protein folding calculations using reduced backbone representations, while the extensions to DOPE illustrate general principles for improving knowledge-based potentials.