A tale of two secondary structure elements: when a beta-hairpin becomes an alpha-helix

Unveiling the dimer/monomer propensities of Smad MH1-DNA complexes

•

MH1 domains of BMP-activated Smads populate a dimer-monomer equilibrium.

•

Swapping loop1 between BMP- and TGFβ Smads turns dimers into monomers and vice versa.

•

BMP-responsive genomic regions have a lower average count of motifs than TGFβ ones.

Smad transcription factors are the main downstream effectors of the Transforming growth factor β superfamily (TGFβ) signalling network. The DNA complexes determined here by X-ray crystallography for the Bone Morphogenetic Proteins (BMP) activated Smad5 and Smad8 proteins reveal that all MH1 domains bind [GGC(GC)|(CG)] motifs similarly, although TGFβ-activated Smad2/3 and Smad4 MH1 domains bind as monomers whereas Smad1/5/8 form helix-swapped dimers. Dimers and monomers are also present in solution, as revealed by NMR. To decipher the characteristics that defined these dimers, we designed chimeric MH1 domains and characterized them using X-ray crystallography. We found that swapping the loop1 between TGFβ- and BMP- activated MH1 domains switches the dimer/monomer propensities. When we scanned the distribution of Smad-bound motifs in ChIP-Seq peaks (Chromatin immunoprecipitation followed by high-throughput sequencing) in Smad-responsive genes, we observed specific site clustering and spacing depending on whether the peaks correspond to BMP- or TGFβ-responsive genes. We also identified significant correlations between site distribution and monomer or dimer propensities. We propose that the MH1 monomer or dimer propensity of Smads contributes to the distinct motif selection genome-wide and together with the MH2 domain association, help define the composition of R-Smad/Smad4 trimeric complexes.

Exploring the sequence fitness landscape of a bridge between protein folds

Most foldable protein sequences adopt only a single native fold. Recent protein design studies have, however, created protein sequences which fold into different structures apon changes of environment, or single point mutation, the best characterized example being the switch between the folds of the GA and GB binding domains of streptococcal protein G. To obtain further insight into the design of sequences which can switch folds, we have used a computational model for the fitness landscape of a single fold, built from the observed sequence variation of protein homologues. We have recently shown that such coevolutionary models can be used to design novel foldable sequences. By appropriately combining two of these models to describe the joint fitness landscape of GA and GB, we are able to describe the propensity of a given sequence for each of the two folds. We have successfully tested the combined model against the known series of designed GA/GB hybrids. Using Monte Carlo simulations on this landscape, we are able to identify pathways of mutations connecting the two folds. In the absence of a requirement for domain stability, the most frequent paths go via sequences in which neither domain is stably folded, reminiscent of the propensity for certain intrinsically disordered proteins to fold into different structures according to context. Even if the folded state is required to be stable, we find that there is nonetheless still a wide range of sequences which are close to the transition region and therefore likely fold switches, consistent with recent estimates that fold switching may be more widespread than had been thought.

While most proteins self-assemble (or “fold”) to a unique three-dimensional structure, a few have been identified that can fold into two distinct structures. These so-called “metamorphic” proteins that can switch folds have attracted a lot of recent interest, and it has been suggested that they may be much more widespread than currently appreciated. We have developed a computational model that captures the propensity of a given protein sequence to fold into either one of two specific structures (GA and GB), in order to investigate which sequences are able to fold to both GA and GB (“switch sequences”), versus just one of them. Our model predicts that there is a large number of switch sequences that could fold into both structures, but also that the most likely such sequences are those for which the folded structures have low stability, in agreement with available experimental data. This also suggests that intrinsically disordered proteins which can fold into different structures on binding may provide an evolutionary path in sequence space between protein folds.

Network Connectivity, Centrality and Fragmentation in the Greek-Key Protein Topology

Understanding and computationally predicting the protein folding process remains one of the most challenging scientific problems and has uniquely garnered the interdisciplinary efforts of researchers from both the biological, chemical, physical and computational disciplines. Previous studies have demonstrated the importance of long-range interactions in guiding the native structure. However, predicting how the native long-range interaction network forms to generate a specific topology from among all other conformations remains unresolved. The present research study conducts an exploratory study to identify amino acids and long-range interactions that have the potential to play a key role in building and maintaining the protein topology. Towards this end, the application of network science is utilized and developed to analyze the structures of a group of proteins that share a common Greek-key topology but differ in sequence, secondary structure and function. We investigate the idea that the residues with high betweeness centrality score are potentially significant in maintaining the protein network and in governing the Greek-key topology. This hypothesis is tested by two different computational methods: through a fragmentation test and by the analysis of diameter impacts. In summary, we find a subset of selected residues in similar geographical positions in all model proteins, which demonstrates the role of these specific residues and regions in governing the Greek-key topology from a network perspective.

Dependence of internal friction on folding mechanism

An outstanding challenge in protein folding is understanding the origin of "internal friction" in folding dynamics, experimentally identified from the dependence of folding rates on solvent viscosity. A possible origin suggested by simulation is the crossing of local torsion barriers. However, it was unclear why internal friction varied from protein to protein or for different folding barriers of the same protein. Using all-atom simulations with variable solvent viscosity, in conjunction with transition-path sampling to obtain reaction rates and analysis via Markov state models, we are able to determine the internal friction in the folding of several peptides and miniproteins. In agreement with experiment, we find that the folding events with greatest internal friction are those that mainly involve helix formation, while hairpin formation exhibits little or no evidence of friction. Via a careful analysis of folding transition paths, we show that internal friction arises when torsion angle changes are an important part of the folding mechanism near the folding free energy barrier. These results suggest an explanation for the variation of internal friction effects from protein to protein and across the energy landscape of the same protein.

Transformation of an alpha-helix peptide into a beta-hairpin induced by addition of a fragment results in creation of a coexisting state

Intrinsic rules of determining the tertiary structure of a protein have been unknown partly because physicochemical factors that contribute to stabilization of a protein structure cannot be represented as a linear combination of local interactions. To clarify the rules on the nonlinear term caused by nonlocal interaction in a protein, we tried to transform a peptide that has a fully helical structure (Target Peptide or TP) into a peptide that has a beta-hairpin structure (Designed Peptide or DP) by adding seven residues to the C terminus of TP. According to analyses of nuclear magnetic resonance measurements, while the beta-hairpin structure is stabilized in some DPs, it is evident that the helical structure observed in TP is also persistent and even extended throughout the length of the molecule. As a result, we have produced a peptide molecule that contains both the alpha-helix and beta-hairpin conformation at an almost equally populated level. The helical structures contained in these DPs were more stable than the helix in TP, suggesting that stabilizing one conformation does not result in destabilizing the other conformation. These DPs can thus be regarded as an isolated peptide version of the chameleon sequence, which has the capability of changing the secondary structure depending on the context of the surrounding environment in a protein structure. The fact that the transformation of one secondary structure caused stabilization of both the original and the induced structure would shed light on the mechanism of protein folding.

Folding mechanisms of proteins with high sequence identity but different folds

The problem of how a protein folds from a linear chain of amino acids to the three-dimensional structure necessary for function is often investigated using proteins with a low degree of sequence identity that adopt different folds. The design of pairs of proteins with a high degree of sequence identity but different folds offers the opportunity for a complementary study; in two highly similar sequences, which residues are the most important in directing folding to a particular structure? Here we use molecular dynamics simulations to characterize the folding-unfolding pathways of a pair of proteins designed by Bryan and co-workers [Alexander, P. A., et al. (2005) Biochemistry 44, 14045-14054; He, Y. N., et al. (2005) Biochemistry 44, 14055-14061]. Despite being 59% identical, the two protein sequences fold to two different structures. The first sequence folds to the alpha+beta protein G structure and the second to the all-alpha-helical protein A structure. We show that the final protein structure is determined early along the folding pathway. In folding to the protein G structure, the single alpha-helix (alpha1) and the beta3-beta4 turn fold early. Formation of the hairpin turn essentially prevents folding to helical structure in this region of the protein. This early structure is then consolidated by formation of long-range hydrophobic interactions between alpha1 and the beta3-beta4 turn. The protein A sequence differs both in the residues that form the beta3-beta4 turn and also in many of the residues that form the early hydrophobic interactions in the protein G structure. Instead, in the protein A sequence, a more hierarchical mechanism is observed, with helices folding before many of the tertiary interactions are formed. We find that small, but critical, sequence differences determine the topology of the protein early along the folding pathway, which help to explain the process by which one fold can evolve into another.

Engineering Diverse Changes in β-Turn Propensities in the N-Terminal β-Hairpin of Ubiquitin Reveals Significant Effects on Stability and Kinetics but a Robust Folding Transition State

Using the N-terminal 17-residue beta-hairpin of ubiquitin as a "host" for mutational studies, we have investigated the influence of the beta-turn sequence on protein stability and folding kinetics by replacing the native G-bulged turn (TLTGK) with more flexible analogues (TG3K and TG5K) and a series of four-residue type I' beta-turn sequences, commonly found in beta-hairpins. Although a statistical analysis of type I' turns demonstrates residue preferences at specific sites, the frequency of occurrence appears to only broadly correlate with experimentally determined protein stabilities. The subsequent engineering of context-dependent non-native tertiary contacts involving turn residues is shown to produce large changes in stability. Relatively few point mutations have been described that probe secondary structure formation in ubiquitin in a manner that is independent of tertiary contacts. To this end, we have used the more rigorous rate-equilibrium free energy relationship (Leffler analysis), rather than the two-point phi value analysis, to show for a family of engineered beta-turn mutants that stability (range of approximately 20 kJ/mol) and folding kinetics (190-fold variation in refolding rate) are linearly correlated (alpha(f) = 0.74 +/- 0.08). The data are consistent with a transition state that is robust with regard to a wide range of statistically favored and disfavored beta-turn mutations and implicate a loosely assembled beta-hairpin as a key template in transition state stabilization with the beta-turn playing a central role.

Engineering Stabilising β-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.

Effects of segment substitution on the structure and stability of immunoglobulin G binding domain of streptococcal protein G

Structural formation of segments plays pivotal roles in protein folding and stability, but how the segment influences the structural ensemble remains elusive. We engineered two hybrid proteins by replacing the central helical segment of immunoglobulin G binding domain of streptococcal protein G with an alpha-helix or beta(2)-strand element of a structural homologue, the immunoglobulin G binding domain of streptococcal protein L. The results show that substitution by the alpha-helical sequence retains a folded structure predominantly with a three-stranded beta-sheet but slightly destabilizes the compact ensemble, while substitution by the beta(2)-strand sequence completely destroys the structural formation. The finding implies that the local segment may influence the tertiary structure and overall stability, and the tertiary interactions may modulate structural formation of the segment, which might be considered when studying protein folding, prediction, design, and engineering.

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.

The effect of an ionic detergent on the natively unfolded beta-dystroglycan ectodomain and on its interaction with alpha-dystroglycan

Dystroglycan (DG) is an adhesion complex, expressed in a wide variety of tissues, formed by an extracellular and a transmembrane subunit, alpha-DG and beta-DG, respectively, interacting noncovalently. Recently, we have shown that the recombinant ectodomain of beta-DG, beta-DG(654-750), behaves as a natively unfolded protein, as it is able to bind the C-terminal domain of alpha-DG, while not displaying a defined structural organization. We monitored the effect of a commonly used denaturing agent, the anionic detergent sodium dodecylsulphate (SDS), on beta-DG(654-750) using a number of biophysical techniques. Very low concentrations of SDS (< or =2 mM) affect both tryptophan fluorescence and circular dichroism of beta-DG, and significantly perturb the interaction with the alpha-DG subunit as shown by solid-phase binding assays and fluorescence titrations in solution. This result confirms, as recently proposed for natively unfolded proteins, that beta-DG(654-750) exists in a native state, which is crucial to fulfill its biological function. Two-dimensional NMR analysis shows that SDS does not induce any evident conformational rearrangement within the ectodomain of beta-DG. Its first 70 amino acids, which show a lower degree of mobility, interact with the detergent, but this does not change the amount of secondary structure, whereas the highly flexible and mobile C-terminal region of beta-DG(654-750) remains largely unaffected, even at a very high SDS concentration (up to 50 mM). Our data indicate that SDS can be used as a useful tool for investigating natively unfolded proteins, and confirm that the beta-DG ectodomain is an interesting model system.

Stability and Folding Kinetics of a Ubiquitin Mutant with a Strong Propensity for Nonnative β-Hairpin Conformation in the Unfolded State

A F45W mutant of yeast ubiquitin has been used as a model system to examine the effects of nonnative local interactions on protein folding and stability. Mutating the native TLTGK G-bulged type I turn in the N-terminal beta-hairpin to NPDG stabilizes a nonnative beta-strand alignment in the isolated peptide fragment. However, NMR structural analysis of the native and mutant proteins shows that the NPDG mutant is forced to adopt the native beta-strand alignment and an unfavorable type I NPDG turn. The mutant is significantly less stable (approximately 9 kJ mol(-1)) and folds 30 times slower than the native sequence, demonstrating that local interactions can modulate protein stability and that attainment of a nativelike beta-hairpin conformation in the transition state ensemble is frustrated by the turn mutations. Surprising, alcoholic cosolvents [5-10% (v/v) TFE] are shown to accelerate the folding rate of the NPDG mutant. We conclude, backed-up by NMR data on the peptide fragments, that even though nonnative states in the denatured ensemble are highly populated and their stability further enhanced in the presence of cosolvents, the simultaneous increase in the proportion of nativelike secondary structure (hairpin or helix), in rapid equilibrium with nonnative states, is sufficient to accelerate the folding process. It is evident that modulating local interactions and increasing nonnative secondary structure propensities can change protein stability and folding kinetics. However, nonlocal contacts formed in the global cooperative folding event appear to determine structural specificity.

Coupled prediction of protein secondary and tertiary structure

The strong coupling between secondary and tertiary structure formation in protein folding is neglected in most structure prediction methods. In this work we investigate the extent to which nonlocal interactions in predicted tertiary structures can be used to improve secondary structure prediction. The architecture of a neural network for secondary structure prediction that utilizes multiple sequence alignments was extended to accept low-resolution nonlocal tertiary structure information as an additional input. By using this modified network, together with tertiary structure information from native structures, the Q3-prediction accuracy is increased by 7-10% on average and by up to 35% in individual cases for independent test data. By using tertiary structure information from models generated with the ROSETTA de novo tertiary structure prediction method, the Q3-prediction accuracy is improved by 4-5% on average for small and medium-sized single-domain proteins. Analysis of proteins with particularly large improvements in secondary structure prediction using tertiary structure information provides insight into the feedback from tertiary to secondary structure.

Mechanism of fast protein folding

An explosion of in vitro experimental data on the folding of proteins has revealed many examples of folding in the millisecond or faster timescale, often occurring in the absence of stable intermediate states. We review experimental methods for measuring fast protein folding kinetics, and then discuss various analytical models used to interpret these data. Finally, we classify general mechanisms that have been proposed to explain fast protein folding into two catagories, heterogeneous and homogeneous, reflecting the nature of the transition state. One heterogeneous mechanism, the diffusion-collision mechanism, can be used to interpret experimental data for a number of proteins.

Importance of secondary structural specificity determinants in protein folding: insertion of a native beta-sheet sequence into an alpha-helical coiled-coil

To examine how a short secondary structural element derived from a native protein folds when in a different protein environment, we inserted an 11-residue beta-sheet segment (cassette) from human immunoglobulin fold, Fab new, into an alpha-helical coiled-coil host protein (cassette holder). This de novo design protein model, the structural cassette mutagenesis (SCM) model, allows us to study protein folding principles involving both short- and long-range interactions that affect secondary structure stability and conformation. In this study, we address whether the insertion of this beta-sheet cassette into the alpha-helical coiled-coil protein would result in conformational change nucleated by the long-range tertiary stabilization of the coiled-coil, therefore overriding the local propensity of the cassette to form beta-sheet, observed in its native immunoglobulin fold. The results showed that not only did the nucleating helices of the coiled-coil on either end of the cassette fail to nucleate the beta-sheet cassette to fold with an alpha-helical conformation, but also the entire chimeric protein became a random coil. We identified two determinants in this cassette that prevented coiled-coil formation: (1) a tandem dipeptide NN motif at the N-terminal of the beta-sheet cassette, and (2) the hydrophilic Ser residue, which would be buried in the hydrophobic core if the coiled-coil structure were to fold. By amino acid substitution of these helix disruptive residues, that is, either the replacement of the NN motif with high helical propensity Ala residues or the substitution of Ser with Leu to enhance hydrophobicity, we were able to convert the random coil chimeric protein into a fully folded alpha-helical coiled-coil. We hypothesized that this NN motif is a "secondary structural specificity determinant" which is very selective for one type of secondary structure and may prevent neighboring residues from adopting an alternate protein fold. These sequences with secondary structural specificity determinants have very strong local propensity to fold into a specific secondary structure and may affect overall protein folding by acting as a folding initiation site.

Structural mimicry of proline kinks: tertiary packing interactions support local structural distortions

Proline residues in the helical segments of soluble and transmembrane proteins have received special attention from both a structural and functional perspective. A feature of these helices is the structural distortion termed "proline-kink", which has been associated with the presence of the proline residue. However, a recent report on the yeast heat-shock transcription factor of Kluyveromyces lactis (HSF_KL) suggests that these proline-associated deformations can be achieved in the absence of proline residues, thus raising the question of the mechanisms responsible for the structural mimicry of proline-related features. In this study, the specific interactions responsible for the distortion were characterized by comparative analysis of the atomic details of the packing interactions that surround the evolutionarily conserved proline-kink in the alpha2 helix of HSF_KL and a set of 39 structurally related proteins that lacked the distortion. The mechanistic details inferred from this analysis were confirmed with molecular dynamics simulations. The study shows that the packing interactions between the alpha2 and alpha1 helices in HSF_KL are responsible for the stabilization of the conserved kink, whether a proline residue that divides the helix into segments is present or not. The proline-kink can facilitate the formation of tertiary packing interactions that would otherwise not be possible. However, it is the ability to establish differential packing interactions for the helix segments, rather than the structural properties of the proline-kink itself, that emerges as the key factor for the characteristic distortion.

Unspecific hydrophobic stabilization of folding transition states

Here we present a method for determining the inference of non-native conformations in the folding of a small domain, alpha-spectrin Src homology 3 domain. This method relies on the preservation of all native interactions after Tyr/Phe exchanges in solvent-exposed, contact-free positions. Minor changes in solvent exposure and free energy of the denatured ensemble are in agreement with the reverse hydrophobic effect, as the Tyr/Phe mutations slightly change the polypeptide hydrophilic/hydrophobic balance. Interestingly, more important Gibbs energy variations are observed in the transition state ensemble (TSE). Considering the small changes induced by the H/OH replacements, the observed energy variations in the TSE are rather notable, but of a magnitude that would remain undetected under regular mutations that alter the folded structure free energy. Hydrophobic residues outside of the folding nucleus contribute to the stability of the TSE in an unspecific nonlinear manner, producing a significant acceleration of both unfolding and refolding rates, with little effect on stability. These results suggest that sectors of the protein transiently reside in non-native areas of the landscape during folding, with implications in the reading of phi values from protein engineering experiments. Contrary to previous proposals, the principle that emerges is that non-native contacts, or conformations, could be beneficial in evolution and design of some fast folding proteins.

Hints of nonhierarchical folding of acidic fibroblast growth factor

We have analyzed by circular dichroism (CD) and proton nuclear magnetic resonance (NMR) the helical propensity of the all-beta protein acidic fibroblast growth factor (aFGF) and two peptides corresponding to beta-strand 8 (beta8 peptide, amino acids 95-107) and the beta-strand 8/turn/beta-strand 9 hairpin (beta8/9 peptide, amino acids 95-114), which has been involved in receptor binding. A secondary structure prediction of aFGF carried out by several procedures labels the 95-104 sequence as predominantly alpha-helical. A titration of aFGF with 2,2,2-trifluoroethanol (TFE) induces a change in the far-UV CD spectrum of the protein giving rise to a prominent alpha-helical shape (22% alpha-helix). The cooperativity of the transition and the moderate TFE concentrations used (midpoint at 24%) suggest that the effect of TFE is specific. Moreover, a titration performed at pH 2 yields a higher amount of alpha-helix (55%) at a smaller TFE concentration. Synthetic peptides containing the beta8 and beta8/9 sequences display a random coil conformation at pH 7 but acquire alpha-helical structure in the presence of TFE, methanol, and SDS micelles. At pH below 3.0 a significant amount (20-30%) of alpha-helical conformation is present in both the beta8 and beta8/9 peptides even in the absence of other solvent additives. The secondary structure of the peptides was determined by proton nuclear magnetic resonance (1H NMR). These results suggest that the 95-114 sequence of aFGF has helical propensity and that the protein may fold nonhierarchically in the early steps of folding, acquiring its final beta-structure by a later interaction with the rest of the polypeptide.

Elucidating protein secondary structures using alpha-carbon recurrence quantifications

Secondary structures of proteins were studied by recurrence quantification analysis (RQA). High-resolution, 3-dimensional coordinates of alpha-carbon atoms comprising a set of 68 proteins were downloaded from the Protein Data Bank. By fine-tuning four recurrence parameters (radius, line, residue, separation), it was possible to establish excellent agreement between percent contribution of alpha-helix and beta-sheet structures determined independently by RQA and that of the DSSP algorithm (Define Secondary Structure of Proteins). These results indicate that there is an equivalency between these two techniques, which are based upon totally different pattern recognition strategies. RQA enhances qualitative contact maps by quantifying the arrangements of recurrent points of alpha carbons close in 3-dimensional space. For example, the radius was systematically increased, moving the analysis beyond local alpha-carbon neighborhoods in order to capture super-secondary and tertiary structures. However, differences between proteins could only be detected within distances up to about 6-11 A, but not higher. This result underscores the complexity of alpha-carbon spacing when super-secondary structures appear at larger distances. Finally, RQA-defined secondary structures were found to be robust against random displacement of alpha carbons upwards of 1 A. This finding has potential import for the dynamic functions of proteins in motion.

Protein folding theory: from lattice to all-atom models

This review focuses on recent advances in understanding protein folding kinetics in the context of nucleation theory. We present basic concepts such as nucleation, folding nucleus, and transition state ensemble and then discuss recent advances and challenges in theoretical understanding of several key aspects of protein folding kinetics. We cover recent topology-based approaches as well as evolutionary studies and molecular dynamics approaches to determine protein folding nucleus and analyze other aspects of folding kinetics. Finally, we briefly discuss successful all-atom Monte-Carlo simulations of protein folding and conclude with a brief outlook for the future.

Understanding the sequence determinants of conformational switching using protein design

An important goal of protein design is to understand the forces that stabilize a particular fold in preference to alternative folds. Here, we describe an extension of earlier studies in which we successfully designed a stable, native-like helical protein that is 50% identical in sequence to a predominantly beta-sheet protein, the B1 domain of Streptococcal IgG-binding protein G. We report the characteristics of a series of variants of our original design that have even higher sequence identity to the B1 domain. Their properties illustrate the extent to which protein stability and conformation can be modulated through careful manipulation of key amino acid residues. Our results have implications for understanding conformational change phenomena of central biological importance and in probing the malleability of the sequence/structure relationship.

Fitting peptides into the RNA world

The structures of several peptide-RNA complexes have been reported in the past year, underscoring the diverse nature of RNA structure and protein interactions. In general, specific peptide conformations are stabilized by the surrounding RNA framework; this is strikingly similar to how peptides are stabilized upon interaction with proteins.

Is folding of beta-lactoglobulin non-hierarchic? Intermediate with native-like beta-sheet and non-native alpha-helix

The refolding of beta-lactoglobulin, a beta-barrel protein consisting of beta strands betaA-betaI and one major helix, is unusual because non-native alpha-helices are formed at the beginning of the process. We studied the refolding kinetics of bovine beta-lactoglobulin A at pH 3 using the stopped-flow circular dichroism and manual H/(2)H exchange pulse labeling coupled with heteronuclear NMR. The protection pattern from the H/(2)H exchange of the native state indicated the presence of a stable hydrophobic core consisting of betaF, betaG and betaH strands. The protection pattern of the kinetic intermediate obtained about one second after initiating the reaction was compared with that of the native state. In this relatively late kinetic intermediate, which still contains some non-native helical structure, the disulfide-bonded beta-hairpin made up of betaG and betaH strands was formed, but the rest of the molecule was fluctuating, where the non-native alpha-helices may reside. Subsequently, the core beta-sheet extends, accompanied by a further alpha-helix to beta-sheet transition. Thus, the refolding of beta-lactoglobulin exhibits two elements: the critical role of the core beta-sheet is consistent with the hierarchic mechanism, whereas the alpha-helix to beta-sheet transition suggests the non-hierarchic mechanism.