Kinetic Studies of the Folding of Heterodimeric Monellin: Evidence for Switching between Alternative Parallel Pathways

Understanding the heterogeneity intrinsic to protein folding

Relating the native fold of a protein to its amino acid sequence remains a fundamental problem in biology. While computer algorithms have demonstrated recently their prowess in predicting what structure a particular amino acid sequence will fold to, an understanding of how and why a specific protein fold is achieved remains elusive. A major challenge is to define the role of conformational heterogeneity during protein folding. Recent experimental studies, utilizing time-resolved FRET, hydrogen-exchange coupled to mass spectrometry, and single-molecule force spectroscopy, often in conjunction with simulation, have begun to reveal how conformational heterogeneity evolves during folding, and whether an intermediate ensemble of defined free energy consists of different sub-populations of molecules that may differ significantly in conformation, energy and entropy.

Differentiating between the sequence of structural events on alternative pathways of folding of a heterodimeric protein

Distinguishing between competing pathways of folding of a protein, on the basis of how they differ in their progress of structure acquisition, remains an important challenge in protein folding studies. A previous study had shown that the heterodimeric protein, double chain monellin (dcMN) switches between alternative folding pathways upon a change in guanidine hydrochloride (GdnHCl) concentration. In the current study, the folding of dcMN has been characterized by the pulsed hydrogen exchange (HX) labeling methodology used in conjunction with mass spectrometry. Quantification of the extent to which folding intermediates accumulate and then disappear with time of folding at both low and high GdnHCl concentrations, where the folding pathways are known to be different, shows that the folding mechanism is describable by a triangular three-state mechanism. Structural characterization of the productive folding intermediates populated on the alternative pathways has enabled the pathways to be differentiated on the basis of the progress of structure acquisition that occurs on them. The intermediates on the two pathways differ in the extent to which the α-helix and the rest of the β-sheet have acquired structure that is protective against HX. The major difference is, however, that β2 has not acquired any protective structure in the intermediate formed on one pathway, but it has acquired significant protective structure in the intermediate formed on the alternative pathway. Hence, the sequence of structural events is different on the two alternative pathways.

Structural Characterization of the Cooperativity of Unfolding of a Heterodimeric Protein using Hydrogen Exchange-Mass Spectrometry

Little is known about how the sequence of structural changes in one chain of a heterodimeric protein is coupled to those in the other chain during protein folding and unfolding reactions, and whether individual secondary structural changes in the two chains occur in one or many coordinated steps. Here, the unfolding mechanism of a small heterodimeric protein, double chain monellin, has been characterized using hydrogen exchange-mass spectrometry. Transient structure opening, which enables HX, was found to be describable by a five state N ↔ I₁ ↔ I₂ ↔ I₃ ↔ U mechanism. Structural changes occur gradually in the first three steps, and cooperatively in the last step. β strands 2, 4 and 5, as well as the α-helix undergo transient unfolding during all three non-cooperative steps, while β1 and the two loops on both sides of the helix undergo transient unfolding during the first two steps. In the absence of GdnHCl, only β3 in chain A of the protein unfolds during the last cooperative step, while in the presence of 1 M GdnHCl, not only β3, but also β2 in chain B unfolds cooperatively. Hence, the extent of cooperative structural change and size of the cooperative unfolding unit increase when the protein is destabilized by denaturant. The naturally evolved two-chain variant of monellin folds and unfolds in a more cooperative manner than does a single chain variant created artificially, suggesting that increasing folding cooperativity, even at the cost of decreasing stability, may be a driving force in the evolution of proteins.

Mapping Distinct Sequences of Structure Formation Differentiating Multiple Folding Pathways of a Small Protein

To determine experimentally how the multiple folding pathways of a protein differ, in the order in which the structural parts are assembled, has been a long-standing challenge. To resolve whether structure formation during folding can progress in multiple ways, the complex folding landscape of monellin has been characterized, structurally and temporally, using the multisite time-resolved FRET methodology. After an initial heterogeneous polypeptide chain collapse, structure formation proceeds on parallel pathways. Kinetic analysis of the population evolution data across various protein segments provides a clear structural distinction between the parallel pathways. The analysis leads to a phenomenological model that describes how and when discrete segments acquire structure independently of each other in different subensembles of protein molecules. When averaged over all molecules, structure formation is seen to progress as α-helix formation, followed by core consolidation, then β-sheet formation, and last end-to-end distance compaction. Parts of the protein that are closer in the primary sequence acquire structure before parts separated by longer sequence.

Intrinsic Disorder in a Well-Folded Globular Protein

The folded structure of the heterodimeric sweet protein monellin mimics single-chain proteins with topology β1-α1-β2-β3-β4-β5 (chain A: β3-β4-β5; chain B: β1-α1-β2). Furthermore, like naturally occurring single-chain proteins of a similar size, monellin folds cooperatively with no detectable intermediates. However, the two monellin chains, A and B, are marginally structured in isolation and fold only upon binding to each other. Thus, monellin presents a unique opportunity to understand the design of intrinsically disordered proteins that fold upon binding. Here, we study the folding of a single-chain variant of monellin (scMn) using simulations of an all heavy-atom structure-based model. These simulations can explain mechanistic details derived from scMn experiments performed using several different structural probes. scMn folds cooperatively in our structure-based simulations, as is also seen in experiments. We find that structure formation near the transition-state ensemble of scMn is not uniformly distributed but is localized to a hairpin-like structure which contains one strand from each chain (β2, β3). Thus, the sequence and the underlying energetics of heterodimeric monellin promote the early formation of the interchain interface (β2-β3). By studying computational scMn mutants whose "interchain" interactions are deleted, we infer that this energy distribution allows the two protein chains to remain largely disordered when this interface is not folded. From these results, we suggest that cutting the protein backbone of a globular protein between residues which lie within its folding nucleus may be one way to construct two disordered fragments which fold upon binding.

Site-specific time-resolved FRET reveals local variations in the unfolding mechanism in an apparently two-state protein unfolding transition

Using multi-site time-resolved FRET, it is shown that equilibrium unfolding of monellin is not only heterogeneous, but that the degree of non-cooperativity differs between the sole α-helix and different parts of the β-sheet.

Homodimeric Escherichia coli Toxin CcdB (Controller of Cell Division or Death B Protein) Folds via Parallel Pathways

The existence of parallel pathways in the folding of proteins seems intuitive, yet remains controversial. We explore the folding kinetics of the homodimeric Escherichia coli toxin CcdB (Controller of Cell Division or Death B protein) using multiple optical probes and approaches. Kinetic studies performed as a function of protein and denaturant concentrations demonstrate that the folding of CcdB is a four-state process. The two intermediates populated during folding are present on parallel pathways. Both form by rapid association of the monomers in a diffusion limited manner and appear to be largely unstructured, as they are silent to the optical probes employed in the current study. The existence of parallel pathways is supported by the insensitivity of the amplitudes of the refolding kinetic phases to the different probes used in the study. More importantly, interrupted refolding studies and ligand binding studies clearly demonstrate that the native state forms in a biexponential manner, implying the presence of at least two pathways. Our studies indicate that the CcdA antitoxin binds only to the folded CcdB dimer and not to any earlier folding intermediates. Thus, despite being part of the same operon, the antitoxin does not appear to modulate the folding pathway of the toxin encoded by the downstream cistron. This study highlights the utility of ligand binding in distinguishing between sequential and parallel pathways in protein folding studies, while also providing insights into molecular interactions during folding in Type II toxin-antitoxin systems.

Molecular Dynamics Driven Design of pH-Stabilized Mutants of MNEI, a Sweet Protein

MNEI is a single chain derivative of monellin, a plant protein that can interact with the human sweet taste receptor, being therefore perceived as sweet. This unusual physiological activity makes MNEI a potential template for the design of new sugar replacers for the food and beverage industry. Unfortunately, applications of MNEI have been so far limited by its intrinsic sensitivity to some pH and temperature conditions, which could occur in industrial processes. Changes in physical parameters can, in fact, lead to irreversible protein denaturation, as well as aggregation and precipitation. It has been previously shown that the correlation between pH and stability in MNEI derives from the presence of a single glutamic residue in a hydrophobic pocket of the protein. We have used molecular dynamics to study the consequences, at the atomic level, of the protonation state of such residue and have identified the network of intramolecular interactions responsible for MNEI stability at acidic pH. Based on this information, we have designed a pH-independent, stabilized mutant of MNEI and confirmed its increased stability by both molecular modeling and experimental techniques.

Force-dependent switch in protein unfolding pathways and transition-state movements

Although it is known that single-domain proteins fold and unfold by parallel pathways, demonstration of this expectation has been difficult to establish in experiments. Unfolding rate, [Formula: see text], as a function of force f, obtained in single-molecule pulling experiments on src SH3 domain, exhibits upward curvature on a [Formula: see text] plot. Similar observations were reported for other proteins for the unfolding rate [Formula: see text]. These findings imply unfolding in these single-domain proteins involves a switch in the pathway as f or [Formula: see text] is increased from a low to a high value. We provide a unified theory demonstrating that if [Formula: see text] as a function of a perturbation (f or [Formula: see text]) exhibits upward curvature then the underlying energy landscape must be strongly multidimensional. Using molecular simulations we provide a structural basis for the switch in the pathways and dramatic shifts in the transition-state ensemble (TSE) in src SH3 domain as f is increased. We show that a single-point mutation shifts the upward curvature in [Formula: see text] to a lower force, thus establishing the malleability of the underlying folding landscape. Our theory, applicable to any perturbation that affects the free energy of the protein linearly, readily explains movement in the TSE in a β-sandwich (I27) protein and single-chain monellin as the denaturant concentration is varied. We predict that in the force range accessible in laser optical tweezer experiments there should be a switch in the unfolding pathways in I27 or its mutants.