Cytochrome b562 folding triggered by electron transfer: approaching the speed limit for formation of a four-helix-bundle protein

Construction of An Artificial Photosynthesis System with A Single CdS QDs-Ferritin Hybrid Molecule

Establishing artificial photosynthesis systems in a simple but effective manner to mitigate the greenhouse effect and address the energy crisis remains challenging. The combination of photocatalysis technology with bioengineering is an emerging field with great potential to construct such artificial photosynthesis systems, but so far, it has barely been explored in this area. Herein, an artificial photocatalysis platform is constructed with high CO2 conversion and H2O splitting capability by integration of CdS QDs into the intra-subunit interface of H-type ferritin (Marsupenaeus japonicus), a natural ferroxidase through protein interface redesign. The crystal structure of the synthesized CdS QDs with engineered ferritin molecules as bio-templates confirmed the design at an atomic level. Notably, upon absorbing desirable visible light (≈420 nm), such a single CdS-ferritin hybrid molecule is able to selectively catalyze the reduction of CO2 into HCOOH (≈90%), meanwhile catalyzing the oxidation of H2O into O2 in an aqueous environment. The O2 production rate reached to 180 µmol g-1 h-1, and the HCOOH output hit almost 800 µmol g-1 h-1. This work advances the utilization of the four-helix bundle structure for crafting artificial photosynthesis systems, facilitating the seamless integration of bioengineering and photocatalysis technology.

Profiling the Heme-Binding Proteomes of Bacteria Using Chemical Proteomics

Heme is a cofactor with myriad roles and essential to almost all living organisms. Beyond classical gas transport and catalytic functions, heme is increasingly appreciated as a tightly controlled signalling molecule regulating protein expression. However, heme acquisition, biosynthesis and regulation is poorly understood beyond a few model organisms, and the heme-binding proteome has not been fully characterised in bacteria. Yet as heme homeostasis is critical for bacterial survival, heme-binding proteins are promising drug targets. Herein we report a chemical proteomics method for global profiling of heme-binding proteins in live cells for the first time. Employing a panel of heme-based clickable and photoaffinity probes enabled the profiling of 32-54 % of the known heme-binding proteomes in Gram-positive and Gram-negative bacteria. This simple-to-implement profiling strategy could be interchangeably applied to different cell types and systems and fuel future research into heme biology.

Directed Self-Assembly of Dimeric Building Blocks into Networklike Protein Origami to Construct Hydrogels

Engineering proteins to construct self-assemblies is of crucial significance not only for understanding the sophisticated living systems but also for fabricating advanced materials with unexplored functions. However, due to the inherent chemical heterogeneity and structural complexity of the protein surface, designing complex protein assemblies in an anisotropic fashion remains challenging. Here, we describe a self-assembly approach to fabricating protein origami with a networklike structure by designing dual noncovalent interactions on the different positions of a single protein building block. With dimeric proteins as building blocks, 1D protein filaments were constructed by the designed metal coordination at key protein interfaces. Subsequently, the network superstructures were created by the cross-linking of the 1D protein filaments at branch point linkages through the second designed π-π stacking interactions. Notably, upon increasing the protein concentration, the formed protein networks convert into hydrogels with reversible, injectable, and self-healing properties, which have the ability to promote bone regeneration. This strategy could be used to fabricate other protein-based materials with unexplored functions.

Computationally Guided Redesign of a Heme-free Cytochrome with Native-like Structure and Stability

Metals can play key roles in stabilizing protein structures, but ensuring their proper incorporation is a challenge when a metalloprotein is overexpressed in a non-native cellular environment. Here, we have used computational protein design tools to redesign cytochrome b₅₆₂ (cyt b₅₆₂), which relies on the binding of its heme cofactor to achieve its proper fold, into a stable, heme-free protein. The resulting protein, ApoCyt, features only four mutations and no metal-ligand or covalent bonds, yet displays improved stability over cyt b₅₆₂. Mutagenesis studies and X-ray crystal structures reveal that the increase in stability is due to the computationally prescribed mutations, which stabilize the protein fold through a combination of hydrophobic packing interactions, hydrogen bonds, and cation-π interactions. Upon installation of the relevant mutations, ApoCyt is capable of assembling into previously reported, cytochrome-based trimeric and tetrameric assemblies, demonstrating that ApoCyt retains the structure and assembly properties of cyt b₅₆₂. The successful design of ApoCyt therefore enables further functional diversification of cytochrome-based assemblies and demonstrates that structural metal cofactors can be replaced by a small number of well-designed, non-covalent interactions.

Reversible structure transformation between protein nanocages and nanorods controlled by small molecules

Constructing different protein nanostructures by using identical building blocks, while realizing their structural transformation in response to external stimuli, remains a challenge. Here, we fabricated protein nanocages and nanorods by using dimeric TmFtn as a building block and reacting with Mg²⁺/(α, L-lysine) with polymerization degrees of 9 (PLL₉) and 15 (PLL₁₅), respectively. Notably, the reversible shape transformation of these two supramolecular protein architectures with different dimensions can be achievable in response to external stimuli.

Protein interface redesign facilitates the transformation of nanocage building blocks to 1D and 2D nanomaterials

Although various artificial protein nanoarchitectures have been constructed, controlling the transformation between different protein assemblies has largely been unexplored. Here, we describe an approach to realize the self-assembly transformation of dimeric building blocks by adjusting their geometric arrangement. Thermotoga maritima ferritin (TmFtn) naturally occurs as a dimer; twelve of these dimers interact with each other in a head-to-side manner to generate 24-meric hollow protein nanocage in the presence of Ca²⁺ or PEG. By tuning two contiguous dimeric proteins to interact in a fully or partially side-by-side fashion through protein interface redesign, we can render the self-assembly transformation of such dimeric building blocks from the protein nanocage to filament, nanorod and nanoribbon in response to multiple external stimuli. We show similar dimeric protein building blocks can generate three kinds of protein materials in a manner that highly resembles natural pentamer building blocks from viral capsids that form different protein assemblies.

Designed Artificial Protein Heterodimers With Coupled Functions Constructed Using Bio-Orthogonal Chemistry

The formation of protein complexes is central to biology, with oligomeric proteins more prevalent than monomers. The coupling of functionally and even structurally distinct protein units can lead to new functional properties not accessible by monomeric proteins alone. While such complexes are driven by evolutionally needs in biology, the ability to link normally functionally and structurally disparate proteins can lead to new emergent properties for use in synthetic biology and the nanosciences. Here we demonstrate how two disparate proteins, the haem binding helical bundle protein cytochrome b 562 and the β-barrel green fluorescent protein can be combined to form a heterodimer linked together by an unnatural triazole linkage. The complex was designed using computational docking approaches to predict compatible interfaces between the two proteins. Models of the complexes where then used to engineer residue coupling sites in each protein to link them together. Genetic code expansion was used to incorporate azide chemistry in cytochrome b 562 and alkyne chemistry in GFP so that a permanent triazole covalent linkage can be made between the two proteins. Two linkage sites with respect to GFP were sampled. Spectral analysis of the new heterodimer revealed that haem binding and fluorescent protein chromophore properties were retained. Functional coupling was confirmed through changes in GFP absorbance and fluorescence, with linkage site determining the extent of communication between the two proteins. We have thus shown here that is possible to design and build heterodimeric proteins that couple structurally and functionally disparate proteins to form a new complex with new functional properties.

Unfolding cytochromes c-b562 and Rd apo b562

We have analyzed the early stages of unfolding of cytochromes c-b₅₆₂ (PDB ID: 2BC5) and Rd apo b₅₆₂ (PDB ID: 1YYJ). Our geometrical approach proceeds from an analysis of the crystal structure reported for each protein. We quantify, residue-by-residue and region-by-region, the spatial and angular changes in the structure as the protein denatures, and quantify differences that result from the seven residues that differ in the two proteins. Using two independent analyses, one based on spatial metrics and the second on angular metrics, we establish the order of unfolding of the five helices in cyt c-b₅₆₂ and the four helices in the apo protein. For the two helices nearest the N-terminal end of both proteins, the ones in the apo protein unfold first. For the two helices nearest the C-terminal end, the interior helix of the apo protein unfolds first, whereas the terminal helix of the holo protein unfolds first. Excluded-volume effects (repulsive interactions) are minimized in turning regions; the overall range in Δ values is Δ = 36.3 Å³ for cyt c-b₅₆₂ and Δ = 36.6 Å³ for the apo protein, whereas the span for all 20 amino acids is Δ = 167.7 Å³. As our work indicates that the interior helix of cytochrome c-b₅₆₂ is the first to fold, we suggest that this helix protects the heme from misligation, consistent with ultrafast folding over a minimally frustrated funneled landscape.

Chaperone action of a cofactor in protein folding

Previously, we have reported that ATP accelerates the folding and unfolding of Escherichia coli glyceraldehyde-3-phosphate dehydrogenase (GAPDH), which is a glycolytic enzyme utilizing NAD⁺ as a cofactor. Because ATP and NAD⁺ share the ADP moiety, we hypothesized that NAD⁺ also accelerates the folding of GAPDH and that the common structural motif between ATP and NAD⁺ is responsible for the chaperone activity. After confirming that NAD⁺ indeed accelerates the folding of GAPDH, we examined the chaperone activity of the structural fragments of NAD⁺ (ADP, AMP, adenosine, and nicotinamide monophosphate). Our finding showed that ADP and AMP significantly speed up the folding of GAPDH, while adenosine and nicotinamide monophosphate do not. ADP and AMP also dramatically speed up the unfolding of GAPDH by selectively stabilizing a transition state in which GAPDH has a partially unfolded conformation. Similar to the previously reported effect of ATP on the equilibrium unfolding of GAPDH, a partially unfolded intermediate also accumulates in the presence of ADP and AMP. Based on the effect of the structural fragments of NAD⁺ on the folding of GAPDH, we identified that AMP is the structural determinant of the chaperone activity of ATP and NAD⁺ . Also, we propose a plausible model to explain how NAD⁺ accelerates the folding of GAPDH through a stepwise development of molecular interactions with the protein.

My journey in academia: things not on the CV

I am a professor at Chalmers University of Technology in Sweden. I trained in chemistry in Sweden but went to the USA for my postdoc. I remained there for 12 years, being faculty at two American universities, before I returned to Sweden for a professorship in the northern city of Umeå. More recently, I returned to my alma mater Chalmers University of Technology in Gothenburg, where I have taken on senior leadership roles. On paper, my career trajectory looks straightforward, but there are many detrimental aspects and lucky coincidences that are not listed on my CV. Life in academia is never easy, and one is never ‘done’. But working in academia is wonderful, as it provides so much freedom and creativity, including being very accommodating towards having kids. Here, I will describe my own personal journey, with the hope of inspiring young women to follow their own path in academia. Yes, there is still bias against women in academia, but change is happening, and the many benefits of being an academic beat such drawbacks.

Switching protein metalloporphyrin binding specificity by design from iron to fluorogenic zinc

Metalloporphyrins play important roles in areas ranging from biology to nanoscience. Using computational design, we converted metalloporphyrin specificity of cytochrome b₅₆₂ from iron to fluorogenic zinc. The new variant had a near total preference for zinc representing a switch in specificity, which greatly enhanced the negligible aqueous fluorescence of free ZnPP in vitro and in vivo.

Copper-binding energetics of amicyanin in different folding states

Amicyanin is a type I copper protein that mediates electron transfer between methylamine dehydrogenase and cytochrome c-_551i for energy production in Paracoccus denitrificans. Although the Met98 axial ligand of amicyanin has been shown to dictate metal selectivity and specificity during protein folding, the mechanism involved in copper-mediated amicyanin folding is unknown. Here, we kinetically and spectroscopically described reaction steps for incorporating copper into fully and less folded apo-amicyanin and established thermodynamic parameters for two amicyanin folding states. The rate constant for the incorporation of copper into fully folded apo-amicyanin at 25 °C was almost 1.5-fold lower than that for the initial phase of copper addition to the less folded apo-amicyanin. However, the rate constant was 10-fold higher than that of the second phase of copper addition to less folded apo-amicyanin at 25 °C. When overall binding energetic parameters (ΔH° and ΔS°) for the incorporation of copper into fully folded apo-amicyanin were measured by the van't Hoff method and isothermal titration calorimetry, the values were more positive than those determined for less folded apo-amicyanin. This indicates that during amicyanin biogenesis, copper rapidly binds to an unfolded apo-amicyanin active site, inducing protein folding and favorably influencing subsequent organization of copper ligands.

Stabilization of protein structure through π-π interaction in the second coordination sphere of pseudoazurin

Noncovalent, weak interactions in the second coordination sphere of the copper active site of Pseudoazurin (PAz) from Achromobacter cycloclastes were examined using a series of Met16X variants. In this study, the differences in protein stability due to the changes in the nature of the 16th amino acid (Met, Phe, Val, Ile) were investigated by electrospray ionization mass spectrometry (ESI-MS) and far-UV circular dichroism (CD) as a result of acid denaturation. The percentage of native states (folded holo forms) of Met16Phe variants was estimated to be 75% at pH 2.9 although the wild-type (WT), Met16Val and Met16Ile PAz, became completely unfolded. The high stability under acidic conditions is correlated with the result of the active site being stabilized by the aromatic substitution of the Met16 residue. The π-π interaction in the second coordination sphere makes a significant contribution to the stability of active site and the protein matrix.

Domain-swapped cytochrome cb562 dimer and its nanocage encapsulating a Zn-SO4 cluster in the internal cavity

Three domain-swapped cytochrome cb₅₆₂ dimers formed a unique cage structure with a Zn–SO₄ cluster inside the cavity.

Protein nanostructures have been gaining in interest, along with developments in new methods for construction of novel nanostructures. We have previously shown that c-type cytochromes and myoglobin form oligomers by domain swapping. Herein, we show that a four-helix bundle protein cyt cb₅₆₂, with the cyt b₅₆₂ heme attached to the protein moiety by two Cys residues insertion, forms a domain-swapped dimer. Dimeric cyt cb₅₆₂ did not dissociate to monomers at 4 °C, whereas dimeric cyt b₅₆₂ dissociated under the same conditions, showing that heme attachment to the protein moiety stabilizes the domain-swapped structure. According to X-ray crystallographic analysis of dimeric cyt cb₅₆₂, the two helices in the N-terminal region of one protomer interacted with the other two helices in the C-terminal region of the other protomer, where Lys51–Asp54 served as a hinge loop. The heme coordination structure of the dimer was similar to that of the monomer. In the crystal, three domain-swapped cyt cb₅₆₂ dimers formed a unique cage structure with a Zn–SO₄ cluster inside the cavity. The Zn–SO₄ cluster consisted of fifteen Zn²⁺ and seven SO₄^2– ions, whereas six additional Zn²⁺ ions were detected inside the cavity. The cage structure was stabilized by coordination of the amino acid side chains of the dimers to the Zn²⁺ ions and connection of two four-helix bundle units through the conformation-adjustable hinge loop. These results show that domain swapping can be applied in the construction of unique protein nanostructures.

Microdroplet fusion mass spectrometry for fast reaction kinetics

We investigated the fusion of high-speed liquid droplets as a way to record the kinetics of liquid-phase chemical reactions on the order of microseconds. Two streams of micrometer-size droplets collide with one another. The droplets that fused (13 μm in diameter) at the intersection of the two streams entered the heated capillary inlet of a mass spectrometer. The mass spectrum was recorded as a function of the distance x between the mass spectrometer inlet and the droplet fusion center. Fused droplet trajectories were imaged with a high-speed camera, revealing that the droplet fusion occurred approximately within a 500-μm radius from the droplet fusion center and both the size and the speed of the fused droplets remained relatively constant as they traveled from the droplet fusion center to the mass spectrometer inlet. Evidence is presented that the reaction effectively stops upon entering the heated inlet of the mass spectrometer. Thus, the reaction time was proportional to x and could be measured and manipulated by controlling the distance x. Kinetic studies were carried out in fused water droplets for acid-induced unfolding of cytochrome c and hydrogen-deuterium exchange in bradykinin. The kinetics of the former revealed the slowing of the unfolding rates at the early stage of the reaction within 50 μs. The hydrogen-deuterium exchange revealed the existence of two distinct populations with fast and slow exchange rates. These studies demonstrated the power of this technique to detect reaction intermediates in fused liquid droplets with microsecond temporal resolution.

Intrachain contact dynamics in unfolded cytochrome cb562

We have investigated intrachain contact dynamics in unfolded cytochrome cb562 by monitoring heme quenching of excited ruthenium photosensitizers covalently bound to residues along the polypeptide. Intrachain diffusion for chemically denatured proteins proceeds on the microsecond time scale with an upper limit of 0.1 μs. The rate constants exhibit a power-law dependence on the number of peptide bonds between the heme and Ru complex. The power-law exponent of -1.5 is consistent with theoretical models for freely jointed Gaussian chains, but its magnitude is smaller than that reported for several synthetic polypeptides. Contact formation within a stable loop was examined in a His63-heme ligated form of the protein under denaturing conditions. Loop formation accelerated contact kinetics for the Ru66 labeling site, owing to reduction in the length of the peptide separating redox sites. For other labeling sites within the stable loop, quenching rates were modestly reduced compared to the open chain polymer.

Single-Molecule Analysis of Cytochrome c Folding by Monitoring the Lifetime of an Attached Fluorescent Probe

Conformational dynamics of proteins are important for function. However, obtaining information about specific conformations is difficult for samples displaying heterogeneity. Here, time-resolved fluorescence resonance energy transfer is used to characterize the folding of single cytochrome c molecules. In particular, measurements of the fluorescence lifetimes of individual cytochrome c molecules labeled with a single dye that is quenched by energy transfer to the heme were used to monitor conformational transitions of the protein under partially denaturing conditions. These studies indicate significantly more conformational heterogeneity than has been described previously. Importantly, the use of a purified singly-labeled sample made a direct comparison to ensemble data possible. The distribution of lifetimes of single-proteins was compared to the distribution of lifetimes determined from analysis of ensemble lifetime fluorescence data. The results show broad agreement between single-molecule and ensemble data, with a similar range of lifetimes. However, the single-molecule data reveal greater conformational heterogeneity.

Cooperativity, local-nonlocal coupling, and nonnative interactions: principles of protein folding from coarse-grained models

Coarse-grained, self-contained polymer models are powerful tools in the study of protein folding. They are also essential to assess predictions from less rigorous theoretical approaches that lack an explicit-chain representation. Here we review advances in coarse-grained modeling of cooperative protein folding, noting in particular that the Levinthal paradox was raised in response to the experimental discovery of two-state-like folding in the late 1960s, rather than to the problem of conformational search per se. Comparisons between theory and experiment indicate a prominent role of desolvation barriers in cooperative folding, which likely emerges generally from a coupling between local conformational preferences and nonlocal packing interactions. Many of these principles have been elucidated by native-centric models, wherein nonnative interactions may be treated perturbatively. We discuss these developments as well as recent applications of coarse-grained chain modeling to knotted proteins and to intrinsically disordered proteins.

Investigation of an anomalously accelerating substitution in the folding of a prototypical two-state protein

The folding rates of two-state single-domain proteins are generally resistant to small-scale changes in amino acid sequence. For example, having surveyed here over 700 single-residue substitutions in 24 well-characterized two-state proteins, we find that the majority (55%) of these substitutions affect folding rates by less than a factor of 2, and that only 9% affect folding rates by more than a factor of 8. Among those substitutions that significantly affect folding rates, we find that accelerating substitutions are an order of magnitude less common than those that decelerate the process. One of the most extreme outliers in this data set, an arginine-to-phenylalanine substitution at position 48 (R48F) of chymotrypsin inhibitor 2 (CI2), accelerates the protein's folding rate by a factor of 36 relative to that of the wild-type protein and is the most accelerating substitution reported to date in a two-state protein. In order to better understand the origins of this anomalous behavior, we have characterized the kinetics of multiple additional substitutions at this position. We find that substitutions at position 48 in CI2 fall into two distinct classes. The first, comprising residues that ablate the charge of the wild-type arginine but retain the hydrophobicity of its alkane chain, accelerate folding by at least 10-fold. The second class, comprising all other residues, produces folding rates within a factor of two of the wild-type rate. A significant positive correlation between hydrophobicity and folding rate across all of the residues we have characterized at this position suggests that the hydrophobic methylene units of the wild-type arginine play a significant role in stabilizing the folding transition state. Likewise, studies of the pH dependence of the histidine substitution indicate a strong correlation between folding rate and charge state. Thus, mutations that ablate the arginine's positive charge while retaining the hydrophobic contacts of its methylene units tend to dramatically accelerate folding. Previous studies have suggested that arginine 48 plays an important functional role in CI2, which may explain why it is highly conserved despite the anomalously large deceleration it produces in the folding of this protein.

Cleavage site specificity and conformational selection in type I collagen degradation

Excessive degradation of type I collagen is associated with a variety of human diseases such as arthritis, tumor metastasis, and atherosclerosis. Methods that further our understanding of collagenolysis may therefore provide insights into the mechanism of several important disorders. Prior experiments suggest that cleavage of collagen in vitro requires intact full-length collagenase, a multidomain protein containing both a catalytic and a hemopexin-like domain. In this work we demonstrate that type I collagen can be degraded at room temperature, a temperature well below the melting temperature of type I collagen, by collagenase deletion mutants that only contain the catalytic domain of the enzyme. Furthermore, these mutant enzymes hydrolyze the same peptide bond that is recognized by the corresponding full-length enzymes. Hence enzyme specificity at room temperature is achieved without the hemopexin-like domain. We demonstrate that these findings can be explained in light of a conformational selection mechanism that dictates that collagenases preferentially recognize and cleave preformed partially unfolded states of collagen.

Guanidine hydrochloride-induced unfolding of the three heme coordination states of the CO-sensing transcription factor, CooA

CooA is a heme-dependent CO-sensing transcription factor that has three observable heme coordination states. There is some evidence that each CooA heme state has a distinct protein conformation; the goal of this study was to characterize these conformations by measuring their structural stabilities through guanidine hydrochloride (GuHCl) denaturation. By studying the denaturation processes of the Fe(III) state of WT CooA and several variants, we were able to characterize independent unfolding processes for each domain of CooA. This information was used to compare the unfolding profiles of various CooA heme activation states [Fe(III), Fe(II), and Fe(II)-CO] to show that the heme coordination state changes the stability of the effector binding domain. A mechanism consistent with the data predicts that all CooA coordination states and variants undergo unfolding of the DNA-binding domain between 2 and 3 M GuHCl with a free energy of unfolding of approximately 17 kJ/mol, while unfolding of the heme domain is variable and dependent on the heme coordination state. The findings support a model in which changes in heme ligation alter the structural stability of the heme domain and dimer interface but do not alter the stability of the DNA-binding domain. These studies provide evidence that the domains of transcription factors are modular and that allosteric signaling occurs through changes in the relative positions of the protein domains without affecting the structure of the DNA-binding region.

Folding energy landscape of cytochrome cb562

Cytochrome cb(562) is a variant of an Escherichia coli four-helix bundle b-type heme protein in which the porphyrin prosthetic group is covalently ligated to the polypeptide near the terminus of helix 4. Studies from other laboratories have shown that the apoprotein folds rapidly without the formation of intermediates, whereas the holoprotein loses heme before native structure can be attained. Time-resolved fluorescence energy transfer (TRFET) measurements of cytochrome cb(562) refolding triggered using an ultrafast continuous-flow mixer (150 micros dead time) reveal that heme attachment to the polypeptide does not interfere with rapid formation of the native structure. Analyses of the TRFET data produce distributions of Trp-59-heme distances in the protein before, during, and after refolding. Characterization of the moments and time evolution of these distributions provides compelling evidence for a refolding mechanism that does not involve significant populations of intermediates. These observations suggest that the cytochrome b(562) folding energy landscape is minimally frustrated and able to tolerate the introduction of substantial perturbations (i.e., the heme prosthetic group) without the formation of deep misfolded traps.

The Ultimate Speed Limit to Protein Folding Is Conformational Searching

More than a dozen proteins are known to be ultrafast folders. In addition to being fast, their kinetics is unusual. Like traditional rate processes, fast folding proteins have activation barriers at low temperatures, but unlike traditional processes, they have negative activation energies at high temperatures. We develop a model of ultrafast folders that joins a macroscopic mass-action model with a microscopic energy landscape description; we call it the Thruway Search Model. We find good agreement with experimental rates and equilibria on 13 ultrafast folders. The observed folding rates are found to be proportional to the number of microscopic folding routes: fast-folding proteins have more parallel microscopic routes on energy landscapes. At high temperatures, where traditional barriers are small, the remaining bottleneck is a search through denatured conformations to find thruway routes to the native state. Negative activation arises because increasing temperature expands the denatured ensemble, broadening the search, slowing the folding to the native state. We find that the upper estimate of the free energy barriers are positive but small, as little as 0.5 kT.

Probing the cytochrome c' folding landscape

The folding kinetics of R. palustris cytochrome c' (cyt c') have been monitored by heme absorption and native Trp72 fluorescence at pH 5. The Trp72 fluorescence burst signal suggests early compaction of the polypeptide ensemble. Analysis of heme transient absorption spectra reveals deviations from two-state behavior, including a prominent slow phase that is accelerated by the prolyl isomerase cyclophilin. A nonnative proline configuration (Pro21) likely interferes with the formation of the helical bundle surrounding the heme.

Folding rate prediction using n-order contact distance for proteins with two- and three-state folding kinetics

It is a challenging task to understand the relationship between sequences and folding rates of proteins. Previous studies are found that one of contact order (CO), long-range order (LRO), total contact distance (TCD), chain topology parameter (CTP), and effective length (Leff) has a significant correlation with folding rate of proteins. In this paper, we introduce a new parameter called n-order contact distance (nOCD) and use it to predict folding rate of proteins with two- and three-state folding kinetics. A good linear correlation between the folding rate logarithm lnkf and nOCD with n=1.2, alpha=0.6 is found for two-state folders (correlation coefficient is -0.809, P-value<0.0001) and n=2.8, alpha=1.5 for three-state folders (correlation coefficient is -0.816, P-value<0.0001). However, this correlation is completely absent for three-state folders with n=1.2, alpha=0.6 (correlation coefficient is 0.0943, P-value=0.661) and for two-state folders with n=2.8, alpha=1.5 (correlation coefficient is -0.235, P-value=0.2116). We also find that the average number of contacts per residue Pm in the interval of m for two-state folders is smaller than that for three-state folders. The probability distribution P(gamma) of residue having gamma pairs of contacts fits a Gaussian distribution for both two- and three-state folders. We observe that the correlations between square radius of gyration S2 and number of residues for two- and three-state folders are both good, and the correlation coefficient is 0.908 and 0.901, and the slope of the fitting line is 1.202 and 0.795, respectively. Maybe three-state folders are more compact than two-state folders. Comparisons with nTCD and nCTP are also made, and it is found that nOCD is the best one in folding rate prediction.

Folding of Desulfovibrio desulfuricans flavodoxin is accelerated by cofactor fly-casting

Folding of cofactor-binding proteins involves ligand binding in addition to polypeptide folding. We here assess the kinetic folding/binding landscape for Desulfovibrio desulfuricans flavodoxin that coordinates an FMN cofactor. The apo-form folds in a two-step process involving a burst-phase intermediate. Studies on Tyr98Ala and Trp60Ala variants reveal that these aromatics-that stack with the FMN in the holo-form-are not participating in the apo-protein folding pathway. However, these residues are essential for FMN interactions with the unfolded protein during refolding of holo-flavodoxin. Unfolding of wild-type holo-flavodoxin is coupled to FMN dissociation whereas for Tyr98Ala and Trp60Ala holo-variants, FMN dissociates before polypeptide unfolding. Both variants refold as apo-proteins before FMN rebinds. In sharp contrast, refolding of unfolded wild-type holo-flavodoxin is over an order of magnitude faster than that of the apo-form, the pathway does not include a burst-phase intermediate, and the speed is independent of FMN excess ratio. These observations demonstrate that FMN binds rapidly to the unfolded polypeptide and guides folding straight to the native state. As this path to functional D. desulfuricans holo-flavodoxin is faster than if the cofactor binds to pre-folded apo-protein, this is one of few examples where molecular recognition via a "fly-casting" mechanism is kinetically favored.

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.

Effects of Heme on the Structure of the Denatured State and Folding Kinetics of Cytochrome b562

Heme-linked proteins, such as cytochromes, are popular subjects for protein folding studies. There is the underlying question of whether the heme affects the structure of the denatured state by cross-linking it and forming other interactions, which would perturb the folding pathway. We have studied wild-type and mutant cytochrome b562 from Escherichia coli, a 106 residue four-alpha-helical bundle. The holo protein apparently refolds with a half-life of 4 micros in its ferrous state. We have analysed the folding of the apo protein using continuous-flow fluorescence as well as stopped-flow fluorescence and CD. The apo protein folded much more slowly with a half-life of 270 micros that was unaffected by the presence of exogenous heme. We examined the nature of the denatured states of both holo and apo proteins by NMR methods over a range of concentrations of guanidine hydrochloride. The starting point for folding of the holo protein in concentrations of denaturant around the denaturation transition was a highly ordered native-like species with heme bound. Fully denatured holo protein at higher concentrations of denaturant consisted of denatured apo protein and free heme. Our results suggest that the very fast folding species of denatured holo protein is in a compact state, whereas the normal folding pathway from fully denatured holo protein consists of the slower folding of the apo protein followed by the binding of heme. These data should be considered in the analysis of folding of heme proteins.

Native geometry and the dynamics of protein folding

In this paper, we investigate the role of native geometry on the kinetics of protein folding based on simple lattice models and Monte Carlo simulations. Results obtained within the scope of the Miyazawa-Jernigan indicate the existence of two dynamical folding regimes depending on the protein chain length. For chains larger than 80 amino acids, the folding performance is sensitive to the native state's conformation. Smaller chains, with less than 80 amino acids, fold via two-state kinetics and exhibit a significant correlation between the contact order parameter and the logarithmic folding times. In particular, chains with N=48 amino acids were found to belong to two broad classes of folding, characterized by different cooperativity, depending on the contact order parameter. Preliminary results based on the Go model show that the effect of long-range contact interaction strength in the folding kinetics is largely dependent on the native state's geometry.

FMN binding and unfolding of Desulfovibrio desulfuricans flavodoxin: "hidden" intermediates at low denaturant concentrations

The flavin mononucleotide (FMN) cofactor in Desulfovibrio desulfuricans flavodoxin stays associated with the polypeptide upon guanidine hydrochloride (GuHCl) induced unfolding. Using isothermal titration calorimetry (ITC), we determined the affinity of FMN for the flavodoxin polypeptide as a function of both urea and GuHCl concentrations (pH 7, 25 degrees C). The FMN affinity for folded and GuHCl-unfolded flavodoxin differs 10-fold, which is in agreement with the difference in thermodynamic stability between the apo- and holo-forms. In contrast, the urea-unfolded protein does not interact with FMN and equilibrium unfolding of holo-flavodoxin in urea results in FMN dissociation prior to polypeptide unfolding. ANS-binding, near-UV circular dichroism (CD), acrylamide quenching and FMN-emission experiments reveal the presence of native-like intermediates, not detected by far-UV CD and aromatic fluorescence detection methods, in low concentrations of both denaturants. Time-resolved experiments show that FMN binding is fastest at GuHCl concentrations where the native-like intermediate species is populated.

The folding of spectrin domains II: phi-value analysis of R16

Studies on the folding of helical proteins have shown a wide range of different mechanisms and highlighted the importance of helical propensity as a factor in determining folding mechanism. Here, we contribute to this interesting field with the protein engineering phi-value analysis of the 16th domain of chicken brain alpha-spectrin, R16. The fortuitous curvature seen in the unfolding arm of the chevron plot allows us to investigate both early and late events in folding. R16 is the first two-state helical protein for which this has been possible.

Role of Cofactors in Folding of the Blue-Copper Protein Azurin

Many proteins in living cells coordinate cofactors, such as metal ions, to attain their activity. Since the cofactors in such cases often can interact with their corresponding unfolded polypeptides in vitro, it is important to unravel how cofactors modulate protein folding. In this review, I will discuss the role of cofactors in folding of the blue-copper protein Pseudomonas aeruginosa azurin. In the case of both copper (Cu(II) and Cu(I)) and zinc (Zn(II)), the metal can bind to unfolded azurin. The residues involved in copper (Cu(II) and Cu(I)) coordination in the unfolded state have been identified as Cys112, His117, and Met121. The affinities of Cu(II), Cu(I), and Zn(II) are all higher for the folded than for the unfolded azurin polypeptide, resulting in metal stabilization of the native state as compared to the stability of apo-azurin. Cu(II), Zn(II), and several apo forms of azurin all fold in two-state kinetic reactions with roughly identical polypeptide-folding speeds. This suggests that the native-state beta-barrel topology, not cofactor interactions or thermodynamic stability, determines azurin's folding barrier. Nonetheless, copper binds much more rapidly (i.e., 4 orders of magnitude) to unfolded azurin than to folded azurin. Therefore, the fastest route to functional azurin is through copper binding before polypeptide folding; this sequence of events may be the relevant biological pathway.

Fast and faster: a designed variant of the B-domain of protein A folds in 3 microsec

We have introduced the mutation glycine 29 to alanine, designed to increase the rate of protein folding, into the B-domain of protein A (BdpA). From NMR lineshape analysis, we find the G29A mutation increases the folding rate constant by threefold; the folding time is 3 microsec. Although wild-type BdpA folds extremely fast, simple-point mutations can still speed up the folding; thus, the folding rate is not evolutionarily maximized. The short folding time of G29A BdpA (the shortest time yet reported) makes it an attractive candidate for an all-atom molecular dynamics simulation that could potentially show a complete folding reaction starting from an extended chain. We also constructed a fluorescent variant of BdpA by mutating phenylalanine 13 to tryptophan, allowing fluorescence-based time-resolved temperature-jump measurements. Temperature jumps and NMR complement each other, and give a very complete picture of the folding kinetics.

The protein folding 'speed limit'

How fast can a protein possibly fold? This question has stimulated experimentalists to seek fast folding proteins and to engineer them to fold even faster. Proteins folding at or near the speed limit are prime candidates for all-atom molecular dynamics simulations. They may also have no free energy barrier, allowing the direct observation of intermediate structures on the pathways from the unfolded to the folded state. Both experimental and theoretical approaches predict a speed limit of approximately N/100micros for a generic N-residue single-domain protein, with alpha proteins folding faster than beta or alphabeta. The predicted limits suggest that most known ultrafast folding proteins can be engineered to fold more than ten times faster.

Quantifying the roughness on the free energy landscape: entropic bottlenecks and protein folding rates

The prediction of protein folding rates and mechanisms is currently of great interest in the protein folding community. A close comparison between theory and experiment in this area is promising to advance our understanding of the physical-chemical principles governing the folding process. The delicate interplay of entropic and energetic/enthalpic factors in the protein free energy regulates the details of this complex reaction. In this article, we propose the use of topological descriptors to quantify the amount of heterogeneity in the configurational entropy contribution to the free energy. We apply the procedure to a set of 16 two-state folding proteins. The results offer a clean and simple theoretical explanation for the experimentally measured folding rates and mechanisms, in terms of the intrinsic entropic roughness along the populated folding routes on the protein free energy landscape.

Using protein folding rates to test protein folding theories

The fastest simple, kinetically two-state protein folds a million times more rapidly than the slowest. Here we review many recent theories of protein folding kinetics in terms of their ability to qualitatively rationalize, if not quantitatively predict, this fundamental experimental observation.

Refolding processes of cytochrome P450cam from ferric and ferrous acid forms to the native conformation. Formations of folding intermediates with non-native heme coordination state

Changes in heme coordination state and protein conformation of cytochrome P450(cam) (P450(cam)), a b-type heme protein, were investigated by employing pH jump experiments coupled with time-resolved optical absorption, fluorescence, circular dichroism, and resonance Raman techniques. We found a partially unfolded form (acid form) of ferric P450(cam) at pH 2.5, in which a Cys(-)-heme coordination bond in the native conformation was ruptured. When the pH was raised to pH 7.5, the acid form refolded to the native conformation through a distinctive intermediate. Formations of similar acid and intermediate forms were also observed for ferrous P450(cam). Both the ferric and ferrous forms of the intermediate were found to have an unidentified axial ligand of the heme at the 6th coordination sphere, which is vacant in the high spin ferric and ferrous forms at the native conformation. For the ferrous form, it was also indicated that the 5th axial ligand is different from the native cysteinate. The folding intermediates identified in this study demonstrate occurrences of non-native coordination state of heme during the refolding processes of the large b-type heme protein, being akin to the well known folding intermediates of cytochromes c, in which c-type heme is covalently attached to a smaller protein.

Folding and form: insights from lattice simulations

Monte Carlo simulations of a Miyazawa-Jernigan lattice-polymer model indicate that, depending on the native structure's geometry, the model exhibits two broad classes of folding mechanisms for two-state folders. Folding to native structures of low contact order is driven by backbone distance and is characterized by a progressive accumulation of structure towards the native fold. By contrast, folding to high contact order targets is dominated by intermediate stage contacts not present in the native fold, yielding a more cooperative folding process.

Ultrafast folding of alpha3D: a de novo designed three-helix bundle protein

Here, we describe the folding/unfolding kinetics of alpha3D, a small designed three-helix bundle. Both IR temperature jump and ultrafast fluorescence mixing methods reveal a single-exponential process consistent with a minimal folding time of 3.2 +/- 1.2 micros (at approximately 50 degrees C), indicating that a protein can fold on the 1- to 5-micros time scale. Furthermore, the single-exponential nature of the relaxation indicates that the prefactor for transition state (TS)-folding models is probably >or=1 (micros)-1 for a protein of this size and topology. Molecular dynamics simulations and IR spectroscopy provide a molecular rationale for the rapid, single-exponential folding of this protein. alpha3D shows a significant bias toward local helical structure in the thermally denatured state. The molecular dynamics-simulated TS ensemble is highly heterogeneous and dynamic, allowing access to the TS via multiple pathways.

The Nature of the Rate-limiting Steps in the Refolding of the Cofactor-dependent Protein Aspartate Aminotransferase

The refolding of mitochondrial aspartate aminotransferase (mAAT; EC 2.6.1.1) has been studied following unfolding in 6 m guanidine hydrochloride for different periods of time. Whereas reactivation of equilibrium-unfolded mAAT is sigmoidal, reactivation of the short term unfolded protein displays a double exponential behavior consistent with the presence of fast and slow refolding species. The amplitude of the fast phase decreases with increasing unfolding times (k approximately 0.75 min(-1) at 20 degrees C) and becomes undetectable at equilibrium unfolding. According to hydrogen exchange and stopped-flow intrinsic fluorescence data, unfolding of mAAT appears to be complete in less than 10 s, but hydrolysis of the Schiff base linking the coenzyme pyridoxal 5'-phosphate (PLP) to the polypeptide is much slower (k approximately 0.08 min(-1)). This implies the existence in short term unfolded samples of unfolded species with PLP still attached. However, since the disappearance of the fast refolding phase is about 10-fold faster than the release of PLP, the fast refolding phase does not correspond to folding of the coenzyme-containing molecules. The fast refolding phase disappears more rapidly in the pyridoxamine and apoenzyme forms of mAAT, both of which lack covalently attached cofactor. Thus, bound PLP increases the kinetic stability of the fast refolding unfolding intermediates. Conversion between fast and slow folding forms also takes place in an early folding intermediate. The presence of cyclophilin has no effect on the reactivation of either equilibrium or short term unfolded mAAT. These results suggest that proline isomerization may not be the only factor determining the slow refolding of this cofactor-dependent protein.

Approaching the speed limit for Greek Key β-barrel formation: transition-state movement tunes folding rate of zinc-substituted azurin

Azurin is a blue-copper protein with a beta-barrel structure of Greek Key topology. In vitro, copper can be substituted with zinc without change in protein structure. We here analyze the kinetic folding behavior of zinc-substituted Pseudomonas aeruginosa azurin. Our findings can be summarized in three key conclusions: first, zinc remains strongly bound to the polypeptide upon unfolding, suggesting that the cofactor may bind to the protein before polypeptide folding in vivo. Second, the semi-logarithmic plot of folding and unfolding rates for zinc-substituted azurin as a function of denaturant concentration exhibits curvature due to a changing transition-state structure. Third, the extrapolated folding speed in water for zinc-substituted azurin is similar to that of other proteins with the same topology, implying that there is a speed limit that can be modulated by stability-driven transition-state movement for formation of beta-barrel structures with Greek Key topology.

Contact order revisited: influence of protein size on the folding rate

Guided by the recent success of empirical model predicting the folding rates of small two-state folding proteins from the relative contact order (CO) of their native structures, by a theoretical model of protein folding that predicts that logarithm of the folding rate decreases with the protein chain length L as L(2/3), and by the finding that the folding rates of multistate folding proteins strongly correlate with their sizes and have very bad correlation with CO, we reexamined the dependence of folding rate on CO and L in attempt to find a structural parameter that determines folding rates for the totality of proteins. We show that the Abs_CO = CO x L, is able to predict rather accurately folding rates for both two-state and multistate folding proteins, as well as short peptides, and that this Abs_CO scales with the protein chain length as L(0.70 +/- 0.07) for the totality of studied single-domain proteins and peptides.