Cooperative folding of a protein mini domain: the peripheral subunit-binding domain of the pyruvate dehydrogenase multienzyme complex

From covalent transition states in chemistry to noncovalent in biology: from β- to Φ-value analysis of protein folding

Solving the mechanism of a chemical reaction requires determining the structures of all the ground states on the pathway and the elusive transition states linking them. 2024 is the centenary of Brønsted's landmark paper that introduced the β-value and structure-activity studies as the only experimental means to infer the structures of transition states. It involves making systematic small changes in the covalent structure of the reactants and analysing changes in activation and equilibrium-free energies. Protein engineering was introduced for an analogous procedure, Φ-value analysis, to analyse the noncovalent interactions in proteins central to biological chemistry. The methodology was developed first by analysing noncovalent interactions in transition states in enzyme catalysis. The mature procedure was then applied to study transition states in the pathway of protein folding - 'part (b) of the protein folding problem'. This review describes the development of Φ-value analysis of transition states and compares and contrasts the interpretation of β- and Φ-values and their limitations. Φ-analysis afforded the first description of transition states in protein folding at the level of individual residues. It revealed the nucleation-condensation folding mechanism of protein domains with the transition state as an expanded, distorted native structure, containing little fully formed secondary structure but many weak tertiary interactions. A spectrum of transition states with various degrees of structural polarisation was then uncovered that spanned from nucleation-condensation to the framework mechanism of fully formed secondary structure. Φ-analysis revealed how movement of the expanded transition state on an energy landscape accommodates the transition from framework to nucleation-condensation mechanisms with a malleability of structure as a unifying feature of folding mechanisms. Such movement follows the rubric of analysis of classical covalent chemical mechanisms that began with Brønsted. Φ-values are used to benchmark computer simulation, and Φ and simulation combine to describe folding pathways at atomic resolution.

Unfolding Transitions of Peripheral Subunit Binding Domains Show Cooperative Behavior

Characterization of native, intermediate, and denatured states is crucial for understanding the factors influencing the stability of proteins. We have carried out molecular dynamics simulations to study the unfolding of three peripheral subunit binding domains (PSBDs): E. coli BBL, Bacillus stearothermophilus E3BD, and human hbSBD, at three different temperatures: 300, 330, and 400 K, and in the presence of two solvents: water and 5 M guanidinium hydrochloride (GndCl) solution. These proteins share similar folds, with two parallel helices, maintained via a hydrophobic core comprising residues from their interconnecting loop. BBL is more sensitive to thermal and chemical denaturation in comparison to hbSBD, and E3BD is the most stable of all of the three proteins. The effect of temperature on the stability of these proteins is more pronounced in "water-only" simulations compared to that in the presence of guanidium hydrochloride in high concentrations. Our results show cooperative unfolding transitions of these proteins, which are triggered by an initial melting of the C-terminal helix H2. The consequent loss of interhelical interactions or native contacts, as observed, leads to the subsequent melting of the N-terminal helix H1.

Isotope-Labeled Aspartate Sidechain as a Non-Perturbing Infrared Probe: Application to Investigate the Dynamics of a Carboxylate Buried Inside a Protein

Because of their negatively charged carboxylates, aspartate and glutamate are frequently found at the active or binding site of proteins. However, studying a specific carboxylate in proteins that contain multiple aspartates and/or glutamates via infrared spectroscopy is difficult due to spectral overlap. We show, herein, that isotopic-labeling of the aspartate sidechain can overcome this limitation as the resultant ¹³C=O asymmetric stretching vibration resides in a transparent region of the protein IR spectrum. Applicability of this site-specific vibrational probe is demonstrated by using it to assess the dynamics of an aspartate ion buried inside a small protein via two-dimensional infrared spectroscopy.

Experimental and Computational Analysis of Protein Stabilization by Gly-to-d-Ala Substitution: A Convolution of Native State and Unfolded State Effects

The rational and predictable enhancement of protein stability is an important goal in protein design. Most efforts target the folded state, however stability is the free energy difference between the folded and unfolded states thus both are suitable targets. Strategies directed at the unfolded state usually seek to decrease chain entropy by introducing cross-links or by replacing glycines. Cross-linking has led to mixed results. Replacement of glycine with an l-amino acid, while reducing the entropy of the unfolded state, can introduce unfavorable steric interactions in the folded state, since glycine is often found in conformations that require a positive φ angle such as helical C-capping motifs or type I' and II″ β-turns. l-Amino acids are strongly disfavored in these conformations, but d-amino acids are not. However, there are few reported examples and conflicting results have been obtained when glycines are replaced with d-Ala. We critically examine the effect of Gly-to-d-Ala substitutions on protein stability using experimental approaches together with molecular dynamics simulations and free energy calculations. The data, together with a survey of high resolution structures, show that the vast majority of proteins can be stabilized by substitution of C-capping glycines with d-Ala. Sites suitable for substitutions can be identified via sequence alignment with a high degree of success. Steric clashes in the native state due to the new side chain are rarely observed, but are likely responsible for the destabilizing or null effect observed for the small subset of Gly-to-d-Ala substitutions which are not stabilizing. Changes in backbone solvation play less of a role. Favorable candidates for d-Ala substitution can be identified using a rapid algorithm based on molecular mechanics.

Estimation of protein folding free energy barriers from calorimetric data by multi-model Bayesian analysis

The realization that folding free energy barriers can be small enough to result in significant population of the species at the barrier top has sprouted in several methods to estimate folding barriers from equilibrium experiments. Some of these approaches are based on fitting the experimental thermogram measured by differential scanning calorimetry (DSC) to a one-dimensional representation of the folding free-energy surface (FES). Different physical models have been used to represent the FES: (1) a Landau quartic polynomial as a function of the total enthalpy, which acts as an order parameter; (2) the projection onto a structural order parameter (i.e. number of native residues or native contacts) of the free energy of all the conformations generated by Ising-like statistical mechanical models; and (3) mean-field models that define conformational entropy and stabilization energy as functions of a continuous local order parameter. The fundamental question that emerges is how can we obtain robust, model-independent estimates of the thermodynamic folding barrier from the analysis of DSC experiments. Here we address this issue by comparing the performance of various FES models in interpreting the thermogram of a protein with a marginal folding barrier. We chose the small α-helical protein PDD, which folds-unfolds in microseconds crossing a free energy barrier previously estimated as ~1 RT. The fits of the PDD thermogram to the various models and assumptions produce FES with a consistently small free energy barrier separating the folded and unfolded ensembles. However, the fits vary in quality as well as in the estimated barrier. Applying Bayesian probabilistic analysis we rank the fit performance using a statistically rigorous criterion that leads to a global estimate of the folding barrier and its precision, which for PDD is 1.3 ± 0.4 kJ mol(-1). This result confirms that PDD folds over a minor barrier consistent with the downhill folding regime. We have further validated the multi-model Bayesian approach through the analysis of two additional protein systems: gpW, a midsize single-domain with α + β topology that also folds in microseconds and has been previously catalogued as a downhill folder, and α-spectrin SH3, a domain of similar size but with a β-barrel fold, slow-folding kinetics and two-state-like thermodynamics. From a general viewpoint, the Bayesian analysis developed here results in a statistically robust, virtually model-independent, method to estimate the thermodynamic free-energy barriers to protein folding from DSC thermograms. Our method appears to be sufficiently accurate to consistently detect small differences in the barrier height, and thus opens up the possibility of characterizing experimentally the changes in thermodynamic folding barriers induced by single-point mutations on proteins within the downhill regime.

The human peripheral subunit-binding domain folds rapidly while overcoming repulsive Coulomb forces

Peripheral subunit binding domains (PSBDs) are integral parts of large multienzyme complexes involved in carbohydrate metabolism. PSBDs facilitate shuttling of prosthetic groups between different catalytic subunits. Their protein surface is characterized by a high density of positive charges required for binding to subunits within the complex. Here, we investigated folding thermodynamics and kinetics of the human PSBD (HSBD) using circular dichroism and tryptophan fluorescence experiments. HSBD was only marginally stable under physiological solvent conditions but folded within microseconds via a barrier-limited apparent two-state transition, analogous to its bacterial homologues. The high positive surface-charge density of HSBD leads to repulsive Coulomb forces that modulate protein stability and folding kinetics, and appear to even induce native-state movement. The electrostatic strain was alleviated at high solution-ionic-strength by Debye-Hückel screening. Differences in ionic-strength dependent characteristics among PSBD homologues could be explained by differences in their surface charge distributions. The findings highlight the trade-off between protein function and stability during protein evolution.

Navigating the downhill protein folding regime via structural homologues

Proteins that fold over free-energy barriers <or= 3RT are classified as downhill folders. This regime is characterized by equilibrium unfolding that is proportionally broader and more complex the lower the folding barrier. Downhill proteins are also expected to fold up in a few microseconds. However, the relationship between rate and equilibrium signatures is affected by other factors such as protein size and folding topology. Here we perform a direct comparison of the kinetics and equilibrium unfolding of two structural homologues: BBL and PDD. BBL folds-unfolds in just approximately 1 micros at 335 K and displays the equilibrium signatures expected for a protein at the bottom of the downhill folding regime. PDD, which has the same 3D structure and size, folds-unfolds approximately 8 times more slowly and, concomitantly, exhibits all the downhill equilibrium signatures to a lesser degree. Our results demonstrate that the equilibrium signatures of downhill folding are proportional to the changes in folding rate once structural and size-scaling effects are factored out. This conclusion has two important implications: (1) it confirms that the quantitative analysis of equilibrium experiments in ultrafast folding proteins does provide direct information about free-energy barriers, a result that is incompatible with the conventional view of protein folding as a highly activated process, and (2) it advocates for equilibrium-kinetic studies of homologous proteins as a powerful tool to navigate the downhill folding regime via comparative analysis. The latter should prove extremely useful for the investigation of sequence, functional, and evolutionary determinants of protein folding barriers.

Direct observation of barrier-limited folding of BBL by single-molecule fluorescence resonance energy transfer

One controversial area in protein folding mechanisms is whether some small, ultra-fast-folding proteins exist in distinct native and denatured state ensembles, separated by an energy barrier, or if there is a continuum of states between native and denatured. In theory, the simplest way of distinguishing between single-state barrierless or "downhill" folding and conventional separate state folding is by single-molecule spectroscopy, which can detect either distinct populations of proteins or a continuum. But, the time resolution of approximately 1 ms of most confocal fluorescence microscopes for single-molecule fluorescence resonance energy transfer (SM-FRET) is longer than that for the structural relaxation of proteins such as BBL, whose mechanism of folding is controversial. We have constructed a highly sensitive confocal fluorescence microscope and measured the distribution of FRET efficiencies of appropriately labeled BBL in time bins of 50 and 200 mus under conditions in which its structural relaxation time is 340 mus or less. The experiments are at the very limits of detection because of signal artefacts from shot noise, photo-bleaching, and other events that broaden signals of individual states so they appear to coalesce. However, with appropriate tuning of the thresholds for detection and length of data collection, we clearly observed 2 distinct states of BBL, with FRET efficiencies corresponding to native and denatured states. The population of each state varied with GdmCl or urea during chemical denaturation transitions corresponding to conventional barrier-limited folding at 279 K and pH 7 and pH 5.8. The folding of BBL is accordingly barrier limited.

Observation of noncooperative folding thermodynamics in simulations of 1BBL

One of the predictions of the energy landscape theory of protein folding is the possibility of barrierless, "downhill" folding under certain conditions. The protein 1BBL has been proposed to fold by such a downhill mechanism, though this is a matter of some dispute. We carried out extensive replica exchange molecular dynamics simulations on 1BBL in explicit solvent to address this controversy and provide a microscopic picture of its folding thermodynamics. Our simulations show two distinct structural transitions in the folding of 1BBL. A low-temperature transition involves a disordering of the protein's tertiary structure without loss of secondary structure. A distinct, higher temperature transition involves the complete loss of secondary structure and dissolution of the hydrophobic core. In contrast, control simulations of the 1BBL homolog E3BD show a single high temperature unfolding transition. Further simulations of 1BBL at high ionic strength show a significant destabilization of helix II but not helix I, suggesting that the apparent folding cooperativity of 1BBL may be highly dependent on experimental conditions. Although our simulations cannot provide definitive evidence of downhill folding in 1BBL, they clearly show evidence of a complex, non-two-state folding process.

Structural bases for the specific interactions between the E2 and E3 components of the Thermus thermophilus 2-oxo acid dehydrogenase complexes

Pyruvate dehydrogenase (PDH), branched-chain 2-oxo acid dehydrogenase (BCDH) and 2-oxoglutarate dehydrogenase (OGDH) are multienzyme complexes that play crucial roles in several common metabolic pathways. These enzymes belong to a family of 2-oxo acid dehydrogenase complexes that contain multiple copies of three different components (E1, E2 and E3). For the Thermus thermophilus enzymes, depending on its substrate specificity (pyruvate, branched-chain 2-oxo acid or 2-oxoglutarate), each complex has distinctive E1 (E1p, E1b or E1o) and E2 (E2p, E2b or E2o) components and one of the two possible E3 components (E3b and E3o). (The suffixes, p, b and o identify their respective enzymes, PDH, BCDH and OGDH.) Our biochemical characterization demonstrates that only three specific E3*E2 complexes can form (E3b*E2p, E3b*E2b and E3o*E2o). X-ray analyses of complexes formed between the E3 components and the peripheral subunit-binding domains (PSBDs), derived from the corresponding E2-binding partners, reveal that E3b interacts with E2p and E2b in essentially the same manner as observed for Geobacillus stearothermophilus E3*E2p, whereas E3o interacts with E2o in a novel fashion. The buried intermolecular surfaces of the E3b*PSBDp/b and E3o*PSBDo complexes differ in size, shape and charge distribution and thus, these differences presumably confer the binding specificities for the complexes.

Rational design, structural and thermodynamic characterization of a hyperstable variant of the villin headpiece helical subdomain

A hyperstable variant of the small independently folded helical subdomain (HP36) derived from the F-actin binding villin headpiece was designed by targeting surface electrostatic interactions and helical propensity. A double mutant N68A, K70M was significantly more stable than wild type. The Tm of wild type in aqueous buffer is 73.0 degrees C, whereas the double mutant did not display a complete unfolding transition. The double mutant could not be completely unfolded even by 10 M urea. In 3 M urea, the Tm of wild type is 54.8 degrees C while that of the N68AK70M double mutant is 73.9 degrees C. Amide H/2H exchange studies show that the pattern of exchange is very similar for wild type and the double mutant. The structures of a K70M single mutant and the double mutant were determined by X-ray crystallography and are identical to that of the wild type. Analytical ultracentrifugation demonstrates that the proteins are monomeric. The hyperstable mutant described here is expected to be useful for folding studies of HP36 because studies of the wild type domain have sometimes been limited by its marginal stability. The results provide direct evidence that naturally occurring miniature protein domains have not been evolutionarily optimized for global stability. The stabilizing effect of this double mutant could not be predicted by sequence analysis because K70 is conserved in the larger intact headpiece for functional reasons.

The transition state for folding of a peripheral subunit-binding domain contains robust and ionic-strength dependent characteristics

The denaturant dependencies of the folding and unfolding kinetics were used to characterize the structure of the transition state for folding of E3BD, a peripheral subunit-binding domain. For the majority of E3BD mutants, the Phi-values calculated at 298 K from the analysis of chevron plots were in good agreement with those previously determined at 325 K using Arrhenius analysis. This agreement further demonstrates the general robustness of Phi-value analyses, since different experiments, methods of denaturation and thermodynamic assumptions were used to determine each set of Phi(F) values. The structure of the transition state for folding was grossly conserved at 298 K and 325 K, with residues in Helix I playing a lesser role in folding than those located in the 3(10) helix, disordered loop and Helix II. However, the energetic contributions of a cluster of basic residues close to the N-terminus and Helix I, which are an integral part of the ligand-binding site, were susceptible to ionic strength effects because of electrostatic strain in native and transition states of E3BD at low ionic strength. We found no evidence of the downhill folding previously proposed for E3BD, even though the conditions employed in this study significantly increased the energetic bias towards the native state.

Ultra-fast Barrier-limited Folding in the Peripheral Subunit-binding Domain Family

We have determined the solution structures, equilibrium properties and ultra-fast folding kinetics for three bacterial homologues of the peripheral subunit-binding domain (PSBD) family. The mesophilic homologue, BBL, was less stable than the thermophilic and hyper-thermophilic variants (E3BD and POB, respectively). The broad unfolding transitions of each PSBD, when probed by different techniques, were essentially superimposable, consistent with co-operative denaturation. Temperature-jump and continuous-flow fluorescence methods were used to measure the folding kinetics for E3BD, POB and BBL. E3BD folded fairly rapidly at 298K (folding half-time approximately 25 micros) and BBL and POB folded even faster (folding half-times approximately 3-5 micros). The variations in equilibrium and kinetic behaviour observed for the PSBD family resembles that of the homeodomain family, where the folding pattern changes from apparent two-state transitions to multi-state kinetics as the denatured state becomes more structured. The faster folding of POB may be a consequence of its higher propensity to form helical structure in the region corresponding to the folding nucleus of E3BD. The ultra-fast folding of BBL appears to be a consequence of residual structure in the denatured ensemble, as with engrailed homeodomain. We discuss issues concerning "one-state", downhill folding, and find no evidence for, and strong evidence against, it occurring in these PSBDs. The shorter construct used previously for BBL was destabilized significantly and the stability further perturbed by the introduction of fluorescent probes. Thermal titrations for 11 side-chains scattered around the protein, when probed by (13)C-NMR experiments, could be fit globally to a common co-operative transition.

Simulation and experiment at high temperatures: ultrafast folding of a thermophilic protein by nucleation-condensation

We used Phi-value analysis to characterise the transition state for folding of a thermophilic protein at the relatively high temperature of 325 K. PhiF values for the folding of the three-helix bundle, peripheral subunit binding domain from Bacillus stearothermophilus (E3BD) were determined by temperature-jump experiments in the absence of chemical denaturants. E3BD folded in microseconds through a highly diffuse transition state. Excellent agreement was observed between experiment and the results from eight (independent) molecular dynamics simulations of unfolding at 373 K. We used a combination of heteronuclear NMR experiments and molecular dynamics simulations to characterise the denatured ensemble, and found that it contained very little persistent, residual structure. However, those regions that adopt helical structure in the native state were found by simulation to be poised for helix formation in the denatured state. These regions also had significant structure in the transition state for folding. The overall folding pathway appears to be nucleation-condensation.

One-state Downhill versus Conventional Protein Folding

Classical protein folding invokes a cooperative transition between distinct thermodynamic states that are individually populated at equilibrium and separated by an energy barrier. It has been proposed, however, that the small protein, BBL, undergoes one-step downhill folding whereby it folds non-cooperatively to its native state without encountering an appreciable energy barrier. Only a single conformational ensemble is populated under given conditions, and so the denatured state ensemble progressively changes into the native structure. A wide dispersion of thermal denaturation midpoints that was observed for an extrinsically labelled fragment of BBL is proposed to be evidence for its one-state, downhill folding, a phenomenon that is also suggested to be functionally important for BBL and its homologues. We found, however, that thermal denaturation of unlabelled wild-type BBL was highly cooperative, with very similar transition midpoints for the melting of secondary and tertiary interactions, as well as for individual residues when monitored by NMR. Similar results were also observed for two other homologues, E3BD and POB. Further, the extrinsic fluorophores perturbed the unfolding energetics of labelled BBL, and complicated its equilibrium behaviour. One-step downhill folding may well occur for some proteins that do not have distinct folded states but not for BBL and its well-folded homologues.

Thermal stability and DNA binding activity of a variant form of the Sso7d protein from the archeon Sulfolobus solfataricus truncated at leucine 54

Sso7d is a 62-residue, basic protein from the hyperthermophilic archaeon Sulfolobus solfataricus. Around neutral pH, it exhibits a denaturation temperature close to 100 degrees C and a non-sequence-specific DNA binding activity. Here, we report the characterization by circular dichroism and fluorescence measurements of a variant form of Sso7d truncated at leucine 54 (L54Delta). It is shown that L54Delta has a folded conformation at neutral pH and that its thermal unfolding is a reversible process, represented well by the two-state N <=> D transition model, with a denaturation temperature of 53 degrees C. Fluorescence titration experiments indicate that L54Delta binds tightly to calf thymus DNA, even though the binding parameters are smaller than those of the wild-type protein. Therefore, the truncation of eight residues at the C-terminus of Sso7d markedly affects the thermal stability of the protein, which nevertheless retains a folded structure and DNA binding activity.

Rapid cooperative two-state folding of a miniature alpha-beta protein and design of a thermostable variant

There is a great deal of interest in developing small stably folded miniature proteins. A limited number of these molecules have been described, however they typically have not been characterized in depth. In particular, almost no detailed studies of the thermodynamics and folding kinetics of these proteins have been reported. Here we describe detailed studies of the thermodynamics and kinetics of folding of a 39 residue mixed alpha-beta protein (NTL9(1-39)) derived from the N-terminal domain of the ribosomal protein L9. The protein folds cooperatively and rapidly in a two-state fashion to a native state typical of those found for normal globular proteins. At pH 5.4 in 20mM sodium acetate, 100mM NaCl the temperature of maximum stability is 6 degrees C, the t(m) is 65.3 degrees C, deltaH degrees (t(m)) is between 24.6 kcalmol(-1) and 26.3 kcalmol(-1), and deltaC(p) degrees is 0.38 kcalmol(-1)deg(-1). The thermodynamic parameters are in the range expected on the basis of per residue values determined from databases of globular proteins. H/2H exchange measurements reveal a set of amides that exchange via global unfolding, exactly as expected for a normal cooperatively folded globular protein. Kinetic measurements show that folding is two-state folding. The folding rate is 640 s(-1) and the value of deltaG degrees calculated from the folding and unfolding rates is in excellent agreement with the equilibrium value. A designed thermostable variant, generated by mutating K12 to M, was characterized and found to have a t(m) of 82 degrees C. Equilibrium and kinetic measurements demonstrate that its folding is cooperative and two-state.

Experimental identification of downhill protein folding

Theory predicts the existence of barrierless protein folding. Without barriers, folding should be noncooperative and the degree of native structure should be coupled to overall protein stability. We investigated the thermal unfolding of the peripheral subunit binding domain from Escherichia coli's 2-oxoglutarate dehydrogenase multienzyme complex (termed BBL) with a combination of spectroscopic techniques and calorimetry. Each technique probed a different feature of protein structure. BBL has a defined three-dimensional structure at low temperatures. However, each technique showed a distinct unfolding transition. Global analysis with a statistical mechanical model identified BBL as a downhill-folding protein. Because of BBL's biological function, we propose that downhill folders may be molecular rheostats, in which effects could be modulated by altering the distribution of an ensemble of structures.

Characterization of large peptide fragments derived from the N-terminal domain of the ribosomal protein L9: definition of the minimum folding motif and characterization of local electrostatic interactions

A set of peptides derived from the N-terminal domain of the ribosomal protein L9 (NTL9) have been characterized in an effort to define the minimum unit of this domain required to fold and to provide model peptides for the analysis of electrostatic interactions in the unfolded state. NTL9 is a 56-residue alpha-beta protein with a beta1-loop-beta2-alpha1-beta3-alpha2 topology. The beta-sheet together with the first helix comprise a simple example of a common supersecondary motif called the split beta-alpha-beta fold. Peptides corresponding to the beta1-loop-beta2 unit are unstructured even when constrained by an introduced disulfide. The pK(a)s of Asp-8 and Glu-17 in these peptides are slightly lower than the values found for shorter peptides but are considerably higher than the values in NTL9. A 34-residue peptide, which represents the beta1-loop-beta2-alpha1 portion of NTL9, is also unstructured. In contrast, a 39-residue peptide corresponding to the entire split beta-alpha-beta motif is folded and monomeric as judged by near- and far-UV CD, two-dimensional NMR, ANS binding experiments, pK(a) measurements, and analytical ultracentrifugation. The fold is very similar to the structure of this region in the intact protein. Thermal and urea unfolding experiments show that it is cooperatively folded with a DeltaG degrees of unfolding of 1.8-2.0 kcal/mol and a T(m) of 58 degrees C. This peptide represents the first demonstration of the independent folding of an isolated split beta-alpha-beta motif, and is one of only four naturally occurring sequences of fewer than 40 residues that has been shown to fold cooperatively in the absence of disulfides or ligand binding.

NMR Structures of 36 and 73-residue Fragments of the Calreticulin P-domain

Calreticulin (CRT) is an abundant, soluble molecular chaperone of the endoplasmic reticulum. Similar to its membrane-bound homolog calnexin (CNX), it is a lectin that promotes the folding of proteins carrying N-linked glycans. Both proteins cooperate with an associated co-chaperone, the thiol-disulfide oxidoreductase ERp57. This enzyme catalyzes the formation of disulfide bonds in CNX and CRT-bound glycoprotein substrates. Previously, we solved the NMR structure of the central proline-rich P-domain of CRT comprising residues 189-288. This structure shows an extended hairpin topology, with three short anti-parallel beta-sheets, three small hydrophobic clusters, and one helical turn at the tip of the hairpin. We further demonstrated that the residues 225-251 at the tip of the CRT P-domain are involved in direct contacts with ERp57. Here, we show that the CRT P-domain fragment CRT(221-256) constitutes an autonomous folding unit, and has a structure highly similar to that of the corresponding region in CRT(189-288). Of the 36 residues present in CRT(221-256), 32 form a well-structured core, making this fragment one of the smallest known natural sequences to form a stable non-helical fold in the absence of disulfide bonds or tightly bound metal ions. CRT(221-256) comprises all the residues of the intact P-domain that were shown to interact with ERp57. Isothermal titration microcalorimetry (ITC) now showed affinity of this fragment for ERp57 similar to that of the intact P-domain, demonstrating that CRT(221-256) may be used as a low molecular mass mimic of CRT for further investigations of the interaction with ERp57. We also solved the NMR structure of the 73-residue fragment CRT(189-261), in which the tip of the hairpin and the first beta-sheet are well structured, but the residues 189-213 are disordered, presumably due to lack of stabilizing interactions across the hairpin.

Identification of key amino acid residues in the assembly of enzymes into the pyruvate dehydrogenase complex of Bacillus stearothermophilus: a kinetic and thermodynamic analysis

Structural studies have shown that electrostatic interactions play a major part in the binding of dihydrolipoyl dehydrogenase (E3) to the peripheral subunit-binding domain (PSBD) of the dihydrolipoyl acyltransferase (E2) in the assembly of the pyruvate dehydrogenase multienzyme complex of Bacillus stearothermophilus. The binding is characterized by a small, unfavorable enthalpy change (deltaH degrees = +2.2 kcal/mol) and a large, positive entropy change (TdeltaS degrees = +14.8 kcal/mol). The contributions of individual surface residues of the PSBD of E2 to its interaction with E3 have been assessed by alanine-scanning mutagenesis, surface plasmon resonance detection, and isothermal titration calorimetry. The mutation R135A in the PSBD gave rise to a significant decrease (120-fold) in the binding affinity; two other mutations (R139A and R156A) were associated with smaller effects. The binding of the R135A mutant to E3 was accompanied by a favorable enthalpy (deltaH degrees = -2.6 kcal/mol) and a less positive entropy change (TdeltaS degrees = +7.2 kcal/mol). The midpoint melting temperature (T(m)) of E3-PSBD complexes was determined by differential scanning calorimetry. The R135A mutation caused a significant decrease (5 degrees C) in the T(m), compared with the wild-type complex. The results reveal the importance of Arg135 of the PSBD as a key residue in the molecular recognition of E3 by E2, and as a major participant in the overall entropy-driven binding process. Further, the effects of mutagenesis on the deltaCp of subunit association illustrate the difficulties in attributing changes in heat capacity to specific classes of interactions.

The role of aromatic residues in the hydrophobic core of the villin headpiece subdomain

Small autonomously folding proteins are of interest as model systems to study protein folding, as the same molecule can be used for both experimental and computational approaches. The question remains as to how well these minimized peptide model systems represent larger native proteins. For example, is the core of a minimized protein tolerant to mutation like larger proteins are? Also, do minimized proteins use special strategies for specifying and stabilizing their folded structure? Here we examine these questions in the 35-residue autonomously folding villin headpiece subdomain (VHP subdomain). Specifically, we focus on a cluster of three conserved phenylalanine (F) residues F47, F51, and F58, that form most of the hydrophobic core. These three residues are oriented such that they may provide stabilizing aromatic-aromatic interactions that could be critical for specifying the fold. Circular dichroism and 1D-NMR spectroscopy show that point mutations that individually replace any of these three residues with leucine were destabilized, but retained the native VHP subdomain fold. In pair-wise replacements, the double mutant that retains F58 can adopt the native fold, while the two double mutants that lack F58 cannot. The folding of the double mutant that retains F58 demonstrates that aromatic-aromatic interactions within the aromatic cluster are not essential for specifying the VHP subdomain fold. The ability of the VHP subdomain to tolerate mutations within its hydrophobic core indicates that the information specifying the three dimensional structure is distributed throughout the sequence, as observed in larger proteins. Thus, the VHP subdomain is a legitimate model for larger, native proteins.

Folding rate prediction using total contact distance

Linear regression analysis found that either contact order (CO) or long-range order (LRO) parameter has a significant correlation with the logarithms of folding rates. This suggests that sequence separation per contact and total number of contacts are both important in determining the rate of folding. Here, the two factors are incorporated into a new parameter, total contact distance (TCD). Using a database of 28 two-state or weakly three-state folding proteins, TCD is found to be the most accurate among the three parameters (CO, LRO, and TCD) in terms of correlation and prediction. It provides even more accurate prediction than the best neural network results with two descriptors (contact order and stability per residue). The improvement is achieved in all three-structural classes (all alpha, beta, and mixed). The accuracy of total contact distance in predicting folding rates is essentially unchanged if "short"-ranged contacts (absolute value of i - j < or = 14) are not included in calculation. Thus, only long-range contacts with a sequence separation of more than 14 residues are important in determining the rate of folding. This is consistent with the results from the long-range order parameter. One of the significant outliers in prediction is found to be associated with the only protein in the database that involves nonlocal disulfide bonds. Removing the protein leads to a correlation coefficient of 0.89 between experimental observed and predicted folding rates in jackknife cross validation. The corresponding values for CO and LRO are 0.71 and 0.80, respectively.

Folding intermediates of a model three-helix bundle protein. Pressure and cold denaturation studies

The stability and equilibrium unfolding of a model three-helix bundle protein, alpha(3)-1, by guanidine hydrochloride (GdnHCl), hydrostatic pressure, and temperature have been investigated. The combined use of these denaturing agents allowed detection of two partially folded states of alpha(3)-1, as monitored by circular dichroism, intrinsic fluorescence emission, and fluorescence of the hydrophobic probe bis-ANS (4,4'-dianilino-1,1'-binaphthyl-5,5'-disulfonic acid). The overall free-energy change for complete unfolding of alpha(3)-1, determined from GdnHCl unfolding data, is +4.6 kcal/mol. The native state is stabilized by -1.4 kcal/mol relative to a partially folded pressure-denatured intermediate (I(1)). Cold denaturation at high pressure gives rise to a second partially (un)folded conformation (I(2)), suggesting a significant contribution of hydrophobic interactions to the stability of alpha(3)-1. The free energy of stabilization of the native-like state relative to I(2) is evaluated to be -2.5 kcal/mol. Bis-ANS binding to the pressure- and cold-denatured states indicates the existence of significant residual hydrophobic structure in the partially (un)folded states of alpha(3)-1. The demonstration of folding intermediates of alpha(3)-1 lends experimental support to a number of recent protein folding simulation studies of other three-helix bundle proteins that predicted the existence of such intermediates. The results are discussed in terms of the significance of de novo designed proteins for protein folding studies.

Swinging arms and swinging domains in multifunctional enzymes: catalytic machines for multistep reactions

Multistep chemical reactions are increasingly seen as important in a growing number of complex biotransformations. Covalently attached prosthetic groups or swinging arms, and their associated protein domains, are essential to the mechanisms of active-site coupling and substrate channeling in a number of the multifunctional enzyme systems responsible. The protein domains, for which the posttranslational machinery in the cell is highly specific, are crucially important, contributing to the processes of molecular recognition that define and protect the substrates and the catalytic intermediates. The domains have novel folds and move by virtue of conformationally flexible linker regions that tether them to other components of their respective multienzyme complexes. Structural and mechanistic imperatives are becoming apparent as the assembly pathways and the coupling of multistep reactions catalyzed by these dauntingly complex molecular machines are unraveled.

Rational modification of protein stability by the mutation of charged surface residues

Continuum methods were used to calculate the electrostatic contributions of charged and polar side chains to the overall stability of a small 41-residue helical protein, the peripheral subunit-binding domain. The results of these calculations suggest several residues that are destabilizing, relative to hydrophobic isosteres. One position was chosen to test the results of these calculations. Arg8 is located on the surface of the protein in a region of positive electrostatic potential. The calculations suggest that Arg8 makes a significant, unfavorable electrostatic contribution to the overall stability. The experiments described in this paper represent the first direct experimental test of the theoretical methods, taking advantage of solid-phase peptide synthesis to incorporate approximately isosteric amino acid substitutions. Arg8 was replaced with norleucine (Nle), an amino acid that is hydrophobic and approximately isosteric, or with alpha-amino adipic acid (Aad), which is also approximately isosteric but oppositely charged. In this manner, it is possible to isolate electrostatic interactions from the effects of hydrophobic and van der Waals interactions. Both Arg8Nle and Arg8Aad are more thermostable than the wild-type sequence, testifying to the validity of the calculations. These replacements led to stability increases at 52.6 degrees C, the T(m) of the wild-type, of 0.86 and 1.08 kcal mol(-)(1), respectively. The stability of Arg8Nle is particularly interesting as a rare case in which replacement of a surface charge with a hydrophobic residue leads to an increase in the stability of the protein.

Submillisecond folding of the peripheral subunit-binding domain

Folding and unfolding rates have been measured for the peripheral subunit-binding domain, a small three-helix protein. The protein folds very fast, with rates too rapid to be measured using traditional stopped-flow techniques. Folding and unfolding rates were measured as a function of temperature using dynamic NMR lineshape analysis. At the lowest temperature at which there is sufficient broadening to measure rates, 41 degrees C, the folding rate is 16,050 s(-1). Thus, the halftime required for folding is 43 micros. At the same temperature, the unfolding rate is 2800 s(-1). Identical rates were measured using resolved resonances from Val16 in the loop and Val21 at the end of the 310-helix. Folding rates have been correlated with protein topology, and this correlation is consistent with the rapid folding of the peripheral subunit-binding domain. The results presented here show that the peripheral subunit-binding domain is the third fastest folding protein for which rates have been estimated. The folding rate is the fastest that has been directly measured and provides further support for the importance of chain topology as a major determinant of folding rates.

Folding propensities of synthetic peptide fragments covering the entire sequence of phage 434 Cro protein

The phage 434 Cro protein, the N-terminal domain of its repressor (R1-69) and that of phage lambda (lambda6-85) constitute a group of small, monomeric, single-domain folding units consisting of five helices with striking structural similarity. The intrinsic helix stabilities in lambda6-85 have been correlated to its rapid folding behavior, and a residual hydrophobic cluster found in R1-69 in 7 M urea has been proposed as a folding initiation site. To understand the early events in the folding of 434 Cro, and for comparison with R1-69 and lambda6-85, we examined the conformational behavior of five peptides covering the entire 434 Cro sequence in water, 40% (by volume) TFE/water, and 7 M urea solutions using CD and NMR. Each peptide corresponds to a helix and adjacent residues as identified in the native 434 Cro NMR and crystal structures. All are soluble and monomeric in the solution conditions examined except for the peptide corresponding to the 434 Cro helix 4, which has low water solubility. Helix formation is observed for the 434 Cro helix 1 and helix 2 peptides in water, for all the peptides in 40% TFE and for none in 7 M urea. NMR data indicate that the helix limits in the peptides are similar to those in the native protein helices. The number of side-chain NOEs in water and TFE correlates with the helix content, and essentially none are observed in 7 M urea for any peptide, except that for helix 5, where a hydrophobic cluster may be present. The low intrinsic folding propensities of the five helices could account for the observed stability and folding behavior of 434 Cro and is, at least qualitatively, in accord with the results of the recently described diffusion-collision model incorporating intrinsic helix propensities.

Nativelike structure and stability in a truncation mutant of a protein minidomain: the peripheral subunit-binding domain

Despite its small size, the peripheral subunit-binding domain from the dihydrolipoamide acetyltransferase component of the Bacillus stearothermophilus pyruvate dehydrogenase multienzyme complex adopts a unique, compact structure. To determine whether the full 43 residue sequence is required for the domain to adopt a stable, nativelike structure, 3 proteins of different lengths were prepared. Psbd41 corresponds to residues 3-43 of the domain, psbd36 spans residues 6-41, and psbd33 comprises residues 7-39. Psbd41 folds in a cooperative, two-state fashion with a Tm of 53 degrees C and a stability at 25 degrees C of 2.2 kcal mol-1. Psbd36 is nearly as stable with a Tm of 48 degrees C and a stability of 1.8 kcal mol-1. Similar m-values and heat capacities suggest that psbd36 and psbd41 bury approximately the same surface area. Minimal differences in CalphaH and NH chemical shifts between psbd41 and psbd36 show that the two sequences adopt the same tertiary fold. On a per residue basis, DeltaH degrees and DeltaC degrees p fall within the range typical for single-domain globular proteins. Psbd33 is significantly less stable. It is not fully folded at 25 degrees C, and at all temperatures it shows broadened NMR lines. ANS titrations provide evidence that this is due to an equilibrium between nativelike and unfolded molecules rather than formation of a molten globule. The fraction of psbd33 molecules which are folded appear to adopt the same structure as the full-length domain. Thus, although more than the 33 residue core is required to form a fully stable native structure, the entire sequence is not required for folding.

Conformational analysis of peptide fragments derived from the peripheral subunit-binding domain from the pyruvate dehydrogenase multienzyme complex of Bacillus stearothermophilus: evidence for nonrandom structure in the unfolded state

There is currently a great deal of interest in the early events in protein folding. Two issues that have generated particular interest are the nature of the unfolded state under native conditions and the role of local interactions in folding. Here, we report the results of a study of a set of peptides derived from a small two-helix protein, the peripheral subunit-binding domain of the pyruvate dehydrogenase multienzyme complex. Five peptides of overlapping sequence were prepared, including sequences corresponding to each of the helices and to the region connecting them. The peptides were characterized by CD and, where possible, nmr. A peptide corresponding to the second helix is between 12 and 17% helical at neutral pH. CD also indicates a lower percentage of helical structure in the peptide corresponding to the first alpha-helix, although the values of the alpha-proton chemical shifts suggest some preference for nonrandom structure. Peptides corresponding to the interhelical loop, which in the full domain contains two overlapping beta-turns and a 5-residue 3(10)-helix, are less structured. There is no significant change in the helicity of any of these peptides with pH. To test for fragment complementation, CD spectra of the two peptides derived from each helix and the long connecting peptide were compared to the spectra of each possible pair, as well as to a mixture containing all three. No increase in structure was observed. We complement our peptide studies by characterizing a point mutant, D34V, which disrupts a critical hydrogen bonding network. This mutant is unable to fold and provides a useful model of the denatured state. The mutant is between 9 and 16% helical as judged by CD. The modest amount of helical structure formed in some of the peptide fragments and in the point mutant suggests that the denatured state of the peripheral subunit binding domain is not completely unstructured. This may contribute to the very rapid folding observed for the intact protein.

Global analysis of the effects of temperature and denaturant on the folding and unfolding kinetics of the N-terminal domain of the protein L9

The folding and unfolding kinetics of the N-terminal domain of the ribosomal protein L9 have been measured at temperatures between 7 and 85 degrees C and between 0 and 6 M guanidine deuterium chloride. Stopped-flow fluorescence was used to measure rates below 55 degrees C and NMR lineshape analysis was used above 55 degrees C. The amplitudes and rate profiles of the stopped-flow fluorescence experiments are consistent with a two-state folding mechanism, and plots of ln(k) versus guanidine deuterium chloride concentration show the classic v-shape indicative of two-state folding. There is no roll over in the plots when the experiments are repeated in the presence of 400 mM sodium sulfate. Temperature and denaturant effects were fit simultaneously to the simple model k=D exp(-DeltaG*/RT) where DeltaG* represents the change in apparent free energy between the transition state and the folded or unfolded state and D represents the maximum possible folding speed. DeltaG* is assumed to vary linearly with denaturant concentration and the Gibbs-Helmholtz equation is used to model stability changes with temperature. Approximately 60% of the surface area buried upon folding is buried in the transition state as evidenced by changes in the heat capacity and m value between the unfolded state and the transition state. The equilibrium thermodynamic parameters, DeltaCp degrees, m and DeltaG degrees, all agree with the values calculated from the kinetic experiments, providing additional evidence that folding is two-state. The folding rates at 0 M guanidine hydrochloride show a non-Arrhenius temperature dependence typical of globular proteins. When the folding rates are examined along constant DeltaG degrees/T contours they display an Arrhenius temperature dependence with a slope of -8600 K. This indicates that for this system, the non-Arrhenius temperature dependence of folding can be accounted for by the anomalous temperature dependence of the interactions which stabilize proteins.