Do proteins with similar folds have similar transition state structures? A diffuse transition state of the 169 residue apoflavodoxin

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.

Protposer: the web server that readily proposes protein stabilizing mutations with high PPV

Protein stability is a requisite for most biotechnological and medical applications of proteins. As natural proteins tend to suffer from a low conformational stability ex vivo, great efforts have been devoted toward increasing their stability through rational design and engineering of appropriate mutations. Unfortunately, even the best currently used predictors fail to compute the stability of protein variants with sufficient accuracy and their usefulness as tools to guide the rational stabilisation of proteins is limited. We present here Protposer, a protein stabilising tool based on a different approach. Instead of quantifying changes in stability, Protposer uses structure- and sequence-based screening modules to nominate candidate mutations for subsequent evaluation by a logistic regression model, carefully trained to avoid overfitting. Thus, Protposer analyses PDB files in search for stabilization opportunities and provides a ranked list of promising mutations with their estimated success rates (eSR), their probabilities of being stabilising by at least 0.5 kcal/mol. The agreement between eSRs and actual positive predictive values (PPV) on external datasets of mutations is excellent. When Protposer is used with its Optimal kappa selection threshold, its PPV is above 0.7. Even with less stringent thresholds, Protposer largely outperforms FoldX, Rosetta and PoPMusiC. Indicating the PDB file of the protein suffices to obtain a ranked list of mutations, their eSRs and hints on the likely source of the stabilization expected. Protposer is a distinct, straightforward and highly successful tool to design protein stabilising mutations, and it is freely available for academic use at http://webapps.bifi.es/the-protposer.

Mutational Analysis of Protein Folding Transition States: Phi Values

The analysis of protein folding reactions by monitoring the kinetic effects of specifically designed single-point mutations, the so-termed phi-value analysis, has been a favorite technique to experimentally probe the mechanisms of protein folding. The idea behind phi-value analysis is that the effects that mutations have on the folding and unfolding rate constants report on the energetic/structural features of the folding transition state ensemble (TSE), which is the highest point in the free energy surface connecting the native and unfolded states, and thus the rate limiting step that ultimately defines the folding mechanism. For single-domain, two-state folding proteins, the general procedure to perform the phi-value analysis of protein folding is relatively simple to implement in the lab. Once the mutations have been produced and purified, the researcher needs to follow a few specific guidelines to perform the experiments and to analyze the data so produced. In this chapter, a step-by-step description of how to measure and interpret the effects induced by site-directed mutations on the folding and unfolding rate constants of a protein of interest is provided. Some possible solutions to the most typical problems that arise when performing phi-value analysis in the lab are also provided.

Mechanism of the small ATP-independent chaperone Spy is substrate specific

ATP-independent chaperones are usually considered to be holdases that rapidly bind to non-native states of substrate proteins and prevent their aggregation. These chaperones are thought to release their substrate proteins prior to their folding. Spy is an ATP-independent chaperone that acts as an aggregation inhibiting holdase but does so by allowing its substrate proteins to fold while they remain continuously chaperone bound, thus acting as a foldase as well. The attributes that allow such dual chaperoning behavior are unclear. Here, we used the topologically complex protein apoflavodoxin to show that the outcome of Spy's action is substrate specific and depends on its relative affinity for different folding states. Tighter binding of Spy to partially unfolded states of apoflavodoxin limits the possibility of folding while bound, converting Spy to a holdase chaperone. Our results highlight the central role of the substrate in determining the mechanism of chaperone action.

Direct examination of the relevance for folding, binding and electron transfer of a conserved protein folding intermediate

Near the minimum free energy basin of proteins where the native ensemble resides, partly unfolded conformations of slightly higher energy can be significantly populated under native conditions. It has been speculated that they play roles in molecular recognition and catalysis, but they might represent contemporary features of the evolutionary process without functional relevance. Obtaining conclusive evidence on these alternatives is difficult because it requires comparing the performance of a given protein when populating and when not populating one such intermediate, in otherwise identical conditions. Wild type apoflavodoxin populates under native conditions a partly unfolded conformation (10% of molecules) whose unstructured region includes the binding sites for the FMN cofactor and for redox partner proteins. We recently engineered a thermostable variant where the intermediate is no longer detectable. Using the wild type and variant, we assess the relevance of the intermediate comparing folding kinetics, cofactor binding kinetics, cofactor affinity, X-ray structure, intrinsic dynamics, redox potential of the apoflavodoxin-cofactor complex (Fld), its affinity for partner protein FNR, and electron transfer rate within the Fld/FNR physiological complex. Our data strongly suggest the intermediate state, conserved in long-chain apoflavodoxins, is not required for the correct assembly of flavodoxin nor does it contribute to shape its electron transfer properties. This analysis can be applied to evaluate other native basin intermediates.

Differences in the mechanical unfolding pathways of apo- and copper-bound azurins

Metalloproteins carry out diverse biological functions including metal transport, electron transfer, and catalysis. At present, the influence of metal cofactors on metalloprotein stability is not well understood. Here, we report the mechanical stability and unfolding pathway of azurin, a cupredoxin family protein with β-barrel topology and type I copper-binding centre. Single-molecule force spectroscopy (SMFS) experiments reveal 2-state and 3-state unfolding pathways for apo-azurin. The intermediate in the 3-state pathway occurs at an unfolding contour length of 7.5 nm from the native state. Steered molecular dynamics (SMD) simulations show that apo-azurin unfolds via a first transition state (TS) where β2Β–β8 and β7–β8 strand pairs rupture to form the intermediate, which subsequently unfolds by the collective rupture of remaining strands. SMFS experiments on holo-azurin exhibit an additional 4-state pathway besides the 2-state and 3-state pathways. The unfolding contour length leading to the first intermediate is 6.7 nm suggesting a sequestration of ~1 nm polypeptide chain length by the copper. SMD simulations reveal atomistic details of the copper sequestration and predict a combined β4–β7 pair and copper coordination sphere rupture to create the third TS in the 4-state pathway. Our systematic studies provide detailed mechanistic insights on modulation of protein mechanical properties by metal-cofactors.

Folding of proteins with a flavodoxin-like architecture

The flavodoxin-like fold is a protein architecture that can be traced back to the universal ancestor of the three kingdoms of life. Many proteins share this α-β parallel topology and hence it is highly relevant to illuminate how they fold. Here, we review experiments and simulations concerning the folding of flavodoxins and CheY-like proteins, which share the flavodoxin-like fold. These polypeptides tend to temporarily misfold during unassisted folding to their functionally active forms. This susceptibility to frustration is caused by the more rapid formation of an α-helix compared to a β-sheet, particularly when a parallel β-sheet is involved. As a result, flavodoxin-like proteins form intermediates that are off-pathway to native protein and several of these species are molten globules (MGs). Experiments suggest that the off-pathway species are of helical nature and that flavodoxin-like proteins have a nonconserved transition state that determines the rate of productive folding. Folding of flavodoxin from Azotobacter vinelandii has been investigated extensively, enabling a schematic construction of its folding energy landscape. It is the only flavodoxin-like protein of which cotranslational folding has been probed. New insights that emphasize differences between in vivo and in vitro folding energy landscapes are emerging: the ribosome modulates MG formation in nascent apoflavodoxin and forces this polypeptide toward the native state.

Contribution to the prediction of the fold code: application to immunoglobulin and flavodoxin cases

Background

Folding nucleus of globular proteins formation starts by the mutual interaction of a group of hydrophobic amino acids whose close contacts allow subsequent formation and stability of the 3D structure. These early steps can be predicted by simulation of the folding process through a Monte Carlo (MC) coarse grain model in a discrete space. We previously defined MIRs (Most Interacting Residues), as the set of residues presenting a large number of non-covalent neighbour interactions during such simulation. MIRs are good candidates to define the minimal number of residues giving rise to a given fold instead of another one, although their proportion is rather high, typically [15-20]% of the sequences. Having in mind experiments with two sequences of very high levels of sequence identity (up to 90%) but different folds, we combined the MIR method, which takes sequence as single input, with the “fuzzy oil drop” (FOD) model that requires a 3D structure, in order to estimate the residues coding for the fold. FOD assumes that a globular protein follows an idealised 3D Gaussian distribution of hydrophobicity density, with the maximum in the centre and minima at the surface of the “drop”. If the actual local density of hydrophobicity around a given amino acid is as high as the ideal one, then this amino acid is assigned to the core of the globular protein, and it is assumed to follow the FOD model. Therefore one obtains a distribution of the amino acids of a protein according to their agreement or rejection with the FOD model.

Results

We compared and combined MIR and FOD methods to define the minimal nucleus, or keystone, of two populated folds: immunoglobulin-like (Ig) and flavodoxins (Flav). The combination of these two approaches defines some positions both predicted as a MIR and assigned as accordant with the FOD model. It is shown here that for these two folds, the intersection of the predicted sets of residues significantly differs from random selection. It reduces the number of selected residues by each individual method and allows a reasonable agreement with experimentally determined key residues coding for the particular fold. In addition, the intersection of the two methods significantly increases the specificity of the prediction, providing a robust set of residues that constitute the folding nucleus.

Diffuse transition state structure for the unfolding of a leucine-rich repeat protein

Tandem-repeat proteins, such as leucine-rich repeats, comprise arrays of small structural motifs that pack in a linear fashion to produce elongated architectures. They lack contacts between residues that are distant in primary sequence, a feature that distinguishes them from the complex topologies of globular proteins. Here we have investigated the unfolding pathway of the leucine-rich repeat domain of the mRNA export protein TAP (TAPLRR) using Φ-value analysis. Whereas most of the tandem-repeat proteins studied to date have been found to unfold via a polarised mechanism in which only a small, localised number of repeats are structured in the transition state, the unfolding mechanism of TAPLRR is more diffuse in nature. In the transition state for unfolding of TAPLRR, three of the four LRRs are highly structured and non-native interactions are formed within the N-terminal α-helical cap and the first LRR. Thus, the α-helical cap plays an important role in which non-native interactions are required to provide a scaffold for the LRRs to pack against in the folding reaction.

Early folding events protect aggregation-prone regions of a β-rich protein

Protein folding and aggregation inevitably compete with one another. This competition is even keener for proteins with frustrated landscapes, such as those rich in β structure. It is interesting that, despite their rugged energy landscapes and high β sheet content, intracellular lipid-binding proteins (iLBPs) appear to successfully avoid aggregation, as they are not implicated in aggregation diseases. In this study, we used a canonical iLBP, cellular retinoic acid-binding protein 1 (CRABP1), to understand better how folding is favored over aggregation. Analysis of folding kinetics of point mutants reveals that the folding pathway of CRABP1 involves early barrel closure. This folding mechanism protects sequences in CRABP1 that comprise cores of aggregates as identified by nuclear magnetic resonance. The amino acid conservation pattern in other iLBPs suggests that early barrel closure may be a general strategy for successful folding and minimization of aggregation. We suggest that folding mechanisms in general may incorporate steps that disfavor aggregation.

Defining the nature of thermal intermediate in 3 state folding proteins: apoflavodoxin, a study case

The early stages of the thermal unfolding of apoflavodoxin have been determined by using atomistic multi microsecond-scale molecular dynamics (MD) simulations complemented with a variety of experimental techniques. Results strongly suggest that the intermediate is reached very early in the thermal unfolding process and that it has the properties of an “activated” form of the native state, where thermal fluctuations in the loops break loop-loop contacts. The unrestrained loops gain then kinetic energy corrupting short secondary structure elements without corrupting the core of the protein. The MD-derived ensembles agree with experimental observables and draw a picture of the intermediate state inconsistent with a well-defined structure and characteristic of a typical partially disordered protein. Our results allow us to speculate that proteins with a well packed core connected by long loops might behave as partially disordered proteins under native conditions, or alternatively behave as three state folders. Small details in the sequence, easily tunable by evolution, can yield to one or the other type of proteins.

A simplistic view of protein structure tends to emphasize the opposition between the native state and the denatured ensemble of unfolded conformations. In addition to these extreme conformations, proteins subjected to a variety of perturbations often populate alternative partly unfolded conformations, some of which are close in energy to the native state and, accordingly, can be populated under native or quasi-native conditions. There is increasing evidence that these “perturbed” conformations participate in protein function or, in some cases, are related to the outcome of folding diseases. We have used the “state of the art” molecular dynamics combined with a variety of experimental techniques to characterize for the first time, to our knowledge, the thermal intermediate of a three-state folding protein (apoflavodoxin). Based on our results we have been able to suggest a general mechanism of thermal unfolding in complex proteins and to determine interesting links between thermal intermediates and partially unfolded proteins.

Why do protein folding rates correlate with metrics of native topology?

For almost 15 years, the experimental correlation between protein folding rates and the contact order parameter has been under scrutiny. Here, we use a simple simulation model combined with a native-centric interaction potential to investigate the physical roots of this empirical observation. We simulate a large set of circular permutants, thus eliminating dependencies of the folding rate on other protein properties (e.g. stability). We show that the rate-contact order correlation is a consequence of the fact that, in high contact order structures, the contact order of the transition state ensemble closely mirrors the contact order of the native state. This happens because, in these structures, the native topology is represented in the transition state through the formation of a network of tertiary interactions that are distinctively long-ranged.

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

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

Design and structure of an equilibrium protein folding intermediate: a hint into dynamical regions of proteins

Partly unfolded protein conformations close to the native state may play important roles in protein function and in protein misfolding. Structural analyses of such conformations which are essential for their fully physicochemical understanding are complicated by their characteristic low populations at equilibrium. We stabilize here with a single mutation the equilibrium intermediate of apoflavodoxin thermal unfolding and determine its solution structure by NMR. It consists of a large native region identical with that observed in the X-ray structure of the wild-type protein plus an unfolded region. Small-angle X-ray scattering analysis indicates that the calculated ensemble of structures is consistent with the actual degree of expansion of the intermediate. The unfolded region encompasses discontinuous sequence segments that cluster in the 3D structure of the native protein forming the FMN cofactor binding loops and the binding site of a variety of partner proteins. Analysis of the apoflavodoxin inner interfaces reveals that those becoming destabilized in the intermediate are more polar than other inner interfaces of the protein. Natively folded proteins contain hydrophobic cores formed by the packing of hydrophobic surfaces, while natively unfolded proteins are rich in polar residues. The structure of the apoflavodoxin thermal intermediate suggests that the regions of natively folded proteins that are easily responsive to thermal activation may contain cores of intermediate hydrophobicity.

Residue-specific analysis of frustration in the folding landscape of repeat beta/alpha protein apoflavodoxin

Flavodoxin adopts the common repeat beta/alpha topology and folds in a complex kinetic reaction with intermediates. To better understand this reaction, we analyzed a set of Desulfovibrio desulfuricans apoflavodoxin variants with point mutations in most secondary structure elements by in vitro and in silico methods. By equilibrium unfolding experiments, we first revealed how different secondary structure elements contribute to overall protein resistance to heat and urea. Next, using stopped-flow mixing coupled with far-UV circular dichroism, we probed how individual residues affect the amount of structure formed in the experimentally detected burst-phase intermediate. Together with in silico folding route analysis of the same point-mutated variants and computation of growth in nucleation size during early folding, computer simulations suggested the presence of two competing folding nuclei at opposite sides of the central beta-strand 3 (i.e., at beta-strands 1 and 4), which cause early topological frustration (i.e., misfolding) in the folding landscape. Particularly, the extent of heterogeneity in folding nuclei growth correlates with the in vitro burst-phase circular dichroism amplitude. In addition, phi-value analysis (in vitro and in silico) of the overall folding barrier to apoflavodoxin's native state revealed that native-like interactions in most of the beta-strands must form in transition state. Our study reveals that an imbalanced competition between the two sides of apoflavodoxin's central beta-sheet directs initial misfolding, while proper alignment on both sides of beta-strand 3 is necessary for productive folding.

Macromolecular crowding modulates folding mechanism of alpha/beta protein apoflavodoxin

Protein dynamics in cells may be different from those in dilute solutions in vitro, because the environment in cells is highly concentrated with other macromolecules. This volume exclusion because of macromolecular crowding is predicted to affect both equilibrium and kinetic processes involving protein conformational changes. To quantify macromolecular crowding effects on protein folding mechanisms, we investigated the folding energy landscape of an alpha/beta protein, apoflavodoxin, in the presence of inert macromolecular crowding agents, using in silico and in vitro approaches. By means of coarse-grained molecular simulations and topology-based potential interactions, we probed the effects of increased volume fractions of crowding agents (phi(c)) as well as of crowding agent geometry (sphere or spherocylinder) at high phi(c). Parallel kinetic folding experiments with purified Desulfovibro desulfuricans apoflavodoxin in vitro were performed in the presence of Ficoll (sphere) and Dextran (spherocylinder) synthetic crowding agents. In conclusion, we identified the in silico crowding conditions that best enhance protein stability, and discovered that upon manipulation of the crowding conditions, folding routes experiencing topological frustrations can be either enhanced or relieved. Our test-tube experiments confirmed that apoflavodoxin's time-resolved folding path is modulated by crowding agent geometry. Macromolecular crowding effects may be a tool for the manipulation of protein-folding and function in living cells.

Folding, stability and shape of proteins in crowded environments: experimental and computational approaches

How the crowded environment inside cells affects folding, stability and structures of proteins is a vital question, since most proteins are made and function inside cells. Here we describe how crowded conditions can be created in vitro and in silico and how we have used this to probe effects on protein properties. We have found that folded forms of proteins become more compact in the presence of macromolecular crowding agents; if the protein is aspherical, the shape also changes (extent dictated by native-state stability and chemical conditions). It was also discovered that the shape of the macromolecular crowding agent modulates the folding mechanism of a protein; in addition, the extent of asphericity of the protein itself is an important factor in defining its folding speed.

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.

Folding domain B of protein A on a dynamically partitioned free energy landscape

The B domain of staphylococcal protein A (BdpA) is a small helical protein that has been studied intensively in kinetics experiments and detailed computer simulations that include explicit water. The simulations indicate that BdpA needs to reorganize in crossing the transition barrier to facilitate folding its C-terminal helix (H3) onto the nucleus formed from helices H1 and H2. This process suggests frustration between two partially ordered forms of the protein, but recent varphi value measurements indicate that the transition structure is relatively constant over a broad range of temperatures. Here we develop a simplistic model to investigate the folding transition in which properties of the free energy landscape can be quantitatively compared with experimental data. The model is a continuation of the Muñoz-Eaton model to include the intermittency of contacts between structured parts of the protein, and the results compare variations in the landscape with denaturant and temperature to varphi value measurements and chevron plots of the kinetic rates. The topography of the model landscape (in particular, the feature of frustration) is consistent with detailed simulations even though variations in the varphi values are close to measured values. The transition barrier is smaller than indicated by the chevron data, but it agrees in order of magnitude with a similar alpha-carbon type of model. Discrepancies with the chevron plots are investigated from the point of view of solvent effects, and an approach is suggested to account for solvent participation in the model.

SIMPLE estimate of the free energy change due to aliphatic mutations: Superior predictions based on first principles

The bioinformatics revolution of the last decade has been instrumental in the development of empirical potentials to quantitatively estimate protein interactions for modeling and design. Although computationally efficient, these potentials hide most of the relevant thermodynamics in 5-to-40 parameters that are fitted against a large experimental database. Here, we revisit this longstanding problem and show that a careful consideration of the change in hydrophobicity, electrostatics, and configurational entropy between the folded and unfolded state of aliphatic point mutations predicts 20-30% less false positives and yields more accurate predictions than any published empirical energy function. This significant improvement is achieved with essentially no free parameters, validating past theoretical and experimental efforts to understand the thermodynamics of protein folding. Our first principle analysis strongly suggests that both the solute-solute van der Waals interactions in the folded state and the electrostatics free energy change of exposed aliphatic mutations are almost completely compensated by similar interactions operating in the unfolded ensemble. Not surprisingly, the problem of properly accounting for the solvent contribution to the free energy of polar and charged group mutations, as well as of mutations that disrupt the protein backbone remains open.