The first of the two reviews that follow is based on Professor Dill's Hans Neurath Award lecture given at the Twelfth Symposium of the Protein Society in July, 1998. Profs. Dill and Nussinov have collaborated to present a comprehensive picture of our current understanding of the protein folding process from the initial stages to the point at which the protein binds its ligands.

This paper surveys the emerging role of statistical mechanics and polymer theory in protein folding. In the polymer perspective, the folding code is more a solvation code than a code of local φψ propensities. The polymer perspective resolves two classic puzzles: (1) the Blind Watchmaker's Paradox that biological proteins could not have originated from random sequences, and (2) Levinthal's Paradox that the folded state of a protein cannot be found by random search. Both paradoxes are traditionally framed in terms of random unguided searches through vast spaces, and vastness is equated with impossibility. But both processes are partly guided. The searches are more akin to balls rolling down funnels than balls rolling aimlessly on flat surfaces. In both cases, the vastness of the search is largely irrelevant to the search time and success. These ideas are captured by energy and fitness landscapes. Energy landscapes give a language for bridging between microscopics and macroscopics, for relating folding kinetics to equilibrium fluctuations, and for developing new and faster computational search strategies.

Folding funnels have been the focus of considerable attention during the last few years. These have mostly been discussed in the general context of the theory of protein folding. Here we extend the utility of the concept of folding funnels, relating them to biological mechanisms and function. In particular, here we describe the shape of the funnels in light of protein synthesis and folding; flexibility, conformational diversity, and binding mechanisms; and the associated binding funnels, illustrating the multiple routes and the range of complexed conformers. Specifically, the walls of the folding funnels, their crevices, and bumps are related to the complexity of protein folding, and hence to sequential vs. nonsequential folding. Whereas the former is more frequently observed in eukaryotic proteins, where the rate of protein synthesis is slower, the latter is more frequent in prokaryotes, with faster translation rates. The bottoms of the funnels reflect the extent of the flexibility of the proteins. Rugged floors imply a range of conformational isomers, which may be close on the energy landscape. Rather than undergoing an induced fit binding mechanism, the conformational ensembles around the rugged bottoms argue that the conformers, which are most complementary to the ligand, will bind to it with the equilibrium shifting in their favor. Furthermore, depending on the extent of the ruggedness, or of the smoothness with only a few minima, we may infer nonspecific, broad range vs. specific binding. In particular, folding and binding are similar processes, with similar underlying principles. Hence, the shape of the folding funnel of the monomer enables making reasonable guesses regarding the shape of the corresponding binding funnel. Proteins having a broad range of binding, such as proteolytic enzymes or relatively nonspecific endonucleases, may be expected to have not only rugged floors in their folding funnels, but their binding funnels will also behave similarly, with a range of complexed conformations. Hence, knowledge of the shape of the folding funnels is biologically very useful. The converse also holds: If kinetic and thermodynamic data are available, hints regarding the role of the protein and its binding selectivity may be obtained. Thus, the utility of the concept of the funnel carries over to the origin of the protein and to its function.

Elicitins, produced by most of the phytopathogenic fungi of the genus Phytophthora, provoke in tobacco both remote leaf necrosis and the induction of a resistance against subsequent attack by various microorganisms. Despite the recent description of the three-dimensional crystal structure of cryptogein (CRY), the molecular basis of the interactions between Phytophthora and plants largely remains unknown. The X-ray crystal structure, refined at 2.1 Å, of a ligand complexed, mutated CRY, K13H, is reported. Analysis of this structure reveals that CRY is able to encapsulate a ligand that induces only a minor conformational change in the protein structure. The ligand has been identified as an ergosterol by gas chromatographic analysis coupled with mass spectrometry analysis. This result is consistent with biochemical data that have shown that elicitins are a distinct class of Sterol Carrier Proteins (SCP). Data presented here provide the first structural description of the pertinent features of the elicitin sterol interaction and permit a reassessment of the importance of both the key residue 13 and the mobility of the omega loop for the accessibility of the sterol to the cavity. The biological implications thereof are discussed. This paper reports the first structure of a SCP/sterol complex.

Two equilibrium intermediates have previously been observed in the urea denaturation of the α subunit of tryptophan synthase (αTS) from Escherichia coli, an eight-stranded β/α barrel protein. In the current study, a series of amino-terminal fragments were characterized to probe the elementary folding units that may be in part responsible for this complex behavior. Stop-codon mutagenesis was used to produce eight fragments ranging in size from 105–214 residues and containing incremental elements of secondary structure. Equilibrium studies by circular dichroism indicate that all of these fragments are capable of adopting secondary structure. All except for the shortest fragment fold cooperatively. The addition of the fourth, sixth, and eighth β-strands leads to distinct increases in structure, cooperativity, and/or stability, suggesting that folding involves the modular assembly of βαβ supersecondary structural elements. One-dimensional NMR titrations at high concentrations of urea, probing the environment around His92, were also performed to test for the presence of residual structure in the fragments. All fragments that contained the first four βα units of structure exhibited a cooperative unfolding transition at high concentrations of urea with significant but reduced stability relative to the full-length protein. These results suggest that the residual structure in αTS requires the participation of hydrophobic residues in multiple β-strands that span the entire sequence.

We carry out a systematic analysis of the correlation between similarity of protein three-dimensional structures and their evolutionary relationships. The structural similarity is quantitatively identified by an all-against-all comparison of the spatial arrangement of secondary structural elements in nonredundant 967 representative proteins, and the evolutionary relationship is judged according to the definition of superfamily in the SCOP database. We find the following symmetry rule: a protein pair that has similar folds but belong to different superfamilies has (with a very rare exception) certain internal symmetry in its common similar folds. Possible reasons behind the symmetry rule are discussed.

The first crystal structure of an inorganic pyrophosphatase (S-PPase) from an archaebacterium, the thermophile Sulfolobus acidocaldarius, has been solved by molecular replacement and refined to an R-factor of 19.7% at 2.7 Å. S-PPase is a D3 homohexameric protein with one Mg2+ per active site in a position similar to, but not identical with, the first activating metal in mesophilic pyrophosphatases (PPase). In mesophilic PPases, Asp65, Asp70, and Asp102 coordinate the Mg2+, while only Asp65 and Asp102 do in S-PPase, and the Mg2+ moves by 0.7 Å. S-PPase may therefore be deactivated at low temperature by mispositioning a key metal ion.

The monomer S-PPase structure is very similar to that of Thermus thermophilus (T-PPase) and Escherichia coli (E-PPase), root-mean-square deviations around 1 Å/Cα. But the hexamer structures of S- and T-PPase are more tightly packed and more similar to each other than they are to that of E-PPase, as shown by the increase in surface area buried upon oligomerization. In T-PPase, Arg116 creates an interlocking ionic network to both twofold and threefold related monomers; S-PPase has hydrophilic interactions to threefold related monomers absent in both E- and T-PPase. In addition, the thermostable PPases have about 7% more hydrogen bonds per monomer than E-PPase, and, especially in S-PPase, additional ionic interactions anchor the C-terminus to the rest of the protein. Thermostability in PPases is thus due to subtle improvements in both monomer and oligomer interactions.

Cytochrome cH is the electron donor to the oxidase in methylotrophic bacteria. Its amino acid sequence suggests that it is a typical Class I cytochrome c, but some features of the sequence indicated that its structure might be of special interest. The structure of oxidized cytochrome cH has been solved to 2.0 Å resolution by X-ray diffraction. It has the classical tertiary structure of the Class 1 cytochromes c but bears a closer gross resemblance to mitochondrial cytochrome c than to the bacterial cytochrome c2. The left-hand side of the haem cleft is unique; in particular, it is highly hydrophobic, the usual water is absent, and the “conserved” Tyr67 is replaced by tryptophan. A number of features of the structure demonstrate that the usual hydrogen bonding network involving water in the haem channel is not essential and that other mechanisms may exist for modulation of redox potentials in this cytochrome.

Principles of protein thermostability have been studied by comparing structures of thermostable proteins with mesophilic counterparts that have a high degree of sequence identity. Two tetrameric NADP(H)-dependent alcohol dehydrogenases, one from Clostridium beijerinckii (CBADH) and the other from Thermoanaerobacter brockii (TBADH), having exceptionally high (75%) sequence identity, differ by 30° in their melting temperatures. The crystal structures of CBADH and TBADH in their holo-enzyme form have been determined at a resolution of 2.05 and 2.5 Å, respectively. Comparison of these two very similar structures (RMS difference in Cα = 0.8 Å) revealed several features that can account for the higher thermal stability of TBADH. These include additional ion pairs, “charged-neutral” hydrogen bonds, and prolines as well as improved stability of α-helices and tighter molecular packing. However, a deeper structural insight, based on the location of stabilizing elements, suggests that enhanced thermal stability of TBADH is due mainly to the strategic placement of structural determinants at positions that strengthen the interface between its subunits. This is also supported by mutational analysis of structural elements at critical locations. Thus, it is the reinforcement of the quaternary structure that is most likely to be a primary factor in preserving enzymatic activity of this oligomeric bacterial ADH at elevated temperatures.

The X-ray structure of the complex formed between yeast 5-aminolaevulinic acid dehydratase (ALAD) and the inhibitor laevulinic acid has been determined at 2.15 Å resolution. The inhibitor binds by forming a Schiff base link with one of the two invariant lysines at the catalytic center: Lys263. It is known that this lysine forms a Schiff base link with substrate bound at the enzyme's so-called P-site. The carboxyl group of laevulinic acid makes hydrogen bonds with the side-chain–OH groups of Tyr329 and Ser290, as well as with the main-chain >NH group of Ser290. The aliphatic moiety of the inhibitor makes hydrophobic interactions with surrounding aromatic residues in the protein including Phe219, which resides in the flap covering the active site. Our analysis strongly suggests that the same interactions will be made by P-side substrate and also indicates that the substrate that binds at the enzyme's A-site will interact with the enzyme's zinc ion bound by three cysteines (133, 135, and 143). Inhibitor binding caused a substantial ordering of the active site flap (residues 217–235), which was largely invisible in the native electron density map and indicates that this highly conserved yet flexible region has a specific role in substrate binding during catalysis.

Extended retro (reversed) peptide sequences have not previously been accommodated within functional proteins. Here, we show that the entire transmembrane portion of the β-barrel of the pore-forming protein α-hemolysin can be formed by retrosequences comprising a total of 175 amino acid residues, 25 contributed by the central sequence of each subunit of the heptameric pore. The properties of wild-type and retro heptamers in planar bilayers are similar. The single-channel conductance of the retro pore is 15% less than that of the wild-type heptamer and its current-voltage relationship denotes close to ohmic behavior, while the wild-type pore is weakly rectifying. Both wild-type and retro pores are very weakly anion selective. These results and the examination of molecular models suggest that β-barrels may be especially accepting of retro sequences compared to other protein folds. Indeed, the ability to form a retro domain could be diagnostic of a β-barrel, explaining, for example, the activity of the retro forms of many membrane-permeabilizing peptides. By contrast with the wild-type subunits, monomeric retro subunits undergo premature assembly in the absence of membranes, most likely because the altered central sequence fails to interact with the remainder of the subunit, thereby initiating assembly. Despite this difficulty, a technique was devised for obtaining heteromeric pores containing both wild-type and retro subunits. Most probably as a consequence of unfavorable interstrand side-chain interactions, the heteromeric pores are less stable than either the wild-type or retro homoheptamers, as judged by the presence of subconductance states in single-channel recordings. Knowledge about the extraordinary plasticity of the transmembrane β-barrel of α-hemolysin will be very useful in the de novo design of functional membrane proteins based on the β-barrel motif.

A sequence motif that is Similar to Ubiquitin (SUb) has been identified in the Saccharomyces cerevisiae ubiquitin-specific protease Ubp6. SUb is conserved in all known Ubp6 homologues from a spectrum of eukaryotic species and is also present in a group of hypothetical proteins of unknown function (Unk1-3) present in sequence databases. An N-terminal deletion mutant of Ubp6 that lacks SUb is still capable of cleaving α-linked ubiquitin fusions, suggesting that SUb forms a separate domain to the catalytic core of Ubp6 and demonstrating that it is not required for in vitro cleavage activity. A homology model of the 78 N-terminal amino acids of human Ubp6, based on the known fold of ubiquitin, is presented. In human Ubp6, SUb shares only 20% sequence identity with ubiquitin. Even weaker similarity occurs between S. cerevisiae SUb and ubiquitin. The homology model supports a ubiquitin-like fold for SUb and suggests that two conserved Lys residues, corresponding to Lys48 and Lys63 of ubiquitin, are functionally important.

The Escherichia coli tryptophan repressor protein (TR) represses the transcription of several genes in response to the concentration of tryptophan in the environment. In the co-crystal structure of TR bound to a DNA fragment containing its target very few direct contacts between TR and the DNA were observed. In contrast, a number of solvent mediated contacts were apparent. NMR solution structures, however, did not resolve any solvent mediated bonds at the complex interface. To probe for the role of water in TR operator recognition, the effect of osmolytes on the interactions between TR and a target oligonucleotide bearing the operator site was examined. In the absence of specific solvent mediated hydrogen bonding interactions between the protein and the DNA, increasing osmolyte concentration is expected to strongly stabilize the TR operator interaction due to the large amount of macromolecular surface area buried upon complexation. The results of our studies indicate that xylose did not alter the binding affinity significantly, while glycerol and PEG had a small stabilizing effect. A study of binding as a function of betaine concentration revealed that this osmolyte at low concentration results in a stabilization of the 1:1 TR/operator complex, but at higher concentrations leads to a switching between binding modes to favor tandem binding. Analysis of the effects of betaine on the 1:1 complex suggest that this osmolyte has about 78% of the expected effect. If one accepts the analysis in terms of the number of water molecules excluded upon complexation, these results suggest that about 75 water molecules remain at the interface of the 1:1 dimer/DNA complex. This value is consistent with the number of water molecules found at the interface in the crystallographically determined structure and supports the notion that interfacial waters play an important thermodynamic role in the specific complexation of one TR dimer with its target DNA. However, the complexity of the effects of betaine and the small or negligible effects of the other osmolytes could also arise from osmolyte induced competition between antagonistic coupled reactions.

The refolding of barstar, the intracellular inhibitor of barnase, is dominated by the slow formation of a cis peptidyl prolyl bond in the native protein. The triple mutant C40/82A P27A in which two cysteine residues and one trans proline were replaced by alanine was used as model system to investigate the kinetics and structural consequences of the trans/cis interconversion of Pro48. One- and two-dimensional real-time NMR spectroscopy was used to follow the trans/cis interconversion after folding was initiated by rapid dilution of the urea denatured protein. Series of 1H, 15N HSQC spectra acquired with and without the addition of peptidyl prolyl isomerase unambiguously revealed the accumulation of a transient trans-Pro48 intermediate within the dead time of the experiment. Subtle chemical shift differences between the native state and the intermediate spectra indicate that the intermediate is predominantly native-like with a local rearrangement in the Pro48 loop and in the β-sheet region including residues Tyr47, Ala82, Thr85, and Val50.

BBA1 is a designed protein that has only 23 residues. It is the smallest protein without disulfide bridges that has a well-defined tertiary structure in solution. We have performed unfolding molecular dynamics simulations on BBA1 and some of its mutants at 300, 330, 360, and 400 K to study their kinetic stability as well as the unfolding mechanism of BBA1. It was shown that the unfolding simulations can provide insights into the forces that stabilize the protein. Packing, hydrophobic interactions, and a salt bridge between Asp12 and Lys16 were found to be important to the protein's stability. The unfolding of BBA1 goes through two major steps: (1) disruption of the hydrophobic core and (2) unfolding of the helix. The β-hairpin remains stable in the unfolding because of the high stability of the type II′ turn connecting the two β-strands.

The X-ray structure of the Escherichia coli aspartate transcarbamoylase with the bisubstrate analog phosphonacetyl-l-aspartate (PALA) bound shows that PALA interacts with Lys84 from an adjacent catalytic chain. To probe the function of Lys84, site-specific mutagenesis was used to convert Lys84 to alanine, threonine, and asparagine. The K84N and K84T enzymes exhibited 0.08 and 0.29% of the activity of the wild-type enzyme, respectively. However, the K84A enzyme retained 12% of the activity of the wild-type enzyme. For each of these enzymes, the affinity for aspartate was reduced 5- to 10-fold, and the affinity for carbamoyl phosphate was reduced 10- to 30-fold. The enzymes K84N and K84T exhibited no appreciable cooperativity, whereas the K84A enzyme exhibited a Hill coefficient of 1.8. The residual cooperativity and enhanced activity of the K84A enzyme suggest that in this enzyme another mechanism functions to restore catalytic activity. Modeling studies as well as molecular dynamics simulations suggest that in the case of only the K84A enzyme, the lysine residue at position 83 can reorient into the active site and complement for the loss of Lys84. This hypothesis was tested by the creation and analysis of the K83A enzyme and a double mutant enzyme (DM) that has both Lys83 and Lys84 replaced by alanine. The DM enzyme has no cooperativity and exhibited 0.18% of wild-type activity, while the K83A enzyme exhibited 61% of wild-type activity. These data suggest that Lys84 is not only catalytically important, but is also essential for binding both substrates and creation of the high-activity, high-affinity active site. Since low-angle X-ray scattering demonstrated that the mutant enzymes can be converted to the R-structural state, the loss of cooperativity must be related to the inability of these mutant enzymes to form the high-activity, high-affinity active site characteristic of the R-functional state of the enzyme.

Urea-induced denaturations of RNase T1 and reduced and carboxyamidated RNase T1 (RTCAM) as a function of temperature were analyzed using the linear extrapolation method, and denaturation m values, ΔCp, ΔH, ΔS, and ΔG quantities were determined. Because both ΔCp and m values are believed to reflect the protein surface area newly exposed on denaturation, the prediction is that the ratio of m values for RNase T1 and RTCAM should equal the ΔCp ratio for the two proteins. This is not the case, for it is found that the m value of RTCAM is 1.5 times that of RNase T1, while the denaturation ΔCp's for the two proteins are identical. The paradox of why the two parameters, m and ΔCp, are not equivalent in their behavior is of importance in the interpretations of their respective molecular-level meanings. It is found that the measured denaturation ΔCp's are consistent with ΔCp's calculated on the basis of empirical relationships between the change in surface area on denaturation (ΔASA), and that the measured m value of RNase T1 agrees with m calculated from empirical data relating m to ΔASA. However, the measured m of RTCAM is so much out of line with its calculated m as to call into question the validity of always equating m with surface area newly exposed on denaturation.

The N-terminal 17 residues of ubiquitin have been shown by 1H NMR to fold autonomously into a β-hairpin structure in aqueous solution. This structure has a specific, native-like register, though side-chain contacts differ in detail from those observed in the intact protein. An autonomously folding hairpin has previously been identified in the case of streptococcal protein G, which is structurally homologous with ubiquitin, but remarkably, the two are not in topologically equivalent positions in the fold. This suggests that the organization of folding may be quite different for proteins sharing similar tertiary structures. Two smaller peptides have also been studied, corresponding to the isolated arms of the N-terminal hairpin of ubiquitin, and significant differences from simple random coil predictions observed in the spectra of these subfragments, suggestive of significant limitation of the backbone conformational space sampled, presumably as a consequence of the strongly β-structure favoring composition of the sequences. This illustrates the ability of local sequence elements to express a propensity for β-structure even in the absence of actual sheet formation. Attempts were made to estimate the population of the folded state of the hairpin, in terms of a simple two-state folding model. Using published “random coil” values to model the unfolded state, and values derived from native ubiquitin for the putative unique, folded state, it was found that the apparent population varied widely for different residues and with different NMR parameters. Use of the spectra of the subfragment peptides to provide a more realistic model of the unfolded state led to better agreement in the estimates that could be obtained from chemical shift and coupling constant measurements, while making it clear that some other approaches to population estimation could not give meaningful results, because of the tendency to populate the β-region of conformational space even in the absence of the hairpin structure.

Hepatitis C virus (HCV) nonstructural protein 3 (NS3) has been shown to possess protease and helicase activities and has also been demonstrated to spontaneously associate with nonstructural protein NS4A (NS4A) to form a stable complex. Previous attempts to produce the NS3/NS4A complex in recombinant baculovirus resulted in a protein complex that aggregated and precipitated in the absence of nonionic detergent and high salt. A single-chain form of the NS3/NS4A complex (His-NS4A21–32–GSGS-NS33–631) was constructed in which the NS4A core peptide is fused to the N-terminus of the NS3 protease domain as previously described (Taremi et al., 1998). This protein contains a histidine tagged NS4A peptide (a.a. 21–32) fused to the full-length NS3 (a.a. 3–631) through a flexible tetra amino acid linker. The recombinant protein was expressed to high levels in Escherichia coli, purified to homogeneity, and examined for NTPase, nucleic acid unwinding, and proteolytic activities. The single-chain recombinant NS3-NS4A protein possesses physiological properties equivalent to those of the NS3/NS4A complex except that this novel construct is stable, soluble and sixfold to sevenfold more active in unwinding duplex RNA. Comparison of the helicase activity of the single-chain recombinant NS3-NS4A with that of the full-length NS3 (without NS4A) and that of the helicase domain alone suggested that the presence of the protease domain and at least the NS4A core peptide are required for optimal unwinding activity.

The structure of a bifunctional 5,10-methylene-tetrahydrofolate dehydrogenase/cyclohydrolase from Escherichia coli has been determined at 2.5 Å resolution in the absence of bound substrates and compared to the NADP-bound structure of the homologous enzyme domains from a trifunctional human synthetase enzyme. Superposition of these structures allows the identification of a highly conserved cluster of basic residues that are appropriately positioned to serve as a binding site for the poly-γ-glutamyl tail of the tetrahydrofolate substrate. Modeling studies and molecular dynamic simulations of bound methylene-tetrahydrofolate and NADP shows that this binding site would allow interaction of the nicotinamide and pterin rings in the dehydrogenase active site. Comparison of these enzymes also indicates differences between their active sites that might allow the development of inhibitors specific to the bacterial target.