Relationship of Leffler (Bronsted) alpha values and protein folding Phi values to position of transition-state structures on reaction coordinates

The Mechanism of Folding of Human Frataxin in Comparison to the Yeast Homologue - Broad Energy Barriers and the General Properties of the Transition State

The funneled energy landscape theory suggests that the folding pathway of homologous proteins should converge at the late stages of folding. In this respect, proteins displaying a broad energy landscape for folding are particularly instructive, allowing inferring both the early, intermediate and late stages of folding. In this paper we explore the folding mechanisms of human frataxin, an essential mitochondrial protein linked to the neurodegenerative disorder Friedreich's ataxia. Building upon previous studies on the yeast homologue, the folding pathway of human frataxin is thoroughly examined, revealing a mechanism implying the presence of a broad energy barrier, reminiscent of the yeast counterpart. Through an extensive site-directed mutagenesis, we employed a Φ -value analysis to map native-like contacts in the folding transition state. The presence of a broad energy barrier facilitated the exploration of such contacts in both early and late folding events. We compared results from yeast and human frataxin providing insights into the impact of native topology on the folding mechanism and elucidating the properties of the underlying free energy landscape. The findings are discussed in the context of the funneled energy landscape theory of protein folding.

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

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

Understanding the Binding Transition State After the Conformational Selection Step: The Second Half of the Molecular Recognition Process Between NS1 of the 1918 Influenza Virus and Host p85β

Biomolecular recognition often involves conformational changes as a prerequisite for binding (i.e., conformational selection) or concurrently with binding (i.e., induced-fit). Recent advances in structural and kinetic approaches have enabled the detailed characterization of protein motions at atomic resolution. However, to fully understand the role of the conformational dynamics in molecular recognition, studies on the binding transition state are needed. Here, we investigate the binding transition state between nonstructural protein 1 (NS1) of the pandemic 1918 influenza A virus and the human p85β subunit of PI3K. 1918 NS1 binds to p85β via conformational selection. We present the free-energy mapping of the transition and bound states of the 1918 NS1:p85β interaction using linear free energy relationship and ϕ-value analyses. We find that the binding transition state of 1918 NS1 and p85β is structurally similar to the bound state with well-defined binding orientation and hydrophobic interactions. Our finding provides a detailed view of how protein motion contributes to the development of intermolecular interactions along the binding reaction coordinate.

Converging experimental and computational views of the knotting mechanism of a small knotted protein

MJ0366 from Methanocaldococcus jannaschii is the smallest topologically knotted protein known to date. 92 residues in length, MJ0366 ties a trefoil (31) knot by threading its C-terminal helix through a buttonhole formed by the remainder of the secondary structure elements. By generating a library of point mutations at positions pertinent to the knot formation, we systematically evaluated the contributions of individual residues to the folding stability and kinetics of MJ0366. The experimental Φ-values were used as restraints to computationally generate an ensemble of conformations that correspond to the transition state of MJ0366, which revealed several nonnative contacts. The importance of these nonnative contacts in stabilizing the transition state of MJ0366 was confirmed by a second round of mutagenesis, which also established the pivotal role of F15 in stapling the network of hydrophobic interactions around the threading C-terminal helix. Our converging experimental and computational results show that, despite the small size, the transition state of MJ0366 is formed at a very late stage of the folding reaction coordinate, following a polarized pathway. Eventually, the formation of extensive native contacts, as well as a number of nonnative ones, leads to the threading of the C-terminal helix that defines the topological knot.

Mapping the energy landscape of protein-ligand binding via linear free energy relationships determined by protein NMR relaxation dispersion

Biochemical signaling is mediated by complexes between macromolecular receptors and their ligands, with the duration of the signal being directly related to the lifetime of the ligand-receptor complex. In the field of drug design, the recognition that drug efficacy in vivo depends on the lifetime of the drug-protein complex has spawned the concept of designing drugs with particular binding kinetics. To advance this field it is critical to investigate how the molecular details of designed ligands might affect the binding kinetics, as well as the equilibrium binding constant. Here we use protein NMR relaxation dispersion to determine linear free energy relationships involving the on- and off-rates and the affinity for a series of congeneric ligands targeting the carbohydrate recognition domain of galectin-3. Using this approach we determine the energy landscape and the position of the transition state along the reaction coordinate of protein-ligand binding. The results show that ligands exhibiting reduced off-rates achieve this by primarily stabilizing the bound state, but do not affect the transition state to any greater extent. The transition state forms early, that is, it is located significantly closer to the free state than to the bound state, suggesting a critical role of desolvation. Furthermore, the data suggest that different subclasses of ligands show different behavior with respect to these characteristics.

Demonstration of Binding Induced Structural Plasticity in a SH2 Domain

SH2 domains are common protein interaction domains able to recognize short aminoacidic sequences presenting a phosphorylated tyrosine (pY). In spite of their fundamental importance for cell physiology there is a lack of information about the mechanism by which these domains recognize and bind their natural ligands. The N-terminal SH2 (N-SH2) domain of PI3K mediates the interaction with different scaffolding proteins and is known to recognize a specific pY-X-X-M consensus sequence. These interactions are at the cross roads of different molecular pathways and play a key role for cell development and division. By combining mutagenesis, chemical kinetics and NMR, here we provide a complete characterization of the interaction between N-SH2 and a peptide mimicking the scaffolding protein Gab2. Our results highlight that N-SH2 is characterized by a remarkable structural plasticity, with the binding reaction being mediated by a diffused structural region and not solely by the residues located in the binding pocket. Furthermore, the analysis of kinetic data allow us to pinpoint an allosteric network involving residues far from the binding pocket involved in specificity. Results are discussed on the light of previous works on the binding properties of SH2 domains.

Understanding the Mechanism of Recognition of Gab2 by the N-SH2 Domain of SHP2

Gab2 is a scaffold protein with a crucial role in colocalizing signaling proteins and it is involved in the regulation of several important molecular pathways. SHP2 is a protein phosphatase that binds, through its two SH2 domains, specific consensus sequences presenting a phosphorylated tyrosine located on the disordered tail of Gab2. To shed light on the details of such a fundamental interaction for the physiology of the cell, we present a complete mutational analysis of the kinetics of binding between the N-SH2 domain of SHP2 and a peptide mimicking a specific region of Gab2. By analyzing kinetic data, we determined structural features of the transition state of the N-SH2 domain binding to Gab2, highlighting a remarkable cooperativity of the binding reaction. Furthermore, comparison of these data with ones previously obtained for another SH2 domain suggests the presence of underlying general features characterizing the binding process of SH2 domains. Data are discussed under the light of previous works on SH2 domains.

Rapid Quantification of Protein-Ligand Binding via 19F NMR Lineshape Analysis

Fluorine incorporation is ideally suited to many NMR techniques, and incorporation of fluorine into proteins and fragment libraries for drug discovery has become increasingly common. Here, we use one-dimensional 19F NMR lineshape analysis to quantify the kinetics and equilibrium thermodynamics for the binding of a fluorine-labeled Src homology 3 (SH3) protein domain to four proline-rich peptides. SH3 domains are one of the largest and most well-characterized families of protein recognition domains and have a multitude of functions in eukaryotic cell signaling. First, we showe that fluorine incorporation into SH3 causes only minor structural changes to both the free and bound states using amide proton temperature coefficients. We then compare the results from lineshape analysis of one-dimensional 19F spectra to those from two-dimensional 1H-15N heteronuclear single quantum coherence spectra. Their agreement demonstrates that one-dimensional 19F lineshape analysis is a robust, low-cost, and fast alternative to traditional heteronuclear single quantum coherence-based experiments. The data show that binding is diffusion limited and indicate that the transition state is highly similar to the free state. We also measured binding as a function of temperature. At equilibrium, binding is enthalpically driven and arises from a highly positive activation enthalpy for association with small entropic contributions. Our results agree with those from studies using different techniques, providing additional evidence for the utility of 19F NMR lineshape analysis, and we anticipate that this analysis will be an effective tool for rapidly characterizing the energetics of protein interactions.

An Early Association between the α-Helix of the TEAD Binding Domain of YAP and TEAD Drives the Formation of the YAP:TEAD Complex

The Hippo pathway is an evolutionarily conserved signaling pathway that is involved in the control of organ size and development. The TEAD transcription factors are the most downstream elements of the Hippo pathway, and their transcriptional activity is regulated via the interaction with different co-regulators such as YAP. The structure of the YAP:TEAD complex shows that YAP binds to TEAD via two distinct secondary structure elements, an α-helix and an Ω-loop, and site-directed mutagenesis experiments revealed that the Ω-loop is the "hot spot" of this interaction. While much is known about how YAP and TEAD interact with each other, little is known about the mechanism leading to the formation of a complex between these two proteins. Here we combine site-directed mutagenesis with pre-steady-state kinetic measurements to show that the association between these proteins follows an apparent one-step binding mechanism. Furthermore, linear free energy relationships and a Φ analysis suggest that binding-induced folding of the YAP α-helix to TEAD occurs independently of and before formation of the Ω-loop interface. Thus, the binding-induced folding of YAP appears not to conform to the concomitant formation of tertiary structure (nucleation-condensation) usually observed for coupled binding and folding reactions. Our findings demonstrate how a mechanism reminiscent of the classical framework (diffusion-collision) mechanism of protein folding may operate in disorder-to-order transitions involving intrinsically disordered proteins.

Templated folding of intrinsically disordered proteins

Much of our current knowledge of biological chemistry is founded in the structure-function relationship, whereby sequence determines structure that determines function. Thus, the discovery that a large fraction of the proteome is intrinsically disordered, while being functional, has revolutionized our understanding of proteins and raised new and interesting questions. Many intrinsically disordered proteins (IDPs) have been determined to undergo a disorder-to-order transition when recognizing their physiological partners, suggesting that their mechanisms of folding are intrinsically different from those observed in globular proteins. However, IDPs also follow some of the classic paradigms established for globular proteins, pointing to important similarities in their behavior. In this review, we compare and contrast the folding mechanisms of globular proteins with the emerging features of binding-induced folding of intrinsically disordered proteins. Specifically, whereas disorder-to-order transitions of intrinsically disordered proteins appear to follow rules of globular protein folding, such as the cooperative nature of the reaction, their folding pathways are remarkably more malleable, due to the heterogeneous nature of their folding nuclei, as probed by analysis of linear free-energy relationship plots. These insights have led to a new model for the disorder-to-order transition in IDPs termed "templated folding," whereby the binding partner dictates distinct structural transitions en route to product, while ensuring a cooperative folding.

Binding induced folding: Lessons from the kinetics of interaction between NTAIL and XD

Intrinsically Disordered Proteins (IDPs) are a class of protein that exert their function despite lacking a well-defined three-dimensional structure, which is sometimes achieved only upon binding to their natural ligands. This feature implies the folding of IDPs to be generally coupled with a binding event, representing an interesting challenge for kinetic studies. In this review, we recapitulate some of the most important findings of IDPs binding-induced folding mechanisms obtained by analyzing their binding kinetics. Furthermore, by focusing on the interaction between the Measles virus N_TAIL protein, a prototypical IDP, and its physiological partner, the X domain, we recapitulate the major theoretical and experimental approaches that were used to describe binding induced folding.

Structural characterization of an on-pathway intermediate and transition state in the folding of the N-terminal SH2 domain from SHP2

Src Homology 2 (SH2) domains are a class of protein domains that present a conserved three-dimensional structure and possess a crucial role in mediating protein-protein interactions. Despite their importance and abundance in the proteome, knowledge about the folding properties of SH2 domain is limited. Here we present an extensive mutational analysis (Φ value analysis) of the folding pathway of the N-SH2 domain of the Src homology region 2 domain-containing phosphatase-2 (SHP2) protein, a 104 residues domain that presents the classical SH2 domain fold (two α-helices flanking a central β-sheet composed of 3-5 antiparallel β-strands), with a fundamental role in mediating the interaction of SHP2 with its substrates and triggering key metabolic pathways in the cell. By analysing folding kinetic data we demonstrated that the folding pathway of N-SH2 presents an obligatory on-pathway intermediate that accumulates during the folding reaction. The production of 24 conservative site-directed variants allowed us to perform a Φ value analysis, by which we could fully characterize the intermediate and the transition state native-like interactions, providing a detailed quantitative analysis of the folding pathway of N-SH2. Results highlight the presence of a hydrophobic nucleus that stabilizes the intermediate, leading to a higher degree of native-like interactions in the transition state. Data are discussed and compared with previous works on SH2 domains.

How Robust Is the Mechanism of Folding-Upon-Binding for an Intrinsically Disordered Protein?

The mechanism of interaction of an intrinsically disordered protein (IDP) with its physiological partner is characterized by a disorder-to-order transition in which a recognition and a binding step take place. Even if the mechanism is quite complex, IDPs tend to bind their partner in a cooperative manner such that it is generally possible to detect experimentally only the disordered unbound state and the structured complex. The interaction between the disordered C-terminal domain of the measles virus nucleoprotein (N_TAIL) and the X domain (XD) of the viral phosphoprotein allows us to detect and quantify the two distinct steps of the overall reaction. Here, we analyze the robustness of the folding of N_TAIL upon binding to XD by measuring the effect on both the folding and binding steps of N_TAIL when the structure of XD is modified. Because it has been shown that wild-type XD is structurally heterogeneous, populating an on-pathway intermediate under native conditions, we investigated the binding to 11 different site-directed variants of N_TAIL of one particular variant of XD (I504A XD) that populates only the native state. Data reveal that the recognition and the folding steps are both affected by the structure of XD, indicating a highly malleable pathway. The experimental results are briefly discussed in the light of previous experiments on other IDPs.

The Cost of Long Catalytic Loops in Folding and Stability of the ALS-Associated Protein SOD1

A conspicuous feature of the amyotrophic lateral sclerosis (ALS)-associated protein SOD1 is that its maturation into a functional enzyme relies on local folding of two disordered loops into a catalytic subdomain. To drive the disorder-to-order transition, the protein employs a single Zn²⁺ ion. The question is then if the entropic penalty of maintaining such disordered loops in the immature apoSOD1 monomer is large enough to explain its unusually low stability, slow folding, and pathological aggregation in ALS. To find out, we determined the effects of systematically altering the SOD1-loop lengths by protein redesign. The results show that the loops destabilize the apoSOD1 monomer by ∼3 kcal/mol, rendering the protein marginally stable and accounting for its aggregation behavior. Yet the effect on the global folding kinetics remains much smaller with a transition-state destabilization of <1 kcal/mol. Notably, this 1/3 transition-state to folded-state stability ratio provides a clear-cut example of the enigmatic disagreement between the Leffler α value from loop-length alterations (typically 1/3) and the "standard" reaction coordinates based on solvent perturbations (typically >2/3). Reconciling the issue, we demonstrate that the disagreement disappears when accounting for the progressive loop shortening that occurs along the folding pathway. The approach assumes a consistent Flory loop entropy scaling factor of c = 1.48 for both equilibrium and kinetic data and has the added benefit of verifying the tertiary interactions of the folding nucleus as determined by phi-value analysis. Thus, SOD1 not only represents a case where evolution of key catalytic function has come with the drawback of a destabilized apo state but also stands out as a well-suited model system for exploring the physicochemical details of protein self-organization.

Kinetic Insights into the Binding between the nSH3 Domain of CrkII and Proline-Rich Motifs in cAbl

The interaction between CrkII and cAbl is implicated in diverse cellular processes. This interaction starts with the binding of the N-terminal Src homology 3 (nSH3) domain of CrkII to the proline-rich motifs of cAbl (PRMs^cAbl). Despite its critical importance, the detailed binding mechanism between the nSH3 domain and PRMs remains elusive. In this study, we used nuclear magnetic resonance Carr-Purcell-Meiboom-Gill relaxation dispersion experiment to study the binding kinetics between the nSH3 domain of CrkII and PRMs^cAbl. Our results highlight that the nSH3 domain binds to three PRMs^cAbl with very high on- and off-rate constants, indicating the transient nature of the binding. To further characterize the binding transition state, we conducted the Eyring and linear free energy relationship analyses using temperature-dependent kinetic data. These data indicate that the binding transition state of the nSH3 domain and PRM is accompanied by small activation enthalpy, owing to partial desolvation of the transition state. These results also highlight the similarity between the transition and free states, in terms of structure and energetics. Although the binding of the nSH3 domain and PRM displays the features consistent with a diffusion-limited process within our experimental conditions, further tests are necessary to determine if the binding is a true diffusion-limited process.

Towards a structural biology of the hydrophobic effect in protein folding

The hydrophobic effect is a major driving force in protein folding. A complete understanding of this effect requires the description of the conformational states of water and protein molecules at different temperatures. Towards this goal, we characterise the cold and hot denatured states of a protein by modelling NMR chemical shifts using restrained molecular dynamics simulations. A detailed analysis of the resulting structures reveals that water molecules in the bulk and at the protein interface form on average the same number of hydrogen bonds. Thus, even if proteins are 'large' particles (in terms of the hydrophobic effect, i.e. larger than 1 nm), because of the presence of complex surface patterns of polar and non-polar residues their behaviour can be compared to that of 'small' particles (i.e. smaller than 1 nm). We thus find that the hot denatured state is more compact and richer in secondary structure than the cold denatured state, since water at lower temperatures can form more hydrogen bonds than at high temperatures. Then, using Φ-value analysis we show that the structural differences between the hot and cold denatured states result in two alternative folding mechanisms. These findings thus illustrate how the analysis of water-protein hydrogen bonds can reveal the molecular origins of protein behaviours associated with the hydrophobic effect.

Molecular Recognition by Templated Folding of an Intrinsically Disordered Protein

Intrinsically disordered proteins often become structured upon interacting with their partners. The mechanism of this ‘folding upon binding’ process, however, has not been fully characterised yet. Here we present a study of the folding of the intrinsically disordered transactivation domain of c-Myb (c-Myb) upon binding its partner KIX. By determining the structure of the folding transition state for the binding of wild-type and three mutational variants of KIX, we found a remarkable plasticity of the folding pathway of c-Myb. To explain this phenomenon, we show that the folding of c-Myb is templated by the structure of KIX. This adaptive folding behaviour, which occurs by heterogeneous nucleation, differs from the robust homogeneous nucleation typically observed for globular proteins. We suggest that this templated folding mechanism may enable intrinsically disordered proteins to achieve specific and reliable binding with multiple partners while avoiding aberrant interactions.

Interplay between partner and ligand facilitates the folding and binding of an intrinsically disordered protein

Protein-protein interactions are at the heart of regulatory and signaling processes in the cell. In many interactions, one or both proteins are disordered before association. However, this disorder in the unbound state does not prevent many of these proteins folding to a well-defined, ordered structure in the bound state. Here we examine a typical system, where a small disordered protein (PUMA, p53 upregulated modulator of apoptosis) folds to an α-helix when bound to a groove on the surface of a folded protein (MCL-1, induced myeloid leukemia cell differentiation protein). We follow the association of these proteins using rapid-mixing stopped flow, and examine how the kinetic behavior is perturbed by denaturant and carefully chosen mutations. We demonstrate the utility of methods developed for the study of monomeric protein folding, including β-Tanford values, Leffler α, Φ-value analysis, and coarse-grained simulations, and propose a self-consistent mechanism for binding. Folding of the disordered protein before binding does not appear to be required and few, if any, specific interactions are required to commit to association. The majority of PUMA folding occurs after the transition state, in the presence of MCL-1. We also examine the role of the side chains of folded MCL-1 that make up the binding groove and find that many favor equilibrium binding but, surprisingly, inhibit the association process.

Hydrophobic interactions between the voltage sensor and pore mediate inactivation in Kv11.1 channels

Kv11.1 channels are critical for the maintenance of a normal heart rhythm. The flow of potassium ions through these channels is controlled by two voltage-regulated gates, termed "activation" and "inactivation," located at opposite ends of the pore. Crucially in Kv11.1 channels, inactivation gating occurs much more rapidly, and over a distinct range of voltages, compared with activation gating. Although it is clear that the fourth transmembrane segments (S4), within each subunit of the tetrameric channel, are important for controlling the opening and closing of the activation gate, their role during inactivation gating is much less clear. Here, we use rate equilibrium free energy relationship (REFER) analysis to probe the contribution of the S4 "voltage-sensor" helix during inactivation of Kv11.1 channels. Contrary to the important role that charged residues play during activation gating, it is the hydrophobic residues (Leu529, Leu530, Leu532, and Val535) that are the key molecular determinants of inactivation gating. Within the context of an interconnected multi-domain model of Kv11.1 inactivation gating, our REFER analysis indicates that the S4 helix and the S4-S5 linker undergo a conformational rearrangement shortly after that of the S5 helix and S5P linker, but before the S6 helix. Combining REFER analysis with double mutant cycle analysis, we provide evidence for a hydrophobic interaction between residues on the S4 and S5 helices. Based on a Kv11.1 channel homology model, we propose that this hydrophobic interaction forms the basis of an intersubunit coupling between the voltage sensor and pore domain that is an important mediator of inactivation gating.

Structure of the transition state for the binding of c-Myb and KIX highlights an unexpected order for a disordered system

A classical dogma of molecular biology dictates that the 3D structure of a protein is necessary for its function. However, a considerable fraction of the human proteome, although functional, does not adopt a defined folded state under physiological conditions. These intrinsically disordered proteins tend to fold upon binding to their partners with a molecular mechanism that is elusive to experimental characterization. Indeed, although many hypotheses have been put forward, the functional role (if any) of disorder in these intrinsically denatured systems is still shrouded in mystery. Here, we characterize the structure of the transition state of the binding-induced folding in the reaction between the KIX domain of the CREB-binding protein and the transactivation domain of c-Myb. The analysis, based on the characterization of a series of conservative site-directed mutants, reveals a very high content of native-like structure in the transition state and indicates that the recognition between KIX and c-Myb is geometrically precise. The implications of our results in the light of previous work on intrinsically unstructured systems are discussed.

Pore helices play a dynamic role as integrators of domain motion during Kv11.1 channel inactivation gating

Proteins that form ion-selective pores in the membrane of cells are integral to many rapid signaling processes, including regulating the rhythm of the heartbeat. In potassium channels, the selectivity filter is critical for both endowing an exquisite selectivity for potassium ions, as well as for controlling the flow of ions through the pore. Subtle rearrangements in the complex hydrogen-bond network that link the selectivity filter to the surrounding pore helices differentiate conducting (open) from nonconducting (inactivated) conformations of the channel. Recent studies suggest that beyond the selectivity filter, inactivation involves widespread rearrangements of the channel protein. Here, we use rate equilibrium free energy relationship analysis to probe the structural changes that occur during selectivity filter gating in Kv11.1 channels, at near atomic resolution. We show that the pore helix plays a crucial dynamic role as a bidirectional interface during selectivity filter gating. We also define the molecular bases of the energetic coupling between the pore helix and outer helix of the pore domain that occurs early in the transition from open to inactivated states, as well as the coupling between the pore helix and inner helix late in the transition. Our data demonstrate that the pore helices are more than just static structural elements supporting the integrity of the selectivity filter; instead they play a crucial dynamic role during selectivity filter gating.

Side-chain interactions form late and cooperatively in the binding reaction between disordered peptides and PDZ domains

Intrinsically disordered proteins are very common and mediate numerous protein-protein and protein-DNA interactions. While it is clear that these interactions are instrumental for the life of the mammalian cell, there is a paucity of data regarding their molecular binding mechanisms. Here we have used short peptides as a model system for intrinsically disordered proteins. Linear free energy relationships based on rate and equilibrium constants for the binding of these peptides to ordered target proteins, PDZ domains, demonstrate that native side-chain interactions form mainly after the rate-limiting barrier for binding and in a cooperative fashion. This finding suggests that these disordered peptides first form a weak encounter complex with non-native interactions. The data do not support the recent notion that the affinities of intrinsically disordered proteins toward their targets are generally governed by their association rate constants. Instead, we observed the opposite for peptide-PDZ interactions, namely, that changes in K(d) correlate with changes in k(off).

Linear rate-equilibrium relations arising from ion channel-bilayer energetic coupling

Linear rate-equilibrium (RE) relations, also known as linear free energy relations, are widely observed in chemical reactions, including protein folding, enzymatic catalysis, and channel gating. Despite the widespread occurrence of linear RE relations, the principles underlying the linear relation between changes in activation and equilibrium energy in macromolecular reactions remain enigmatic. When examining amphiphile regulation of gramicidin channel gating in lipid bilayers, we noted that the gating process could be described by a linear RE relation with a simple geometric interpretation. This description is possible because the gating process provides a well-understood reaction, in which structural changes in a bilayer-embedded model protein can be studied at the single-molecule level. It is thus possible to obtain quantitative information about the energetics of the reaction transition state and its position on a spatial coordinate. It turns out that the linear RE relation for the gramicidin monomer-dimer reaction can be understood, and the quantitative relation between changes in activation energy and equilibrium energy can be interpreted, by considering the effects of amphiphiles on the changes in bilayer elastic energy associated with channel gating. We are not aware that a similar simple mechanistic explanation of a linear RE relation has been provided for a chemical reaction in a macromolecule. RE relations generally should be useful for examining how amphiphile-induced changes in bilayer properties modulate membrane protein folding and function, and for distinguishing between direct (e.g., due to binding) and indirect (bilayer-mediated) effects.

Protein folding mechanisms studied by pulsed oxidative labeling and mass spectrometry

Deciphering the mechanisms of protein folding remains a considerable challenge. In this review we discuss the application of pulsed oxidative labeling for tracking protein structural changes in a time-resolved fashion. Exposure to a microsecond OH pulse at selected time points during folding induces the oxidation of solvent-accessible side chains, whereas buried residues are protected. Oxidative modifications can be detected by mass spectrometry. Folding is associated with dramatic accessibility changes, and therefore this method can provide detailed mechanistic insights. Solvent accessibility patterns are complementary to H/D exchange investigations, which report on the extent of hydrogen bonding. This review highlights the application of pulsed OH labeling to soluble proteins as well as membrane proteins.

Mapping the sequence of conformational changes underlying selectivity filter gating in the K(v)11.1 potassium channel

The potassium channel selectivity filter both discriminates between K(+) and sodium ions and contributes to gating of ion flow. Static structures of conducting (open) and nonconducting (inactivated) conformations of this filter are known; however, the sequence of protein rearrangements that connect these two states is not. We show that closure of the selectivity filter gate in the human K(v)11.1 K(+) channel (also known as hERG, for ether-a-go-go-related gene), a key regulator of the rhythm of the heartbeat, is initiated by K(+) exit, followed in sequence by conformational rearrangements of the pore domain outer helix, extracellular turret region, voltage sensor domain, intracellular domains and pore domain inner helix. In contrast to the simple wave-like sequence of events proposed for opening of ligand-gated ion channels, a complex spatial and temporal sequence of widespread domain motions connect the open and inactivated states of the K(v)11.1 K(+) channel.

Relationships between structural dynamics and functional kinetics in oligomeric membrane receptors

Recent efforts to broaden understanding of the molecular mechanisms of membrane receptors in signal transduction make use of rate-equilibrium free-energy relationships (REFERs), previously applied to chemical reactions, enzyme kinetics, and protein folding. For oligomeric membrane receptors, we distinguish between a), the Leffler parameter α_L, to characterize the global transition state for the interconversion between conformations; and b), the Fersht parameter, ϕ_F, to assign the degree of progression of individual residue positions at the transition state. For both α_L and ϕ_F, insights are achieved by using harmonic energy profiles to reflect the dynamic nature of proteins, as illustrated with single-channel results reported for normal and mutant nicotinic receptors. We also describe new applications of α_L based on published results. For large-conductance calcium-activated potassium channels, data are satisfactorily fit with an α_L value of 0.65, in accord with REFERs. In contrast, results reported for the flip conformational state of glycine and nicotinic receptors are in disaccord with REFERs, since they yield α_L values outside the usual limits of 0–1. Concerning published ϕ_F values underlying the conformational wave hypothesis for nicotinic receptors, we note that interpretations may be complicated by variations in the width of harmonic energy profiles.

Time-dependent changes in side-chain solvent accessibility during cytochrome c folding probed by pulsed oxidative labeling and mass spectrometry

The current work employs a novel approach for characterizing structural changes during the refolding of acid-denatured cytochrome c (cyt c). At various time points (ranging from 10 ms to 5 min) after a pH jump from 2 to 7, the protein is exposed to a microsecond hydroxyl radical (.OH) pulse that induces oxidative labeling of solvent-exposed side chains. Most of the covalent modifications appear as +16-Da adducts that are readily detectable by mass spectrometry. The overall extent of labeling decreases as folding proceeds, reflecting dramatic changes in the accessibility of numerous residues. Peptide mapping and tandem mass spectrometry reveal that the side chains of C14, C17, H33, F46, Y48, W59, M65, Y67, Y74, M80, I81, and Y97 are among the dominant sites of oxidation. Temporal changes in the accessibility of these residues are consistent with docking of the N- and C-terminal helices as early as 10 ms. However, structural reorganization at the helix interface takes place up to at least 1 s. Initial misligation of the heme iron by H33 leads to distal crowding, giving rise to low solvent accessibility of the displaced (native) M80 ligand and the adjacent I81. W59 retains a surprisingly high level of accessibility long into the folding process, indicating the presence of packing defects in the hydrophobically collapsed core. Overall, the results of this work are consistent with previous hydrogen/deuterium exchange studies that proposed a foldon-mediated mechanism. The structural data obtained by .OH labeling monitor the packing and burial of side chains, whereas hydrogen/deuterium exchange primarily monitors the formation of secondary structure elements. Hence, the two approaches yield complementary information. Considering the very short time scale of pulsed oxidative labeling, an extension of the approach used here to sub-millisecond folding studies should be feasible.

Relationship between nucleotide binding and ion channel gating in cystic fibrosis transmembrane conductance regulator

We have employed rate-equilibrium free energy relationship (REFER) analysis to characterize the dynamic events involved in the allosteric regulation of cystic fibrosis transmembrane conductance regulator (CFTR) function. A wide range of different hydrolysable and poorly hydrolysable nucleoside triphosphates were used to elucidate the role of ATP hydrolysis in CFTR function. The linearity of the REFER plots and Phi values near unity for all ligands tested implies that CFTR channel gating is a reversible thermally driven process with all structural reorganization in the binding site(s) completed prior to channel opening. This is consistent with the requirement for nucleotide binding for channel opening. However, the channel structural transition from the open to the closed state occurs independently of any events in the binding sites. Similar results were obtained on substitution of amino acids at coupling joints between both nucleotide binding domains (NBD) and cytoplasmic loops (CL) in opposite halves of the protein, indicating that any structural reorganization there also had occurred in the channel closed state. The fact that fractional Phi values were not observed in either of these distant sites suggests that there may not be a deterministic 'lever-arm' mechanism acting between nucleotide binding sites and the channel gate. These findings favour a stochastic coupling between binding and gating in which all structural transitions are thermally driven processes. We speculate that increase of channel open state probability is due to reduction of the number of the closed state configurations available after physical interaction between ligand bound NBDs and the channel.

Minimal folding pathways for coarse-grained biopolymer fragments

The minimal folding pathway or trajectory for a biopolymer can be defined as the transformation that minimizes the total distance traveled between a folded and an unfolded structure. This involves generalizing the usual Euclidean distance from points to one-dimensional objects such as a polymer. We apply this distance here to find minimal folding pathways for several candidate protein fragments, including the helix, the beta-hairpin, and a nonplanar structure where chain noncrossing is important. Comparing the distances traveled with root mean-squared distance and mean root-squared distance, we show that chain noncrossing can have large effects on the kinetic proximity of apparently similar conformations. Structures that are aligned to the beta-hairpin by minimizing mean root-squared distance, a quantity that closely approximates the true distance for long chains, show globally different orientation than structures aligned by minimizing root mean-squared distance.

Combined kinetic and thermodynamic analysis of alpha-helical membrane protein unfolding

The analytical toolkit developed for investigations into water-soluble protein folding has yet to be applied in earnest to membrane proteins. A major problem is the difficulty in collecting kinetic data, which are crucial to understanding any reaction. Here, we combine kinetic and thermodynamic studies of the reversible unfolding of an alpha-helical membrane protein to provide a definitive value for the reaction free energy and a means to probe the transition state. Our analyses show that the major unfolding step in the SDS-induced denaturation of bacteriorhodopsin involves a reduction in alpha-helical structure and proceeds with a large free-energy change; both our equilibrium and kinetic measurements predict that the free energy of unfolding in the absence of denaturant is +20 kcal.mol(-1), with an associated m-value of 25 kcal.mol(-1). The rate of unfolding in the absence of denaturant, k(u)(H(2)O), is surprisingly very slow ( approximately 10(-15) s(-1)). The kinetics also give information on the transition state for this major unfolding step, with a value for beta (m(f)/[m(f) + m(u)]) of approximately 0.1, indicating that the transition state is close to the unfolded state. We thus present a basis for mapping the structural and energetic properties of membrane protein folding by mutagenesis and classical kinetics.

Conformational dynamics of the alphaM3 transmembrane helix during acetylcholine receptor channel gating

Muscle acetylcholine receptors are synaptic ion channels that "gate" between closed- and open-channel conformations. We used Phi-value analysis to probe the transition state of the diliganded gating reaction with regard to residues in the M3, membrane-spanning helix of the muscle acetylcholine receptor alpha-subunit. Phi (a fraction between 1 and 0) parameterizes the extent to which a mutation changes the opening versus the closing rate constant and, for a linear reaction mechanism, the higher the Phi-value, the "earlier" the gating motion. In the upper half of alphaM3 the gating motions of all five tested residues were temporally correlated (Phi approximately 0.30) and serve to link structural changes occurring at the middle of the M2, pore-lining helix with those occurring at the interface of the extracellular and transmembrane domains. alphaM3 belongs to a complex and diverse set of synchronously moving parts that change structure relatively late in the channel-opening process. The propagation of the gating Brownian conformational cascade has a complex spatial distribution in the transmembrane domain.

Prediction of folding transition-state position (βT) of small, two-state proteins from local secondary structure content

Folding kinetics of proteins is governed by the free energy and position of transition states. But attempts to predict the position of folding transition state on reaction pathway from protein structure have been met with only limited success, unlike the folding-rate prediction. Here, we find that the folding transition-state position is related to the secondary structure content of native two-state proteins. We present a simple method for predicting the transition-state position from their alpha-helix, turn and polyproline secondary structures. The method achieves 81% correlation with experiment over 24 small, two-state proteins, suggesting that the local secondary structure content, especially for content of alpha-helix, is a determinant of the solvent accessibility of the transition state ensemble and size of folding nucleus.

Conformational entropy of alanine versus glycine in protein denatured states

The presence of a solvent-exposed alanine residue stabilizes a helix by 0.4-2 kcal.mol(-1) relative to glycine. Various factors have been suggested to account for the differences in helical propensity, from the higher conformational freedom of glycine sequences in the unfolded state to hydrophobic and van der Waals' stabilization of the alanine side chain in the helical state. We have performed all-atom molecular dynamics simulations with explicit solvent and exhaustive sampling of model peptides to address the backbone conformational entropy difference between Ala and Gly in the denatured state. The mutation of Ala to Gly leads to an increase in conformational entropy equivalent to approximately 0.4 kcal.mol(-1) in a fully flexible denatured, that is, unfolded, state. But, this energy is closely counterbalanced by the (measured) difference in free energy of transfer of the glycine and alanine side chains from the vapor phase to water so that the unfolded alanine- and glycine-containing peptides are approximately isoenergetic. The helix-stabilizing propensity of Ala relative to Gly thus mainly results from more favorable interactions of Ala in the folded helical structure. The small difference in energetics in the denatured states means that the Phi-values derived from Ala --> Gly scanning of helices are a very good measure of the extent of formation of structure in proteins with little residual structure in the denatured state.

Folding transition-state and denatured-state ensembles of FSD-1 from folding and unfolding simulations

Characterization of the folding transition-state ensemble and the denatured-state ensemble is an important step toward a full elucidation of protein folding mechanisms. We report herein an investigation of the free-energy landscape of FSD-1 protein by a total of four sets of folding and unfolding molecular dynamics simulations with explicit solvent. The transition-state ensemble was initially identified from unfolding simulations at 500 K and was verified by simulations at 300 K starting from the ensemble structures. The denatured-state ensemble and the early-stage folding were studied by a combination of unfolding simulations at 500 K and folding simulations at 300 K starting from the extended conformation. A common feature of the transition-state ensemble was the substantial formation of the native secondary structures, including both the alpha-helix and beta-sheet, with partial exposure of the hydrophobic core in the solvent. Both the native and non-native secondary structures were observed in the denatured-state ensemble and early-stage folding, consistent with the smooth experimental melting curve. Interestingly, the contact orders of the transition-state ensemble structures were similar to that of the native structure and were notably lower than those of the compact structures found in early-stage folding, implying that chain and topological entropy might play significant roles in protein folding. Implications for FSD-1 folding mechanisms and the rate-limiting step are discussed. Analyses further revealed interesting non-native interactions in the denatured-state ensemble and early-stage folding and the possibility that destabilization of these interactions could help to enhance the stability and folding rate of the protein.

Agonist-activated ion channels

This paper looks at ion channels as an example of the pharmacologist's stock in trade, the action of an agonist on a receptor to produce a response. Looked at in this way, ion channels have been helpful because they are still the only system which is simple enough for quantitative investigation of transduction mechanisms. A short history is given of attempts to elucidate what happens between the time when agonist first binds, and the time when the channel opens.

The contribution of the residues from the main hydrophobic core of ribonuclease A to its pressure-folding transition state

The role of hydrophobic interactions established by the residues that belong to the main hydrophobic core of ribonuclease A in its pressure-folding transition state was investigated using the Phi-value method. The folding kinetics was studied using pressure-jump techniques both in the pressurization and depressurization directions. The ratio between the folding activation volume and the reaction volume (beta p-value), which is an index of the compactness or degree of solvation of the transition state, was calculated. All the positions analyzed presented fractional Phi f-values, and the lowest were those corresponding to the most critical positions for the ribonuclease A stability. The structure of the transition state of the hydrophobic core of ribonuclease A, from the point of view of formed interactions, is a relatively, uniformly expanded form of the folded structure with a mean Phi f-value of 0.43. This places it halfway between the folded and unfolded states. On the other hand, for the variants, the average of beta p-values is 0.4, suggesting a transition state that is 40% native-like. Altogether the results suggest that the pressure-folding transition state of ribonuclease A looks like a collapsed globule with some secondary structure and a weakened hydrophobic core. A good correlation was found between the Phi f-values and the Deltabeta p-values. Although the nature of the transition state inferred from pressure-induced folding studies and the results of the protein engineering method have been reported to be consistent for other proteins, to the best of our knowledge this is the first direct comparison using a set of mutants.

Pressure as a tool to study protein-unfolding/refolding processes: The case of ribonuclease A

This paper gives an overview of the application of high-pressure to study the folding/unfolding processes of proteins using Ribonuclease A as a model protein. A particular focus is the study of pressure-equilibrium unfolding and folding kinetics using variants and the information obtained by comparing these with the wild-type enzyme.

Determination of an ensemble of structures representing the intermediate state of the bacterial immunity protein Im7

We present a detailed structural characterization of the intermediate state populated during the folding and unfolding of the bacterial immunity protein Im7. We achieve this result by incorporating a variety of experimental data available for this species in molecular dynamics simulations. First, we define the structure of the exchange-competent intermediate state of Im7 by using equilibrium hydrogen-exchange protection factors. Second, we use this ensemble to predict Phi-values and compare the results with the experimentally determined Phi-values of the kinetic refolding intermediate. Third, we predict chemical-shift measurements and compare them with the measured chemical shifts of a mutational variant of Im7 for which the kinetic folding intermediate is the most stable state populated at equilibrium. Remarkably, we found that the properties of the latter two species are predicted with high accuracy from the exchange-competent intermediate that we determined, suggesting that these three states are characterized by a similar architecture in which helices I, II, and IV are aligned in a native-like, but reorganized, manner. Furthermore, the structural ensemble that we obtained enabled us to rationalize the results of tryptophan fluorescence experiments in the WT protein and a series of mutational variants. The results show that the integration of diverse sets of experimental data at relatively low structural resolution is a powerful approach that can provide insights into the structural organization of this conformationally heterogeneous three-helix intermediate with unprecedented detail and highlight the importance of both native and non-native interactions in stabilizing its structure.

A critical assessment of the topomer search model of protein folding using a continuum explicit-chain model with extensive conformational sampling

Recently, a series of closely related theoretical constructs termed the "topomer search model" (TSM) has been proposed for the folding mechanism of small, single-domain proteins. A basic assumption of the proposed scenarios is that the rate-limiting step in folding is an essentially unbiased, diffusive search for a conformational state called the native topomer defined by an overall native-like topological pattern. Successes in correlating TSM-predicted folding rates with that of real proteins have been interpreted as experimental support for the model. To better delineate the physics entailed, key TSM concepts are examined here using extensive Langevin dynamics simulations of continuum C(alpha) chain models. The theoretical native topomers of four experimentally well-studied two-state proteins are characterized. Consistent with the TSM perspective, we found that the sizes of the native topomers increase with experimental folding rate. However, a careful determination of the corresponding probabilities that the native topomers are populated during a random search fails to reproduce the previously predicted folding rates. Instead, our results indicate that an unbiased TSM search for the native topomer amounts to a Levinthal-like process that would take an impossibly long average time to complete. Furthermore, intraprotein contacts in all four native topomers considered exhibit no apparent correlation with the experimental phi-values determined from the folding kinetics of these proteins. Thus, the present findings suggest that certain basic, generic yet essential energetic features in protein folding are not accounted for by TSM scenarios to date.

Determination of the folding transition states of barnase by using PhiI-value-restrained simulations validated by double mutant PhiIJ-values

The protein barnase folds from the denatured state into its native conformation via a high-energy intermediate. Using PhiI-values determined experimentally from single-point mutations as restraints in all-atom molecular dynamics simulations, we have determined ensembles of structures corresponding to the transition states for the formation of the folding intermediate and its conversion into the native state. We have also introduced a stringent validation of the approach used to calculate such structures by predicting interaction PhiIJ-values determined experimentally from a series of double-mutant cycles. The ensembles that we have obtained illustrate the heterogeneity in the nucleation-condensation process by which barnase folds. Obligatory steps of this process include the sequential formation of two folding nuclei, which correspond to the two main hydrophobic cores of the protein. Nonobligatory steps include the elongation of the strand beta1 and the formation of the helix alpha2. The results confirm that the use of experimental observables such as PhiI-values as restraints in molecular dynamics simulations is a powerful general strategy to characterize the relatively heterogeneous structural ensembles that populate nonnative regions of the energy landscapes of proteins.

Kinetic folding mechanism of PDZ2 from PTP-BL

PDZ domains represent a large family of protein-interaction modules associated with a variety of unrelated proteins with different functions. We report a complete characterization of the kinetic folding mechanism of a fluorescent variant of PDZ2 from PTP-BL, investigated under a variety of different experimental conditions. For this purpose, we engineered a fluorescent variant of this protein Y43W (called pseudo-wild-type, pWT43). The results suggest the presence of a high-energy intermediate in the folding of PDZ2, as revealed by a pronounced non-linear dependence of the unfolding rate constant on denaturant concentration. Such an intermediate may or may not be detectable depending on the experimental conditions, giving rise to apparent two-state folding under stabilizing conditions (e.g. in the presence of sodium sulfate). Interestingly, even under these conditions, three-state folding can be restored by selectively destabilizing the native-like rate-limiting barrier by one specific mutation (V44A). Finally, we show that data taken on pWT43 under different experimental conditions (e.g. different pH values from 2.1 to 8.0 or in the presence of a stabilizing salt) and also data on a site-directed conservative mutant can be rationalized in terms of a simple reaction scheme involving a single set of intermediates and transition states.