Slowing down downhill folding: a three-probe study

Chemical interactions modulate λ6-85 stability in cells

Because steric crowding is most effective when the crowding agent is similar in size to the molecule that it acts upon and the average macromolecule inside cells is much larger than a small protein or peptide, steric crowding is not predicted to affect their folding inside cells. On the other hand, chemical interactions should perturb in-cell structure and stability because they arise from interactions between the surface of the small protein or peptide and its environment. Indeed, previous in vitro measurements of the λ-repressor fragment, λ6-85 , in crowding matrices comprised of Ficoll or protein crowders support these predictions. Here, we directly quantify the in-cell stability of λ6-85 and distinguish the contribution of steric crowding and chemical interactions to its stability. Using a FRET-labeled λ6-85 construct, we find that the fragment is stabilized by 5°C in-cells compared to in vitro. We demonstrate that this stabilization cannot be explained by steric crowding because, as anticipated, Ficoll has no effect on λ6-85 stability. We find that the in-cell stabilization arises from chemical interactions, mimicked in vitro by mammalian protein extraction reagent (M-PER™). Comparison between FRET values in-cell and in Ficoll confirms that U-2 OS cytosolic crowding is reproduced at macromolecule concentrations of 15% w/v. Our measurements validate the cytomimetic of 15% Ficoll and 20% M-PER™ that we previously developed for protein and RNA folding studies. However, because the in-cell stability of λ6-85 is reproduced by 20% v/v M-PER™ alone, we predict that this simplified mixture could be a useful tool to predict the in-cell behaviors of other small proteins and peptides.

Threading single proteins through pores to compare their energy landscapes

Protein function correlates with its structural dynamics. While theoretical approaches to studying protein energy landscapes are well developed, experimental methods that enable probing these landscapes of proteins remain challenging. We used solid-state nanopores to study the translocation behavior of three mutants of a helix bundle protein and quantified the number of energetically accessible conformational states for each mutant. We found that a slower-folding mutant with access to more conformational states translocates faster than a faster-folding mutant with a smaller number of accessible states, suggesting that ease of folding and ease of translocation are at odds in this case.

Translocation of proteins is correlated with structural fluctuations that access conformational states higher in free energy than the folded state. We use electric fields at the solid-state nanopore to control the relative free energy and occupancy of different protein conformational states at the single-molecule level. The change in occupancy of different protein conformations as a function of electric field gives rise to shifts in the measured distributions of ionic current blockades and residence times. We probe the statistics of the ionic current blockades and residence times for three mutants of the λ-repressor family in order to determine the number of accessible conformational states of each mutant and evaluate the ruggedness of their free energy landscapes. Translocation becomes faster at higher electric fields when additional flexible conformations are available for threading through the pore. At the same time, folding rates are not correlated with ease of translocation; a slow-folding mutant with a low-lying intermediate state translocates faster than a faster-folding two-state mutant. Such behavior allows us to distinguish among protein mutants by selecting for the degree of current blockade and residence time at the pore. Based on these findings, we present a simple free energy model that explains the complementary relationship between folding equilibrium constants and translocation rates.

Co-translational folding of α-helical proteins: structural studies of intermediate-length variants of the λ repressor

Nascent polypeptide chains fold cotranslationally, but the atomic‐level details of this process remain unknown. Here, we report crystallographic, de novo modeling, and spectroscopic studies of intermediate‐length variants of the λ repressor N‐terminal domain. Although the ranges of helical regions of the half‐length variant were almost identical to those of the full‐length protein, the relative orientations of these helices in the intermediate‐length variants differed. Our results suggest that cotranslational folding of the λ repressor initially forms a helical structure with a transient conformation, as in the case of a molten globule state. This conformation subsequently matures during the course of protein synthesis.

Database

Structural data are available in the PDB under the accession numbers http://www.rcsb.org/pdb/search/structidSearch.do?structureId=5ZCA and http://www.rcsb.org/pdb/search/structidSearch.do?structureId=3WOA.

Globular Protein Folding In Vitro and In Vivo

In vitro, computational, and theoretical studies of protein folding have converged to paint a rich and complex energy landscape. This landscape is sensitively modulated by environmental conditions and subject to evolutionary pressure on protein function. Of these environments, none is more complex than the cell itself, where proteins function in the cytosol, in membranes, and in different compartments. A wide variety of kinetic and thermodynamics experiments, ranging from single-molecule studies to jump kinetics and from nuclear magnetic resonance to imaging on the microscope, have elucidated how protein energy landscapes facilitate folding and how they are subject to evolutionary constraints and environmental perturbation. Here we review some recent developments in the field and refer the reader to some original work and additional reviews that cover this broad topic in protein science.

Cooperative folding near the downhill limit determined with amino acid resolution by hydrogen exchange

The relationship between folding cooperativity and downhill, or barrier-free, folding of proteins under highly stabilizing conditions remains an unresolved topic, especially for proteins such as λ-repressor that fold on the microsecond timescale. Under aqueous conditions where downhill folding is most likely to occur, we measure the stability of multiple H bonds, using hydrogen exchange (HX) in a λYA variant that is suggested to be an incipient downhill folder having an extrapolated folding rate constant of 2 × 10(5) s(-1) and a stability of 7.4 kcal·mol(-1) at 298 K. At least one H bond on each of the three largest helices (α1, α3, and α4) breaks during a common unfolding event that reflects global denaturation. The use of HX enables us to both examine folding under highly stabilizing, native-like conditions and probe the pretransition state region for stable species without the need to initiate the folding reaction. The equivalence of the stability determined at zero and high denaturant indicates that any residual denatured state structure minimally affects the stability even under native conditions. Using our ψ analysis method along with mutational ϕ analysis, we find that the three aforementioned helices are all present in the folding transition state. Hence, the free energy surface has a sufficiently high barrier separating the denatured and native states that folding appears cooperative even under extremely stable and fast folding conditions.

Observation of complete pressure-jump protein refolding in molecular dynamics simulation and experiment

Density is an easily adjusted variable in molecular dynamics (MD) simulations. Thus, pressure-jump (P-jump)-induced protein refolding, if it could be made fast enough, would be ideally suited for comparison with MD. Although pressure denaturation perturbs secondary structure less than temperature denaturation, protein refolding after a fast P-jump is not necessarily faster than that after a temperature jump. Recent P-jump refolding experiments on the helix bundle λ-repressor have shown evidence of a <3 μs burst phase, but also of a ~1.5 ms "slow" phase of refolding, attributed to non-native helical structure frustrating microsecond refolding. Here we show that a λ-repressor mutant is nonetheless capable of refolding in a single explicit solvent MD trajectory in about 19 μs, indicating that the burst phase observed in experiments on the same mutant could produce native protein. The simulation reveals that after about 18.5 μs of conformational sampling, the productive structural rearrangement to the native state does not occur in a single swift step but is spread out over a brief series of helix and loop rearrangements that take about 0.9 μs. Our results support the molecular time scale inferred for λ-repressor from near-downhill folding experiments, where transition-state population can be seen experimentally, and also agrees with the transition-state transit time observed in slower folding proteins by single-molecule spectroscopy.

Fast protein folding kinetics

Fast-folding proteins have been a major focus of computational and experimental study because they are accessible to both techniques: they are small and fast enough to be reasonably simulated with current computational power, but have dynamics slow enough to be observed with specially developed experimental techniques. This coupled study of fast-folding proteins has provided insight into the mechanisms, which allow some proteins to find their native conformation well <1 ms and has uncovered examples of theoretically predicted phenomena such as downhill folding. The study of fast folders also informs our understanding of even 'slow' folding processes: fast folders are small; relatively simple protein domains and the principles that govern their folding also govern the folding of more complex systems. This review summarizes the major theoretical and experimental techniques used to study fast-folding proteins and provides an overview of the major findings of fast-folding research. Finally, we examine the themes that have emerged from studying fast folders and briefly summarize their application to protein folding in general, as well as some work that is left to do.

Multiscaled exploration of coupled folding and binding of an intrinsically disordered molecular recognition element in measles virus nucleoprotein

Numerous relatively short regions within intrinsically disordered proteins (IDPs) serve as molecular recognition elements (MoREs). They fold into ordered structures upon binding to their partner molecules. Currently, there is still a lack of in-depth understanding of how coupled binding and folding occurs in MoREs. Here, we quantified the unbound ensembles of the α-MoRE within the intrinsically disordered C-terminal domain of the measles virus nucleoprotein. We developed a multiscaled approach by combining a physics-based and an atomic hybrid model to decipher the mechanism by which the α-MoRE interacts with the X domain of the measles virus phosphoprotein. Our multiscaled approach led to remarkable qualitative and quantitative agreements between the theoretical predictions and experimental results (e.g., chemical shifts). We found that the free α-MoRE rapidly interconverts between multiple discrete partially helical conformations and the unfolded state, in accordance with the experimental observations. We quantified the underlying global folding-binding landscape. This leads to a synergistic mechanism in which the recognition event proceeds via (minor) conformational selection, followed by (major) induced folding. We also provided evidence that the α-MoRE is a compact molten globule-like IDP and behaves as a downhill folder in the induced folding process. We further provided a theoretical explanation for the inherent connections between "downhill folding," "molten globule," and "intrinsic disorder" in IDP-related systems. Particularly, we proposed that binding and unbinding of IDPs proceed in a stepwise way through a "kinetic divide-and-conquer" strategy that confers them high specificity without high affinity.

Transient helical structure during PI3K and Fyn SH3 domain folding

A growing list of proteins, including the β-sheet-rich SH3 domain, is known to transiently populate a compact α-helical intermediate before settling into the native structure. Examples have been discovered in cryogenic solvent as well as by pressure jumps. Earlier studies of λ repressor mutants showed that transient states with excess helix are robust in an all-α protein. Here we extend a previous study of src SH3 domain to two new SH3 sequences, phosphatidylinositol 3-kinase (PI3K) and a Fyn mutant, to see how robust such helix-rich transients are to sequence variations in this β-sheet fold. We quantify helical structure by circular dichroism (CD), protein compactness by small-angle X-ray scattering (SAXS), and transient helical populations by cryo-stopped-flow CD. Our results show that transient compact helix-rich intermediates are easily accessible on the folding landscape of different SH3 domains. In molecular dynamics simulations, force field errors are often blamed for transient non-native structure. We suggest that experimental examples of very fast α-rich transient misfolding could become a more subtle test for further force field improvements than observation of the native state alone.

Microsecond folding experiments and simulations: a match is made

For the past two decades, protein folding experiments have been speeding up from the second or millisecond time scale to the microsecond time scale, and full-atom simulations have been extended from the nanosecond to the microsecond and even millisecond time scale. Where the two meet, it is now possible to compare results directly, allowing force fields to be validated and refined, and allowing experimental data to be interpreted in atomistic detail. In this perspective we compare recent experiments and simulations on the microsecond time scale, pointing out the progress that has been made in determining native structures from physics-based simulations, refining experiments and simulations to provide more quantitative underlying mechanisms, and tackling the problems of multiple reaction coordinates, downhill folding, and complex underlying structure of unfolded or misfolded states.

Structural study of hNck2 SH3 domain protein in solution by circular dichroism and X-ray solution scattering

We have done conformational study of hNck2 SH3 domain by means of far-ultraviolet (far-UV) circular dichroism (CD) and X-ray solution scattering (XSS). The results indicated that the following: (1) hNck2 SH3 domain protein exhibited concentration dependent monomer-dimer transition at neutral pH, while the secondary structure of this protein was independent of the protein concentration. (2) The hNck2 SH3 domain also exhibited pH dependent monomer-dimer transition. This monomer-dimer transition was accompanied with helix-β transition of the secondary structural change. Moreover, the acid-induced conformation, which was previously studied by Liu and Song by CD and nuclear magnetic resonance (NMR), was found to be not compact, but the conformation of the protein at acidic pH was similar to the cold denatured state (C-state) reported by Yamada et al. for equine β-lactoglobulin. We calculated that a structure of the equilibrium helix-rich intermediate of the hNck2 SH3 domain by DAMMIF program.

Integrative structural modeling with small angle X-ray scattering profiles

Recent technological advances enabled high-throughput collection of Small Angle X-ray Scattering (SAXS) profiles of biological macromolecules. Thus, computational methods for integrating SAXS profiles into structural modeling are needed more than ever. Here, we review specifically the use of SAXS profiles for the structural modeling of proteins, nucleic acids, and their complexes. First, the approaches for computing theoretical SAXS profiles from structures are presented. Second, computational methods for predicting protein structures, dynamics of proteins in solution, and assembly structures are covered. Third, we discuss the use of SAXS profiles in integrative structure modeling approaches that depend simultaneously on several data types.

Unifying model for two-state and downhill protein folding

A protein-folding model is proposed at the amino acid level, in which the folding process is divided into two successive stages: the rate-determining step, dominated by the "stochastic interactions"of solvent molecules, and the rapid phase, dominated by the "order interactions"among atoms in polypeptide. The master equation approach is used to investigate the folding kinetics, and an analytical treatment of the master equation yields a simple three-parameter expression for folding time. It is found that both two-state and downhill protein-folding kinetics can be described by a unifying model.

Time resolved SAXS and RNA folding

Small angle X-ray scattering provides low resolution structural information about macromolecules in solution. When coupled with rapid mixing methods, SAXS reports time-dependent conformational changes of RNA induced by the addition of Mg(2+) to trigger folding. Thus time-resolved SAXS provides unique information about the global or overall structures of transient intermediates populated during folding. Notably, SAXS provides information about the earliest folding events, which can evade detection by other methods.

alpha-Helical burst on the folding pathway of FHA domains from Rad53 and Ki67

We investigated refolding processes of beta-sheeted protein FHA domains (FHA1 domain of Rad53 and Ki67 FHA domain) by cryo-stopped-flow (SF) method combined with far-ultraviolet (far-UV) circular dichroism (CD, the average secondary structure content) and small angle X-ray scattering (SAXS, measuring the radius of gyration). In case of FHA1 domain of Rad53, no detectable time course was observed except the initial burst on its refolding process at 4 degrees C, suggesting that the FHA1 domain of Rad53 was already refolded to its native state within the dead time of the SF apparatus and the rate of the refolding is too fast to be observed at this temperature. In contrast, there was an observable alpha-helical burst at -15 degrees C and -20 degrees C in the presence of 45% ethylene glycol (EGOH) by CD-SF. Besides, the radius of gyration (Rg) of the burst phase intermediate at -20 degrees C shows the intermediate is already compact, and the compaction process was accompanied with the decrease of alpha-helical content at the same temperature. In case of Ki67 FHA domain, ellipticity change at 222 nm was observed on its refolding pathway at -28 degrees C in the presence of 45% EGOH and 2 mM DTT, indicating that Ki67 FHA domain also takes non-native alpha-helix-rich intermediate on its folding pathway. Time-resolved SAXS experiment was done. As the signal/noise ratio is low, we could not observe the time-dependent signal change through the time course. However, the initial Rg value was obtained as 18.2 +/- 0.5 A, which is much smaller than the unfolded Rg value (26.5 +/- 1.2 A), and is slightly larger than the native one (15.9 +/- 1.8 A). These results suggest that Ki67 FHA domain also forms compact non-native alpha-helix-rich intermediate before refolding to its native beta-structure on the refolding pathway. These results are in good agreement with other beta-proteins, such as bovine beta-lactoglobulin (BLG), src SH3 domain proteins. It seems the alpha-helical burst phases appear on the folding pathway of beta-sandwiched proteins.

A survey of lambda repressor fragments from two-state to downhill folding

We survey the two-state to downhill folding transition by examining 20 lambda(6-85)* mutants that cover a wide range of stabilities and folding rates. We investigated four new lambda(6-85)* mutants designed to fold especially rapidly. Two were engineered using the core remodeling of Lim and Sauer, and two were engineered using Ferreiro et al.'s frustratometer. These proteins have probe-dependent melting temperatures as high as 80 degrees C and exhibit a fast molecular phase with the characteristic temperature dependence of the amplitude expected for downhill folding. The survey reveals a correlation between melting temperature and downhill folding previously observed for the beta-sheet protein WW domain. A simple model explains this correlation and predicts the melting temperature at which downhill folding becomes possible. An X-ray crystal structure with a 1.64-A resolution of a fast-folding mutant fragment shows regions of enhanced rigidity compared to the full wild-type protein.