Cofactor effects on the protein folding reaction: acceleration of alpha-lactalbumin refolding by metal ions

About 30% of proteins require cofactors for their proper folding. The effects of cofactors on the folding reaction have been investigated with alpha-lactalbumin as a model protein and metal ions as cofactors. Metal ions accelerate the refolding of alpha-lactalbumin by lessening the energy barrier between the molten globule state and the transition state, mainly by decreasing the difference of entropy between the two states. These effects are linked to metal ion binding to the protein in the native state. Hence, relationships between the metal affinities for the intermediate states and those for the native state are observed. Some residual specificity for the calcium ion is still observed in the molten globule state, this specificity getting closer in the transition state to that of the native state. The comparison between kinetic and steady-state data in association with the Phi value method indicates the binding of the metal ions on the unfolded state of alpha-lactalbumin. Altogether, these results provide insight into cofactor effects on protein folding. They also suggest new possibilities to investigate the presence of residual native structures in the unfolded state of protein and the effects of such structures on the protein folding reaction and on protein stability.     

Extending the folding nucleus of ubiquitin with an independently folding beta-hairpin finger: hurdles to rapid folding arising from the stabilisation of local interactions

The N-terminal beta-hairpin sequence of ubiquitin has been implicated as a folding nucleation site. To extend and stabilise the ubiquitin folding nucleus, we have inserted an autonomously folding 14-residue peptide sequence beta4 which in isolation forms a highly populated beta-hairpin (>70%) stabilised by local interactions. NMR structural analysis of the ubiquitin mutant (Ubeta4) shows that the hairpin finger is fully structured and stabilises ubiquitin by approximately 8kJmol(-1). Protein engineering and kinetic (phi(F)-value) analysis of a series of Ubeta4 mutants shows that the hairpin extension of Ubeta4 is also significantly populated in the transition state (phi(F)-values >0.7) and has the effect of templating the formation of native contacts in the folding nucleus of ubiquitin. However, at low denaturant concentrations the chevron plot of Ubeta4 shows a small deviation from linearity (roll-over effect), indicative of the population of a compact collapsed state, which appears to arise from over-stabilisation of local interactions. Destabilising mutations within the native hairpin sequence and within the engineered hairpin extension, but not elsewhere, eliminate this non-linearity and restore apparent two-state behaviour. The pitfall to stabilising local interactions is to present hurdles to the rapid and efficient folding of small proteins down a smooth folding funnel by trapping partially folded or misfolded states that must unfold or rearrange before refolding.     

Disulfide bonds and protein folding

The applications of disulfide-bond chemistry to studies of protein folding, structure, and stability are reviewed and illustrated with bovine pancreatic ribonuclease A (RNase A). After surveying the general properties and advantages of disulfide-bond studies, we illustrate the mechanism of reductive unfolding with RNase A, and discuss its application to probing structural fluctuations in folded proteins. The oxidative folding of RNase A is then described, focusing on the role of structure formation in the regeneration of the native disulfide bonds. The development of structure and conformational order in the disulfide intermediates during oxidative folding is characterized. Partially folded disulfide species are not observed, indicating that disulfide-coupled folding is highly cooperative. Contrary to the predictions of "rugged funnel" models of protein folding, misfolded disulfide species are also not observed despite the potentially stabilizing effect of many nonnative disulfide bonds. The mechanism of regenerating the native disulfide bonds suggests an analogous scenario for conformational folding. Finally, engineered covalent cross-links may be used to assay for the association of protein segments in the folding transition state, as illustrated with RNase A.     

Interpretation of protein folding psi values

The structural characterization of transition states is essential for understanding the mechanism of protein folding. Analyzing the effect of mutations on protein stability and folding kinetics in phi-value analysis is commonly used to gain information about the presence of side-chain interactions in transition states. Recently, specific binding of ligands to engineered binding sites was applied to monitor the formation of local structures in transition states (psi analysis). A surprising result from psi analysis was the presence of parallel folding pathways in all reported studies and a major discrepancy between phi and psi values measured in the same protein. Here, we show that psi values cannot be analyzed in the same way as other rate-equilibrium free energy relationships due to the involvement of bimolecular reactions that may have different dissociation constants for the native, unfolded and transition state. As a consequence, psi values reflect the relative binding energy (kappa) of the transition state only for the extreme values of kappa=0 or kappa=1. In all other cases, non-linear rate-equilibrium free-energy relationships (Leffler plots) are observed. This apparently indicates the presence of parallel folding pathways even if folding occurs over a single homogeneous transition state. Consequently, the results from Leffler plots do not yield information about the structural properties of the transition state. This explains the lack of agreement between results from psi analysis and other methods used to characterize protein folding transition states. We further show that the same considerations apply for the analysis of the effect of pH on protein folding.     

Effects of tetracycline and spectinomycin on the tertiary structure of ribosomal RNA in the Escherichia coli 30 S ribosomal subunit.

Structural analysis of the 16 S rRNA in the 30 S subunit and 70 S ribosome in the presence of ribosome-specific antibiotics was performed to determine whether they produced rRNA structural changes that might provide further insight to their action. An UV cross-linking procedure that determines the pattern and frequency of intramolecular 16 S RNA cross-links was used to detect differences reflecting structural changes. Tetracycline and spectinomycin have specific effects detected by this assay. The presence of tetracycline inhibits the cross-link C967xC1400 completely, increases the frequency of cross-link C1402x1501 twofold, and decreases the cross-link G894xU244 by one-half without affecting other cross-links. Spectinomycin reduces the frequency of the cross-link C934xU1345 by 60% without affecting cross-linking at other sites. The structural changes occur at concentrations at which the antibiotics exert their inhibitory effects. For spectinomycin, the apparent binding site and the affected cross-linking site are distant in the secondary structure but are close in tertiary structure in several recent models, indicating a localized effect. For tetracycline, the apparent binding sites are significantly separated in both the secondary and the three-dimensional structures, suggesting a more regional effect.     

Crystallographic study of a site-specifically cross-linked protein complex with a genetically incorporated photoreactive amino acid

The benzophenone photophore is widely used to photo-cross-link macromolecules. Recent developments in genetic code expansion have allowed the biosynthesis of proteins with p-benzoyl-L-phenylalanine (pBpa) at defined sites, for covalent bonding with interacting proteins. However, the structure of a photo-cross-linked protein complex had not been revealed, and thus neither the actual structure of the "photobridge" in a complex nor the influence of this covalent bridge on the overall complex structure was known. In this study, we determine the crystal structure of the cross-linked complex of the liver oncoprotein gankyrin and the C-terminal domain of S6 proteasomal protein (S6C), at 2.05 ? resolution. First, the photoreactive amino acid was separately incorporated into gankyrin at 16 sites on the protein surface, and two variants that efficiently formed a covalent bond with S6C were found. The yield of one of the cross-linked products, with pBpa in place of Arg85 in gankyrin, was maximized for crystallization via optimization of the duration of complex exposure to 365 nm light. The structure revealed that the carbonyl group of the benzophenone of pBpa85 formed a covalent bond exclusively with the Cγ atom of Glu356 in S6C, showing the high selectivity of formation of cross-links by pBpa. In addition, the cross-linked structure exhibited little structural distortion from the native complex structure. Our results demonstrated that cross-linking with site-specifically incorporated pBpa preserves the native binding mode and is useful for probing protein-protein interactions.     

In vitro protein folding by E. coli ribosome: unfolded protein splitting 70S to interact with 50S subunit

Folding of unfolded protein on Escherichia coli 70S ribosome is accompanied by rapid dissociation of the ribosome into 50S and 30S subunits. The dissociation rate of 70S ribosome with unfolded protein is much faster than that caused by combined effect of translation and polypeptide release factors known to be involved in the dissociation of ribosome into subunits. The protein then reaches a “folding competent” state on 50S and is released to take up native conformation by itself. Release before attaining the folding competent state or prevention of release by cross-linking it with ribosome, would not allow the protein to get back to its native conformation.     

The folding transition state of protein L is extensive with nonnative interactions (and not small and polarized)

Progress in understanding protein folding relies heavily upon an interplay between experiment and theory. In particular, readily interpretable experimental data that can be meaningfully compared to simulations are required. According to standard mutational ? analysis, the transition state for Protein L contains only a single hairpin. However, we demonstrate here using ψ analysis with engineered metal ion binding sites that the transition state is extensive, containing the entire four-stranded β sheet. Underreporting of the structural content of the transition state by ? analysis also occurs for acyl phosphatase [Pandit, A. D., Jha, A., Freed, K. F. &; Sosnick, T. R., (2006). Small proteins fold through transition states with native-like topologies. J. Mol. Biol. 361, 755–770], ubiquitin [Sosnick, T. R., Dothager, R. S. &; Krantz, B. A., (2004). Differences in the folding transition state of ubiquitin indicated by ? and ψ analyses. Proc. Natl Acad. Sci. USA 101, 17377–17382] and BdpA [Baxa, M., Freed, K. F. &; Sosnick, T. R., (2008). Quantifying the structural requirements of the folding transition state of protein A and other systems. J. Mol. Biol. 381, 1362–1381]. The carboxy-terminal hairpin in the transition state of Protein L is found to be nonnative, a significant result that agrees with our Protein Data Bank-based backbone sampling and all-atom simulations. The nonnative character partially explains the failure of accepted experimental and native-centric computational approaches to adequately describe the transition state. Hence, caution is required even when an apparent agreement exists between experiment and theory, thus highlighting the importance of having alternative methods for characterizing transition states.     

Small proteins fold through transition states with native-like topologies

The folding pathway of common-type acyl phosphatase (ctAcP) is characterized using psi-analysis, which identifies specific chain-chain contacts using bi-histidine (biHis) metal-ion binding sites. In the transition state ensemble (TSE), the majority of the protein is structured with a near-native topology, only lacking one beta-strand and an alpha-helix. psi-Values are zero or unity for all sites except one at the amino terminus of helix H2. This fractional psi-value remains unchanged when three metal ions of differing coordination geometries are used, indicating this end of the helix experiences microscopic heterogeneity through fraying in the TSE. Ubiquitin, the other globular protein characterized using psi-analysis, also exhibits a single consensus TSE structure. Hence, the TSE of both proteins have converged to a single configuration, albeit one that contains some fraying at the periphery. Models of the TSE of both proteins are created using all-atom Langevin dynamics simulations using distance constraints derived from the experimental psi-values. For both proteins, the relative contact order of the TS models is approximately 80% of the native value. This shared value viewed in the context of the known correlation between contact order and folding rates, suggests that other proteins will have a similarly high fraction of the native contact order. This constraint greatly limits the range of possible configurations at the rate-limiting step.     

Snapshots of a dynamic folding nucleus in zinc-substituted Pseudomonas aeruginosa azurin

Zinc-substituted Pseudomonas aeruginosa azurin folds in two-state equilibrium and kinetic reactions. In the unfolded state, the zinc ion remains bound to the unfolded polypeptide via two native-state ligands (His117 and Cys112). The significantly curved Chevron plot for zinc-substituted azurin was earlier ascribed to movement of the folding-transition state. At low concentrations of denaturant, the transition state occurs early in the folding reaction (low Tanford beta-value), whereas at high-denaturant concentration, it moves closer to the native structure (high Tanford beta-value). Here, we use this movement to track the formation and growth of zinc-substituted azurin's folding nucleus with atomic resolution using protein engineering. The average phi (phi) value for 17 positions (covering all secondary-structure elements) goes from 0.25 in 0 M GuHCl (beta approximately 0.46) to 0.76 in 4 M GuHCl (beta approximately 0.86); a phi-value of 1 or 0 indicates native-like or unfolded-like interactions, respectively. Analysis of individual phi-values reveals a delocalized nucleus where structure condenses around a leading density centered on Leu50 in the core. The diffuse moving transition state for zinc-substituted azurin is in sharp contrast to the fixed polarized folding nucleus observed for apo-azurin. The dramatic difference in apparent kinetic behavior for the two forms of azurin can be rationalized as a minor alteration on a common free-energy profile that exhibits a broad activation barrier.     

Characterization of the molten globule conformation of V26A ubiquitin by far-UV circular dichroic spectroscopy and amide hydrogen/deuterium exchange

The molten globular conformation of V26A ubiquitin (valine to alanine mutation at residue 26) was studied by nuclear magnetic resonance spectroscopy in conjunction with amide hydrogen/deuterium exchange. Most of the amide protons that are involved in the native secondary structures were observed to be protected in the molten globule state with the protection factors from 1.2 to 6.7. These protection factors are about 2 to 6 orders of magnitude smaller than those of the native state. These observations indicate that V26A molten globule has native-like backbone structure with marginal stability. The comparison of amide protection factors of V26A ubiquitin molten globule state with those of initial collapsed state of the wild type ubiquitin suggests that V26A ubiquitin molten globule state is located close to unfolded state in the folding reaction coordinate. It is considered that V26A ubiquitin molten globule is useful model to study early events in protein folding reaction.     

Unfolding pathway of CotA-laccase and the role of copper on the prevention of refolding through aggregation of the unfolded state

Copper is a redox-active metal and the main player in electron transfer reactions occurring in multicopper oxidases. The role of copper in the unfolding pathway and refolding of the multicopper oxidase CotA laccase in vitro was solved using double-jump stopped-flow experiments. Unfolding of apo- and holo-CotA was described as a three-state process with accumulation of an intermediate in between the native and unfolded state. Copper stabilizes the native holo-CotA but also the intermediate state showing that copper is still bound to this state. Also, copper binds to unfolded holo-CotA in a non-native coordination promoting CotA aggregation and preventing refolding to the native structure. These results gather information on unfolding/folding pathways of multicopper oxidases and show that copper incorporation in vivo should be a tight controlled process as copper binding to the unfolded state under native conditions promotes protein aggregation.     

Engineering diverse changes in beta-turn propensities in the N-terminal beta-hairpin of ubiquitin reveals significant effects on stability and kinetics but a robust folding transition state

Using the N-terminal 17-residue beta-hairpin of ubiquitin as a "host" for mutational studies, we have investigated the influence of the beta-turn sequence on protein stability and folding kinetics by replacing the native G-bulged turn (TLTGK) with more flexible analogues (TG3K and TG5K) and a series of four-residue type I' beta-turn sequences, commonly found in beta-hairpins. Although a statistical analysis of type I' turns demonstrates residue preferences at specific sites, the frequency of occurrence appears to only broadly correlate with experimentally determined protein stabilities. The subsequent engineering of context-dependent non-native tertiary contacts involving turn residues is shown to produce large changes in stability. Relatively few point mutations have been described that probe secondary structure formation in ubiquitin in a manner that is independent of tertiary contacts. To this end, we have used the more rigorous rate-equilibrium free energy relationship (Leffler analysis), rather than the two-point phi value analysis, to show for a family of engineered beta-turn mutants that stability (range of approximately 20 kJ/mol) and folding kinetics (190-fold variation in refolding rate) are linearly correlated (alpha(f) = 0.74 +/- 0.08). The data are consistent with a transition state that is robust with regard to a wide range of statistically favored and disfavored beta-turn mutations and implicate a loosely assembled beta-hairpin as a key template in transition state stabilization with the beta-turn playing a central role.     

Autonomous folding of a peptide corresponding to the N-terminal beta-hairpin from ubiquitin.

The N-terminal 17 residues of ubiquitin have been shown by 1H NMR to fold autonomously into a beta-hairpin structure in aqueous solution. This structure has a specific, native-like register, though side-chain contacts differ in detail from those observed in the intact protein. An autonomously folding hairpin has previously been identified in the case of streptococcal protein G, which is structurally homologous with ubiquitin, but remarkably, the two are not in topologically equivalent positions in the fold. This suggests that the organization of folding may be quite different for proteins sharing similar tertiary structures. Two smaller peptides have also been studied, corresponding to the isolated arms of the N-terminal hairpin of ubiquitin, and significant differences from simple random coil predictions observed in the spectra of these subfragments, suggestive of significant limitation of the backbone conformational space sampled, presumably as a consequence of the strongly beta-structure favoring composition of the sequences. This illustrates the ability of local sequence elements to express a propensity for beta-structure even in the absence of actual sheet formation. Attempts were made to estimate the population of the folded state of the hairpin, in terms of a simple two-state folding model. Using published "random coil" values to model the unfolded state, and values derived from native ubiquitin for the putative unique, folded state, it was found that the apparent population varied widely for different residues and with different NMR parameters. Use of the spectra of the subfragment peptides to provide a more realistic model of the unfolded state led to better agreement in the estimates that could be obtained from chemical shift and coupling constant measurements, while making it clear that some other approaches to population estimation could not give meaningful results, because of the tendency to populate the beta-region of conformational space even in the absence of the hairpin structure.     

Effects of the difference in the unfolded-state ensemble on the folding of Escherichia coli dihydrofolate reductase

The unfolded state of a protein is an ensemble of a large number of conformations ranging from fully extended to compact structures. To investigate the effects of the difference in the unfolded-state ensemble on protein folding, we have studied the structure, stability, and folding of "circular" dihydrofolate reductase (DHFR) from Escherichia coli in which the N and C-terminal regions are cross-linked by a disulfide bond, and compared the results with those of disulfide-reduced "linear" DHFR. Equilibrium studies by circular dichroism, difference absorption spectra, solution X-ray scattering, and size-exclusion chromatography show that whereas the native structures of both proteins are essentially the same, the unfolded state of circular DHFR adopts more compact conformations than the unfolded state of the linear form, even with the absence of secondary structure. Circular DHFR is more stable than linear DHFR, which may be due to the decrease in the conformational entropy of the unfolded state as a result of circularization. Kinetic refolding measurements by stopped-flow circular dichroism and fluorescence show that under the native conditions both proteins accumulate a burst-phase intermediate having the same structures and both fold by the same complex folding mechanism with the same folding rates. Thus, the effects of the difference in the unfolded state of circular and linear DHFRs on the refolding reaction are not observed after the formation of the intermediate. This suggests that for the proteins with close termini in the native structure, early compaction of a protein molecule to form a specific folding intermediate with the N and C-terminal regions in close proximity is a crucial event in folding. If there is an enhancement in the folding reflecting the reduction in the breadth of the unfolded-state ensemble for circular DHFR, this acceleration must occur in the sub-millisecond time-range.     

Kinetic analysis of folding and unfolding the 56 amino acid IgG-binding domain of streptococcal protein G.

The 56 amino acid B domain of protein G (GB) is a stable globular folding unit with no disulfide cross-links. The physical properties of GB offer extraordinary flexibility for evaluating the energetics of the folding reaction. The protein is monomeric and very soluble in both folded and unfolded forms. The folding reaction has been previously examined by differential scanning calorimetry (Alexander et al., 1992) and found to exhibit two-state unfolding behavior over a wide pH range with an unfolding transition near 90 degrees C (GB1) at neutral pH. Here, the kinetics of folding and unfolding two naturally occurring versions of GB have been measured using stopped-flow mixing methods and analyzed according to transition-state theory. GB contains no prolines, and the kinetics of folding and unfolding can be fit to a single, first-order rate constant over the temperature range of 5-35 degrees C. The major thermodynamic changes going from the unfolded state to the transition state are (1) a large decrease in heat capacity (delta Cp), indicating that the transition state is compact and solvent inaccessible relative to the unfolded state; (2) a large loss of entropy; and (3) a small increase in enthalpy. The most surprising feature of the folding of GB compared to that of previously studied proteins is that its folding approximates a rapid diffusion controlled process with little increase in enthalpy going from the unfolded to the transition state.     

Hydrophobic core variant ubiquitin forms a molten globule conformation at acidic pH

The conformational properties of hydrophobic core variant ubiquitin (Val26 to Ala mutation) in an acidic solution were studied. The intrinsic tryptophan fluorescence emission spectrum, far-UV and near-UV circular dichroic spectra, the fluorescence emission spectrum of 8-anilinonaphthalene-1-sulfonic acid in the presence of V26A ubiquitin, and urea-induced unfolding measurements indicate this variant ubiquitin to be in the partially folded molten globule conformation in solution at pH 2. The folding kinetics from molten globule to the native state was nearly identical to those from the unfolded state to the native state. This observation suggests that the equilibrium molten globule state of hydrophobic core variant ubiquitin is an on-pathway folding intermediate.     

The folding pathway of ubiquitin from all-atom molecular dynamics simulations

The folding (unfolding) pathway of ubiquitin is probed using all-atom molecular dynamics simulations. We dissect the folding pathway using two techniques: first, we probe the folding pathway of ubiquitin by calculating the evolution of structural properties over time and second, we identify the rate determining transition state for folding. The structural properties that we look at are hydrophobic solvent accessible surface area (SASA) and Calpha-root-mean-square deviation (rmsd). These properties on their own tell us relatively little about the folding pathway of ubiquitin; however, when plotted against each other, they become powerful tools for dissecting ubiquitin's folding mechanism. Plots of Calpha-rmsd against SASA serve as a phase space trajectories for the folding of ubiquitin. In this study, these plots show that ubiquitin folds to the native state via the population of an intermediate state. This is shown by an initial hydrophobic collapse phase followed by a second phase of secondary structure arrangement. Analysis of the structure of the intermediate state shows that it is a collapsed species with very little secondary structure. In reconciling these observations with recent experimental data, the transition that we observe in our simulations from the unfolded state (U) to the intermediate state (I) most likely occurs in the dead-time of the stopped flow instrument. The folding pathway of ubiquitin is probed further by identification of the rate-determining transition state for folding. The method used for this is essential dynamics, which utilizes a principal component analysis (PCA) on the atomic fluctuations throughout the simulation. The five transition state structures identified in silico are in good agreement with the experimentally determined transition state. The calculation of phi-values from the structures generated in the simulations is also carried out and it shows a good correlation with the experimentally measured values. An initial analysis of the denatured state shows that it is compact with fluctuating regions of nonnative secondary structure. It is found that the compactness in the denatured state is due to the burial of some hydrophobic residues. We conclude by looking at a correlation between folding kinetics and residual structure in the denatured state. A hierarchical folding mechanism is then proposed for ubiquitin.     

Effects of soman inhibition and of structural differences on cholinesterase molecular dynamics: a neutron scattering study

Incoherent elastic neutron scattering experiments on members of the cholinesterase family were carried out to investigate how molecular dynamics is affected by covalent inhibitor binding and by differences in primary and quaternary structure. Tetrameric native and soman-inhibited human butyrylcholinesterase (HuBChE) as well as native dimeric Drosophila melanogaster acetylcholinesterase (DmAChE) hydrated protein powders were examined. Atomic mean-square displacements (MSDs) were found to be identical for native HuBChE and for DmAChE in the whole temperature range examined, leading to the conclusion that differences in activity and substrate specificity are not reflected by a global modification of subnanosecond molecular dynamics. MSDs of native and soman-inhibited HuBChE were identical below the thermal denaturation temperature of the native enzyme, indicating a common mean free-energy surface. Denaturation of the native enzyme is reflected by a relative increase of MSDs consistent with entropic stabilization of the unfolded state. The results suggest that the stabilization of HuBChE phosphorylated by soman is due to an increase in free energy of the unfolded state due to a decrease in entropy.