Is protein folding hierarchic? II. Folding intermediates and transition states

Kinetic Analysis of Guanidine Hydrochloride Inactivation of β-Galactosidase in the Presence of Galactose

Inactivation of purified β-Galactosidase was done with GdnHCl in the absence and presence of varying [galactose] at 50°C and at pH 4.5. Lineweaver-Burk plots of initial velocity data, in the presence and absence of guanidine hydrochloride (GdnHCl) and galactose, were used to determine the relevant K_m and V_max values, with p-nitrophenyl β-D-galactopyranoside (pNPG) as substrate, S. Plots of ln([P]_∞ − [P]_t) against time in the presence of GdnHCl yielded the inactivation rate constant, A. Plots of A versus [S] at different galactose concentrations were straight lines that became increasingly less steep as the [galactose] increased, showing that A was dependent on [S]. Slopes and intercepts of the 1/[P]_∞ versus 1/[S] yielded k₊₀ and k'₊₀, the microscopic rate constants for the free enzyme and the enzyme-substrate complex, respectively. Plots of k₊₀ and k'₊₀ versus [galactose] showed that galactose protected the free enzyme as well as the enzyme-substrate complex (only at the lowest and highest [galactose]) against GdnHCl inactivation. In the absence of galactose, GdnHCl exhibited some degree of non-competitive inhibition. In the presence of GdnHCl, galactose exhibited competitive inhibition at the lower [galactose] of 5 mM which changed to non-competitive as the [galactose] increased. The implications of our findings are further discussed.

An instrument for fast acquisition of fluorescence decay curves at picosecond resolution designed for "double kinetics" experiments: application to fluorescence resonance excitation energy transfer study of protein folding

The information obtained by studying fluorescence decay of labeled biopolymers is a major resource for understanding the dynamics of their conformations and interactions. The lifetime of the excited states of probes attached to macromolecules is in the nanosecond time regime, and hence, a series of snapshot decay curves of such probes might - in principle - yield details of fast changes of ensembles of labeled molecules down to sub-microsecond time resolution. Hence, a major current challenge is the development of instruments for the low noise detection of fluorescence decay curves within the shortest possible time intervals. Here, we report the development of an instrument, picosecond double kinetics apparatus, that enables recording of multiple fluorescence decay curves with picosecond excitation pulses over wide spectral range during microsecond data collection for each curve. The design is based on recording and averaging multiphoton pulses of fluorescence decay using a fast 13 GHz oscilloscope during microsecond time intervals at selected time points over the course of a chemical reaction or conformational transition. We tested this instrument in a double kinetics experiment using reference probes (N-acetyl-tryptophanamide). Very low stochastic noise level was attained, and reliable multi-parameter analysis such as derivation of distance distributions from time resolved FRET (fluorescence resonance excitation energy transfer) measurements was achieved. The advantage of the pulse recording and averaging approach used here relative to double kinetics methods based on the established time correlated single photon counting method, is that in the pulse recording approach, averaging of substantially fewer kinetic experiments is sufficient for obtaining the data. This results in a major reduction in the consumption of labeled samples, which in many cases, enables the performance of important experiments that were not previously feasible.

Resolution of the effects induced by W → F substitutions on the conformation and dynamics of the amyloid-forming apomyoglobin mutant W7FW14F

Myoglobin is an alpha-helical globular protein containing two highly conserved tryptophanyl residues at positions 7 and 14 in the N-terminal region. The simultaneous substitution of the two residues increases the susceptibility of the polypeptide chain to misfold, causing amyloid aggregation under physiological condition, i.e., neutral pH and room temperature. The role played by tryptophanyl residues in driving the folding process has been investigated by examining three mutated apomyoglobins, i.e., W7F, W14F, and the amyloid-forming mutant W7FW14F, by an integrated approach based on far-ultraviolet (UV) circular dichroism (CD) analysis, fluorescence spectroscopy, and complementary proteolysis. Particular attention has been devoted to examine the conformational and dynamic properties of the equilibrium intermediate formed at pH 4.0, since it represents the early organized structure from which the native fold originates. The results show that the W → F substitutions at position 7 and 14 differently affect the structural organization of the AGH subdomain of apomyoglobin. The combined effect of the two substitutions in the double mutant impairs the formation of native-like contacts and favors interchain interactions, leading to protein aggregation and amyloid formation.

Impact of ion binding on poly-L-lysine (un)folding energy landscape and kinetics

We utilize T-jump UV resonance Raman spectroscopy (UVRR) to study the impact of ion binding on the equilibrium energy landscape and on (un)folding kinetics of poly-L-lysine (PLL). We observe that the relaxation rates of the folded conformations (including π-helix (bulge), pure α-helix, and turns) of PLL are slower than those of short alanine-based peptides. The PLL pure α-helix folding time is similar to that of short alanine-based peptides. We for the first time have directly observed that turn conformations are α-helix and π-helix (bulge) unfolding intermediates. ClO(4)(-) binding to the Lys side chain -NH(3)(+) groups and the peptide backbone slows the α-helix unfolding rate compared to that in pure water, but little impacts the folding rate, resulting in an increased α-helix stability. ClO(4)(-) binding significantly increases the PLL unfolding activation barrier but little impacts the folding barrier. Thus, the PLL folding coordinate(s) differs from the unfolding coordinate(s). The-π helix (bulge) unfolding and folding coordinates do not directly go through the α-helix energy well. Our results clearly demonstrate that PLL (un)folding is not a two-state process.

Molecular chaperones as therapeutic targets to counteract proteostasis defects

The health of cells is preserved by the levels and correct folding states of the proteome, which is generated and maintained by the proteostasis network, an integrated biological system consisting of several cytoprotective and degradative pathways. Indeed, the health conditions of the proteostasis network is a fundamental prerequisite to life as the inability to cope with the mismanagement of protein folding arising from genetic, epigenetic, and micro-environment stress appears to trigger a whole spectrum of unrelated diseases. Here we describe the potential functional role of the proteostasis network in tumor biology and in conformational diseases debating on how the signaling branches of this biological system may be manipulated to develop more efficacious and selective therapeutic strategies. We discuss the dual strategy of these processes in modulating the folding activity of molecular chaperones in order to counteract the antithetic proteostasis deficiencies occurring in cancer and loss/gain of function diseases. Finally, we provide perspectives on how to improve the outcome of these disorders by taking advantage of proteostasis modeling.

β-sheet topology prediction with high precision and recall for β and mixed α/β proteins

The prediction of the correct -sheet topology for pure and mixed proteins is a critical intermediate step toward the three dimensional protein structure prediction. The predicted beta sheet topology provides distance constraints between sequentially separated residues, which reduces the three dimensional search space for a protein structure prediction algorithm. Here, we present a novel mixed integer linear optimization based framework for the prediction of -sheet topology in and mixed proteins. The objective is to maximize the total strand-to-strand contact potential of the protein. A large number of physical constraints are applied to provide biologically meaningful topology results. The formulation permits the creation of a rank-ordered list of preferred -sheet arrangements. Finally, the generated topologies are re-ranked using a fully atomistic approach involving torsion angle dynamics and clustering. For a large, non-redundant data set of 2102 and mixed proteins with at least 3 strands taken from the PDB, the proposed approach provides the top 5 solutions with average precision and recall greater than 78%. Consistent results are obtained in the -sheet topology prediction for blind targets provided during the CASP8 and CASP9 experiments, as well as for actual and predicted secondary structures. The -sheet topology prediction algorithm, BeST, is available to the scientific community at http://selene.princeton.edu/BeST/.

Identification of helix capping and b-turn motifs from NMR chemical shifts

We present an empirical method for identification of distinct structural motifs in proteins on the basis of experimentally determined backbone and (13)C(β) chemical shifts. Elements identified include the N-terminal and C-terminal helix capping motifs and five types of β-turns: I, II, I', II' and VIII. Using a database of proteins of known structure, the NMR chemical shifts, together with the PDB-extracted amino acid preference of the helix capping and β-turn motifs are used as input data for training an artificial neural network algorithm, which outputs the statistical probability of finding each motif at any given position in the protein. The trained neural networks, contained in the MICS (motif identification from chemical shifts) program, also provide a confidence level for each of their predictions, and values ranging from ca 0.7-0.9 for the Matthews correlation coefficient of its predictions far exceed those attainable by sequence analysis. MICS is anticipated to be useful both in the conventional NMR structure determination process and for enhancing on-going efforts to determine protein structures solely on the basis of chemical shift information, where it can aid in identifying protein database fragments suitable for use in building such structures.

Exploring the energy landscapes of protein folding simulations with Bayesian computation

Nested sampling is a Bayesian sampling technique developed to explore probability distributions localized in an exponentially small area of the parameter space. The algorithm provides both posterior samples and an estimate of the evidence (marginal likelihood) of the model. The nested sampling algorithm also provides an efficient way to calculate free energies and the expectation value of thermodynamic observables at any temperature, through a simple post processing of the output. Previous applications of the algorithm have yielded large efficiency gains over other sampling techniques, including parallel tempering. In this article, we describe a parallel implementation of the nested sampling algorithm and its application to the problem of protein folding in a Gō-like force field of empirical potentials that were designed to stabilize secondary structure elements in room-temperature simulations. We demonstrate the method by conducting folding simulations on a number of small proteins that are commonly used for testing protein-folding procedures. A topological analysis of the posterior samples is performed to produce energy landscape charts, which give a high-level description of the potential energy surface for the protein folding simulations. These charts provide qualitative insights into both the folding process and the nature of the model and force field used.

Interconnection of salt-induced hydrophobic compaction and secondary structure formation depends on solution conditions: revisiting early events of protein folding at single molecule resolution

What happens in the early stage of protein folding remains an interesting unsolved problem. Rapid kinetics measurements with cytochrome c using submillisecond continuous flow mixing devices suggest simultaneous formation of a compact collapsed state and secondary structure. These data seem to indicate that collapse formation is guided by specific short and long range interactions (heteropolymer collapse). A contrasting interpretation also has been proposed, which suggests that the collapse formation is rapid, nonspecific, and a trivial solvent related compaction, which could as well be observed by a homopolymer (homopolymer collapse). We address this controversy using fluorescence correlation spectroscopy (FCS), which enables us to monitor the salt-induced compaction accompanying collapse formation and the associated time constant directly at single molecule resolution. In addition, we follow the formation of secondary structure using far UV CD. The data presented here suggest that both these models (homopolymer and heteropolymer) could be applicable depending on the solution conditions. For example, the formation of secondary structure and compact state is not simultaneous in aqueous buffer. In aqueous buffer, formation of the compact state occurs through a two-state co-operative transition following heteropolymer formalism, whereas secondary structure formation takes place gradually. In contrast, in the presence of urea, a compaction of the protein radius occurs gradually over an extended range of salt concentration following homopolymer formalism. The salt-induced compaction and the formation of secondary structure take place simultaneously in the presence of urea.

COVARIATION ANALYSIS OF LOCAL AMINO ACID SEQUENCES IN RECURRENT PROTEIN LOCAL STRUCTURES

Local structural information is supposed to be frequently encoded in local amino acid sequences. Previous research only indicated that some local structure positions have specific residue preferences in some particular local structures. However, correlated pairwise replacements for interacting residues in recurrent local structural motifs from unrelated proteins have not been studied systematically. We introduced a new method fusing statistical covariation analysis and local structure-based alignment. Systematic analysis of structure-based multiple alignments of recurrent local structures from unrelated proteins in representative subset of Protein Databank indicates that covarying residue pairs with statistical significance exist in local structural motifs, in particular β-turns and helix caps. These residue pairs are mostly linked through polar functional groups with direct or indirect hydrogen bonding. Hydrophobic interaction is also a major factor in constraining pairwise amino acid residue replacement in recurrent local structures. We also found correlated residue pairs that are not clearly linked with through-space interactions. The physical constrains underlying these covariations are less clear. Overall, covarying residue pairs with statistical significance exist in local structures from unrelated proteins. The existence of sequence covariations in local structural motifs from unrelated proteins indicates that many relics of local relations are still retained in the tertiary structures after protein folding. It supports the notion that some local structural information is encoded in local sequences and the local structural codes could play important roles in determining native state protein folding topology.

The Levinthal paradox of the interactome

The central biological question of the 21st century is: how does a viable cell emerge from the bewildering combinatorial complexity of its molecular components? Here, we estimate the combinatorics of self-assembling the protein constituents of a yeast cell, a number so vast that the functional interactome could only have emerged by iterative hierarchic assembly of its component sub-assemblies. A protein can undergo both reversible denaturation and hierarchic self-assembly spontaneously, but a functioning interactome must expend energy to achieve viability. Consequently, it is implausible that a completely "denatured" cell could be reversibly renatured spontaneously, like a protein. Instead, new cells are generated by the division of pre-existing cells, an unbroken chain of renewal tracking back through contingent conditions and evolving responses to the origin of life on the prebiotic earth. We surmise that this non-deterministic temporal continuum could not be reconstructed de novo under present conditions.

Probing molecular kinetics with Markov models: metastable states, transition pathways and spectroscopic observables

Markov (state) models (MSMs) have attracted a lot of interest recently as they (1) can probe long-term molecular kinetics based on short-time simulations, (2) offer a way to analyze great amounts of simulation data with relatively little subjectivity of the analyst, (3) provide insight into microscopic quantities such as the ensemble of transition pathways, and (4) allow simulation data to be reconciled with measurement data in a rigorous and explicit way. Here we sketch our current perspective of Markov models and explain in short their theoretical basis and assumptions. We describe transition path theory which allows the entire ensemble of protein folding pathways to be investigated and that combines naturally with Markov models. Experimental observations can be naturally linked to Markov models with the dynamical fingerprint theory, by which experimentally observable timescales can be equipped with an understanding of the structural rearrangement processes that take place at these timescales. The concepts of this paper are illustrated by a simple kinetic model of protein folding.

Search for identical octapeptides in unrelated proteins: Structural plasticity revisited

Since proteins are dynamic in nature, they can alter their local structure in response to changes in their environment factors such as temperature, pH, phosphorylation, and binding of other small molecules. These conformational changes are extremely important for the correct folding and functioning of proteins. There are also a number of diseases associated with protein conformational change such as amyloid diseases. To stimulate research into the above factors which specify one conformation over another, different theoretical models have been proposed and tested against sequence similar distant structure protein fragments. In order to simplify the computational complexity of identifying conformational changes in proteins, various local sequence search algorithms were employed and the structural plasticity in unrelated proteins was examined by various research groups. In the present work, we revisit the mechanism of structural plasticity in unrelated proteins with increased number of structures in Protein Data Bank by comparing identical octapeptides in unrelated proteins with dictionary of protein secondary structure extracted from existing experimental data. Our goal is to bring out the influence of hydrophobic residues, hydrophilic residues, flanking residues, difference in secondary structural propensities of surrounding residues, difference in phi-psi angles and local and nonlocal interactions in identical octapeptides adopting different conformations. Also we have used surrounding hydrophobicity, environment dependent interaction energy, atomic mean force potential, structural unit contacts and difference profiles models to explore the factors which cause structural plasticity. The results discussed here may provide insights into protein folding, design and function.

Dry molten globule intermediates and the mechanism of protein unfolding

New experimental results show that either gain or loss of close packing can be observed as a discrete step in protein folding or unfolding reactions. This finding poses a significant challenge to the conventional two-state model of protein folding. Results of interest involve dry molten globule (DMG) intermediates, an expanded form of the protein that lacks appreciable solvent. When an unfolding protein expands to the DMG state, side chains unlock and gain conformational entropy, while liquid-like van der Waals interactions persist. Four unrelated proteins are now known to form DMGs as the first step of unfolding, suggesting that such an intermediate may well be commonplace in both folding and unfolding. Data from the literature show that peptide amide protons are protected in the DMG, indicating that backbone structure is intact despite loss of side-chain close packing. Other complementary evidence shows that secondary structure formation provides a major source of compaction during folding. In our model, the major free-energy barrier separating unfolded from native states usually occurs during the transition between the unfolded state and the DMG. The absence of close packing at this barrier provides an explanation for why phi-values, derived from a Brønsted-Leffler plot, depend primarily on structure at the mutational site and not on specific side-chain interactions. The conventional two-state folding model breaks down when there are DMG intermediates, a realization that has major implications for future experimental work on the mechanism of protein folding.

Elucidating Peptide and Protein Structure and Dynamics: UV Resonance Raman Spectroscopy

UV resonance Raman spectroscopy (UVRR) is a powerful method that has the requisite selectivity and sensitivity to incisively monitor biomolecular structure and dynamics in solution. In this perspective, we highlight applications of UVRR for studying peptide and protein structure and the dynamics of protein and peptide folding. UVRR spectral monitors of protein secondary structure, such as the Amide III(3) band and the C(α)-H band frequencies and intensities can be used to determine Ramachandran Ψ angle distributions for peptide bonds. These incisive, quantitative glimpses into conformation can be combined with kinetic T-jump methodologies to monitor the dynamics of biomolecular conformational transitions. The resulting UVRR structural insight is impressive in that it allows differentiation of, for example, different α-helix-like states that enable differentiating π- and 3(10)- states from pure α-helices. These approaches can be used to determine the Gibbs free energy landscape of individual peptide bonds along the most important protein (un)folding coordinate. Future work will find spectral monitors that probe peptide bond activation barriers that control protein (un)folding mechanisms. In addition, UVRR studies of sidechain vibrations will probe the role of side chains in determining protein secondary, tertiary and quaternary structures.

Stabilization of proteins for storage

Following isolation and purification, it is often necessary to store proteins and peptides for extended periods of time before performing detailed biophysical, enzymatic, and structural proteomics. Therefore, it is essential that the pure target protein maintain its original biological (or functional) behavior over an extended period of storage which may range from weeks to years. Protein pharmaceuticals must remain viable following extensive shipping and storage, and they must remain devoid of all possible inactivation processes. The shelf life of a protein depends on both the intrinsic nature of the protein and the storage conditions. Proteins (especially enzymes) must be stored at an appropriate temperature and pH range and frequently in the presence of concentrated (approximately 1 M) glycerol, sucrose, or a similar substance, for the proteins to retain activity and prevent aggregation. This article discusses the major causes of protein inactivation and describes a range of measures that can be adopted to maintain the stability and solubility of proteins.

Theoretical perspectives on protein folding

Understanding how monomeric proteins fold under in vitro conditions is crucial to describing their functions in the cellular context. Significant advances in theory and experiments have resulted in a conceptual framework for describing the folding mechanisms of globular proteins. The sizes of proteins in the denatured and folded states, cooperativity of the folding transition, dispersions in the melting temperatures at the residue level, and timescales of folding are, to a large extent, determined by N, the number of residues. The intricate details of folding as a function of denaturant concentration can be predicted by using a novel coarse-grained molecular transfer model. By watching one molecule fold at a time, using single-molecule methods, investigators have established the validity of the theoretically anticipated heterogeneity in the folding routes and the N-dependent timescales for the three stages in the approach to the native state. Despite the successes of theory, of which only a few examples are documented here, we conclude that much remains to be done to solve the protein folding problem in the broadest sense.

Helix-coil energetics for helix formers and breakers reflect context and temperature: mutants of helically robust, guest-sensitive homopeptide hosts

The natural amino acids are primarily helix breakers at the low assignment temperatures characteristic of many studies, but recent genomic analyses of thermophilic proteins suggest that at high temperatures, some breakers may become strong helix formers. Moreover, the breaker/former inventory has not been previously characterized at the physiologically relevant temperature of 37 degrees C. The versatility of 13C==O NMR chemical shifts as helicity reporters allows construction of two mutant peptide series, tailored to expand the range of temperature assignments for helical propensities and derived from the core hosts tL-Ala9XxxAla9-tL and tL-AlaNva4XxxNva4Ala9-tL, Nva=norvaline. For three limiting guests Xxx, the helix former Nva and the breakers Gly and Pro, we report wXxx[T] assignments at seven temperatures from 2 to 80 degrees C, validating our reasoning and paving the way for assignment of a definitive wXxx[T] data-base.

Site-specific phosphorylation induces functionally active conformation in the intrinsically disordered N-terminal activation function (AF1) domain of the glucocorticoid receptor

Intrinsically disordered (ID) regions are disproportionately higher in cell signaling proteins and are predicted to have much larger frequency of phosphorylation sites than ordered regions, suggesting an important role in their regulatory capacity. In this study, we show that AF1, an ID activation domain of the glucocorticoid receptor (GR), adopts a functionally folded conformation due to its site-specific phosphorylation by p38 mitogen-activated protein kinase, which is involved in apoptotic and gene-inductive events initiated by the GR. Further, we show that site-specific phosphorylation-induced secondary and tertiary structure formation specifically facilitates AF1's interaction with critical coregulatory proteins and subsequently its transcriptional activity. These data demonstrate a mechanism through which ID activation domain of the steroid receptors and other similar transcription factors may adopt a functionally active conformation under physiological conditions.

Insights into protein aggregation by NMR characterization of insoluble SH3 mutants solubilized in salt-free water

Protein aggregation in vivo has been extensively associated with a large spectrum of human diseases. On the other hand, mechanistic insights into protein aggregation in vitro were incomplete due to the inability in solubilizing insoluble proteins for high-resolution biophysical investigations. However, a new avenue may be opened up by our recent discovery that previously-thought insoluble proteins can in fact be solubilized in salt-free water. Here we use this approach to study the NMR structural and dynamic properties of an insoluble SH3 mutant with a naturally-occurring insertion of Val22 at the tip of the diverging turn. The obtained results reveal: 1) regardless of whether the residue is Val, Ala, Asp or Arg, the insertion will render the first hNck2 SH3 domain to be insoluble in buffers. Nevertheless, all four mutants could be solubilized in salt-free water and appear to be largely unfolded as evident from their CD and NMR HSQC spectra. 2) Comparison of the chemical shift deviations reveals that while in V22-SH3 the second helical region is similarly populated as in the wild-type SH3 at pH 2.0, the first helical region is largely unformed. 3) In V22-SH3, many non-native medium-range NOEs manifest to define non-native helical conformations. In the meanwhile a small group of native-like long-range NOEs still persists, indicating the existence of a rudimentary native-like tertiary topology. 4) Although overall, V22-SH3 has significantly increased backbone motions on the ps-ns time scale, some regions still own restricted backbone motions as revealed by analyzing ¹⁵N relaxation data. Our study not only leads to the establishment of the first high-resolution structural and dynamic picture for an insoluble protein, but also shed more light on the molecular events for the nonhierarchical folding mechanism. Furthermore, a general mechanism is also proposed for in vivo protein aggregation triggered by the genetic mutation and posttranslational modification.

Salt dependence of an alpha-helical peptide folding energy landscapes

We used CD, UV resonance Raman spectroscopy, and molecular dynamics simulation to examine the impact of salts on the conformational equilibria and the Ramachandran Psi angle (un)folding Gibbs free energy landscape coordinate of a mainly polyalanine alpha-helical peptide, AP of sequence AAAAA(AAARA)(3)A. NaClO(4) stabilizes alpha-helical-like conformations more than does NaCl, which stabilizes more than Na(2)SO(4) at identical ionic strengths. This alpha-helix stabilization ordering is the reverse of the Hofmeister series of anions in their ability to disorder water hydrogen bonding. Much of the NaClO(4) alpha-helix stabilization results from ClO(4)(-) association with the AP terminal -NH(3)(+) groups and Arg side chains. ClO(4)(-) stabilizes 3(10)-helix conformations but destabilizes turn conformations. The decreased Cl(-) and SO(4)(2-) AP alpha-helix stabilization probably results from a decreased association with the Arg and terminal -NH(3)(+) groups. Cl(-) is expected to have a smaller binding affinity and thus stabilizes alpha-helical conformations intermediately between NaClO(4) and Na(2)SO(4). Electrostatic screening stabilizes pi-bulge conformations.

Context-independent, temperature-dependent helical propensities for amino acid residues

Assigned from data sets measured in water at 2, 25, and 60 degrees C containing (13)C=O NMR chemical shifts and [theta](222) ellipticities, helical propensities are reported for the 20 genetically coded amino acids, as well as for norvaline and norleucine. These have been introduced by chemical synthesis at central sites within length-optimized, spaced, solubilized Ala(19) hosts. The resulting polyalanine-derived, quantitative propensity sets express for each residue its temperature-dependent but context-independent tendency to forego a coil state and join a preexisting helical conformation. At 2 degrees C their rank ordering is: P << G < H < C, T, N < S < Y, F, W < V, D < K < Q < I < R, M < L < E < A; at 60 degrees C the rank becomes: H, P < G < C < R, K < T, Y, F < N, V < S < Q < W, D < I, M < E < A < L. The DeltaDeltaG values, kcal/mol, relative to alanine, for the cluster T, N, S, Y, F, W, V, D, Q, imply that at 2 degrees C all are strong breakers: DeltaDeltaG(mean) = +0.63 +/- 0.11, but at 60 degrees C their breaking tendencies are dramatically attenuated and converge toward the mean: DeltaDeltaG(mean) = +0.25 +/- 0.07. Accurate modeling of helix-rich proteins found in thermophiles, mesophiles, and organisms that flourish near 0 degrees C thus requires appropriately matched propensity sets. Comparisons are offered between the temperature-dependent propensity assignments of this study and those previously assigned by the Scheraga group; the special problems that attend propensity assignments for charged residues are illustrated by lysine guest data; and comparisons of errors in helicity assignments from shifts and ellipticity data show that the former provide superior precision and accuracy.

Mechanism of formation of the C-terminal beta-hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus. I. Importance of hydrophobic interactions in stabilization of beta-hairpin structure

We previously studied a 16-amino acid-residue fragment of the C-terminal beta-hairpin of the B3 domain (residues 46-61), [IG(46-61)] of the immunoglobulin binding protein G from Streptoccocus, and found that hydrophobic interactions and the turn region play an important role in stabilizing the structure. Based on these results, we carried out systematic structural studies of peptides derived from the sequence of IG (46-61) by systematically shortening the peptide by one residue at a time from both the C- and the N-terminus. To determine the structure and stability of two resulting 12- and 14-amino acid-residue peptides, IG(48-59) and IG(47-60), respectively, we carried out circular dichroism, NMR, and calorimetric studies of these peptides in pure water. Our results show that IG(48-59) possesses organized three-dimensional structure stabilized by hydrophobic interactions (Tyr50-Phe57 and Trp48-Val59) at T = 283 and 305 K. At T = 313 K, the structure breaks down because of increased chain entropy, but the turn region is preserved in the same position observed for the structure of the whole protein. The breakdown of structure occurs near the melting temperature of this peptide (T(m) = 310 K) measured by differential scanning calorimetry (DSC). The melting temperature of IG(47-60) determined by DSC is T(m) = 330 K and its structure is similar to that of the native beta-hairpin at all (lower) temperatures examined (283-313 K). Both of these truncated sequences are conserved in all known amino acid sequences of the B domains of the immunoglobulin binding protein G from bacteria. Thus, this study contributes to an understanding of the mechanism of folding of this whole family of proteins, and provides information about the mechanism of formation and stabilization of a beta-hairpin structural element.

Mechanism of formation of the C-terminal beta-hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus. II. Interplay of local backbone conformational dynamics and long-range hydrophobic interactions in hairpin formation

Two peptides, corresponding to the turn region of the C-terminal beta-hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus, consisting of residues 51-56 [IG(51-56)] and 50-57 [IG(50-57)], respectively, were studied by circular dichroism and NMR spectroscopy at various temperatures and by differential scanning calorimetry. Our results show that the part of the sequence corresponding to the beta-turn in the native structure (DDATKT) of the B3 domain forms bent conformations similar to those observed in the native protein. The formation of a turn is observed for both peptides in a broad range of temperatures (T = 283-323 K), which confirms the conclusion drawn from our previous studies of longer sequences from the C-terminal beta-hairpin of the B3 domain of the immunoglobulin binding protein G (16, 14, and 12 residues), that the DDATKT sequence forms a nucleation site for formation of the beta-hairpin structure of peptides corresponding to the C-terminal part of all the B domains of the immunoglobulin binding protein G. We also show and discuss the role of long-range hydrophobic interactions as well as local conformational properties of polypeptide chains in the mechanism of formation of the beta-hairpin structure.

Reconstruction and stability of secondary structure elements in the context of protein structure prediction

Efficient and accurate reconstruction of secondary structure elements in the context of protein structure prediction is the major focus of this work. We present a novel approach capable of reconstructing alpha-helices and beta-sheets in atomic detail. The method is based on Metropolis Monte Carlo simulations in a force field of empirical potentials that are designed to stabilize secondary structure elements in room-temperature simulations. Particular attention is paid to lateral side-chain interactions in beta-sheets and between the turns of alpha-helices, as well as backbone hydrogen bonding. The force constants are optimized using contrastive divergence, a novel machine learning technique, from a data set of known structures. Using this approach, we demonstrate the applicability of the framework to the problem of reconstructing the overall protein fold for a number of commonly studied small proteins, based on only predicted secondary structure and contact map. For protein G and chymotrypsin inhibitor 2, we are able to reconstruct the secondary structure elements in atomic detail and the overall protein folds with a root mean-square deviation of <10 A. For cold-shock protein and the SH3 domain, we accurately reproduce the secondary structure elements and the topology of the 5-stranded beta-sheets, but not the barrel structure. The importance of high-quality secondary structure and contact map prediction is discussed.

Amino acid insertion reveals a necessary three-helical intermediate in the folding pathway of the colicin E7 immunity protein Im7

The small (87-residue) α-helical protein Im7 (an inhibitor protein for colicin E7 that provides immunity to cells producing colicin E7) folds via a three-state mechanism involving an on-pathway intermediate. This kinetic intermediate contains three of four native helices that are oriented in a non-native manner so as to minimise exposed hydrophobic surface area at this point in folding. The short (6-residue) helix III has been shown to be unstructured in the intermediate ensemble and does not dock onto the developing hydrophobic core until after the rate-limiting transition state has been traversed. After helix III has docked, it adopts an α-helical secondary structure, and the side chains of residues within this region provide contacts that are crucial to native-state stability. In order to probe further the role of helix III in the folding mechanism of Im7, we created a variant that contains an eight-amino-acid polyalanine-like helix stabilised by a Glu-Arg salt bridge and an Asn-Pro-Gly capping motif, juxtaposed C-terminal to the natural 6-residue helix III. The effect of this insertion on the structure of the native protein and its folding mechanism were studied using NMR and ϕ-value analysis, respectively. The results reveal a robust native structure that is not perturbed by the presence of the extended helix III. Mutational analysis performed to probe the folding mechanism of the redesigned protein revealed a conserved mechanism involving the canonical three-helical intermediate. The results suggest that folding via a three-helical species stabilised by both native and non-native interactions is an essential feature of Im7 folding, independent of the helical propensity of helix III.

Analysis of protein contacts into Protein Units

Three-dimensional structures of proteins are the support of their biological functions. Their folds are maintained by inter-residue interactions which are one of the main focuses to understand the mechanisms of protein folding and stability. Furthermore, protein structures can be composed of single or multiple functional domains that can fold and function independently. Hence, dividing a protein into domains is useful for obtaining an accurate structure and function determination. In previous studies, we enlightened protein contact properties according to different definitions and developed a novel methodology named Protein Peeling. Within protein structures, Protein Peeling characterizes small successive compact units along the sequence called protein units (PUs). The cutting done by Protein Peeling maximizes the number of contacts within the PUs and minimizes the number of contacts between them. This method is so a relevant tool in the context of the protein folding research and particularly regarding the hierarchical model proposed by George Rose. Here, we accurately analyze the PUs at different levels of cutting, using a non-redundant protein databank. Distribution of PU sizes, number of PUs or their accessibility are screened to determine their common and different features. Moreover, we highlight the preferential amino acid interactions inside and between PUs. Our results show that PUs are clearly an intermediate level between secondary structures and protein structural domains.

Improving strand pairing prediction through exploring folding cooperativity

The topology of beta-sheets is defined by the pattern of hydrogen-bonded strand pairing. Therefore, predicting hydrogen bonded strand partners is a fundamental step towards predicting beta-sheet topology. At the same time, finding the correct partners is very difficult due to long range interactions involved in strand pairing. Additionally, patterns of amino acids involved, in beta-sheet formations are very general and therefore difficult to use for computational recognition of specific contacts between strands. In this work, we report a new strand pairing algorithm. To address above mentioned difficulties, our algorithm attempts to mimic elements of the folding process. Namely, in addition to ensuring that the predicted hydrogen bonded strand pairs satisfy basic global consistency constraints, it takes into account hypothetical folding pathways. Consistently with this view, introducing hydrogen bonds between a pair of strands changes the probabilities of forming hydrogen bonds between other pairs of strand. We demonstrate that this approach provides an improvement over previously proposed algorithms. We also compare the performance of this method to that of a global optimization algorithm that poses the problem as integer linear programming optimization problem and solves it using ILOG CPLEX package.

pH-induced equilibrium unfolding of apomyoglobin: substitutions at conserved Trp14 and Met131 and non-conserved Val17 positions

A number of residues in globins family are well conserved but are not directly involved in the primary oxygen-carrying function of these proteins. A possible role for these conserved, non-functional residues has been suggested in promoting a rapid and correct folding process to the native tertiary structure. To test this hypothesis, we have studied pH-induced equilibrium unfolding of mutant apomyoglobins with substitutions of the conserved residues Trp14 and Met131, which are not involved in the function of myoglobin, by various amino acids. This allowed estimating their impact on the stability of various conformational states of the proteins and selecting conditions for a folding kinetics study. The results obtained from circular dichroism, tryptophan fluorescence, and differential scanning microcalorimetry for these mutant proteins were compared with those for the wild type protein and for a mutant with the non-conserved Val17 substituted by Ala. In the native folded state, all of the mutant apoproteins have a compact globular structure, but are destabilized in comparison to the wild type protein. The pH-induced denaturation of the mutant proteins occurs through the formation of a molten globule-like intermediate similar to that of the wild type protein. Thermodynamic parameters for all of the proteins were calculated using the three state model. Stability of equilibrium intermediates at pH ~4.0 was shown to be slightly affected by the mutations. Thus, all of the above substitutions influence the stability of the native state of these proteins. The cooperativity of conformational transitions and the exposed to solvent protein surface were also changed, but not for the substitution at Val17.

Multiple routes and structural heterogeneity in protein folding

Experimental studies show that many proteins fold along sequential pathways defined by folding intermediates. An intermediate may not always be a single population of molecules but may consist of subpopulations that differ in their average structure. These subpopulations are likely to fold via independent pathways. Parallel folding and unfolding pathways appear to arise because of structural heterogeneity. For some proteins, the folding pathways can effectively switch either because different subpopulations of an intermediate get populated under different folding conditions, or because intermediates on otherwise hidden pathways get stabilized, leading to their utilization becoming discernible, or because mutations stabilize different substructures. Therefore, the same protein may fold via different pathways in different folding conditions. Multiple folding pathways make folding robust, and evolution is likely to have selected for this robustness to ensure that a protein will fold under the varying conditions prevalent in different cellular contexts.

Prediction of protein local structures and folding fragments based on building-block library

In recent years, protein structure prediction using local structure information has made great progress. In this study, a novel and effective method is developed to predict the local structure and the folding fragments of proteins. First, the proteins with known structures are split into fragments. Second, these fragments, represented by dihedrals, are clustered to produce the building blocks (BBs). Third, an efficient machine learning method is used to predict the local structures of proteins from sequence profiles. Finally, a bi-gram model, trained by an iterated algorithm, is introduced to simulate the interactions of these BBs. For test proteins, the building-block lattice is constructed, which contains all the folding fragments of the proteins. The local structures and the optimal fragments are then obtained by the dynamic programming algorithm. The experiment is performed on a subset of the PDB database with sequence identity less than 25%. The results show that the performance of the method is better than the method that uses only sequence information. When multiple paths are returned, the average classification accuracy of local structures is 72.27% and the average prediction accuracy of local structures is 67.72%, which is a significant improvement in comparison with previous studies. The method can predict not only the local structures but also the folding fragments of proteins. This work is helpful for the ab initio protein structure prediction and especially, the understanding of the folding process of proteins.

Extrinsic fluorescent dyes as tools for protein characterization

Noncovalent, extrinsic fluorescent dyes are applied in various fields of protein analysis, e.g. to characterize folding intermediates, measure surface hydrophobicity, and detect aggregation or fibrillation. The main underlying mechanisms, which explain the fluorescence properties of many extrinsic dyes, are solvent relaxation processes and (twisted) intramolecular charge transfer reactions, which are affected by the environment and by interactions of the dyes with proteins. In recent time, the use of extrinsic fluorescent dyes such as ANS, Bis-ANS, Nile Red, Thioflavin T and others has increased, because of their versatility, sensitivity and suitability for high-throughput screening. The intention of this review is to give an overview of available extrinsic dyes, explain their spectral properties, and show illustrative examples of their various applications in protein characterization.

Thermal unfolding of Escherichia coli trigger factor studied by ultra-sensitive differential scanning calorimetry

Temperature-induced unfolding of Escherichia coli trigger factor (TF) and its domain truncation mutants, NM and MC, were studied by ultra-sensitive differential scanning calorimetry (UC-DSC). Detailed thermodynamic analysis showed that thermal induced unfolding of TF and MC involves population of dimeric intermediates. In contrast, the thermal unfolding of the NM mutant involves population of only monomeric states. Covalent cross-linking experiments confirmed the presence of dimeric intermediates during thermal unfolding of TF and MC. These data not only suggest that the dimeric form of TF is extremely resistant to thermal unfolding, but also provide further evidence that the C-terminal domain of TF plays a vital role in forming and stabilizing the dimeric structure of the TF molecule. Since TF is the first molecular chaperone that nascent polypeptides encounter in eubacteria, the stable dimeric intermediates of TF populated during thermal denaturation might be important in responding to stress damage to the cell, such as heat shock.