Effect In Protein Folding Research Articles

Effects of Topology and Sequence in Protein Folding Linked via Conformational Fluctuations

Experiments have compared the folding of proteins with different amino acid sequences but the same basic structure, or fold. Results indicate that folding is robust to sequence variations for proteins with some nonlocal folds, such as all-β, whereas the folding of more local, all-α proteins typically exhibits a stronger sequence dependence. Here, we use a coarse-grained model to systematically study how variations in sequence perturb the folding energy landscapes of three model sequences with 3α, 4β + α, and β-barrel folds, respectively. These three proteins exhibit folding features in line with experiments, including expected rank order in the cooperativity of the folding transition and stability-dependent shifts in the location of the free-energy barrier to folding. Using a generalized-ensemble simulation approach, we determine the thermodynamics of around 2000 sequence variants representing all possible hydrophobic or polar single- and double-point mutations. From an analysis of the subset of stability-neutral mutations, we find that folding is perturbed in a topology-dependent manner, with the β-barrel protein being the most robust. Our analysis shows, in particular, that the magnitude of mutational perturbations of the transition state is controlled in part by the size or “width” of the underlying conformational ensemble. This result suggests that the mutational robustness of the folding of the β-barrel protein is underpinned by its conformationally restricted transition state ensemble, revealing a link between sequence and topological effects in protein folding.

Nonnative contact effects in protein folding.

A comprehensive understanding of protein folding includes the knowledge of the formation of individual secondary structures, tertiary structure, and the effects of non-native contacts on these folding events. The measurement of these microscopic events has been posing challenges for experiment and molecular simulation. In this work, we performed enhanced sampling MD simulations for three proteins (NTL9, NuG2b, and CspA) and analyzed minimum free energy paths on multi-dimensional free energy landscapes to explore the underlying folding mechanisms. Consistencies can be seen between the present simulations and the existing experiments as well as other MD simulations. Quantitative analysis reveals the nucleation-condensation folding mechanism indicating the concurrent build-up of secondary and tertiary structures for the three proteins and gives the detailed formation sequence of individual native secondary structure elements. More importantly, nonnative contacts are generally observed among the proteins, creating a nonnative environment to affect the folding of individual secondary structure elements. A general tendency is that the secondary structure element(s) where the maximal nonnative contacts are observed have the largest formation free-energy barrier(s), corresponding to the rate-limiting step(s) of the folding for proteins that follow the nucleation-condensation mechanism. In summary, while native contacts determine the folding mechanism and pathway, non-native contacts play an important role in determining the protein folding thermodynamics by influencing the free energies of individual secondary structure element formation.

Shielding effect in protein folding

One of the most important interactions responsible for protein folding and stability are hydrogen bonds between peptide groups. There is a constant competition between the water molecules and peptide groups in a hydrogen bond formation. Also side-chains take part in this process by reducing hydration of peptide group (shielding effect) that promotes the protein folding. In this paper, a new approach to take into account a shielding effect is presented. A modification of the energy function is derived and incorporated into the UNited RESidue (UNRES) force field. Canonical Molecular Dynamics and Replica Exchange Molecular Dynamics with UNRES force field is applied to study the influence of this effect on protein structure, folding kinetics and free energy landscapes. The results of test calculations suggest that even small contribution of this effect into energy function changes force field behavior as well as speeds up the folding process significantly.

Towards a structural biology of the hydrophobic effect in protein folding.

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

Accounting for a mirror-image conformation as a subtle effect in protein folding.

By using local (free-energy profiles along the amino acid sequence and (13)C(α) chemical shifts) and global (principal component) analyses to examine the molecular dynamics of protein-folding trajectories, generated with the coarse-grained united-residue force field, for the B domain of staphylococcal protein A, we are able to (i) provide the main reason for formation of the mirror-image conformation of this protein, namely, a slow formation of the second loop and part of the third helix (Asp29-Asn35), caused by the presence of multiple local conformational states in this portion of the protein; (ii) show that formation of the mirror-image topology is a subtle effect resulting from local interactions; (iii) provide a mechanism for how protein A overcomes the barrier between the metastable mirror-image state and the native state; and (iv) offer a plausible reason to explain why protein A does not remain in the metastable mirror-image state even though the mirror-image and native conformations are at least energetically compatible.

The Rise and Fall of the Hydrophobic Effect in Protein Folding and Protein-Protein Association, and Molecular Recognition

In the beginning everything was explained in Biochemistry in terms of hydrogen-bonds (HB). Then, the devastating blow, known as the HB-inventory argument came; hydrogen bonding with water molecules compete with intramolecular hydrogen-bonds. As a result, the HBs paradigm fell from grace. The void created was immediately filled by Kauzmann’s idea of hydrophobic (HφO) effect which reigned supreme in biochemical literature for over 50 years (1960-2010). Cracks in the HB-inventory argument on one hand, and doubts about the adequacy of Kauzmann’s model for the HφO effect, have led to a comeback of the HBs, along with a host of new hydrophilic (HφI) effects. The HφO effects lost much of its power - which it never really had - in explaining protein folding and protein-protein association. Instead, the more powerful and richer repertoire of HφI effects took over the reins. The interactions also offered simple and straightforward answers to the problems of protein folding, and protein-protein association.

Ion Specificity in α-Helical Folding Kinetics

The influence of the salts KCl, NaCl, and NaI at molar concentrations on the α-helical folding kinetics of the alanine-based oligopeptide Ace-AEAAAKEAAAKA-Nme is investigated by means of (explicit-water) molecular dynamics simulations and a diffusional analysis. The mean first passage times for folding and unfolding are found to be highly salt-specific. In particular, the folding times increase about 1 order of magnitude for the sodium salts. The drastic slowing can be traced to long-lived, compact configurations of the partially folded peptide, in which sodium ions are tightly bound by several carbonyl and carboxylate groups. This multiple trapping leads to a nonexponential residence time distribution of the cations in the first solvation shell of the peptide. The analysis of α-helical folding in the framework of diffusion in a reduced (one-dimensional) free energy landscape further shows that the salt not only specifically modifies equilibrium properties but also induces kinetic barriers due to individual ion binding. In the sodium salts, for instance, the peptide's configurational mobility (or "diffusivity") can decrease about 1 order of magnitude. This study demonstrates the highly specific action of ions and highlights the intimate coupling of intramolecular friction and solvent effects in protein folding.

Alpha helical crossovers favor right‐handed supersecondary structures by kinetic trapping: The phone cord effect in protein folding

The remarkable predominance of right-handedness in beta-alpha-beta helical crossovers has been previously explained in terms of thermodynamic stability and kinetic accessibility, but a different kinetic trapping mechanism may also play a role. If the beta-sheet contacts are made before the crossover helix is fully formed, and if the backbone angles of the folding helix follows the energetic pathway of least resistance, then the helix would impart a torque on the ends of the two strands. Such a torque would tear apart a left-handed conformation but hold together a right-handed one. Right-handed helical crossovers predominate even in all-alpha proteins, where previous explanations based on the preferred twist of the beta sheet do not apply. Using simple molecular simulations, we can reproduce the right-handed preference in beta-alpha-beta units, without imposing specific beta-strand geometry. The new kinetic trapping mechanism is dubbed the "phone cord effect" because it is reminiscent of the way a helical phone cord forms superhelices to relieve torsional stress. Kinetic trapping explains the presence of a right-handed superhelical preference in alpha helical crossovers and provides a possible folding mechanism for knotted proteins.

Hydrodynamic Selection of the Kinetic Pathway of a Polymer Coil-Globule Transition

Recently, the role of hydrodynamic interactions in the selection of a kinetic pathway for phase transitions has attracted considerable attention. Here we study this problem numerically by taking as an example a coil-globule transition of a single polymer, which is a prototype model of protein folding. When a swollen polymer collapses into a globule state, hydrodynamic interactions accelerate the transition. We find, on the other hand, that when a rather compact polymer collapses into the same final state, hydrodynamic interactions decelerate the transition due to a slow squeezing process of the solvent. We reveal that the degree of the initial enhancement of anisotropy of the polymer configuration determines whether hydrodynamic interactions accelerate or decelerate the collapsing dynamics. We also discuss the possible relevance of squeezing flow effects in protein folding.

ESBRI: A web server for evaluating salt bridges in proteins

Salt bridges can play important roles in protein structure and function and have stabilizing and destabilizing effects in protein folding. ESBRI is a software available as web tool which analyses the salt bridges in a protein structure, starting from the atomic coordinates. In the case of protein complexes, the salt bridges between protein chains can be evaluated, as well as those among specific charged amino acids and the different protein subunits, in order to obtain useful information regard the protein-protein interaction. The service is available at the URL: http://bioinformatica.isa.cnr.it/ESBRI/

Simulating Virtual Protein Cα Traces with Applications

We propose a simple procedure for generating virtual protein C(alpha) traces. One of the key ingredients of our method, to build a three-dimensional structure from a random sequence of amino acids, is to work directly on torsional angles of the chain which we sample from a von Mises distribution. With simple modeling of the hydrophobic effect in protein folding, the procedure produces compact and globular structures. Some characteristics of real proteins (i.e., compactness and globularity) are well mimicked by this procedure. These virtual traces are used to assess algorithms for matching protein structures or functional sites.

An analytical study of the interplay between geometrical and energetic effects in protein folding

Analytical studies have several advantages for an understanding of the mechanisms of protein folding such as the interplay between geometrical and energetic effects. In this paper, we introduce a Gaussian filament with a C(alpha) structure-based (Go) potential as a new theoretical scheme based on a Hamiltonian approach. This model takes into account geometrical information in a realistic fashion without the need of phenomenological descriptions. In order to make this model more appropriate for comparison with protein folding simulations and experiments, we introduce a many-body interaction into the potential term to enhance cooperativity. We apply our new analytical model to a beta-hairpin-type peptide and compare our results with a molecular dynamics simulation of a structure-based model.

Ionic-Strength-Dependent Effects in Protein Folding: Analysis of Rate Equilibrium Free-Energy Relationships and Their Interpretation

The traditional approach to studying protein folding involves applying a perturbation, usually denaturant or mutation, and determining the effect upon the free energy of folding, DeltaG0, and the activation free energy, DeltaG(not equal). Data collected as a function of the perturbation can be used to construct rate equilibrium free-energy relationships, which report on the development of interactions in the transition state for folding. We examine the use of the ionic-strength-dependent rate equilibrium free-energy relationship in protein folding using the N-terminal domain of L9, a small alpha-beta protein, as a model system. Folding is two-state for the range of ionic strength examined, 0.045-1.52 M. The plot of DeltaG(not equal) versus DeltaG0 is linear (r2= 0.918), with a slope equal to 0.45. The relatively low value of the slope indicates that the ionic-strength-dependent interactions are modestly developed in the transition state. The slope is, however, greater than that of a plot of DeltaG(not equal) versus DeltaG0 constructed by varying pH, thus demonstrating directly that ionic-strength-dependent studies probe more than simple electrostatic interactions. Potential transition movement was probed by analysis of the denaturant, ionic strength cross-interaction parameters. The values are small but nonzero and positive, suggesting a small shift of the transition state toward the native state as the protein is destabilized, i.e., Hammond behavior. The complications that arise in the interpretation of ionic-strength-dependent rate equilibrium free-energy relationships are discussed, and it is concluded that the ionic-strength-dependent studies do not provide a reliable indicator of the role of electrostatic interactions. Complications include incomplete screening of electrostatic interactions, specific ion binding, Hofmeister effects, and the potential presence of electrostatic interactions in the denatured state ensemble.

Deducing the Energetic Cost of Protein Folding in Zinc Finger Proteins Using Designed Metallopeptides

Zinc finger transcription factors represent the largest single class of metalloproteins in the human genome. Binding of Zn(II) to their canonical Cys4, Cys3His1, or Cys2His2 sites results in metal-induced protein folding events required to achieve their proper structure for biological activity. The thermodynamic contribution of Zn(II) in each of these coordination spheres toward protein folding is poorly understood because of the coupled nature of the metal-ligand and protein-protein interactions. Using an unstructured peptide scaffold, GGG, we have employed fluorimetry, potentiometry, and calorimetry to determine the thermodynamics of Zn(II) binding to the Cys4, Cys3His1, and Cys2His2 ligand sets with minimal interference from protein folding effects. The data show that Zn(II) complexation is entropy driven and modulated by proton release. The formation constants for Zn(II)-GGG with a Cys4, Cys3His1, or Cys2His2 site are 5.6 x 10(16), 1.5 x 10(15), or 2.5 x 10(13) M(-1), respectively. Thus, the Zn(II)-Cys4, Zn(II)-Cys3His1, and Zn(II)-Cys2His2 interactions can provide up to 22.8, 20.7, and 18.3 kcal/mol, respectively, in driving force for protein stabilization, folding, and/or assembly at pH values above the ligand pKa values. While the contributions from the three coordination motifs differ by 4.5 kcal/mol in Zn(II) affinity at pH 9.0, they are equivalent at physiological pH, DeltaG = -16.8 kcal/mol or a Ka = 2.0 x 10(12) M(-1). Calorimetric data show that this is due to proton-based enthalpy-entropy compensation between the favorable entropic term from proton release and the unfavorable enthalpic term due to thiol deprotonation. Since protein folding effects have been minimized in the GGG scaffold, these peptides possess nearly the tightest Zn(II) affinities possible for their coordination motifs. The Zn(II) affinities in each coordination motif are compared between the GGG scaffold and natural zinc finger proteins to determine the free energy required to fold the latter. Several proteins have identical Zn(II) affinities to GGG. That is, little, if any, of their Zn(II) binding energy is required to fold the protein, whereas some have affinities weakened by up to 5.7 kcal/mol; i.e., the Zn(II) binding energy is being used to fold the protein.

Solvation and desolvation effects in protein folding: native flexibility, kinetic cooperativity and enthalpic barriers under isostability conditions

As different parts of a protein chain approach one another during folding, they are expected to encounter desolvation barriers before optimal packing is achieved. This impediment originates from the water molecule's finite size, which entails a net energetic cost for water exclusion when the formation of compensating close intraprotein contacts is not yet complete. Based on recent advances, we extend our exploration of these microscopic elementary desolvation barriers' roles in the emergence of generic properties of protein folding. Using continuum Gō-like Cα chain models of chymotrypsin inhibitor 2 (CI2) and barnase as examples, we underscore that elementary desolvation barriers between a protein's constituent groups can significantly reduce native conformational fluctuations relative to model predictions that neglected these barriers. An increasing height of elementary desolvation barriers leads to thermodynamically more cooperative folding/unfolding transitions (i.e., higher overall empirical folding barriers) and higher degrees of kinetic cooperativity as manifested by more linear rate-stability relationships under constant temperature. Applying a spatially non-uniform thermodynamic parametrization we recently introduced for the pairwise Cα potentials of mean force, the present barnase model further illustrates that desolvation is a probable physical underpinning for the experimentally observed high intrinsic enthalpic folding barrier under isostability conditions.

Modeling the Interplay between Geometrical and Energetic Effects in Protein Folding

A theoretical framework is constructed with the aid of a free-energy functional method that is capable of describing the interplay between geometrical and energetic effects on protein folding. In this paper, we generalize a free-energy functional model based on polymer theory to make it more appropriate for comparison with protein folding simulations and experiments. This generalization is made by introducing cooperativity into the configurational entropy and the internal energy. Modifications to configurational entropy enable the model to account for the loop-loop interactions, a contribution neglected in the original model. Modifications to the internal energy introduce many-body corrections, which are needed to establish quantitative contact to simulations as well as experimental observations. To demonstrate the efficiency of the modified analytical model, we compare our results with C(alpha) structure-based (Go) model simulations of chymotrypsin inhibitor II and the SH3 domain of src.

Heat capacity effects in protein folding and ligand binding: a re-evaluation of the role of water in biomolecular thermodynamics

Large “anomalous” heat capacity (Δ C p) effects are a common feature of the thermodynamics of biomolecular interactions in aqueous solution and, as a result of the improved facility for direct calorimetric measurements, there is a growing body of experimental data for such effects in protein folding, protein–protein and protein–ligand interactions. Conventionally such heat capacity effects have been ascribed to hydrophobic interactions, and there are some remarkably convincing demonstrations of the usefulness of this concept. Nonetheless, there is also increasing evidence that hydrophobic interactions are not the only possible source of such effects. Here we re-evaluate the possible contributions of other interactions to the heat capacity changes to be expected for cooperative biomolecular folding and binding processes, with particular reference to the role of hydrogen bonding and solvent water interactions. Simple models based on the hydrogen-bonding propensity of water as a function of temperature give quantitative estimates of Δ C p that compare well with experimental observations for both protein folding and ligand binding. The thermodynamic contribution of bound waters in protein complexes is also estimated. The prediction from simple lattice models is that trapping of water in a complex should give more exothermic binding (ΔΔ H−6 to −12 kJ mol −1) with lower entropy (ΔΔ S 0≈−11 J mol −1 K −1) and more negative Δ C p (by about −75 J mol −1 K −1) per water molecule. More generally, it is clear that significant Δ C p effects are to be expected for any macromolecular process involving a multiplicity of cooperative weak interactions of whatever kind.

Sequence-space selection of cooperative model proteins

We introduce a cooperative lattice model to account for specific many-body effects in protein folding. Cooperativity arises as local solvent environments are determined by the large-scale organization of the chain. In contrast to an uncorrelated lattice, a selection pressure in sequence space arises: foldable sequences endowed with sharp equilibrium conformation ensembles are scarce, indicating that a few sequences are able to exclude the solvent where it matters to warrant the integrity of the structures they form in solution.

Hammond Behavior versus Ground State Effects in Protein Folding: Evidence for Narrow Free Energy Barriers and Residual Structure in Unfolded States

Apparent transition state movement upon mutation or changes in solvent conditions is frequently observed in protein folding and is often interpreted in terms of Hammond behavior. This led to the conclusion that barrier regions in protein folding are broad maxima on the free energy landscape. Here, we use the concept of self-interaction and cross-interaction parameters to test experimental data of 21 well-characterized proteins for Hammond behavior. This allows us to characterize the origin of transition state movements along different reaction coordinates. Only one of the 21 proteins shows a small but coherent transition state movement in agreement with the Hammond postulate. In most proteins the structure of the transition state is insensitive to changes in protein stability. The apparent change in the position of the transition state upon mutation, which is frequently observed in φ-value analysis, is in most cases due to ground-state effects caused by structural changes in the unfolded state. This argues for significant residual structure in unfolded polypeptide chains of many proteins. Disruption of these residual interactions by mutation often leads to decreased folding rates, which implies that these interactions are still present in the transition state. The failure to detect Hammond behavior shows that the free energy barriers encountered by a folding polypeptide chain are generally rather narrow and robust maxima for all experimentally explorable reaction coordinates.

The Hydrophobic Effect: A New Insight from Cold Denaturation and a Two-State Water Structure

Herein we provide a new insight into the hydrophobic effect in protein folding. Our proposition explains the molecular basis of cold denaturation, and of intermediate states in heat and their absence in cold denaturation. The exposure of non-polar surface reduces the entropy and enthalpy of the system, at low and at high temperatures. At low temperatures the favorable reduction in enthalpy overcomes the unfavorable reduction in entropy, leading to cold denaturation. At high temperatures, folding/unfolding is a two-step process: in the first, the entropy gain leads to hydrophobic collapse, in the second, the reduction in enthalpy due to protein-protein interactions leads to the native state. The different entropy and enthalpy contributions to the Gibbs energy change at each step at high, and at low, temperatures can be conveniently explained by a two-state model of the water structure. The model provides a clear view of the dominant factors in protein folding and stability. Consequently, it appears to provide a microscopic view of the hydrophobic effect and is consistently linked to macroscopic thermodynamic parameters.