Exploring subdomain cooperativity in T4 lysozyme I: structural and energetic studies of a circular permutant and protein fragment

Protein engineering using circular permutation - structure, function, stability, and applications

Protein engineering is important for creating novel variants from natural proteins, enabling a wide range of applications. Approaches such as rational design and directed evolution are routinely used to make new protein variants. Computational tools like de novo design can introduce new protein folds. Expanding the amino acid repertoire to include unnatural amino acids with non-canonical side chains in vitro by native chemical ligation and in vivo via codon expansion methods broadens sequence and structural possibilities. Circular permutation (CP) is an invaluable approach to redesigning a protein by rearranging the amino acid sequence, where the connectivity of the secondary structural elements is altered without changing the overall structure of the protein. Artificial CP proteins (CPs) are employed in various applications such as biocatalysis, sensing of small molecules by fluorescence, genome editing, ligand-binding protein switches, and optogenetic engineering. Many studies have shown that CP can lead to either reduced or enhanced stability or catalytic efficiency. The effects of CP on a protein's energy landscape cannot be predicted a priori. Thus, it is important to understand how CP can affect the thermodynamic and kinetic stability of a protein. In this review, we discuss the discovery and advancement of techniques to create protein CP, and existing reviews on CP. We delve into the plethora of biological applications for designed CP proteins. We subsequently discuss the experimental and computational reports on the effects of CP on the thermodynamic and kinetic stabilities of proteins of various topologies. An understanding of the various aspects of CP will allow the reader to design robust CP proteins for their specific purposes.

Circular permutation at azurin's active site slows down its folding

Circular permutation (CP) is a technique by which the primary sequence of a protein is rearranged to create new termini. The connectivity of the protein is altered but the overall protein structure generally remains unperturbed. Understanding the effect of CP can help design robust proteins for numerous applications such as in genetic engineering, optoelectronics, and improving catalytic activity. Studies on different protein topologies showed that CP usually affects protein stability as well as unfolding rates. Though a significant number of proteins contain metals or other cofactors, reports of metalloprotein CPs are rare. Thus, we chose a bacterial metalloprotein, azurin, and its CP within the metal-binding site (cpF114). We studied the stabilities, folding, and unfolding rates of apo- and Zn2+-bound CP azurin using fluorescence and circular dichroism. The introduced CP had destabilizing effects on the protein. Also, the folding of the Zn2+-CP protein was much slower than that of the Zn2+-WT or apo-protein. We compared this study to our previously reported azurin-cpN42, where we had observed an equilibrium and kinetic intermediate. cpF114 exhibits an apparent two-state equilibrium unfolding but has an off-pathway kinetic intermediate. Our study hinted at CP as a method to modify the energy landscape of proteins to alter their folding pathways. WT azurin, being a faster folder, may have evolved to optimize the folding rate of metal-bound protein compared to its CPs, albeit all of them have the same structure and function. Our study underscores that protein sequence and protein termini positions are crucial for metalloproteins. TOC Figure. (Top) Zn2+-azurin WT structure (PDB code: 1E67) and 2-D topology diagram of Zn2+-cpF114 azurin. (Bottom) Cartoon diagram representing folding (red arrows) and unfolding (blue arrows) of apo- and Zn2+- WT and cpF114 azurins. The width of the arrows represents the rate of the corresponding processes.

Local Xenon-Protein Interaction Produces Global Conformational Change and Allosteric Inhibition in Lysozyme

Noble gases have well-established biological effects, yet their molecular mechanisms remain poorly understood. Here, we investigated, both experimentally and computationally, the molecular modes of xenon (Xe) action in bacteriophage T4 lysozyme (T4L). By combining indirect gassing methods with a colorimetric lysozyme activity assay, a reversible, Xe-specific (20 ± 3)% inhibition effect was observed. Accelerated molecular dynamic simulations revealed that Xe exerts allosteric inhibition on the protein by expanding a C-terminal hydrophobic cavity. Xe-induced cavity expansion results in global conformational changes, with long-range transduction distorting the active site where peptidoglycan binds. Interestingly, the peptide substrate binding site that enables lysozyme specificity does not change conformation. Two T4L mutants designed to reshape the C-terminal Xe cavity established a correlation between cavity expansion and enzyme inhibition. This work also highlights the use of Xe flooding simulations to identify new cryptic binding pockets. These results enrich our understanding of Xe-protein interactions at the molecular level and inspire further biochemical investigations with noble gases.

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.

Interpretation of Single-Molecule Force Experiments on Proteins Using Normal Mode Analysis

Single-molecule force spectroscopy experiments allow protein folding and unfolding to be explored using mechanical force. Probably the most informative technique for interpreting the results of these experiments at the structural level makes use of steered molecular dynamics (MD) simulations, which can explicitly model the protein under load. Unfortunately, this technique is computationally expensive for many of the most interesting biological molecules. Here, we find that normal mode analysis (NMA), a significantly cheaper technique from a computational perspective, allows at least some of the insights provided by MD simulation to be gathered. We apply this technique to three non-homologous proteins that were previously studied by force spectroscopy: T4 lysozyme (T4L), Hsp70 and the glucocorticoid receptor domain (GCR). The NMA results for T4L and Hsp70 are compared with steered MD simulations conducted previously, and we find that we can recover the main results. For the GCR, which did not undergo MD simulation, our approach identifies substructures that correlate with experimentally identified unfolding intermediates. Overall, we find that NMA can make a valuable addition to the analysis toolkit for the structural analysis of single-molecule force experiments on proteins.

Resolving dynamics and function of transient states in single enzyme molecules

We use a hybrid fluorescence spectroscopic toolkit to monitor T4 Lysozyme (T4L) in action by unraveling the kinetic and dynamic interplay of the conformational states. In particular, by combining single-molecule and ensemble multiparameter fluorescence detection, EPR spectroscopy, mutagenesis, and FRET-positioning and screening, and other biochemical and biophysical tools, we characterize three short-lived conformational states over the ns-ms timescale. The use of 33 FRET-derived distance sets, to screen available T4L structures, reveal that T4L in solution mainly adopts the known open and closed states in exchange at 4 µs. A newly found minor state, undisclosed by, at present, more than 500 crystal structures of T4L and sampled at 230 µs, may be actively involved in the product release step in catalysis. The presented fluorescence spectroscopic toolkit will likely accelerate the development of dynamic structural biology by identifying transient conformational states that are highly abundant in biology and critical in enzymatic reactions.

T4 Lysozyme (T4L) is a model protein whose structure is extensively studied. Here the authors combine single-molecule and ensemble FRET measurements, FRET-positioning and screening and EPR spectroscopy to study the structural dynamics of T4L and describe its conformational landscape during the catalytic cycle by an extended Michaelis–Menten mechanism and identify an excited conformational state of the enzyme.

A Double-Resonance CEST Experiment To Study Multistate Protein Conformational Exchange: An Application to Protein Folding

Despite the importance of protein dynamics to function, studying exchange between multiple conformational states remains a challenge because sparsely populated states are invisible to conventional techniques. CEST NMR experiments can detect minor states with lifetimes between 5 and 200 ms populated to a level of just ∼1%. However, CEST often cannot provide the exchange mechanism for processes involving three or more states, leaving the role of the detected minor states unknown. Here a double-resonance CEST experiment to determine the kinetics of multistate exchange is presented. The approach that involves irradiating resonances from two minor states simultaneously is used to study the exchange of T4 lysozyme (T4L) between the dominant native state and two minor states, the unfolded state and a second minor state (B), each populated to only ∼4%. Regular CEST does not provide the folding mechanism, but double-resonance CEST clearly shows that T4L can fold directly without going through B.

Investigating the Effect of Chain Connectivity on the Folding of a Beta-Sheet Protein On and Off the Ribosome

Determining the relationship between protein folding pathways on and off the ribosome remains an important area of investigation in biology. Studies on isolated domains have shown that alteration of the separation of residues in a polypeptide chain, while maintaining their spatial contacts, may affect protein stability and folding pathway. Due to the vectorial emergence of the polypeptide chain from the ribosome, chain connectivity may have an important influence upon cotranslational folding. Using MATH, an all β-sandwich domain, we investigate whether the connectivity of residues and secondary structure elements is a key determinant of when cotranslational folding can occur on the ribosome. From Φ-value analysis, we show that the most structured region of the transition state for folding in MATH includes the N and C terminal strands, which are located adjacent to each other in the structure. However, arrest peptide force-profile assays show that wild-type MATH is able to fold cotranslationally, while some C-terminal residues remain sequestered in the ribosome, even when destabilized by 2–3 kcal mol⁻¹. We show that, while this pattern of Φ-values is retained in two circular permutants in our studies of the isolated domains, one of these permutants can fold only when fully emerged from the ribosome. We propose that in the case of MATH, onset of cotranslational folding is determined by the ability to form a sufficiently stable folding nucleus involving both β-sheets, rather than by the location of the terminal strands in the ribosome tunnel.

•

Adjacent N and C terminal strands are most structured region in the transition state.

•

Two circular permutants retain the same folding pathway as wild-type MATH.

•

On the ribosome, early emergence of terminal strands does not promote earlier folding.

•

Formation of both β-sheets is energetically critical for folding on the ribosome.

•

Folding pathway minimizes formation of partly structured states prone to mis-folding.

Modulating long-range energetics via helix stabilization: A case study using T4 lysozyme

Cooperative protein folding requires distant regions of a protein to interact and provide mutual stabilization. The mechanism of this long-distance coupling remains poorly understood. Here, we use T4 lysozyme (T4L*) as a model to investigate long-range communications across two subdomains of a globular protein. T4L* is composed of two structurally distinct subdomains, although it behaves in a two-state manner at equilibrium. The subdomains of T4L* are connected via two topological connections: the N-terminal helix that is structurally part of the C-terminal subdomain (the A-helix) and a long helix that spans both subdomains (the C-helix). To understand the role that the C-helix plays in cooperative folding, we analyzed a circularly permuted version of T4L* (CP13*), whose subdomains are connected only by the C-helix. We demonstrate that when isolated as individual fragments, both subdomains of CP13* can fold autonomously into marginally stable conformations. The energetics of the N-terminal subdomain depend on the formation of a salt bridge known to be important for stability in the full-length protein. We show that the energetic contribution of the salt bridge to the stability of the N-terminal fragment increases when the C-helix is stabilized, such as occurs upon folding of the C-terminal subdomain. These results suggest a model where long-range energetic coupling is mediated by helix stabilization and not specific tertiary interactions.

Effect of circular permutations on transient partial unfolding in proteins

Under native conditions, proteins can undergo transient partial unfolding, which may cause proteins to misfold or aggregate. A change in sequence connectivity by circular permutation may affect the energetics of transient partial unfolding in proteins without altering the three-dimensional structures. Using Escherichia coli dihydrofolate reductase (DHFR) as a model system, we investigated how circular permutation affects transient partial unfolding in proteins. We constructed three circular permutants, CP18, CP37, and CP87, with the new N-termini at residue 18, 37, and 87, respectively, and probed transient partial unfolding by native-state proteolysis. The new termini in CP18, CP37, and CP87 are within, near, and distal to the Met20 loop, which is known to be dynamic and also part of the region that undergoes transient unfolding in wild-type DHFR. The stabilities of both native and partially unfolded forms of CP18 are similar to those of wild-type DHFR, suggesting that the influence of introducing new termini in a dynamic region to the protein is minimal. CP37 has a significantly more accessible partially unfolded form than wild-type DHFR, demonstrating that introducing new termini near a dynamic region may promote transient partial unfolding. CP87 has significantly destabilized native and partially unfolded forms, confirming that modification of the folded region in a partially unfolded form destabilizes the partially unfolded form similar to the native form. Our findings provide valuable guidelines to control transient partial unfolding in designing circular permutants in proteins.

Network representation of conformational transitions between hidden intermediates of Rd-apocytochrome b562

The folding kinetics of Rd-apocytochrome b562 is two-state, but native-state hydrogen exchange experiments show that there are discrete partially unfolded (PUF) structures in equilibrium with the native state. These PUF structures are called hidden intermediates because they are not detected in kinetic experiments and they exist after the rate-limiting step. Structures of the mimics of hidden intermediates of Rd-apocytochrome b562 are resolved by NMR. Based upon their relative stability and structural features, the folding mechanism was proposed to follow a specific pathway (unfolded → rate-limiting transition state → PUF1 → PUF2 → native). Investigating the roles of equilibrium PUF structures in folding kinetics and their interrelationship not only deepens our understanding of the details of folding mechanism but also provides guides in protein design and prevention of misfolding. We performed molecular dynamics simulations starting from a hidden intermediate and the native state of Rd-apocytochrome b562 in explicit solvent, for a total of 37.18 μs mainly with Anton. We validated our simulations by detailed comparison with experimental data and other computations. We have verified that we sampled the post rate-limiting transition state region only. Markov state model was used to analyze the simulation results. We replace the specific pathway model with a network model. Transition-path theory was employed to calculate the net effective flux from the most unfolded state towards the most folded state in the network. The proposed sequential folding pathway via PUF1 then more stable, more native-like PUF2 is one of the routes in our network, but it is not dominant. The dominant path visits PUF2 without going through PUF1. There is also a route from PUF1 directly to the most folded state in the network without visiting PUF2. Our results indicate that the PUF states are not necessarily sequential in the folding. The major routes predicted in our network are testable by future experiments such as single molecule experiment.

Protein folding and unfolding under force

The recent revolution in optics and instrumentation has enabled the study of protein folding using extremely low mechanical forces as the denaturant. This exciting development has led to the observation of the protein folding process at single molecule resolution and its response to mechanical force. Here, we describe the principles and experimental details of force spectroscopy on proteins, with a focus on the optical tweezers instrument. Several recent results will be discussed to highlight the importance of this technique in addressing a variety of questions in the protein folding field.

The folding unit of phosphofructokinase-2 as defined by the biophysical properties of a monomeric mutant

Escherichia coli phosphofructokinase-2 (Pfk-2) is an obligate homodimer that follows a highly cooperative three-state folding mechanism N2 ↔ 2I ↔ 2U. The strong coupling between dissociation and unfolding is a consequence of the structural features of its interface: a bimolecular domain formed by intertwining of the small domain of each subunit into a flattened β-barrel. Although isolated monomers of E. coli Pfk-2 have been observed by modification of the environment (changes in temperature, addition of chaotropic agents), no isolated subunits in native conditions have been obtained. Based on in silico estimations of the change in free energy and the local energetic frustration upon binding, we engineered a single-point mutant to destabilize the interface of Pfk-2. This mutant, L93A, is an inactive monomer at protein concentrations below 30 μM, as determined by analytical ultracentrifugation, dynamic light scattering, size exclusion chromatography, small-angle x-ray scattering, and enzyme kinetics. Active dimer formation can be induced by increasing the protein concentration and by addition of its substrate fructose-6-phosphate. Chemical and thermal unfolding of the L93A monomer followed by circular dichroism and dynamic light scattering suggest that it unfolds noncooperatively and that the isolated subunit is partially unstructured and marginally stable. The detailed structural features of the L93A monomer and the F6P-induced dimer were ascertained by high-resolution hydrogen/deuterium exchange mass spectrometry. Our results show that the isolated subunit has overall higher solvent accessibility than the native dimer, with the exception of residues 240-309. These residues correspond to most of the β-meander module and show the same extent of deuterium uptake as the native dimer. Our results support the idea that the hydrophobic core of the isolated monomer of Pfk-2 is solvent-penetrated in native conditions and that the β-meander module is not affected by monomerizing mutations.

Probing the folded state and mechanical unfolding pathways of T4 lysozyme using all-atom and coarse-grained molecular simulation

The Bacteriophage T4 Lysozyme (T4L) is a prototype modular protein comprised of an N-terminal and a C-domain domain, which was extensively studied to understand the folding/unfolding mechanism of modular proteins. To offer detailed structural and dynamic insights to the folded-state stability and the mechanical unfolding behaviors of T4L, we have performed extensive equilibrium and steered molecular dynamics simulations of both the wild-type (WT) and a circular permutation (CP) variant of T4L using all-atom and coarse-grained force fields. Our all-atom and coarse-grained simulations of the folded state have consistently found greater stability of the C-domain than the N-domain in isolation, which is in agreement with past thermostatic studies of T4L. While the all-atom simulation cannot fully explain the mechanical unfolding behaviors of the WT and the CP variant observed in an optical tweezers study, the coarse-grained simulations based on the Go model or a modified elastic network model (mENM) are in qualitative agreement with the experimental finding of greater unfolding cooperativity in the WT than the CP variant. Interestingly, the two coarse-grained models predict different structural mechanisms for the observed change in cooperativity between the WT and the CP variant--while the Go model predicts minor modification of the unfolding pathways by circular permutation (i.e., preserving the general order that the N-domain unfolds before the C-domain), the mENM predicts a dramatic change in unfolding pathways (e.g., different order of N/C-domain unfolding in the WT and the CP variant). Based on our simulations, we have analyzed the limitations of and the key differences between these models and offered testable predictions for future experiments to resolve the structural mechanism for cooperative folding/unfolding of T4L.

The response of Greek key proteins to changes in connectivity depends on the nature of their secondary structure

What governs the balance between connectivity and topology in regulating the mechanism of protein folding? We use circular permutation to vary the order of the helices in the all-α Greek key protein FADD (Fas-associated death domain) to investigate this question. Unlike all-β Greek key proteins, where changes in the order of secondary structure cause a shift in the folding nucleus, the position of the nucleus in FADD is unchanged, even when permutation reduces the complexity significantly. We suggest that this is because local helical contacts are so dominant that permutation has little effect on the entropic cost of forming the folding nucleus whereas, in all-β Greek key proteins, all interactions in the nucleus are long range. Thus, the type of secondary structure modulates the sensitivity of proteins to changes in connectivity.

Role of cavities and hydration in the pressure unfolding of T4 lysozyme

It is well known that high hydrostatic pressures can induce the unfolding of proteins. The physical underpinnings of this phenomenon have been investigated extensively but remain controversial. Changes in solvation energetics have been commonly proposed as a driving force for pressure-induced unfolding. Recently, the elimination of void volumes in the native folded state has been argued to be the principal determinant. Here we use the cavity-containing L99A mutant of T4 lysozyme to examine the pressure-induced destabilization of this multidomain protein by using solution NMR spectroscopy. The cavity-containing C-terminal domain completely unfolds at moderate pressures, whereas the N-terminal domain remains largely structured to pressures as high as 2.5 kbar. The sensitivity to pressure is suppressed by the binding of benzene to the hydrophobic cavity. These results contrast to the pseudo-WT protein, which has a residual cavity volume very similar to that of the L99A-benzene complex but shows extensive subglobal reorganizations with pressure. Encapsulation of the L99A mutant in the aqueous nanoscale core of a reverse micelle is used to examine the hydration of the hydrophobic cavity. The confined space effect of encapsulation suppresses the pressure-induced unfolding transition and allows observation of the filling of the cavity with water at elevated pressures. This indicates that hydration of the hydrophobic cavity is more energetically unfavorable than global unfolding. Overall, these observations point to a range of cooperativity and energetics within the T4 lysozyme molecule and illuminate the fact that small changes in physical parameters can significantly alter the pressure sensitivity of proteins.

Modulation of a protein free-energy landscape by circular permutation

Circular permutations usually retain the native structure and function of a protein while inevitably perturbing its folding dynamics. By using simulations with a structure-based model and a rigorous methodology to determine free-energy surfaces from trajectories, we evaluate the effect of a circular permutation on the free-energy landscape of the protein T4 lysozyme. We observe changes which, although subtle, largely affect the cooperativity between the two subdomains. Such a change in cooperativity has been previously experimentally observed and recently also characterized using single molecule optical tweezers and the Crooks relation. The free-energy landscapes show that both the wild type and circular permutant have an on-pathway intermediate, previously experimentally characterized, in which one of the subdomains is completely formed. The landscapes, however, differ in the position of the rate-limiting step for folding, which occurs before the intermediate in the wild type and after in the circular permutant. This shift of transition state explains the observed change in the cooperativity. The underlying free-energy landscape thus provides a microscopic description of the folding dynamics and the connection between circular permutation and the loss of cooperativity experimentally observed.

Remodeling of the folding free energy landscape of staphylococcal nuclease by cavity-creating mutations

The folding of staphylococcal nuclease (SNase) is known to proceed via a major intermediate in which the central OB subdomain is folded and the C-terminal helical subdomain is disordered. To identify the structural and energetic determinants of this folding free energy landscape, we have examined in detail, using high-pressure NMR, the consequences of cavity creating mutations in each of the two subdomains of an ultrastable SNase, Δ+PHS. The stabilizing mutations of Δ+PHS enhanced the population of the major folding intermediate. Cavity creation in two different regions of the Δ+PHS reference protein, despite equivalent effects on global stability, had very distinct consequences on the complexity of the folding free energy landscape. The L125A substitution in the C-terminal helix of Δ+PHS slightly suppressed the major intermediate and promoted an additional excited state involving disorder in the N-terminus, but otherwise decreased landscape heterogeneity with respect to the Δ+PHS background protein. The I92A substitution, located in the hydrophobic OB-fold core, had a much more profound effect, resulting in a significant increase in the number of intermediate states and implicating the entire protein structure. Denaturant (GuHCl) had very subtle and specific effects on the landscape, suppressing some states and favoring others, depending upon the mutational context. These results demonstrate that disrupting interactions in a region of the protein with highly cooperative, unfrustrated folding has very profound effects on the roughness of the folding landscape, whereas the effects are less pronounced for an energetically equivalent substitution in an already frustrated region.

Engineering a model protein cavity to catalyze the Kemp elimination

Synthetic cavitands and protein cavities have been widely studied as models for ligand recognition. Here we investigate the Met102 → His substitution in the artificial L99A cavity in T4 lysozyme as a Kemp eliminase. The resulting enzyme had k(cat)/K(M) = 0.43 M(-1) s(-1) and a (k(cat)/K(M))/k(uncat) = 10(7) at pH 5.0. The crystal structure of this enzyme was determined at 1.30 Å, as were the structures of four complexes of substrate and product analogs. The absence of ordered waters or hydrogen bonding interactions, and the presence of a common catalytic base (His102) in an otherwise hydrophobic, buried cavity, facilitated detailed analysis of the reaction mechanism and its optimization. Subsequent substitutions increased eliminase activity by an additional four-fold. As activity-enhancing substitutions were engineered into the cavity, protein stability decreased, consistent with the stability-function trade-off hypothesis. This and related model cavities may provide templates for studying protein design principles in radically simplified environments.

The ribosome modulates nascent protein folding

Proteins are synthesized by the ribosome and generally must fold to become functionally active. Although it is commonly assumed that the ribosome affects the folding process, this idea has been extremely difficult to demonstrate. We have developed an experimental system to investigate the folding of single ribosome-bound stalled nascent polypeptides with optical tweezers. In T4 lysozyme, synthesized in a reconstituted in vitro translation system, the ribosome slows the formation of stable tertiary interactions and the attainment of the native state relative to the free protein. Incomplete T4 lysozyme polypeptides misfold and aggregate when free in solution, but they remain folding-competent near the ribosomal surface. Altogether, our results suggest that the ribosome not only decodes the genetic information and synthesizes polypeptides, but also promotes efficient de novo attainment of the native state.

High-pressure EPR reveals conformational equilibria and volumetric properties of spin-labeled proteins

Identifying equilibrium conformational exchange and characterizing conformational substates is essential for elucidating mechanisms of function in proteins. Site-directed spin labeling has previously been employed to detect conformational changes triggered by some event, but verifying conformational exchange at equilibrium is more challenging. Conformational exchange (microsecond-millisecond) is slow on the EPR time scale, and this proves to be an advantage in directly revealing the presence of multiple substates as distinguishable components in the EPR spectrum, allowing the direct determination of equilibrium constants and free energy differences. However, rotameric exchange of the spin label side chain can also give rise to multiple components in the EPR spectrum. Using spin-labeled mutants of T4 lysozyme, it is shown that high-pressure EPR can be used to: (i) demonstrate equilibrium between spectrally resolved states, (ii) aid in distinguishing conformational from rotameric exchange as the origin of the resolved states, and (iii) determine the relative partial molar volume (ΔV°) and isothermal compressibility (Δβ(T)) of conformational substates in two-component equilibria from the pressure dependence of the equilibrium constant. These volumetric properties provide insight into the structure of the substates. Finally, the pressure dependence of internal side-chain motion is interpreted in terms of volume fluctuations on the nanosecond time scale, the magnitude of which may reflect local backbone flexibility.

Multifrequency electron spin resonance study of the dynamics of spin labeled T4 lysozyme

An extensive set of electron spin resonance spectra was obtained over a wide range of frequencies (9, 95, 170, and 240 GHz) and temperatures (2 to 32 degrees C) to explore the dynamic modes of nitroxide-labeled T4 lysozyme in solution. A commonly used nitroxide side chain (R1), or a methylated analogue with hindered internal motion (R2), was substituted for the native side chain at solvent-exposed helical sites, 72 or 131. The spectra at all four frequencies were simultaneously fit with the slowly relaxing local structure (SRLS) model. Good fits were achieved at all the temperatures. Two principle dynamic modes are included in the SRLS model, the global tumbling of the protein and the internal motion consisting of backbone fluctuations and side chain isomerizations. Three distinct spectral components were required for R1 and two for R2 to account for the spectra at all temperatures. One is a highly ordered and slow motional component, which is observed in the spectra of both R1 and R2; it may correspond to conformers stabilized by interaction with the protein surface. The fraction of this component decreases with increasing temperature and is more populated in the R2 spectra, possibly arising from stronger interaction of the nitroxide ring with the protein surface due to the additional methyl group. The other two components of R1 and the second component of R2 are characterized by fast anisotropic diffusion and relatively low ordering, most likely corresponding to conformers having little or no interactions with nearby residues. Ficoll of different concentrations was added to increase the solution viscosity, thereby slowing down the global tumbling of the protein. A significant effect of Ficoll on the internal motion of an immobilized component was apparent in R2 but not in R1. The ability of such multifrequency studies to separate the effects of faster internal modes of motion from slower overall motions is clearly demonstrated, and its utility in future studies is considered.

The folding cooperativity of a protein is controlled by its chain topology

The three-dimensional structures of proteins often show a modular architecture comprised of discrete structural regions or domains. Cooperative communication between these regions is important for catalysis, regulation and efficient folding; lack of coupling has been implicated in the formation of fibrils and other misfolding pathologies. How different structural regions of a protein communicate and contribute to a protein's overall energetics and folding, however, is still poorly understood. Here we use a single-molecule optical tweezers approach to induce the selective unfolding of particular regions of T4 lysozyme and monitor the effect on other regions not directly acted on by force. We investigate how the topological organization of a protein (the order of structural elements along the sequence) affects the coupling and folding cooperativity between its domains. To probe the status of the regions not directly subjected to force, we determine the free energy changes during mechanical unfolding using Crooks' fluctuation theorem. We pull on topological variants (circular permutants) and find that the topological organization of the polypeptide chain critically determines the folding cooperativity between domains and thus what parts of the folding/unfolding landscape are explored. We speculate that proteins may have evolved to select certain topologies that increase coupling between regions to avoid areas of the landscape that lead to kinetic trapping and misfolding.

Engineering an artificial zymogen by alternate frame protein folding

Alternate frame folding (AFF) is a novel mechanism by which allostery can be introduced into a protein where none may have existed previously. We employ this technology to convert the cytotoxic ribonuclease barnase into an artificial zymogen that is activated by HIV-1 protease. The AFF modification entails partial duplication of the polypeptide chain and mutation of a key catalytic residue in one of the duplicated segments. The resulting molecule can fold in one of two "frames" to yield the wild-type structure or a circularly permuted form in which the positions of the N- and C-termini are exchanged with a surface loop. It cannot take on both structures simultaneously because each competes for a shared amino acid sequence. An HIV-1 protease recognition sequence is inserted into one of the surface loops in the nonpermuted frame, and cleavage induces a shift from the nonpermuted fold to the permuted fold. Using the AFF mechanism, we were able to suppress k(cat)/K(M) by 250-fold in the proenzyme relative to wild-type barnase. HIV-1 protease cleavage subsequently increases k(cat)/K(M) by 130-fold. AFF is significant because it is general and can in principle be used to control activity of many enzymes, including those whose functions are not regulated by any existing mechanism.

Subdomain competition, cooperativity, and topological frustration in the folding of CheY

The folding of multidomain proteins often proceeds in a hierarchical fashion with individual domains folding independent of one another. A large single-domain protein, however, can consist of multiple modules whose folding may be autonomous or interdependent in ways that are unclear. We used coarse-grained simulations to explore the folding landscape of the two-subdomain bacterial response regulator CheY. Thermodynamic and kinetic characterization shows the landscape to be highly analogous to the four-state landscape reported for another two-subdomain protein, T4 lysozyme. An on-pathway intermediate structured in the more stable nucleating subdomain was observed, as were transient states frustrated in off-pathway contacts prematurely structured in the weaker subdomain. Local unfolding, or backtracking, was observed in the frustrated state before the native conformation could be reached. Nonproductive frustration was attributable to competition for van der Waals contacts between the two subdomains. In an accompanying article, stopped-flow kinetic measurements support an off-pathway burst-phase intermediate, seemingly consistent with our prediction of early frustration in the folding landscape of CheY. Comparison of the folding mechanisms for CheY, T4 lysozyme, and interleukin-1 beta leads us to postulate that subdomain competition is a general feature of large single-domain proteins with multiple folding modules.

Atomic force microscopy reveals parallel mechanical unfolding pathways of T4 lysozyme: evidence for a kinetic partitioning mechanism

Kinetic partitioning is predicted to be a general mechanism for proteins to fold into their well defined native three-dimensional structure from unfolded states following multiple folding pathways. However, experimental evidence supporting this mechanism is still limited. By using single-molecule atomic force microscopy, here we report experimental evidence supporting the kinetic partitioning mechanism for mechanical unfolding of T4 lysozyme, a small protein composed of two subdomains. We observed that on stretching from its N and C termini, T4 lysozyme unfolds by multiple distinct unfolding pathways: the majority of T4 lysozymes unfold in an all-or-none fashion by overcoming a dominant unfolding kinetic barrier; and a small fraction of T4 lysozymes unfold in three-state fashion involving unfolding intermediate states. The three-state unfolding pathways do not follow well defined routes, instead they display variability and diversity in individual unfolding pathways. The unfolding intermediate states are local energy minima along the mechanical unfolding pathways and are likely to result from the residual structures present in the two subdomains after crossing the main unfolding barrier. These results provide direct evidence for the kinetic partitioning of the mechanical unfolding pathways of T4 lysozyme, and the complex unfolding behaviors reflect the stochastic nature of kinetic barrier rupture in mechanical unfolding processes. Our results demonstrate that single-molecule atomic force microscopy is an ideal tool to investigate the folding/unfolding dynamics of complex multimodule proteins that are otherwise difficult to study using traditional methods.

Exploring subdomain cooperativity in T4 lysozyme II: uncovering the C-terminal subdomain as a hidden intermediate in the kinetic folding pathway

Intermediates along a protein's folding pathway can play an important role in its biology. Previous kinetics studies have revealed an early folding intermediate for T4 lysozyme, a small, well-characterized protein composed of an N-terminal and a C-terminal subdomain. Pulse-labeling hydrogen exchange studies suggest that residues from both subdomains contribute to the structure of this intermediate. On the other hand, equilibrium native state hydrogen experiments have revealed a high-energy, partially unfolded form of the protein that has an unstructured N-terminal subdomain and a structured C-terminal subdomain. To resolve this discrepancy between kinetics and equilibrium data, we performed detailed kinetics analyses of the folding and unfolding pathways of T4 lysozyme, as well as several point mutants and large-scale variants. The data support the argument for the presence of two distinct intermediates, one present on each side of the rate-limiting transition state barrier. The effects of circular permutation and site-specific mutations in the wild-type and circular permutant background, as well as a fragment containing just the C-terminal subdomain, support a model for the unfolding intermediate with an unfolded N-terminal and a folded C-terminal subdomain. Our results suggest that the partially unfolded form identified by native state hydrogen exchange resides on the folded side of the rate-limiting transition state and is, therefore, under most conditions, a "hidden" intermediate.