Backbone dynamics, fast folding, and secondary structure formation in helical proteins and peptides

We discuss the construction of a simple, off-lattice model protein with a comparatively detailed representation of the protein backbone, and use it to address some general aspects of the folding kinetics of a small helical protein and two peptide fragments. The model makes use of an associative memory hamiltonian to smoothly interpolate between the limits of a native contact only, or Go, potential and a statistical pair potential derived from a database of known structures. We have observed qualitatively different behavior in these two limits. In the Go limit, we see apparently barrier-less folding. As we increase the roughness of the model energy landscape, we can observe the emergence of the characteristic activated temperature dependence previously seen in lattice studies and analytical theories. We are also able to study the dependence of the folding kinetics on local interactions such as hydrogen bonds, and we discuss the implications of these results for the formation of secondary structure at intermediate stages of the folding reaction.     

Deciphering the timescales and mechanisms of protein folding using minimal off-lattice models.

Considerable insights into the mechanisms and timescales of protein folding have been obtained from detailed studies of minimal off-lattice models. These models are coarse-grained representations of polypeptide chains. Many novel predictions of the mechanisms and timescales of the folding of proteins have been made using simulations of off-lattice models. The concepts derived from these simulations have been used to analyze the recent experiments and simulations of proteins and peptides.     

Sequence-based study of two related proteins with different folding behaviors

Z(SPA-1) is an engineered protein that binds to its parent, the three-helix-bundle Z domain of staphylococcal protein A. Uncomplexed Z(SPA-1) shows a reduced helix content and a melting behavior that is less cooperative, compared with the wild-type Z domain. Here we show that the difference in folding behavior between these two sequences can be partly understood in terms of an off-lattice model with 5-6 atoms per amino acid and a minimalistic potential, in which folding is driven by backbone hydrogen bonding and effective hydrophobic attraction.     

The effects of nonnative interactions on protein folding rates: theory and simulation

Proteins are minimally frustrated polymers. However, for realistic protein models, nonnative interactions must be taken into account. In this paper, we analyze the effect of nonnative interactions on the folding rate and on the folding free energy barrier. We present an analytic theory to account for the modification on the free energy landscape upon introduction of nonnative contacts, added as a perturbation to the strong native interactions driving folding. Our theory predicts a rate-enhancement regime at fixed temperature, under the introduction of weak, nonnative interactions. We have thoroughly tested this theoretical prediction with simulations of a coarse-grained protein model, by using an off-lattice C(alpha)model of the src-SH3 domain. The strong agreement between results from simulations and theory confirm the nontrivial result that a relatively small amount of nonnative interaction energy can actually assist the folding to the native structure.     

Coarse-grained models of protein folding: toy models or predictive tools?

Coarse-grained models are emerging as a practical alternative to all-atom simulations for the characterization of protein folding mechanisms over long time scales. While a decade ago minimalist toy models were mainly designed to test general hypotheses on the principles regulating protein folding, the latest coarse-grained models are increasingly realistic and can be used to characterize quantitatively the detailed folding mechanism of specific proteins. The ability of such models to reproduce the essential features of folding dynamics suggests that each single atomic degree of freedom is not by itself particularly relevant to folding and supports a statistical mechanical approach to characterize folding transitions. When combined with more refined models and with experimental studies, the systematic investigation of protein systems and complexes using coarse-grained models can advance our theoretical understanding of the actual organizing principles that emerge from the complex network of interactions among protein atomic constituents.     

Dependence of folding dynamics and structural stability on the location of a hydrophobic pair in beta-hairpins

We study the dependence of folding time, nucleation site, and stability of a model beta-hairpin on the location of a cross-strand hydrophobic pair, using a coarse-grained off-lattice model with the aid of Monte Carlo simulations. Our simulations have produced 6500 independent folding trajectories dynamically, forming the basis for extensive statistical analysis. Four folding pathways, zipping-out, middle-out, zipping-in, and reptation, have been closely monitored and discussed in all seven sequences studied. A hydrophobic pair placed near the beta-turn or in the middle section effectively speed up folding; a hydrophobic pair placed close to the terminal ends or next to the beta-turn encourages stability of the entire chain.     

Solvent effects on the conformational transition of a model polyalanine peptide

We have investigated the folding of polyalanine by combining discontinuous molecular dynamics simulation with our newly developed off-lattice intermediate-resolution protein model. The thermodynamics of a system containing a single Ac-KA(14)K-NH(2) molecule has been explored by using the replica exchange simulation method to map out the conformational transitions as a function of temperature. We have also explored the influence of solvent type on the folding process by varying the relative strength of the side-chain's hydrophobic interactions and backbone hydrogen bonding interactions. The peptide in our simulations tends to mimic real polyalanine in that it can exist in three distinct structural states: alpha-helix, beta-structures (including beta-hairpin and beta-sheet-like structures), and random coil, depending upon the solvent conditions. At low values of the hydrophobic interaction strength between nonpolar side-chains, the polyalanine peptide undergoes a relatively sharp transition between an alpha-helical conformation at low temperatures and a random-coil conformation at high temperatures. As the hydrophobic interaction strength increases, this transition shifts to higher temperatures. Increasing the hydrophobic interaction strength even further induces a second transition to a beta-hairpin, resulting in an alpha-helical conformation at low temperatures, a beta-hairpin at intermediate temperatures, and a random coil at high temperatures. At very high values of the hydrophobic interaction strength, polyalanines become beta-hairpins and beta-sheet-like structures at low temperatures and random coils at high temperatures. This study of the folding of a single polyalanine-based peptide sets the stage for a study of polyalanine aggregation in a forthcoming paper.     

Folding of bundles of alpha-helices in solution, membranes, and adsorbed overlayers

We propose a coarse-grained lattice model for Monte Carlo simulations of folding of proteins consisting of several alpha-helices. A chain representing a protein is considered to contain A and B monomers forming relatively stiff A subchains, mimicking helices, and flexible B links between these subchains, respectively. Using this model, we simulate (1) folding of four-helix proteins in solution; (2) folding of membrane proteins containing one, two, or four helices; and (3) refolding of four-helix proteins adsorbed at the liquid-solid interface. For these cases, we show typical scenarios of protein folding and refolding and study the dependence of the folding time on the chain length. Combining the latter results with those already available in the literature, we discuss the relative rates of folding of proteins belonging to different classes.     

Orientational potentials extracted from protein structures improve native fold recognition

We develop coarse-grained, distance- and orientation-dependent statistical potentials from the growing protein structural databases. For protein structural classes (alpha, beta, and alpha/beta), a substantial number of backbone-backbone and backbone-side-chain contacts stabilize the native folds. By taking into account the importance of backbone interactions with a virtual backbone interaction center as the 21st anisotropic site, we construct a 21 x 21 interaction scheme. The new potentials are studied using spherical harmonics analysis (SHA) and a smooth, continuous version is constructed using spherical harmonic synthesis (SHS). Our approach has the following advantages: (1) The smooth, continuous form of the resulting potentials is more realistic and presents significant advantages for computational simulations, and (2) with SHS, the potential values can be computed efficiently for arbitrary coordinates, requiring only the knowledge of a few spherical harmonic coefficients. The performance of the new orientation-dependent potentials was tested using a standard database of decoy structures. The results show that the ability of the new orientation-dependent potentials to recognize native protein folds from a set of decoy structures is strongly enhanced by the inclusion of anisotropic backbone interaction centers. The anisotropic potentials can be used to develop realistic coarse-grained simulations of proteins, with direct applications to protein design, folding, and aggregation.     

Protein folding pathways and state transitions described by classical equations of motion of an elastic network model

Protein topology defined by the matrix of residue contacts has proved to be a fruitful basis for the study of protein dynamics. The widely implemented coarse-grained elastic network model of backbone fluctuations has been used to describe crystallographic temperature factors, allosteric couplings, and some aspects of the folding pathway. In the present study, we develop a model of protein dynamics based on the classical equations of motion of a damped network model (DNM) that describes the folding path from a completely unfolded state to the native conformation through a single-well potential derived purely from the native conformation. The kinetic energy gained through the collapse of the protein chain is dissipated through a friction term in the equations of motion that models the water bath. This approach is completely general and sufficiently fast that it can be applied to large proteins. Folding pathways for various proteins of different classes are described and shown to correlate with experimental observations and molecular dynamics and Monte Carlo simulations. Allosteric transitions between alternative protein structures are also modeled within the DNM through an asymmetric double-well potential.     

Folding of the protein domain hbSBD

The folding of the alpha-helix domain hbSBD of the mammalian mitochondrial branched-chain alpha-ketoacid dehydrogenase complex is studied by the circular dichroism technique in absence of urea. Thermal denaturation is used to evaluate various thermodynamic parameters defining the equilibrium unfolding, which is well described by the two-state model with the folding temperature T(F) = 317.8 +/- 1.95 K and the enthalpy change DeltaH(G) = 19.67 +/- 2.67 kcal/mol. The folding is also studied numerically using the off-lattice coarse-grained Go model and the Langevin dynamics. The obtained results, including the population of the native basin, the free-energy landscape as a function of the number of native contacts, and the folding kinetics, also suggest that the hbSBD domain is a two-state folder. These results are consistent with the biological function of hbSBD in branched-chain alpha-ketoacid dehydrogenase.     

Emergence of highly designable protein-backbone conformations in an off-lattice model

Despite the variety of protein sizes, shapes, and backbone configurations found in nature, the design of novel protein folds remains an open problem. Within simple lattice models it has been shown that all structures are not equally suitable for design. Rather, certain structures are distinguished by unusually high designability: the number of amino acid sequences for which they represent the unique lowest energy state; sequences associated with such structures possess both robustness to mutation and thermodynamic stability. Here we report that highly designable backbone conformations also emerge in a realistic off-lattice model. The highly designable conformations of a chain of 23 amino acids are identified and found to be remarkably insensitive to model parameters. Although some of these conformations correspond closely to known natural protein folds, such as the zinc finger and the helix-turn-helix motifs, others do not resemble known folds and may be candidates for novel fold design.     

Folding Trp-cage to NMR resolution native structure using a coarse-grained protein model

We develop a coarse-grained protein model with a simplified amino acid interaction potential. Using this model, we perform discrete molecular dynamics folding simulations of a small 20-residue protein--Trp-cage--from a fully extended conformation. We demonstrate the ability of the Trp-cage model to consistently reach conformations within 2-angstroms backbone root-mean-square distance from the corresponding NMR structures. The minimum root-mean-square distance of Trp-cage conformations in simulations can be <1 angstroms. Our findings suggest that, at least in the case of Trp-cage, a detailed all-atom protein model with a molecular mechanics force field is not necessary to reach the native state of a protein. Our results also suggest that the success of folding Trp-cage in our simulations and in the reported all-atom molecular mechanics simulation studies may be mainly due to the special stabilizing features specific to this miniprotein.     

Toward minimalist models of larger proteins: a ubiquitin-like protein

Our recently developed off-lattice bead model capable of simulating protein structures with mixed alpha/beta content has been extended to model the folding of a ubiquitin-like protein and provides a means for examining the more complex kinetics involved in the folding of larger proteins. Using trajectories generated from constant-temperature Langevin dynamics simulations and sampling with the multiple multi-histogram method over five-order parameters, we are able to characterize the free energy landscape for folding and find evidence for folding through compact intermediates. Our model reproduces the observation that the C-terminus loop structure in ubiquitin is the last to fold in the folding process and most likely plays a spectator role in the folding kinetics. The possibility of a productive metastable intermediate along the folding pathway consisting of collapsed states with no secondary structure, and of intermediates or transition structures involving secondary structural elements occurring early in the sequence, is also supported by our model. The kinetics of folding remain multi-exponential below the folding temperature, with glass-like kinetics appearing at T/T(f) approximately 0.86. This new physicochemical model, designed to be predictive, helps validate the value of modeling protein folding at this level of detail for genomic-scale studies, and motivates further studies of other protein topologies and the impact of more complex energy functions, such as the addition of solvation forces.     

Structural transitions of confined model proteins: molecular dynamics simulation and experimental validation

Proteins fold in a confined space not only in vivo, i.e., folding assisted by molecular chaperons and chaperonins in a crowded cellular medium, but also in vitro as in production of recombinant proteins. Despite extensive work on protein folding in bulk, little is known about how and to what extent the thermodynamics and kinetics of protein folding are altered by confinement. In this work, we use a Gō-like off-lattice model to investigate the folding and stability of an all beta-sheet protein in spherical cages of different sizes and surface hydrophobicity. We find whereas extreme confinement inhibits correct folding, a hydrophilic cage stabilizes the protein due to restriction of the unfolded configurations. In a hydrophobic cage, however, strong attraction from the cage surface destabilizes the confined protein because of competition between self-aggregation and adsorption of hydrophobic residues. We show that the kinetics of protein collapse and folding is strongly correlated with both the cage size and the surface hydrophobicity. It is demonstrated that a cage of moderate size and hydrophobicity optimizes both the folding yield and kinetics of structural transitions. To support the simulation results, we have also investigated the refolding of hen-egg lysozyme in the presence of cetyltrimethylammoniumbromide (CTAB) surfactants that provide an effective confinement of the proteins by micellization. The influence of the surfactant hydrophobicity on the structural and biological activity of the protein is determined with circular dichroism spectrum, fluorescence emission spectrum, and biological activity assay. It is shown that, as predicted by coarse-grained simulations, CTAB micelles facilitate the collapse of denatured lysozyme, whereas the addition of beta-cyclodextrin-grafted-PNIPAAm, a weakly hydrophobic stripper, dissociates CTAB micelles and promotes the conformational rearrangement and thereby gives an improved recovery of lysozyme activity.