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

Thermodynamics of protein folding using a modified Wako-Saitô-Muñoz-Eaton model

Herein, we propose a modified version of the Wako-Saitô-Muñoz-Eaton (WSME) model. The proposed model introduces an empirical temperature parameter for the hypothetical structural units (i.e., foldons) in proteins to include site-dependent thermodynamic behavior. The thermodynamics for both our proposed model and the original WSME model were investigated. For a system with beta-hairpin topology, a mathematical treatment (contact-pair treatment) to facilitate the calculation of its partition function was developed. The results show that the proposed model provides better insight into the site-dependent thermodynamic behavior of the system, compared with the original WSME model. From this site-dependent point of view, the relationship between probe-dependent experimental results and model's thermodynamic predictions can be explained. The model allows for suggesting a general principle to identify foldon behavior. We also find that the backbone hydrogen bonds may play a role of structural constraints in modulating the cooperative system. Thus, our study may contribute to the understanding of the fundamental principles for the thermodynamics of protein folding.

Replica exchange molecular dynamics of the thermodynamics of fibril growth of Alzheimer's Aβ42 peptide

The growth of amyloid fibrils is studied by replica exchange molecular dynamics in an implicit solvent. Our data indicate that extremely long simulation times (at least a few hundred ns) are necessary to study the thermodynamics of fibril elongation in detail. However some aspects of the aggregation process are already accessible on the time scales available in the present study. A peak in the specific heat indicates a docking temperature of T(dock) ≈ 320 K. Irreversible locking requires lower temperatures with the locking temperature estimated as T(lock) ≈ 280 K. In our simulation the fibril grows from both sides with the C-terminal of the incoming monomer attaching to the C-terminal of the peptides in the fibril forming a β-sheet on the fibril edge. Our simulation indicates that the C-terminal is crucial for aggregation.

Internal and environmental effects on folding and dimerization of the Alzheimer's β amyloid peptide

Amyloid deposits are a hallmark of many diseases. In the case of Alzheimer's disease a turn between 21Ala and 30Ala, stabilized by a salt bridge between 22Glu/23Asp and 28Lys, may nucleate folding and aggregation of the Aβ peptide. In the present paper we test this hypothesis by studying how salt bridge and turn formation vary with intrinsic and environmental changes, and how these changes effect folding and aggregation of the Aβ peptide.

Folding mechanism of β-hairpin trpzip2: heterogeneity, transition state and folding pathways

We review the studies on the folding mechanism of the β-hairpin tryptophan zipper 2 (trpzip2) and present some additional computational results to refine the picture of folding heterogeneity and pathways. We show that trpzip2 can have a two-state or a multi-state folding pattern, depending on whether it folds within the native basin or through local state basins on the high-dimensional free energy surface; Trpzip2 can fold along different pathways according to the packing order of tryptophan pairs. We also point out some important problems related to the folding mechanism of trpzip2 that still need clarification, e.g., a wide distribution of the computed conformations for the transition state ensemble.

Replica exchange molecular dynamics simulations of coarse-grained proteins in implicit solvent

Current approaches aimed at determining the free energy surface of all-atom medium-size proteins in explicit solvent are slow and are not sufficient to converge to equilibrium properties. To ensure a proper sampling of the configurational space, it is preferable to use reduced representations such as implicit solvent and/or coarse-grained protein models, which are much lighter computationally. Each model must be verified, however, to ensure that it can recover experimental structures and thermodynamics. Here we test the coarse-grained implicit solvent OPEP model with replica exchange molecular dynamics (REMD) on six peptides ranging in length from 10 to 28 residues: two alanine-based peptides, the second beta-hairpin from protein G, the Trp-cage and zinc-finger motif, and a dimer of a coiled coil peptide. We show that REMD-OPEP recovers the proper thermodynamics of the systems studied, with accurate structural description of the beta-hairpin and Trp-cage peptides (within 1-2 A from experiments). The light computational burden of REMD-OPEP, which enables us to generate many hundred nanoseconds at each temperature and fully assess convergence to equilibrium ensemble, opens the door to the determination of the free energy surface of larger proteins and assemblies.

Folding mechanisms of individual beta-hairpins in a Go model of Pin1 WW domain by all-atom molecular dynamics simulations

This paper examines the folding mechanism of an individual beta-hairpin in the presence of other hairpins by using an off-lattice model of a small triple-stranded antiparallel beta-sheet protein, Pin1 WW domain. The turn zipper model and the hydrophobic collapse model originally developed for a single beta-hairpin in literature is confirmed to be useful in describing beta-hairpins in model Pin1 WW domain. We find that the mechanism for folding a specific hairpin is independent of whether it folds first or second, but the formation process are significantly dependent on temperature. More specifically, beta1-beta2 hairpin folds via the turn zipper model at a low temperature and the hydrophobic collapse model at a high temperature, while the folding of beta2-beta3 hairpin follows the turn zipper model at both temperatures. The change in folding mechanisms is interpreted by the interplay between contact stability (enthalpy) and loop lengths (entropy), the effect of which is temperature dependent.

The complex folding pathways of protein A suggest a multiple-funnelled energy landscape

Folding proteins into their native states requires the formation of both secondary and tertiary structures. Many questions remain, however, as to whether these form into a precise order, and various pictures have been proposed that place the emphasis on the first or the second level of structure in describing folding. One of the favorite test models for studying this question is the B domain of protein A, which has been characterized by numerous experiments and simulations. Using the activation-relaxation technique coupled with a generic energy model (optimized potential for efficient peptide structure prediction), we generate more than 50 folding trajectories for this 60-residue protein. While the folding pathways to the native state are fully consistent with the funnel-like description of the free energy landscape, we find a wide range of mechanisms in which secondary and tertiary structures form in various orders. Our nonbiased simulations also reveal the presence of a significant number of non-native beta and alpha conformations both on and off pathway, including the visit, for a non-negligible fraction of trajectories, of fully ordered structures resembling the native state of nonhomologous proteins.

Minimum model for the alpha-helix-beta-hairpin transition in proteins

We present a minimal model for proteins, which is able to capture the structural conversion between the alpha-helix and beta-hairpin. In most regimes of the parameter space, the model produces a stable structure at a low temperature; in a few limited regimes of the parameter space, the model displays an beta-hairpin transition as the physical conditions vary. These variations include a perturbation on hydrogen bonding propensity at the middle of the modeled chain, or the change of the hydrophobicity of a designated pair along the chain. Using Monte Carlo simulations, we demonstrate the structural conversion by means of state diagrams, heat capacity maps, and free energy maps.

Aggregating the amyloid Abeta(11-25) peptide into a four-stranded beta-sheet structure

We present a detailed analysis of the structural properties of one monomer of Abeta(11-25) as well as of the aggregation mechanisms for four chains of Abeta(11-25) using the activation-relaxation technique coupled with a generic energy potential. Starting from a random distribution of these four chains, we find that the system assembles rapidly into a random globular state that evolves into three- and four-stranded antiparallel beta-sheets. The aggregation process is considerably accelerated by the presence of preformed dimers. We also find that the reptation mechanism already identified in shorter peptides plays a significant role here in allowing the structure to reorganize without having to fully dissociate.

Coarse-grained protein molecular dynamics simulations

A limiting factor in biological science is the time-scale gap between experimental and computational trajectories. At this point, all-atom explicit solvent molecular dynamics (MD) are clearly too expensive to explore long-range protein motions and extract accurate thermodynamics of proteins in isolated or multimeric forms. To reach the appropriate time scale, we must then resort to coarse graining. Here we couple the coarse-grained OPEP model, which has already been used with activated methods, to MD simulations. Two test cases are studied: the stability of three proteins around their experimental structures and the aggregation mechanisms of the Alzheimer's Abeta16-22 peptides. We find that coarse-grained isolated proteins are stable at room temperature within 50 ns time scale. Based on two 220 ns trajectories starting from disordered chains, we find that four Abeta16-22 peptides can form a three-stranded beta sheet. We also demonstrate that the reptation move of one chain over the others, first observed using the activation-relaxation technique, is a kinetically important mechanism during aggregation. These results show that MD-OPEP is a particularly appropriate tool to study qualitatively the dynamics of long biological processes and the thermodynamics of molecular assemblies.

Evaluation of coarse grained models for hydrogen bonds in proteins

Backbone hydrogen bonds contribute very importantly to the stability of proteins and therefore they must be appropriately represented in protein folding simulations. Simple models are frequently used in theoretical approaches to this process, but their simplifications are often confronted with the need to be true to the physics of the interactions. Here we study the effects of different levels of coarse graining in the modeling of backbone hydrogen bonds. We study three different models taken from the bibliography in a twofold fashion. First, we calculate the hydrogen bonds in 2gb1, an (alpha + beta)-protein, and see how different backbone representations and potentials can mimic the effects of real hydrogen bonds both in helices and sheets. Second, we use an evolutionary method for protein fragment assembly to locate the global energy minimum for a set of small beta-proteins with these models. This way, we assess the effects of coarse graining in hydrogen bonding models and show what can be expected from them when used in simulation experiments.

Molecular dynamics simulations of folding processes of a β-hairpin in an implicit solvent

Computer simulations of beta-hairpin folding are relatively difficult, especially those based on the explicit water model. This greatly limits the complete analysis and understanding of their folding mechanisms. In this paper, we use the generalized Born/solvent accessible implicit solvent model to simulate the folding processes of a nine-residue beta-hairpin. We find that the beta-hairpin can fold into its native structure very easily, even using the traditional molecular dynamics method. This allows us to extract 21 complete folding events and investigate the folding process sufficiently. Our results show that there exist four most stable states on the free energy landscape of the short peptide, one native state and three intermediates. We find that two of the non-native stable states have almost the same potential energy as the native state but with lower entropy. This suggests that the native state can be stabilized entropically. Furthermore, we find that the folding processes of this peptide have common features: to fold into its native state, the peptide undergoes a continuous collapsing-extending-recollapsing process to adjust the positions of the side chains in order to form the native middle inter-strand hydrogen bonds. The formations of these bonds are the key step of the folding process. Once these bonds are formed, the peptide can fold into the native state quickly.

Model systems for β-hairpins and β-sheets

beta-Sheets and alpha-helices are the two principal secondary structures in proteins. However, our understanding of beta-sheet structure lags behind that of alpha-helices, largely because, until recently, there was no model system to study the beta-sheet secondary structure in isolation. With the development of well-folded beta-hairpins, this is changing rapidly. Recent advances include: increased understanding of the relative contributions of turn, strand and sidechain interactions to beta-hairpin and beta-sheet stability, with the role of aromatic residues as a common subtheme; experimental and theoretical kinetic and thermodynamic studies of beta-hairpin and beta-sheet folding; de novo protein design, including all-beta structures, mixed alpha/beta motifs and switchable systems; and the creation of functional beta-hairpins.