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

Sequential Translocation of Polypeptides across the Bacterial Outer Membrane through the Trimeric Autotransporter Pathway

Trimeric autotransporter adhesins (TAAs) are a family of bacterial outer membrane (OM) proteins that are comprised of three identical subunits. Each subunit contains an N-terminal extracellular ("passenger") domain and a short C-terminal segment that contributes four β strands to a single 12-stranded β barrel. The mechanism by which the passenger domains are translocated across the OM and the energetics of the translocation reaction are poorly understood. To address these issues, we examined the secretion of modified versions of the passenger domain of UpaG, a TAA produced by Escherichia coli CFT073. Using the SpyTag-SpyCatcher system to probe passenger domain localization, we found that both intrinsically disordered polypeptides fused to the UpaG passenger domain and artificially disulfide-bonded polypeptides were secreted effectively but relatively slowly. Surprisingly, we also found that in some cases, the three nonnative passenger domain segments associated with a single trimer were secreted sequentially. Photo-cross-linking experiments indicated that incompletely assembled UpaG derivatives remained bound to the barrel assembly machinery (Bam) complex until all three passenger domains were fully secreted. Taken together, our results strongly suggest that the secretion of polypeptides through the TAA pathway is coordinated with the assembly of the β barrel domain and that the folding of passenger domains in the extracellular space maximizes the rate of secretion. Furthermore, our work provides evidence for an unprecedented sequential mode of protein translocation, at least under specific experimental conditions.IMPORTANCE Trimeric autotransporter adhesins (TAAs) are specialized bacterial outer membrane proteins consisting of three identical subunits. TAAs contain large extracellular domains that trimerize and promote virulence, but the mechanism by which they are secreted is poorly understood. We found that the extracellular domains of a native TAA were secreted rapidly but that disordered and artificially folded polypeptides fused to native passenger domains were secreted in a slow, sequential fashion. Our results strongly suggest that the efficient secretion of native extracellular domains is driven by their trimerization following export but that alternative energy sources can be harnessed to secrete nonnative polypeptides. Furthermore, we obtained evidence that TAA extracellular domains are secreted before the assembly of the linked membrane spanning domain is completed.

Exposing the Nucleation Site in α-Helix Folding: A Joint Experimental and Simulation Study

One of the fundamental events in protein folding is α-helix formation, which involves sequential development of a series of helical hydrogen bonds between the backbone C═O group of residues i and the -NH group of residues i + 4. While we now know a great deal about α-helix folding dynamics, a key question that remains to be answered is where the productive helical nucleation event occurs. Statistically, a helical nucleus (or the first helical hydrogen-bond) can form anywhere within the peptide sequence in question; however, the one that leads to productive folding may only form at a preferred location. This consideration is based on the fact that the α-helical structure is inherently asymmetric, due to the specific alignment of the helical hydrogen bonds. While this hypothesis is plausible, validating it is challenging because there is not an experimental observable that can be used to directly pinpoint the location of the productive nucleation process. Therefore, in this study we combine several techniques, including peptide cross-linking, laser-induced temperature-jump infrared spectroscopy, and molecular dynamics simulations, to tackle this challenge. Taken together, our experimental and simulation results support an α-helix folding mechanism wherein the productive nucleus is formed at the N-terminus, which propagates toward the C-terminal end of the peptide to yield the folded structure. In addition, our results show that incorporation of a cross-linker can lead to formation of differently folded conformations, underscoring the need for all-atom simulations to quantitatively assess the proposed cross-linking design.

Unified understanding of folding and binding mechanisms of globular and intrinsically disordered proteins

Extensive experimental and theoretical studies have advanced our understanding of the mechanisms of folding and binding of globular proteins, and coupled folding and binding of intrinsically disordered proteins (IDPs). The forces responsible for conformational changes and binding are common in both proteins; however, these mechanisms have been separately discussed. Here, we attempt to integrate the mechanisms of coupled folding and binding of IDPs, folding of small and multi-subdomain proteins, folding of multimeric proteins, and ligand binding of globular proteins in terms of conformational selection and induced-fit mechanisms as well as the nucleation–condensation mechanism that is intermediate between them. Accumulating evidence has shown that both the rate of conformational change and apparent rate of binding between interacting elements can determine reaction mechanisms. Coupled folding and binding of IDPs occurs mainly by induced-fit because of the slow folding in the free form, while ligand binding of globular proteins occurs mainly by conformational selection because of rapid conformational change. Protein folding can be regarded as the binding of intramolecular segments accompanied by secondary structure formation. Multi-subdomain proteins fold mainly by the induced-fit (hydrophobic collapse) mechanism, as the connection of interacting segments enhances the binding (compaction) rate. Fewer hydrophobic residues in small proteins reduce the intramolecular binding rate, resulting in the nucleation–condensation mechanism. Thus, the folding and binding of globular proteins and IDPs obey the same general principle, suggesting that the coarse-grained, statistical mechanical model of protein folding is promising for a unified theoretical description of all mechanisms.

Flexibility vs Preorganization: Direct Comparison of Binding Kinetics for a Disordered Peptide and Its Exact Preorganized Analogues

Many intrinsically disordered proteins, which are prevalent in nature, fold only upon binding their structured partner proteins. Such proteins have been hypothesized to have a kinetic advantage over their folded, preorganized analogues in binding their partner proteins. Here we determined the effects of ligand preorganization on the k_on for a biomedically important system: an intrinsically disordered p53 peptide ligand and the MDM2 protein receptor. Based on direct simulations of binding pathways, computed k_on values for fully disordered and preorganized p53 peptide analogues were within error of each other, indicating little if any kinetic advantage to being disordered or preorganized for binding the MDM2 protein. We also examined the effects of increasing the concentration of MDM2 on the extent to which its mechanism of binding to the p53 peptide is induced fit vs conformational selection. Results predict that the mechanism is solely induced fit if the unfolded state of the peptide is more stable than its folded state; otherwise, the mechanism shifts from being dominated by conformational selection at low MDM2 concentration to induced fit at high MDM2 concentration. Taken together, our results are relevant to any protein binding process that involves a disordered peptide of a similar length that forms a single α-helix upon binding a partner protein. Such disorder-to-helix transitions are common among protein interactions of disordered proteins and are therefore of fundamental biological interest.

Microscopic nucleation and propagation rates of an alanine-based α-helix

An infrared temperature-jump (T-jump) study by Huang et al. (Proc. Natl. Acad. Sci. U. S. A., 2002, 99, 2788-2793) showed that the conformational relaxation kinetics of an alanine-based α-helical peptide depend not only on the final temperature (Tf) but also on the initial temperature (Ti) when Tf is fixed. Their finding indicates that the folding free energy landscape of this peptide is non-two-state like, allowing for the population of conformational ensembles with different helical lengths and relaxation times in the temperature range of the experiment. Because α-helix folding involves two fundamental events, nucleation and propagation, the results of Huang et al. thus present a unique opportunity to determine their rate constants - a long-sought goal in the study of the helix-coil transition dynamics. Herein, we capitalize on this notion and develop a coarse-grained kinetic model to globally fit the thermal unfolding curve and T-jump kinetic traces of this peptide. Using this strategy, we are able to explicitly determine the microscopic rate constants of the kinetic steps encountered in the nucleation and propagation processes. Our results reveal that the time taken to form one α-helical turn (i.e., an α-helical segment with one helical hydrogen bond) is about 315 ns, whereas the time taken to elongate this nucleus by one residue (or backbone unit) is 5.9 ns, depending on the position of the residue.

Effect of helix length on the stability of the Lac repressor antiparallel coiled coil

The helix length dependence of the stability of antiparallel four-chain coiled coils is investigated using eight synthetic peptides (Lac21-Lac28) whose sequences are derived from the tetramerization domain of the Lac repressor protein. Previous studies using analytical ultracentrifugation sedimentation equilibrium experiments to characterize Lac21 and Lac28 justifies the use of a two state model to describe the unfolding behavior of these two peptides. Using circular dichroism spectropolarimetry as a measure of tetramer assembly, both chemical and thermal denaturation experiments were carried out to determine thermodynamic parameters. We found that the hydrophobic core residues provide the greatest impact on stability and, as a consequence, must reorganize the register of the antiparallel helices to accommodate the burial of the nonpolar amino acids. Addition of noncore residues appears to have only a minor effect on stability, and in some cases, show a slight destabilization.

Dihedral angle entropy measures for intrinsically disordered proteins

Protein stability is based on a delicate balance between energetic and entropic factors. Intrinsically disordered proteins (IDPs) interacting with a folded partner protein in the act of binding can order the IDP to form the correct functional interface by decrease in the overall free energy. In this work, we evaluate the part of the entropic cost of ordering an IDP arising from their dihedral states. The IDP studied is a leucine zipper dimer that we simulate with molecular dynamics and find that it does show disorder in six phi and psi dihedral angles of the N terminal sequence of one monomer. Essential to ascertain is the degree of disorder in the IDP, and we do so by considering the entire, discretized probability distribution function of N dihedrals with M conformers per dihedral. A compositional clustering method is introduced, whereby the NS = N(M) states are formed from the Cartesian product of each dihedral's conformational space. Clustering is carried out with a version of a k-means algorithm that accounts for the circular nature of dihedral angles. For the 12 dihedrals each found to have three conformers, among the resulting 531441 states, their populations show that the first 100 (500) most populated states account for ∼65% (∼90%) of the entire population, indicating that there are strong dependencies among the dihedrals' conformations. These state populations are used to evaluate a Kullback-Leibler divergence entropy measure and obtain the dihedral configurational entropy S. At 300 K, TS ∼ 3 kcal/mol, showing that IDP entropy, while roughly half that would be expected from independently distributed dihedrals, can be a decisive contributor to the free energy of this IDP binding and ordering.

Sequence-resolved free energy profiles of stress-bearing vimentin intermediate filaments

Intermediate filaments (IFs) are key to the mechanical strength of metazoan cells. Their basic building blocks are dimeric coiled coils mediating hierarchical assembly of the full-length filaments. Here we use single-molecule force spectroscopy by optical tweezers to assess the folding and stability of coil 2B of the model IF protein vimentin. The coiled coil was unzipped from its N and C termini. When pulling from the C terminus, we observed that the coiled coil was resistant to force owing to the high stability of the C-terminal region. Pulling from the N terminus revealed that the N-terminal half is considerably less stable. The mechanical pulling assay is a unique tool to study and control seed formation and structure propagation of the coiled coil. We then used rigorous theory-based deconvolution for a model-free extraction of the energy landscape and local stability profiles. The data obtained from the two distinct pulling directions complement each other and reveal a tripartite stability of the coiled coil: a labile N-terminal half, followed by a medium stability section and a highly stable region at the far C-terminal end. The different stability regions provide important insight into the mechanics of IF assembly.

Slow folding-unfolding kinetics of an octameric β-peptide bundle

β-Peptide foldamers offer attractive frameworks for examining the effect of backbone flexibility on the dynamics of protein folding. Herein, we study the folding-unfolding kinetics of a β-peptide, Acid-1Y,1 which folds in aqueous solution into an octameric bundle of peptides in a conformation known as the 14-helix. Acid-1Y is comprised exclusively of β-amino acids, which differ from α-amino acids by the addition of a single methylene into the backbone. We aim to understand how the additional degree of freedom and increased backbone flexibility in the β-amino acid affect folding dynamics and to measure folding rates of this octameric β-peptide. Previously, we found that the T-jump induced relaxation kinetics of a monomeric β-peptide that forms a monomeric 14-helix occurred on the nanosecond time scale2 and were noticeably slower than a similar alanine-based α-helical peptide.3 Additionally, in comparison to similar α-helices, the relaxation rates showed a weaker dependence on temperature. Here, we find that the T-jump induced relaxation kinetics of the octameric β-peptide occurs on an even slower time scale (minutes) and the unfolding relaxation rates show a large dependence on temperature. These differences indicate that folding energy landscapes of β-peptide secondary and quaternary structure are markedly distinct from one another and also from their α-helical counterparts.

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.

Simplified protein models: predicting folding pathways and structure using amino acid sequences

We demonstrate the ability of simultaneously determining a protein's folding pathway and structure using a properly formulated model without prior knowledge of the native structure. Our model employs a natural coordinate system for describing proteins and a search strategy inspired by the observation that real proteins fold in a sequential fashion by incrementally stabilizing nativelike substructures or "foldons." Comparable folding pathways and structures are obtained for the twelve proteins recently studied using atomistic molecular dynamics simulations [K. Lindorff-Larsen, S. Piana, R. O. Dror, D. E. Shaw, Science 334, 517 (2011)], with our calculations running several orders of magnitude faster. We find that nativelike propensities in the unfolded state do not necessarily determine the order of structure formation, a departure from a major conclusion of the molecular dynamics study. Instead, our results support a more expansive view wherein intrinsic local structural propensities may be enhanced or overridden in the folding process by environmental context. The success of our search strategy validates it as an expedient mechanism for folding both in silico and in vivo.

The loop hypothesis: contribution of early formed specific non-local interactions to the determination of protein folding pathways

The extremely fast and efficient folding transition (in seconds) of globular proteins led to the search for some unifying principles embedded in the physics of the folding polypeptides. Most of the proposed mechanisms highlight the role of local interactions that stabilize secondary structure elements or a folding nucleus as the starting point of the folding pathways, i.e., a "bottom-up" mechanism. Non-local interactions were assumed either to stabilize the nucleus or lead to the later steps of coalescence of the secondary structure elements. An alternative mechanism was proposed, an "up-down" mechanism in which it was assumed that folding starts with the formation of very few non-local interactions which form closed long loops at the initiation of folding. The possible biological advantage of this mechanism, the "loop hypothesis", is that the hydrophobic collapse is associated with ordered compactization which reduces the chance for degradation and misfolding. In the present review the experiments, simulations and theoretical consideration that either directly or indirectly support this mechanism are summarized. It is argued that experiments monitoring the time-dependent development of the formation of specifically targeted early-formed sub-domain structural elements, either long loops or secondary structure elements, are necessary. This can be achieved by the time-resolved FRET-based "double kinetics" method in combination with mutational studies. Yet, attempts to improve the time resolution of the folding initiation should be extended down to the sub-microsecond time regime in order to design experiments that would resolve the classes of proteins which first fold by local or non-local interactions.

Folding and binding of an intrinsically disordered protein: fast, but not 'diffusion-limited'

Coupled folding and binding of intrinsically disordered proteins (IDPs) is prevalent in biology. As the first step toward understanding the mechanism of binding, it is important to know if a reaction is 'diffusion-limited' as, if this speed limit is reached, the association must proceed through an induced fit mechanism. Here, we use a model system where the 'BH3 region' of PUMA, an IDP, forms a single, contiguous α-helix upon binding the folded protein Mcl-1. Using stopped-flow techniques, we systematically compare the rate constant for association (k(+)) under a number of solvent conditions and temperatures. We show that our system is not 'diffusion-limited', despite having a k(+) in the often-quoted 'diffusion-limited' regime (10(5)-10(6) M(-1) s(-1) at high ionic strength) and displaying an inverse dependence on solvent viscosity. These standard tests, developed for folded protein-protein interactions, are not appropriate for reactions where one protein is disordered.

Mutational effect of structural parameters on coiled-coil stability of proteins

Understanding the parameters that influence the melting temperature of coiled-coils (CC) and their stability is very important. We have analyzed 45 CC mutants of DNA binding protein, electron transport protein, hydrolase, oxidoreductase, and transcription factors. Many mutants have been observed at Tm = 40 °C–60 °C with ΔS = 9–11 kcal/°C mol, ΔG = −400 to −450 kcal/mol, and Keq = 0.98–1.03. The multiple regression analysis of Tm reveals that influences of thermodynamic parameters are strong (R = 0.97); chemical parameters are moderate (R = 0.63); and the geometrical parameters are negligible (R = 0.19). The combination of all these three parameters exhibits a little higher influence on Tm (R = 0.98). From the analysis, it has been concluded that the thermodynamic parameters alone are very important in stability studies on protein coil mutants. Besides, the derived regression model would have been useful for the reliable prediction of the melting temperature of coil mutants.

Slow, reversible, coupled folding and binding of the spectrin tetramerization domain

Many intrinsically disordered proteins (IDPs) are significantly unstructured under physiological conditions. A number of these IDPs have been shown to undergo coupled folding and binding reactions whereby they can gain structure upon association with an appropriate partner protein. In general, these systems display weaker binding affinities than do systems with association between completely structured domains, with micromolar K(d) values appearing typical. One such system is the association between α- and β-spectrin, where two partially structured, incomplete domains associate to form a fully structured, three-helix bundle, the spectrin tetramerization domain. Here, we use this model system to demonstrate a method for fitting association and dissociation kinetic traces where, using typical biophysical concentrations, the association reactions are expected to be highly reversible. We elucidate the unusually slow, two-state kinetics of spectrin assembly in solution. The advantages of studying kinetics in this regime include the potential for gaining equilibrium constants as well as rate constants, and for performing experiments with low protein concentrations. We suggest that this approach would be particularly appropriate for high-throughput mutational analysis of two-state reversible binding processes.

Simulations of temperature and salt concentration effects on bZIP, a basic region leucine zipper

Basic region leucine zipper (bZIP) transcription factors are dimeric proteins that recognize DNA. The monomers consist of a leucine zipper subdomain responsible for dimerization and a highly basic DNA recognition subdomain. Twelve explicit solvent molecular dynamics (MD) trajectories were run on the GCN4 bZIP transcriptional factor in the absence of DNA at three temperatures and two ion concentrations (0 mM with Cl(-) ions to neutralize the bZIP and 200 mM with additional Na(+) and Cl(-) ions) to probe the conformational ensemble that the basic region samples. In most trajectories, the basic region exhibits an alligator-jaw-like opening and closing (only one monomer moves), versus scissor-like motion, by a mainly rigid body, hinge motion centered on three "fork" residues that span the basic region to the coiled coil. In this motion, the α-helical character of the basic region monomers is mostly maintained. A broad range of distances is accessed, consistent with the absence of particular interactions for the basic region monomers. In two of the trajectories, the basic region monomers "collapse" to form a stable state. The coiled coil, leucine zipper subdomain is very stable for all of the trajectories. Ion solvation of the charged residue side chains is transient, on the scale of a few picoseconds. There is no evidence for persistent specific ion salt bridges to charged residues. For 0 mM, only certain basic region positively charged residues are substantially Cl(-) ion salt bridged. For 200 mM, in addition, some basic region negatively (positively) charged residues are salt bridged to Na(+) (Cl(-)) ions. The different ion solvation patterns at the two ion concentrations are not greatly temperature sensitive, and the conformational sampling found in the MD is remarkably unperturbed by ion concentration and/or temperature.

Structural basis of RNA binding by leucine zipper GCN4

Recently, we showed that leucine zipper (LZ) motifs of basic leucine zipper (bZIP) transcription factors GCN4 and c-Jun are capable of catalyzing degradation of RNA (Nikolaev et al., PLoS ONE 2010; 5:e10765). This observation is intriguing given the tight regulation of RNA turnover control and the antiquity of bZIP transcription factors. To support further mechanistic studies, herein, we elucidated RNA binding interface of the GCN4 leucine zipper motif from yeast. Solution NMR experiments showed that the LZ-RNA interaction interface is located in the first two heptads of LZ moiety, and that only the dimeric (coiled coil) LZ conformation is capable of binding RNA. Site-directed mutagenesis of the LZ-GCN4 RNA binding interface showed that substrate binding is facilitated by lysine and arginine side chains, and that at least one nucleophilic residue is located in proximity to the RNA phosphate backbone. Further studies in the context of full-length bZIP factors are envisaged to address the biological relevance of LZ RNase activity.

A hamiltonian replica exchange method for building protein-protein interfaces applied to a leucine zipper

Leucine zippers consist of alpha helical monomers dimerized (or oligomerized) into alpha superhelical structures known as coiled coils. Forming the correct interface of a dimer from its monomers requires an exploration of configuration space focused on the side chains of one monomer that must interdigitate with sites on the other monomer. The aim of this work is to generate good interfaces in short simulations starting from separated monomers. Methods are developed to accomplish this goal based on an extension of a previously introduced [Su and Cukier, J. Phys. Chem. B 113, 9595, (2009)] hamiltonian temperature replica exchange method (HTREM), which scales the hamiltonian in both potential and kinetic energies that was used for the simulation of dimer melting curves. The new method, HTREM_MS (MS designates mean square), focused on interface formation, adds restraints to the hamiltonians for all but the physical system, which is characterized by the normal molecular dynamics force field at the desired temperature. The restraints in the nonphysical systems serve to prevent the monomers from separating too far, and have the dual aims of enhancing the sampling of close in configurations and breaking unwanted correlations in the restrained systems. The method is applied to a 31-residue truncation of the 33-residue leucine zipper (GCN4-p1) of the yeast transcriptional activator GCN4. The monomers are initially separated by a distance that is beyond their capture length. HTREM simulations show that the monomers oscillate between dimerlike and monomerlike configurations, but do not form a stable interface. HTREM_MS simulations result in the dimer interface being faithfully reconstructed on a 2 ns time scale. A small number of systems (one physical and two restrained with modified potentials and higher effective temperatures) are sufficient. An in silico mutant that should not dimerize because it lacks charged residues that provide electrostatic stabilization of the dimer does not with HTREM_MS, giving confidence in the method. The interface formation time scale is sufficiently short that using HTREM_MS as a screening tool to validate leucine zipper design methods may be feasible.

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.

Essential role of coiled coils for aggregation and activity of Q/N-rich prions and PolyQ proteins

The functional switch of glutamine/asparagine (Q/N)-rich prions and the neurotoxicity of polyQ-expanded proteins involve complex aggregation-prone structural transitions, commonly presumed to be forming β sheets. By analyzing sequences of interaction partners of these proteins, we discovered a recurrent presence of coiled-coil domains both in the partners and in segments that flank or overlap Q/N-rich and polyQ domains. Since coiled coils can mediate protein interactions and multimerization, we studied their possible involvement in Q/N-rich and polyQ aggregations. Using circular dichroism and chemical crosslinking, we found that Q/N-rich and polyQ peptides form α-helical coiled coils in vitro and assemble into multimers. Using structure-guided mutagenesis, we found that coiled-coil domains modulate in vivo properties of two Q/N-rich prions and polyQ-expanded huntingtin. Mutations that disrupt coiled coils impair aggregation and activity, whereas mutations that enhance coiled-coil propensity promote aggregation. These findings support a coiled-coil model for the functional switch of Q/N-rich prions and for the pathogenesis of polyQ-expansion diseases.

The folding of single domain proteins--have we reached a consensus?

Rather than stressing the most recent advances in the field, this review highlights the fundamental topics where disagreement remains and where adequate experimental data are lacking. These topics include properties of the denatured state and the role of residual structure, the nature of the fundamental steps and barriers, the extent of pathway heterogeneity and non-native interactions, recent comparisons between theory and experiment, and finally, dynamical properties of the folding reaction.

Designing the folding mechanics of coiled coils

Naturally occurring coiled coils are often not homogeneous throughout their entire structure but rather interrupted by sequence discontinuities and non-coiled-coil-forming subsegments. We apply atomic force microscopy to locally probe the mechanical folding/unfolding process of a well-understood model coiled coil when unstructured subsegments with different sizes are added. We find that the refolding force decreases from 7.8 pN with increasing size of the added unstructured subsegment, while the unfolding properties of the model coiled coil remain unchanged. We show that this behavior results from the increased size of the nucleation seed which has to form before further coiled-coil folding can proceed. Since the nucleation seed size is linked to the width of the energetic folding barrier, we are able to directly measure the dependence of folding forces on the barrier width. Our results allow the design of coiled coils with designated refolding forces by simply adjusting the nucleation seed size.

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.

Temperature and chemical denaturant dependence of forced unfolding of titin I27

Single-molecule force measurement opens a new door for investigating detailed biomolecular interactions and their thermodynamic properties by pulling molecules apart while monitoring the force exerted on them. Recent advances in the nonequilibrium work theorem allows one to determine the free-energy landscapes of these events. Such information is valuable for understanding processes such as protein and RNA folding and receptor-ligand binding. Here, we used force as a physical parameter under the traditional chemical and temperature denaturing environment to alter the protein folding energy landscape and compared the change in the unfolding free-energy barrier of the I27 domain of human cardiac titin. We found that the trends in protein unfolding free-energy barriers are consistent for single-molecule force measurements and bulk chemical and temperature studies. The results suggest that the information from single-molecule pulling experiments are meaningful and useful for understanding the mechanism of folding of titin I27.

Molecular dynamics guided study of salt bridge length dependence in both fluorinated and non-fluorinated parallel dimeric coiled-coils

The alpha-helical coiled-coil is one of the most common oligomerization motifs found in both native and engineered proteins. To better understand the stability and dynamics of the coiled-coil motifs, including those modified by fluorination, several fluorinated and nonfluorinated parallel dimeric coiled-coil protein structures were designed and modeled. We also attempt to investigate how changing the length and geometry of the important stabilizing salt bridges influences the coiled-coil protein structure. Molecular dynamics (MD) and free energy simulations with AMBER used a particle mesh Ewald treatment of the electrostatics in explicit TIP3P solvent with balanced force field treatments. Preliminary studies with legacy force fields (ff94, ff96, and ff99) show a profound instability of the coiled-coil structures in short MD simulation. Significantly, better behavior is evident with the more balanced ff99SB and ff03 protein force fields. Overall, the results suggest that the coiled-coil structures can readily accommodate the larger acidic arginine or S-2,7-diaminoheptanedoic acid mutants in the salt bridge, whereas substitution of the smaller L-ornithine residue leads to rapid disruption of the coiled-coil structure on the MD simulation time scale. This structural distortion of the secondary structure allows both the formation of large hydration pockets proximal to the charged groups and within the hydrophobic core. Moreover, the increased structural fluctuations and movement lead to a decrease in the water occupancy lifetimes in the hydration pockets. In contrast, analysis of the hydration in the stable dimeric coiled-coils shows high occupancy water sites along the backbone residues with no water occupancy in the hydrophobic core, although transitory water interactions with the salt bridge residues are evident. The simulations of the fluorinated coiled-coils suggest that in some cases fluorination electrostatically stabilizes the intermolecular coiled-coil salt bridges. Structural analyses also reveal different side chain rotamer preferences for leucine when compared with 5,5,5,5',5',5'-hexafluoroleucine mutants. These observed differences in the side chain rotamer populations suggest differential changes in the side chain conformational entropy upon coiled-coil formation when the protein is fluorinated. The free energy of hydration of the isolated 5,5,5,5',5',5'-hexafluoroleucine amino acid is calculated to be 1.1 kcal/mol less stable than leucine; this hydrophobic penalty in the monomer may provide a driving force for coiled-coil dimer formation. Estimation of the ellipticity at 222 nm from a series of snapshots from the MD simulations with DicroCalc shows distinct increases in the ellipticity when the coiled-coil is fluorinated, which suggests that the helicity in the folded coiled-coils is greater when fluorinated.

Quantifying the structural requirements of the folding transition state of protein A and other systems

The B-domain of protein A is a small three-helix bundle that has been the subject of considerable experimental and theoretical investigation. Nevertheless, a unified view of the structure of the transition-state ensemble (TSE) is still lacking. To characterize the TSE of this surprisingly challenging protein, we apply a combination of psi analysis (which probes the role of specific side-chain to side-chain contacts) and kinetic H/D amide isotope effects (which measures hydrogen-bond content), building upon previous studies using mutational phi analysis (which probes the energetic influence of side-chain substitutions). The second helix is folded in the TSE, while helix formation appears just at the carboxy and amino termini of the first and third helices, respectively. The experimental data suggest a homogenous yet plastic TS with a native-like topology. This study generalizes our earlier conclusion, based on two larger alpha/beta proteins, that the TSEs of most small proteins achieve approximately 70% of their native state's relative contact order. This high percentage limits the degree of possible TS heterogeneity and requires a reevaluation of the structural content of the TSE of other proteins, especially when they are characterized as small or polarized.

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.

Molecular basis of coiled-coil formation

Coiled coils have attracted considerable interest as design templates in a wide range of applications. Successful coiled-coil design strategies therefore require a detailed understanding of coiled-coil folding. One common feature shared by coiled coils is the presence of a short autonomous helical folding unit, termed "trigger sequence," that is indispensable for folding. Detailed knowledge of trigger sequences at the molecular level is thus key to a general understanding of coiled-coil formation. Using a multidisciplinary approach, we identify and characterize here the molecular determinants that specify the helical conformation of the monomeric early folding intermediate of the GCN4 coiled coil. We demonstrate that a network of hydrogen-bonding and electrostatic interactions stabilize the trigger-sequence helix. This network is rearranged in the final dimeric coiled-coil structure, and its destabilization significantly slows down GCN4 leucine zipper folding. Our findings provide a general explanation for the molecular mechanism of coiled-coil formation.

All-atom ab initio folding of a diverse set of proteins

Natural proteins fold to a unique, thermodynamically dominant state. Modeling of the folding process and prediction of the native fold of proteins are two major unsolved problems in biophysics. Here, we show successful all-atom ab initio folding of a representative diverse set of proteins by using a minimalist transferable-energy model that consists of two-body atom-atom interactions, hydrogen bonding, and a local sequence-energy term that models sequence-specific chain stiffness. Starting from a random coil, the native-like structure was observed during replica exchange Monte Carlo (REMC) simulation for most proteins regardless of their structural classes; the lowest energy structure was close to native-in the range of 2-6 A root-mean-square deviation (rmsd). Our results demonstrate that the successful folding of a protein chain to its native state is governed by only a few crucial energetic terms.

Secondary structure provides a template for the folding of nearby polypeptides

Although protein structures are primarily encoded by their sequences, they are also critically dependent on environmental factors such as solvents and interactions with other molecules. Here we investigate how the folding-energy landscape of a short peptide is altered by interactions with another peptide, by performing atomistic replica-exchange molecular dynamics simulations of polyalanines in various environments. We analyzed the free-energy landscapes of Ala7 and Ala8 in isolation, near an alpha-helix template, and near a beta-strand template. The isolated Ala7 and Ala8 at 270 K were mainly in polyproline II helix conformations and in equilibrium between the alpha-helix and polyproline II helix, respectively, in harmony with the experiment. Interestingly, we found remarkably strong secondary-structure "templating"; namely, the alpha-helix template enhanced alpha-helix conformation and the beta-strand template induced beta-strand conformation in the simulated Ala8. The alpha-helix template lowered the nearby dielectric constant, which strengthened hydrogen bonds in the simulated Ala8, leading to alpha-helix stabilization. The beta-strand template provided hydrogen bond positions to the simulated Ala8, sharply inducing beta-strand structure. With or without templates, the energy landscape of Ala8 is always funnel-like and centered at the alpha-helix conformation, whereas entropic contribution disfavors the alpha-helix, leading to subtle competition. Secondary-structure templating may play a critical role in protein conformation dynamics in the cellular environment.

Helix stability and hydrophobicity in the folding mechanism of the bacterial immunity protein Im9

Recent models suggest that the mechanism of protein folding is determined by the balance between the stability of secondary structural elements and the hydrophobicity of the sequence. Here we determine the role of these factors in the folding kinetics of Im9* by altering the secondary structure propensity or hydrophobicity of helices I, II or IV by the substitution of residues at solvent exposed sites. The folding kinetics of each variant were measured at pH 7.0 and 10 degrees C, under which conditions wild-type Im9* folds with two-state kinetics. We show that increasing the helicity of these sequences in regions known to be structured in the folding intermediate of Im7*, switches the folding of Im9* from a two- to three-state mechanism. By contrast, increasing the hydrophobicity of helices I or IV has no effect on the kinetic folding mechanism. Interestingly, however, increasing the hydrophobicity of solvent-exposed residues in helix II stabilizes the folding intermediate and the rate-limiting transition state, consistent with the view that this helix makes significant non-native interactions during folding. The results highlight the generic importance of intermediates in folding and show that such species can be populated by increasing helical propensity or by stabilizing inter-helix contacts through non-native interactions.