Peripheral Subunit-binding Domain Research Articles

Dimerization of a 5-kDa domain defines the architecture of the 5-MDa gammaproteobacterial pyruvate dehydrogenase complex.

The Escherichia coli pyruvate dehydrogenase complex (PDHc) is a ~5 MDa assembly of the catalytic subunits pyruvate dehydrogenase (E1), dihydrolipoamide acetyltransferase (E2), and dihydrolipoamide dehydrogenase (E3). The PDHc core is a cubic complex of eight E2 homotrimers. Homodimers of the peripheral subunits E1 and E3 associate with the core by binding to the peripheral subunit binding domain (PSBD) of E2. Previous reports indicated that 12 E1 dimers and 6 E3 dimers bind to the 24-meric E2 core. Using an assembly arrested E2 homotrimer (E23), we show that two of the three PSBDs in the E23 dimerize, that each PSBD dimer cooperatively binds two E1 dimers, and that E3 dimers only bind to the unpaired PSBD in E23. This mechanism is preserved in wild-type PDHc, with an E1 dimer:E2 monomer:E3 dimer stoichiometry of 16:24:8. The conserved PSBD dimer interface indicates that PSBD dimerization is the previously unrecognized architectural determinant of gammaproteobacterial PDHc megacomplexes.

Studying micro to millisecond protein dynamics using simple amide 15N CEST experiments supplemented with major-state R2 and visible peak-position constraints.

Over the last decade amide 15N CEST experiments have emerged as a popular tool to study protein dynamics that involves exchange between a 'visible' major state and sparsely populated 'invisible' minor states. Although initially introduced to study exchange between states that are in slow exchange with each other (typical exchange rates of, 10 to 400s-1), they are now used to study interconversion between states on the intermediate to fast exchange timescale while still using low to moderate (5 to 350Hz) 'saturating' B1 fields. The 15N CEST experiment is very sensitive to exchange as the exchange delay TEX can be quite long (~0.5s) allowing for a large number of exchange events to occur making it a very powerful tool to detect minor sates populated ([Formula: see text]) to as low as 1%. When systems are in fast exchange and the 15N CEST data has to be described using a model that contains exchange, the exchange parameters are often poorly defined because the [Formula: see text] versus [Formula: see text] and [Formula: see text] versus exchange rate ([Formula: see text]) plots can be quite flat with shallow or no minima and the analysis of such 15N CEST data can lead to wrong estimates of the exchange parameters due to the presence of 'spurious' minima. Here we show that the inclusion of experimentally derived constraints on the intrinsic transverse relaxation rates and the inclusion of visible state peak-positions during the analysis of amide 15N CEST data acquired with moderate B1 values (~50 to ~350Hz) results in convincing minima in the [Formula: see text] versus [Formula: see text] and the [Formula: see text] versus [Formula: see text] plots even when exchange occurs on the 100μs timescale. The utility of this strategy is demonstrated on the fast-folding Bacillus stearothermophilus peripheral subunit binding domain that folds with a rate constant ~104s-1. Here the analysis of 15N CEST data alone results in [Formula: see text] versus [Formula: see text] and [Formula: see text] versus [Formula: see text] plots that contain shallow minima, but the inclusion of visible-state peak positions and restraints on the intrinsic transverse relaxation rates of both states during the analysis of the 15N CEST data results in pronounced minima in the [Formula: see text] versus [Formula: see text] and [Formula: see text] versus [Formula: see text] plots and precise exchange parameters even in the fast exchange regime ([Formula: see text]~5). Using this strategy we find that the folding rate constant of PSBD is invariant (~10,500s-1) from 33.2 to 42.9°C while the unfolding rates (~70 to ~500s-1) and unfolded state populations (~0.7 to ~4.3%) increase with temperature. The results presented here show that protein dynamics occurring on the 10 to 104s-1 timescale can be studied using amide 15N CEST experiments.

MRPS36 provides a structural link in the eukaryotic 2-oxoglutarate dehydrogenase complex.

The tricarboxylic acid cycle is the central pathway of energy production in eukaryotic cells and plays a key part in aerobic respiration throughout all kingdoms of life. One of the pivotal enzymes in this cycle is 2-oxoglutarate dehydrogenase complex (OGDHC), which generates NADH by oxidative decarboxylation of 2-oxoglutarate to succinyl-CoA. OGDHC is a megadalton protein complex originally thought to be assembled from three catalytically active subunits (E1o, E2o, E3). In fungi and animals, however, the protein MRPS36 has more recently been proposed as a putative additional component. Based on extensive cross-linking mass spectrometry data supported by phylogenetic analyses, we provide evidence that MRPS36 is an important member of the eukaryotic OGDHC, with no prokaryotic orthologues. Comparative sequence analysis and computational structure predictions reveal that, in contrast with bacteria and archaea, eukaryotic E2o does not contain the peripheral subunit-binding domain (PSBD), for which we propose that MRPS36 evolved as an E3 adaptor protein, functionally replacing the PSBD. We further provide a refined structural model of the complete eukaryotic OGDHC of approximately 3.45 MDa with novel mechanistic insights.

Structure of the native pyruvate dehydrogenase complex reveals the mechanism of substrate insertion

The pyruvate dehydrogenase complex (PDHc) links glycolysis to the citric acid cycle by converting pyruvate into acetyl-coenzyme A. PDHc encompasses three enzymatically active subunits, namely pyruvate dehydrogenase, dihydrolipoyl transacetylase, and dihydrolipoyl dehydrogenase. Dihydrolipoyl transacetylase is a multidomain protein comprising a varying number of lipoyl domains, a peripheral subunit-binding domain, and a catalytic domain. It forms the structural core of the complex, provides binding sites for the other enzymes, and shuffles reaction intermediates between the active sites through covalently bound lipoyl domains. The molecular mechanism by which this shuttling occurs has remained elusive. Here, we report a cryo-EM reconstruction of the native E. coli dihydrolipoyl transacetylase core in a resting state. This structure provides molecular details of the assembly of the core and reveals how the lipoyl domains interact with the core at the active site.

Characterization of lysine acetylation in the peripheral subunit-binding domain of the E2 subunit of the pyruvate dehydrogenase-2-oxoglutarate dehydrogenase hybrid complex from Corynebacterium glutamicum.

In Corynebacterium glutamicum, pyruvate dehydrogenase (PDH) and 2-oxoglutarate dehydrogenase (ODH) form a unique hybrid complex in which CgE1p and CgE1o are associated with the CgE2-CgE3 subcomplex. We analyzed the role of a lysine acetylation site in the peripheral subunit-binding domain of CgE2 in PDH and ODH functions. Acetylation-mimic substitution at Lys391 of CgE2 severely reduced the interaction of CgE2 with CgE1p and CgE3, but not with CgE1o, indicating the critical role of this residue in the assembly of CgE1p and CgE3 into the complex. It also suggested that Lys391 acetylation inhibited the binding of CgE1p and CgE3 to CgE2, thereby affecting PDH and ODH activities. Interestingly, the CgE2-K391R variant strain showed increased l-glutamate production and reduced pyruvate accumulation. Kinetic analysis suggested that the increased affinity of the K391R variant toward pyruvate might be advantageous for l-glutamate production.

The N-terminal 1-55 residues domain of pyruvate dehydrogenase from Escherichia coli assembles as a dimer in solution.

The pyruvate dehydrogenase complex (PDHc) from Escherichia coli is a large protein complex consisting of multiple copies of the pyruvate dehydrogenase (E1ec), dihydrolipoamide acetyltransferase (E2ec) and dihydrolipoamide dehydrogenase (E3ec). The N-terminal domain (NTD, residues 1-55) of E1ec plays a critical role in the interaction between E1ec and E2ec and the whole PDHc activity. Using circular dichroism, size-exclusion chromatography and dynamic light scattering spectroscopy, we show that the NTD of E1ec presents dimeric assembly under physiological condition. Pull-down and isothermal titration calorimetry binding assays revealed that the E2ec peripheral subunit-binding domain (PSBD) forms a very stable complex with the NTD, indicating the isolated NTD functionally interacts with PSBD and the truncated E1ec (E1ec∆NTD) does not interact with PSBD. These findings are important to understand the mechanism of PDHc and other thiamine-based multi-component enzymes.

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.

Measuring the signs of the methyl 1H chemical shift differences between major and 'invisible' minor protein conformational states using methyl 1H multi-quantum spectroscopy.

Carr-Purcell-Meiboom-Gill (CPMG) type relaxation dispersion experiments are now routinely used to characterise protein conformational dynamics that occurs on the μs to millisecond (ms) timescale between a visible major state and 'invisible' minor states. The exchange rate(s) ([Formula: see text]), population(s) of the minor state(s) and the absolute value of the chemical shift difference [Formula: see text] (ppm) between different exchanging states can be extracted from the CPMG data. However the sign of [Formula: see text] that is required to reconstruct the spectrum of the 'invisible' minor state(s) cannot be obtained from CPMG data alone. Building upon the recently developed triple quantum (TQ) methyl [Formula: see text] CPMG experiment (Yuwen in Angew Chem 55:11490-11494, 2016) we have developed pulse sequences that use carbon detection to generate and evolve single quantum (SQ), double quantum (DQ) and TQ coherences from methyl protons in the indirect dimension to measure the chemical exchange-induced shifts of the SQ, DQ and TQ coherences from which the sign of [Formula: see text] is readily obtained for two state exchange. Further a combined analysis of the CPMG data and the difference in exchange induced shifts between the SQ and DQ resonances and between the SQ and TQ resonances improves the estimates of exchange parameters like the population of the minor state. We demonstrate the use of these experiments on two proteins undergoing exchange: (1) the ~18kDa cavity mutant of T4 Lysozyme ([Formula: see text]) and (2) the [Formula: see text]kDa Peripheral Sub-unit Binding Domain (PSBD) from the acetyl transferase of Bacillus stearothermophilus ([Formula: see text]).

Novel Binding Motif and New Flexibility Revealed by Structural Analyses of a Pyruvate Dehydrogenase-Dihydrolipoyl Acetyltransferase Subcomplex from the Escherichia coli Pyruvate Dehydrogenase Multienzyme Complex

The Escherichia coli pyruvate dehydrogenase multienzyme complex contains multiple copies of three enzymatic components, E1p, E2p, and E3, that sequentially carry out distinct steps in the overall reaction converting pyruvate to acetyl-CoA. Efficient functioning requires the enzymatic components to assemble into a large complex, the integrity of which is maintained by tethering of the displaced, peripheral E1p and E3 components to the E2p core through non-covalent binding. We here report the crystal structure of a subcomplex between E1p and an E2p didomain containing a hybrid lipoyl domain along with the peripheral subunit-binding domain responsible for tethering to the core. In the structure, a region at the N terminus of each subunit in the E1p homodimer previously unseen due to crystallographic disorder was observed, revealing a new folding motif involved in E1p-E2p didomain interactions, and an additional, unexpected, flexibility was discovered in the E1p-E2p didomain subcomplex, both of which probably have consequences in the overall multienzyme complex assembly. This represents the first structure of an E1p-E2p didomain subcomplex involving a homodimeric E1p, and the results may be applicable to a large range of complexes with homodimeric E1 components. Results of HD exchange mass spectrometric experiments using the intact, wild type 3-lipoyl E2p and E1p are consistent with the crystallographic data obtained from the E1p-E2p didomain subcomplex as well as with other biochemical and NMR data reported from our groups, confirming that our findings are applicable to the entire E1p-E2p assembly.

Structure and Function of the Catalytic Domain of the Dihydrolipoyl Acetyltransferase Component in Escherichia coli Pyruvate Dehydrogenase Complex

The Escherichia coli pyruvate dehydrogenase complex (PDHc) catalyzing conversion of pyruvate to acetyl-CoA comprises three components: E1p, E2p, and E3. The E2p is the five-domain core component, consisting of three tandem lipoyl domains (LDs), a peripheral subunit binding domain (PSBD), and a catalytic domain (E2pCD). Herein are reported the following. 1) The x-ray structure of E2pCD revealed both intra- and intertrimer interactions, similar to those reported for other E2pCDs. 2) Reconstitution of recombinant LD and E2pCD with E1p and E3p into PDHc could maintain at least 6.4% activity (NADH production), confirming the functional competence of the E2pCD and active center coupling among E1p, LD, E2pCD, and E3 even in the absence of PSBD and of a covalent link between domains within E2p. 3) Direct acetyl transfer between LD and coenzyme A catalyzed by E2pCD was observed with a rate constant of 199 s(-1), comparable with the rate of NADH production in the PDHc reaction. Hence, neither reductive acetylation of E2p nor acetyl transfer within E2p is rate-limiting. 4) An unprecedented finding is that although no interaction could be detected between E1p and E2pCD by itself, a domain-induced interaction was identified on E1p active centers upon assembly with E2p and C-terminally truncated E2p proteins by hydrogen/deuterium exchange mass spectrometry. The inclusion of each additional domain of E2p strengthened the interaction with E1p, and the interaction was strongest with intact E2p. E2p domain-induced changes at the E1p active site were also manifested by the appearance of a circular dichroism band characteristic of the canonical 4'-aminopyrimidine tautomer of bound thiamin diphosphate (AP).

Insight to the Interaction of the Dihydrolipoamide Acetyltransferase (E2) Core with the Peripheral Components in the Escherichia coli Pyruvate Dehydrogenase Complex via Multifaceted Structural Approaches

Multifaceted structural approaches were undertaken to investigate interaction of the E2 component with E3 and E1 components from the Escherichia coli pyruvate dehydrogenase multienzyme complex (PDHc), as a representative of the PDHc from Gram-negative bacteria. The crystal structure of E3 at 2.5 Å resolution reveals similarity to other E3 structures and was an important starting point for understanding interaction surfaces between E3 and E2. Biochemical studies revealed that R129E-E2 and R150E-E2 substitutions in the peripheral subunit-binding domain (PSBD) of E2 greatly diminished PDHc activity, affected interactions with E3 and E1 components, and affected reductive acetylation of E2. Because crystal structures are unavailable for any complete E2-containing complexes, peptide-specific hydrogen/deuterium exchange mass spectrometry was used to identify loci of interactions between 3-lipoyl E2 and E3. Two peptides from the PSBD, including Arg-129, and three peptides from E3 displayed statistically significant reductions in deuterium uptake resulting from interaction between E3 and E2. Of the peptides identified on E3, two were from the catalytic site, and the third was from the interface domain, which for all known E3 structures is believed to interact with the PSBD. NMR clearly demonstrates that there is no change in the lipoyl domain structure on complexation with E3. This is the first instance where the entire wild-type E2 component was employed to understand interactions with E3. A model for PSBD-E3 binding was independently constructed and found to be consistent with the importance of Arg-129, as well as revealing other electrostatic interactions likely stabilizing this complex.

Investigation of Protein‐Protein Interactions in the Pyruvate Dehydrogenase Complex from Escherichia coli by Hydrogen/Deuterium Exchange Mass Spectrometry

The pyruvate dehydrogenase multienzyme complex (PDHc) from Escherichia coli is a molecular machine consisting of multiple copies of three enzyme components: pyruvate dehydrogenase (E1p), dihydrolipoyl acetyltransferase (E2p), and dihydrolipoyl dehydrogenase (E3p), which cooperate with each other to perform the oxidative decarboxylation of pyruvate to acetyl‐CoA. To probe the protein‐protein interactions and conformational changes between these components, we used hydrogen/deuterium exchange mass spectrometry (DXMS). Based on our study, the regions of interactions are mapped as follows: 1. The interaction of E2p with the active regions of E1p that were missing from previous X‐ray crystallography (an inner loop formed by residues 401–413, outer loop formed by residues 541–557 and N‐terminal residues 1–55); 2. The binding loci on peripheral‐subunit binding domain (PSBD) of E2p that interact with E1p and E3p; 3. The interaction of the E3p redox‐active region with E1p and E2p. This DXMS method has proved to be fast, accurate and reliable in characterizing the interaction of components in PDHc. The authors thank NIHGM‐050380 for financial support.

Interchain Acetyl Transfer in the E2 Component of Bacterial Pyruvate Dehydrogenase Suggests a Model with Different Roles for Each Chain in a Trimer of the Homooligomeric Component

The bacterial pyruvate dehydrogenase complex carries out conversion of pyruvate to acetyl-coenzyme A with the assistance of thiamin diphosphate (ThDP), several other cofactors, and three principal protein components, E1-E3, each present in multiple copies. The E2 component forms the core of the complexes, each copy consisting of variable numbers of lipoyl domains (LDs, lipoic acid covalently amidated at a lysine residue), peripheral subunit binding domains (PSBDs), and catalytic (or core) domains (CDs). The reaction starts with a ThDP-dependent decarboxylation on E1 to an enamine/C2α̃ carbanion, followed by oxidation and acetyl transfer to form S-acetyldihydrolipoamide E2, and then transfer of this acetyl group from the LD to coenzyme A on the CD. The dihydrolipoamide E2 is finally reoxidized by the E3 component. This report investigates whether the acetyl group is passed from the LD to the CD in an intra- or interchain reaction. Using an Escherichia coli E2 component having a single LD, two types of constructs were prepared: one with a Lys to Ala substitution in the LD at the Lys carrying the lipoic acid, making E2 incompetent toward post-translational ligation of lipoic acid and, hence, toward reductive acetylation, and the other in which the His believed to catalyze the transthiolacetylation in the CD is substituted with A or C, the absence of His rendering it incompetent toward acetyl-CoA formation. Both kinetic evidence and mass spectrometric evidence support interchain transfer of the acetyl groups, providing a novel model for the presence of multiples of three chains in all E2 components, and their assembly in bacterial enzymes.

Effects of Two Solvent Conditions on the Free Energy Landscape of the BBL Peripheral Subunit Binding Domain

BBL is a small independently folding domain with two main parallel helices. The experiment of C(α) secondary shifts has shown that changing the pH from ~7 to ~5 results in the reduced helicity at the C-terminus of helix 2. Combining constant pH molecular dynamics with replica exchange, we sampled the protein conformation space and protonation states extensively under a neutral pH condition and an acidic condition. Our results reveal that the solvent conditions influence the free energy landscape. Under the neutral pH condition, the denatured state and the native state are well separated. The condition of the acidic pH reshapes the free energy surface, leading to a broadly populated denatured-state basin and a low free energy barrier between the denatured state and the native state. The acidic pH shifts the equilibrium between the denatured state and the native state in favor of the denatured state. Caution must be used to interpret experimental data under the acidic condition because the contribution of the denatured state is significant. Our simulation results are supported by the fact that the calculated chemical shifts are in good agreement with the experiment data.

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.

Nuclear Magenetic Resonance Evidence for the Role of the Flexible Regions of the E1 Component of the Pyruvate Dehydrogenase Complex from Gram Negative Bacteria

Sequence‐specific NMR assignments were made for six residues in the N‐terminal 1–55 region, and one in each of the inner (401–413) and outer (541–557) active center loops of the E1 component of E. coli pyruvate dehydrogenase complex (E1ec, 2x100 kDa), flexible regions not seen in the x‐ray structure. The 2D HSQC spectrum revealed 103 resonances of the 109 expected in the three regions. Site‐specific substitutions and appropriate labeling patterns assisted the assignments. Resonances assigned to the N‐terminal region diminished in size/were absent in the spectrum in the presence of E2ec1‐190 comprising the lipoyl domain, peripheral subunit binding domain and hinge regions, indicating that the N‐terminal region became organized. Assignment of backbone NH resonances for residues G402 and G542 revealed dynamic features of the loops. These resonances, but not those assigned to the N‐terminal region, were no longer seen on saturation of the enzyme with a stable pre‐decarboxylation analogue C2alpha‐phosphonolactylThDP (PLThDP), indicating that the dynamic loop regions became ordered. The results provide the first structural support to the notion that the N‐terminal region of the E1ec of this entire class of bacterial pyruvate dehydrogenase complexes is responsible for binding the E2 component specifically, and dynamic motions of loops of E1ec are restricted by PLThDP. Supported by NIH‐GM‐050380.

Nuclear Magnetic Resonance Evidence for the Role of the Flexible Regions of the E1 Component of the Pyruvate Dehydrogenase Complex from Gram-negative Bacteria

Most bacterial pyruvate dehydrogenase complexes from either gram-positive or gram-negative bacteria have E1 components with an alpha(2) homodimeric quaternary structure. In a sequel to our previous publications, we present the first NMR study on the flexible regions of the E1 component from Escherichia coli and its biological relevance. We report sequence-specific NMR assignments for 6 residues in the N-terminal 1-55 region and for a glycine in each of the two mobile active center loops of the E1 component, a 200-kDa homodimer. This was accomplished by using site-specific substitutions and appropriate labeling patterns along with a peptide with the sequence corresponding to the N-terminal 1-35 amino acids of the E1 component. To study the functions of these mobile regions, we also examined the spectra in the presence of (a) a reaction intermediate analog known to affect the mobility of the active center loops, (b) an E2 component construct consisting of a lipoyl domain and peripheral subunit binding domain, and (c) a peptide corresponding to the amino acid sequence of the E2 peripheral subunit binding domain. Deductions from the NMR studies are in excellent agreement with our functional finding, providing a clear indication that the N-terminal region of the E1 interacts with the E2 peripheral subunit binding domain and that this interaction precedes reductive acetylation. The results provide the first structural support to the notion that the N-terminal region of the E1 component of this entire class of bacterial pyruvate dehydrogenase complexes is responsible for binding the E2 component.

Cooperative Folding Kinetics of BBL Protein and Peripheral Subunit-Binding Domain Homologues

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

Interplaying roles of native topology and chain length in marginally cooperative and noncooperative folding of small protein fragments

AbstractCoarse‐grained chain simulations were used to study fragments of two homologous proteins of the peripheral subunit‐binding domain (PSBD) family, Bacillus stearothermophilus PSBD (E3BD) and Escherichia coli 2‐oxo‐glutarate dehydrogenase PSBD (BBL). To ascertain a robust rank order of folding cooperativity, native‐centric intraprotein interactions were modeled by (i) a common Gō‐like potential, and (ii) native‐centric potentials with desolvation barriers or (iii) many‐body terms. Homologous proteins can possess substantially different folding cooperativity. Consistent with experiment, our calculations indicated that E3BD fragments fold more cooperatively than BBL fragments of approximately the same chain length. For a given fragment, native contacts deduced from Protein Data Bank structures can vary significantly depending on the number of residues that the structure encompasses in addition to those of the fragment itself, resulting in variation in model folding cooperativity predicted using different native contact sets for the same fragment. This observation underscores that folding cooperativity of these fragments can be extremely sensitive to change in chain length. Thus, a ∼40‐residue protein fragment's folding cooperativity, or lack thereof, does not necessarily imply essentially identical behaviors for super‐ or sub‐fragments with only several residues more, or less, than the given fragment. Ramifications for experimental investigations of downhill folding are discussed. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009

The Folding Mechanism of BBL: Plasticity of Transition-State Structure Observed within an Ultrafast Folding Protein Family

Studies on members of protein families with similar structures but divergent sequences provide insights into the effects of sequence composition on the mechanism of folding. Members of the peripheral subunit-binding domain (PSBD) family fold ultrafast and approach the smallest size for cooperatively folding proteins. Φ-Value analysis of the PSBDs E3BD and POB reveals folding via nucleation–condensation through structurally very similar, polarized transition states. Here, we present a Φ-value analysis of the family member BBL and found that it also folds by a nucleation–condensation mechanism. The mean Φ values of BBL, E3BD, and POB were near identical, indicating similar fractions of non-covalent interactions being formed in the transition state. Despite the overall conservation of folding mechanism in this protein family, however, the pattern of Φ values determined for BBL revealed a larger dispersion of the folding nucleus across the entire structure, and the transition state was less polarized. The observed plasticity of transition-state structure can be rationalized by the different helix-forming propensities of PSBD sequences. The very strong helix propensity in the first helix of BBL, relative to E3BD and POB, appears to recruit more structure formation in that helix in the transition state at the expense of weaker interactions in the second helix. Differences in sequence composition can modulate transition-state structure of even the smallest natural protein domains.