The folding landscape of Streptomyces griseus protease B reveals the energetic costs and benefits associated with evolving kinetic stability

Like most extracellular bacterial proteases, Streptomyces griseus protease B (SGPB) and alpha-lytic protease (alphaLP) are synthesized with covalently attached pro regions necessary for their folding. In this article, we characterize the folding free energy landscape of SGPB and compare it to the folding landscapes of alphaLP and trypsin, a mammalian homolog that folds independently of its zymogen peptide. In contrast to the thermodynamically stable native state of trypsin, SGPB and alphaLP fold to native states that are thermodynamically marginally stable or unstable, respectively. Instead, their apparent stability arises kinetically, from unfolding free energy barriers that are both large and highly cooperative. The unique unfolding transitions of SGPB and alphaLP extend their functional lifetimes under highly degradatory conditions beyond that seen for trypsin; however, the penalty for evolving kinetic stability is remarkably large in that each factor of 2.4-8 in protease resistance is accompanied by a cost of ~10(5) in the spontaneous folding rate and ~5-9 kcal/mole in thermodynamic stability. These penalties have been overcome by the coevolution of increasingly effective pro regions to facilitate folding. Despite these costs, kinetic stability appears to be a potent mechanism for developing native-state properties that maximize protease longevity.     

Interdependent folding of the N- and C-terminal domains defines the cooperative folding of alpha-lytic protease

Alpha-lytic protease (alphaLP) serves as an important model in achieving a quantitative and physical understanding of protein folding reactions. Synthesized as a pro-protease, alphaLP belongs to an interesting class of proteins that require pro regions to facilitate their proper folding. alphaLP's pro region (Pro) acts as a potent folding catalyst for the protease, accelerating alphaLP folding to its native conformation nearly 10(10)-fold. Structural and mutational studies suggested that Pro's considerable foldase activity is directed toward structuring the alphaLP C-terminal domain (CalphaLP), a seemingly folding-impaired domain, which is believed to contribute significantly to the high-energy folding and unfolding transition states of alphaLP. Pro-mediated nucleation of alphaLP folding within CalphaLP was hypothesized to subsequently enable the alphaLP N-terminal domain (NalphaLP) to dock and fold, completing the formation of native protease. In this paper, we find that ternary folding reactions of Pro and noncovalent NalphaLP and CalphaLP domains are unaffected by the order in which the components are added or by the relative concentrations of the alphaLP domains, indicating that neither discrete CalphaLP structuring nor docking of the two alphaLP domains is involved in the folding transition state. Instead, the rate-limiting step of these folding reactions appears to be a slow and concerted rearrangement of the NalphaLP and CalphaLP domains to form active protease. This cooperative and interdependent folding of both protease domains defines the large alphaLP folding barrier and is an apparent extension of the highly cooperative alphaLP unfolding transition that imparts the protease with remarkable kinetic stability and functional longevity.     

alpha-lytic protease precursor: characterization of a structured folding intermediate

The bacterial alpha-lytic protease (alphaLP) is synthesized as a precursor containing a large N-terminal pro region (Pro) transiently required for correct folding of the protease [Silen, J. L., and Agard, D. A. (1989) Nature 341, 462-464]. Upon folding, the precursor is autocatalyticly cleaved to yield a tight-binding inhibitory complex of the pro region and the fully folded protease (Pro/alphaLP). An in vitro purification and refolding protocol has been developed for production of the disulfide-bonded precursor. A combination of spectroscopic approaches have been used to compare the structure and stability of the precursor with either the Pro/alphaLP complex or isolated Pro. The precursor and complex have significant similarities in secondary structure but some differences in tertiary structure, as well as a dramatic difference in stability. Correlations with isolated Pro suggest that the pro region part of the precursor is fully folded and acts to stabilize and structure the alphaLP region. Precursor folding is shown to be biphasic with the fast phase matching the rate of pro region folding. Further, the rate-limiting step in oxidative folding is formation of the disulfide bonds and autocatalytic processing occurs rapidly thereafter. These studies suggests a model in which the pro region folds first and catalyzes folding of the protease domain, forming the active site and finally causing autocatalytic cleavage of the bond separating pro region and protease. This last processing step is critical as it allows the protease N-terminus to rearrange, providing the majority of net stabilization of the product Pro/alphaLP complex.     

The pro region N-terminal domain provides specific interactions required for catalysis of alpha-lytic protease folding

The extracellular bacterial protease, alpha-lytic protease (alphaLP), is synthesized with a large, two-domain pro region (Pro) that catalyzes the folding of the protease to its native conformation. In the absence of its Pro folding catalyst, alphaLP encounters a very large folding barrier (DeltaG = 30 kcal mol(-1)) that effectively prevents the protease from folding (t(1/2) of folding = 1800 years). Although homology data, mutational studies, and structural analysis of the Pro.alphaLP complex suggested that the Pro C-terminal domain (Pro C-domain) serves as the minimum "foldase" unit responsible for folding catalysis, we find that the Pro N-terminal domain (Pro N-domain) is absolutely required for alphaLP folding. Detailed kinetic analysis of Pro N-domain point mutants and a complete N-domain deletion reveal that the Pro N-domain both provides direct interactions with alphaLP that stabilize the folding transition state and confers stability to the Pro C-domain. The Pro N- and C-domains make conflicting demands upon native alphaLP binding that are alleviated in the optimized interface of the folding transition state complex. From these studies, it appears that the extremely high alphaLP folding barrier necessitates the presence of both the Pro domains; however, alphaLP homologues with less demanding folding barriers may not require both domains, thus possibly explaining the wide variation in the pro region size of related pro-proteases.     

The folding landscape of an alpha-lytic protease variant reveals the role of a conserved beta-hairpin in the development of kinetic stability

Most secreted bacterial proteases, including alpha-lytic protease (alphaLP), are synthesized with covalently attached pro regions necessary for their folding. The alphaLP folding landscape revealed that its pro region, a potent folding catalyst, is required to circumvent an extremely large folding free energy of activation that appears to be a consequence of its unique unfolding transition. Remarkably, the alphaLP native state is thermodynamically unstable; a large unfolding free energy barrier is solely responsible for the persistence of its native state. Although alphaLP folding is well characterized, the structural origins of its remarkable folding mechanism remain unclear. A conserved beta-hairpin in the C-terminal domain was identified as a structural element whose formation and positioning may contribute to the large folding free energy barrier. In this article, we characterize the folding of an alphaLP variant with a more favorable beta-hairpin turn conformation (alphaLP(beta-turn)). Indeed, alphaLP(beta-turn) pro region-catalyzed folding is faster than that for alphaLP. However, instead of accelerating spontaneous folding, alphaLP(beta-turn) actually unfolds more slowly than alphaLP. Our data support a model where the beta-hairpin is formed early, but its packing with a loop in the N-terminal domain happens late in the folding reaction. This tight packing at the domain interface enhances the kinetic stability of alphaLP(beta-turn), to nearly the same degree as the change between alphaLP and a faster folding homolog. However, alphaLP(beta-turn) has impaired proteolytic activity that negates the beneficial folding properties of this variant. This study demonstrates the evolutionary limitations imposed by the simultaneous optimization of folding and functional properties.     

Two energetically disparate folding pathways of alpha-lytic protease share a single transition state

The Lysobacter enzymogenes alpha-lytic protease (alphaLP) is synthesized with a 166 amino acid pro region (Pro) that catalyzes the folding of the 198 amino acid protease into its native conformation. An extraordinary feature of this system is the very high energy barrier (DeltaG = 30 kcal mol-1) that effectively prevents alphaLP from folding in the absence of Pro (t1/2 = 1800 years). A pair of mutations has been isolated in the protease that completely suppresses the catalytic defect incurred in Pro by truncation of its last three amino acids. These mutations also accelerate the folding of alphaLP in the absence of Pro by 400-fold. An energetic analysis of the two folding reactions indicates that the mutations stabilize the transition states of both the catalyzed and uncatalyzed folding reactions by 3 kcal mol-1. This finding points to a single transition state for these two distinct and energetically disparate folding pathways, and raises the possibility that all alphaLP folding pathways share the same transition state.     

Mesophile versus thermophile: insights into the structural mechanisms of kinetic stability

Obtaining detailed knowledge of folding intermediate and transition state (TS) structures is critical for understanding protein folding mechanisms. Comparisons between proteins adapted to survive extreme temperatures with their mesophilic homologs are likely to provide valuable information on the interactions relevant to the unfolding transition. For kinetically stable proteins such as alpha-lytic protease (alphaLP) and its family members, their large free energy barrier to unfolding is central to their biological function. To gain new insights into the mechanisms that underlie kinetic stability, we have determined the structure and high temperature unfolding kinetics of a thermophilic homolog, Thermobifida fusca protease A (TFPA). These studies led to the identification of a specific structural element bridging the N and C-terminal domains of the protease (the "domain bridge") proposed to be associated with the enhanced high temperature kinetic stability in TFPA. Mutagenesis experiments exchanging the TFPA domain bridge into alphaLP validate this hypothesis and illustrate key structural details that contribute to TFPA's increased kinetic thermostability. These results lead to an updated model for the unfolding transition state structure for this important class of proteases in which domain bridge undocking and unfolding occurs at or before the TS. The domain bridge appears to be a structural element that can modulate the degree of kinetic stability of the different members of this class of proteases.     

Determinants of Affinity and Proteolytic Stability in Interactions of Kunitz Family Protease Inhibitors with Mesotrypsin

An important functional property of protein protease inhibitors is their stability to proteolysis. Mesotrypsin is a human trypsin that has been implicated in the proteolytic inactivation of several protein protease inhibitors. We have found that bovine pancreatic trypsin inhibitor (BPTI), a Kunitz protease inhibitor, inhibits mesotrypsin very weakly and is slowly proteolyzed, whereas, despite close sequence and structural homology, the Kunitz protease inhibitor domain of the amyloid precursor protein (APPI) binds to mesotrypsin 100 times more tightly and is cleaved 300 times more rapidly. To define features responsible for these differences, we have assessed the binding and cleavage by mesotrypsin of APPI and BPTI reciprocally mutated at two nonidentical residues that make direct contact with the enzyme. We find that Arg at P₁ (versus Lys) favors both tighter binding and more rapid cleavage, whereas Met (versus Arg) at P′₂ favors tighter binding but has minimal effect on cleavage. Surprisingly, we find that the APPI scaffold greatly enhances proteolytic cleavage rates, independently of the binding loop. We draw thermodynamic additivity cycles analyzing the interdependence of P₁ and P′₂ substitutions and scaffold differences, finding multiple instances in which the contributions of these features are nonadditive. We also report the crystal structure of the mesotrypsin·APPI complex, in which we find that the binding loop of APPI displays evidence of increased mobility compared with BPTI. Our data suggest that the enhanced vulnerability of APPI to mesotrypsin cleavage may derive from sequence differences in the scaffold that propagate increased flexibility and mobility to the binding loop.     

Effect of signal peptide on the stability and folding kinetics of maltose binding protein

While the role of the signal sequence in targeting proteins to specific subcellular compartments is well characterized, there are fewer studies that characterize its effects on the stability and folding kinetics of the protein. We report a detailed characterization of the folding kinetics and thermodynamic stabilities of maltose binding protein (MBP) and its precursor form, preMBP. Isothermal GdmCl and urea denaturation as a function of temperature and thermal denaturation studies have been carried out to compare stabilities of the two proteins. preMBP was found to be destabilized by about 2-6 kcal/mol (20-40%) with respect to MBP. Rapid cleavage of the signal peptide by various proteases shows that the signal peptide is accessible in the native form of preMBP. The observed rate constant of the major slow phase in folding was decreased 5-fold in preMBP relative to MBP. The rate constants of unfolding were similar at 25 degrees C, but preMBP also exhibited a large burst phase change in unfolding that was absent in MBP. At 10 degrees C, preMBP exhibited a higher unfolding rate than MBP as well as a large burst phase. The appreciable destabilization of MBP by signal peptide is functionally relevant, because it enhances the likelihood of finding the protein in an unfolded translocation-competent form and may influence the interactions of the protein with the translocation machinery. Destabilization is likely to result from favorable interactions between the hydrophobic signal peptide and other hydrophobic regions that are exposed in the unfolded state.     

Towards a consistent modeling of protein thermodynamic and kinetic cooperativity: how applicable is the transition state picture to folding and unfolding?

To what extent do general features of folding/unfolding kinetics of small globular proteins follow from their thermodynamic properties? To address this question, we investigate a new simplified protein chain model that embodies a cooperative interplay between local conformational preferences and hydrophobic burial. The present four-helix-bundle 55mer model exhibits protein-like calorimetric two-state cooperativity. It rationalizes native-state hydrogen exchange observations. Our analysis indicates that a coherent, self-consistent physical account of both the thermodynamic and kinetic properties of the model leads naturally to the concept of a native state ensemble that encompasses considerable conformational fluctuations. Such a multiple-conformation native state is seen to involve conformational states similar to those revealed by native-state hydrogen exchange. Many of these conformational states are predicted to lie below native baselines commonly used in interpreting calorimetric data. Folding and unfolding kinetics are studied under a range of intrachain interaction strengths as in experimental chevron plots. Kinetically determined transition midpoints match well with their thermodynamic counterparts. Kinetic relaxations are found to be essentially single-exponential over an extended range of model interaction strengths. This includes the entire unfolding regime and a significant part of a folding regime with a chevron rollover, as has been observed for real proteins that fold with non-two-state kinetics. The transition state picture of protein folding and unfolding is evaluated by comparing thermodynamic free energy profiles with actual kinetic rates. These analyses suggest that some chevron rollovers may arise from an internal frictional effect that increasingly impedes chain motions with more native conditions, rather than being caused by discrete deadtime folding intermediates or shifts of the transition state peak as previously posited.     

Oxidative modification of Escherichia coli glutamine synthetase. Decreases in the thermodynamic stability of protein structure and specific changes in the active site conformation.

Metal catalyzed oxidation of specific amino acid residues has been proposed to be an important physiological mechanism of marking proteins for proteolytic degradation. After initial oxidative inactivation of dodecameric Escherichia coli glutamine synthetase (GS), the integrity of the GS active site and protein structure was assessed by monitoring ATP binding, observing a susceptibility of GS to tryptic cleavage, and comparative thermodynamic analysis. The tryptic cleavage rates of an active site linked central loop were significantly accelerated for the oxidized conformer. This tryptic cleavage was essentially prevented in the presence of glutamate for native GS but not for the oxidized conformer. The integrity of the ATP binding site in the oxidized GS was substantially altered as indicated by the reduction in fluorescence enhancement associated with ATP binding. Decreases in the free energies of quaternary protein structure and subunit interactions due to oxidative modification were determined by temperature and urea induced unfolding equilibrium measurements. Comparative thermal stability measurements of a partial unfolding transition indicated that the loss in stabilization free energy for the oxidized GS conformer was 1.3 kcal/mol dodecamer. Under alkaline conditions, the urea-induced disruption of quaternary and tertiary structures of oxidized and native GS were examined. This comparative analysis revealed that the free energies of the subunit interactions and unfolding of the dissociated monomers for oxidized GS were decreased by 1.5 and 1.7 kcal/mol, respectively. Our results suggest that small free energy decreases in GS protein structural stability of only 1-2 kcal/mol may be responsible for the selective proteolytic turnover of the oxidized GS.     

On the unfolding of alpha-lytic protease and the role of the pro region

Molecular dynamics simulations of alpha-lytic protease (alphaLP) alone and complexed with its pro region (PRO) are performed to understand the origin of its high unfolding (and folding) barrier when it is alone and how the pro region lowers this barrier. At room temperature, alphaLP exhibits lower dynamic fluctuations than alpha-chymotrypsin. Simulation of PRO alone led to reorientation of its N terminal helix and collapse to a more compact state. A model for the uncleaved proenzyme was built and found to be stable in the time scale of the simulations. Energetic analysis suggests that the origin of strain in the uncleaved proenzyme compared with the cleaved complex is in the intramolecular backbone electrostatic interactions of the cleaved strand. In high temperature simulations, the interaction of the long beta hairpin of the enzyme with the C terminal beta sheet of PRO is among the most stable in the complex and a likely "nucleation site" for folding. In the course of unfolding, the C terminal tail of PRO is sometimes observed to intervene between the long hairpin and the aspartate loop of the enzyme, perhaps thereby lowering the energy barrier for separation of the two hairpins. Tighter interactions at the interface between the enzyme and its pro region are also occasionally observed, providing an additional mechanism for unfolding catalysis. Simulations of a mutant enzyme where the buried ion pair residues R102 and D142 were replaced by W and L, respectively, did not display any distinguishable behavior compared with the wild type.     

Kinetic barriers to the folding of horse cytochrome C in the reduced state

To determine the kinetic barrier in the folding of horse cytochrome c, a CO-liganded derivative of cytochrome c, called carbonmonoxycytochrome c, has been prepared by exploiting the thermodynamic reversibility of ferrocytochrome c unfolding induced by guanidinium hydrochloride (GdnHCl), pH 7. The CO binding properties of unfolded ferrocytochrome c, studied by 13C NMR and optical spectroscopy, are remarkably similar to those of native myoglobin and isolated chains of human hemoglobin. Equilibrium unfolding transitions of ferrocytochrome c in the presence and the absence of CO observed by both excitation energy transfer from the lone tryptophan to the ferrous heme and far-UV circular dichroism (CD) indicate no accumulation of structural intermediates to a detectable level. Values of thermodynamic parameters obtained by two-state analysis of fluorescence transitions are DeltaG(H2O) = 11.65(+/-1.13) kcal x mol(-1) and C(m) = 3.9(+/-0.1) M GdnHCl in the presence of CO, and DeltaG(H2O)=19.3(+/-0.5) kcal x mol(-1) and C(m) = 5.1(+/-0.1) M GdnHCl in the absence of CO, indicating destabilization of ferrocytochrome c by approximately 7.65 kcal x mol(-1) due to CO binding. The native states of ferrocytochrome c and carbonmonoxycytochrome c are nearly identical in terms of structure and conformation except for the Fe2+-M80 --> Fe2+-CO replacement. Folding and unfolding kinetics as a function of GdnHCl, studied by stopped-flow fluorescence, are significantly different for the two proteins. Both refold fast, but carbonmonoxycytochrome c refolds 2-fold faster (tau = 1092 micros at 10 degrees C) than ferrocytochrome c. Linear extrapolation of the folding rates to the ordinate of the chevron plot projects this value of tau to 407 micros. The unfolding rate of the former in water, estimated by extrapolation, is faster by more than 10 orders of magnitude. Significant differences are also observed in rate-denaturant gradients in the chevron. Formation and disruption of the Fe2+-M80 coordination contact clearly impose high-energy kinetic barriers to folding and unfolding of ferrocytochrome c. The unfolding barrier due to the Fe2+-M80 bond provides sufficient kinetic stability to the native state of ferrocytochrome c to perform its physiological function as an electron donor.     

Weak cooperativity in the core causes a switch in folding mechanism between two proteins of the cks family

The human protein ckshs1 (cks1) is a 79 residue alpha/beta protein with low thermodynamic and kinetic stability. Its folding mechanism was probed by mutation at sites throughout the structure. Many of the mutations caused changes in the slope of the unfolding arm of the chevron plot. The effects can be rationalised in terms of either transition-state movement or native-state "breathing", and in either case, the magnitude of the effect enables the sequence of events in the folding reaction to be determined. Those sites that fold early exhibit a small perturbation, whilst those sites that fold late exhibit a large perturbation. The results show that cks1 folds sequential pairs of beta-strands first; beta1/beta2 and beta3/beta4. Subsequently, these pairs pack against each other and onto the alpha-helical region to form the core. The folding process of cks1 contrasts with that of the homologue, suc1. The 113 residue suc1 has the same beta-sheet core structure but, additionally, two large insertions that confer much greater thermodynamic and kinetic stability. The more extensive network of tertiary interactions in suc1 provides sufficient enthalpic gain to overcome the entropic cost of forming the core and thus tips the balance in favour of non-local interactions: the non-local, central beta-strand pair, beta2/beta4, forms first and the periphery strands pack on later. Moreover, the greater cooperativity of the core of suc1 protects its folding from perturbation and consequently the slope of the unfolding arm of the chevron plot is much less sensitive to mutation.     

Partially unfolded forms and non-two-state folding of a beta-sandwich: FHA domain from Arabidopsis receptor kinase-associated protein phosphatase

FHA domains adopt a beta-sandwich fold with 11 strands. The first evidence of partially unfolded forms of a beta-sandwich is derived from native-state hydrogen exchange (NHX) of the forkhead-associated (FHA) domain from kinase-associated protein phosphatase from Arabidopsis. The folding kinetics of this FHA domain indicate that EX2 behavior prevails at pH 6.3. In the chevron plot, rollover in the folding arm and bends in the unfolding arm suggest folding intermediates. NHX of this FHA domain suggests a core of six most stable beta-strands and two loops, characterized by rare global unfolding events. Flanking this stable core are beta-strands and recognition loops with less stability, termed subglobal motifs. These suggest partially unfolded forms (near-native intermediates) with two levels of stability. The spatial separation of the subglobal motifs on the flanks suggests possible parallelism in their folding as additional beta-strands align with the stable core of six strands. Intermediates may contribute to differences in stabilities and m-values suggested by NHX or kinetics relative to chemical denaturation. Residual structure in the unfolded regime is suggested by superprotection of beta-strand 6 and by GdmCl-dependence of adjustments in amide NMR spectra and residual optical signal. The global folding stability depends strongly on pH, with at least 3 kcal/mol more stability at pH 7.3 than at pH 6.3. This FHA domain is hypothesized to fold progressively with initial hydrophobic collapse of its stable six-stranded core followed by addition of less stable flanking beta-strands and ordering of recognition loops.     

Folding pathway mediated by an intramolecular chaperone: dissecting conformational changes coincident with autoprocessing and the role of Ca(2+) in subtilisin maturation.

Subtilisin is produced as a precursor that requires its N-terminal propeptide to chaperone the folding of its protease domain. Once folded, subtilisin adopts a remarkably stable conformation, which has been attributed to a high affinity Ca(2+) binding site. We investigated the role of the metal ligand in the maturation of pro-subtilisin, a process that involves folding, autoprocessing and partial degradation. Our results establish that although Ca(2+) ions can stabilize the protease domain, the folding and autoprocessing of pro-subtilisin take place independent of Ca(2+) ion. We demonstrate that the stabilizing effect of calcium is observed only after the completion of autoprocessing and that the metal ion appears to be responsible for shifting the folding equilibrium towards the native conformation in both mature subtilisin and the autoprocessed propeptide:subtilisin complex. Furthermore, the addition of active subtilisin to unautoprocessed pro-subtilisin in trans does not facilitate precursor maturation, but rather promotes rapid autodegradation. The primary cleavage site that initiates this autodegradation is at Gln19 in the N-terminus of mature subtilisin. This corresponds to the loop that links alpha-helix-2 and beta-strand-1 in mature subtilisin and has indirect effects on the formation of the Ca(2+) binding site. Our results show that the N-terminus of mature subtilisin undergoes rearrangement subsequent to propeptide autoprocessing. Since this structural change enhances the proteolytic stability of the precursor, our results suggest that the autoprocessing reaction must be completed before the release of active subtilisin in order to maximize folding efficiency.     

Rapid folding and unfolding of Apaf-1 CARD

Caspase recruitment domains (CARDs) are members of the death domain superfamily and contain six antiparallel helices in an alpha-helical Greek key topology. We have examined the equilibrium and kinetic folding of the CARD of Apaf-1 (apoptotic protease activating factor 1), which consists of 97 amino acid residues, at pH 6 and pH 8. The results showed that an apparent two state equilibrium mechanism is not adequate to describe the folding of Apaf-1 CARD at either pH, suggesting the presence of intermediates in equilibrium unfolding. Interestingly, the results showed that the secondary structure is less stable than the tertiary structure, based on the transition mid-points for unfolding. Single mixing and sequential mixing stopped-flow studies showed that Apaf-1 CARD folds and unfolds rapidly and suggest a folding mechanism that contains parallel channels with two unfolded conformations folding to the native conformation. Kinetic simulations show that a slow folding phase is described by a third conformation in the unfolded ensemble that interconverts with one or both unfolded species. Overall, the native ensemble is formed rapidly upon refolding. This is in contrast to other CARDs in which folding appears to be dominated by formation of kinetic traps.     

The Folding Free-Energy Surface of HIV-1 Protease: Insights into the Thermodynamic Basis for Resistance to Inhibitors

Spontaneous mutations at numerous sites distant from the active site of human immunodeficiency virus type 1 protease enable resistance to inhibitors while retaining enzymatic activity. As a benchmark for probing the effects of these mutations on the conformational adaptability of this dimeric β-barrel protein, the folding free-energy surface of a pseudo-wild-type variant, HIV-PR^?, was determined by a combination of equilibrium and kinetic experiments on the urea-induced unfolding/refolding reactions. The equilibrium unfolding reaction was well described by a two-state model involving only the native dimeric form and the unfolded monomer. The global analysis of the kinetic folding mechanism reveals the presence of a fully folded monomeric intermediate that associates to form the native dimeric structure. Independent analysis of a stable monomeric version of the protease demonstrated that a small-amplitude fluorescence phase in refolding and unfolding, not included in the global analysis of the dimeric protein, reflects the presence of a transient intermediate in the monomer folding reaction. The partially folded and fully folded monomers are only marginally stable with respect to the unfolded state, and the dimerization reaction provides a modest driving force at micromolar concentrations of protein. The thermodynamic properties of this system are such that mutations can readily shift the equilibrium from the dimeric native state towards weakly folded states that have a lower affinity for inhibitors but that could be induced to bind to their target proteolytic sites. Presumably, subsequent secondary mutations increase the stability of the native dimeric state in these variants and, thereby, optimize the catalytic properties of the resistant human immunodeficiency virus type 1 protease.     

Understanding the folding and stability of a designed WW domain protein with replica exchange molecular dynamics simulations

WW domain proteins are usually regarded as simple models for understanding the folding mechanism of β-sheet. CC45 is an artificial protein that is capable of folding into the same structure as WW domain. In this article, the replica exchange molecular dynamics simulations are performed to investigate the folding mechanism of CC45. The analysis of thermal stability shows that β-hairpin 1 is more stable than β-hairpin 2 during the unfolding process. Free energy analysis shows that the unfolding of this protein substantially proceeds through solvating the smaller β-hairpin 2, followed by the unfolding of β-hairpin 1. We further propose the unfolding process of CC45 and the folding mechanism of two β-hairpins. These results are similar to the previous folding studies of formin binding protein 28 (FBP28). Compared with FBP28, it is found that CC45 has more aromatic residues in N-terminal loop, and these residues contact with C-terminal loop to form the outer hydrophobic core, which increases the stability of CC45. Knowledge about the stability and folding behaviour of CC45 may help in understanding the folding mechanisms of the β-sheet and in designing new WW domains.     

Loop and stem dynamics during RNA hairpin folding and unfolding

2-Aminopurine (2AP) is a fluorescent adenine analog that probes mainly base stacking in nucleic acids. We labeled the loop or the stem of the RNA hairpin gacUACGguc with 2AP to study folding thermodynamics and kinetics at both loci. Thermal melts and fast laser temperature jumps detected by 2AP fluorescence monitored the stability and folding/unfolding kinetics. The observed thermodynamic and kinetic traces of the stem and loop mutants, though strikingly different at a first glance, can be fitted to the same free-energy landscape. The differences between the two probe locations arise because base stacking decreases upon unfolding in the stem, whereas it increases in the loop. We conclude that 2AP is a conservative adenine substitution for mapping out the contributions of different RNA structural elements to the overall folding process. Molecular dynamics (MD) totaling 0.6 μsec were performed to look at the conformations populated by the RNA at different temperatures. The combined experimental data, and MD simulations lead us to propose a minimal four-state free-energy landscape for the RNA hairpin. Analysis of this landscape shows that a sequential folding model is a good approximation for the full folding dynamics. The frayed state formed initially from the native state is a heterogeneous ensemble of structures whose stem is frayed either from the end or from the loop.