Does alpha-helix folding necessarily provide an energy source for the protein-lipid binding?

Lipid-induced alpha-helix folding, which occurs in many lipid surface-binding proteins and peptides such as apolipoproteins and synucleins, has been proposed to provide an energy source for protein-lipid interactions. We propose that in a system comprised of a phospholipid surface and a small polypeptide that is unfolded in solution and binds reversibly to lipid surface, helical folding involves expenditure of free energy as compared to a similar polypeptide that is alpha-helical in solution. This is a consequence of the entropic cost of helix folding that is illustrated in a simple thermodynamic model and exemplifies the general "key-into-lock" paradigm of protein-ligand binding. Even though this simple model does not explicitly address the protein-induced lipid re-arrangement and may not directly apply to large proteins that undergo significant tertiary structural changes upon lipid binding, it suggests that the notion of helix folding as an energy source for lipid binding should be treated with caution.     

Exploring structures in protein folding funnels with free energy functionals: the transition state ensemble

We use free energy functionals that account for the partial ordering of residues in the transition state ensemble to characterize the free energy surfaces for fast folding proteins. We concentrate on chymotrypsin inhibitor and lambda-repressor. We show how the explicit cooperativity that can arise from many body forces, such as side-chain ordering or hydrophobic surface burial, determines the crossover from folding with a large delocalized nucleus and the specific small classical nucleus of the type envisioned in nucleation growth scenarios. We compare the structural correlations present in the transition state ensemble obtained from free energy functionals with those inferred from experiment using extrathermodynamic free energy relations for folding time obtained via protein engineering kinetics experiments. We also use the free energy functionals to examine both the size of barriers and multidimensional representations of the free energy profiles in order to address the question of appropriate reaction coordinates for folding.     

A heteropolymer model study for the mechanism of protein folding

A heteropolymer model of randomly self-interacting chains in two dimensions is studied with numerical simulations in order to elucidate the folding mechanism of protein. We find that the model occasionally shows folding propensity depending on the sequence of random numbers given to the chain. We study the thermodynamic and kinematic roles in the folding mechanism by grouping the local energy minima found in the simulations into clusters according to the similarity of their conformations. It is suggested that the local minima to which some heteropolymers show a folding tendency are always the lowest energy states of the energy spectrum within a cluster, though which cluster is selected depends on the sequence. For the eight random sequences we study, we find that the energy gap between the ground state and excited states is little correlated with folding or nonfolding. We rather find that folding propensities are correlated with the global structure of the average energy surface, implying a dominant kinetic role in the folding mechanism, although thermal factors cannot be ignored as the mechanism of choosing the ground state within a cluster of states connected by small deformations. We suggest that a hierarchical cluster structure plays an important role in selecting a unique folded state out of the huge number of local minima of heteropolymers. © 1997 John Wiley & Sons, Inc.     

Hierarchic finite level energy landscape model: to describe the refolding kinetics of phosphoglycerate kinase

One of the most intriguing predictions of energy landscape models is the existence of non-exponential protein folding kinetics caused by hierarchical structures in the landscapes. Here we provide the strongest evidence so far of such hierarchy and determine the time constants and weights of the kinetic components of the suggested hierarchic energy landscape. To our knowledge, the idea of hierarchical folding energy barriers has never been tested over such a broad timescale. Refolding of yeast phosphoglycerate kinase was initiated from the guanidine-unfolded state by stopped-flow or manual mixing and monitored by tryptophan fluorescence from 1 ms to 15 min. The strategy to build a model that describes folding of yeast phosphoglycerate kinase was to start from the simplest paradigm and modify it stepwise to the necessary minimal extent after repeated comparisons with the experiments. We made no a priori assumptions about the folding landscape. The result was a hierarchic finite level landscape model that quantitatively describes the refolding of yeast phosphoglycerate kinase from 1 ms to 15 min. The early steps of the folding process happen in the upper region of the landscape, where the surface has a hierarchic structure. This leads to stretched kinetics in the early phase of the folding. The lower region of the energy landscape is dominated by a trap that reflects the accumulation of molten globule intermediate state. From this intermediate, the protein can reach the global energy minimum corresponding to the native state through a cross-barrier folding step.     

A new method to measure cortical growth in the developing brain

Folding of the cerebral cortex is a critical phase of brain development in higher mammals but the biomechanics of folding remain incompletely understood. During folding, the growth of the cortical surface is heterogeneous and anisotropic. We developed and applied a new technique to measure spatial and directional variations in surface growth from longitudinal magnetic resonance imaging (MRI) studies of a single animal or human subject. MRI provides high resolution 3D image volumes of the brain at different stages of development. Surface representations of the cerebral cortex are obtained by segmentation of these volumes. Estimation of local surface growth between two times requires establishment of a point-to-point correspondence ("registration") between surfaces measured at those times. Here we present a novel approach for the registration of two surfaces in which an energy function is minimized by solving a partial differential equation on a spherical surface. The energy function includes a strain-energy term due to distortion and an "error energy" term due to mismatch between surface features. This algorithm, implemented with the finite element method, brings surface features into approximate alignment while minimizing deformation in regions without explicit matching criteria. The method was validated by application to three simulated test cases and applied to characterize growth of the ferret cortex during folding. Cortical surfaces were created from MRI data acquired in vivo at 14 days, 21 days, and 28 days of life. Deformation gradient and Lagrangian strain tensors describe the kinematics of growth over this interval. These quantitative results illuminate the spatial, temporal, and directional patterns of growth during cortical folding.     

The Energy Landscape,Folding Pathways and the Kinetics of a Knotted Protein

The folding pathway and rate coefficients of the folding of a knotted protein are calculated for a potential energy function with minimal energetic frustration. A kinetic transition network is constructed using the discrete path sampling approach, and the resulting potential energy surface is visualized by constructing disconnectivity graphs. Owing to topological constraints, the low-lying portion of the landscape consists of three distinct regions, corresponding to the native knotted state and to configurations where either the N or C terminus is not yet folded into the knot. The fastest folding pathways from denatured states exhibit early formation of the N terminus portion of the knot and a rate-determining step where the C terminus is incorporated. The low-lying minima with the N terminus knotted and the C terminus free therefore constitute an off-pathway intermediate for this model. The insertion of both the N and C termini into the knot occurs late in the folding process, creating large energy barriers that are the rate limiting steps in the folding process. When compared to other protein folding proteins of a similar length, this system folds over six orders of magnitude more slowly.     

Protein energy landscape exploration with structure-based models

Exploring the multi-dimensional energy landscape of a large protein in detail is a computational challenge. Such investigations may include analysis of multiple folding pathways, rate constants for important conformational transitions, locating intermediate states populated during folding, estimating energetic and entropic barriers that separate populated basins, and visualising a high-dimensional surface. The complexity of the landscape can be simplified through coarse-grained structure-based models (SBMs). These widely used coarse-grained representations of proteins provide a minimalist approximation to the free energy landscape, which subsumes the folding behaviour of many single-domain proteins. Here we describe the combination of SBMs with discrete path sampling (DPS), and show how this approach can provide details of the landscape and folding pathways. Combining SBMs and DPS provides an efficient framework for sampling the protein free energy landscape and for calculating various kinetic and thermodynamic quantities.     

154 Protein folding and binding funnels: a common driving force and a common mechanism

Under the free energy landscape theory, both the protein-folding and protein–ligand binding processes are driven by the decrease in total Gibbs free energy of the protein-solvent or protein–ligand-solvent system, which involves the non-complementary changes between the entropy and enthalpy, ultimately leading to a global free energy minimization of these thermodynamic systems (Ji & Liu, 2011; Liu et al., 2012; Yang, Ji & Liu, 2012). In the case of protein folding, the lowering of the system free energy coupled with the gradual reduction in conformational degree of freedom of the folding intermediates determines that the shape of the free energy landscape for protein folding must be funnel-like (Dill & Chan, 1997), rather than non-funneled shapes (Ben-Naim, 2012). In the funnel-like free energy landscape, protein folding can be viewed as going down the hill via multiple parallel routes from a vast majority of individual non-native states on surface outside the funnel to the native states located around the bottom of the funnel. The first stage of folding, i.e. the rapid hydrophobic collapse process, is driven by the solvent entropy maximization. Concretely, the water molecules squeeze and sequestrate the hydrophobic amino acid side chains within the interior of the folding intermediates while exposing the polar and electrostatically charged side chains on the intermediate surface so as to minimize the solvent-accessible surface area of the solute and thus, the minimal contacts between the folding intermediates and the water molecules. This will maximize the entropy of the solvent, thus contributing substantially to lowering of the system free energy due to an absolute advantage of the solvent in both quantity and mass (Yang, Ji & Liu, 2012). The resulting molten globule states (Ohgushi & Wada, 1983), within which a few transient secondary structural components and tertiary contacts have been formed but many native contacts or close residue–residue interactions has yet to form, need to be further sculptured into the native states. This is a relatively slow “bottleneck” process because the competitive interactions between protein residues within the folding intermediates and between residues and water molecules may repeat many rounds to accumulate a large enough number of stable noncovalent bonds capable of counteracting the conformational entropy loss of the intermediates, thus putting this bottleneck stage under the enthalpy control (i.e. negative enthalpy change), contributing further to the lowering of the system free energy. Although the protein–ligand association occurs around the rugged bottom of the free energy landscape, the exclusion of water from the binding interfaces and the formation of noncovalent bonds between the two partners can still lower the system free energy. In conjunction with the loss of the rotational and translational degrees of freedom of the two partners as well as the loss of the conformational entropy of the protein, these processes could merge, downwards expand, and further narrow the free energy wells within which the protein–ligand binding process takes place, thereby making them look like a funnel, which we term the binding funnel. In this funnel, the free energy downhill process follows a similar paradigm to the protein-folding process. For example, if the initial collisions/contacts occur between the properly complementary interfaces of the protein and ligand, a large amount of water molecules (which usually form a water network around the solute surface) will be displaced to suit the need for maximizing the solvent entropy. This process is similar to that of the hydrophobic collapse during protein folding, resulting in a loosely associated protein–ligand complex that needs also to be further adapted into a tight complex, i.e. the second step which is mainly driven by the negative enthalpy change through intermolecular competitive interactions to gradually accumulate the noncovalent bonds and ultimately, to stabilize the complex at a tightly bound state. Taken together, we conclude that whether in the protein-folding or in the protein–ligand binding process, both the entropy-driven first step and the enthalpy-driven second step contribute to the lowering of the system free energy, resulting in the funnel-like folding or binding free energy landscape.     

Single-molecule unfolding force distributions reveal a funnel-shaped energy landscape

The protein folding process is described as diffusion on a high-dimensional energy landscape. Experimental data showing details of the underlying energy surface are essential to understanding folding. So far in single-molecule mechanical unfolding experiments a simplified model assuming a force-independent transition state has been used to extract such information. Here we show that this so-called Bell model, although fitting well to force velocity data, fails to reproduce full unfolding force distributions. We show that by applying Kramers' diffusion model, we were able to reconstruct a detailed funnel-like curvature of the underlying energy landscape and establish full agreement with the data. We demonstrate that obtaining spatially resolved details of the unfolding energy landscape from mechanical single-molecule protein unfolding experiments requires models that go beyond the Bell model.     

Statistical mechanics of protein folding by exhaustive enumeration

It is hard to construct theories for the folding of globular proteins because they are large and complicated molecules having enormous numbers of nonnative conformations and having native states that are complicated to describe. Statistical mechanical theories of protein folding are constructed around major simplifying assumptions about the energy as a function of conformation and/or simplifications of the representation of the polypeptide chain, such as one point per residue on a cubic lattice. It is not clear how the results of these theories are affected by their various simplifications. Here we take a very different simplification approach where the chain is accurately represented and the energy of each conformation is calculated by a not unreasonable empirical function. However, the set of amino acid sequences and allowed conformations is so restricted that it becomes computationally feasible to examine them all. Hence we are able to calculate melting curves for thermal denaturation as well as the detailed kinetic pathway of refolding. Such calculations are based on a novel representation of the conformations as points in an abstract 12-dimensional Euclidean conformation space. Fast folding sequences have relatively high melting temperatures, native structures with relatively low energies, small kinetic barriers between local minima, and relatively many conformations in the global energy minimum's watershed. In contrast to other folding theories, these models show no necessary relationship between fast folding and an overall funnel shape to the energy surface, or a large energy gap between the native and the lowest nonnative structure, or the depth of the native energy minimum compared to the roughness of the energy landscape. Proteins 32:425–437, 1998. © 1998 Wiley-Liss, Inc.     

Use of a quantitative structure-property relationship to design larger model proteins that fold rapidly

A quantitative structure-property relationship (QSPR) was used to design model protein sequences that fold repeatedly and relatively rapidly to stable target structures. The specific model was a 125-residue heteropolymer chain subject to Monte Carlo dynamics on a simple cubic lattice. The QSPR was derived from an analysis of a database of 200 sequences by a statistical method that uses a genetic algorithm to select the sequence attributes that are most important for folding and a neural network to determine the corresponding functional dependence of folding ability on the chosen attributes. The QSPR depends on the number of anti-parallel sheet contacts, the energy gap between the native state and quasi-continuous part of the spectrum and the total energy of the contacts between surface residues. Two Monte Carlo procedures were used in series to optimize both the target structures and the sequences. We generated 20 fully optimized sequences and 60 partially optimized control sequences and tested each for its ability to fold in dynamic MC simulations. Although sequences in which either the number of anti-parallel sheet contacts or the energy of the surface residues is non-optimal are capable of folding almost as well as fully optimized ones, sequences in which only the energy gap is optimized fold markedly more slowly. Implications of the results for the design of proteins are discussed.     

Shape of the free energy barriers for protein folding probed by multiple perturbation analysis

The characterization of the free energy barriers has been a major goal in studies on the mechanism of protein folding. Testing the effect of mutations or denaturants on protein folding reactions revealed that transition state movement is rare, suggesting that folding barriers are robust and narrow maxima on the free energy landscape. Here we demonstrate that the application of multiple perturbations allows the observation of small transition state movements that escape detection in single perturbation experiments. We used tendamistat as a model protein to test the broadness of the free energy barriers. Tendamistat folds over two consecutive transition states and through a high-energy intermediate. Measuring the combined effect of temperature and denaturant on the position of the transition state in the wild-type protein and in several mutants revealed that the early transition state shows significant transition state movement. Its accessible surface area state becomes more native-like with destabilization of the native state by temperature. To the same extent, the entropy of the early transition state becomes more native-like with increasing denaturant concentration, in accordance with Hammond behavior. The position of the late transition state, in contrast, is much less sensitive to the applied perturbations. These results suggest that the barriers in protein folding become increasingly narrow as the folding polypeptide chain approaches the native state.     

Folding is not required for bilayer insertion: replica exchange simulations of an alpha-helical peptide with an explicit lipid bilayer

We implement the replica exchange molecular dynamics algorithm to study the interactions of a model peptide (WALP-16) with an explicitly represented DPPC membrane bilayer. We observe the spontaneous, unbiased insertion of WALP-16 into the DPPC bilayer and its folding into an alpha-helix with a transbilayer orientation. The free energy surface suggests that the insertion of the peptide into the DPPC bilayer precedes secondary structure formation. Although the peptide has some propensity to form a partially helical structure in the interfacial region of the DPPC/water system, this state is not a productive intermediate but rather an off-pathway trap for WALP-16 insertion. Equilibrium simulations show that the observed insertion/folding pathway mirrors the potential of mean force (PMF). Calculation of the enthalpic and entropic contributions to this PMF show that the surface bound conformation of WALP-16 is significantly lower in energy than other conformations, and that the insertion of WALP-16 into the bilayer without regular secondary structure is enthalpically unfavorable by 5-10 kcal/mol/residue. The observed insertion/folding pathway disagrees with the dominant conceptual model, which is that a surface-bound helix is an obligatory intermediate for the insertion of alpha-helical peptides into lipid bilayers. In our simulations, the observed insertion/folding pathway is favored because of a large (>100 kcal/mol) increase in system entropy that occurs when the unstructured WALP-16 peptide enters the lipid bilayer interior. The insertion/folding pathway that is lowest in free energy depends sensitively on the near cancellation of large enthalpic and entropic terms. This suggests the possibility that intrinsic membrane peptides may have a diversity of insertion/folding behaviors depending on the exact system of peptide and lipid under consideration.     

Probe-dependent and nonexponential relaxation kinetics: unreliable signatures of downhill protein folding

The theoretical suggestion that some proteins may encounter no significant free energy barriers during their folding raises an important question: What experimental signature does this downhill folding produce? Several authors have argued that nonexponential (and especially stretched exponential) or probe‐dependent kinetics represent useful experimental signatures of a downhill free energy surface. Here we examine more closely the connection between unusual kinetics and downhill energy surfaces. Simulation of diffusive relaxation dynamics on a variety of generically downhill, two‐dimensional free energy surfaces shows that these surfaces do not necessarily produce either probe‐dependent or significantly nonexponential kinetics. Conversely, we find that two‐dimensional surfaces with significant (>3k_B T) energy barriers can readily give rise to nonexponential and probe‐dependent kinetics. These results show that downhill folding does not constitute a necessary or sufficient condition for nonexponential and/or probe‐dependent folding kinetics. One cannot easily prove or disprove that a protein folds downhill simply by studying its relaxation kinetics. Proteins 2007. © 2007 Wiley‐Liss, Inc.     

Use of a potential of mean force to analyze free energy contributions in protein folding.

A method for calculation of the free energy of residues as a function of residue burial is proposed. The method is based on the potential of mean force, with a reaction coordinate expressed by residue burial. Residue burials are calculated from high-resolution protein structures. The largest individual contributions to the free energy of a residue are found to be due to the hydrophobic interactions of the nonpolar atoms, interactions of the main chain polar atoms, and interactions of the charged groups of residues Arg and Lys. The contribution to the free energy of folding due to the uncharged side chain polar atoms is small. The contribution to the free energy of folding due to the main chain polar atoms is favorable for partially buried residues and less favorable or unfavorable for fully buried residues. Comparison of the accessible surface areas of proteins and model spheres shows that proteins deviate considerably from a spherical shape and that the deviations increase with the size of a protein. The implications of these results for protein folding are also discussed.     

The C-terminal Tails of the Bacterial Chaperonin GroEL Stimulate Protein Folding by Directly Altering the Conformation of a Substrate Protein

Many essential cellular proteins fold only with the assistance of chaperonin machines like the GroEL-GroES system of Escherichia coli. However, the mechanistic details of assisted protein folding by GroEL-GroES remain the subject of ongoing debate. We previously demonstrated that GroEL-GroES enhances the productive folding of a kinetically trapped substrate protein through unfolding, where both binding energy and the energy of ATP hydrolysis are used to disrupt the inhibitory misfolded states. Here, we show that the intrinsically disordered yet highly conserved C-terminal sequence of the GroEL subunits directly contributes to substrate protein unfolding. Interactions between the C terminus and the non-native substrate protein alter the binding position of the substrate protein on the GroEL apical surface. The C-terminal tails also impact the conformational state of the substrate protein during capture and encapsulation on the GroEL ring. Importantly, removal of the C termini results in slower overall folding, reducing the fraction of the substrate protein that commits quickly to a productive folding pathway and slowing several kinetically distinct folding transitions that occur inside the GroEL-GroES cavity. The conserved C-terminal tails of GroEL are thus important for protein folding from the beginning to the end of the chaperonin reaction cycle.     

The folding energy landscape of the dimerization domain of Escherichia coli Trp repressor: a joint experimental and theoretical investigation

Enhanced structural insights into the folding energy landscape of the N-terminal dimerization domain of Escherichia coli tryptophan repressor, [2-66]2 TR, were obtained from a combined experimental and theoretical analysis of its equilibrium folding reaction. Previous studies have shown that the three intertwined helices in [2-66]2 TR are sufficient to drive the formation of a stable dimer for the full-length protein, [2-107]2 TR. The monomeric and dimeric folding intermediates that appear during the folding reactions of [2-66]2 TR have counterparts in the folding mechanism of the full-length protein. The equilibrium unfolding energy surface on which the folding and dimerization reactions occur for [2-66]2 TR was examined with a combination of native-state hydrogen exchange analysis, pepsin digestion and matrix-assisted laser/desorption mass spectrometry performed at several concentrations of protein and denaturant. Peptides corresponding to all three helices in [2-66]2 TR show multi-layered protection patterns consistent with the relative stabilities of the dimeric and monomeric folding intermediates. The observation of protection exceeding that offered by the dimeric intermediate in segments from all three helices implies that a segment-swapping mechanism may be operative in the monomeric intermediate. Protection greater than that expected from the global stability for a single amide hydrogen in a peptide from the C-helix possibly and another from the A-helix may reflect non-random structure, possibly a precursor for segment swapping, in the urea-denatured state. Native topology-based model simulations that correspond to a funnel energy landscape capture both the monomeric and dimeric intermediates suggested by the HX MS data and provide a rationale for the progressive acquisition of secondary structure in their conformational ensembles.