Asymmetric effect of domain interactions on the kinetics of folding in yeast phosphoglycerate kinase

Protein Folding Stability and Kinetics in Alginate Hydrogels

Proteins are commonly encapsulated in alginate gels for drug delivery and tissue-engineering applications. However, there is limited knowledge of how encapsulation impacts intrinsic protein properties such as folding stability or unfolding kinetics. Here, we use fast relaxation imaging (FReI) to image protein unfolding in situ in alginate hydrogels after applying a temperature jump. Based on changes in the Förster resonance energy transfer (FRET) response of FRET-labeled phosphoglycerate kinase (PGK), we report the quantitative impact of multiple alginate hydrogel concentrations on protein stability and folding dynamics. The gels stabilize PGK by increasing its melting temperature up to 18.4 °C, and the stabilization follows a nonmonotonic dependence on the alginate density. In situ kinetic measurements also reveal that PGK deviates more from two-state folding behavior in denser gels and that the gel decreases the unfolding rate and accelerates the folding rate of PGK, compared to buffer. Phi-value analysis suggests that the folding transition state of an encapsulated protein is structurally similar to that of folded protein. This work reveals both beneficial and negative impacts of gel encapsulation on protein folding, as well as potential mechanisms contributing to altered stability.

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.

Protein folding in vitro and in the cell: From a solitary journey to a team effort

Correct protein folding is essential for the health and function of living organisms. Yet, it is not well understood how unfolded proteins reach their native state and avoid aggregation, especially within the cellular milieu. Some proteins, especially small, single-domain and apparent two-state folders, successfully attain their native state upon dilution from denaturant. Yet, many more proteins undergo misfolding and aggregation during this process, in a concentration-dependent fashion. Once formed, native and aggregated states are often kinetically trapped relative to each other. Hence, the early stages of protein life are absolutely critical for proper kinetic channeling to the folded state and for long-term solubility and function. This review summarizes current knowledge on protein folding/aggregation mechanisms in buffered solution and within the bacterial cell, highlighting early stages. Remarkably, teamwork between nascent chain, ribosome, trigger factor and Hsp70 molecular chaperones enables all proteins to overcome aggregation propensities and reach a long-lived bioactive state.

Stabilization and Kinetics of an Adsorbed Protein Depends on the Poly(N-isopropylacrylamide) Grafting Density

The solubility transition at the lower critical solution temperature (LCST, 32 °C) of poly(N-isopropylacrylamide) (PNIPAM) is widely used as a thermal switch to rapidly and reversibly capture and release proteins and cells. It is generally assumed that proteins adsorbed to PNIPAM above the LCST are unaffected by polymer interactions. Here we show that the folding stability of the enzyme phosphoglycerate kinase (PGK) is increased by interactions with end-grafted PNIPAM films above the LCST. We systematically compare two protein mutants with different stabilities. The stabilization mirrors the degree of protein adsorption under grafting conditions studied previously. Maximum stabilization occurs when proteins adsorb to low density, collapsed polymer "mushrooms". In the denser polymer "brush" regime, protein stabilization decreases back to a value indistinguishable from the bulk solution, consistent with low protein adsorption on dense, collapsed brushes. The temperature-dependent kinetics measured by Fast Relaxation Imaging reveals that PNIPAM does not affect the overall folding/unfolding mechanism. Based on the different stabilizations of two mutants and the relaxation kinetics, we hypothesize that the polymer acts mainly by increasing the conformational entropy of the folded protein by interacting with the protein surface and less by crowding the unfolded state of PGK.

Applying text mining to identify relevant literature in food science: Cold denaturation as a case study

In a research environment characterized by the five V's of big data, volume, velocity, variety, value, and veracity, the need to develop tools that quickly screen a large number of publications into relevant work is an increasing area of concern, and the data-rich food industry is no exception. Here, a combination of latent Dirichlet allocation and food keyword searches were employed to analyze and filter a dataset of 6102 publications about cold denaturation. After using the Python toolkit generated in this work, the approach yielded 22 topics that provide background and insight on the direction of research in this field, as well as identified the publications in this dataset which are most pertinent to the food industry with precision and recall of 0.419 and 0.949, respectively. Precision is related to the relevance of a paper in the filtered dataset and the recall represents papers which were not identified in the screening method. Lastly, gaps in the literature based on keyword trends are identified to improve the knowledge base of cold denaturation as it relates to the food industry. This approach is generalizable to any similarly organized dataset, and the code is available upon request. Practical Application: A common problem in research is that when you are an expert in one field, learning about another field is difficult, because you may lack the vocabulary and background needed to read cutting edge literature from a new discipline. The Python toolkit developed in this research can be applied by any researcher that is new to a field to identify what the key literature is, what topics they should familiarize themselves with, and what the current trends are in the field. Using this structure, researchers can greatly speed up how they identify new areas to research and find new projects.

Probing the Effects of Local Frustration in the Folding of a Multidomain Protein

Our current knowledge of protein folding is primarily based on experimental data obtained from isolated domains. In fact, because of their complexity, multidomain proteins have been elusive to the experimental characterization. Thus, the folding of a domain in isolation is generally assumed to resemble what should be observed for more complex structural architectures. Here we compared the folding mechanism of a protein domain in isolation and in the context of its supramodular multidomain structure. By carrying out an extensive mutational analysis we illustrate that while the early events of folding are malleable and influenced by the absence/presence of the neighboring structures, the late events appear to be more robust. These effects may be explained by analyzing the local frustration patterns of the domain, providing critical support for the funneled energy landscape theory of protein folding, and highlighting the role of protein frustration in sculpting the early events of the reaction.

Protein folding and surface interaction phase diagrams in vitro and in cells

Protein stability is subject to environmental perturbations such as pressure and crowding, as well as sticking to other macromolecules and quinary structure. Thus, the environment inside and outside the cell plays a key role in how proteins fold, interact, and function on the scale from a few molecules to macroscopic ensembles. This review discusses three aspects of protein phase diagrams: first, the relevance of phase diagrams to protein folding and function in vitro and in cells; next, how the evolution of protein surfaces impacts on interaction phase diagrams; and finally, how phase separation plays a role on much larger length-scales than individual proteins or oligomers, when liquid phase-separated regions form to assist protein function and cell homeostasis.

Investigating the trade-off between folding and function in a multidomain Y-family DNA polymerase

The way in which multidomain proteins fold has been a puzzling question for decades. Until now, the mechanisms and functions of domain interactions involved in multidomain protein folding have been obscure. Here, we develop structure-based models to investigate the folding and DNA-binding processes of the multidomain Y-family DNA polymerase IV (DPO4). We uncover shifts in the folding mechanism among ordered domain-wise folding, backtracking folding, and cooperative folding, modulated by interdomain interactions. These lead to ‘U-shaped’ DPO4 folding kinetics. We characterize the effects of interdomain flexibility on the promotion of DPO4–DNA (un)binding, which probably contributes to the ability of DPO4 to bypass DNA lesions, which is a known biological role of Y-family polymerases. We suggest that the native topology of DPO4 leads to a trade-off between fast, stable folding and tight functional DNA binding. Our approach provides an effective way to quantitatively correlate the roles of protein interactions in conformational dynamics at the multidomain level.

Mapping Multiple Distances in a Multidomain Protein for the Identification of Folding Intermediates

The investigation and understanding of the folding mechanism of multidomain proteins is still a challenge in structural biology. The use of single-molecule Förster resonance energy transfer offers a unique tool to map conformational changes within the protein structure. Here, we present a study following denaturant-induced unfolding transitions of yeast phosphoglycerate kinase by mapping several inter- and intradomain distances of this two-domain protein, exhibiting a quite heterogeneous behavior. On the one hand, the development of the interdomain distance during the unfolding transition suggests a classical two-state unfolding behavior. On the other hand, the behavior of some intradomain distances indicates the formation of a compact and transient molten globule intermediate state. Furthermore, different intradomain distances measured within the same domain show pronounced differences in their unfolding behavior, underlining the fact that the choice of dye attachment positions within the polypeptide chain has a substantial impact on which unfolding properties are observed by single-molecule Förster resonance energy transfer measurements. Our results suggest that, to fully characterize the complex folding and unfolding mechanism of multidomain proteins, it is necessary to monitor multiple intra- and interdomain distances because a single reporter can lead to a misleading, partial, or oversimplified interpretation.

Unraveling the Mechanical Unfolding Pathways of a Multidomain Protein: Phosphoglycerate Kinase

Phosphoglycerate kinase (PGK) is a highly conserved enzyme that is crucial for glycolysis. PGK is a monomeric protein composed of two similar domains and has been the focus of many studies for investigating interdomain interactions within the native state and during folding. Previous studies used traditional biophysical methods (such as circular dichroism, tryptophan fluorescence, and NMR) to measure signals over a large ensemble of molecules, which made it difficult to observe transient changes in stability or structure during unfolding and refolding of single molecules. Here, we unfold single molecules of PGK using atomic force spectroscopy and steered molecular dynamic computer simulations to examine the conformational dynamics of PGK during its unfolding process. Our results show that after the initial forced separation of its domains, yeast PGK (yPGK) does not follow a single mechanical unfolding pathway; instead, it stochastically follows two distinct pathways: unfolding from the N-terminal domain or unfolding from the C-terminal domain. The truncated yPGK N-terminal domain unfolds via a transient intermediate, whereas the structurally similar isolated C-terminal domain has no detectable intermediates throughout its mechanical unfolding process. The N-terminal domain in the full-length yPGK displays a strong unfolding intermediate 13% of the time, whereas the truncated domain (yPGKNT) transitions through the intermediate 81% of the time. This effect indicates that the mechanical properties of yPGK cannot be simply deduced from the mechanical properties of its constituents. We also find that Escherichia coli PGK is significantly less mechanically stable as compared to yPGK, contrary to bulk unfolding measurements. Our results support the growing body of observations that the folding behavior of multidomain proteins is difficult to predict based solely on the studies of isolated domains.

Same but not alike: Structure, flexibility and energetics of domains in multi-domain proteins are influenced by the presence of other domains

The majority of the proteins encoded in the genomes of eukaryotes contain more than one domain. Reasons for high prevalence of multi-domain proteins in various organisms have been attributed to higher stability and functional and folding advantages over single-domain proteins. Despite these advantages, many proteins are composed of only one domain while their homologous domains are part of multi-domain proteins. In the study presented here, differences in the properties of protein domains in single-domain and multi-domain systems and their influence on functions are discussed. We studied 20 pairs of identical protein domains, which were crystallized in two forms (a) tethered to other proteins domains and (b) tethered to fewer protein domains than (a) or not tethered to any protein domain. Results suggest that tethering of domains in multi-domain proteins influences the structural, dynamic and energetic properties of the constituent protein domains. 50% of the protein domain pairs show significant structural deviations while 90% of the protein domain pairs show differences in dynamics and 12% of the residues show differences in the energetics. To gain further insights on the influence of tethering on the function of the domains, 4 pairs of homologous protein domains, where one of them is a full-length single-domain protein and the other protein domain is a part of a multi-domain protein, were studied. Analyses showed that identical and structurally equivalent functional residues show differential dynamics in homologous protein domains; though comparable dynamics between in-silico generated chimera protein and multi-domain proteins were observed. From these observations, the differences observed in the functions of homologous proteins could be attributed to the presence of tethered domain. Overall, we conclude that tethered domains in multi-domain proteins not only provide stability or folding advantages but also influence pathways resulting in differences in function or regulatory properties.

High prevalence of multi-domain proteins in proteomes has been attributed to higher stability and functional and folding advantages of the multi-domain proteins. Influence of tethering of domains on the overall properties of proteins has been well studied but its influence on the properties of the constituent domains is largely unaddressed. Here, we investigate the influence of tethering of domains in multi-domain proteins on the structural, dynamics and energetics properties of the constituent domains and its implications on the functions of proteins. To this end, comparative analyses were carried out for identical protein domains crystallized in tethered and untethered forms. Also, comparative analyses of single-domain proteins and their homologous multi-domain proteins were performed. The analyses suggest that tethering influences the structural, dynamic and energetic properties of constituent protein domains. Our observations hint at regulation of protein domains by tethered domains in multi-domain systems, which may manifest at the differential function observed between single-domain and homologous multi-domain proteins.

On-chip measurements of protein unfolding from direct observations of micron-scale diffusion

The unfolding process of BSA in solution as a function of pH was studied by microfluidic diffusional sizing device.

Investigations of protein folding, unfolding and stability are critical for the understanding of the molecular basis of biological structure and function. We describe here a microfluidic approach to probe the unfolding of unlabelled protein molecules in microliter volumes. We achieve this objective through the use of a microfluidic platform, which allows the changes in molecular diffusivity upon folding and unfolding to be detected directly. We illustrate this approach by monitoring the unfolding of bovine serum albumin in solution as a function of pH. These results show the viability of probing protein stability on chip in small volumes.

Protein Assembly and Building Blocks: Beyond the Limits of the LEGO Brick Metaphor

Proteins, like other biomolecules, have a modular and hierarchical structure. Various building blocks are used to construct proteins of high structural complexity and diverse functionality. In multidomain proteins, for example, domains are fused to each other in different combinations to achieve different functions. Although the LEGO brick metaphor is justified as a means of simplifying the complexity of three-dimensional protein structures, several fundamental properties (such as allostery or the induced-fit mechanism) make deviation from it necessary to respect the plasticity, softness, and cross-talk that are essential to protein function. In this work, we illustrate recently reported protein behavior in multidomain proteins that deviates from the LEGO brick analogy. While earlier studies showed that a protein domain is often unaffected by being fused to another domain or becomes more stable following the formation of a new interface between the tethered domains, destabilization due to tethering has been reported for several systems. We illustrate that tethering may sometimes result in a multidomain protein behaving as "less than the sum of its parts". We survey these cases for which structure additivity does not guarantee thermodynamic additivity. Protein destabilization due to fusion to other domains may be linked in some cases to biological function and should be taken into account when designing large assemblies.

Identification of cancer fusion drivers using network fusion centrality

SUMMARY

Gene fusions are being discovered at an increasing rate using massively parallel sequencing technologies. Prioritization of cancer fusion drivers for validation cannot be performed using traditional single-gene based methods because fusions involve portions of two partner genes. To address this problem, we propose a novel network analysis method called fusion centrality that is specifically tailored for prioritizing gene fusions. We first propose a domain-based fusion model built on the theory of exon/domain shuffling. The model leads to a hypothesis that a fusion is more likely to be an oncogenic driver if its partner genes act like hubs in a network because the fusion mutation can deregulate normal functions of many other genes and their pathways. The hypothesis is supported by the observation that for most known cancer fusion genes, at least one of the fusion partners appears to be a hub in a network, and even for many fusions both partners appear to be hubs. Based on this model, we construct fusion centrality, a multi-gene-based network metric, and use it to score fusion drivers. We show that the fusion centrality outperforms other single gene-based methods. Specifically, the method successfully predicts most of 38 newly discovered fusions that had validated oncogenic importance. To our best knowledge, this is the first network-based approach for identifying fusion drivers.

AVAILABILITY

Matlab code implementing the fusion centrality method is available upon request from the corresponding authors.

Role of domain interactions in the collective motion of phosphoglycerate kinase

Protein function is governed by the underlying conformational dynamics of the molecule. The experimental and theoretical work leading to contemporary understanding of enzyme dynamics was mostly restricted to the large-scale movements of single-domain proteins. Collective movements resulting from a regulatory interplay between protein domains is often crucial for enzymatic activity. It is not clear, however, how our knowledge could be extended to describe collective near-equilibrium motions of multidomain enzymes. We examined the effect of domain interactions on the low temperature near equilibrium dynamics of the native state, using phosphoglycerate kinase as model protein. We measured thermal activation of tryptophan phosphorescence quenching to explore millisecond-range protein motions. The two protein domains of phosphoglycerate kinase correspond to two dynamic units, but interdomain interactions link the motion of the two domains. The effect of the interdomain interactions on the activation of motions in the individual domains is asymmetric. As the temperature of the frozen protein is increased from the cryogenic, motions of the N domain are activated first. This is a partial activation, however, and the full dynamics of the domain becomes activated only after the activation of the C domain.

Comparing the folding and misfolding energy landscapes of phosphoglycerate kinase

Partitioning of polypeptides between protein folding and amyloid formation is of outstanding pathophysiological importance. Using yeast phosphoglycerate kinase as model, here we identify the features of the energy landscape that decide the fate of the protein: folding or amyloidogenesis. Structure formation was initiated from the acid-unfolded state, and monitored by fluorescence from 10 ms to 20 days. Solvent conditions were gradually shifted between folding and amyloidogenesis, and the properties of the energy landscape governing structure formation were reconstructed. A gradual transition of the energy landscape between folding and amyloid formation was observed. In the early steps of both folding and misfolding, the protein searches through a hierarchically structured energy landscape to form a molten globule in a few seconds. Depending on the conditions, this intermediate either folds to the native state in a few minutes, or forms amyloid fibers in several days. As conditions are changed from folding to misfolding, the barrier separating the molten globule and native states increases, although the barrier to the amyloid does not change. In the meantime, the native state also becomes more unstable and the amyloid more stable. We conclude that the lower region of the energy landscape determines the final protein structure.

Protein stability and folding kinetics in the nucleus and endoplasmic reticulum of eucaryotic cells

We measure the stability and folding relaxation rate of phosphoglycerate kinase (PGK) Förster resonance energy transfer (FRET) constructs localized in the nucleus or in the endoplasmic reticulum (ER) of eukaryotic cells. PGK has a more compact native state in the cellular compartments than in aqueous solution. Its native FRET signature is similar to that previously observed in a carbohydrate-crowding matrix, consistent with crowding being responsible for the compact native state of PGK in the cell. PGK folds through multiple states in vitro, but its folding kinetics is more two-state-like in the ER, so the folding mechanism can be modified by intracellular compartments. The nucleus increases PGK stability and folding rate over the cytoplasm and ER, even though the density of crowders in the nucleus is no greater than in the ER or cytoplasm. Nuclear folding kinetics (and to a lesser extent, thermodynamics) vary less from cell to cell than in the cytoplasm or ER, indicating a more homogeneous crowding and chemical environment in the nucleus.

What lessons can be learned from studying the folding of homologous proteins?

The studies of the folding of structurally related proteins have proved to be a very important tool for investigating protein folding. Here we review some of the insights that have been gained from such studies. Our highlighted studies show just how such an investigation should be designed and emphasise the importance of the synergy between experiment and theory. We also stress the importance of choosing the right system carefully, exploiting the excellent structural and sequence databases at our disposal.

Studying the folding of multidomain proteins

There have been relatively few detailed comprehensive studies of the folding of protein domains (or modules) in the context of their natural covalently linked neighbors. This is despite the fact that a significant proportion of the proteome consists of multidomain proteins. In this review we highlight some key experimental investigations of the folding of multidomain proteins to draw attention to the difficulties that can arise in analyzing such systems. The evidence suggests that interdomain interactions can significantly affect stability, folding, and unfolding rates. However, preliminary studies suggest that folding pathways are unaffected-to this extent domains can be truly considered to be independent folding units. Nonetheless, it is clear that interactions between domains cannot be ignored, in particular when considering the effects of mutations.

Predicting protein folding pathways at the mesoscopic level based on native interactions between secondary structure elements

Background

Since experimental determination of protein folding pathways remains difficult, computational techniques are often used to simulate protein folding. Most current techniques to predict protein folding pathways are computationally intensive and are suitable only for small proteins.

Results

By assuming that the native structure of a protein is known and representing each intermediate conformation as a collection of fully folded structures in which each of them contains a set of interacting secondary structure elements, we show that it is possible to significantly reduce the conformation space while still being able to predict the most energetically favorable folding pathway of large proteins with hundreds of residues at the mesoscopic level, including the pig muscle phosphoglycerate kinase with 416 residues. The model is detailed enough to distinguish between different folding pathways of structurally very similar proteins, including the streptococcal protein G and the peptostreptococcal protein L. The model is also able to recognize the differences between the folding pathways of protein G and its two structurally similar variants NuG1 and NuG2, which are even harder to distinguish. We show that this strategy can produce accurate predictions on many other proteins with experimentally determined intermediate folding states.

Conclusion

Our technique is efficient enough to predict folding pathways for both large and small proteins at the mesoscopic level. Such a strategy is often the only feasible choice for large proteins. A software program implementing this strategy (SSFold) is available at .

The folding pathway of a single domain in a multidomain protein is not affected by its neighbouring domain

Domains are the structural, functional, and evolutionary components of proteins. Most folding studies to date have concentrated on the folding of single domains, but more than 70% of human proteins contain more than one domain, and interdomain interactions can affect both the stability and the folding kinetics. Whether the folding pathway is altered by interdomain interactions is not yet known. Here we investigated the effect of a folded neighbouring domain on the folding pathway of spectrin R16 (the 16th α-helical repeat from chicken brain α-spectrin) by using the two-domain construct R1516. The R16 folds faster and unfolds more slowly in the presence of its folded neighbour R15 (the 15th α-helical repeat from chicken brain α-spectrin). An extensive Φ-value analysis of the R16 domain in R1516 was completed to compare the transition state of the R16 domain alone with that of the R16 domain in a multidomain construct. The results indicate that the folding pathways are the same. This result validates the current approach of breaking up larger proteins into domains for the study of protein folding pathways.

Distinguishing specific and nonspecific interdomain interactions in multidomain proteins

Multidomain proteins account for over two-thirds of the eukaryotic genome. Although there have been extensive studies into the biophysical properties of isolated domains, few have investigated how the domains interact. Spectrin is a well-characterized multidomain protein with domains linked in tandem array by contiguous helices. Several of these domains have been shown to be stabilized by their neighbors. Until now, this stabilization has been attributed to specific interactions between the natural neighbors, however we have recently observed that nonnatural neighboring domains can also induce a significant amount of stabilization. Here we investigate this nonnative stabilizing effect. We created spectrin-titin domain pairs of both spectrin R16 and R17 with a single titin I27 domain at either the N- or the C-terminus and found that spectrin domains are significantly stabilized, through slowed unfolding, by nonnative interactions at the C-terminus only. Of particular importance, we show that specific interactions between natural folded neighbors at either terminus confer even greater stability by additionally increasing the folding rate constants. We demonstrate that it is possible to distinguish between natural stabilizing interactions and nonspecific stabilizing effects through examination of the kinetics of well chosen mutant proteins. This work adds to the complexity of studying multidomain proteins.

Comparison of Proteolytic Susceptibility in Phosphoglycerate Kinases from Yeast and E. coli: Modulation of Conformational Ensembles Without Altering Structure or Stability

Escherichia coli phosphoglycerate kinase (PGK) is resistant to proteolytic cleavage while the yeast homolog from Saccharomyces cerevisiae is not. We have explored the biophysical basis of this surprising difference. The sequences of these homologs are 39% identical and 56% similar. Determination of the crystal structure for the E. coli protein and comparison to the previously solved yeast structure reveals that the two proteins have extremely similar tertiary structures, and their global stabilities determined by equilibrium denaturation are also very similar. The extrapolated unfolding rate of E. coli PGK is, however, 10(5) slower than that of the yeast homolog. This surprisingly large difference in unfolding rates appears to arise from a divergence in the extent of cooperativity between the two structural domains (the N and C-domains) that make up these kinases. This is supported by: (1) the C-domain of E. coli PGK cannot be expressed or fold independently of the N-domain, while both domains of the yeast protein fold in isolation into stable structures and (2) the energetics and kinetics of the proteolytically sensitive state of E. coli PGK match those for global unfolding. This suggests that proteolysis occurs from the globally unfolded state of E. coli PGK, while the characteristics defining the yeast homolog suggest that proteolysis occurs upon unfolding of only the C-domain, with the N-domain remaining folded and consequently resistant to cleavage.

The folding and evolution of multidomain proteins

Analyses of genomes show that more than 70% of eukaryotic proteins are composed of multiple domains. However, most studies of protein folding focus on individual domains and do not consider how interactions between domains might affect folding. Here, we address this by analysing the three-dimensional structures of multidomain proteins that have been characterized experimentally and observe that where the interface is small and loosely packed, or unstructured, the folding of the domains is independent. Furthermore, recent studies indicate that multidomain proteins have evolved mechanisms to minimize the problems of interdomain misfolding.

The influence of interdomain interactions on the intradomain motions in yeast phosphoglycerate kinase: a molecular dynamics study

A 3-ns molecular dynamics simulation in explicit solvent was performed to examine the inter- and intradomain motions of the two-domain enzyme yeast phosphoglycerate kinase without the presence of substrates. To elucidate contributions from individual domains, simulations were carried out on the complete enzyme as well as on each isolated domain. The enzyme is known to undergo a hinge-bending type of motion as it cycles from an open to a closed conformation to allow the phosphoryl transfer occur. Analysis of the correlation of atomic movements during the simulations confirms hinge bending in the nanosecond timescale: the two domains of the complete enzyme exhibit rigid body motions anticorrelated with respect to each other. The correlation of the intradomain motions of both domains converges, yielding a distinct correlation map in the enzyme. In the isolated domain simulations-in which interdomain interactions cannot occur-the correlation of domain motions no longer converges and shows a very small correlation during the same simulation time. This result points to the importance of interdomain contacts in the overall dynamics of the protein. The secondary structure elements responsible for interdomain contacts are also discussed.

Apparent cooperativity in the folding of multidomain proteins depends on the relative rates of folding of the constituent domains

Approximately 75% of eukaryotic proteins contain more than one so-called independently folding domain. However, there have been relatively few systematic studies to investigate the effect of interdomain interactions on protein stability and fewer still on folding kinetics. We present the folding of pairs of three-helix bundle spectrin domains as a paradigm to indicate how complex such an analysis can be. Equilibrium studies show an increase in denaturant concentration required to unfold the domains with only a single unfolding transition; however, in some cases, this is not accompanied by the increase in m value, which would be expected if the protein is a truly cooperative, all-or-none system. We analyze the complex kinetics of spectrin domain pairs, both wild-type and carefully selected mutants. By comparing these pairs, we are able to demonstrate that equilibrium data alone are insufficient to describe the folding of multidomain proteins and to quantify the effects that one domain can have on its neighbor.

Hierarchic finite level energy landscape model: to describe the refolding kinetics of phosphoglycerate kinase

One of the most intriguing predictions of energy landscape models is the existence of non-exponential protein folding kinetics caused by hierarchical structures in the landscapes. Here we provide the strongest evidence so far of such hierarchy and determine the time constants and weights of the kinetic components of the suggested hierarchic energy landscape. To our knowledge, the idea of hierarchical folding energy barriers has never been tested over such a broad timescale. Refolding of yeast phosphoglycerate kinase was initiated from the guanidine-unfolded state by stopped-flow or manual mixing and monitored by tryptophan fluorescence from 1 ms to 15 min. The strategy to build a model that describes folding of yeast phosphoglycerate kinase was to start from the simplest paradigm and modify it stepwise to the necessary minimal extent after repeated comparisons with the experiments. We made no a priori assumptions about the folding landscape. The result was a hierarchic finite level landscape model that quantitatively describes the refolding of yeast phosphoglycerate kinase from 1 ms to 15 min. The early steps of the folding process happen in the upper region of the landscape, where the surface has a hierarchic structure. This leads to stretched kinetics in the early phase of the folding. The lower region of the energy landscape is dominated by a trap that reflects the accumulation of molten globule intermediate state. From this intermediate, the protein can reach the global energy minimum corresponding to the native state through a cross-barrier folding step.