Local versus global fold switching in protein evolution: insight from a three-letter continuous model

Simulations of a protein fold switch reveal crowding-induced population shifts driven by disordered regions

Macromolecular crowding effects on globular proteins, which usually adopt a single stable fold, have been widely studied. However, little is known about crowding effects on fold-switching proteins, which reversibly switch between distinct folds. Here we study the mutationally driven switch between the folds of GA and GB, the two 56-amino acid binding domains of protein G, using a structure-based dual-basin model. We show that, in the absence of crowders, the fold populations PA and PB can be controlled by the strengths of contacts in the two folds, κA and κB. A population balance, PA ≈ PB, is obtained for κB/κA = 0.92. The resulting model protein is subject to crowding at different packing fractions, ϕc. We find that crowding increases the GB population and reduces the GA population, reaching PB/PA ≈ 4 at ϕc = 0.44. We analyze the ϕc-dependence of the crowding-induced GA-to-GB switch using scaled particle theory, which provides a qualitative, but not quantitative, fit of our data, suggesting effects beyond a spherical description of the folds. We show that the terminal regions of the protein chain, which are intrinsically disordered only in GA, play a dominant role in the response of the fold switch to crowding effects.

Crowding-induced protein destabilization in the absence of soft attractions

It is generally assumed that volume exclusion by macromolecular crowders universally stabilizes the native states of proteins and destabilization suggests soft attractions between crowders and protein. Here we show that proteins can be destabilized even by crowders that are purely repulsive. With a coarse-grained sequence-based model, we study the folding thermodynamics of two sequences with different native folds, a helical hairpin and a β-barrel, in a range of crowder volume fractions, φ_c. We find that the native state, N, remains structurally unchanged under crowded conditions, while the size of the unfolded state, U, decreases monotonically with φ_c. Hence, for all φ_c>0, U is entropically disfavored relative to N. This entropy-centric view holds for the helical hairpin protein, which is stabilized under all crowded conditions as quantified by changes in either the folding midpoint temperature, T_m, or the free energy of folding. We find, however, that the β-barrel protein is destabilized under low-T, low-φ_c conditions. This destabilization can be understood from two characteristics of its folding: 1) a relatively compact U at T<T_m, such that U is only weakly disfavored entropically by the crowders; and 2) a transient, compact, and relatively low-energy nonnative state that has a maximum population of only a few percent at φ_c=0, but increasing monotonically with φ_c. Overall, protein destabilization driven by hard-core effects appears possible when a compaction of U leads to even a modest population of compact nonnative states that are energetically competitive with N.

The C-terminal domain of transcription factor RfaH: Folding, fold switching and energy landscape

We simulate the folding and fold switching of the C-terminal domain (CTD) of the transcription factor RfaH using an all-atom physics-based model augmented with a dual-basin structure-based potential energy term. We show that this hybrid model captures the essential thermodynamic behavior of this metamorphic domain, that is, a change in the global free energy minimum from an α-helical hairpin to a 5-stranded β-barrel upon the dissociation of the CTD from the rest of the protein. Using Monte Carlo sampling techniques, we then analyze the energy landscape of the CTD in terms of progress variables for folding toward the two folds. We find that, below the folding transition, the energy landscape is characterized by a single, dominant funnel to the native β-barrel structure. The absence of a deep funnel to the α-helical hairpin state reflects a negligible population of this fold for the isolated CTD. We observe, however, a higher α-helix structure content in the unfolded state compared to results from a similar but fold switch-incompetent version of our model. Moreover, in folding simulations started from an extended chain conformation we find transiently formed α-helical structure, occurring early in the process and disappearing as the chain progresses toward the thermally stable β-barrel state.

Effects of Topology and Sequence in Protein Folding Linked via Conformational Fluctuations

Experiments have compared the folding of proteins with different amino acid sequences but the same basic structure, or fold. Results indicate that folding is robust to sequence variations for proteins with some nonlocal folds, such as all-β, whereas the folding of more local, all-α proteins typically exhibits a stronger sequence dependence. Here, we use a coarse-grained model to systematically study how variations in sequence perturb the folding energy landscapes of three model sequences with 3α, 4β + α, and β-barrel folds, respectively. These three proteins exhibit folding features in line with experiments, including expected rank order in the cooperativity of the folding transition and stability-dependent shifts in the location of the free-energy barrier to folding. Using a generalized-ensemble simulation approach, we determine the thermodynamics of around 2000 sequence variants representing all possible hydrophobic or polar single- and double-point mutations. From an analysis of the subset of stability-neutral mutations, we find that folding is perturbed in a topology-dependent manner, with the β-barrel protein being the most robust. Our analysis shows, in particular, that the magnitude of mutational perturbations of the transition state is controlled in part by the size or "width" of the underlying conformational ensemble. This result suggests that the mutational robustness of the folding of the β-barrel protein is underpinned by its conformationally restricted transition state ensemble, revealing a link between sequence and topological effects in protein folding.

Dimension conversion and scaling of disordered protein chains

To extract protein dimension and energetics information from single-molecule fluorescence resonance energy transfer spectroscopy (smFRET) data, it is essential to establish the relationship between the distributions of the radius of gyration (Rg) and the end-to-end (donor-to-acceptor) distance (Ree). Here, we performed a coarse-grained molecular dynamics simulation to obtain a conformational ensemble of denatured proteins and intrinsically disordered proteins. For any disordered chain with fixed length, there is an excellent linear correlation between the average values of Rg and Ree under various solvent conditions, but the relationship deviates from the prediction of a Gaussian chain. A modified conversion formula was proposed to analyze smFRET data. The formula reduces the discrepancy between the results obtained from FRET and small-angle X-ray scattering (SAXS). The scaling law in a coil-globule transition process was examined where a significant finite-size effect was revealed, i.e., the scaling exponent may exceed the theoretical critical boundary [1/3, 3/5] and the prefactor changes notably during the transition. The Sanchez chain model was also tested and it was shown that the mean-field approximation works well for expanded chains.