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

Dilute phase oligomerization can oppose phase separation and modulate material properties of a ribonucleoprotein condensate

SignificanceA large subclass of biomolecular condensates are linked to RNA regulation and are known as ribonucleoprotein (RNP) bodies. While extensive work has identified driving forces for biomolecular condensate formation, relatively little is known about forces that oppose assembly. Here, using a fungal RNP protein, Whi3, we show that a portion of its intrinsically disordered, glutamine-rich region modulates phase separation by forming transient alpha helical structures that promote the assembly of dilute phase oligomers. These oligomers detour Whi3 proteins from condensates, thereby impacting the driving forces for phase separation, the protein-to-RNA ratio in condensates, and the material properties of condensates. Our findings show how nanoscale conformational and oligomerization equilibria can influence mesoscale phase equilibria.

Hierarchical cascades of instability govern the mechanics of coiled coils: helix unfolding precedes coil unzipping

Coiled coils are a fundamental emergent motif in proteins found in structural biomaterials, consisting of α-helical secondary structures wrapped in a supercoil. A fundamental question regarding the thermal and mechanical stability of coiled coils in extreme environments is the sequence of events leading to the disassembly of individual oligomers from the universal coiled-coil motifs. To shed light on this phenomenon, here we report atomistic simulations of a trimeric coiled coil in an explicit water solvent and investigate the mechanisms underlying helix unfolding and coil unzipping in the assembly. We employ advanced sampling techniques involving steered molecular dynamics and metadynamics simulations to obtain the free-energy landscapes of single-strand unfolding and unzipping in a three-stranded assembly. Our comparative analysis of the free-energy landscapes of instability pathways shows that coil unzipping is a sequential process involving multiple intermediates. At each intermediate state, one heptad repeat of the coiled coil first unfolds and then unzips due to the loss of contacts with the hydrophobic core. This observation suggests that helix unfolding facilitates the initiation of coiled-coil disassembly, which is confirmed by our 2D metadynamics simulations showing that unzipping of one strand requires less energy in the unfolded state compared with the folded state. Our results explain recent experimental findings and lay the groundwork for studying the hierarchical molecular mechanisms that underpin the thermomechanical stability/instability of coiled coils and similar protein assemblies.

pH-jump induced leucine zipper folding beyond the diffusion limit

The folding of a pH-sensitive leucine zipper, that is, a GCN4 mutant containing eight glutamic acid residues, has been investigated. A pH-jump induced by a caged proton (o-nitrobenzaldehyde, oNBA) is employed to initiate the process, and time-resolved IR spectroscopy of the amide I band is used to probe it. The experiment has been carefully designed to minimize the buffer capacity of the sample solution so that a large pH jump can be achieved, leading to a transition from a completely unfolded to a completely folded state with a single laser shot. In order to eliminate the otherwise rate-limiting diffusion-controlled step of the association of two peptides, they have been covalently linked. The results for the folding kinetics of the cross-linked peptide are compared with those of an unlinked peptide, which reveals a detailed picture of the folding mechanism. That is, folding occurs in two steps, one on an ∼1-2 μs time scale leading to a partially folded α-helix even in the monomeric case and a second one leading to the final coiled-coil structure on distinctively different time scales of ∼30 μs for the cross-linked peptide and ∼200 μs for the unlinked peptide. By varying the initial pH, it is found that the folding mechanism is consistent with a thermodynamic two-state model, despite the fact that a transient intermediate is observed in the kinetic experiment.

Emergence of structure through protein-protein interactions and pH changes in dually predicted coiled-coil and disordered regions of centrosomal proteins

Human centrosomal proteins show a significant, 3.5 fold, bias to be both unstructured and coiled-coils with respect to generic human proteins, based on results from state of the art bioinformatics tools. We hypothesize that this bias means that these proteins adopt an ensemble of disordered and partially helical conformations, with the latter becoming stabilized when these proteins form complexes. Characterization of the structural properties of 13 peptides from 10 different centrosomal proteins ranging in size from 20 to 61 residues by biophysical methods led us to confirm our hypothesis in most cases. Interestingly, the secondary structure adopted by most of these peptides becomes stabilized at acidic pH and it is concentration dependent. For two of them, PIK3R1(453-513) and BRCA1(1253-1273), we observed not only the stabilization of helical structure through self-association, but also the presence of β-structures linked to the formation of high molecular weight oligomers. These oligomers are the predominant forms detected by CD, but unobservable by liquid state NMR. BRCA1(1397-1424) and MAP3K11(396-441) populate helical structures that can also self-associate at pH3 through oligomeric species. Four peptides, derived from three proteins, namely CCNA2(103-123), BRCA1(1253-1273), BRCA1(1397-1424) and PIK3R1(453-513), can form intermolecular associations that are concomitant with alpha or beta structure stabilization. The self-phosphorylation previously described for the kinase NEK2 did not lead to any stabilization in the peptide's structure of NEK2(303-333), NEK2(341-361), and NEK2(410-430). Based on these results, obtained from a series of peptides derived from a significant number of different centrosomal proteins, we propose that conformational polymorphism, modulated by intermolecular interactions is a general property of centrosomal proteins.

Conformational dynamics is more important than helical propensity for the folding of the all α-helical protein Im7

Im7 folds via an on-pathway intermediate that contains three of the four native α-helices. The missing helix, helix III, is the shortest and its failure to be formed until late in the pathway is related to frustration in the structure. Im7H3M3, a 94-residue variant of the 87-residue Im7 in which helix III is the longest of the four native helices, also folds via an intermediate. To investigate the structural basis for this we calculated the frustration in the structure of Im7H3M3 and used NMR to investigate its dynamics. We found that the native state of Im7H3M3 is highly frustrated and in equilibrium with an intermediate state that lacks helix III, similar to Im7. Model-free analysis identified residues with chemical exchange contributions to their relaxation that aligned with the residues predicted to have highly frustrated interactions, also like Im7. Finally, we determined properties of urea-denatured Im7H3M3 and identified four clusters of interacting residues that corresponded to the α-helices of the native protein. In Im7 the cluster sizes were related to the lengths of the α-helices with cluster III being the smallest but in Im7H3M3 cluster III was also the smallest, despite this region forming the longest helix in the native state. These results suggest that the conformational properties of the urea-denatured states promote formation of a three-helix intermediate in which the residues that form helix III remain non-helical. Thus it appears that features of the native structure are formed early in folding linked to collapse of the unfolded state.

Nonequilibrium dynamics of helix reorganization observed by transient 2D IR spectroscopy

The relaxation of helical structures very close to equilibrium is observed via transient 2D IR spectroscopy. An initial distribution of synthetically distorted helices having an unnatural bridge linking the 10th and 12th residues of an alanine-rich α-helix is released to evolve into the equilibrium distribution of α-helix conformations. The bridge constrains the structure to be slightly displaced from the full α-helix equilibrium near these residues, yet the peptide is not unfolded completely. The release is accomplished by a subpicosecond pulse of UV irradiation. The resulting 2D IR signals are used to obtain snapshots of the ∼100-ps helical conformational reorganization of the distorted dihedral angle and distance between amide units at chemical bond length-scale resolution. The decay rates of the angle between the dipoles, dihedral angles, and distance autocorrelations obtained from molecular dynamics simulations support the experiments, providing evidence that the final helix collapse conforms to linear response theory.

Unique features of the anti-parallel, heterodimeric coiled-coil interaction between methyl-cytosine binding domain 2 (MBD2) homologues and GATA zinc finger domain containing 2A (GATAD2A/p66α)

The methyl-cytosine binding domain 2 (MBD2)-nucleosome remodeling and deacetylase (NuRD) complex recognizes methylated DNA and silences expression of associated genes through histone deacetylase and nucleosome remodeling functions. Our previous structural work demonstrated that a coiled-coil interaction between MBD2 and GATA zinc finger domain containing 2A (GATAD2A/p66α) proteins recruits the chromodomain helicase DNA-binding protein (CHD4/Mi2β) to the NuRD complex and is necessary for MBD2-mediated DNA methylation-dependent gene silencing in vivo (Gnanapragasam, M. N., Scarsdale, J. N., Amaya, M. L., Webb, H. D., Desai, M. A., Walavalkar, N. M., Wang, S. Z., Zu Zhu, S., Ginder, G. D., and Williams, D. C., Jr. (2011) p66α-MBD2 coiled-coil interaction and recruitment of Mi-2 are critical for globin gene silencing by the MBD2-NuRD complex. Proc. Natl. Acad. Sci. U.S.A. 108, 7487-7492). The p66α-MBD2 interaction differs from most coiled-coils studied to date by forming an anti-parallel heterodimeric complex between two peptides that are largely monomeric in isolation. To further characterize unique features of this complex that drive heterodimeric specificity and high affinity binding, we carried out biophysical analyses of MBD2 and the related homologues MBD3, MBD3-like protein 1 (MBD3L1), and MBD3-like protein 2 (MBD3L2) as well as specific mutations that modify charge-charge interactions and helical propensity of the coiled-coil domains. Analytical ultracentrifugation analyses show that the individual peptides remain monomeric in isolation even at 300 μM in concentration for MBD2. Circular dichroism analyses demonstrate a direct correlation between helical content of the coiled-coil domains in isolation and binding affinity for p66α. Furthermore, complementary electrostatic surface potentials and inherent helical content of each peptide are necessary to maintain high-affinity association. These factors lead to a binding affinity hierarchy of p66α for the different MBD2 homologues (MBD2 ≈ MBD3 > MBD3L1 ≈ MBD3L2) and suggest a hierarchical regulatory model in tissue and life cycle stage-specific silencing by NuRD complexes.

Dynamics and mechanisms of coupled protein folding and binding reactions

Protein folding coupled to binding of a specific ligand is frequently observed in biological processes. In recent years numerous studies have addressed the structural properties of the unfolded proteins in the absence of their ligands. Surprisingly few time-resolved investigations on coupled folding and binding reactions have been published up to date and the dynamics and kinetic mechanisms of these processes are still only poorly understood. Especially, it is still unsolved for most systems which conformation of the protein is recognized by the ligand (conformational selection vs. folding-after-binding) and whether the ligand influences the folding kinetics. Here we review experimental methods, kinetic models and time-resolved experimental studies of coupled folding and binding reactions.

The effects of pK(a) tuning on the thermodynamics and kinetics of folding: design of a solvent-shielded carboxylate pair at the a-position of a coiled-coil

The tuning of the pK(a) of ionizable residues plays a critical role in various protein functions, such as ligand-binding, catalysis, and allostery. Proteins harness the free energy of folding to position ionizable groups in highly specific environments that strongly affect their pK(a) values. To investigate the interplay among protein folding kinetics, thermodynamics, and pK(a) modulation, we introduced a pair of Asp residues at neighboring interior positions of a coiled-coil. A single Asp residue was replaced for an Asn side chain at the a-position of the coiled-coil from GCN4, which was also crosslinked at the C-terminus via a flexible disulfide bond. The thermodynamic and kinetic stability of the system was measured by circular dichroism and stopped-flow fluorescence as a function of pH and concentration of guanidine HCl. Both sets of data are consistent with a two-state equilibrium between fully folded and unfolded forms. Distinct pK(a) values of 6.3 and 5.35 are assigned to the first and second protonation of the Asp pair; together they represent an energetic difference of 5 kcal/mol relative to the protonation of two Asp residues with unperturbed pK(a) values. Analysis of the rate data as a function of pH and denaturant concentration allowed calculation of the kinetic constants for the conformational transitions of the peptide with the Asp residues in the doubly protonated, singly protonated, and unprotonated forms. The doubly and singly protonated forms fold rapidly, and a ϕ-value analysis shows that their contribution to folding occurs subsequent to the transition state ensemble for folding. By contrast, the doubly charged state shows a reduced rate of folding and a ϕ-value near 0.5 indicative of a repulsive interaction, and possibly also heterogeneity in the transition state ensemble.

Electrostatic contacts in the activator protein-1 coiled coil enhance stability predominantly by decreasing the unfolding rate

The hypothesis is tested that Jun-Fos activator protein-1 coiled coil interactions are dominated during late folding events by the formation of intricate intermolecular electrostatic contacts. A previously derived cJun-FosW was used as a template as it is a highly stable relative of the wild-type cJun-cFos coiled coil protein (thermal melting temperature = 63 degrees C versus 16 degrees C), allowing kinetic folding data to be readily extracted. An electrostatic mutant, cJun(R)-FosW(E), was created to generate six Arg-Glu interactions at e-g'+1 positions between cJun(R) and FosW(E), and investigations into how their contribution to stability is manifested in the folding pathway were undertaken. The evidence now strongly indicates that the formation of interhelical electrostatic contacts exert their effect predominantly on the coiled coil unfolding/dissociation rate. This has major implications for future antagonist design whereby kinetic rules could be applied to increase the residency time of the antagonist-peptide complex, and therefore significantly increase the efficacy of the antagonist.

The importance of being flexible: the case of basic region leucine zipper transcriptional regulators

Large volumes of protein sequence and structure data acquired by proteomic studies led to the development of computational bioinformatic techniques that made possible the functional annotation and structural characterization of proteins based on their primary structure. It has become evident from genome-wide analyses that many proteins in eukaryotic cells are either completely disordered or contain long unstructured regions that are crucial for their biological functions. The content of disorder increases with evolution indicating a possibly important role of disorder in the regulation of cellular systems. Transcription factors are no exception and several proteins of this class have recently been characterized as premolten/molten globules. Yet, mammalian cells rely on these proteins to control expression of their 30,000 or so genes. Basic region:leucine zipper (bZIP) DNA-binding proteins constitute a major class of eukaryotic transcriptional regulators. This review discusses how conformational flexibility "built" into the amino acid sequence allows bZIP proteins to interact with a large number of diverse molecular partners and to accomplish their manifold cellular tasks in a strictly regulated and coordinated manner.

Mutational analysis of the stability of the H2A and H2B histone monomers

The eukaryotic histone heterodimer H2A-H2B folds through an obligatory dimeric intermediate that forms in a nearly diffusion-limited association reaction in the stopped-flow dead time. It is unclear whether there is partial folding of the isolated monomers before association. To address the possible contributions of structure in the monomers to the rapid association, we characterized H2A and H2B monomers in the absence of their heterodimeric partner. By far-UV circular dichroism, the H2A and H2B monomers are 15% and 31% helical, respectively--significantly less than observed in X-ray crystal structures. Acrylamide quenching of the intrinsic Tyr fluorescence was indicative of tertiary structure. The H2A and H2B monomers exhibit free energies of unfolding of 2.5 and 2.9 kcal mol(-1), respectively; at 10 microM, the sum of the stability of the monomers is approximately 60% of the stability of the native dimer. The helical content, stability, and m values indicate that H2B has a more stable, compact structure than H2A. The monomer m values are larger than expected for the extended histone fold motif, suggesting that the monomers adopt an overly collapsed structure. Stopped-flow refolding-initiated from urea-denatured monomers or the partially folded monomers populated at low denaturant concentrations-yielded essentially identical rates, indicating that monomer folding is productive in the rapid association and folding of the heterodimer. A series of Ala and Gly mutations were introduced into H2A and H2B to probe the importance of helix propensity on the structure and stability of the monomers. The mutational studies show that the central alpha-helix of the histone fold, which makes extensive intermonomer contacts, is structured in H2B but only partially folded in H2A.

Conformational stability and folding mechanisms of dimeric proteins

The folding of multisubunit proteins is of tremendous biological significance since the large majority of proteins exist as protein-protein complexes. Extensive experimental and computational studies have provided fundamental insights into the principles of folding of small monomeric proteins. Recently, important advances have been made in extending folding studies to multisubunit proteins, in particular homodimeric proteins. This review summarizes the equilibrium and kinetic theory and models underlying the quantitative analysis of dimeric protein folding using chemical denaturation, as well as the experimental results that have been obtained. Although various principles identified for monomer folding also apply to the folding of dimeric proteins, the effects of subunit association can manifest in complex ways, and are frequently overlooked. Changes in molecularity typically give rise to very different overall folding behaviour than is observed for monomeric proteins. The results obtained for dimers have provided key insights pertinent to understanding biological assembly and regulation of multisubunit proteins. These advances have set the stage for future advances in folding involving protein-protein interactions for natural multisubunit proteins and unnatural assemblies involved in disease.

Microsecond melting of a folding intermediate in a coiled-coil peptide, monitored by T-jump/UV Raman spectroscopy

A truncated version of the GCN4 coiled-coil peptide has been studied by ultraviolet resonance Raman (UVRR) spectroscopy with 197 nm excitation, where amide modes are optimally enhanced. Although the CD melting curve could be satisfactorily described with a two-state transition having Tm = 30 degrees C, singular value decomposition of the UVRR data yielded three principal components, whose temperature dependence implicates an intermediate form between the folded and unfolded forms, with formation and melting temperatures of 10 and 40 degrees C. Two alpha-helical amide III bands, at 1340 and 1300 cm(-1), melted out selectively at 10 and 40 degrees C, respectively, and are assigned to hydrated and unhydrated helical regions. The hydrated regions are proposed to be melted in the intermediate form, while the unhydrated regions are intact. Time-resolved UVRR spectra following laser-induced temperature jumps revealed two relaxations, with time constants of 0.2 and 15 mus. These are suggested to reflect the melting times of hydrated and unhydrated helices. The unhydrated helical region may be associated with a 14-residue "trigger" sequence that has been identified in the C-terminal half of GCN4. Dehydration of helices may be a key step in the folding of coiled-coils.

Improved Stability of the Jun-Fos Activator Protein-1 Coiled Coil Motif

Two c-Jun leucine zipper variants that bind with high affinity to c-Fos have been selected using semirational design combined with protein-fragment complementation assays (JunW) or phage display selection (JunW(Ph1)). Enriched winners differ from each other in only two of ten semi-randomized positions, with DeltaT(m) values of 28 degrees C and 37 degrees C over wild-type. cFos-JunW, cFos-JunW(Ph1), and two intermediate mutants (cFos-JunW(Q21R) and cFos-JunW(E23K)) display biphasic kinetics in the folding direction, indicating the existence of a folding intermediate. The first reaction phase is fast and concentration-dependent, showing that the intermediate is readily populated and dimeric. The second phase is independent of concentration and is exponential. In contrast, in the unfolding direction, all molecules display two-state kinetics. Collectively this implies a transition state between unfolded helices and dimeric intermediate that is readily traversed in both directions. We demonstrate that the added stability of cFos-JunW(Ph1) relative to cFos-JunW is achieved via a combination of kinetic rate changes; cFos-JunW(E23K) has an increased initial dimerization rate, prior to the major transition state barrier while cFos-JunW(Q21R) displays a decreased unfolding rate. The former implies that improved hydrophobic burial and helix-stabilizing mutations exert their effect on the initial, rapid, monomer-collision event. In contrast, electrostatic interactions exert their effect late in the folding pathway. Although our focus is the leucine zipper region of the oncogenic transcription factor Activator Protein-1, coiled coils are ubiquitous and highly specific in their recognition of partners. Consequently, generating kinetic-based rules to predict and engineer their stability is of major significance in peptide-based drug design and nano-biotechnology.

Computational Study of a Single Surface-Immobilized Two-Stranded Coiled-Coil Polypeptide

The equilibrium structure and dynamics of a two-stranded coiled-coil polypeptide are investigated via Langevin dynamics simulations. An off-lattice model of the polypeptide chain is employed, which gives rise to a well-defined helical dimer native state and two-state folding kinetics. The behavior of the freely diffusing and surface-immobilized polypeptide is studied under different surface and denaturation conditions. The effect of surface immobilization on the distributions of structural and dynamical properties is considered in detail. The relationship between the simulation results and recent single-molecule fluorescence resonance energy transfer experiments performed on the two-stranded coiled-coil from the yeast transcription factor GCN4 (Jia et al. Chem. Phys. 1999, 247, 69; Talaga et al. Proc. Natl. Acad. Sci. U.S.A. 2000, 97, 13021) is discussed.

Site-specific experiments on folding/unfolding of Jun coiled coils: thermodynamic and kinetic parameters from spin inversion transfer nuclear magnetic resonance at leucine-18

The 32-residue leucine zipper subsequence, called here Jun-lz, associates in benign media to form a parallel two-stranded coiled coil. Studies are reported of its thermal unfolding/folding transition by circular dichroism (CD) on samples of natural isotopic abundance and by both equilibrium and spin inversion transfer (SIT) nuclear magnetic resonance (NMR) on samples labeled at the leucine-18 alpha-carbon with 99% 13C. The data cover a wide range of temperature and concentration, and show that Jun-lz unfolds below room temperature, being far less stable than some other leucine zippers such as GCN4. 13C-NMR shows two well-separated resonances. We ascribe the upfield one to 13C spins on unfolded single chains and the downfield one to 13C spins on coiled-coil dimers. Their relative intensities provide a measure of the unfolding equilibrium constant. In SIT NMR, the recovery of the equilibrium magnetization after one resonance is inverted is modulated in part by the unfolding and folding rate constants, which are accessible from the data. Global Bayesian analysis of the equilibrium and SIT NMR data provide values for the standard enthalpy, entropy, and heat capacity of unfolding, and show the latter to be unusually large. The CD results are compatible with the NMR findings. Global Bayesian analysis of the SIT NMR data yields the corresponding activation parameters for unfolding and folding. The results show that both reaction directions are activated processes. Activation for unfolding is entropy driven, enthalpy opposed. Activation for folding is strongly enthalpy opposed and somewhat entropy opposed, falsifying the idea that the barrier for folding is solely due to a purely entropic search for properly registered partners. The activation heat capacity is much larger for folding, so almost the entire overall change is due to the folding direction. This latter finding, if it applies to GCN4 leucine zippers, clears up an extant apparent disagreement between folding rate constants for GCN4 as determined by chevron analysis and NMR in differing temperature regimes.

Stabilization of bzip peptides through incorporation of fluorinated aliphatic residues

Two fluorinated amino acids, 5,5,5-trifluoroisoleucine (5TFI) and (2S,3R)-4,4,4-trifluorovaline (4TFV), which have been shown to serve as isoleucine surrogates in protein synthesis in Escherichia coli, have been incorporated in vivo into basic leucine zipper (bzip) peptides derived from GCN4. The extents of residue-specific incorporation of 5TFI and 4TFV were 90 and 88 %, respectively, of the encoded isoleucine residues, as evidenced by MALDI mass spectrometry and amino acid analysis. Both circular dichroism and equilibrium sedimentation studies of the fluorinated bzip peptides indicated preservation of secondary and higher-order protein structure. Thermal-denaturation experiments showed an increase of 27 degrees C in melting temperature when isoleucine was replaced by 5TFI. However, the T(m) of the peptide containing 4TFV was increased by only 4 degrees C over that of the peptide containing valine. Similar trends were observed in chemical denaturation studies in which DeltaDeltaG(unfold) in water was determined to be 2.1 or 0.3 kcal mol(-1) upon incorporation of 5TFI or 4TFV, respectively. When the fluorinated peptides were tested for DNA binding, both their affinity and specificity were similar to those of the respective hydrogenated peptides. These results suggest that fluorinated amino acids, even when introduced into the same positions, can have markedly different effects on the physical properties of proteins, while having little impact on secondary and higher-order structure.

Manipulation of hydrogel assembly and growth factor delivery via the use of peptide-polysaccharide interactions

The design of materials in which assembly, mechanical response, and biological properties are controlled by protein-polysaccharide interactions could provide materials that mimic the biological environment and find use in the delivery of growth factors. In the investigations reported here, a heparin-binding, coiled-coil peptide, PF4 ZIP, was employed to mediate the assembly of heparinized polymers. The heparin-binding affinity of this peptide was compared with that of other heparin-binding peptides (HBP) via heparin-sepharose chromatography and surface plasmon resonance (SPR) experiments. Results from these experiments indicate that PF4 ZIP demonstrates a higher heparin-binding affinity and heparin association rate when compared to the heparin-binding domains of antithrombin III (ATIII) and heparin-interacting protein (HIP). Viscoelastic hydrogels were formed upon the association of PF4 ZIP-functionalized star poly(ethylene glycol) (PEG-PF4 ZIP) with low-molecular-weight heparin-functionalized star PEG (PEG-LMWH). The viscoelastic properties of the hydrogels can be altered via variations in the ratio of LMWH to PF4 ZIP. bFGF release from these gels have also been investigated. Comparison of the bFGF release profiles with the hydrogel erosion profiles indicates that bFGF delivery from this class of hydrogels is mainly an erosion-controlled process and the rates of bFGF release can be modulated via choice of HBP or via variations in the mechanical properties of the hydrogels. Manipulation of hydrogel physical properties and erosion profiles will provide novel materials for controlled growth factor delivery and other biomedical applications.

Phi-analysis of the folding of the B domain of protein A using multiple optical probes

We examined the co-operativity of ultra-fast folding of a protein and whether the Phi-value analysis of its transition state depended on the location of the optical probe. We incorporated in turn a tryptophan residue into each of the three helices of the B domain of Protein A. Each Trp mutant of the three-helix bundle protein was used as a pseudo-wild-type parent for Phi-analysis in which the intrinsic Trp fluorescence probed the formation of each helix during the transition state. Apart from local effects in the immediate vicinity of the probe, the three separate sets of Phi-values were in excellent agreement, demonstrating the overall co-operativity of folding and the robustness of the Phi-analysis. The transition state of folding of Protein A contains the second helix being well formed with many stabilizing tertiary hydrophobic interactions. In contrast, the first and the third helices are more poorly structured in the transition state. The mechanism of folding thus involves the concurrent formation of secondary and tertiary interactions, and is towards the nucleation-condensation extreme in the nucleation-condensation-framework continuum of mechanism, with helix 2 being the nucleus. We provide an error analysis of Phi-values derived purely from the kinetics of two-state chevron plots.

T-jump infrared study of the folding mechanism of coiled-coil GCN4-p1

Partially folded intermediates have been frequently observed in equilibrium and kinetic protein folding studies. However, folding intermediates that exist at the native side of the rate-limiting step are rather difficult to study because they often evade detection by conventional folding kinetic methods. Here, we demonstrated that a laser-induced temperature-jump method can potentially be used to identify the existence of such post-transition or hidden intermediates. Specifically, we studied two cross-linked variants of GCN4-p1 coiled-coil. The GCN4 leucine zipper has been studied extensively and most of these studies have regarded it as a two-state folder. Our static circular dichroism and infrared data also indicate that the thermal unfolding of these two monomeric coiled-coils can be adequately described by an apparent two-state model. However, their temperature-jump-induced relaxation kinetics exhibit non-monoexponential behavior, dependent upon sequence and temperature. Taken together, our results support a folding mechanism wherein at least one folding intermediate populates behind the main rate-limiting step.

The Structure of the Major Transition State for Folding of an FF Domain from Experiment and Simulation

We have analysed the transition state of folding of the four-helix FF domain from HYPA/FBP11 by high-resolution experiment and simulation as part of a continuing effort to understand the principles of folding and the refinement of predictive methods. The major transition state for folding was subjected to a Phi-value analysis utilising 50 mutants. The transition state contained a nucleus for folding centred around the end of helix 1 (H1) and the beginning of helix 2 (H2). Secondary structure in this region was fully formed (PhiF=0.9-1) and tertiary interactions were well developed. Interactions in the distal part of the native structure were weak (PhiF=0-0.2). The hydrophobic core and other parts of the protein displayed intermediate Phi-values, suggesting that interactions coalesce as the end of H1 and beginning of H2 are in the process of being formed. The distribution of Phi-values resembled that of barnase, which folds via an intermediate, rather than that of CI2 which folds by a concerted nucleation-condensation mechanism. The overall picture of the transition state structure identified in molecular dynamics simulations is in qualitative agreement, with the turn connecting H1 and H2 being formed, a loosened core, and H4 partially unfolded and detached from the core. There are some differences in the details and interpretation of specific Phi-values.

The folding of spectrin domains II: phi-value analysis of R16

Studies on the folding of helical proteins have shown a wide range of different mechanisms and highlighted the importance of helical propensity as a factor in determining folding mechanism. Here, we contribute to this interesting field with the protein engineering phi-value analysis of the 16th domain of chicken brain alpha-spectrin, R16. The fortuitous curvature seen in the unfolding arm of the chevron plot allows us to investigate both early and late events in folding. R16 is the first two-state helical protein for which this has been possible.

Differences in the folding transition state of ubiquitin indicated by phi and psi analyses

We compare the folding transition state (TS) of ubiquitin previously identified by using psi analysis to that determined by using analysis. Both methods attempt to identify interactions and their relative populations at the rate-limiting step for folding. The TS ensemble derived from psi analysis has an extensive native-like chain topology, with a four-stranded beta-sheet network and a portion of the major helix. According to analysis, however, the TS is much smaller and more polarized, with only a local helix/hairpin motif. We find that structured regions can have values far from unity, the canonical value for such sites, because of structural relaxation of the TS. Consequently, these sites may be incorrectly interpreted as contributing little to the structure of the TS. These results stress the need for caution when interpreting and drawing conclusions from analysis alone and highlight the need for more specific tools for examining the structure and energetics of the TS ensemble.

Fast folding of a helical protein initiated by the collision of unstructured chains

To examine whether helix formation necessarily precedes chain collision, we have measured the folding of a fully helical coiled coil that has been specially engineered to have negligible intrinsic helical propensity but high overall stability. The folding rate approaches the diffusion-limited value and is much faster than possible if folding is contingent on precollision helix formation. Therefore, the collision of two unstructured chains is the initial step of the dominant kinetic pathway, whereas helicity exerts its influence only at a later step. Folding from an unstructured encounter complex may be efficient and robust, which has implications for any biological process that couples folding to binding.

Helix formation and the unfolded state of a 52-residue helical protein

A growing class of proteins in biological processes has been found to be unfolded on isolation under normal solution conditions. We have used NMR spectroscopy to characterize the structural and dynamic properties of the unfolded and partially folded states of a 52-residue alanine-rich protein (Ala-14) at temperatures from -5 degrees C to 40 degrees C. At 40 degrees C, alanine residues in Ala-14 adopt phi and psi angles, consistent with a significant ensemble population of polyproline II conformation. Analysis of relaxation rates in the protein reveals that a series of residues, Gln 35-Ala 36-Ala 37-Lys 38-Asp 39-Asp 40-Ala 41-Ala 42, displays slow motional dynamics at both -5 degrees C and 40 degrees C. Temperature-dependent chemical shift changes indicate that this region is the site of helix initiation. The remaining N-terminal residues become increasingly dynamic as they extend from the nucleation site. The C terminus remains dynamic and changes less with temperature, indicating it is relatively unstructured. Ala-14 provides a high-resolution portrait of the unfolded state and the process of helix nucleation and propagation in the absence of tertiary contacts, information that bears on early events in protein folding.

Clustering of large hydrophobes in the hydrophobic core of two-stranded alpha-helical coiled-coils controls protein folding and stability

The de novo design and biophysical characterization of two 60-residue peptides that dimerize to fold as parallel coiled-coils with different hydrophobic core clustering is described. Our goal was to investigate whether designing coiled-coils with identical hydrophobicity but with different hydrophobic clustering of non-polar core residues (each contained 6 Leu, 3 Ile, and 7 Ala residues in the hydrophobic core) would affect helical content and protein stability. The disulfide-bridged P3 and P2 differed dramatically in alpha-helical structure in benign conditions. P3 with three hydrophobic clusters was 98% alpha-helical, whereas P2 was only 65% alpha-helical. The stability profiles of these two analogs were compared, and the enthalpy and heat capacity changes upon denaturation were determined by measuring the temperature dependence by circular dichroism spectroscopy and confirmed by differential scanning calorimetry. The results showed that P3 assembled into a stable alpha-helical two-stranded coiled-coil and exhibited a native protein-like cooperative two-state transition in thermal melting, chemical denaturation, and calorimetry experiments. Although both peptides have identical inherent hydrophobicity (the hydrophobic burial of identical non-polar residues in equivalent heptad coiled-coil positions), we found that the context dependence of an additional hydrophobic cluster dramatically increased stability of P3 (Delta Tm approximately equal to 18 degrees C and Delta[urea](1/2) approximately equal to 1.5 M) as compared with P2. These results suggested that hydrophobic clustering significantly stabilized the coiled-coil structure and may explain how long fibrous proteins like tropomyosin maintain chain integrity while accommodating polar or charged residues in regions of the protein hydrophobic core.

Hg(II) binding to a weakly associated coiled coil nucleates an encoded metalloprotein fold: a kinetic analysis

A detailed kinetic analysis of metal encapsulation by a de novo-designed protein is described. The kinetic mechanism of Hg(II) encapsulation in the three-stranded coiled coil formed by the peptide CH(3)CO-G LKALEEK CKALEEK LKALEEK G-NH(2) (Baby L9C) is derived by global analysis. The mechanism involves rapid initial collapse of two peptides by Hg(II) forming Hg(Baby L9C(-H))(2) with a linear thiolato Hg(II) bound to the cysteine sulfur atoms. Here, Baby L9C(-H) denotes Baby L9C with the cysteine thiol deprotonated. Addition of the third peptide, forming the three-stranded coiled coil, is the rate-determining step and results in an intermediate state involving two separate species. One of the species, termed the properly folded intermediate, undergoes rapid deprotonation of the third cysteine thiol, yielding the desired three-stranded coiled coil with an encapsulated trigonal thiolato Hg(II). The other species, termed the misfolded intermediate, rearranges in an experimentally distinguishable step to the properly folded intermediate. The order of the reaction involving the addition of the third peptide with respect to the concentration of Baby L9C indicates that addition of the third helix only occurs through reaction of Hg(Baby L9C(-H))(2) and Baby L9C that is unassociated with a coiled coil. Temperature dependence of the reaction afforded activation parameters for both the addition of the third helix (deltaH = 20(2) kcalmol; deltaS= 40(5) calmol K) and the rearrangement of the misfolded intermediate steps (deltaH = 23(2) kcalmol; deltaS= 27(5) calmol K). The mechanism is discussed with regard to metalloprotein folding and metalloprotein design.

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.

Unfolding of a leucine zipper is not a simple two-state transition

Temperature-induced unfolding of the leucine zipper, an alpha-helical, double-stranded, coiled-coil, was studied by circular dichroism spectroscopy, spectrofluorimetry and heat capacity scanning calorimetry. It is shown that this process does not represent a simple two-state transition, as previously believed, but consists of several stages. The first transition starts at the very beginning of heating from 0 degrees C and proceeds with significant heat absorption and decrease of ellipticity. This transition does not depend on the concentration of protein and is sensitive to modification of the N terminus; it is therefore associated with unfolding or fraying of this part of the leucine zipper. The second transition takes place at a considerably higher temperature; it is more pronounced than the first one and does not depend on the concentration of protein, i.e. it is unimolecular. This transition is sensitive to modification of both termini of the leucine zipper and affects the optical properties of a tryptophan residue placed in the central part of the zipper. It therefore involves the whole dimer but does not result in its dissociation, presumably being associated with some repacking of the coiled-coil. This second transition is followed at higher temperatures by the concentration-dependent cooperative unfolding/dissociation of the two strands. The enthalpy and entropy of the temperature-induced structural changes of the leucine zipper that take place before its cooperative unfolding/dissociation comprises almost 40% of the total enthalpy and entropy of unfolding of the completely folded coiled-coil, the state in which it appears to be below 0 degrees C. Comparison of the total enthalpy of leucine zipper unfolding with that of a single-stranded alpha-helix shows that their temperature-dependence correlates with the extent of intramolecular non-polar contacts and allows an assessment of the enthalpy of hydrogen bonding in alpha-helices, which appears to be about 3.3kJmol(-1) at 20 degrees C.

pH-induced folding of an apoptotic coiled coil

Par-4 is a 38-kD protein pivotal to the apoptotic pathways of various cell types, most notably prostate cells and neurons, where it has been linked to prostate cancer and various neurodegenerative disorders including Alzheimer's and Huntington's diseases and HIV encephalitis. The C-terminal region of Par-4 is responsible for homodimerization and the ability of Par-4 to interact with proposed effector molecules. In this study, we show that the C-terminal 47 residues of Par-4 are natively unfolded at physiological pH and temperature. Evidence is rapidly accumulating that natively unfolded proteins play an important role in various cellular functions and signaling pathways, and that folding can often be induced on complexation with effector molecules or alteration of environment. Here we use primarily CD studies to show that changes in the environment, particularly pH and temperature, can induce the Par-4 C terminus to form a self-associated coiled coil.

Dissociation and unfolding of GCN4 leucine zipper in the presence of sodium dodecyl sulfate

The dissociation and unfolding behavior of the GCN4 leucine zipper has been studied using SDS titration. Circular dichroism (CD) spectra showed that the alpha-helix content of the leucine zipper (20 microM) decreased during the sodium dodecyl sulfate (SDS) titration. However, the alpha-helix content of the leucine zipper still remained significant in the presence of 1 mM SDS, with little change detected when the SDS concentration further increased to 2 mM. The dimer dissociation of the leucine zipper is also a co-operative process during SDS titration; with no dimer remaining when SDS concentration reached 1 mM, as shown by electrophoresis and the the theta(222)/theta(208) ratio. Our results indicate that SDS efficiently induces leucine zipper dimer dissociation with the monomers still partially folded. The experimental results provide important evidence for the previous model that partial helix formation precedes dimerization in coiled coil folding.

Unfolding of apomyoglobin examined by synchrotron footprinting

A new method to examine the structure and stability of proteins using footprinting is applied to examine the unfolding of apomyoglobin. Unlike previous cleavage based footprinting methods, synchrotron X-ray protein footprinting is based on a quantitative determination of the extent and the site of radiolytic modification of amino acid side chains, analyzed using mass spectrometry. The amino acids most susceptible to radiolytic oxidation (cysteine, methionine, phenylalanine, tyrosine, tryptophan, histidine, proline, and leucine) serve as convenient probes of protein structure to monitor changes in solvent accessibility. To determine if the technique can measure quantitative properties of proteins relevant to structure and function, we examined the equilibrium unfolding of apomyoglobin in urea and compared the results to data derived from fluorescence studies under the same conditions.

Helix-stabilized Fv (hsFv) antibody fragments: substituting the constant domains of a Fab fragment for a heterodimeric coiled-coil domain

Antibody Fv fragments would in principle be useful for a variety of biotechnological applications because of their small size and the possibility to produce them in relatively large amounts in recombinant form; however, their limited stability is a drawback. To solve this problem, both domains are usually fused via a peptide linker to form a single-chain Fv (scFv) fragment, but in some cases this leads to a dimerization. We present an alternative format for stabilizing antibody Fv fragments. The C(H)1 and C(L) domain of the Fab fragment were replaced with a heterodimeric coiled coil (WinZip-A2B1), which had previously been selected using a protein-fragment complementation assay in Escherichia coli. This new antibody format was termed helix-stabilized Fv fragment (hsFv), and was compared to the corresponding Fv, Fab and single-chain Fv format. Bacterial growth and expression of the hsFv was significantly improved compared to the Fab fragment. The hsFv fragment formed a heterodimer of heavy and light chain with the expected molecular mass, also under conditions where the scFv fragment was predominantly dimeric. The hsFv fragment was significantly more stable than the Fv fragment, and nearly as stable as the scFv fragment under the conditions used (80 nM protein concentration). Thus, the format of a helix-stabilized Fv (hsFv) fragment can be a useful alternative to existing recombinant antibody formats, especially in cases where poor expression of Fab fragments or multimerization of scFv fragments is a problem.

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.

Molecular characterization of calmodulin trapping by calcium/calmodulin-dependent protein kinase II

Autophosphorylation of alpha-Ca(2+)/calmodulin-dependent protein kinase II (CaM kinase II) at Thr(286) results in calmodulin (CaM) trapping, a >10,000-fold decrease in the dissociation rate of CaM from the enzyme. Here we present the first site-directed mutagenesis study on the dissociation of the high affinity complex between CaM and full-length CaM kinase II. We measured dissociation kinetics of CaM and CaM kinase II proteins by using a fluorescently modified CaM that is sensitive to binding to target proteins. In low [Ca(2+)], the phosphorylated mutant kinase F293A and the CaM mutant E120A/M124A exhibited deficient trapping compared with wild-type. In high [Ca(2+)], the CaM mutations E120A, M124A, and E120A/M124A and the CaM kinase II mutations F293A, F293E, N294A, N294P, and R297E increased dissociation rate constants by factors ranging from 2.3 to 116. We have also identified residues in CaM and CaM kinase II that interact in the trapped state by mutant cycle-based analysis, which suggests that interactions between Phe(293) in the kinase and Glu(120) and Met(124) in CaM specifically stabilize the trapped CaM-CaM kinase II complex. Our studies further show that Phe(293) and Asn(294) in CaM kinase II play dual roles, because they likely destabilize the low affinity state of CaM complexed to unphosphorylated kinase but stabilize the trapped state of CaM bound to phosphorylated kinase.

An intrahelical salt bridge within the trigger site stabilizes the GCN4 leucine zipper

We previously reported that a helical trigger segment within the GCN4 leucine zipper monomer is indispensable for the formation of its parallel two-stranded coiled coil. Here, we demonstrate that the intrinsic secondary structure of the trigger site is largely stabilized by an intrahelical salt bridge. Removal of this surface salt bridge by a single amino acid mutation induced only minor changes in the backbone structure of the GCN4 leucine zipper dimer as verified by nuclear magnetic resonance. The mutation, however, substantially destabilized the dimeric structure. These findings support the proposed hierarchic folding mechanism of the GCN4 coiled coil in which local helix formation within the trigger segment precedes dimerization.

Temperature dependence of the folding and unfolding kinetics of the GCN4 leucine zipper via 13C(alpha)-NMR

Studies by one-dimensional NMR are reported on the interconversion of folded and unfolded forms of the GCN4 leucine zipper in neutral saline buffer. The peptide bears 99% 13C(alpha) labels at three sites: V9, L12, and G31. Time-domain 13C(alpha)-NMR spectra are interpreted by global Bayesian lineshape analysis to extract the rate constants for both unfolding and folding as functions of temperature in the range 47-71 degrees C. The data are well fit by the assumption that the same rate constants apply at each labeled site, confirming that only two conformational states need be considered. Results show that 1) both processes require a free energy of activation; 2) unfolding is kinetically enthalpy-opposed and entropy-driven, while folding is the opposite; and 3) the transition state dimer ensemble averages approximately 40% helical. The activation parameters for unfolding, derived from NMR data at the elevated temperatures where both conformations are populated, lead to estimates of the rate constant at low temperatures (5-15 degrees C) that agree with extant values determined by stopped-flow CD via dilution from denaturing media. However, the corresponding estimated values for the folding rate constant are larger by two to three orders of magnitude than those obtained by stopped flow. We propose that this apparent disagreement is caused by the necessity, in the stopped-flow experiment, for initiation of new helices as the highly denaturant-unfolded molecule adjusts to the newly created benign solvent conditions. This must reduce the success rate of collisions in producing the folded molecule. In the NMR determinations, however, the unfolded chains always have a small, but essential, helix content that makes such initiation unnecessary. Support for this hypothesis is adduced from recent extant experiments on the helix-coil transition in single-chain helical peptides and from demonstration that the folding rate constants for coiled coils, as obtained by stopped flow, are influenced by the nature of the denaturant used.

Dynamics and folding of single two-stranded coiled-coil peptides studied by fluorescent energy transfer confocal microscopy

We report single-molecule measurements on the folding and unfolding conformational equilibrium distributions and dynamics of a disulfide crosslinked version of the two-stranded coiled coil from GCN4. The peptide has a fluorescent donor and acceptor at the N termini of its two chains and a Cys disulfide near its C terminus. Thus, folding brings the two N termini of the two chains close together, resulting in an enhancement of fluorescent resonant energy transfer. End-to-end distance distributions have thus been characterized under conditions where the peptide is nearly fully folded (0 M urea), unfolded (7.4 M urea), and in dynamic exchange between folded and unfolded states (3.0 M urea). The distributions have been compared for the peptide freely diffusing in solution and deposited onto aminopropyl silanized glass. As the urea concentration is increased, the mean end-to-end distance shifts to longer distances both in free solution and on the modified surface. The widths of these distributions indicate that the molecules are undergoing millisecond conformational fluctuations. Under all three conditions, these fluctuations gave nonexponential correlations on 1- to 100-ms time scale. A component of the correlation decay that was sensitive to the concentration of urea corresponded to that measured by bulk relaxation kinetics. The trajectories provided effective intramolecular diffusion coefficients as a function of the end-to-end distances for the folded and unfolded states. Single-molecule folding studies provide information concerning the distributions of conformational states in the folded, unfolded, and dynamically interconverting states.

Folding of a three-stranded coiled coil

Coiled coils consist of two or more amphipathic a-helices wrapped around each other to form a superhelical structure stabilized at the interhelical interface by hydrophobic residues spaced in a repeating 3-4 sequence pattern. Dimeric coiled coils have been shown to often form in a single step reaction in which association and folding of peptide chains are tightly coupled. Here, we ask whether such a simple folding mechanism may also apply to the formation of a three-stranded coiled coil. The designed 29-residue peptide LZ16A was shown previously to be in a concentration-dependent equilibrium between unfolded monomer (M), folded dimer (D), and folded trimer (T). We show by time-resolved fluorescence change experiments that folding of LZ16A to D and T can be described by 2M (k1)<==>(k(-1)) D and M + D (k2)<==>(k(-2)) T. The following rate constants were determined (25 degrees C, pH 7): k1 = 7.8 x 10(4) M(-1) s(-1), k(-1) = 0.015 s(-1), k2 = 6.5 x 10(5) M(-1) s(-1), and k(-2) = 1.1 s(-1). In a separate experiment, equilibrium binding constants were determined from the change with concentration of the far-ultraviolet circular dichroism spectrum of LZ16A and were in good agreement with the kinetic rate constants according to K(D) = k1/2k(-1) and K(T) = k2/k(-2). Furthermore, pulsed hydrogen-exchange experiments indicated that only unfolded M and folded D and T were significantly populated during folding. The results are compatible with a two-step reaction in which a subpopulation of association competent (e.g., partly helical) monomers associate to dimeric and trimeric coiled coils.

A breakdown of symmetry in the folding transition state of protein L

The 62 residue IgG binding domain of protein L consists of a central alpha-helix packed on a four-stranded beta-sheet formed by N and C-terminal beta-hairpins. The overall topology of the protein is quite symmetric: the beta-hairpins have similar lengths and make very similar interactions with the central helix. Characterization of the effects of 70 point mutations distributed throughout the protein on the kinetics of folding and unfolding reveals that this symmetry is completely broken during folding; the first beta-hairpin is largely structured while the second beta-hairpin and helix are largely disrupted in the folding transition state ensemble. The results are not consistent with a "hydrophobic core first" picture of protein folding; the first beta-hairpin appears to be at least as ordered at the rate limiting step in folding as the hydrophobic core.