The role of helix formation in the folding of a fully alpha-helical coiled coil

Preformed secondary structure drives the association reaction of GCN4-p1, a model coiled-coil system

The structure of the transition state for the rate-limiting step in the folding and association of the homodimeric coiled-coil peptide GCN4-p1, was probed by mutational analysis. A series of quadruple amino acid replacements that spanned the helix propensity scale were made at the four external f positions in the heptad repeat. Equilibrium and kinetic circular dichroism studies demonstrate that both the stability and the unfolding and refolding rate constants vary with helix propensity but also reflect interactions of the altered side-chains with their local environments. Pairwise replacements and fragment studies show that the two C-terminal heptads are the likely source of the nucleating helices. Helix-helix recognition between preformed elements of secondary structure plays an important role in this fundamental folding reaction.

Minimal aggregate size and minimal fusion unit for the first fusion pore of influenza hemagglutinin-mediated membrane fusion

The data of Melikyan et al. (J. Gen. Physiol. 106:783, 1995) for the time required for the first measurable step of fusion, the formation of the first flickering conductivity pore between influenza hemagglutinin (HA) expressing cells and planar bilayers, has been analyzed using a new mass action kinetic model. The analysis incorporates a rigorous distinction between the minimum number of HA trimers aggregated at the nascent fusion site (which is denoted the minimal aggregate size) and the number of those trimers that must to undergo a slow essential conformational change before the first fusion pore could form (which is denoted the minimal fusion unit). At least eight (and likely more) HA trimers aggregated at the nascent fusion site. Remarkably, of these eight (or more) HAs, only two or three must undergo the essential conformational change slowly before the first fusion pore can form. Whether the conformational change of these first two or three HAs are sufficient for the first fusion pore to form or whether the remaining HAs within the aggregate must rapidly transform in a cooperative manner cannot be determined kinetically. Remarkably, the fitted halftime for the essential HA conformational change is roughly 10(4) s, which is two orders of magnitude slower than the observed halftime for fusion. This is because the HAs refold with distributed kinetics and because the conductance assay monitored the very first aggregate to succeed in forming a first fusion pore from an ensemble of hundreds or thousands (depending upon the cell line) of fusogenic HA aggregates within the area of apposition between the cell and the planar bilayer. Furthermore, the average rate constant for this essential conformational change was at least 10(7) times slower than expected for a simple coiled coil conformational change, suggesting that there is either a high free energy barrier to fusion and/or very many nonfusogenic conformations in the refolding landscape. Current models for HA-mediated fusion are examined in light of these new constraints on the early structure and evolution of the nascent fusion site. None completely comply with the data.

Rate and equilibrium constants for protein unfolding and refolding determined by hydrogen exchange-mass spectrometry

Studies of protein unfolding and refolding may help us understand the more general problem of protein folding. Recent studies from this laboratory demonstrated that the unfolding and refolding of a large protein, rabbit muscle aldolase (M(r) 157 kDa), can be studied by combining amide hydrogen exchange and mass spectrometry. Results of these studies indicated that aldolase has three unfolding domains which likely unfold sequentially. Urea was used to increase the populations of partially unfolded states which were labeled with deuterium following a brief exposure to D(2)O. Electrospray ionization mass spectra of both the intact protein and its peptic fragments had multiple envelopes of isotope peaks from which the populations of unfolded forms were determined. The present study extends the previous investigations to include different urea concentrations and kinetic modeling of data taken as the system approaches equilibrium. Analysis of these results gives rate and equilibrium constants describing the unfolding and refolding processes characteristic of aldolase destabilized in urea. The change in solvent-accessible surface, which has been used as a reaction coordinate for protein folding, is estimated from the dependence of the equilibrium and rate constants on the concentration of urea.

Restraint-driven formation of alpha-helical coiled coils in molecular dynamics simulations

The alpha-helical coiled coil motif is among the first characterized and widely found architecture of protein structures. We report here a fast and reliable approach of simulated annealing molecular dynamics (SA/MD) for predicting the three-dimensional structures of various alpha-helical coiled coils of heptad repeat. One key element of our simulation involves a geometric restraint requiring residues occupying the first and fourth positions of the heptad to orient to the angle of their respective statistical average derived from a survey of coiled-coil structures deposited in the Protein Data Bank. Another is the incorporation of subunit rotation and inversion operations for generating symmetrized protein assemblies during the dynamics simulations. The procedure is fully automated and can be applied to different oligomerization states of identical subunits, as well as both parallel and antiparallel arrangements. Despite simplicity, the formation of five coiled-coil prototype systems driven by the restraint-based SA/MD approach shows that the level of prediction accuracy achieved previously by more elaborate procedures can be retained. The present work thus provides validation of a simulation approach that can be employed to utilize a wide variety of knowledge-based geometric restraints for structural prediction of symmetrical or pseudo-symmetrical protein systems.

Transition state heterogeneity in GCN4 coiled coil folding studied by using multisite mutations and crosslinking

We have investigated the folding behavior of dimeric and covalently crosslinked versions of the 33-residue alpha-helical GCN4-p1 coiled coil derived from the leucine zipper region of the transcriptional activator GCN4. The effects of multisite substitutions indicate that folding occurs along multiple routes with nucleation sites located throughout the protein. The similarity in activation energies of the different routes together with an analysis of intrinsic helical propensities indicate that minimal helix is present before a productive collision of the two chains. However, approximately one-third to one-half of the total helical structure is formed in the postcollision transition state ensemble. For the crosslinked, monomeric version, folding occurs along a single robust pathway. Here, the region nearest the crosslink, with the least helical propensity, is structured in the transition state whereas the region farthest from the tether, with the most propensity, is completely unstructured. Hence, the existence of transition state heterogeneity and the selection of folding routes critically depend on chain topology.

De novo simulations of the folding thermodynamics of the GCN4 leucine zipper

Entropy Sampling Monte Carlo (ESMC) simulations were carried out to study the thermodynamics of the folding transition in the GCN4 leucine zipper (GCN4-lz) in the context of a reduced model. Using the calculated partition functions for the monomer and dimer, and taking into account the equilibrium between the monomer and dimer, the average helix content of the GCN4-lz was computed over a range of temperatures and chain concentrations. The predicted helix contents for the native and denatured states of GCN4-lz agree with the experimental values. Similar to experimental results, our helix content versus temperature curves show a small linear decline in helix content with an increase in temperature in the native region. This is followed by a sharp transition to the denatured state. van't Hoff analysis of the helix content versus temperature curves indicates that the folding transition can be described using a two-state model. This indicates that knowledge-based potentials can be used to describe the properties of the folded and unfolded states of proteins.

Reinterpretation of GCN4-p1 folding kinetics: partial helix formation precedes dimerization in coiled coil folding

The folding of coiled coil peptides has traditionally been interpreted in terms of native dimer and unfolded monomers. Calculations using AGADIR and experimental studies of fragments suggest that the monomers of the coiled coil peptide, GCN4-p1, contain significant residual helical structure. A simple model based on diffusion-collision theory predicts not only the measured folding rate within an order of magnitude, but also predicts remarkably well the effect of alanine to glyXcine mutations. We suggest that intrinsic helix stability is a major determinant of the folding rate of the GCN4 coiled coil.

Effects of varying the local propensity to form secondary structure on the stability and folding kinetics of a rapid folding mixed alpha/beta protein: characterization of a truncation mutant of the N-terminal domain of the ribosomal protein L9

The N-terminal domain of the ribosomal protein L9 forms a split betaalphabeta structure with a long C-terminal helix. The folding transitions of a 56 residue version of this protein have previously been characterized, here we report the results of a study of a truncation mutant corresponding to residues 1-51. The 51 residue protein adopts the same fold as the 56 residue protein as judged by CD and two-dimensional NMR, but it is less stable as judged by chemical and thermal denaturation experiments. Studies with synthetic peptides demonstrate that the C-terminal helix of the 51 residue version has very little propensity to fold in isolation in contrast to the C-terminal helix of the 56 residue variant. The folding rates of the two proteins, as measured by stopped-flow fluorescence, are essentially identical, indicating that formation of local structure in the C-terminal helix is not involved in the rate-limiting step of folding.

A specific transition state for S-peptide combining with folded S-protein and then refolding

We measured the folding and unfolding kinetics of mutants for a simple protein folding reaction to characterize the structure of the transition state. Fluorescently labeled S-peptide analogues combine with S-protein to form ribonuclease S analogues: initially, S-peptide is disordered whereas S-protein is folded. The fluorescent probe provides a convenient spectroscopic probe for the reaction. The association rate constant, kon, and the dissociation rate constant, koff, were both determined for two sets of mutants. The dissociation rate constant is measured by adding an excess of unlabeled S-peptide analogue to a labeled complex (RNaseS*). This strategy allows kon and koff to be measured under identical conditions so that microscopic reversibility applies and the transition state is the same for unfolding and refolding. The first set of mutants tests the role of the alpha-helix in the transition state. Solvent-exposed residues Ala-6 and Gln-11 in the alpha-helix of native RNaseS were replaced by the helix destabilizing residues glycine or proline. A plot of log kon vs. log Kd for this series of mutants is linear over a very wide range, with a slope of -0.3, indicating that almost all of the molecules fold via a transition state involving the helix. A second set of mutants tests the role of side chains in the transition state. Three side chains were investigated: Phe-8, His-12, and Met-13, which are known to be important for binding S-peptide to S-protein and which also contribute strongly to the stability of RNaseS*. Only the side chain of Phe-8 contributes significantly, however, to the stability of the transition state. The results provide a remarkably clear description of a folding transition state.

Long-time dynamics of Met-enkephalin: comparison of theory with Brownian dynamics simulations

A recent theory for the long time dynamics of flexible chain molecules is applied for the first time to a peptide of biological importance, the neurotransmitter met-enkephalin. The dynamics of met-enkephalin is considerably more complicated than that of the previously studied glycine oligomers; met-enkephalin contains the interesting motions of phenyl groups and of side chains relative to the backbone, motions that are present in general flexible peptides. The theory extends the generalized Rouse (GR) model used to study the dynamics of polymers by providing a systematic procedure for including the contributions from the memory function matrices neglected in the GR theory. The new method describes the dynamics by time correlation functions instead of individual trajectories. These correlation functions are analytically expressed in terms of a set of equilibrium averages and the eigenvalues and eigenfunctions of the diffusion operator. The predictions of the theory are compared with Brownian dynamics (BD) simulations, so that both theory and simulation use identical potential functions and solvent models. The theory thus contains no adjustable parameters. Inclusion of the memory function contributions profoundly affects the dynamics. The theory produces very good agreement with the BD simulations for the global motions of met-enkephalin. It also correctly predicts the long-time relaxation rate for local motions.

Conformational analysis of peptide fragments derived from the peripheral subunit-binding domain from the pyruvate dehydrogenase multienzyme complex of Bacillus stearothermophilus: evidence for nonrandom structure in the unfolded state

There is currently a great deal of interest in the early events in protein folding. Two issues that have generated particular interest are the nature of the unfolded state under native conditions and the role of local interactions in folding. Here, we report the results of a study of a set of peptides derived from a small two-helix protein, the peripheral subunit-binding domain of the pyruvate dehydrogenase multienzyme complex. Five peptides of overlapping sequence were prepared, including sequences corresponding to each of the helices and to the region connecting them. The peptides were characterized by CD and, where possible, nmr. A peptide corresponding to the second helix is between 12 and 17% helical at neutral pH. CD also indicates a lower percentage of helical structure in the peptide corresponding to the first alpha-helix, although the values of the alpha-proton chemical shifts suggest some preference for nonrandom structure. Peptides corresponding to the interhelical loop, which in the full domain contains two overlapping beta-turns and a 5-residue 3(10)-helix, are less structured. There is no significant change in the helicity of any of these peptides with pH. To test for fragment complementation, CD spectra of the two peptides derived from each helix and the long connecting peptide were compared to the spectra of each possible pair, as well as to a mixture containing all three. No increase in structure was observed. We complement our peptide studies by characterizing a point mutant, D34V, which disrupts a critical hydrogen bonding network. This mutant is unable to fold and provides a useful model of the denatured state. The mutant is between 9 and 16% helical as judged by CD. The modest amount of helical structure formed in some of the peptide fragments and in the point mutant suggests that the denatured state of the peripheral subunit binding domain is not completely unstructured. This may contribute to the very rapid folding observed for the intact protein.

The single helix in protein L is largely disrupted at the rate-limiting step in folding

To investigate the role of helix formation in the folding of protein L, a 62 residue alpha/beta protein, we studied the consequences of both single and multiple mutations in the helix on the kinetics of folding. A triple mutant with 11 additional carbon atoms in core residues in the amino-terminal portion of the helix folded substantially faster than wild type, suggesting that hydrophobic association with residues elsewhere in the protein occurs at the rate-limiting step in folding. However, helix-destabilizing mutations had little effect on the rate of folding; in particular, a triple glycine substitution on the solvent-exposed side of the helix increased the unfolding rate 56-fold while reducing the folding rate less than threefold. Thus, in contrast to the predictions of models of folding involving the coalescence of well-formed secondary structure elements, the single helix in protein L appears to be largely disrupted at the rate-limiting step in folding and unfolding.

An autonomous folding unit mediates the assembly of two-stranded coiled coils

Subunit oligomerization of many proteins is mediated by coiled-coil domains. Although the basic features contributing to the thermodynamic stability of coiled coils are well understood, the mechanistic details of their assembly have not yet been dissected. Here we report a 13-residue sequence pattern that occurs with limited sequence variations in many two-stranded coiled coils and that is absolutely required for the assembly of the Dictyostelium discoideum actin-bundling protein cortexillin I and the yeast transcriptional activator GCN4. The functional relationship between coiled-coil "trigger" sequences was manifested by replacing the intrinsic trigger motif of GCN4 with the related sequence from cortexillin I. We demonstrate that these trigger sequences represent autonomous helical folding units that, in contrast to arbitrarily chosen heptad repeats, can mediate coiled-coil formation. Aside from being of general interest for protein folding, trigger motifs should be of particular importance in the protein de novo design.

Supercoiled protein motifs: the collagen triple-helix and the alpha-helical coiled coil

The collagen triple-helix and the alpha-helical coiled coil represent the two basic supercoiled multistranded protein motifs. Originally they were characterized in fibrous proteins, but have been found more recently in a number of other proteins containing rod-shaped domains. Coiled-coil domains are responsible for the oligomerization of proteins, as well as other specific functions, while the triple-helix domains associate to form supramolecular structures and bind a variety of ligands. Both structures were originally solved by fiber diffraction, and recent crystallographic studies on small proteins and peptide models have confirmed the structure and provided molecular details. The differences in the molecular conformations of these two motifs and the interactions stabilizing these conformations are discussed. The molecular structures of both motifs constrain the amino acid sequence to recognizable patterns, requiring the (Gly-X-Y)n repeating sequence for the collagen triple-helix and a less stringent heptad repeat requirement (h-x-x-h-x-x-x)n for the coiled-coil domains, where h represents hydrophobic residues. The features and roles of these supercoiled domains in proteins are considered when they are found adjacent in the same protein.

Amide protection in an early folding intermediate of cytochrome c

BACKGROUND

For many proteins, compact states appear long before the rate-limiting step in the formation of the native structure. A key issue is whether the initial collapse of the chain is driven by random or more specific hydrophobic interactions.

RESULTS

Hydrogen-exchange labeling coupled with NMR was used to monitor the formation of stable hydrogen-bonded and solvent-excluded structure in horse cytochrome c (cyt c). Protection was measured using a hydrogen exchange/folding competition protocol at variable pH and short competition time (2 ms). Protection factors of threefold to eightfold were observed in all three alpha helices of cyt c, whereas other regions showed no significant protection. This suggests that the compact states that are present contain segments of marginally stable hydrogen-bonded structure. When the intermediate(s) are destabilized, only amide protons from Cys14, Ala15 and His18 show significant protection, indicating a region of persistent residual structure near the covalently bound heme group in the unfolded protein. Fluorescence-detected stopped-flow studies showed that the maximum protection factor in the early intermediate is consistent with its unfolding equilibrium constant.

CONCLUSIONS

Together with previous fluorescence and CD results, the observed pattern of amide protection is consistent with the early formation of an alpha-helical core domain in an ensemble of compact states, indicating that efficient folding is facilitated by stepwise acquisition of native structural elements. These specific early interactions are established on the sub-millisecond time scale, prior to the rate-limiting step for folding.

Protein folding: matching theory and experiment

The impact of folding funnels and folding simulations on the way experimentalists interpret results is examined. The image of the transition state has changed from a unique species that has a strained configuration, with a correspondingly high free energy, to a more ordinary folding intermediate, whose balance between limited conformational entropy and stabilizing contacts places it at the top of the free energy barrier. Evidence for a broad transition barrier comes from studies showing that mutations can change the position of the barrier. The main controversial issue now is whether populated folding intermediates are productive on-pathway intermediates or dead-end traps. Direct experimental evidence is needed. Theories suggesting that populated intermediates are trapped in a glasslike state are usually based on mechanisms which imply that trapping would only be extremely short-lived (e.g., nanoseconds) in water at 25 degrees C. There seems to be little experimental evidence for long-lived trapping in monomers, if folding aggregates are excluded. On the other hand, there is good evidence for kinetic trapping in dimers. alpha-Helix formation is currently the fastest known process in protein folding, and incipient helices are present at the start of folding. Fast helix formation has the effect of narrowing drastically the choice of folding routes. Thus helix formation can direct folding. It changes the folding metaphor from pouring liquid down a folding funnel to a train leaving a switchyard with only a few choices of exit tracks.

A distinct 14 residue site triggers coiled-coil formation in cortexillin I

We have investigated the process of the assembly of the Dictyostelium discoideum cortexillin I oligomerization domain (Ir) into a tightly packed, two-stranded, parallel coiled-coil structure using a variety of recombinant polypeptide chain fragments. The structures of these Ir fragments were analyzed by circular dichroism spectroscopy, analytical ultracentrifugation and electron microscopy. Deletion mapping identified a distinct 14 residue site within the Ir coiled coil, Arg311-Asp324, which was absolutely necessary for dimer formation, indicating that heptad repeats alone are not sufficient for stable coiled-coil formation. Moreover, deletion of the six N-terminal heptad repeats of Ir led to the formation of a four- rather than a two-helix structure, suggesting that the full-length cortexillin I coiled-coil domain behaves as a cooperative folding unit. Most interestingly, a 16 residue peptide containing the distinct coiled-coil 'trigger' site Arg311-Asp324 yielded approximately 30% helix formation as monomer, in aqueous solution. pH titration and NaCl screening experiments revealed that the peptide's helicity depends strongly on pH and ionic strength, indicating that electrostatic interactions by charged side chains within the peptide are critical in stabilizing its monomer helix. Taken together, these findings demonstrate that Arg311-Asp324 behaves as an autonomous helical folding unit and that this distinct Ir segment controls the process of coiled-coil formation of cortexillin I.

Diffusion-controlled DNA recognition by an unfolded, monomeric bZIP transcription factor

Basic leucine zipper (bZIP) transcription factors are dimers that recognize mainly palindromic DNA sites. It has been assumed that bZIP factors have to form a dimer in order to bind to their target DNA. We find that DNA binding of both monomeric and dimeric bZIP transcription factor GCN4 is diffusion-limited and that, therefore, the rate of dimerization of the bZIP domain does not affect the rate of DNA recognition and GCN4 need not dimerize in order to bind to its specific DNA site. The results have implications for the mechanism by which bZIP transcription factors find their target sites for transcriptional regulation.

Contrasting roles for symmetrically disposed beta-turns in the folding of a small protein

To investigate the role of turns in protein folding, we have characterized the effects of combinatorial and site-directed mutations in the two beta-turns of peptostreptococcal protein L on folding thermodynamics and kinetics. Sequences of folded variants recovered from combinatorial libraries using a phase display selection method were considerably more variable in the second turn than in the first turn. These combinatorial mutants as well as strategically placed point mutants in the two turns had a similar range of thermodynamic stabilities, but strikingly different folding kinetics. A glycine to alanine substitution in the second beta-turn increased the rate of unfolding more than tenfold but had little effect on the rate of folding, while mutation of a symmetrically disposed glycine residue in the first turn had little effect on unfolding but slowed the rate of folding nearly tenfold. These results demonstrate that the role of beta-turns in protein folding is strongly context-dependent, and suggests that the first turn is formed and the second turn disrupted in the folding transition state.

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.

Ultrafast signals in protein folding and the polypeptide contracted state

To test the significance of ultrafast protein folding signals (<<1 msec), we studied cytochrome c (Cyt c) and two Cyt c fragments with major C-terminal segments deleted. The fragments remain unfolded under all conditions and so could be used to define the unfolded baselines for protein fluorescence and circular dichroism (CD) as a function of denaturant concentration. When diluted from high to low denaturant in kinetic folding experiments, the fragments readjust to their new baseline values in a "burst phase" within the mixing dead time. The fragment burst phase reflects a contraction of the polypeptide from a more extended unfolded condition at high denaturant to a more contracted unfolded condition in the poorer, low denaturant solvent. Holo Cyt c exhibits fluorescence and CD burst phase signals that are essentially identical to the fragment signals over the whole range of final denaturant concentrations, evidently reflecting the same solvent-dependent, relatively nonspecific contraction and not the formation of a specific folding intermediate. The significance of fast folding signals in Cyt c and other proteins is discussed in relation to the hypothesis of an initial rate-limiting search-nucleation-collapse step in protein folding [Sosnick, T. R., Mayne, L. & Englander, S. W. (1996) Proteins Struct. Funct. Genet. 24, 413-426].

Modeling transmembrane helical oligomers

Recently, methods for the analysis and design of water-soluble, oligomeric bundles of alpha helices, including coiled coils, have reached a high level of sophistication. These same methods may now be applied to transmembrane helical bundles. Studies of the transmembrane domains of glycophorin, phospholamban, and the M2 protein from influenza A virus exemplify this general approach.

Thermal unfolding in a GCN4-like leucine zipper: 13C alpha NMR chemical shifts and local unfolding curves

13C alpha chemical shifts and site-specific unfolding curves are reported for 12 sites on a 33-residue, GCN4-like leucine zipper peptide (GCN4-lzK), ranging over most of the chain and sampling most heptad positions. Data were derived from NMR spectra of nine synthetic, isosequential peptides bearing 99% 13C alpha at sites selected to avoid spectral overlap in each peptide. At each site, separate resonances appear for unfolded and folded forms, and most sites show resonances for two folded forms near room temperature. The observed chemical shifts suggest that 1) urea-unfolded GCN4-lzK chains are randomly coiled; 2) thermally unfolded chains include significant transient structure, except at the ends; 3) the coiled-coli structure in the folded chains is atypical near the C-terminus; 4) only those interior sites surrounded by canonical interchain salt bridges fail to show two folded forms. Local unfolding curves, obtained from integrated resonance intensities, show that 1) sites differ in structure content and in melting temperature, so the equilibrium population must comprise more than two molecular conformations; 2) there is significant end-fraying, even at the lowest temperatures, but thermal unfolding is not a progressive unwinding from the ends; 3) residues 9-16 are in the lowest melting region; 4) heptad position does not dictate stability; 5) significant unfolding occurs below room temperature, so the shallow, linear decline in backbone CD seen there has conformational significance. It seems that only a relatively complex array of conformational states could underlie these findings.

Structure of bacteriophage T4 fibritin: a segmented coiled coil and the role of the C-terminal domain

BACKGROUND

Oligomeric coiled-coil motifs are found in numerous protein structures; among them is fibritin, a structural protein of bacteriophage T4, which belongs to a class of chaperones that catalyze a specific phage-assembly process. Fibritin promotes the assembly of the long tail fibers and their subsequent attachment to the tail baseplate; it is also a sensing device that controls the retraction of the long tail fibers in adverse environments and, thus, prevents infection. The structure of fibritin had been predicted from sequence and biochemical analyses to be mainly a triple-helical coiled coil. The determination of its structure at atomic resolution was expected to give insights into the assembly process and biological function of fibritin, and the properties of modified coiled-coil structures in general.

RESULTS

The three-dimensional structure of fibritin E, a deletion mutant of wild-type fibritin, was determined to 2.2 A resolution by X-ray crystallography. Three identical subunits of 119 amino acid residues form a trimeric parallel coiled-coil domain and a small globular C-terminal domain about a crystallographic threefold axis. The coiled-coil domain is divided into three segments that are separated by insertion loops. The C-terminal domain, which consists of 30 residues from each subunit, contains a beta-propeller-like structure with a hydrophobic interior.

CONCLUSIONS

The residues within the C-terminal domain make extensive hydrophobic and some polar intersubunit interactions. This is consistent with the C-terminal domain being important for the correct assembly of fibritin, as shown earlier by mutational studies. Tight interactions between the C-terminal residues of adjacent subunits counteract the latent instability that is suggested by the structural properties of the coiled-coil segments. Trimerization is likely to begin with the formation of the C-terminal domain which subsequently initiates the assembly of the coiled coil. The interplay between the stabilizing effect of the C-terminal domain and the labile coiled-coil domain may be essential for the fibritin function and for the correct functioning of many other alpha-fibrous proteins.

Folding funnels and energy landscapes of larger proteins within the capillarity approximation

The characterization of protein-folding kinetics with increasing chain length under various thermodynamic conditions is addressed using the capillarity picture in which distinct spatial regions of the protein are imagined to be folded or trapped and separated by interfaces. The quantitative capillarity theory is based on the nucleation theory of first-order transitions and the droplet analysis of glasses and random magnets. The concepts of folding funnels and rugged energy landscapes are shown to be applicable in the large size limit just as for smaller proteins. An ideal asymptotic free-energy profile as a function of a reaction coordinate measuring progress down the funnel is shown to be quite broad. This renders traditional transition state theory generally inapplicable but allows a diffusive picture with a transition-state region to be used. The analysis unifies several scaling arguments proposed earlier. The importance of fluctuational fine structure both to the free-energy profile and to the glassy dynamics is highlighted. The fluctuation effects lead to a very broad trapping-time distribution. Considerations necessary for understanding the crossover between the mean field and capillarity pictures of the energy landscapes are discussed. A variety of mechanisms that may roughen the interfaces and may lead to a complex structure of the transition-state ensemble are proposed.

Protein folding: the endgame

The last stage of protein folding, the "endgame," involves the ordering of amino acid side-chains into a well defined and closely packed configuration. We review a number of topics related to this process. We first describe how the observed packing in protein crystal structures is measured. Such measurements show that the protein interior is packed exceptionally tightly, more so than the protein surface or surrounding solvent and even more efficiently than crystals of simple organic molecules. In vitro protein folding experiments also show that the protein is close-packed in solution and that the tight packing and intercalation of side-chains is a final and essential step in the folding pathway. These experimental observations, in turn, suggest that a folded protein structure can be described as a kind of three-dimensional jigsaw puzzle and that predicting side-chain packing is possible in the sense of solving this puzzle. The major difficulty that must be overcome in predicting side-chain packing is a combinatorial "explosion" in the number of possible configurations. There has been much recent progress towards overcoming this problem, and we survey a variety of the approaches. These approaches differ principally in whether they use ab initio (physical) or more knowledge-based methods, how they divide up and search conformational space, and how they evaluate candidate configurations (using scoring functions). The accuracy of side-chain prediction depends crucially on the (assumed) positioning of the main-chain. Methods for predicting main-chain conformation are, in a sense, not as developed as that for side-chains. We conclude by surveying these methods. As with side-chain prediction, there are a great variety of approaches, which differ in how they divide up and search space and in how they score candidate conformations.

Modeling protein folding: the beauty and power of simplicity

It is argued that simplified models capture key features of protein stability and folding, whereas more detailed models may be more appropriate for protein structure prediction. A brief overview of experimental and theoretical results is presented that corroborates these points. I argue that statistical models capture the key principle of protein stability-cooperativity- and therefore provide a reasonable estimate of protein free energy whereas more detailed but less physically transparent calculations fail to do so. I also explain that the previously published claim that simple models give predictions that are inconsistent with experiments on polypeptide block-copolymers is based on incomplete analysis of such experiment.