Tuning the Heterogeneous Early Folding Dynamics of Phosphoglycerate Kinase

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.

Proteins: "Boil 'Em, Mash 'Em, Stick 'Em in a Stew"

Cells of the vast majority of organisms are subject to temperature, pressure, pH, ionic strength, and other stresses. We discuss these effects in the light of protein folding and protein interactions in vitro, in complex environments, in cells, and in vivo. Protein phase diagrams provide a way of organizing different structural ensembles that occur under stress and how one can move among ensembles. Experiments that perturb biomolecules in vitro or in cells by stressing them have revealed much about the underlying forces that are competing to control protein stability, folding, and function. Two phenomena that emerge and serve to broadly classify effects of the cellular environment are crowding (mainly due to repulsive forces) and sticking (mainly due to attractive forces). The interior of cells is closely balanced between these emergent effects, and stress can tip the balance one way or the other. The free energy scale involved is small but significant on the scale of the "on/off switches" that control signaling in cells or of protein-protein association with a favorable function such as increased enzyme processivity. Quantitative tools from biophysical chemistry will play an important role in elucidating the world of crowding and sticking under stress.

Critical phenomena in the temperature-pressure-crowding phase diagram of a protein

In the cell, proteins fold and perform complex functions through global structural rearrangements. Function requires a protein to be at the brink of stability to be susceptible to small environmental fluctuations, yet stable enough to maintain structural integrity. These apparently conflicting behaviors are exhibited by systems near a critical point, where distinct phases merge-a concept beyond previous studies indicating proteins have a well-defined folded/unfolded phase boundary in the pressure-temperature plane. Here, by modeling the protein phosphoglycerate kinase (PGK) on the temperature (T), pressure (P), and crowding volume-fraction (ϕ) phase diagram, we demonstrate a critical transition where phases merge, and PGK exhibits large structural fluctuations. Above the critical point, the difference between the intermediate and unfolded phases disappears. When ϕ increases, the critical point moves to lower T _c. We verify the calculations with experiments mapping the T-P-ϕ space, which likewise reveal a critical point at 305 K and 170 MPa that moves to lower T _c as ϕ increases. Crowding places PGK near a critical line in its natural parameter space, where large conformational changes can occur without costly free energy barriers. Specific structures are proposed for each phase based on simulation.

Heat shock-induced chaperoning by Hsp70 is enabled in-cell

Recent work has shown that weak protein-protein interactions are susceptible to the cellular milieu. One case in point is the binding of heat shock proteins (Hsps) to substrate proteins in cells under stress. Upregulation of the Hsp70 chaperone machinery at elevated temperature was discovered in the 1960s, and more recent studies have shown that ATPase activity in one Hsp70 domain is essential for control of substrate binding by the other Hsp70 domain. Although there are several denaturant-based assays of Hsp70 activity, reports of ATP-dependent binding of Hsp70 to a globular protein substrate under heat shock are scarce. Here we show that binding of heat-inducible Hsp70 to phosphoglycerate kinase (PGK) is remarkably different in vitro compared to in-cell. We use fluorescent-labeled mHsp70 and ePGK, and begin by showing that mHsp70 passes the standard β-galactosidase assay, and that it does not self-aggregate until 50°C in presence of ATP. Yet during denaturant refolding or during in vitro heat shock, mHsp70 shows only ATP-independent non-specific sticking to ePGK, as evidenced by nearly identical results with an ATPase activity-deficient K71M mutant of Hsp70 as a control. Addition of Hsp40 (co-factor) or Ficoll (crowder) does not reduce non-specific sticking, but cell lysate does. Therefore, Hsp70 does not act as an ATP-dependent chaperone on its substrate PGK in vitro. In contrast, we observe only specific ATP-dependent binding of mHsp70 to ePGK in mammalian cells, when compared to the inactive Hsp70 K71M mutant. We hypothesize that enhanced in-cell activity is not due to an unknown co-factor, but simply to a favorable shift in binding equilibrium caused by the combination of crowding and osmolyte/macromolecular interactions present in the cell. One candidate mechanism for such a favorable shift in binding equilibrium is the proven ability of Hsp70 to bind near-native states of substrate proteins in vitro. We show evidence for early onset of binding in-cell. Our results suggest that Hsp70 binds PGK preemptively, prior to its full unfolding transition, thus stabilizing it against further unfolding. We propose a "preemptive holdase" mechanism for Hsp70-substrate binding. Given our result for PGK, more proteins than one might think based on in vitro assays may be chaperoned by Hsp70 in vivo. The cellular environment thus plays an important role in maintaining proper Hsp70 function.

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.

Labeling for Quantitative Comparison of Imaging Measurements in Vitro and in Cells

Qualitative imaging of biomolecular localization and distribution inside cells has revolutionized cell biology. Most of these powerful techniques require modifications to the target biomolecule. Over the past 10 years, these techniques have been extended to quantitative measurements, from in-cell protein folding rates to complex dissociation constants to RNA lifetimes. Such measurements can be affected even when a target molecule is just mildly perturbed by its labels. Here, the impact of labeling on protein (and RNA) structure, stability, and function in cells is discussed via practical examples from the recent literature. General guidelines for selecting and validating modification sites are provided to bring the best from cell biology and imaging to quantitative biophysical experiments inside cells.

Femtosecond to Millisecond Dynamics of Light Induced Allostery in the Avena sativa LOV Domain

The rational engineering of photosensor proteins underpins the field of optogenetics, in which light is used for spatiotemporal control of cell signaling. Optogenetic elements function by converting electronic excitation of an embedded chromophore into structural changes on the microseconds to seconds time scale, which then modulate the activity of output domains responsible for biological signaling. Using time-resolved vibrational spectroscopy coupled with isotope labeling, we have mapped the structural evolution of the LOV2 domain of the flavin binding phototropin Avena sativa (AsLOV2) over 10 decades of time, reporting structural dynamics between 100 fs and 1 ms after optical excitation. The transient vibrational spectra contain contributions from both the flavin chromophore and the surrounding protein matrix. These contributions are resolved and assigned through the study of four different isotopically labeled samples. High signal-to-noise data permit the detailed analysis of kinetics associated with the light activated structural evolution. A pathway for the photocycle consistent with the data is proposed. The earliest events occur in the flavin binding pocket, where a subpicosecond perturbation of the protein matrix occurs. In this perturbed environment, the previously characterized reaction between triplet state isoalloxazine and an adjacent cysteine leads to formation of the adduct state; this step is shown to exhibit dispersive kinetics. This reaction promotes coupling of the optical excitation to successive time-dependent structural changes, initially in the β-sheet and then α-helix regions of the AsLOV2 domain, which ultimately gives rise to Jα-helix unfolding, yielding the signaling state. This model is tested through point mutagenesis, elucidating in particular the key mediating role played by Q513.

ReAsH as a Quantitative Probe of In-Cell Protein Dynamics

The tetracysteine (tc) tag/biarsenical dye system (FlAsH or ReAsH) promises to combine the flexibility of fluorescent protein tags with the small size of dye labels, allowing in-cell study of target proteins that are perturbed by large protein tags. Quantitative thermodynamic and kinetic studies in-cell using FlAsH and ReAsH have been hampered by methodological complexities presented by the fluorescence properties of the tag-dye complex probed by either Förster resonance energy transfer (FRET) or direct excitation. We label the model protein phosphoglycerate kinase (PGK) with AcGFP1 and ReAsH for direct comparison with AcGFP1/mCherry-labeled PGK. We find that fast relaxation imaging (FReI), combining millisecond temperature jump kinetics with fluorescence microscopy detection, circumvents many of the difficulties encountered working with the ReAsH system, allowing us to obtain quantitative FRET measurements of protein stability and kinetics both in vitro and in cells. We also demonstrate the to us surprising result that fluorescence from directly excited, unburied ReAsH at the C-terminus of the model protein also reports on folding in vitro and in cells. Comparing the ReAsH-labeled protein to a construct labeled with two fluorescent protein tags allows us to evaluate how a bulkier protein tag affects protein dynamics in cells and in vitro. We find that the average folding rate in the cell is closer to the in vitro rate with the smaller tag, highlighting the effect of tags on quantitative in-cell measurements.

Dodine as a transparent protein denaturant for circular dichroism and infrared studies

The fungicide dodine combines the cooperative denaturation properties of guanidine with the mM denaturation activity of SDS. It was previously tested only on two small model proteins. Here we show that it can be used as a chemical denaturant for phosphoglycerate kinase (PGK), a much larger two-domain enzyme. In addition to its properties as a chemical denaturant, dodine facilitates thermal denaturation of PGK, and we show for the first time that it also facilitates pressure denaturation of a protein. Much higher quality circular dichroism and amide I' infrared spectra of PGK can be obtained in dodine than in guanidine, opening the possibility for use of dodine as a denaturant when UV or IR detection is desirable. One caution is that dodine denaturation, like other detergent-based denaturants, is less reversible than guanidine denaturation.

The Effect of Fluorescent Protein Tags on Phosphoglycerate Kinase Stability Is Nonadditive

It is frequently assumed that fluorescent protein tags used in biological imaging experiments are minimally perturbing to their host protein. As in-cell experiments become more quantitative and measure rates and equilibrium constants, rather than just "on-off" activity or the presence of a protein, it becomes more important to understand such perturbations. One criterion for a protein modification to be a perturbation is additivity of two perturbations (a linear effect on the protein free energy). Here we show that adding fluorescent protein tags to a host protein in vitro has a large nonadditive effect on its folding free energy. We compare an unlabeled, three singly labeled, and a doubly labeled enzyme (phosphoglycerate kinase). We propose two mechanisms for nonadditivity. In the "quinary interaction" mechanism, two tags interact transiently with one another, relieving the host protein from unfavorable tag-protein interactions. In the "crowding" mechanism, adding two tags provides the minimal crowding necessary to overcome destabilizing interactions of individual tags with the host protein. Both of these mechanisms affect protein stability in cells; we show here that they must also be considered for tagged proteins used for reference in vitro.

Periodic and stochastic thermal modulation of protein folding kinetics

Chemical reactions are usually observed either by relaxation of a bulk sample after applying a sudden external perturbation, or by intrinsic fluctuations of a few molecules. Here we show that the two ideas can be combined to measure protein folding kinetics, either by periodic thermal modulation, or by creating artificial thermal noise that greatly exceeds natural thermal fluctuations. We study the folding reaction of the enzyme phosphoglycerate kinase driven by periodic temperature waveforms. As the temperature waveform unfolds and refolds the protein, its fluorescence color changes due to FRET (Förster resonant Energy Transfer) of two donor/acceptor fluorophores labeling the protein. We adapt a simple model of periodically driven kinetics that nicely fits the data at all temperatures and driving frequencies: The phase shifts of the periodic donor and acceptor fluorescence signals as a function of driving frequency reveal reaction rates. We also drive the reaction with stochastic temperature waveforms that produce thermal fluctuations much greater than natural fluctuations in the bulk. Such artificial thermal noise allows the recovery of weak underlying signals due to protein folding kinetics. This opens up the possibility for future detection of a stochastic resonance for protein folding subject to noise with controllable amplitude.

Protein Denaturation on p-T Axes--Thermodynamics and Analysis

Proteins are essential players in the vast majority of molecular level life processes. Since their structure is in most cases substantial for their correct function, study of their structural changes attracted great interest in the past decades. The three dimensional structure of proteins is influenced by several factors including temperature, pH, presence of chaotropic and cosmotropic agents, or presence of denaturants. Although pressure is an equally important thermodynamic parameter as temperature, pressure studies are considerably less frequent in the literature, probably due to the technical difficulties associated to the pressure studies. Although the first steps in the high-pressure protein study have been done 100 years ago with Bridgman's ground breaking work, the field was silent until the modern spectroscopic techniques allowed the characterization of the protein structural changes, while the protein was under pressure. Recently a number of proteins were studied under pressure, and complete pressure-temperature phase diagrams were determined for several of them. This review summarizes the thermodynamic background of the typical elliptic p-T phase diagram, its limitations and the possible reasons for deviations of the experimental diagrams from the theoretical one. Finally we show some examples of experimentally determined pressure-temperature phase diagrams.

Single molecule conformational memory extraction: p5ab RNA hairpin

Extracting kinetic models from single molecule data is an important route to mechanistic insight in biophysics, chemistry, and biology. Data collected from force spectroscopy can probe discrete hops of a single molecule between different conformational states. Model extraction from such data is a challenging inverse problem because single molecule data are noisy and rich in structure. Standard modeling methods normally assume (i) a prespecified number of discrete states and (ii) that transitions between states are Markovian. The data set is then fit to this predetermined model to find a handful of rates describing the transitions between states. We show that it is unnecessary to assume either (i) or (ii) and focus our analysis on the zipping/unzipping transitions of an RNA hairpin. The key is in starting with a very broad class of non-Markov models in order to let the data guide us toward the best model from this very broad class. Our method suggests that there exists a folding intermediate for the P5ab RNA hairpin whose zipping/unzipping is monitored by force spectroscopy experiments. This intermediate would not have been resolved if a Markov model had been assumed from the onset. We compare the merits of our method with those of others.

Temporal variation of a protein folding energy landscape in the cell

Chemical reaction rate coefficients and free energies are usually time-independent quantities. Protein folding in vitro is one such reaction with a fixed energy landscape. However, in the milieu of the cell, the energy landscape can be modulated in space and time by fluctuations in the intracellular environment such as cytoskeletal rearrangements, changes in biomolecule concentrations, and large scale cellular reorganization. We studied the time dependence of the folding landscape of a FRET-labeled enzyme, yeast phosphoglycerate kinase (PGK-FRET). Living U2OS cells served as our test tube, and the mammalian cell cycle, a process strictly regulated in time, served as our clock. We found that both the rate of folding and the thermodynamic stability of PGK-FRET are cell cycle-dependent. We also assayed folding rates of PGK-FRET in spatial proximity to and far away from mitotic chromosomes. Our results show that expedited folding in DNA-rich regions cannot account for the faster rate of PGK-FRET folding in mitotic cells.

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.

Extracting conformational memory from single-molecule kinetic data

Single-molecule data often come in the form of stochastic time trajectories. A key question is how to extract an underlying kinetic model from the data. A traditional approach is to assume some discrete state model, that is, a model topology, and to assume that transitions between states are Markovian. The transition rates are then selected according to which ones best fit the data. However, in experiments, each apparent state can be a broad ensemble of states or can be hiding multiple interconverting states. Here, we describe a more general approach called the non-Markov memory kernel (NMMK) method. The idea is to begin with a very broad class of non-Markov models and to let the data directly select for the best possible model. To do so, we adapt an image reconstruction approach that is grounded in maximum entropy. The NMMK method is not limited to discrete state models for the data; it yields a unique model given the data, it gives error bars for the model, and it does not assume Markov dynamics. Furthermore, NMMK is less wasteful of data by letting the entire data set determine the model. When the data warrants, the NMMK gives a memory kernel that is Markovian. We highlight, by numerical example, how conformational memory extracted using this method can be translated into useful mechanistic insight.

Temperature dependence of protein folding kinetics in living cells

We measure the stability and folding rate of a mutant of the enzyme phosphoglycerate kinase (PGK) inside bone tissue cells as a function of temperature from 38 to 48 °C. To facilitate measurement in individual living cells, we developed a rapid laser temperature stepping method capable of measuring complete thermal melts and kinetic traces in about two min. We find that this method yields improved thermal melts compared to heating a sample chamber or microscope stage. By comparing results for six cells with in vitro data, we show that the protein is stabilized by about 6 kJ/mole in the cytoplasm, but the temperature dependence of folding kinetics is similar to in vitro. The main difference is a slightly steeper temperature dependence of the folding rate in some cells that can be rationalized in terms of temperature-dependent crowding, local viscosity, or hydrophobicity. The observed rate coefficients can be fitted within measurement uncertainty by an effective two-state model, even though PGK folds by a multistate mechanism. We validate the effective two-state model with a three-state free energy landscape of PGK to illustrate that the effective fitting parameters can represent a more complex underlying free energy landscape.

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.

Kinetic analysis of protein aggregation monitored by real-time 2D solid-state NMR spectroscopy

It is shown that real-time 2D solid-state NMR can be used to obtain kinetic and structural information about the process of protein aggregation. In addition to the incorporation of kinetic information involving intermediate states, this approach can offer atom-specific resolution for all detectable species. The analysis was carried out using experimental data obtained during aggregation of the 10.4 kDa Crh protein, which has been shown to involve a partially unfolded intermediate state prior to aggregation. Based on a single real-time 2D (13)C-(13)C transition spectrum, kinetic information about the refolding and aggregation step could be extracted. In addition, structural rearrangements associated with refolding are estimated and several different aggregation scenarios were compared to the experimental data.

Structure, function, and folding of phosphoglycerate kinase are strongly perturbed by macromolecular crowding

We combine experiment and computer simulation to show how macromolecular crowding dramatically affects the structure, function, and folding landscape of phosphoglycerate kinase (PGK). Fluorescence labeling shows that compact states of yeast PGK are populated as the amount of crowding agents (Ficoll 70) increases. Coarse-grained molecular simulations reveal three compact ensembles: C (crystal structure), CC (collapsed crystal), and Sph (spherical compact). With an adjustment for viscosity, crowded wild-type PGK and fluorescent PGK are about 15 times or more active in 200 mg/ml Ficoll than in aqueous solution. Our results suggest a previously undescribed solution to the classic problem of how the ADP and diphosphoglycerate binding sites of PGK come together to make ATP: Rather than undergoing a hinge motion, the ADP and substrate sites are already located in proximity under crowded conditions that mimic the in vivo conditions under which the enzyme actually operates. We also examine T-jump unfolding of PGK as a function of crowding experimentally. We uncover a nonmonotonic folding relaxation time vs. Ficoll concentration. Theory and modeling explain why an optimum concentration exists for fastest folding. Below the optimum, folding slows down because the unfolded state is stabilized relative to the transition state. Above the optimum, folding slows down because of increased viscosity.

Thermodynamics and kinetics of the pressure unfolding of phosphoglycerate kinase

Due to the relationship between compressibility and volume fluctuations, high-pressure studies provide vital insight into protein dynamics and function. Most high-pressure experiments were performed on small and fast folding proteins or model peptides. Here we show that a detailed kinetic study is necessary to extract reliable information from the high-pressure-induced structural conversion of large, slowly folding proteins. The pressure-jump unfolding kinetics of yeast phosphoglycerate kinase was recorded at pressures between 50 and 150 MPa. The time dependence of the conformational state of the protein was followed by tryptophan fluorescence measurements from 30 s to 2 h. The observed changes were described by a three-state model, and the volume change and the activation volume as well as the midpoint pressure of the transitions between the folded, intermediate, and unfolded states were determined. An interesting feature of the pressure unfolding of phosphoglycerate kinase was that the unfolding process speeds up with increasing pressure, which is the consequence of negative activation volumes for the folded --> intermediate, intermediate --> unfolded, and unfolded --> intermediate transitions.

Protein dynamics: from molecules, to interactions, to biology

Proteins have a remarkably rich diversity of dynamical behaviors, and the articles in this issue of the International Journal of Molecular Sciences are a testament to that fact. From the picosecond motions of single sidechains probed by NMR or fluorescence spectroscopy, to aggregation processes at interfaces that take months, all time scales play a role. Proteins are functional molecules, so by their nature they always interact with their environment. This environment includes water, other biomolecules, or larger cellular structures. In a sense, it also includes the protein molecule itself: proteins are large enough to fold and interact with themselves. These interactions have been honed by evolution to produce behaviors completely different from those of random polymers.

Folding time distributions as an approach to protein folding kinetics

A 27-residue lattice heteropolymer subject to Monte Carlo dynamics on a simple cubic lattice is studied over a range of temperatures. Folding time distributions are used to obtain information concerning the details of folding kinetics. The results are compared with those from methods based on mean force surfaces expressed in terms of a reduced set of variables and on a disconnectivity graph for the same system. A detailed analysis of the folding trajectories is given, and the importance of dead-end traps in determining the folding time is demonstrated. We show that the calculated folding kinetics can be modeled by a system of kinetic equations, with the essential rate constants determined from the Monte Carlo simulations and the resulting folding time distributions. The kinetic equations make possible an analysis of the variation of the importance of different channels with temperature. In particular, we show that the presence of intermediates may be masked in the folding time distributions, with the mean folding time being independent of the height of the barrier between the intermediates and collapsed globule state of the system. This and other results demonstrate that care has to be used in interpreting experimental folding data in terms of the underlying kinetics. Correspondingly, simulations are shown to have to satisfy certain requirements to obtain proper sampling of the dead-end traps.

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.

Mechanistic Basis for the Action of New Cephalosporin Antibiotics Effective against Methicillin- and Vancomycin-resistant Staphylococcus aureus

Emergence of methicillin-resistant Staphylococcus aureus (MRSA) has created challenges in treatment of nosocomial infections. The recent clinical emergence of vancomycin-resistant MRSA is a new disconcerting chapter in the evolution of these strains. S. aureus normally produces four PBPs, which are susceptible to modification by beta-lactam antibiotics, an event that leads to bacterial death. The gene product of mecA from MRSA is a penicillin-binding protein (PBP) designated PBP 2a. PBP 2a is refractory to the action of all commercially available beta-lactam antibiotics. Furthermore, PBP 2a is capable of taking over the functions of the other PBPs of S. aureus in the face of the challenge by beta-lactam antibiotics. Three cephalosporins (compounds 1-3) have been studied herein, which show antibacterial activities against MRSA, including the clinically important vancomycin-resistant strains. These cephalosporins exhibit substantially smaller dissociation constants for the preacylation complex compared with the case of typical cephalosporins, but their pseudo-second-order rate constants for encounter with PBP 2a (k(2)/K(s)) are not very large (< or =200 m(-1) s(-1)). It is documented herein that these cephalosporins facilitate a conformational change in PBP 2a, a process that is enhanced in the presence of a synthetic surrogate for cell wall, resulting in increases in the k(2)/K(s) parameter and in more facile enzyme inhibition. These findings argue that the novel cephalosporins are able to co-opt interactions between PBP 2a and the cell wall in gaining access to the active site in the inhibition process, a set of events that leads to effective inhibition of PBP 2a and the attendant killing of the MRSA strains.

Domain interactions direct misfolding and amyloid formation of yeast phosphoglycerate kinase

There are proteins that are built of two structural domains and are deposited full-length in amyloid plaques formed in various diseases. In spite of the known differences in the mechanisms of folding of single- and multidomain proteins, no published studies can be found that address the role of the domain-domain interactions during misfolding and amyloid formation. By the discovery of the role of domain-domain interactions, here we provide important insight in the submolecular mechanism of amyloid formation. A model system based on yeast phosphoglycerate kinase was designed. This system includes the wild-type yeast phosphoglycerate kinase and single-tryptophan mutants of the individual N and C terminal domains and the complete protein. Electron microscopic measurements proved that amyloid fibrils grow from all mutants under identical conditions as for the wild-type protein. Misfolding and amyloid formation was followed in stopped-flow and manual mixing experiments on the 1 ms to 4 days timescale. Tryptophan fluorescence was used for selective detection of conformational changes accompanying the formation of the amyloidogenic intermediates and the growth of amyloid fibrils. The interactions between the polypeptide chains of the two domains direct the misfolding process from the early steps to the amyloid formation, and influence the final structure. The kinetics of misfolding is different for the individual domains, pointing to the significance of the amino acid sequence. Misfolding of the domains within the complete protein is synchronized indicating that domain-domain interactions direct the misfolding and amyloid formation mechanism.

Correlation between conformational stability of the ternary enzyme-substrate complex and domain closure of 3-phosphoglycerate kinase

3-phosphoglycerate kinase (PGK) is a typical two-domain hinge-bending enzyme with a well-structured interdomain region. The mechanism of domain-domain interaction and its regulation by substrate binding is not yet fully understood. Here the existence of strong cooperativity between the two domains was demonstrated by following heat transitions of pig muscle and yeast PGKs using differential scanning microcalorimetry and fluorimetry. Two mutants of yeast PGK containing a single tryptophan fluorophore either in the N- or in the C-terminal domain were also studied. The coincidence of the calorimetric and fluorimetric heat transitions in all cases indicated simultaneous, highly cooperative unfolding of the two domains. This cooperativity is preserved in the presence of substrates: 3-phosphoglycerate bound to the N domain or the nucleotide (MgADP, MgATP) bound to the C domain increased the structural stability of the whole molecule. A structural explanation of domain-domain interaction is suggested by analysis of the atomic contacts in 12 different PGK crystal structures. Well-defined backbone and side-chain H bonds, and hydrophobic and electrostatic interactions between side chains of conserved residues are proposed to be responsible for domain-domain communication. Upon binding of each substrate newly formed molecular contacts are identified that firstly explain the order of the increased heat stability in the various binary complexes, and secondly describe the possible route of transmission of the substrate-induced conformational effects from one domain to the other. The largest stability is characteristic of the native ternary complex and is abolished in the case of a chemically modified inactive form of PGK, the domain closure of which was previously shown to be prevented [Sinev MA, Razgulyaev OI, Vas M, Timchenko AA & Ptitsyn OB (1989) Eur J Biochem180, 61-66]. Thus, conformational stability correlates with domain closure that requires simultaneous binding of both substrates.

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

The aim of this work is to shed more light on the effect of domain-domain interactions on the kinetics and the pathway of protein folding. A model protein system consisting of several single-tryptophan variants of the two-domain yeast phosphoglycerate kinase (PGK) and its individual domains was studied. Refolding was initiated from the guanidine-unfolded state by stopped-flow or manual mixing and monitored by tryptophan fluorescence from 1 msec to 1000 sec. Denaturant titrations of both individual domains showed apparent two-state unfolding transitions. Refolding kinetics of the individual domains from different denaturant concentrations, however, revealed the presence of intermediate structures during titration for both domains. Refolding of the same domains within the complete protein showed that domain-domain interactions direct the folding of both domains, but in an asymmetric way. Folding of the N domain was already altered within 1 msec, while detectable changes in the folding of the C domain occurred only 60-100 msec after initiating refolding. All mutants showed a hyperfluorescent kinetic intermediate. Both the disappearance of this intermediate and the completion of the folding were significantly faster in the individual N domain than in the complete protein. On the contrary, folding of the individual C domain was slower than in the complete protein. The presence of the C domain directs the refolding of the N domain along a completely different pathway than that of the individual N domain, while folding of the individual C domain follows the same path as within the complete protein.

Downhill protein folding: evolution meets physics

Proteins can be redesigned to fold downhill on a free energy surface characterized by only a few coordinates, confirming a principal prediction of the 'energy-landscape' model. Nonetheless, natural proteins have small but significant barriers. Spectroscopy and kinetics reveal potential biological causes for activation barriers during protein folding: evolution against protein aggregation and for protein function.

Kinetics are probe-dependent during downhill folding of an engineered lambda6-85 protein

The Y22W/Q33Y/G46,48A mutant of the protein lambda6-85 folds in a few microseconds at room temperature. We find that its folding kinetics are probe-dependent under a strong bias toward the native state, a new signature for downhill folding. The IR- and fluorescence-detected relaxation time scales converge when the native bias is removed by raising the temperature, recovering activated two-state folding. Langevin dynamics simulations on one- and 2D free energy surfaces tunable from two-state to downhill folding reproduce the difference between the IR and fluorescence experiments, as well as the temperature and viscosity trends. In addition, the 2D surface reproduces the stretched exponential dynamics that we fit to the glucose solution experimental data at short times. Nonexponential dynamics at <10 micros is a signature either for local free energy minima along the reaction coordinate ("longitudinal roughness"), or for folding on a higher-dimensional free energy surface ("transverse roughness").

Rate-temperature relationships in lambda-repressor fragment lambda 6-85 folding

Two classes of lambda(6-85) mutants (those richer in alanine, and those richer in glycine) have very similar slopes in an Arrhenius plot of the unfolding rates but very different temperature dependencies of the folding rates. Temperature-dependent interactions (e.g., hydrophobicity) play a large role in the initial stages of folding but not in the initial stages of unfolding of lambda(6-85). Placement of the transition state in terms of its surface exposure and entropy shows that at least two reaction coordinates are required to describe folding of all mutants over the full temperature range. The unusual Arrhenius plots of the very fastest mutant provide an additional kinetic signature for downhill folding.

Multiple probes reveal a native-like intermediate during low-temperature refolding of ubiquitin

We investigate the refolding of ubiquitin Phe45Trp/Ile61Ala (Ub(*)I61A) in a low-temperature, high-viscosity buffer, where folding is slowed so that apparent two-state and three-state mechanisms are readily distinguishable. Ub(*)I61A forms a compact ensemble rapidly (as judged from stopped-flow, small-angle X-ray scattering) with a secondary structure signature similar to that of the native state (as judged from stopped-flow circular dichroism from 215 nm to 250 nm), but the fluorescence signature still resembles the guanidinium-denatured state. The compact ensemble forms over a range of solvent and temperature conditions. The native fluorescence signature, which requires the tryptophan residue to be packed tightly, is acquired at least 500 times more slowly. Molecular dynamics simulations at 495 K show no contraction of the backbone in ethylene glycol buffer compared to pure aqueous buffer, and no significant effect on the local backbone structure of the unfolded protein. Only at higher simulation temperature does a backbone contraction appear. Thus, it appears unlikely that the aqueous ethylene glycol buffer fundamentally changes the folding mechanism of ubiquitin. We suggest that ubiquitin forms a compact ensemble with native-like secondary structure, but without tight packing, long before the native state.

Folding lambda-repressor at its speed limit

We show that the five-helix bundle lambda(6-85) can be engineered and solvent-tuned to make the transition from activated two-state folding to downhill folding. The transition manifests itself as the appearance of additional dynamics faster than the activated kinetics, followed by the disappearance of the activated kinetics when the bias toward the native state is increased. Our fastest value of 1 micros for the "speed" limit of lambda(6-85) is measured at low concentrations of a denaturant that smooths the free-energy surface. Complete disappearance of the activated phase is obtained in stabilizing glucose buffer. Langevin dynamics on a rough free-energy surface with variable bias toward the native state provides a robust and quantitative description of the transition from activated to downhill folding. Based on our simulation, we estimate the residual energetic frustration of lambda(6-85) to be delta(2) G approximately 0.64 k(2)T(2). We show that lambda(6-86), as well as very fast folding proteins or folding intermediates estimated to lie near the speed limit, provide a better rate-topology correlation than proteins with larger energetic frustration. A limit of beta > or = 0.7 on any stretching of lambda(6-85) barrier-free dynamics suggests that a low-dimensional free-energy surface is sufficient to describe folding.

Heterogeneous Folding of the trpzip Hairpin: Full Atom Simulation and Experiment

The beta-hairpin trpzip2 can be tuned continuously from a two-state folder to folding on a rough energy landscape without a dominant refolding barrier. At high denaturant concentration, this extremely stable peptide exhibits a single apparent "two-state" transition temperature when monitored by different spectroscopic techniques. However, under optimal folding conditions the hairpin undergoes an unusual folding process with three clusters of melting transitions ranging from 15 degrees C to 160 degrees C, as monitored by 12 different experimental and computational observables. We explain this behavior in terms of a rough free energy landscape of the unfolded peptide caused by multiple tryptophan interactions and alternative backbone conformations. The landscape is mapped out by potentials of mean force derived from replica-exchange molecular dynamics simulations. Implications for deducing cooperativity from denaturant titrations, for the origin of folding cooperativity, and for the folding of thermophilic proteins are pointed out. trpzip is an excellent small tunable model system for the glass-like folding transitions predicted by landscape theory.