Slow Folding Kinetics Research Articles

Heterotrimeric collagen helix with high specificity of assembly results in a rapid rate of folding.

The most abundant natural collagens form heterotrimeric triple helices. Synthetic mimics of collagen heterotrimers have been found to fold slowly, even compared to the already slow rates of homotrimeric helices. These prolonged folding rates are not understood. Here we compare the stabilities, specificities and folding rates of three heterotrimeric collagen mimics designed through a computationally assisted approach. The crystal structure of one ABC-type heterotrimer verified a well-controlled composition and register and elucidated the geometry of pairwise cation-π and axial and lateral salt bridges in the assembly. This collagen heterotrimer folds much faster (hours versus days) than comparable, well-designed systems. Circular dichroism and NMR data suggest the folding is frustrated by unproductive, competing heterotrimer species and these species must unwind before refolding into the thermodynamically favoured assembly. The heterotrimeric collagen folding rate is inhibited by the introduction of preformed competing triple-helical assemblies, which suggests that slow heterotrimer folding kinetics are dominated by the frustration of the energy landscape caused by competing triple helices.

Landscape Zooming toward the Prediction of RNA Cotranscriptional Folding.

RNA molecules fold as they are transcribed. Cotranscriptional folding of RNA plays a critical role in RNA functions in vivo. Present computational strategies focus on simulations where large structural changes may not be completely sampled. Here, we describe an alternative approach to predicting cotranscriptional RNA folding by zooming in and out of the RNA folding energy landscape. By classifying the RNA structural ensemble into "partitions" based on long, stable helices, we zoom out of the landscape and predict the overall slow folding kinetics from the interpartition kinetic network, and for each interpartition transition, we zoom in on the landscape to simulate the kinetics. Applications of the model to the 117-nucleotide E. coli SRP RNA and the 59-nucleotide HIV-1 TAR RNA show agreements with the experimental data and new structural and kinetic insights into biologically significant conformational switches and pathways for these important systems. This approach, by zooming in/out of an RNA folding landscape at different resolutions, might allow us to treat large RNAs in vivo with transcriptional pause, transcription speed, and other in vivo effects.

New findings on the d(TGGGAG) sequence: Surprising anti-HIV-1 activity

The biological relevance of tetramolecular G-quadruplexes especially as anti-HIV agents has been extensively reported in the literature over the last years. In the light of our recent results regarding the slow G-quadruplex folding kinetics of ODNs based on d(TGGGAG) sequence, here we report a systematic anti-HIV screening to investigate the impact of the G-quadruplex folding on their anti-HIV activity. In particular, varying the single stranded concentrations of ODNs, it has been tested a pool of ODN sample solutions with different G-quadruplex concentrations. The anti-HIV assays have been designed favouring the limited kinetics involved in the tetramolecular G4-association based on the d(TGGGAG) sequence. Aiming to determine the stoichiometry of G-quadruplex structures in the same experimental conditions of the anti-HIV assays, a native gel electrophoresis was performed. The gel confirmed the G-quadruplex formation for almost all sample solutions while showing the formation of high order G4 structures for the more concentrated ODNs solutions. The most significant result is the discovery of a potent anti-HIV activity of the G-quadruplex formed by the natural d(TGGGAG) sequence (IC50 = 14 nM) that, until now, has been reported to be completely inactive against HIV infection.

‘Z-DNA like’ fragments in RNA: a recurring structural motif with implications for folding, RNA/protein recognition and immune response

Since the work of Alexander Rich, who solved the first Z-DNA crystal structure, we have known that d(CpG) steps can adopt a particular structure that leads to forming left-handed helices. However, it is still largely unrecognized that other sequences can adopt ‘left-handed’ conformations in DNA and RNA, in double as well as single stranded contexts. These ‘Z-like’ steps involve the coexistence of several rare structural features: a C2’-endo puckering, a syn nucleotide and a lone pair–π stacking between a ribose O4’ atom and a nucleobase. This particular arrangement induces a conformational stress in the RNA backbone, which limits the occurrence of Z-like steps to ≈0.1% of all dinucleotide steps in the PDB. Here, we report over 600 instances of Z-like steps, which are located within r(UNCG) tetraloops but also in small and large RNAs including riboswitches, ribozymes and ribosomes. Given their complexity, Z-like steps are probably associated with slow folding kinetics and once formed could lock a fold through the formation of unique long-range contacts. Proteins involved in immunologic response also specifically recognize/induce these peculiar folds. Thus, characterizing the conformational features of these motifs could be a key to understanding the immune response at a structural level.

Folding Dynamics of B1 Domain of Protein G using Single Molecule Force Spectroscopy

The B1 domain of protein G is a well-studied model system of protein folding that has a mixed secondary structure, high stability and relatively slow folding kinetics. However recent measurements on the sub-millisecond timescale show significant complexity before final folding step at ∼1ms which depends on how folding is observed. Since protein folding is still a poorly understood process, different denaturants have been put into experiment to better understand the folding process of various proteins. Applying force is one natural method for denaturing proteins in order to observe the process of refolding in detail, but the time resolution of standard instruments typically make investigation of fast folding proteins difficult. In this study, two alkyne-containing unnatural amino acid residues were added to both N and C termini of GB1 and connected via click chemistry to double stranded DNA handles. The protein-DNA chimeras were then suspended between a pair of polystyrene beads held in high-resolution dual optical traps and reversible folding/unfolding events were observed both during force-extension pulling/relaxation experiments as well as under constant force feedback conditions. Results presented here show that since GB1 structure undergoes rapid switching between intermediate states of varied extension, thus a two-state model is not appropriate for understanding full folding path of this domain.

Slow and Bimolecular Folding of a De Novo Designed Monomeric Protein DS119

De novo protein design offers a unique means to test and advance our understanding of how proteins fold. However, most current design methods are native structure eccentric and folding kinetics has rarely been considered in the design process. Here, we show that a de novo designed mini-protein DS119, which folds into a βαβ structure, exhibits unusually slow and concentration-dependent folding kinetics. For example, the folding time for 50 μM of DS119 was estimated to be ∼2 s. Stopped-flow fluorescence resonance energy transfer experiments further suggested that its folding was likely facilitated by a transient dimerization process. Taken together, these results highlight the need for consideration of the entire folding energy landscape in de novo protein design and provide evidence suggesting nonnative interactions can play a key role in protein folding.

Identity of Hinge Residues Defines Stability and Kinetics of Spectrin Tetramer Interaction

Spectrin is an elongated structural protein which forms, through inter-subunit interaction, a two-dimensional lattice-work of support for the cytoskeleton and membrane of a cell. The head-to-head interaction of the α- and β subunits, the “tetramer interaction,” is of particular interest, as here two partially unstructured domains undergo a coupled binding and folding reaction to yield a new three-helix bundled domain. In mammals, spectrin tetramers are found in two isoforms, SpI and SpII. SpI, the erythrocyte spectrin, forms relatively weak tetramers that exhibit slow folding and unfolding kinetics, believed to contribute to the intrinsic flexibility and durability of the red blood cell. In contrast SpII, the non-erythrocyte spectrin, rapidly associates to forms tetramers three orders of magnitude more stable than SpI. Intriguingly, the choice of β-subunit isoform does not dramatically alter the properties of the tetramer, indicating that difference in α-subunit isoform confers the majority of increased stability and association rate. Previous work has attributed this effect to a single residue (G46 in SpIα versus R37 in SpIIα) linking the N-terminal partially structured tail to the fully folded first domain of the α-subunit. In order to investigate the importance of helical propensity and charge character to the stability and association kinetics of the spectrin tetramer, we have made a series of point mutants in SpIα and SpIIα. Our early data suggest that increasing the helical propensity of residue 46 in SpIα by mutation to alanine accelerates the association kinetics, but does not alter the thermodynamic stability of the complex, while mutation to arginine dramatically increases the stability and folding kinetics.

Using an E. coli Type 1 secretion system to secrete the mammalian, intracellular protein IFABP in its active form

A biotechnological production of proteins through protein secretion systems might be superior to the conventional cytoplasmic production, because of the absence of large amounts of proteases present in the extracellular space and the ease of purification or downstream processing. However, secretion of proteins is still a trial-and-error approach and many proteins fail to be secreted. Recently, a study of a Type 1 secretion system revealed that the folding rate of the passenger protein dictates secretion efficiency. Here, the well-known MalE failed to be secreted when fused to a C-terminal fragment of the natural substrate haemolysin A. In contrast, slow-folding mutants of MalE were secreted in high yields. However, MalE is a bacterial protein that is targeted to the periplasmic space of E. coli and possesses the intrinsic capability to cross a membrane. Therefore, we applied the same approach for another eukaryotic protein that resides in the cytoplasm. As an example, we chose the intestinal fatty acid binding protein (IFABP) and highlight the universal potential of this Type 1 secretion system to secrete proteins with slow-folding kinetics (here the G121V mutant). Finally, a one-step purification protocol was established yielding 1mg of pure IFABP G121V per liter culture supernatant. Moreover, secreted IFABP G121V was shown to reach a folded state, which is biologically active.

Estimation of protein folding free energy barriers from calorimetric data by multi-model Bayesian analysis

The realization that folding free energy barriers can be small enough to result in significant population of the species at the barrier top has sprouted in several methods to estimate folding barriers from equilibrium experiments. Some of these approaches are based on fitting the experimental thermogram measured by differential scanning calorimetry (DSC) to a one-dimensional representation of the folding free-energy surface (FES). Different physical models have been used to represent the FES: (1) a Landau quartic polynomial as a function of the total enthalpy, which acts as an order parameter; (2) the projection onto a structural order parameter (i.e. number of native residues or native contacts) of the free energy of all the conformations generated by Ising-like statistical mechanical models; and (3) mean-field models that define conformational entropy and stabilization energy as functions of a continuous local order parameter. The fundamental question that emerges is how can we obtain robust, model-independent estimates of the thermodynamic folding barrier from the analysis of DSC experiments. Here we address this issue by comparing the performance of various FES models in interpreting the thermogram of a protein with a marginal folding barrier. We chose the small α-helical protein PDD, which folds-unfolds in microseconds crossing a free energy barrier previously estimated as ~1 RT. The fits of the PDD thermogram to the various models and assumptions produce FES with a consistently small free energy barrier separating the folded and unfolded ensembles. However, the fits vary in quality as well as in the estimated barrier. Applying Bayesian probabilistic analysis we rank the fit performance using a statistically rigorous criterion that leads to a global estimate of the folding barrier and its precision, which for PDD is 1.3 ± 0.4 kJ mol(-1). This result confirms that PDD folds over a minor barrier consistent with the downhill folding regime. We have further validated the multi-model Bayesian approach through the analysis of two additional protein systems: gpW, a midsize single-domain with α + β topology that also folds in microseconds and has been previously catalogued as a downhill folder, and α-spectrin SH3, a domain of similar size but with a β-barrel fold, slow-folding kinetics and two-state-like thermodynamics. From a general viewpoint, the Bayesian analysis developed here results in a statistically robust, virtually model-independent, method to estimate the thermodynamic free-energy barriers to protein folding from DSC thermograms. Our method appears to be sufficiently accurate to consistently detect small differences in the barrier height, and thus opens up the possibility of characterizing experimentally the changes in thermodynamic folding barriers induced by single-point mutations on proteins within the downhill regime.

Folding kinetics and thermodynamics of Pseudomonas syringae effector protein AvrPto provide insight into translocation via the type III secretion system

In order to infect their hosts, many Gram-negative bacteria translocate agents of infection, called effector proteins, through the type III secretion system (TTSS) into the host cytoplasm. This process is thought to require at least partial unfolding of these agents, raising the question of how an effector protein might unfold to enable its translocation and then refold once it reaches the host cytoplasm. AvrPto is a well-studied effector protein of Pseudomonas syringae pv tomato. The presence of a readily observed unfolded population of AvrPto in aqueous solution and the lack of a known secretion chaperone make it ideal for studying the kinetic and thermodynamic characteristics that facilitate translocation. Application of Nzz exchange spectroscopy revealed a global, two-state folding equilibrium with 16% unfolded population, a folding rate of 1.8 s(-1), and an unfolding rate of 0.33 s(-1) at pH 6.1. TrAvrPto stability increases with increasing pH, with only 2% unfolded population observed at pH 7.0. The R(1) relaxation of TrAvrPto, which is sensitive to both the global anisotropy of folded TrAvrPto and slow exchange between folded and unfolded conformations, provided independent verification of the global kinetic rate constants. Given the acidic apoplast in which the pathogen resides and the more basic host cytoplasm, these results offer an intriguing mechanism by which the pH dependence of stability and slow folding kinetics of AvrPto would allow efficient translocation of the unfolded form through the TTSS and refolding into its functional folded form once inside the host.

Aberrant Folding of a Mutant Stat5b Causes Growth Hormone Insensitivity and Proteasomal Dysfunction

A predicted alanine to proline substitution in Stat5b that results in profound short stature, growth hormone insensitivity, and immunodeficiency represents the first natural mutation of this transcription factor in a human. To understand the mechanisms responsible for these pathophysiological abnormalities, we have studied the biochemical and biophysical properties of the mutant Stat5b molecule. In a cellular reconstitution model growth hormone robustly stimulated tyrosine phosphorylation and transcriptional activity of wild-type Stat5b while Stat5bA630P was minimally modified and did not promote reporter gene expression. Steady state levels of Stat5bWT were approximately 3-fold higher than Stat5bA630P in cell extracts prepared with nonionic detergents. Although initial rates of biosynthesis of both proteins were similar, pulse-chase experiments established that the apparent half-life of newly synthesized soluble Stat5bA630P was <15% of Stat5bWT (3.5 h versus >24 h). Stat5bA630P accumulated in cells primarily in cytoplasmic inclusion bodies. Structural analysis of the isolated SH2 domain containing the A630P mutation showed that it resembled the wild-type SH2 segment but that it exhibited reduced thermodynamic stability and slower folding kinetics, displayed an increased hydrophobic surface, and was prone to aggregation in solution. Our results are compatible with a model in which Stat5bA630P is an inactive transcription factor by virtue of its aberrant folding and diminished solubility triggered by a misfolded SH2 domain. The potential for aggregation and formation of cytoplasmic inclusions raises the possibility that Stat5bA630P could produce additional defects through inhibition of proteasome function.

Energetics of Substrate Binding and Catalysis by Class 1 (Glycosylhydrolase Family 47) α-Mannosidases Involved in N-Glycan Processing and Endoplasmic Reticulum Quality Control

Nascent glycoproteins are subject to quality control in the lumen of the endoplasmic reticulum (ER) where they can either be effectively folded with the aid of a collection of ER chaperones or they can be targeted for disposal in a process known as ER-associated degradation. Initiation of the ER disposal process involves selective trimming of N-glycans by ER alpha-mannosidase I and subsequent recognition by the ER degradation-enhancing alpha-mannosidase-like protein family of lectins, both members of glycosylhydrolase family 47. The kinetics and energetics of substrate binding and catalysis by members of this family were investigated here by the analysis of wild type and mutant forms of human ER alpha-mannosidase I. The contributions of several amino acid residues and an enzyme-associated Ca(2+) ion to substrate binding and catalysis were demonstrated by a combination of surface plasmon resonance and enzyme kinetic analyses. One mutant, E330Q, shown previously to alter general acid function within the catalytic site, resulted in an enzyme that possessed increased glycan binding affinity but compromised glycan hydrolysis. This mutant protein was used in a series of glycan binding studies with a library of mannose-containing ligands to examine the energetics of Man(9)GlcNAc(2) substrate interactions. These studies provide a framework for understanding the nature of the unusual substrate interactions within the family 47 mannosidases involved in glycan maturation and ER-associated glycoprotein degradation.

Identification of proline residues responsible for the slow folding kinetics in pectate lyase C by mutagenesis.

The folding mechanism of pectate lyase C (pelC) involves two slow phases that have been attributed to proline isomerization. To have a more detailed and complete understanding of the folding mechanism, experiments have been carried out to identify the prolyl-peptide bonds responsible for the slow kinetics. Site-directed mutagenesis has been used to mutate each of the prolines in pelC to alanine or valine. It has been determined that isomerization of the Leu219-Pro220 peptide bond is responsible for the slowest folding phase observed. The mutant P220A shows kinetic behavior that is identical to the wild-type protein except that the 46-s phase is eliminated. The Leu219-Pro220 peptide bond is cis in the native enzyme. An analysis of the free energy of unfolding of this mutant indicates that the mutation destabilizes the protein by about 4 kcal/mol. However, it appears that the major refolding pathways are unaltered. Further mutations were carried out in order to assign the peptide bond responsible for the 21-s folding phase in pelC. Mutation of the remaining 11 prolines, which are trans in the native enzyme, resulted in no significant changes in the kinetic folding behavior. The conclusion from these experiments is that the 21-s phase involves isomerization of more than one prolyl-peptide bond with similar activation energies.

Folding Kinetics of the All-β-sheet Protein Human Basic Fibroblast Growth Factor, a Structural Homolog of Interleukin-1β

The refolding and unfolding kinetics of the all-beta-sheet protein human basic fibroblast growth factor (hFGF-2) were studied by fluorescence spectroscopy. The kinetics of the unfolding transition are monophasic. The refolding reaction at high and low guanidinium chloride (GdmCl) concentrations is best described by mono- and biphasic folding, respectively. Refolding and unfolding of hFGF-2 (155 amino acids) is very slow compared with other non-disulfide-bonded monomeric proteins of similar size. For example, the rate constant for unfolding at 4.5 mol.liter(-1) GdmCl is 0.006 s(-1), and the refolding rate constants at 0.4 mol.liter(-1) GdmCl are 0.01 s(-1) and 0.0009 s(-1) (15 degrees C, pH 7.0). A characterization of the thermodynamic nature of the folding process using transition state theory revealed that the slow refolding is almost exclusively controlled by entropic factors, namely the strong loss of conformational freedom during refolding. The rate of the slow unfolding kinetics is mainly (and at low denaturant concentrations exclusively) controlled by the large positive change in enthalpy. hFGF-2 shows similar slow folding kinetics to that of its structural homolog interleukin-1beta. Since both proteins show very little sequence identity, it is suggested that their slow folding kinetics are determined by the complex beta-sheet arrangement of the native molecules.

Slow unfolding and refolding kinetics of the mesophilic Rop wild-type protein in the transition range.

We describe the guanidinium hydrochloride induced folding kinetics of the four-helix-bundle protein Rop wild-type (wt) under equilibrium conditions at three temperatures. The choice of appropriate denaturant conditions inside the transition range permitted, in combination with equilibrium transition curves, the determination of both unfolding and refolding rate constants. The ratio of the rate constants at zero denaturant concentration provided equilibrium constants and standard free energy changes that are in good agreement with values obtained in previous differential scanning calorimetry studies. The DeltaG0D values for 19, 25 and 40 degrees C calculated from the present kinetic studies are, respectively, 66.8, 70.8 and 57.2 kJ.mol-1. The unfolding reactions are extremely slow under these conditions. Equilibrium was reached only after 18, 12 and 6 days at 19, 25 and 40 degrees C. These results demonstrate that for Rop wt high stability correlates with slow folding kinetics.

Chaperonin-facilitated protein folding: optimization of rate and yield by an iterative annealing mechanism.

We develop a heuristic model for chaperonin-facilitated protein folding, the iterative annealing mechanism, based on theoretical descriptions of "rugged" conformational free energy landscapes for protein folding, and on experimental evidence that (i) folding proceeds by a nucleation mechanism whereby correct and incorrect nucleation lead to fast and slow folding kinetics, respectively, and (ii) chaperonins optimize the rate and yield of protein folding by an active ATP-dependent process. The chaperonins GroEL and GroES catalyze the folding of ribulose bisphosphate carboxylase at a rate proportional to the GroEL concentration. Kinetically trapped folding-incompetent conformers of ribulose bisphosphate carboxylase are converted to the native state in a reaction involving multiple rounds of quantized ATP hydrolysis by GroEL. We propose that chaperonins optimize protein folding by an iterative annealing mechanism; they repeatedly bind kinetically trapped conformers, randomly disrupt their structure, and release them in less folded states, allowing substrate proteins multiple opportunities to find pathways leading to the most thermodynamically stable state. By this mechanism, chaperonins greatly expand the range of environmental conditions in which folding to the native state is possible. We suggest that the development of this device for optimizing protein folding was an early and significant evolutionary event.