Direct observation of an ensemble of stable collapsed states in the mechanical folding of ubiquitin

Statistical theories of protein folding have long predicted plausible mechanisms for reducing the vast conformational space through distinct ensembles of structures. However, these predictions have remained untested by bulk techniques, because the conformational diversity of folding molecules has been experimentally unapproachable. Owing to recent advances in single molecule force-clamp spectroscopy, we are now able to probe the structure and dynamics of the small protein ubiquitin by measuring its length and mechanical stability during each stage of folding. Here, we discover that upon hydrophobic collapse, the protein rapidly selects a subset of minimum energy structures that are mechanically weak and essential precursors of the native fold. From this much reduced ensemble, the native state is acquired through a barrier-limited transition. Our results support the validity of statistical mechanics models in describing the folding of a small protein on biological timescales.     

Making connections between ultrafast protein folding kinetics and molecular dynamics simulations

Determining the rate of forming the truly folded conformation of ultrafast folding proteins is an important issue for both experiments and simulations. The double-norleucine mutant of the 35-residue villin subdomain is the focus of recent computer simulations with atomistic molecular dynamics because it is currently the fastest folding protein. The folding kinetics of this protein have been measured in laser temperature-jump experiments using tryptophan fluorescence as a probe of overall folding. The conclusion from the simulations, however, is that the rate determined by fluorescence is significantly larger than the rate of overall folding. We have therefore employed an independent experimental method to determine the folding rate. The decay of the tryptophan triplet-state in photoselection experiments was used to monitor the change in the unfolded population for a sequence of the villin subdomain with one amino acid difference from that of the laser temperature-jump experiments, but with almost identical equilibrium properties. Folding times obtained in a two-state analysis of the results from the two methods at denaturant concentrations varying from 1.5-6.0 M guanidinium chloride are in excellent agreement, with an average difference of only 20%. Polynomial extrapolation of all the data to zero denaturant yields a folding time of 220 (+100,-70) ns at 283 K, suggesting that under these conditions the barrier between folded and unfolded states has effectively disappeared--the so-called "downhill scenario."     

Calcium-dependent folding of single calmodulin molecules

Calmodulin is the primary calcium binding protein in living cells. Its function and structure depend strongly on calcium concentration. We used single molecule force spectroscopy by optical tweezers to study the folding of calmodulin in the physiologically relevant range. We find that full-length calmodulin switches from a rich and complex folding behavior at high calcium to a simple folding pathway at apo conditions. Using truncation mutants, we studied the individual domains separately. Folding and stability of the individual domains differ significantly at low calcium concentrations. With increasing calcium, the folding rate constants increase while unfolding rate constants decrease. The complete kinetic as well as energetic behavior of both domains could be modeled using a calcium-dependent three-pathway model. We find that the dominant folding pathway at high calcium concentrations proceeds via a transition state capable of binding one calcium ion. The folding of calmodulin seems to be designed to occur fast robustly over a large range of calcium concentrations and hence energetic stabilities.     

Direct observation of multiple misfolding pathways in a single prion protein molecule

Protein misfolding is a ubiquitous phenomenon associated with a wide range of diseases. Single-molecule approaches offer a powerful tool for deciphering the mechanisms of misfolding by measuring the conformational fluctuations of a protein with high sensitivity. We applied single-molecule force spectroscopy to observe directly the misfolding of the prion protein PrP, a protein notable for having an infectious misfolded state that is able to propagate by recruiting natively folded PrP. By measuring folding trajectories of single PrP molecules held under tension in a high-resolution optical trap, we found that the native folding pathway involves only two states, without evidence for partially folded intermediates that have been proposed to mediate misfolding. Instead, frequent but fleeting transitions were observed into off-pathway intermediates. Three different misfolding pathways were detected, all starting from the unfolded state. Remarkably, the misfolding rate was even higher than the rate for native folding. A mutant PrP with higher aggregation propensity showed increased occupancy of some of the misfolded states, suggesting these states may act as intermediates during aggregation. These measurements of individual misfolding trajectories demonstrate the power of single-molecule approaches for characterizing misfolding directly by mapping out nonnative folding pathways.     

Early events in the folding of four-helix-bundle heme proteins

Topologically homologous four-helix-bundle heme proteins exhibit striking diversity in their refolding kinetics. Cytochrome b562 has been reported to fold on a sub-millisecond time scale, whereas cytochrome c' refolding requires 10 s or more to complete. Heme dissociation in cytochrome b562 interferes with studies of folding kinetics, so a variant of cytochrome b562 (cytochrome c-b562) with a covalent c-type linkage to the heme has been expressed in Escherichia coli. Early events in the electron transfer-triggered folding of Fe(II)-cytochrome c-b562, along with those of Fe(II)-cytochrome c556, have been examined by using time-resolved absorption spectroscopy. Coordination of S(Met) to Fe(II) occurs within 10 mus after reduction of the denatured Fe(III)-cytochromes, and shortly thereafter (100 micros) the heme spectra are indistinguishable from those of the folded proteins. Under denaturing conditions, carbon monoxide binds to the Fe(II)-hemes in approximately 15 ms. By contrast, CO binding cannot compete with refolding in the Fe(II)-cytochromes, thereby confirming that the polypeptide encapsulates the heme in <10 ms. We suggest that Fe-S(Met) ligation facilitates refolding in these four-helix-bundle heme proteins by reducing the conformational freedom of the polypeptide chain.     

Full distance-resolved folding energy landscape of one single protein molecule

Kinetic bulk and single molecule folding experiments characterize barrier properties but the shape of folding landscapes between barrier top and native state is difficult to access. Here, we directly extract the full free energy landscape of a single molecule of the GCN4 leucine zipper using dual beam optical tweezers. To this end, we use deconvolution force spectroscopy to follow an individual molecule’s trajectory with high temporal and spatial resolution. We find a heterogeneous energy landscape of the GCN4 leucine zipper domain. The energy profile is divided into two stable C-terminal heptad repeats and two less stable repeats at the N-terminus. Energies and transition barrier positions were confirmed by single molecule kinetic analysis. We anticipate that deconvolution sampling is a powerful tool for the model-free investigation of protein energy landscapes.     

Direct visualization reveals dynamics of a transient intermediate during protein assembly

Interactions between proteins underlie numerous biological functions. Theoretical work suggests that protein interactions initiate with formation of transient intermediates that subsequently relax to specific, stable complexes. However, the nature and roles of these transient intermediates have remained elusive. Here, we characterized the global structure, dynamics, and stability of a transient, on-pathway intermediate during complex assembly between the Signal Recognition Particle (SRP) and its receptor. We show that this intermediate has overlapping but distinct interaction interfaces from that of the final complex, and it is stabilized by long-range electrostatic interactions. A wide distribution of conformations is explored by the intermediate; this distribution becomes more restricted in the final complex and is further regulated by the cargo of SRP. These results suggest a funnel-shaped energy landscape for protein interactions, and they provide a framework for understanding the role of transient intermediates in protein assembly and biological regulation.     

Millisecond spatiotemporal dynamics of FRET biosensors by the pair correlation function and the phasor approach to FLIM

Here we present a fluctuation-based approach to biosensor Förster resonance energy transfer (FRET) detection that can measure the molecular flow and signaling activity of proteins in live cells. By simultaneous use of the phasor approach to fluorescence lifetime imaging microscopy (FLIM) and cross–pair correlation function (pCF) analysis along a line scanned in milliseconds, we detect the spatial localization of Rho GTPase activity (biosensor FRET signal) as well as the diffusive route adopted by this active population. In particular we find, for Rac1 and RhoA, distinct gradients of activation (FLIM-FRET) and a molecular flow pattern (pCF analysis) that explains the observed polarized GTPase activity. This multiplexed approach to biosensor FRET detection serves as a unique tool for dissection of the mechanism(s) by which key signaling proteins are spatially and temporally coordinated.     

Stable folding core in the folding transition state of an alpha-helical integral membrane protein

Defining the structural features of a transition state is important in understanding a folding reaction. Here, we use Φ-value and double mutant analyses to probe the folding transition state of the membrane protein bacteriorhodopsin. We focus on the final C-terminal helix, helix G, of this seven transmembrane helical protein. Φ-values could be derived for 12 amino acid residues in helix G, most of which have low or intermediate values, suggesting that native structure is disrupted at these amino acid positions in the transition state. Notably, a cluster of residues between E204 and M209 all have Φ-values close to zero. Disruption of helix G is further confirmed by a low Φ-value of 0.2 between residues T170 on helix F and S226 on helix G, suggesting the absence of a native hydrogen bond between helices F and G. Φ-values for paired mutations involved in four interhelical hydrogen bonds revealed that all but one of these bonds is absent in the transition state. The unstructured helix G contrasts with Φ-values along helix B that are generally high, implying native structure in helix B in the transition state. Thus helix B seems to constitute part of a stable folding nucleus while the consolidation of helix G is a relatively late folding event. Polarization of secondary structure correlates with sequence position, with a structured helix B near the N terminus contrasting with an unstructured C-terminal helix G.     

The transition state for folding of an outer membrane protein

Inspired by the seminal work of Anfinsen, investigations of the folding of small water-soluble proteins have culminated in detailed insights into how these molecules attain and stabilize their native folds. In contrast, despite their overwhelming importance in biology, progress in understanding the folding and stability of membrane proteins remains relatively limited. Here we use mutational analysis to describe the transition state involved in the reversible folding of the β-barrel membrane protein PhoPQ-activated gene P (PagP) from a highly disordered state in 10 M urea to a native protein embedded in a lipid bilayer. Analysis of the equilibrium stability and unfolding kinetics of 19 variants that span all eight β-strands of this 163-residue protein revealed that the transition-state structure is a highly polarized, partly formed β-barrel. The results provide unique and detailed insights into the transition-state structure for β-barrel membrane protein folding into a lipid bilayer and are consistent with a model for outer membrane protein folding via a tilted insertion mechanism.     

Solvation in protein (un)folding of melittin tetramer–monomer transition

Protein structural integrity and flexibility are intimately tied to solvation. Here, we examine the effect that changes in bulk and local solvent properties have on protein structure and stability. We observe the change in solvation of an unfolding of the protein model, melittin, in the presence of a denaturant, trifluoroethanol. The peptide system displays a well defined transition in that the tetramer unfolds without disrupting the secondary or tertiary structure. In the absence of local structural perturbation, we are able to reveal exclusively the role of solvation dynamics in protein structure stabilization and the (un)folding pathway. A sudden retardation in solvent dynamics, which is coupled to the change in protein structure, is observed at a critical trifluoroethanol concentration. The large amplitude conformational changes are regulated by the local solvent hydrophobicity and bulk solvent viscosity.     

Role of structural determinants in folding of the sandwich-like protein Pseudomonas aeruginosa azurin

An invariant substructure that forms two interlocked pairs of neighboring beta-strands occurs in essentially all known sandwich-like proteins. Eight conserved positions in these strands were recently shown to act as structural determinants. To test whether the residues at these invariant positions are conserved for mechanistic (i.e., part of folding nucleus) or energetic (i.e., governing native-state stability) reasons, we characterized the folding behavior of eight point-mutated variants of the sandwich-like protein Pseudomonas aeruginosa apo-azurin. We find a simple relationship among the conserved positions: half of the residues form native-like interactions in the folding transition state, whereas the others do not participate in the folding nucleus but govern high native-state stability. Thus, evolutionary preservation of these specific positions gives both mechanistic and energetic advantages to members of the sandwich-like protein family.     

Structure and folding of a designed knotted protein

A very small number of natural proteins have folded configurations in which the polypeptide backbone is knotted. Relatively little is known about the folding energy landscapes of such proteins, or how they have evolved. We explore those questions here by designing a unique knotted protein structure. Biophysical characterization and X-ray crystal structure determination show that the designed protein folds to the intended configuration, tying itself in a knot in the process, and that it folds reversibly. The protein folds to its native, knotted configuration approximately 20 times more slowly than a control protein, which was designed to have a similar tertiary structure but to be unknotted. Preliminary kinetic experiments suggest a complicated folding mechanism, providing opportunities for further characterization. The findings illustrate a situation where a protein is able to successfully traverse a complex folding energy landscape, though the amino acid sequence of the protein has not been subjected to evolutionary pressure for that ability. The success of the design strategy--connecting two monomers of an intertwined homodimer into a single protein chain--supports a model for evolution of knotted structures via gene duplication.     

Direct observation of fast protein conformational switching

Folded proteins can exist in multiple conformational substates. Each substate reflects a local minimum on the free-energy landscape with a distinct structure. By using ultrafast 2D-IR vibrational echo chemical-exchange spectroscopy, conformational switching between two well defined substates of a myoglobin mutant is observed on the approximately 50-ps time scale. The conformational dynamics are directly measured through the growth of cross peaks in the 2D-IR spectra of CO bound to the heme active site. The conformational switching involves motion of the distal histidine/E helix that changes the location of the imidazole side group of the histidine. The exchange between substates changes the frequency of the CO, which is detected by the time dependence of the 2D-IR vibrational echo spectrum. These results demonstrate that interconversion between protein conformational substates can occur on very fast time scales. The implications for larger structural changes that occur on much longer time scales are discussed.     

Redox cycling and kinetic analysis of single molecules of solution-phase nitrite reductase

Single-molecule measurements are a valuable tool for revealing details of enzyme mechanisms by enabling observation of unsynchronized behavior. However, this approach often requires immobilizing the enzyme on a substrate, a process which may alter enzyme behavior. We apply a microfluidic trapping device to allow, for the first time, prolonged solution-phase measurement of single enzymes in solution. Individual redox events are observed for single molecules of a blue nitrite reductase and are used to extract the microscopic kinetic parameters of the proposed catalytic cycle. Changes in parameters as a function of substrate concentration are consistent with a random sequential substrate binding mechanism.     

Cooperative folding kinetics of BBL protein and peripheral subunit-binding domain homologues

Recent experiments claiming that Naf-BBL protein follows a global downhill folding raised an important controversy as to the folding mechanism of fast-folding proteins. Under the global downhill folding scenario, not only do proteins undergo a gradual folding, but folding events along the continuous folding pathway also could be mapped out from the equilibrium denaturation experiment. Based on the exact calculation using a free energy landscape, relaxation eigenmodes from a master equation, and Monte Carlo simulation of an extended Muñoz–Eaton model that incorporates multiscale-heterogeneous pairwise interactions between amino acids, here we show that the very nature of a two-state cooperative transition such as a bimodal distribution from an exact free energy landscape and biphasic relaxation kinetics manifest in the thermodynamics and folding–unfolding kinetics of BBL and peripheral subunit-binding domain homologues. Our results provide an unequivocal resolution to the fundamental controversy related to the global downhill folding scheme, whose applicability to other proteins should be critically reexamined.     

Mapping dynamic protein interactions in MAP kinase signaling using live-cell fluorescence fluctuation spectroscopy and imaging

Fluorescence correlation spectroscopy (FCS), fluorescence cross-correlation spectroscopy (FCCS), and photon counting histograms (PCH) are fluctuation methods that emerged recently as potentially useful tools for obtaining parameters of molecular dynamics, interactions, and oligomerization in vivo. Here, we report the successful implementation of FCS, FCCS, and PCH in live yeast cells using fluorescent protein-tagged proteins expressed from their native chromosomal loci, examining cytosolic dynamics and interactions among components of the mitogen activated protein kinase (MAPK) cascade, a widely occurring signaling motif, in response to mating pheromone. FCS analysis detailed the diffusion characteristics and mobile concentrations of MAPK proteins. FCCS analysis using EGFP and mCherry-tagged protein pairs observed the interactions of Ste7 (MAPK kinase) with the MAPKs, Fus3 or Kss1, and of the scaffold protein, Ste5, with Ste7 and Ste11 (MAPK kinase kinase) in the cytosol, providing in vivo constants of their binding equilibrium. The interaction of Ste5 with Fus3 in the cytosol was below the limit of detection, suggesting a weak interaction, if it exists, with K(d) >400-500 nM. Using PCH, we show that cytosolic Ste5 were mostly monomers. Artificial dimerization of Ste5, as confirmed by PCH, using a dimerizing tag, stimulated the interaction between Ste5 and Fus3. Native Ste5 was found to bind Fus3 preferentially at the cortex in pheromone-treated cells, as detected by fluorescence resonance energy transfer (FRET). These results provide a quantitative spatial map of MAPK complexes in vivo and directly support the model that membrane association and regulation of the Ste5 scaffold are critical steps in MAPK activation.     

Direct single-molecule observation of a protein living in two opposed native structures

Biological activity in proteins requires them to share the energy landscape for folding and global conformational motions, 2 key determinants of function. Although most structural studies to date have focused on fluctuations around a single structural basin, we directly observe the coexistence of 2 symmetrically opposed conformations for a mutant of the Rop-homodimer (Repressor of Primer) in single-molecule fluorescence resonance energy transfer (smFRET) measurements. We find that mild denaturing conditions can affect the sensitive balance between the conformations, generating an equilibrium ensemble consisting of 2 equally occupied structural basins. Despite the need for large-scale conformational rearrangement, both native structures are dynamically and reversibly adopted for the same paired molecules without separation of the constituent monomers. Such an ability of some proteins or protein complexes to switch between conformations by thermal fluctuations and/or minor environmental changes could be central to their ability to control biological function.     

Molecular chaperone function of Mia40 triggers consecutive induced folding steps of the substrate in mitochondrial protein import

Several proteins of the mitochondrial intermembrane space are targeted by internal targeting signals. A class of such proteins with α-helical hairpin structure bridged by two intramolecular disulfides is trapped by a Mia40-dependent oxidative process. Here, we describe the oxidative folding mechanism underpinning this process by an exhaustive structural characterization of the protein in all stages and as a complex with Mia40. Two consecutive induced folding steps are at the basis of the protein-trapping process. In the first one, Mia40 functions as a molecular chaperone assisting α-helical folding of the internal targeting signal of the substrate. Subsequently, in a Mia40-independent manner, folding of the second substrate helix is induced by the folded targeting signal functioning as a folding scaffold. The Mia40-induced folding pathway provides a proof of principle for the general concept that internal targeting signals may operate as a folding nucleus upon compartment-specific activation.