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

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 400 s-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 350 Hz) 'saturating' B1 fields. The 15N CEST experiment is very sensitive to exchange as the exchange delay TEX can be quite long (~0.5 s) 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 ~350 Hz) 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 ~104 s-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,500 s-1) from 33.2 to 42.9 °C while the unfolding rates (~70 to ~500 s-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 104 s-1 timescale can be studied using amide 15N CEST experiments.

GPU-Accelerated All-Atom Particle-Mesh Ewald Continuous Constant pH Molecular Dynamics in Amber

Constant pH molecular dynamics (MD) simulations sample protonation states on the fly according to the conformational environment and user specified pH conditions; however, the current accuracy is limited due to the use of implicit-solvent models or a hybrid solvent scheme. Here, we report the first GPU-accelerated implementation, parametrization, and validation of the all-atom continuous constant pH MD (CpHMD) method with particle-mesh Ewald (PME) electrostatics in the Amber22 pmemd.cuda engine. The titration parameters for Asp, Glu, His, Cys, and Lys were derived for the CHARMM c22 and Amber ff14sb and ff19sb force fields. We then evaluated the PME-CpHMD method using the asynchronous pH replica-exchange titration simulations with the c22 force field for six benchmark proteins, including BBL, hen egg white lysozyme (HEWL), staphylococcal nuclease (SNase), thioredoxin, ribonuclease A (RNaseA), and human muscle creatine kinase (HMCK). The root-mean-square deviation from the experimental pK_a's of Asp, Glu, His, and Cys is 0.76 pH units, and the Pearson's correlation coefficient for the pK_a shifts with respect to model values is 0.80. We demonstrated that a finite-size correction or much enlarged simulation box size can remove a systematic error of the calculated pK_a's and improve agreement with experiment. Importantly, the simulations captured the relevant biology in several challenging cases, e.g., the titration order of the catalytic dyad Glu35/Asp52 in HEWL and the coupled residues Asp19/Asp21 in SNase, the large pK_a upshift of the deeply buried catalytic Asp26 in thioredoxin, and the large pK_a downshift of the deeply buried catalytic Cys283 in HMCK. We anticipate that PME-CpHMD will offer proper pH control to improve the accuracies of MD simulations and enable mechanistic studies of proton-coupled dynamical processes that are ubiquitous in biology but remain poorly understood due to the lack of experimental tools and limitation of current MD simulations.

Structures and comparison of endogenous 2-oxoglutarate and pyruvate dehydrogenase complexes from bovine kidney

The α-keto acid dehydrogenase complex family catalyzes the essential oxidative decarboxylation of α-keto acids to yield acyl-CoA and NADH. Despite performing the same overarching reaction, members of the family have different component structures and structural organization between each other and across phylogenetic species. While native structures of α-keto acid dehydrogenase complexes from bacteria and fungi became available recently, the atomic structure and organization of their mammalian counterparts in native states remain unknown. Here, we report the cryo-electron microscopy structures of the endogenous cubic 2-oxoglutarate dehydrogenase complex (OGDC) and icosahedral pyruvate dehydrogenase complex (PDC) cores from bovine kidney determined at resolutions of 3.5 Å and 3.8 Å, respectively. The structures of multiple proteins were reconstructed from a single lysate sample, allowing direct structural comparison without the concerns of differences arising from sample preparation and structure determination. Although native and recombinant E2 core scaffold structures are similar, the native structures are decorated with their peripheral E1 and E3 subunits. Asymmetric sub-particle reconstructions support heterogeneity in the arrangements of these peripheral subunits. In addition, despite sharing a similar monomeric fold, OGDC and PDC E2 cores have distinct interdomain and intertrimer interactions, which suggests a means of modulating self-assembly to mitigate heterologous binding between mismatched E2 species. The lipoyl moiety lies near a mobile gatekeeper within the interdomain active site of OGDC E2 and PDC E2. Analysis of the twofold related intertrimer interface identified secondary structural differences and chemical interactions between icosahedral and cubic geometries of the core. Taken together, our study provides a direct structural comparison of OGDC and PDC from the same source and offers new insights into determinants of interdomain interactions and of architecture diversity among α-keto acid dehydrogenase complexes.

In Silico Modeling of the Influence of Environment on Amyloid Folding Using FOD-M Model

The role of the environment in amyloid formation based on the fuzzy oil drop model (FOD) is discussed here. This model assumes that the hydrophobicity distribution within a globular protein is consistent with a 3D Gaussian (3DG) distribution. Such a distribution is interpreted as the idealized effect of the presence of a polar solvent-water. A chain with a sequence of amino acids (which are bipolar molecules) determined by evolution recreates a micelle-like structure with varying accuracy. The membrane, which is a specific environment with opposite characteristics to the polar aquatic environment, directs the hydrophobic residues towards the surface. The modification of the FOD model to the FOD-M form takes into account the specificity of the cell membrane. It consists in "inverting" the 3DG distribution (complementing the Gaussian distribution), which expresses the exposure of hydrophobic residues on the surface. It turns out that the influence of the environment for any protein (soluble or membrane-anchored) is the result of a consensus factor expressing the participation of the polar environment and the "inverted" environment. The ratio between the proportion of the aqueous and the "reversed" environment turns out to be a characteristic property of a given protein, including amyloid protein in particular. The structure of amyloid proteins has been characterized in the context of prion, intrinsically disordered, and other non-complexing proteins to cover a wider spectrum of molecules with the given characteristics based on the FOD-M model.

Protein vibrations and their localization behaviour. A numerical scaling analysis

Using a classical force field, we investigate the localization properties of protein normal modes. For a set of eighteen proteins that cover five classes of increasing size, we compute the participation ratio as a measure of the spatial extent of protein vibrations. In this scaling analysis, we find extended low-frequency far-infrared and Terahertz modes, in contrast to localized high-frequency near-infrared vibrations. These regimes are separated by a broad crossover around a wave number of 260 cm^-1. Biophysical and biochemical implications are discussed, and the vibrational localization properties are compared to those of amorphous solids.

Downhill, Ultrafast and Fast Folding Proteins Revised

Research on the protein folding problem differentiates the protein folding process with respect to the duration of this process. The current structure encoded in sequence dogma seems to be clearly justified, especially in the case of proteins referred to as fast-folding, ultra-fast-folding or downhill. In the present work, an attempt to determine the characteristics of this group of proteins using fuzzy oil drop model is undertaken. According to the fuzzy oil drop model, a protein is a specific micelle composed of bi-polar molecules such as amino acids. Protein folding is regarded as a spherical micelle formation process. The presence of covalent peptide bonds between amino acids eliminates the possibility of free mutual arrangement of neighbors. An example would be the construction of co-micelles composed of more than one type of bipolar molecules. In the case of fast folding proteins, the amino acid sequence represents the optimal bipolarity system to generate a spherical micelle. In order to achieve the native form, it is enough to have an external force field provided by the water environment which directs the folding process towards the generation of a centric hydrophobic core. The influence of the external field can be expressed using the 3D Gaussian function which is a mathematical model of the folding process orientation towards the concentration of hydrophobic residues in the center with polar residues exposed on the surface. The set of proteins under study reveals a hydrophobicity distribution compatible with a 3D Gaussian distribution, taken as representing an idealized micelle-like distribution. The structure of the present hydrophobic core is also discussed in relation to the distribution of hydrophobic residues in a partially unfolded form.

How the dyes affect folding of small proteins in single-molecule FRET experiments: A simulation study

A key question in the application of the single-molecule Förster resonance energy transfer (smFRET) technique to study protein folding is how the dyes affect the protein behavior. Understanding of these effects is particularly important for small proteins, for which the dyes, along with their linkers, can be comparable in size (mass) with the protein. Using a coarse-grained model, we simulated folding of BBL protein and two of its FRET constructs. The obtained results suggest that even for small proteins, such as the 45-residue BBL, the appearance of the excluded volume in the protein conformation space due to the presence of dyes does not change the overall picture of folding. At the same time, some deviations from folding of the original protein are observed, in particular, the FRET constructs fold considerably slower than the original protein because the protein collapse in the initial state of folding is slowed down due to the protein loading with relatively massive dyes.

A Hierarchical Approach to Predict Conformation-Dependent Histidine Protonation States in Stable and Flexible Proteins

Solution acidity measured by pH is an important environmental factor that affects protein structure. It influences the protonation state of protein residues, which in turn may be coupled to protein conformational changes, unfolding, and ligand binding. It remains difficult to compute and measure the p K_a of individual residues, as well as to relate them to pH-dependent protein transitions. This paper presents a hierarchical approach to compute the p K_a of individual protonatable residues, specifically histidines, coupled with underlying structural changes of a protein. A fast and efficient free energy perturbation (FEP) algorithm has also been developed utilizing a fast implementation of standard molecular dynamics (MD) algorithms. Specifically, a CUDA version of the AMBER MD engine is used in this paper. Eight histidine p K_a's are computed in a diverse set of pH stable proteins to demonstrate the proposed approach's utility and assess the predictive quality of the AMBER FF99SB force field. A reference molecule is carefully selected and tested for convergence. A hierarchical approach is used to model p K_a's of the six histidine residues of the diphtheria toxin translocation domain (DTT), which exhibits a diverse ensemble of individual conformations and pH-dependent unfolding. The hierarchical approach consists of first sampling equilibrium conformational ensembles of a protein with protonated and neutral histidine residues via long equilibrium MD simulations (Flores-Canales, J. C.; et al. bioRxiv, 2019, 572040). A clustering method is then used to identify sampled protein conformations, and p K_a's of histidines in each protein conformation are computed. Finally, an ensemble averaging formalism is developed to compute weighted average histidine p K_a's. These can be compared with an apparent experimentally measured p K_a of the DTT protein and thus allows us to propose a mechanism of pH-dependent unfolding of the DTT protein.

Facile measurement of protein stability and folding kinetics using a nano differential scanning fluorimeter

With advancements in high-throughput generation of phenotypic data on mutant proteins, it has become important to individually characterize different proteins or their variants rapidly and with minimal sample consumption. We have made use of a nano differential scanning fluorimetric device, from NanoTemper technologies, to rapidly carry out isothermal chemical denaturation and measure folding/unfolding kinetics of proteins and compared these to corresponding data obtained from conventional spectrofluorimetry. We show that using sample volumes 10-50-fold lower than with conventional fluorimetric techniques, one can rapidly and accurately measure thermodynamic and kinetic stability, as well as folding/unfolding kinetics. This method also facilitates characterization of proteins that are difficult to express and purify.

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.

Modeling of Multicolor Single-Molecule Förster Resonance Energy-Transfer Experiments on Protein Folding

Using a coarse-grained, C_α-model of BBL protein, a multicolor single-molecule Förster resonance energy transfer (FRET) experiment is modeled. Three fluorophores are introduced, which, for simplicity, are associated with C_α beads. Two fluorophores are placed at the ends of protein chain and the third one at the middle of the chain. The free-energy surfaces (FESs) depending on the interfluorophore distances and on the FRET efficiencies corresponding to these distances have been constructed and compared with the FESs depending on the conventional collective variables, such as the fraction of native contacts and radius of gyration. It has been found that multicolor experiments can successfully resolve all essential BBL states that are revealed by the conventional FESs. The resolution of these states with the FRET-efficiency histogram is found to be successful if the energy transfer is measured between the fluorophores at the BBL ends. We also show that, although the present model construct of BBL is very simple, it captures some characteristic features of the single-molecule FRET experiments, such as the pattern of the FRET-efficiency histograms and their evolution with the denaturant concentration.

Generalized Born Based Continuous Constant pH Molecular Dynamics in Amber: Implementation, Benchmarking and Analysis

Solution pH plays an important role in structure and dynamics of biomolecular systems; however, pH effects cannot be accurately accounted for in conventional molecular dynamics simulations based on fixed protonation states. Continuous constant pH molecular dynamics (CpHMD) based on the λ-dynamics framework calculates protonation states on the fly during dynamical simulation at a specified pH condition. Here we report the CPU-based implementation of the CpHMD method based on the GBNeck2 generalized Born (GB) implicit-solvent model in the pmemd engine of the Amber molecular dynamics package. The performance of the method was tested using pH replica-exchange titration simulations of Asp, Glu and His side chains in 4 miniproteins and 7 enzymes with experimentally known p K_a's, some of which are significantly shifted from the model values. The added computational cost due to CpHMD titration ranges from 11 to 33% for the data set and scales roughly linearly as the ratio between the titrable sites and number of solute atoms. Comparison of the experimental and calculated p K_a's using 2 ns per replica sampling yielded a mean unsigned error of 0.70, a root-mean-squared error of 0.91, and a linear correlation coefficient of 0.79. Though this level of accuracy is similar to the GBSW-based CpHMD in CHARMM, in contrast to the latter, the current implementation was able to reproduce the experimental orders of the p K_a's of the coupled carboxylic dyads. We quantified the sampling errors, which revealed that prolonged simulation is needed to converge p K_a's of several titratable groups involved in salt-bridge-like interactions or deeply buried in the protein interior. Our benchmark data demonstrate that GBNeck2-CpHMD is an attractive tool for protein p K_a predictions.

pH-Dependent cooperativity and existence of a dry molten globule in the folding of a miniprotein BBL

Solution pH plays an important role in protein dynamics, stability, and folding; however, detailed mechanisms remain poorly understood. Here we use continuous constant pH molecular dynamics in explicit solvent with pH replica exchange to describe the pH profile of the folding cooperativity of a miniprotein BBL, which has drawn intense debate in the past. Our data reconciled the two opposing hypotheses (downhill vs. two-state) and uncovered a sparsely populated unfolding intermediate. As pH is lowered from 7 to 5, the folding barrier vanishes. As pH continues to decrease, the unfolding barrier lowers and denaturation is triggered by the protonation of Asp162, consistent with experimental evidence. Interestingly, unfolding proceeded via an intermediate, with intact secondary structure and a compact, unlocked hydrophobic core shielded from solvent, lending support to the recent hypothesis of a universal dry molten globule in protein folding. Our work demonstrates that constant pH molecular dynamics is a unique tool for testing this and other hypotheses to advance the knowledge in protein dynamics, stability, and folding.

Roles of conformational disorder and downhill folding in modulating protein-DNA recognition

Transcription factors are thought to efficiently search for their target DNA site via a combination of conventional 3D diffusion and 1D diffusion along the DNA molecule mediated by non-specific electrostatic interactions. This process requires the DNA-binding protein to quickly exchange between a search competent and a target recognition mode, but little is known as to how these two binding modes are encoded in the conformational properties of the protein. Here, we investigate this issue on the engrailed homeodomain (EngHD), a DNA-binding domain that folds ultrafast and exhibits a complex conformational behavior consistent with the downhill folding scenario. We explore the interplay between folding and DNA recognition using a coarse-grained computational model that allows us to manipulate the folding properties of the protein and monitor its non-specific and specific binding to DNA. We find that conformational disorder increases the search efficiency of EngHD by promoting a fast gliding search mode in addition to sliding. When gliding, EngHD remains loosely bound to DNA moving linearly along its length. A partially disordered EngHD also binds more dynamically to the target site, reducing the half-life of the specific complex via a spring-loaded mechanism. These findings apply to all conditions leading to partial disorder. However, we also find that at physiologically relevant temperatures EngHD is well folded and can only obtain the conformational flexibility required to accelerate 1D diffusion when it folds/unfolds within the downhill scenario (crossing a marginal free energy barrier). In addition, the conformational flexibility of native downhill EngHD enables its fast reconfiguration to lock into the specific binding site upon arrival, thereby affording finer control of the on- and off-rates of the specific complex. Our results provide key mechanistic insights into how DNA-binding domains optimize specific DNA recognition through the control of their conformational dynamics and folding mechanism.

Observations on AMBER Force Field Performance under the Conditions of Low pH and High Salt Concentrations

Molecular dynamics simulations have become an important tool for the study of structures, dynamics, and functions of biomolecules. Time scales and force field accuracy are key factors for the reliability of these simulations. With the advancement of computational platforms and simulation technologies, all-atom simulations of proteins in explicitly represented aqueous solutions can reach as long as milliseconds. However, there are indications of minor force field flaws in literature. In this work we present our observations on force field accuracy under uncommon conditions. We performed molecular dynamics simulations of BBL refolding in aqueous solutions of acidic pH and high salt concentrations (up to 6 M) with the AMBER99SB-ILDN force field for a microsecond time scale. The reliability of the simulations relies on the accuracy of the physical models of protein, water, and ions. Our simulations show the same trend as experiments: higher salt concentration facilities refolding. However, we have observed the presence of β-sheet in the native helical region as well as α-helix and β-sheet in the native loop region. Some of the nonnative secondary structures are even more stable than native helices. Aside from the secondary structure issue under the uncommon conditions, we have found that the rigidity of glycine dihedral angles in the loop region leads to low root-mean-square deviations but large topological differences from the native structure. Whether this is due to a force field deficiency or not needs further investigations. Recently developed ion parameters exhibit evident liquid features even in the 6 M LiCl solution. However, cation-anion interactions in TIP3P water still seem too strong, leading to high fractions of contact ion pairs. We do not find any specific ion-binding motif, thus we conclude that the effect of salt is a nonspecific electrostatic screening in our simulations. Our observations on the AMBER force field performance under acidic conditions and high salt concentrations show that simulations under extreme conditions can provide informative tests of physical models.

Kinetic and thermodynamic studies reveal chemokine homologues CC11 and CC24 with an almost identical tertiary structure have different folding pathways

BACKGROUND

Proteins with low sequence identity but almost identical tertiary structure and function have been valuable to uncover the relationship between sequence, tertiary structure, folding mechanism and functions. Two homologous chemokines, CCL11 and CCL24, with low sequence identity but similar tertiary structure and function, provide an excellent model system for respective studies.

RESULTS

The kinetics and thermodynamics of the two homologous chemokines were systematically characterized. Despite their similar tertiary structures, CCL11 and CCL24 show different thermodynamic stability in guanidine hydrochloride titration, with D_50% = 2.20 M and 4.96 M, respectively. The kinetics curves clearly show two phases in the folding/unfolding processes of both CCL11 and CCL24, which suggests the existence of an intermediate state in their folding/unfolding processes. The folding pathway of both CCL11 and CCL24 could be well described using a folding model with an on-pathway folding intermediate. However, the folding kinetics and stability of the intermediate state of CCL11 and CCL24 are obviously different.

CONCLUSION

Our results suggest homologous proteins with low sequence identity can display almost identical tertiary structure, but very different folding mechanisms, which applies to homologues in the chemokine protein family, extending the general applicability of the above observation.

How cooperative are protein folding and unfolding transitions?

A thermodynamically and kinetically simple picture of protein folding envisages only two states, native (N) and unfolded (U), separated by a single activation free energy barrier, and interconverting by cooperative two-state transitions. The folding/unfolding transitions of many proteins occur, however, in multiple discrete steps associated with the formation of intermediates, which is indicative of reduced cooperativity. Furthermore, much advancement in experimental and computational approaches has demonstrated entirely non-cooperative (gradual) transitions via a continuum of states and a multitude of small energetic barriers between the N and U states of some proteins. These findings have been instrumental towards providing a structural rationale for cooperative versus noncooperative transitions, based on the coupling between interaction networks in proteins. The cooperativity inherent in a folding/unfolding reaction appears to be context dependent, and can be tuned via experimental conditions which change the stabilities of N and U. The evolution of cooperativity in protein folding transitions is linked closely to the evolution of function as well as the aggregation propensity of the protein. A large activation energy barrier in a fully cooperative transition can provide the kinetic control required to prevent the accumulation of partially unfolded forms, which may promote aggregation. Nevertheless, increasing evidence for barrier-less "downhill" folding, as well as for continuous "uphill" unfolding transitions, indicate that gradual non-cooperative processes may be ubiquitous features on the free energy landscape of protein folding.

β-Structure within the Denatured State of the Helical Protein Domain BBL

Protein denatured states are the origin of both healthy and toxic conformational species. Denatured states of ultrafast folding proteins are of interest in mechanistic studies because they are energetically close to the kinetic bottleneck of folding. However, their transient nature makes them elusive to experiment. Here, we generated the denatured state of the helical domain BBL that is poised to fold in microseconds by a single-point mutation and combined circular dichroism spectroscopy, single-molecule fluorescence fluctuation analysis, and computer simulation to characterize its structure and dynamics. Circular dichroism showed a largely unfolded ensemble with marginal helix but significant β-sheet content. Main-chain structure and dynamics were unaffected by side-chain interactions that stabilize the native state, as revealed by site-directed mutagenesis and nanosecond loop closure kinetics probed by fluorescence correlation spectroscopy. Replica-exchange and constant-temperature molecular dynamics simulations showed a highly collapsed, hydrogen-bonded denatured state containing turn and β-sheet structure and few nucleating helices in an otherwise unfolded ensemble. An irregular β-hairpin element that connects helices in the native fold was poised to be formed. The surprising observation of β-structure in regions that form helices in the native state is reconciled by a generic low-energy pathway from the northwest quadrant of Ramachandran space to the helical basin present under folding conditions, proposed recently. Our results show that, indeed, rapid nucleation of helix emanates from β-structure formed early within a collapsed ensemble of unfolded conformers.

Protein folding of the SAP domain, a naturally occurring two-helix bundle

•

Thol SAP domain is one of the smallest model protein folding domains.

•

SAP domain folds through a diffuse transition state in which helix 1 is most formed.

•

Native state stability is dominated by contacts formed after the transition state.

The SAP domain from the Saccharomyces cerevisiae Tho1 protein is comprised of just two helices and a hydrophobic core and is one of the smallest proteins whose folding has been characterised. Φ-value analysis revealed that Tho1 SAP folds through a transition state where helix 1 is the most extensively formed element of secondary structure and flickering native-like core contacts from Leu35 are also present. The contacts that contribute most to native state stability of Tho1 SAP are not formed in the transition state.

Interaction Networks in Protein Folding via Atomic-Resolution Experiments and Long-Time-Scale Molecular Dynamics Simulations

The integration of atomic-resolution experimental and computational methods offers the potential for elucidating key aspects of protein folding that are not revealed by either approach alone. Here, we combine equilibrium NMR measurements of thermal unfolding and long molecular dynamics simulations to investigate the folding of gpW, a protein with two-state-like, fast folding dynamics and cooperative equilibrium unfolding behavior. Experiments and simulations expose a remarkably complex pattern of structural changes that occur at the atomic level and from which the detailed network of residue-residue couplings associated with cooperative folding emerges. Such thermodynamic residue-residue couplings appear to be linked to the order of mechanistically significant events that take place during the folding process. Our results on gpW indicate that the methods employed in this study are likely to prove broadly applicable to the fine analysis of folding mechanisms in fast folding proteins.

Protein unfolding rates correlate as strongly as folding rates with native structure

Although the folding rates of proteins have been studied extensively, both experimentally and theoretically, and many native state topological parameters have been proposed to correlate with or predict these rates, unfolding rates have received much less attention. Moreover, unfolding rates have generally been thought either to not relate to native topology in the same manner as folding rates, perhaps depending on different topological parameters, or to be more difficult to predict. Using a dataset of 108 proteins including two-state and multistate folders, we find that both unfolding and folding rates correlate strongly, and comparably well, with well-established measures of native topology, the absolute contact order and the long range order, with correlation coefficient values of 0.75 or higher. In addition, compared to folding rates, the absolute values of unfolding rates vary more strongly with native topology, have a larger range of values, and correlate better with thermodynamic stability. Similar trends are observed for subsets of different protein structural classes. Taken together, these results suggest that choosing a scaffold for protein engineering may require a compromise between a simple topology that will fold sufficiently quickly but also unfold quickly, and a complex topology that will unfold slowly and hence have kinetic stability, but fold slowly. These observations, together with the established role of kinetic stability in determining resistance to thermal and chemical denaturation as well as proteases, have important implications for understanding fundamental aspects of protein unfolding and folding and for protein engineering and design.

Elucidation of the interaction loci of the human pyruvate dehydrogenase complex E2·E3BP core with pyruvate dehydrogenase kinase 1 and kinase 2 by H/D exchange mass spectrometry and nuclear magnetic resonance

The human pyruvate dehydrogenase complex (PDC) comprises three principal catalytic components for its mission: E1, E2, and E3. The core of the complex is a strong subcomplex between E2 and an E3-binding protein (E3BP). The PDC is subject to regulation at E1 by serine phosphorylation by four kinases (PDK1–4), an inactivation reversed by the action of two phosphatases (PDP1 and -2). We report H/D exchange mass spectrometric (HDX-MS) and nuclear magnetic resonance (NMR) studies in the first attempt to define the interaction loci between PDK1 and PDK2 with the intact E2·E3BP core and their C-terminally truncated proteins. While the three lipoyl domains (L1 and L2 on E2 and L3 on E3BP) lend themselves to NMR studies and determination of interaction maps with PDK1 and PDK2 at the individual residue level, HDX-MS allowed studies of interaction loci on both partners in the complexes, PDKs, and other regions of the E2·E3BP core, as well, at the peptide level. HDX-MS suggested that the intact E2·E3BP core enhances the binding specificity of L2 for PDK2 over PDK1, while NMR studies detected lipoyl domain residues unique to interaction with PDK1 and PDK2. The E2·E3BP core induced more changes on PDKs than any C-terminally truncated protein, with clear evidence of greater plasticity of PDK1 than of PDK2. The effect of L1L2S paralleled HDX-MS results obtained with the intact E2·E3BP core; hence, L1L2S is an excellent candidate with which to define interaction loci with these two PDKs. Surprisingly, L3S′ induced moderate interaction with both PDKs according to both methods.

Novel microscale approaches for easy, rapid determination of protein stability in academic and commercial settings

Chemical denaturant titrations can be used to accurately determine protein stability. However, data acquisition is typically labour intensive, has low throughput and is difficult to automate. These factors, combined with high protein consumption, have limited the adoption of chemical denaturant titrations in commercial settings. Thermal denaturation assays can be automated, sometimes with very high throughput. However, thermal denaturation assays are incompatible with proteins that aggregate at high temperatures and large extrapolation of stability parameters to physiological temperatures can introduce significant uncertainties. We used capillary-based instruments to measure chemical denaturant titrations by intrinsic fluorescence and microscale thermophoresis. This allowed higher throughput, consumed several hundred-fold less protein than conventional, cuvette-based methods yet maintained the high quality of the conventional approaches. We also established efficient strategies for automated, direct determination of protein stability at a range of temperatures via chemical denaturation, which has utility for characterising stability for proteins that are difficult to purify in high yield. This approach may also have merit for proteins that irreversibly denature or aggregate in classical thermal denaturation assays. We also developed procedures for affinity ranking of protein–ligand interactions from ligand-induced changes in chemical denaturation data, and proved the principle for this by correctly ranking the affinity of previously unreported peptide–PDZ domain interactions. The increased throughput, automation and low protein consumption of protein stability determinations afforded by using capillary-based methods to measure denaturant titrations, can help to revolutionise protein research. We believe that the strategies reported are likely to find wide applications in academia, biotherapeutic formulation and drug discovery programmes.

•

Chemical denaturant titrations are slow, lengthy and consume lots of protein sample.

•

We developed very fast methods for measuring chemical denaturant titrations.

•

Automated titrations can be measured in minutes using only μl sample volumes.

•

Ligand interactions rapidly screened and ranked via changes in protein stability

•

These advances make chemical denaturation more suitable for commercial research.

The folding of a family of three-helix bundle proteins: spectrin R15 has a robust folding nucleus, unlike its homologous neighbours

Three homologous spectrin domains have remarkably different folding characteristics. We have previously shown that the slow-folding R16 and R17 spectrin domains can be altered to resemble the fast folding R15, in terms of speed of folding (and unfolding), landscape roughness and folding mechanism, simply by substituting five residues in the core. Here we show that, by contrast, R15 cannot be engineered to resemble R16 and R17. It is possible to engineer a slow-folding version of R15, but our analysis shows that this protein neither has a rougher energy landscape nor does change its folding mechanism. Quite remarkably, R15 appears to be a rare example of a protein with a folding nucleus that does not change in position or in size when its folding nucleus is disrupted. Thus, while two members of this protein family are remarkably plastic, the third has apparently a restricted folding landscape.

Fast protein folding kinetics

Fast-folding proteins have been a major focus of computational and experimental study because they are accessible to both techniques: they are small and fast enough to be reasonably simulated with current computational power, but have dynamics slow enough to be observed with specially developed experimental techniques. This coupled study of fast-folding proteins has provided insight into the mechanisms, which allow some proteins to find their native conformation well <1 ms and has uncovered examples of theoretically predicted phenomena such as downhill folding. The study of fast folders also informs our understanding of even 'slow' folding processes: fast folders are small; relatively simple protein domains and the principles that govern their folding also govern the folding of more complex systems. This review summarizes the major theoretical and experimental techniques used to study fast-folding proteins and provides an overview of the major findings of fast-folding research. Finally, we examine the themes that have emerged from studying fast folders and briefly summarize their application to protein folding in general, as well as some work that is left to do.

Spatial ranges of driving forces are a key determinant of protein folding cooperativity and rate diversity

The physical basis of two-state-like folding transitions and the tremendous diversity in folding rates is elucidated by directly simulating the folding kinetics of 52 representative proteins. Relative to the results from a common modeling approach, the diversity of the simulated folding rates can be increased from ~10(2.1) to the experimental ~10(6.0) by a modest decrease in the spatial range of the attractive potential. The required theoretical range is consistent with desolvation physics and is notably much more permissive than that needed for two-state-like homopolymer collapse.

Nuclear magnetic resonance approaches in the study of 2-oxo acid dehydrogenase multienzyme complexes--a literature review

The 2-oxoacid dehydrogenase complexes (ODHc) consist of multiple copies of three enzyme components: E1, a 2-oxoacid decarboxylase; E2, dihydrolipoyl acyl-transferase; and E3, dihydrolipoyl dehydrogenase, that together catalyze the oxidative decarboxylation of 2-oxoacids, in the presence of thiamin diphosphate (ThDP), coenzyme A (CoA), Mg²⁺ and NAD⁺, to generate CO₂, NADH and the corresponding acyl-CoA. The structural scaffold of the complex is provided by E2, with E1 and E3 bound around the periphery. The three principal members of the family are pyruvate dehydrogenase (PDHc), 2-oxoglutarate dehydrogenase (OGDHc) and branched-chain 2-oxo acid dehydrogenase (BCKDHc). In this review, we report application of NMR-based approaches to both mechanistic and structural issues concerning these complexes. These studies revealed the nature and reactivity of transient intermediates on the enzymatic pathway and provided site-specific information on the architecture and binding specificity of the domain interfaces using solubilized truncated domain constructs of the multi-domain E2 component in its interactions with the E1 and E3 components. Where studied, NMR has also provided information about mobile loops and the possible relationship of mobility and catalysis.

Exact partition function zeros of the Wako-Saitô-Muñoz-Eaton protein model

I compute exact partition function zeros of the Wako-Saitô-Muñoz-Eaton model for various secondary structural elements and for two proteins, 1BBL and 1I6C, by using both analytic and numerical methods. Two-state and barrierless downhill folding transitions can be distinguished by a gap in the distribution of zeros at the positive real axis.

Using VIPT-jump to distinguish between different folding mechanisms: application to BBL and a Trpzip

Protein folding involves a large number of sequential molecular steps or conformational substates. Thus, experimental characterization of the underlying folding energy landscape for any given protein is difficult. Herein, we present a new method that can be used to determine the major characteristics of the folding energy landscape in question, e.g., to distinguish between activated and barrierless downhill folding scenarios. This method is based on the idea that the conformational relaxation kinetics of different folding mechanisms at a given final condition will show different dependences on the initial condition. We show, using both simulation and experiment, that it is possible to differentiate between disparate kinetic folding models by comparing temperature jump (T-jump) relaxation traces obtained with a fixed final temperature and varied initial temperatures, which effectively varies the initial potential (VIP) of the system of interest. We apply this method (hereafter refer to as VIPT-jump) to two model systems, tryptophan zipper (Trpzip)-2c and BBL, and our results show that BBL exhibits characteristics of barrierless downhill folding, whereas Trpzip-2c folding encounters a free energy barrier. In addition, using the T-jump data of BBL we are able to provide, via Langevin dynamics simulations, a realistic estimate of its conformational diffusion coefficient.

Mutational effects on the folding dynamics of a minimized hairpin

The fold stabilities and folding dynamics of a series of mutants of a model hairpin, KTW-NPATGK-WTE (HP7), are reported. The parent system and the corresponding DPATGK loop species display submicrosecond folding time constants. The mutational studies revealed that ultrafast folding requires both some prestructuring of the loop and a favorable interaction between the chain termini in the transition state. In the case of YY-DPETGT-WY, another submicrosecond folding species [Davis, C. M., Xiao, S., Raleigh, D. P., and Dyer, R. B. (2012) J. Am. Chem. Soc. 134, 14476-14482], a hydrophobic cluster provides the latter. In the case of HP7, the Coulombic interaction between the terminal NH3(+) and CO2(-) units provides this; a C-terminal Glu to amidated Ala mutation results in a 5-fold retardation of the folding rate. The effects of mutations within the reversing loop indicate the balance between loop flexibility (favoring fast conformational searching) and turn formation in the unfolded state is a major factor in determining the folding dynamics. The -NAAAKX- loops examined display no detectable turn formation propensity in other hairpin constructs but do result in stable analogues of HP7. Peptide KTW-NAAAKK-WTE displays the same fold stability as HP7, but both the folding and unfolding time constants are greater by a factor of 20.

Low folding cooperativity of HP35 revealed by single-molecule force spectroscopy and molecular dynamics simulation

Some small proteins, such as HP35, fold at submicrosecond timescale with low folding cooperativity. Although these proteins have been extensively investigated, still relatively little is known about their folding mechanism. Here, using single-molecule force spectroscopy and steered molecule dynamics simulation, we study the unfolding of HP35 under external force. Our results show that HP35 unfolds at extremely low forces without a well-defined unfolding transition state. Subsequently, we probe the structure of unfolded HP35 using the persistence length obtained in the force spectroscopy. We found that the persistence length of unfolded HP35 is around 0.72 nm, >40% longer than typical unstructured proteins, suggesting that there are a significant amount of residual secondary structures in the unfolded HP35. Molecular dynamics simulation further confirmed this finding and revealed that many native contacts are preserved in HP35, even its two ends have been extended up to 8 nm. Our results therefore suggest that retaining a significant amount of secondary structures in the unfolded state of HP35 may be an efficient way to reduce the entropic cost for the formation of tertiary structure and increase the folding speed, although the folding cooperativity is compromised. Moreover, we anticipate that the methods we used in this work can be extended to the study of other proteins with complex folding behaviors and even intrinsically disordered ones.

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.

Observation of two families of folding pathways of BBL

BBL is an independent folding domain of a large multienzyme complex, 2-oxoglutarate dehydrogenase. The folding mechanism of BBL is under debate between the views of noncooperative downhill-type and classical two-state. Extensive replica exchange molecular dynamics simulations of BBL in explicit solvent have shown some non-two-state behaviors despite no definitive evidence of downhill folding. In this work, we postprocess the replica exchange data using our roadmap-based MaxFlux reaction path algorithm to reveal atomically detailed folding pathways. A connected graph is used to organize and visualize the folding pathways initiated from random coils. High structural and transition heterogeneity is seen in the early stage of folding. Two main parallel folding pathways emerge in the later stage; one path shows that tertiary contact and helix formation develop at different stages of folding, whereas the other path exhibits concurrence of secondary and tertiary structure formation to some extent. Because the native state of BBL is sensitive to experimental conditions, we speculate that the relative predominance of the two pathways may vary with the protein construct and solvent conditions, possibly leading to the seeming discrepancy of experimental results. Our roadmap-based reaction path algorithm is a general tool to extract path information from replica exchange.

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.

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

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

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.

Folding simulations of a de novo designed protein with a betaalphabeta fold

betaalphabeta structural motifs are commonly used building blocks in protein structures containing parallel beta-sheets. However, to our knowledge, no stand-alone betaalphabeta structure has been observed in nature to date. Recently, for the first time that we know of, a small protein with an independent betaalphabeta structure (DS119) was successfully designed in our laboratory. To understand the folding mechanism of DS119, in the study described here, we carried out all-atom molecular dynamics and coarse-grained simulations to investigate its folding pathways and energy landscape. From all-atom simulations, we successfully observed the folding event and got a stable folded structure with a minimal root mean-square deviation of 2.6 A with respect to the NMR structure. The folding process can be described as a fast collapse phase followed by rapid formation of the central helix, and then slow formation of a parallel beta-sheet. By using a native-centric Gō-like model, the cooperativity of the system was characterized in terms of the calorimetric criterion, sigmoidal transitions, conformation distribution shifts, and free-energy profiles. DS119 was found to be an incipient downhill folder that folds more cooperatively than a downhill folder, but less cooperatively than a two-state folder. This may reflect the balance between the two structural elements of DS119: the rapidly formed alpha-helix and the slowly formed parallel beta-sheet. Folding times estimated from both the all-atom simulations and the coarse-grained model were at microsecond level, making DS119 another fast folder. Compared to fast folders reported previously, DS119 is, to the best of our knowledge, the first that exhibits a parallel beta-sheet.

What lessons can be learned from studying the folding of homologous proteins?

The studies of the folding of structurally related proteins have proved to be a very important tool for investigating protein folding. Here we review some of the insights that have been gained from such studies. Our highlighted studies show just how such an investigation should be designed and emphasise the importance of the synergy between experiment and theory. We also stress the importance of choosing the right system carefully, exploiting the excellent structural and sequence databases at our disposal.

The folding energy landscape and free energy excitations of cytochrome c

The covalently bound heme cofactor plays a dominant role in the folding of cytochrome c. Because of the complicated inorganic chemistry of the heme, some might consider the folding of cytochrome c to be a special case, following principles different from those used to describe the folding of proteins without cofactors. Recent investigations, however, demonstrate that common models describing folding for many proteins work well for cytochrome c when heme is explicitly introduced, generally providing results that agree with experimental observations. In this Account, we first discuss results from simple native structure-based models. These models include attractive interactions between nonadjacent residues only if they are present in the crystal structure at pH 7. Because attractive nonnative contacts are not included in native structure-based models, their energy landscapes can be described as "perfectly funneled". In other words, native structure-based models are energetically guided towards the native state and contain no energetic traps that would hinder folding. Energetic traps are denoted sources of "frustration", which cause specific transient intermediates to be populated. Native structure-based models do, however, include repulsion between residues due to excluded volume. Nonenergetic traps can therefore exist if the chain, which cannot cross over itself, must partially unfold so that folding can proceed. The ability of native structure-based models to capture this kind of motion is partly responsible for their successful predictions of folding pathways for many types of proteins. Models without frustration describe the sequence of folding events for cytochrome c well (as inferred from hydrogen-exchange experiments), thereby justifying their use as a starting point. At low pH, the experimentally observed folding sequence of cytochrome c deviates from that at pH 7 and from models with perfectly funneled energy landscapes. Here, alternate folding pathways are a result of "chemical frustration". This frustration arises because some regions of the protein are destabilized more than others due to the heterogeneous distribution of titratable residues that are protonated at low pH. Beginning with native structure-based terms, we construct more complex models by adding chemical frustration. These more complex models only modestly perturb the energy landscape, which remains, overall, well funneled. These perturbed models can accurately describe how alternative folding pathways are used at low pH. At alkaline pH, cytochrome c populates distinctly different structural ensembles. For instance, lysine residues are deprotonated and compete for the heme ligation site. The same models that can describe folding at low pH also predict well the structures and relative stabilities of intermediates populated at alkaline pH. The success of models based on funneled energy landscapes suggest that cytochrome c folding is driven primarily by native contacts. The presence of heme appears to add chemical complexity to the folding process, but it does not require fundamental modification of the general principles used to describe folding. Moreover, its added complexity provides a valuable means of probing the folding energy landscape in greater detail than is possible with simpler systems.

Phi-value analysis for ultrafast folding proteins by NMR relaxation dispersion

Proteins that fold rapidly, on the (sub-) microsecond time scale, offer the prospect of direct comparison between experimental data and molecular dynamics simulations. However, experimental studies for such proteins frequently are hindered because folding rates are too fast to measure using conventional stopped-flow methods. To overcome this impediment, NMR spin relaxation dispersion experiments are used to quantify mutational effects on kinetics (DeltaDeltaG(o)), stability (DeltaDeltaG(o)), and phi-values (DeltaDeltaG(dagger)/DeltaDeltaG(o)) for proteins exhibiting chemical exchange line broadening that is fast on the NMR chemical shift time scale. The accuracy of phi-value analysis is enhanced because mutational effects on denatured or intermediate states can be detected through changes in line broadening. The transition and intermediate states of the villin headpiece domain, HP67, are characterized in varying solvent conditions to validate the method.

Engineering a two-helix bundle protein for folding studies

The SAP domain from the Saccharomyces cerevisiae THO1 protein contains a hydrophobic core and just two α-helices. It could provide a system for studying protein folding that bridges the gap between studies on isolated helices and those on larger protein domains. We have engineered the SAP domain for protein folding studies by inserting a tryptophan residue into the hydrophobic core (L31W) and solved its structure. The helical regions had a backbone root mean-squared deviation of 0.9 Å from those of wild type. The mutation L31W destabilised wild type by 0.8 ± 0.1 kcal mol⁻¹. The mutant folded in a reversible, apparent two-state manner with a microscopic folding rate constant of around 3700 s⁻¹ and is suitable for extended studies of folding.

Understanding protein folding cooperativity based on topological consideration

The folding cooperativity is an important issue of protein folding dynamics. Since the native topology plays a significant role in determining the folding behavior of proteins, we believe that it also has close relationship with the folding cooperativity. In the present work, we perform simulations on proteins Naf-BBL, QNND-BBL, CI2, and SH3 with the Gō model and compare their different folding behaviors. By analyzing the weak cooperative folding of protein Naf-BBL in detail, we found that the folding of Naf-BBL shows relatively weak thermodynamic coupling between residues, and such weak coupling is found mainly between the nonlocal native contacts. This finding complements our understandings on the source of barrierless folding of Naf-BBL and promotes us to analyze the topological origins of the poor thermodynamic coupling of Naf-BBL. Then, we further extend our analysis to other two-state and multistate proteins. Based on the considerations of the thermodynamic coupling and kinetic coupling, we conclude that the fraction of scattered native contacts, the difference in loop entropy of contacts, and the long range relative contact order are the major topological factors that influence the folding cooperativity. The combination of these three tertiary structural features shows significant correlations with the folding types of proteins. Moreover, we also discuss the topological factors related to downhill folding. Finally, the generic role of tertiary structure in determining the folding cooperativity is summarized.

Direct observation of ultrafast folding and denatured state dynamics in single protein molecules

Single-molecule fluorescence resonance energy transfer (smFRET) experiments are extremely useful in studying protein folding but are generally limited to time scales of greater than approximately 100 micros and distances greater than approximately 2 nm. We used single-molecule fluorescence quenching by photoinduced electron transfer, detecting short-range events, in combination with fluorescence correlation spectroscopy (PET-FCS) to investigate folding dynamics of the small binding domain BBL with nanosecond time resolution. The kinetics of folding appeared as a 10-micros decay in the autocorrelation function, resulting from stochastic fluctuations between denatured and native conformations of individual molecules. The observed rate constants were probe independent and in excellent agreement with values derived from conventional temperature-jump (T-jump) measurements. A submicrosecond relaxation was detected in PET-FCS data that reported on the kinetics of intrachain contact formation within the thermally denatured state. We engineered a mutant of BBL that was denatured under the reaction conditions that favored folding of the parent wild type ("D(phys)"). D(phys) had the same kinetic signature as the thermally denatured state and revealed segmental diffusion with a time constant of intrachain contact formation of 500 ns. This time constant was more than 10 times faster than folding and in the range estimated to be the "speed limit" of folding. D(phys) exhibited significant deviations from a random coil. The solvent viscosity and temperature dependence of intrachain diffusion showed that chain motions were slaved by the presence of intramolecular interactions. PET-FCS in combination with protein engineering is a powerful approach to study the early events and mechanism of ultrafast protein folding.

Native like structure in the unfolded state of the villin headpiece helical subdomain, an ultrafast folding protein

The villin headpiece subdomain, HP36, is the smallest naturally occurring protein that folds cooperatively. Its small size, rapid folding, and simple three-helix topology have made it an extremely popular system for computational studies of protein folding. The role of unfolded state structure in rapid folding is an area of active investigation, but relatively little is known about the properties of unfolded states under native conditions. A peptide fragment, HP21, which contains the first and second helices of HP36 has been shown to be a good model for structure in the unfolded state of the intact domain but a detailed description of the conformational propensities of HP21 is lacking and the balance between native and nonnative interactions is not known. A series of three-dimensional NMR experiments were performed on (13)C, (15)N-labeled HP21 to investigate in detail its conformational propensities. Analysis of (13)C(alpha), (13)C(beta), (13)CO chemical shifts, Deltadelta(13)C(alpha) - Deltadelta(13)C(beta) secondary shifts, the secondary structure propensity scores, NOEs, (15)N R(2) values and comparison of experimental chemical shifts with those of HP36 and with chemical shifts calculated using the SHIFTS and SHIFTX programs all indicate that there is significant native like structure in the HP21 ensemble, and thus by implication in the unfolded state of HP36.

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.