Kinetic traps in the folding of beta alpha-repeat proteins: CheY initially misfolds before accessing the native conformation

A Hydrophobic Core Stabilizes the Residual Structure in the RRM2 Intermediate State of the ALS-linked Protein TDP-43

Folding intermediates mediate both protein folding and the misfolding and aggregation observed in human diseases, including amyotrophic lateral sclerosis (ALS), and are prime targets for therapeutic interventions. In this study, we identified the core nucleus of structure for a folding intermediate in the second RNA recognition motif (RRM2) of the ALS-linked RNA-binding protein, TDP-43 (TAR DNA-binding protein-43), using a combination of experimental and computational approaches. Urea equilibrium unfolding studies revealed that the RRM2 intermediate state consists of collapsed residual secondary structure localized to the N-terminal half of RRM2, while the C-terminus is largely disordered. Steered molecular dynamics simulations and mutagenesis studies yielded key stabilizing hydrophobic contacts that, when mutated to alanine, severely disrupt the overall fold of RRM2. In combination, these findings suggest a role for this RRM intermediate in normal TDP-43 function as well as serving as a template for misfolding and aggregation through the low stability and non-native secondary structure.

Thermodynamic and kinetic analysis of the LAO binding protein and its isolated domains reveal non-additivity in stability, folding and function

Substrate-binding proteins (SBPs) are used by organisms from the three domains of life for transport and signalling. SBPs are composed of two domains that collectively trap ligands with high affinity and selectivity. To explore the role of the domains and the integrity of the hinge region between them in the function and conformation of SBPs, here, we describe the ligand binding, conformational stability and folding kinetics of the Lysine Arginine Ornithine (LAO) binding protein from Salmonella thiphimurium and constructs corresponding to its two independent domains. LAO is a class II SBP formed by a continuous and a discontinuous domain. Contrary to the expected behaviour based on their connectivity, the discontinuous domain shows a stable native-like structure that binds l-arginine with moderate affinity, whereas the continuous domain is barely stable and shows no detectable ligand binding. Regarding folding kinetics, studies of the entire protein revealed the presence of at least two intermediates. While the unfolding and refolding of the continuous domain exhibited only a single intermediate and simpler and faster kinetics than LAO, the folding mechanism of the discontinuous domain was complex and involved multiple intermediates. These findings suggest that in the complete protein the continuous domain nucleates folding and that its presence funnels the folding of the discontinuous domain avoiding nonproductive interactions. The strong dependence of the function, stability and folding pathway of the lobes on their covalent association is most likely the result of the coevolution of both domains as a single unit.

How soluble misfolded proteins bypass chaperones at the molecular level

Subpopulations of soluble, misfolded proteins can bypass chaperones within cells. The extent of this phenomenon and how it happens at the molecular level are unknown. Through a meta-analysis of the experimental literature we find that in all quantitative protein refolding studies there is always a subpopulation of soluble but misfolded protein that does not fold in the presence of one or more chaperones, and can take days or longer to do so. Thus, some misfolded subpopulations commonly bypass chaperones. Using multi-scale simulation models we observe that the misfolded structures that bypass various chaperones can do so because their structures are highly native like, leading to a situation where chaperones do not distinguish between the folded and near-native-misfolded states. More broadly, these results provide a mechanism by which long-time scale changes in protein structure and function can persist in cells because some misfolded states can bypass components of the proteostasis machinery.

Contrasting effect of ficoll on apo and holo forms of bacterial chemotaxis protein Y: Selective destabilization of the conformationally altered holo form

Chemotaxis Y (CheY), upon metal binding, displays a drastic alteration in its structure and stability. This premise prompted us to study the effect of crowding on the two conformationally distinct states of the same test protein. A comparative analysis on the structure and thermal stability in the presence and absence of the macromolecular crowder, ficoll, and its monomeric unit, sucrose, revealed a contrasting effect of ficoll on the apo and holo forms. In the presence of ficoll while the thermal stability (T_m) of the apo form is enhanced, the thermal stability of the holo form is reduced. The selective lowering of T_m for the holo form in the combined presence of ficoll and sucrose and not in sucrose alone suggests that the contrasting effect is due to the macromolecular nature of ficoll. Since metal-protein interaction remains unperturbed in the presence of ficoll and Mg²⁺ sequestration is ruled out in a systematic manner the alternative possibility for the exclusive reduction in the thermal stability of the holo form is the ficoll-induced modulation of the relative population of apo and holo forms of CheY.

The Non-native Disulfide-Bond-Containing Landscape Orthogonal to the Oxidative Protein-Folding Trajectory: A Necessary Evil?

Oxidative protein folding describes the process by which disulfide-bond-containing proteins mature from their ribosomal, fully reduced and unfolded, origins. Over the past 40 years, a number of exemplar proteins including bovine pancreatic ribonuclease A (RNaseA), bovine pancreatic trypsin inhibitor (BPTI), and hen egg-white lysozyme (HEWL), among others, have provided rich insight into the nature of the intermolecular interactions that drive the formation of the native, biologically active fold. In this Review Article, we revisit the oxidative folding process of RNase A with a focus on reconciling the role of non-native disulfide-bond-containing species that populate the oxidative folding landscape. Toward gaining such an understanding, we project the regeneration pathway onto a Cartesian coordinate system. This helps not only to recognize the magnitude of the seemingly "fruitless", non-native disulfide-bond-containing species that lie orthogonal to the "native-protein-forming" reaction progress but also to reconcile a role for their existence in the regenerative trajectory. Finally, we superimpose the folding funnel onto the regeneration trajectory to draw parallels between oxidative folders and conformational folders (proteins that lack disulfide bonds). The overall objective is to provide the reader with a semi-quantitative description of oxidative protein folding and the barriers to successful regeneration while underscoring a role of seemingly fruitless intermediates.

Crowding by Poly(ethylene glycol) Destabilizes Chemotaxis Protein Y (CheY)

The prevailing understanding of various aspects of biochemical processes, including folding, stability, intermolecular interactions, and the binding of metals, substrates, and inhibitors, is derived from studies carried out under dilute and homogeneous conditions devoid of a crowding-related environment. The effect of crowding-induced modulation on the structure and stability of native and magnesium-dependent Chemotaxis Y (CheY), a bacterial signaling protein, was probed in the presence and absence of poly(ethylene glycol) (PEG). A combined analysis from circular dichroism, intrinsic and extrinsic fluorescence, and tryptophan fluorescence lifetime changes indicates that PEG perturbs the structure but leaves the thermal stability largely unchanged. Intriguingly, while the stability of the protein is enhanced in the presence of magnesium under dilute buffer conditions, PEG-induced crowding leads to reduced thermal stability in the presence of magnesium. Nuclear magnetic resonance (NMR) chemical shift perturbations and resonance broadening for a subset of residues indicate that PEG interacts specifically with a subset of hydrophilic and hydrophobic residues found predominantly in α helices, β strands, and in the vicinity of the metal-binding region. Thus, PEG prompted conformational perturbation, presumably provides a different situation for magnesium interaction, thereby perturbing the magnesium-prompted stability. In summary, our results highlight the dominance of enthalpic contributions between PEG and CheY via both hydrophilic and hydrophobic interactions, which can subtly affect the conformation, modulating the metal-protein interaction and stability, implying that in the context of cellular situation, structure, stability, and magnesium binding thermodynamics of CheY may be different from those measured in dilute solution.

The Folding Pathway of 6aJL2

One-third of the reported cases of light chain amyloidosis are related to the germ line λ6 family; remarkably, healthy individuals express this type of protein in just 2% of the peripheral blood and bone marrow B-cells. The appearance of the disease has been related to the inherent properties of this protein family. A recombinant representative model for λ6 proteins called 6aJL2 containing the amino acid sequence encoded by the 6a and JL2 germ line genes was previously designed and synthesized to study the properties of this family. Previous work on 6aJL2 suggested a simple two-state folding model at 25 °C; no intermediate could be identified either by kinetics or by fluorescence and circular dichroism equilibrium studies, although the presence of an intermediate that is populated at ∼2.4 M urea was suggested by size exclusion chromatography. In this study we employed classic equilibrium and kinetic experiments and analysis to elucidate the detailed folding mechanism of this protein. We identify species that are kinetically accessible and/or are populated at equilibrium. We describe the presence of intermediate and native-like species and propose a five-species folding mechanism at 25 °C at short incubation times, similar to and consistent with those observed in other proteins of this fold. The formation of intermediates in the mechanism of 6aJL2 is faster than that proposed for a V_κ light chain, which could be an important distinction in the amyloidogenic potential of both germ lines.

Mechanism of the small ATP-independent chaperone Spy is substrate specific

ATP-independent chaperones are usually considered to be holdases that rapidly bind to non-native states of substrate proteins and prevent their aggregation. These chaperones are thought to release their substrate proteins prior to their folding. Spy is an ATP-independent chaperone that acts as an aggregation inhibiting holdase but does so by allowing its substrate proteins to fold while they remain continuously chaperone bound, thus acting as a foldase as well. The attributes that allow such dual chaperoning behavior are unclear. Here, we used the topologically complex protein apoflavodoxin to show that the outcome of Spy's action is substrate specific and depends on its relative affinity for different folding states. Tighter binding of Spy to partially unfolded states of apoflavodoxin limits the possibility of folding while bound, converting Spy to a holdase chaperone. Our results highlight the central role of the substrate in determining the mechanism of chaperone action.

Frustration and folding of a TIM barrel protein

Triosephosphate isomerase (TIM) barrel proteins have not only a conserved architecture that supports a myriad of enzymatic functions, but also a conserved folding mechanism that involves on- and off-pathway intermediates. Although experiments have proven to be invaluable in defining the folding free-energy surface, they provide only a limited understanding of the structures of the partially folded states that appear during folding. Coarse-grained simulations employing native centric models are capable of sampling the entire energy landscape of TIM barrels and offer the possibility of a molecular-level understanding of the readout from sequence to structure. We have combined sequence-sensitive native centric simulations with small-angle X-ray scattering and time-resolved Förster resonance energy transfer to monitor the formation of structure in an intermediate in the Sulfolobus solfataricus indole-3-glycerol phosphate synthase TIM barrel that appears within 50 μs and must at least partially unfold to achieve productive folding. Simulations reveal the presence of a major and 2 minor folding channels not detected in experiments. Frustration in folding, i.e., backtracking in native contacts, is observed in the major channel at the initial stage of folding, as well as late in folding in a minor channel before the appearance of the native conformation. Similarities in global and pairwise dimensions of the early intermediate, the formation of structure in the central region that spreads progressively toward each terminus, and a similar rate-limiting step in the closing of the β-barrel underscore the value of combining simulation and experiment to unravel complex folding mechanisms at the molecular level.

Spontaneous refolding of the large multidomain protein malate synthase G proceeds through misfolding traps

Most protein folding studies until now focus on single domain or truncated proteins. Although great insights in the folding of such systems has been accumulated, very little is known regarding the proteins containing multiple domains. It has been shown that the high stability of domains, in conjunction with inter-domain interactions, manifests as a frustrated energy landscape, causing complexity in the global folding pathway. However, multidomain proteins despite containing independently foldable, loosely cooperative sections can fold into native states with amazing speed and accuracy. To understand the complexity in mechanism, studies were conducted previously on the multidomain protein malate synthase G (MSG), an enzyme of the glyoxylate pathway with four distinct and adjacent domains. It was shown that the protein refolds to a functionally active intermediate state at a fast rate, which slowly produces the native state. Although experiments decoded the nature of the intermediate, a full description of the folding pathway was not elucidated. In this study, we use a battery of biophysical techniques to examine the protein's folding pathway. By using multiprobe kinetics studies and comparison with the equilibrium behavior of protein against urea, we demonstrate that the unfolded polypeptide undergoes conformational compaction to a misfolded intermediate within milliseconds of refolding. The misfolded product appears to be stabilized under moderate denaturant concentrations. Further folding of the protein produces a stable intermediate, which undergoes partial unfolding-assisted large segmental rearrangements to achieve the native state. This study reveals an evolved folding pathway of the multidomain protein MSG, which involves surpassing the multiple misfolding traps during refolding.

Structural Rearrangement upon Fragmentation of the Stability Core of the ALS-Linked Protein TDP-43

Amyotrophic lateral sclerosis (ALS) is the most common adult degenerative motor neuron disease. Experimental evidence indicates a direct role of transactive-response DNA-binding protein 43 (TDP-43) in the pathology of ALS and other neurodegenerative diseases. TDP-43 has been identified as a major component of cytoplasmic inclusions in patients with sporadic ALS; however, the molecular basis of the disease mechanism is not yet fully understood. Fragmentation within the second RNA recognition motif (RRM2) of TDP-43 has been observed in patient tissues and may play a role in the formation of aggregates in disease. To determine the structural and dynamical changes resulting from the truncation that could lead to aggregation and toxicity, we performed molecular dynamics simulations of the full-length RRM2 domain (the stability core of TDP-43) and of a truncated variant (where residues 189-207 are deleted to mimic a site of cleavage within RRM2 found in ALS patients). Our simulations show heterogeneous structural reorganization and decreased stability of the truncated RRM2 domain compared to the full-length domain, consistent with previous experimental results. The decreased stability and structural reorganization in the truncated RRM2 result in a higher probability of protein-protein interactions through altered electrostatic surface charges and increased accessibility of hydrophobic residues (including the nuclear export sequence), providing a rationale for the increased cytoplasmic aggregation of RRM2 fragments seen in sporadic ALS patients.

Folding of proteins with a flavodoxin-like architecture

The flavodoxin-like fold is a protein architecture that can be traced back to the universal ancestor of the three kingdoms of life. Many proteins share this α-β parallel topology and hence it is highly relevant to illuminate how they fold. Here, we review experiments and simulations concerning the folding of flavodoxins and CheY-like proteins, which share the flavodoxin-like fold. These polypeptides tend to temporarily misfold during unassisted folding to their functionally active forms. This susceptibility to frustration is caused by the more rapid formation of an α-helix compared to a β-sheet, particularly when a parallel β-sheet is involved. As a result, flavodoxin-like proteins form intermediates that are off-pathway to native protein and several of these species are molten globules (MGs). Experiments suggest that the off-pathway species are of helical nature and that flavodoxin-like proteins have a nonconserved transition state that determines the rate of productive folding. Folding of flavodoxin from Azotobacter vinelandii has been investigated extensively, enabling a schematic construction of its folding energy landscape. It is the only flavodoxin-like protein of which cotranslational folding has been probed. New insights that emphasize differences between in vivo and in vitro folding energy landscapes are emerging: the ribosome modulates MG formation in nascent apoflavodoxin and forces this polypeptide toward the native state.

Atomistic structural ensemble refinement reveals non-native structure stabilizes a sub-millisecond folding intermediate of CheY

The dynamics of globular proteins can be described in terms of transitions between a folded native state and less-populated intermediates, or excited states, which can play critical roles in both protein folding and function. Excited states are by definition transient species, and therefore are difficult to characterize using current experimental techniques. Here, we report an atomistic model of the excited state ensemble of a stabilized mutant of an extensively studied flavodoxin fold protein CheY. We employed a hybrid simulation and experimental approach in which an aggregate 42 milliseconds of all-atom molecular dynamics were used as an informative prior for the structure of the excited state ensemble. This prior was then refined against small-angle X-ray scattering (SAXS) data employing an established method (EROS). The most striking feature of the resulting excited state ensemble was an unstructured N-terminus stabilized by non-native contacts in a conformation that is topologically simpler than the native state. Using these results, we then predict incisive single molecule FRET experiments as a means of model validation. This study demonstrates the paradigm of uniting simulation and experiment in a statistical model to study the structure of protein excited states and rationally design validating experiments.

Using the folding landscapes of proteins to understand protein function

Proteins fold on a biologically-relevant timescale because of a funnel-shaped energy landscape. This landscape is sculpted through evolution by selecting amino-acid sequences that stabilize native interactions while suppressing stable non-native interactions that occur during folding. However, there is strong evolutionary selection for functional residues and these cannot be chosen to optimize folding. Their presence impacts the folding energy landscape in a variety of ways. Here, we survey the effects of functional residues on folding by providing several examples. We then review how such effects can be detected computationally and be used as assays for protein function. Overall, an understanding of how functional residues modulate folding should provide insights into the design of natural proteins and their homeostasis.

Clusters of isoleucine, leucine, and valine side chains define cores of stability in high-energy states of globular proteins: Sequence determinants of structure and stability

Measurements of protection against exchange of main chain amide hydrogens (NH) with solvent hydrogens in globular proteins have provided remarkable insights into the structures of rare high-energy states that populate their folding free-energy surfaces. Lacking, however, has been a unifying theory that rationalizes these high-energy states in terms of the structures and sequences of their resident proteins. The Branched Aliphatic Side Chain (BASiC) hypothesis has been developed to explain the observed patterns of protection in a pair of TIM barrel proteins. This hypothesis supposes that the side chains of isoleucine, leucine, and valine (ILV) residues often form large hydrophobic clusters that very effectively impede the penetration of water to their underlying hydrogen bond networks and, thereby, enhance the protection against solvent exchange. The linkage between the secondary and tertiary structures enables these ILV clusters to serve as cores of stability in high-energy partially folded states. Statistically significant correlations between the locations of large ILV clusters in native conformations and strong protection against exchange for a variety of motifs reported in the literature support the generality of the BASiC hypothesis. The results also illustrate the necessity to elaborate this simple hypothesis to account for the roles of adjacent hydrocarbon moieties in defining stability cores of partially folded states along folding reaction coordinates.

Gradual Folding of an Off-Pathway Molten Globule Detected at the Single-Molecule Level

Molten globules (MGs) are compact, partially folded intermediates that are transiently present during folding of many proteins. These intermediates reside on or off the folding pathway to native protein. Conformational evolution during folding of off-pathway MGs is largely unexplored. Here, we characterize the denaturant-dependent structure of apoflavodoxin's off-pathway MG. Using single-molecule fluorescence resonance energy transfer (smFRET), we follow conversion of unfolded species into MG down to denaturant concentrations that favor formation of native protein. Under strongly denaturing conditions, fluorescence resonance energy transfer histograms show a single peak, arising from unfolded protein. The smFRET efficiency distribution shifts to higher value upon decreasing denaturant concentration because the MG folds. Strikingly, upon approaching native conditions, the fluorescence resonance energy transfer efficiency of the MG rises above that of native protein. Thus, smFRET exposes the misfolded nature of apoflavodoxin's off-pathway MG. We show that conversion of unfolded into MG protein is a gradual, second-order-like process that simultaneously involves separate regions within the polypeptide.

Non-native structure appears in microseconds during the folding of E. coli RNase H

The folding pathway of Escherichia coli RNase H is one of the best experimentally characterized for any protein. In spite of this, spectroscopic studies have never captured the earliest events. Using continuous-flow microfluidic mixing, we have now observed the first several milliseconds of folding by monitoring the tryptophan fluorescence lifetime (60 μs dead time). Two folding intermediates are observed, the second of which is the previously characterized I(core) millisecond intermediate. The new earlier intermediate is likely on-pathway and appears to have long-range non-native structure, providing a rare example of such non-native structure formation in a folding pathway. The tryptophan fluorescence lifetimes also suggest a deviation from native packing in the second intermediate, I(core). Similar results from a fragment of RNase H demonstrate that only half of the protein is significantly involved in this early structure formation. These studies give us a view of the formation of tertiary structure on the folding pathway, which complements previous hydrogen-exchange studies that monitored only secondary structure and observed sequential native structure formation. Our results provide detailed folding information on both a timescale and a size-scale accessible to all-atom molecular dynamics simulations of protein folding.

Modulation of frustration in folding by sequence permutation

Folding of globular proteins can be envisioned as the contraction of a random coil unfolded state toward the native state on an energy surface rough with local minima trapping frustrated species. These substructures impede productive folding and can serve as nucleation sites for aggregation reactions. However, little is known about the relationship between frustration and its underlying sequence determinants. Chemotaxis response regulator Y (CheY), a 129-amino acid bacterial protein, has been shown previously to populate an off-pathway kinetic trap in the microsecond time range. The frustration has been ascribed to premature docking of the N- and C-terminal subdomains or, alternatively, to the formation of an unproductive local-in-sequence cluster of branched aliphatic side chains, isoleucine, leucine, and valine (ILV). The roles of the subdomains and ILV clusters in frustration were tested by altering the sequence connectivity using circular permutations. Surprisingly, the stability and buried surface area of the intermediate could be increased or decreased depending on the location of the termini. Comparison with the results of small-angle X-ray-scattering experiments and simulations points to the accelerated formation of a more compact, on-pathway species for the more stable intermediate. The effect of chain connectivity in modulating the structures and stabilities of the early kinetic traps in CheY is better understood in terms of the ILV cluster model. However, the subdomain model captures the requirement for an intact N-terminal domain to access the native conformation. Chain entropy and aliphatic-rich sequences play crucial roles in biasing the early events leading to frustration in the folding of CheY.

Folding of the RNA recognition motif (RRM) domains of the amyotrophic lateral sclerosis (ALS)-linked protein TDP-43 reveals an intermediate state

Pathological alteration of TDP-43 (TAR DNA-binding protein-43), a protein involved in various RNA-mediated processes, is a hallmark feature of the neurodegenerative diseases amyotrophic lateral sclerosis and frontotemporal lobar degeneration. Fragments of TDP-43, composed of the second RNA recognition motif (RRM2) and the disordered C terminus, have been observed in cytoplasmic inclusions in sporadic amyotrophic lateral sclerosis cases, suggesting that conformational changes involving RRM2 together with the disordered C terminus play a role in aggregation and toxicity. The biophysical data collected by CD and fluorescence spectroscopies reveal a three-state equilibrium unfolding model for RRM2, with a partially folded intermediate state that is not observed in RRM1. Strikingly, a portion of RRM2 beginning at position 208, which mimics a cleavage site observed in patient tissues, increases the population of this intermediate state. Mutually stabilizing interactions between the domains in the tethered RRM1 and RRM2 construct reduce the population of the intermediate state and enhance DNA/RNA binding. Despite the high sequence homology of the two domains, a network of large hydrophobic residues in RRM2 provides a possible explanation for the increased stability of RRM2 compared with RRM1. The cluster analysis suggests that the intermediate state may play a functional role by enhancing access to the nuclear export signal contained within its sequence. The intermediate state may also serve as a molecular hazard linking productive folding and function with pathological misfolding and aggregation that may contribute to disease.

Balancing bond, nonbond, and gō-like terms in coarse grain simulations of conformational dynamics

Characterization of the protein conformational landscape remains a challenging problem, whether it concerns elucidating folding mechanisms, predicting native structures or modeling functional transitions. Coarse-grained molecular dynamics simulation methods enable exhaustive sampling of the energetic landscape at resolutions of biological interest. The general utility of structure-based models is reviewed along with their differing levels of approximation. Simple Gō models incorporate attractive native interactions and repulsive nonnative contacts, resulting in an ideal smooth landscape. Non-Gō coarse-grained models reduce the parameter set as needed but do not include bias to any desired native structure. While non-Gō models have achieved limited success in protein coarse-graining, they can be combined with native structured-based potentials to create a balanced and powerful force field. Recent applications of such Gō-like models have yielded insight into complex folding mechanisms and conformational transitions in large macromolecules. The accuracy and usefulness of reduced representations are also revealed to be a function of the mathematical treatment of the intrinsic bonded topology.

Interplay between drying and stability of a TIM barrel protein: a combined simulation-experimental study

Recent molecular dynamics simulations have suggested important roles for nanoscale dewetting in the stability, function, and folding dynamics of proteins. Using a synergistic simulation-experimental approach on the αTS TIM barrel protein, we validated this hypothesis by revealing the occurrence of drying inside hydrophobic amino acid clusters and its manifestation in experimental measures of protein stability and structure. Cavities created within three clusters of branched aliphatic amino acids [isoleucine, leucine, and valine (ILV) clusters] were found to experience strong water density fluctuations or intermittent dewetting transitions in simulations. Individually substituting 10 residues in the large ILV cluster at the N-terminus with less hydrophobic alanines showed a weakening or diminishing effect on dewetting that depended on the site of the mutation. Our simulations also demonstrated that replacement of buried leucines with isosteric, polar asparagines enhanced the wetting of the N- and C-terminal clusters. The experimental results on the stability, secondary structure, and compactness of the native and intermediate states for the asparagine variants are consistent with the preferential drying of the large N-terminal cluster in the intermediate. By contrast, the region encompassing the small C-terminal cluster experiences only partial drying in the intermediate, and its structure and stability are unaffected by the asparagine substitution. Surprisingly, the structural distortions required to accommodate the replacement of leucine by asparagine in the N-terminal cluster revealed the existence of alternative stable folds in the native basin. This combined simulation-experimental study demonstrates the critical role of drying within hydrophobic ILV clusters in the folding and stability of the αTS TIM barrel.

Clusters of branched aliphatic side chains serve as cores of stability in the native state of the HisF TIM barrel protein

Imidazole-3-glycerol phosphate synthase is a heterodimeric allosteric enzyme that catalyzes consecutive reactions in imidazole biosynthesis through its HisF and HisH subunits. The unusually slow unfolding reaction of the isolated HisF TIM barrel domain from the thermophilic bacteria, Thermotoga maritima, enabled an NMR-based site-specific analysis of the main-chain hydrogen bonds that stabilize its native conformation. Very strong protection against exchange with solvent deuterium in the native state was found in a subset of buried positions in α-helices and pervasively in the underlying β-strands associated with a pair of large clusters of isoleucine, leucine and valine (ILV) side chains located in the α7(βα)8(βα)1-2 and α2(βα)3-6β7 segments of the (βα)8 barrel. The most densely packed region of the large cluster, α3(βα)4-6β7, correlates closely with the core of stability previously observed in computational, protein engineering and NMR dynamics studies, demonstrating a key role for this cluster in determining the thermodynamic and structural properties of the native state of HisF. When considered with the results of previous studies where ILV clusters were found to stabilize the hydrogen-bonded networks in folding intermediates for other TIM barrel proteins, it appears that clusters of branched aliphatic side chains can serve as cores of stability across the entire folding reaction coordinate of one of the most common motifs in biology.

Native-state heterogeneity of β(2)-microglobulin as revealed by kinetic folding and real-time NMR experiments

The kinetic folding of β(2)-microglobulin from the acid-denatured state was investigated by interrupted-unfolding and interrupted-refolding experiments using stopped-flow double-jump techniques. In the interrupted unfolding, we first unfolded the protein by a pH jump from pH7.5 to pH2.0, and the kinetic refolding assay was carried out by the reverse pH jump by monitoring tryptophan fluorescence. Similarly, in the interrupted refolding, we first refolded the protein by a pH jump from pH2.0 to pH7.5 and used a guanidine hydrochloride (GdnHCl) concentration jump as well as the reverse pH jump as unfolding assays. Based on these experiments, the folding is represented by a parallel-pathway model, in which the molecule with the correct Pro32 cis isomer refolds rapidly with a rate constant of 5-6 s(-1), while the molecule with the Pro32 trans isomer refolds more slowly (pH7.5 and 25°C). At the last step of folding, the native-like trans conformer produced on the latter pathway isomerizes very slowly (0.001-0.002 s(-1)) into the native cis conformer. In the GdnHCl-induced unfolding assays in the interrupted refolding, the native-like trans conformer unfolded remarkably faster than the native cis conformer, and the direct GdnHCl-induced unfolding was also biphasic, indicating that the native-like trans conformer is populated at a significant level under the native condition. The one-dimensional NMR and the real-time NMR experiments of refolding further indicated that the population of the trans conformer increases up to 7-9% under a more physiological condition (pH7.5 and 37°C).

Illuminating the off-pathway nature of the molten globule folding intermediate of an α-β parallel protein

Partially folded protein species transiently form during folding of most proteins. Often, these species are molten globules, which may be on- or off-pathway to the native state. Molten globules are ensembles of interconverting protein conformers that have a substantial amount of secondary structure, but lack virtually all tertiary side-chain packing characteristics of natively folded proteins. Due to solvent-exposed hydrophobic groups, molten globules are prone to aggregation, which can have detrimental effects on organisms. The molten globule observed during folding of the 179-residue apoflavodoxin from Azotobacter vinelandii is off-pathway, as it has to unfold before native protein can form. Here, we study folding of apoflavodoxin and characterize its molten globule using fluorescence spectroscopy and Förster Resonance Energy Transfer (FRET). Apoflavodoxin is site-specifically labeled with fluorescent donor and acceptor dyes, utilizing dye-inaccessibility of Cys69 in cofactor-bound protein. Donor (i.e., Alexa Fluor 488) is covalently attached to Cys69 in all apoflavodoxin variants used. Acceptor (i.e., Alexa Fluor 568) is coupled to Cys1, Cys131 and Cys178, respectively. Our FRET data show that apoflavodoxin's molten globule forms in a non-cooperative manner and that its N-terminal 69 residues fold last. In addition, striking conformational differences between molten globule and native protein are revealed, because the inter-label distances sampled in the 111-residue C-terminal segment of the molten globule are shorter than observed for native apoflavodoxin. Thus, FRET sheds light on the off-pathway nature of the molten globule during folding of an α-β parallel protein.

Folding mechanism of an extremely thermostable (βα)(8)-barrel enzyme: a high kinetic barrier protects the protein from denaturation

HisF, the cyclase subunit of imidazole glycerol phosphate synthase (ImGPS) from Thermotoga maritima, is an extremely thermostable (βα)(8)-barrel protein. We elucidated the unfolding and refolding mechanism of HisF. Its unfolding transition is reversible and adequately described by the two-state model, but 6 weeks is necessary to reach equilibrium (at 25 °C). During refolding, initially a burst-phase off-pathway intermediate is formed. The subsequent productive folding occurs in two kinetic phases with time constants of ~3 and ~20 s. They reflect a sequential process via an on-pathway intermediate, as revealed by stopped-flow double-mixing experiments. The final step leads to native HisF, which associates with the glutaminase subunit HisH to form the functional ImGPS complex. The conversion of the on-pathway intermediate to the native protein results in a 10(6)-fold increase of the time constant for unfolding from 89 ms to 35 h (at 4.0 M GdmCl) and thus establishes a high energy barrier to denaturation. We conclude that the extra stability of HisF is used for kinetic protection against unfolding. In its refolding mechanism, HisF resembles other (βα)(8)-barrel proteins.

Environment-dependent long-range structural distortion in a temperature-sensitive point mutant

Extensive environment-dependent rearrangement of the helix-turn-helix DNA recognition region and adjacent L-tryptophan binding pocket is reported in the crystal structure of dimeric E. coli trp aporepressor with point mutation Leu75Phe. In one of two subunits, the eight residues immediately C-terminal to the mutation are shifted forward in helical register by three positions, and the five following residues form an extrahelical loop accommodating the register shift. In contrast, the second subunit has wildtype-like conformation, as do both subunits in an isomorphous wildtype control structure. Treated together as an ensemble pair, the distorted and wildtype-like conformations of the mutant apoprotein agree more fully than either conformation alone with previously reported NOE measurements, and account more completely for its diverse biochemical and biophysical properties. The register-shifted segment Ile79-Ala80-Thr81-Ile82-Thr83 is helical in both conformations despite low helical propensity, suggesting an important structural role for the steric constraints imposed by β-branched residues in helical conformation.

Osmolyte-induced folding of an intrinsically disordered protein: folding mechanism in the absence of ligand

Understanding the interconversion between thermodynamically distinguishable states present in a protein folding pathway provides not only the kinetics and energetics of protein folding but also insights into the functional roles of these states in biological systems. The protein component of the bacterial RNase P holoenzyme from Bacillus subtilis (P protein) was previously shown to be unfolded in the absence of its cognate RNA or other anionic ligands. P protein was used in this study as a model system to explore general features of intrinsically disordered protein (IDP) folding mechanisms. The use of trimethylamine N-oxide (TMAO), an osmolyte that stabilizes the unliganded folded form of the protein, enabled us to study the folding process of P protein in the absence of ligand. Transient stopped-flow kinetic traces at various final TMAO concentrations exhibited multiphasic kinetics. Equilibrium "cotitration" experiments were performed using both TMAO and urea during the titration to produce a urea-TMAO titration surface of P protein. Both kinetic and equilibrium studies show evidence of a previously undetected intermediate state in the P protein folding process. The intermediate state is significantly populated, and the folding rate constants are relatively slow compared to those of intrinsically folded proteins similar in size and topology. The experiments and analysis described serve as a useful example for mechanistic folding studies of other IDPs.

Topological frustration in beta alpha-repeat proteins: sequence diversity modulates the conserved folding mechanisms of alpha/beta/alpha sandwich proteins

The thermodynamic hypothesis of Anfinsen postulates that structures and stabilities of globular proteins are determined by their amino acid sequences. Chain topology, however, is known to influence the folding reaction, in that motifs with a preponderance of local interactions typically fold more rapidly than those with a larger fraction of nonlocal interactions. Together, the topology and sequence can modulate the energy landscape and influence the rate at which the protein folds to the native conformation. To explore the relationship of sequence and topology in the folding of beta alpha-repeat proteins, which are dominated by local interactions, we performed a combined experimental and simulation analysis on two members of the flavodoxin-like, alpha/beta/alpha sandwich fold. Spo0F and the N-terminal receiver domain of NtrC (NT-NtrC) have similar topologies but low sequence identity, enabling a test of the effects of sequence on folding. Experimental results demonstrated that both response-regulator proteins fold via parallel channels through highly structured submillisecond intermediates before accessing their cis prolyl peptide bond-containing native conformations. Global analysis of the experimental results preferentially places these intermediates off the productive folding pathway. Sequence-sensitive Gō-model simulations conclude that frustration in the folding in Spo0F, corresponding to the appearance of the off-pathway intermediate, reflects competition for intra-subdomain van der Waals contacts between its N- and C-terminal subdomains. The extent of transient, premature structure appears to correlate with the number of isoleucine, leucine, and valine (ILV) side chains that form a large sequence-local cluster involving the central beta-sheet and helices alpha2, alpha 3, and alpha 4. The failure to detect the off-pathway species in the simulations of NT-NtrC may reflect the reduced number of ILV side chains in its corresponding hydrophobic cluster. The location of the hydrophobic clusters in the structure may also be related to the differing functional properties of these response regulators. Comparison with the results of previous experimental and simulation analyses on the homologous CheY argues that prematurely folded unproductive intermediates are a common property of the beta alpha-repeat motif.

Interrupted hydrogen/deuterium exchange reveals the stable core of the remarkably helical molten globule of alpha-beta parallel protein flavodoxin

Kinetic intermediates that appear early during protein folding often resemble the relatively stable molten globule intermediates formed by several proteins under mildly denaturing conditions. Molten globules have a substantial amount of secondary structure but lack virtually all tertiary side-chain packing characteristics of natively folded proteins. Due to exposed hydrophobic groups, molten globules are prone to aggregation, which can have detrimental effects on organisms. The molten globule that is observed during folding of alpha-beta parallel flavodoxin from Azotobacter vinelandii is a remarkably non-native species. This folding intermediate is helical and contains no beta-sheet and is kinetically off-pathway to the native state. It can be trapped under native-like conditions by substituting residue Phe(44) for Tyr(44). To characterize this species at the residue level, in this study, use is made of interrupted hydrogen/deuterium exchange detected by NMR spectroscopy. In the molten globule of flavodoxin, the helical region comprising residues Leu(110)-Val(125) is shown to be better protected against exchange than the other ordered parts of the folding intermediate. This helical region is better buried than the other helices, causing its context-dependent stabilization against unfolding. Residues Leu(110)-Val(125) thus form the stable core of the helical molten globule of alpha-beta parallel flavodoxin, which is almost entirely structured. Non-native docking of helices in the molten globule of flavodoxin prevents formation of the parallel beta-sheet of native flavodoxin. Hence, to produce native alpha-beta parallel protein molecules, the off-pathway species needs to unfold.

Residue-specific analysis of frustration in the folding landscape of repeat beta/alpha protein apoflavodoxin

Flavodoxin adopts the common repeat beta/alpha topology and folds in a complex kinetic reaction with intermediates. To better understand this reaction, we analyzed a set of Desulfovibrio desulfuricans apoflavodoxin variants with point mutations in most secondary structure elements by in vitro and in silico methods. By equilibrium unfolding experiments, we first revealed how different secondary structure elements contribute to overall protein resistance to heat and urea. Next, using stopped-flow mixing coupled with far-UV circular dichroism, we probed how individual residues affect the amount of structure formed in the experimentally detected burst-phase intermediate. Together with in silico folding route analysis of the same point-mutated variants and computation of growth in nucleation size during early folding, computer simulations suggested the presence of two competing folding nuclei at opposite sides of the central beta-strand 3 (i.e., at beta-strands 1 and 4), which cause early topological frustration (i.e., misfolding) in the folding landscape. Particularly, the extent of heterogeneity in folding nuclei growth correlates with the in vitro burst-phase circular dichroism amplitude. In addition, phi-value analysis (in vitro and in silico) of the overall folding barrier to apoflavodoxin's native state revealed that native-like interactions in most of the beta-strands must form in transition state. Our study reveals that an imbalanced competition between the two sides of apoflavodoxin's central beta-sheet directs initial misfolding, while proper alignment on both sides of beta-strand 3 is necessary for productive folding.

Modular control of cross-oligomerization: analysis of superstabilized Hsp90 homodimers in vivo

Homo-oligomeric proteins fulfill numerous functions in all cells. The ability to co-express subunits of these proteins that preferentially self-assemble without cross-oligomerizing provides for controlled experiments to analyze the function of mutant homo-oligomers in vivo. Hsp90 is a dimeric chaperone involved in the maturation of many kinases and steroid hormone receptors. We observed that co-expression of different Hsp90 subunits in Saccharomyces cerevisiae caused unpredictable synthetic growth defects due to cross-dimerization. We engineered superstabilized Hsp90 dimers that resisted cross-dimerization with endogenous Hsp90 and alleviated the synthetic growth defect. Superstabilized Hsp90 dimers supported robust growth of S. cerevisiae, indicating that dissociation of Hsp90 dimers could be hindered without compromising essential function. We utilized superstabilized dimers to analyze the activity of ATPase mutant homodimers in a temperature-sensitive yeast background where elevated temperature inactivated all other Hsp90 species. We found that ATP binding and hydrolysis by Hsp90 are both required for the efficient maturation of glucocorticoid receptor and v-Src, confirming the critical role of ATP hydrolysis in the maturation of steroid hormone receptors and kinases in vivo.

Non-native hydrophobic interactions detected in unfolded apoflavodoxin by paramagnetic relaxation enhancement

Transient structures in unfolded proteins are important in elucidating the molecular details of initiation of protein folding. Recently, native and non-native secondary structure have been discovered in unfolded A. vinelandii flavodoxin. These structured elements transiently interact and subsequently form the ordered core of an off-pathway folding intermediate, which is extensively formed during folding of this α–β parallel protein. Here, site-directed spin-labelling and paramagnetic relaxation enhancement are used to investigate long-range interactions in unfolded apoflavodoxin. For this purpose, glutamine-48, which resides in a non-native α-helix of unfolded apoflavodoxin, is replaced by cysteine. This replacement enables covalent attachment of nitroxide spin-labels MTSL and CMTSL. Substitution of Gln-48 by Cys-48 destabilises native apoflavodoxin and reduces flexibility of the ordered regions in unfolded apoflavodoxin in 3.4 M GuHCl, because of increased hydrophobic interactions in the unfolded protein. Here, we report that in the study of the conformational and dynamic properties of unfolded proteins interpretation of spin-label data can be complicated. The covalently attached spin-label to Cys-48 (or Cys-69 of wild-type apoflavodoxin) perturbs the unfolded protein, because hydrophobic interactions occur between the label and hydrophobic patches of unfolded apoflavodoxin. Concomitant hydrophobic free energy changes of the unfolded protein (and possibly of the off-pathway intermediate) reduce the stability of native spin-labelled protein against unfolding. In addition, attachment of MTSL or CMTSL to Cys-48 induces the presence of distinct states in unfolded apoflavodoxin. Despite these difficulties, the spin-label data obtained here show that non-native contacts exist between transiently ordered structured elements in unfolded apoflavodoxin.

Topological switching between an alpha-beta parallel protein and a remarkably helical molten globule

Partially folded protein species transiently exist during folding of most proteins. Often these species are molten globules, which may be on- or off-pathway to native protein. Molten globules have a substantial amount of secondary structure but lack virtually all the tertiary side-chain packing characteristic of natively folded proteins. These ensembles of interconverting conformers are prone to aggregation and potentially play a role in numerous devastating pathologies, and thus attract considerable attention. The molten globule that is observed during folding of apoflavodoxin from Azotobacter vinelandii is off-pathway, as it has to unfold before native protein can be formed. Here we report that this species can be trapped under nativelike conditions by substituting amino acid residue F44 by Y44, allowing spectroscopic characterization of its conformation. Whereas native apoflavodoxin contains a parallel beta-sheet surrounded by alpha-helices (i.e., the flavodoxin-like or alpha-beta parallel topology), it is shown that the molten globule has a totally different topology: it is helical and contains no beta-sheet. The presence of this remarkably nonnative species shows that single polypeptide sequences can code for distinct folds that swap upon changing conditions. Topological switching between unrelated protein structures is likely a general phenomenon in the protein structure universe.

Noncooperative Formation of the off-pathway molten globule during folding of the alpha-beta parallel protein apoflavodoxin

During folding of many proteins, molten globules are formed. These partially folded forms of proteins have a substantial amount of secondary structure but lack virtually all tertiary side-chain packing characteristic of native structures. Molten globules are ensembles of interconverting conformers and are prone to aggregation, which can have detrimental effects on organisms. Consequently, molten globules attract considerable attention. The molten globule that is observed during folding of flavodoxin from Azotobacter vinelandii is a kinetically off-pathway species, as it has to unfold before the native state of the protein can be formed. This intermediate contains helices and can be populated at equilibrium using guanidinium hydrochloride as denaturant, allowing the use of NMR spectroscopy to follow molten globule formation at the residue level. Here, we track changes in chemical shifts of backbone amides, as well as disappearance of resonances of unfolded apoflavodoxin, upon decreasing denaturant concentration. Analysis of the data shows that structure formation within virtually all parts of the unfolded protein precedes folding to the molten globule state. This folding transition is noncooperative and involves a series of distinct transitions. Four structured elements in unfolded apoflavodoxin transiently interact and subsequently form the ordered core of the molten globule. Although hydrophobic, tryptophan side chains are not involved in the latter process. This ordered core is gradually extended upon decreasing denaturant concentration, but part of apoflavodoxin's molten globule remains random coil in the denaturant range investigated. The results presented here, together with those reported on the molten globule of alpha-lactalbumin, show that helical molten globules apparently fold in a noncooperative manner.

Insights from coarse-grained Gō models for protein folding and dynamics

Exploring the landscape of large scale conformational changes such as protein folding at atomistic detail poses a considerable computational challenge. Coarse-grained representations of the peptide chain have therefore been developed and over the last decade have proved extremely valuable. These include topology-based Gō models, which constitute a smooth and funnel-like approximation to the folding landscape. We review the many variations of the Gō model that have been employed to yield insight into folding mechanisms. Their success has been interpreted as a consequence of the dominant role of the native topology in folding. The role of local contact density in determining protein dynamics is also discussed and is used to explain the ability of Gō-like models to capture sequence effects in folding and elucidate conformational transitions.

Subdomain competition, cooperativity, and topological frustration in the folding of CheY

The folding of multidomain proteins often proceeds in a hierarchical fashion with individual domains folding independent of one another. A large single-domain protein, however, can consist of multiple modules whose folding may be autonomous or interdependent in ways that are unclear. We used coarse-grained simulations to explore the folding landscape of the two-subdomain bacterial response regulator CheY. Thermodynamic and kinetic characterization shows the landscape to be highly analogous to the four-state landscape reported for another two-subdomain protein, T4 lysozyme. An on-pathway intermediate structured in the more stable nucleating subdomain was observed, as were transient states frustrated in off-pathway contacts prematurely structured in the weaker subdomain. Local unfolding, or backtracking, was observed in the frustrated state before the native conformation could be reached. Nonproductive frustration was attributable to competition for van der Waals contacts between the two subdomains. In an accompanying article, stopped-flow kinetic measurements support an off-pathway burst-phase intermediate, seemingly consistent with our prediction of early frustration in the folding landscape of CheY. Comparison of the folding mechanisms for CheY, T4 lysozyme, and interleukin-1 beta leads us to postulate that subdomain competition is a general feature of large single-domain proteins with multiple folding modules.