The largest protein observed to fold by two-state kinetic mechanism does not obey contact-order correlation

To probe the folding-energy landscape for a very large protein, we used Borrelia burgdorferi VlsE as a model. VlsE is a single-domain, predominantly alpha-helical protein with 341 residues. Remarkably, time-resolved folding and unfolding processes for VlsE follow two-state behavior. VlsE is by far the largest protein characterized that folds by a two-state kinetic mechanism. Thus, the common rule of thumb, that proteins larger than 110 residues fold by complex, multistate kinetic mechanisms, must be used with caution. In contrast with smaller helical proteins, the folding speed in water for VlsE is slow (5 +/- 2 s-1, pH 7, 20 degrees C) and does not agree (by 4 orders of magnitude in different directions) with the speeds predicted on the basis of native-state contact order and the topomer-search model. It is therefore questionable if the barrier height for folding is defined by gross topology for large two-state folders.     

Prediction of protein folding rates from primary sequences using hybrid sequence representation

The ability to predict protein folding rates constitutes an important step in understanding the overall folding mechanisms. Although many of the prediction methods are structure based, successful predictions can also be obtained from the sequence. We developed a novel method called prediction of protein folding rates (PPFR), for the prediction of protein folding rates from protein sequences. PPFR implements a linear regression model for each of the mainstream folding dynamics including two-, multi-, and mixed-state proteins. The proposed method provides predictions characterized by strong correlations with the experimental folding rates, which equal 0.87 for the two- and multistate proteins and 0.82 for the mixed-state proteins, when evaluated with out-of-sample jackknife test. Based on in-sample and out-of-sample tests, the PPFR's predictions are shown to be better than most of other sequence only and structure-based predictors and complementary to the predictions of the most recent sequence-based QRSM method. We show that simultaneous incorporation of several characteristics, including the sequence, physiochemical properties of residues, and predicted secondary structure provides improved quality. This hybridized prediction model was analyzed to reveal the complementary factors that can be used in tandem to predict folding rates. We show that bigger proteins require more time for folding, higher helical and coil content and the presence of Phe, Asn, and Gln may accelerate the folding process, the inclusion of Ile, Val, Thr, and Ser may slow down the folding process, and for the two-state proteins increased beta-strand content may decelerate the folding process. Finally, PPFR provides strong correlation when predicting sequences with low similarity.     

Importance of native-state topology for determining the folding rate of two-state proteins

Understanding the relationship between amino acid sequences and folding rate of proteins is a challenging task similar to protein folding problem. In this work, we have analyzed the relative importance of protein sequence and structure for predicting the protein folding rates in terms of amino acid properties and contact distances, respectively. We found that the parameters derived with protein sequence (physical-chemical, energetic, and conformational properties of amino acid residues) show very weak correlation (|r| < 0.39) with folding rates of 28 two-state proteins, indicating that the sequence information alone is not sufficient to understand the folding rates of two-state proteins. However, the maximum positive correlation obtained for the properties, number of medium-range contacts, and alpha-helical tendency reveals the importance of local interactions to initiate protein folding. On the other hand, a remarkable correlation (r varies from -0.74 to -0.88) has been obtained between structural parameters (contact order, long-range order, and total contact distance) and protein folding rates. Further, we found that the secondary structure content and solvent accessibility play a marginal role in determining the folding rates of two-state proteins. Multiple regression analysis carried out with the combination of three properties, beta-strand tendency, enthalpy change, and total contact distance improved the correlation to 0.92 with protein folding rates. The relative importance of existing methods along with multiple-regression model proposed in this work will be discussed. Our results demonstrate that the native-state topology is the major determinant for the folding rates of two-state proteins.     

Statistical mechanical treatment of protein conformation. I. Conformational properties of amino acids in proteins.

A statistical mechanical (one-dimensional Ising model) treatment, based on the dominance of short-range interactions, is developed in this series of papers; it is intended as an improvement over empirical prediction schemes for obtaining approximate initial conformations of proteins (to be used to try to deduce the native conformation by subsequent energy minimization). In the present paper, the statistical weights for a two-state model (alpha-helical and other conformations) and for a three-state model (alpha-helical, extended, and other conformations) are evaluated from x-ray data on 16 native proteins. The method for evaluating the statistical weights is presented. Asymmetric alpha-helical nucleation parameters are also evaluated for the 20 naturally occurring amino acids. On the basis of these statistical weights, the conformational properties of the twenty naturally occurring amino acids are discussed. The statistical weights evaluated from x-ray data are also discussed in comparison with experimental results on the helix--coil transition in polyamino acids in solution. The predominant role of short-range interactions, and some possible long-range effects in determining the statistical weights, are discussed in conjunction with the mechanism of protein folding.     

Improved theoretical description of protein folding kinetics from rotations in the phase space of relevant order parameters

A method is introduced to construct a better approximation for the reaction coordinate for protein folding from known order parameters. The folding of a two-state off-lattice alpha helical Go-type protein is studied using molecular dynamics simulations. Folding times are computed directly from simulation, as well as theoretically using an equation derived by considering Brownian-type dynamics for the putative reaction coordinate. Theoretical estimates of the folding time using the number of native contacts (Qn) as a reaction coordinate were seen to differ quite significantly from the true folding time of the protein. By considering the properties of the bimodal free energy surface of this protein as a function of Qn and another relevant coordinate for folding Q (the total number of contacts), we show that by introducing a rotation in the phase space of the order parameters Q and Qn, we can construct a new reaction coordinate q that leads to a fivefold improvement in the estimate of the folding rate. This new coordinate q, resulting from the rotation, lies along the line connecting the unfolded and folded ensemble minima of the free energy map plotted as a function of the original order parameters Q and Qn. Possible reasons for the remaining discrepancy between the folding time computed theoretically and from folding simulations are discussed.     

Mechanically untying a protein slipknot: multiple pathways revealed by force spectroscopy and steered molecular dynamics simulations

Protein structure is highly diverse when considering a wide range of protein types, helping to give rise to the multitude of functions that proteins perform. In particular, certain proteins are known to adopt a knotted or slipknotted fold. How such proteins undergo mechanical unfolding was investigated utilizing a combination of single molecule atomic force microscopy (AFM), protein engineering, and steered molecular dynamics (SMD) simulations to show the mechanical unfolding mechanism of the slipknotted protein AFV3-109. Our results reveal that the mechanical unfolding of AFV3-109 can proceed via multiple parallel unfolding pathways that all cause the protein slipknot to untie and the polypeptide chain to completely extend. These distinct unfolding pathways proceed via either a two- or three-state unfolding process involving the formation of a well-defined, stable intermediate state. SMD simulations predict the same contour length increments for different unfolding pathways as single molecule AFM results, thus providing a plausible molecular mechanism for the mechanical unfolding of AFV3-109. These SMD simulations also reveal that two-state unfolding is initiated from both the N- and C-termini, while three-state unfolding is initiated only from the C-terminus. In both pathways, the protein slipknot was untied during unfolding, and no tightened slipknot conformation was observed. Detailed analysis revealed that interactions between key structural elements lock the knotting loop in place, preventing it from shrinking and the formation of a tightened slipknot conformation. Our results demonstrate the bifurcation of the mechanical unfolding pathway of AFV3-109 and point to the generality of a kinetic partitioning mechanism for protein folding/unfolding.     

Protein folding kinetics: barrier effects in chemical and thermal denaturation experiments

Recent experimental work on fast protein folding brings about an intriguing paradox. Microsecond-folding proteins are supposed to fold near or at the folding speed limit (downhill folding), but yet their folding behavior seems to comply with classical two-state analyses, which imply the crossing of high free energy barriers. However, close inspection of chemical and thermal denaturation kinetic experiments in fast-folding proteins reveals systematic deviations from two-state behavior. Using a simple one-dimensional free energy surface approach we find that such deviations are indeed diagnostic of marginal folding barriers. Furthermore, the quantitative analysis of available fast-kinetic data indicates that many microsecond-folding proteins fold downhill in native conditions. All of these proteins are then promising candidates for an atom-by-atom analysis of protein folding using nuclear magnetic resonance.1 We also find that the diffusion coefficient for protein folding is strongly temperature dependent, corresponding to an activation energy of approximately 1 kJ.mol-1 per protein residue. As a consequence, the folding speed limit at room temperature is about an order of magnitude slower than the approximately 1 micros estimates from high-temperature T-jump experiments. Our analysis is quantitatively consistent with the available thermodynamic and kinetic data on slow two-state folding proteins and provides a straightforward explanation for the apparent fast-folding paradox.     

Two-state folding observed in individual protein molecules

The folding dynamics of small proteins are often described in terms of a simple two-state kinetic model. Within this notion, the behavior of individual molecules is expected to be stochastic, with a protein molecule residing in either the unfolded or the folded state for extended periods of time, with intermittent rapid jumps across the free energy barrier. However, a direct observation of this bistable behavior has not been made to date. Rather, previous reports of folding trajectories of individual proteins have shown an unexpected degree of complexity. This raises the question whether the simple kinetic properties derived from classical experiments on large ensembles of molecules are reflected in the folding paths taken by individual proteins. Here we report single-molecule folding/unfolding trajectories observed by fluorescence resonance energy transfer for a protein that meets all criteria of a two state-system. The trajectories, measured on molecules immobilized in lipid vesicles, demonstrate the anticipated bistable behavior, with steplike transitions between folded and unfolded conformations. They further allow us to put an upper bound on the barrier crossing time.     

Analysis and prediction of protein folding rates using quadratic response surface models

Understanding the relationship between amino acid sequences and folding rates of proteins is an important task in computational and molecular biology. In this work, we have systematically analyzed the composition of amino acid residues for proteins with different ranges of folding rates. We observed that the polar residues, Asn, Gln, Ser, and Lys, are dominant in fast folding proteins whereas the hydrophobic residues, Ala, Cys, Gly, and Leu, prefer to be in slow folding proteins. Further, we have developed a method based on quadratic response surface models for predicting the folding rates of 77 two- and three-state proteins. Our method showed a correlation of 0.90 between experimental and predicted protein folding rates using leave-one-out cross-validation method. The classification of proteins based on structural class improved the correlation to 0.98 and it is 0.99, 0.98, and 0.96, respectively, for all-alpha, all-beta, and mixed class proteins. In addition, we have utilized Baysean classification theory for discriminating two- and three-state proteins, which showed an accuracy of 90%. We have developed a web server for predicting protein folding rates and it is available at http://bioinformatics.myweb.hinet.net/foldrate.htm.     

Statistical mechanical treatment of protein conformation. III. Prediction of protein conformation based on a three-state model.

The method proposed for the evaluation of statistical weights in paper I, and the three-state model [alpha-helical (alpha), extended (epsilon), and other (c) states] formulated in paper II, have been used to develop a procedure to predict the backbone conformations of proteins, based on the concept of the predominant role played by shortrange interactions in determining protein conformation. Conformational probability profiles, in which the probabilities of formation of three consecutive alpha-helical conformations (triad) and of four consecutive extended conformations (tetrad) have been defined relative to their average values over the whole molecule, are calculated for 19 proteins, of which 16 had been used in paper I to evaluate the set of statistical weights of the 20 naturally occurring amino acids. By comparing these conformational probability profiles to experimental x-ray observations, the following results have been obtained: 80% of the alpha-helical regions and 72% of the extended conformational regions have been predicted correctly for the 19 proteins. The percentage of residues predicted correctly is in the range of 53 to 90% for the alpha-helical conformation and in the range of 63 to 88% for the extended conformation for the 19 proteins in the two-state models [alpha-helical (alpha) and other (c) states, and extended (epsilon) and other (c) states]. In the three-state model, the percentage of residues predicted correctly is in the range of 47% to 77 for 19 proteins. These results suggest that the assumption of the dominance of short-range interactions, on which the predictive scheme is based, is a reasonable one. The present predictive method is compared with that of other authors.     

Detailed unfolding and folding of gaseous ubiquitin ions characterized by electron capture dissociation

The unfolding enthalpy of the native state of ubiquitin in solution is 5 to 8 times that of its gaseous ions, as determined by electron capture dissociation (ECD) mass spectrometry. Although two-state folding occurs in solution, the three-state gaseous process proposed for this by Clemmer and co-workers based on ion mobility data is supported in general by ECD mass spectra, including relative product yields, distinct Delta H(unfolding) values between states, site-specific melting temperatures, and folding kinetics indicating a cooperative process. ECD also confirms that the 13+ ions represent separate conformers, possibly with side-chain solvated alpha-helical structures. However, the ECD data on the noncovalent bonding in the 5+ to 13+ ions, determined overall in 69 of the 75 interresidue sites, shows that thermal unfolding proceeds via a diversity of intermediates whose conformational characteristics also depend on charge site locations. As occurs with increased acidity in solution, adding 6 protons to the 5+ ions completely destroys their tertiary noncovalent bonding. However, solvation of the newly protonated sites to the backbone instead increases the stability of the secondary structure (possibly an alpha-helix) of these gaseous ions, while in solution these new sites aid denaturation by solvation in the aqueous medium. Extensive ion equilibration can lead to even more compact and diverse conformers. The three-state unfolding of gaseous ubiquitin appears to involve ensembles of individual chain conformations in a "folding funnel" of parallel reaction paths. This also provides a further caution for characterizing solution conformers from their gas-phase behavior.     

Folding kinetics of a naturally occurring helical peptide: implication of the folding speed limit of helical proteins

The folding mechanism and dynamics of a helical protein may strongly depend on how quickly its constituent alpha-helices can fold independently. Thus, our understanding of the protein folding problem may be greatly enhanced by a systematic survey of the folding rates of individual alpha-helical segments derived from their parent proteins. As a first step, we have studied the relaxation kinetics of the central helix (L9:41-74) of the ribosomal protein L9 from the bacterium Bacillus stearothermophilus , in response to a temperature-jump ( T-jump) using infrared spectroscopy. L9:41-74 has been shown to exhibit unusually high helicity in aqueous solution due to a series of side chain-side chain interactions, most of which are electrostatic in nature, while still remaining monomeric over a wide concentration range. Thus, this peptide represents an excellent model system not only for examining how the folding rate of naturally occurring helices differs from that of the widely studied alanine-based peptides, but also for estimating the folding speed limit of (small) helical proteins. Our results show that the T-jump induced relaxation rate of L9:41-74 is significantly slower than that of alanine-based peptides. For example, at 11 degrees C its relaxation time constant is about 2 micros, roughly seven times slower than that of SPE(5), an alanine-rich peptide of similar chain length. In addition, our results show that the folding rate of a truncated version of L9:41-74 is even slower. Taken together, these results suggest that individual alpha-helical segments in proteins may fold on a time scale that is significantly slower than the folding time of alanine-based peptides. Furthermore, we argue that the relaxation rate of L9:41-74 measured between 8 and 45 degrees C provides a realistic estimate of the ultimate folding rate of (small) helical proteins over this temperature range.     

Dual wavelength parametric test of two-state models for circular dichroism spectra of helical polypeptides: anomalous dichroic properties of alanine-rich peptides

A two-state helix-coil model underlies all calculations of fractional helicities FH from CD spectra of helical polypeptides. The presence of an isodichroic point near 203 nm is widely assumed to validate this model, but is shown here to provide inadequate validation for alanine-rich peptides. A parametric correlation with constant slope B between CD ellipticities at a pair of wavelengths is introduced as a more rigorous two-state test. Correlations of temperature-dependent [theta](222) vs [theta](208) values are reported for a variety of peptides. Constant slopes B are observed for literature CD data obtained from fragments of helical proteins and dimeric helical coiled coils, but parametric correlations of CD data for alanine-rich peptides consistently exhibit anomalous concave upward curvature, characterized by local slopes that are linearly temperature dependent. Low-temperature CD studies together with parametric correlations at a series of wavelengths demonstrate that the curvature anomaly is maximal at 222 nm and localized in the 215-230 nm wavelength region. Precedented structural variation of the phi, psi dihedral angles of the alpha-helix is suggested as a possible explanation. For the important case of alanine-rich peptides, experiments are proposed that may yield temperature corrections for [theta](222) and permit reliable calculations of FH from [theta](222) values.     

An 15N NMR spin relaxation dispersion study of the folding of a pair of engineered mutants of apocytochrome b562

15N relaxation dispersion NMR spectroscopy has been used to study exchange dynamics in a pair of mutants of Rd-apocyt b562, a redesigned four-helix-bundle protein. An analysis of the relaxation data over a range of temperatures establishes that exchange in both proteins is best modeled as two-state and that it derives from the folding/unfolding transition. These results are in accord with predictions based on the reaction coordinate for the folding of the protein determined from native-state hydrogen exchange data [Chu, R.; Pei, W.; Takei, J.; Bai, Y. Biochemistry 2002, 41, 7998-8003]. The kinetics and thermodynamics of the folding transition have been characterized in detail. Although only a narrow range of temperatures could be examined, it is clear that the folding rate temperature profile is distinctly non-Arrhenius for both mutants, with the folding barrier for at least one of them entropic.     

Estimation of electronic coupling in pi-stacked donor-bridge-acceptor systems: correction of the two-state model

Comparison of donor-acceptor electronic couplings calculated within two-state and three-state models suggests that the two-state treatment can provide unreliable estimates of V(da) because of neglecting the multistate effects. We show that in most cases accurate values of the electronic coupling in a pi stack, where donor and acceptor are separated by a bridging unit, can be obtained as V(da) = (E(2)-E(1))mu(12)R(da) + (2E(3)-E(1)-E(2))2mu(13)mu(23)R(da) (2), where E(1), E(2), and E(3) are adiabatic energies of the ground, charge-transfer, and bridge states, respectively, mu(ij) is the transition dipole moments between the states i and j, and R(da) is the distance between the planes of donor and acceptor. In this expression based on the generalized Mulliken-Hush approach, the first term corresponds to the coupling derived within a two-state model, whereas the second term is the superexchange correction accounting for the bridge effect. The formula is extended to bridges consisting of several subunits. The influence of the donor-acceptor energy mismatch on the excess charge distribution, adiabatic dipole and transition moments, and electronic couplings is examined. A diagnostic is developed to determine whether the two-state approach can be applied. Based on numerical results, we showed that the superexchange correction considerably improves estimates of the donor-acceptor coupling derived within a two-state approach. In most cases when the two-state scheme fails, the formula gives reliable results which are in good agreement (within 5%) with the data of the three-state generalized Mulliken-Hush model.     

Slow Unfolded-State Structuring in Acyl-CoA Binding Protein Folding Revealed by Simulation and Experiment

Protein folding is a fundamental process in biology, key to understanding many human diseases. Experimentally, proteins often appear to fold via simple two- or three-state mechanisms involving mainly native-state interactions, yet recent network models built from atomistic simulations of small proteins suggest the existence of many possible metastable states and folding pathways. We reconcile these two pictures in a combined experimental and simulation study of acyl-coenzyme A binding protein (ACBP), a two-state folder (folding time ～10 ms) exhibiting residual unfolded-state structure, and a putative early folding intermediate. Using single-molecule FRET in conjunction with side-chain mutagenesis, we first demonstrate that the denatured state of ACBP at near-zero denaturant is unusually compact and enriched in long-range structure that can be perturbed by discrete hydrophobic core mutations. We then employ ultrafast laminar-flow mixing experiments to study the folding kinetics of ACBP on the microsecond time scale. These studies, along with Trp-Cys quenching measurements of unfolded-state dynamics, suggest that unfolded-state structure forms on a surprisingly slow (～100 μs) time scale, and that sequence mutations strikingly perturb both time-resolved and equilibrium smFRET measurements in a similar way. A Markov state model (MSM) of the ACBP folding reaction, constructed from over 30 ms of molecular dynamics trajectory data, predicts a complex network of metastable stables, residual unfolded-state structure, and kinetics consistent with experiment but no well-defined intermediate preceding the main folding barrier. Taken together, these experimental and simulation results suggest that the previously characterized fast kinetic phase is not due to formation of a barrier-limited intermediate but rather to a more heterogeneous and slow acquisition of unfolded-state structure.     

Statistical mechanical treatment of protein conformation. II. A three-state model for specific-sequence copolymers of amino acids.

A one-dimensional three-state Ising model [involving alpha-helical (alpha), extended (epsilon), and coil (or other) (c) states] for specific-sequence copolymers of amino acids ahs been formulated in order to treat the conformational states of proteins. This model involves four parameters (wh,iota, vh, iota, v episilon, iota, and uc, iota), and requires a 4 X 4 matrix for generating statistical weights. Some problems in applying this model to a specific-sequence copolymer of amino acids are discussed. A nearest-neighbor approximation for treating this three-state model is also formulated; it requires a 3 X 3 matrix, in which the same four parameters appear, but (as with the 4 X 4 matrix treatment) only three parameters (wh, uh, and v epsilon) are required if relative statistical weights are used. The relationship between the present three-state model (3 X 3 matrix treatment) and models of the helix--coil transition is discussed. Then, the three-state model (3 X 3 matrix treatment) is incorporated into an earlier (Tanaka--Scheraga) model of the helix-coil transition, in which asymmetric nucleation of helical sequences is taken into account. A method for calculating molecular averages and conformational-sequence probabilities, P(iota/eta/(rho)), i.e., the probability of finding a sequence of eta residues in a specific conformational state (rho), starting at the iotath position of the chain, is described. Two alternative methods for calculating P(iota/eta/(rho)), that can be applied to a model involving any number of states, are proposed and presented; one is the direct matrix-multiplication method, and the other uses a first-order a priori probability and a conditional probability. In this paper, these calculations are performed with the nearest-neighbor model, and without the feature of asymmetric nucleation. Finally, it is indicated how the three-state model and the methods for computing P(iota/eta/(rho)) can be applied to predict protein conformation.     

Insertion and assembly of membrane proteins via simulation

Interactions of lipids are central to the folding and stability of membrane proteins. Coarse-grained molecular dynamics simulations have been used to reveal the mechanisms of self-assembly of protein/membrane and protein/detergent complexes for representatives of two classes of membrane protein, namely, glycophorin (a simple alpha-helical bundle) and OmpA (a beta-barrel). The accuracy of the coarse-grained simulations is established via comparison with the equivalent atomistic simulations of self-assembly of protein/detergent micelles. The simulation of OmpA/bilayer self-assembly reveals how a folded outer membrane protein can be inserted in a bilayer. The glycophorin/bilayer simulation supports the two-state model of membrane folding, in which transmembrane helix insertion precedes dimer self-assembly within a bilayer. The simulations also suggest that a dynamic equilibrium exists between the glycophorin helix monomer and dimer within a bilayer. The simulated glycophorin helix dimer is remarkably close in structure to that revealed by NMR. Thus, coarse-grained methods may help to define mechanisms of membrane protein (re)folding and will prove suitable for simulation of larger scale dynamic rearrangements of biological membranes.