Make Illumina part of your DNA  Sign up for PNAS Online eTocs
Link: Info for AuthorsLink: Editorial BoardLink: AboutLink: SubscribeLink: AdvertiseLink: ContactLink: Sitemap Link: PNAS Home
Proceedings of the National Academy of Sciences
Link: Current Issue "" Link: Archives "" Link: Online Submission ""  Link: Advanced Search

Published online on February 20, 2001, 10.1073/pnas.041611998
PNAS | February 27, 2001 | vol. 98 | no. 5 | 2188-2193


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a colleague
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My File Cabinet
Right arrow Download to citation manager
Right arrow Request Copyright Permission
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via ISI Web of Science (55)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Levy, Y.
Right arrow Articles by Becker, O. M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Levy, Y.
Right arrow Articles by Becker, O. M.
Social Bookmarking
 Add to CiteULike   Add to Complore   Add to Connotea   Add to Del.icio.us   Add to Digg  
What's this?

 Previous Article  | Table of Contents |  Next Article 

Chemistry
Solvent effects on the energy landscapes and folding kinetics of polyalanine

Yaakov Levy, Joshua Jortner*, and Oren M. Beckerdagger

Department of Chemical Physics, School of Chemistry, Tel Aviv University, Ramat Aviv, Tel Aviv 69978, Israel

Contributed by Joshua Jortner, December 22, 2000


    Abstract
Top
Abstract
Introduction
Methods
Conclusion
References

The effect of a solvation on the thermodynamics and kinetics of polyalanine (Ala12) is explored on the basis of its energy landscapes in vacuum and in an aqueous solution. Both energy landscapes are characterized by two basins, one associated with alpha -helical structures and the other with coil and beta -structures of the peptide. In both environments, the basin that corresponds to the alpha -helical structure is considerably narrower than the basin corresponding to the beta -state, reflecting their different contributions to the entropy of the peptide. In vacuum, the alpha -helical state of Ala12 constitutes the native state, in agreement with common helical propensity scales, whereas in the aqueous medium, the alpha -helical state is destabilized, and the beta -state becomes the native state. Thus solvation has a dramatic effect on the energy landscape of this peptide, resulting in an inverted stability of the two states. Different folding and unfolding time scales for Ala12 in hydrophilic and hydrophobic chemical environments are caused by the higher entropy of the native state in water relative to vacuum. The concept of a helical propensity has to be extended to incorporate environmental solvent effects.


    Introduction
Top
Abstract
Introduction
Methods
Conclusion
References

The chemical environment exerts a fundamental influence on the structure, thermodynamics, and dynamics of polypeptides. Particularly, the solvent may affect the dynamics and structure of the polypeptide and consequently alter its function, as has been demonstrated in a variety of biologically important phenomena, ranging from the rate of oxygen uptake in myoglobin to the stabilization of opposite- charged side-chain pairs at the surface of proteins (1, 2). The variance of the properties of a polypeptide in different solvents depends on the nature of the solvent-polypeptide intermolecular interactions, which lead to a rich repertoire of phenomena. These interactions involve the effect of organic solvents on the destabilization of the hydrophobic core and the exposure of side chains, as well as the opposite effects of aqueous solvents on the protein structures favoring the hydrophilic protein surface and the hydrophobic core (1, 3). Theoretical and computational evidence (4-6) for medium effects on polypeptide structures has accumulated. Following the Zimm-Bragg theory of the helix-coil equilibrium (7), it has been established that short polypeptides should not form helices in water. Indeed, numerous studies report that the relative tendency for helix formation in water is low at physiological temperatures (6, 8).

The structure of polyalanine peptides is of considerable interest as, in general, regardless of specific chemical environments, the commonly reported secondary structure propensity scales for amino acids (9-11) rank alanine as having the highest alpha -helical propensity. However, experimental (12-14) and computational (4-6, 15-23) studies showed that polyalanines tend to adopt random-coil conformations in aqueous solution. The ambiguity of the helical propensities, even for alanine, which is known as an excellent helix former, may indicate that these helical propensity scales do not reflect the intrinsic properties of individual residues irrespective of the environment, and that solvent effects have to be taken into account.

The thermodynamics and kinetics of a protein are determined by its energy landscape (24, 25). The relation between energy landscapes and kinetics of complex systems, e.g., clusters (26, 27) and proteins (28), is being explored in terms of a master equation. Previous studies by Y.L. and O.M.B. showed how structural constraints (29, 30) and several specific point mutations (31) affect the topography and topology of the energy landscape of peptides and proteins. Medium effects on the energy landscape of peptides will elucidate the features of solvation on the structure, thermodynamics, and dynamics of these complex systems.

In this paper, we report the observation of dramatically different energy landscapes of a peptide (dodecaalanine, Ala12) in vacuum and in water, which reflects the different environmental chemical properties of its kinetics and thermodynamic stability in hydrophobic and hydrophilic solvents. Polyalanine was chosen to explore the effect of solvent on the energy landscape, because previous experimental (12-14) and computational (4-6, 15-23) studies demonstrate its structural sensitivity to the solvent environment. Two energy landscapes of Ala12 are presented. The first corresponds to the peptide conformational space in vacuum and represents the peptide energy landscape in a hydrophobic medium. The second corresponds to the conformational space in an implicit water solvent and represents the peptide landscape in aqueous solution. A solvent-induced switchover from the native alpha -helical state in vacuum to the native beta -state in aqueous solution is manifested, whereas the implications of the environmentally inverted stability for the thermodynamics and dynamics of this peptide are explored.


    Methods
Top
Abstract
Introduction
Methods
Conclusion
References

Reconstruction of an energy landscape relies on collecting a large sample of conformations, minimizing each of them to the nearest locally stable minima, and then using these local minima to characterize the landscape. To obtain an overview of the molecular energy landscape accessible to the polyalanine at physiological temperatures, we sample the conformation space of Ala12 in each environment by using a high-temperature molecular dynamics trajectory. To ensure that the conformational sample will cover the entire conformation space available to the peptide, we started the high-temperature molecular dynamics trajectory from two distinct Ala12 conformers: an ideal alpha -helix and a beta -hairpin.

Technically, in each environment the sampling procedure includes two 1-ns molecular dynamics trajectories [performed with the CHARMM program (32)], one at 400 K and one at 500 K, for each of the two starting structures. Conformations are sampled every 4 ps along the high-temperature trajectories, resulting in a total of 1,000 conformations of Ala12 in each environment. The protocols of the simulations are similar to those used in previous work (29, 30). The simulations for Ala12 in water used the EEF1 implicit solvation model to take into account the solvent effect (33, 34). The EEF1 expresses the solvation free energy as a sum over group contributions, where the solvation free energy of each group is corrected for screening by the surrounding groups. In addition to the sampling simulations, two 350-K trajectories in vacuum, as well as two 350-K trajectories in water, were performed to study the folding and unfolding kinetics of the system in these environments. The folding trajectory started from a beta -hairpin conformation (0% alpha -helical and 83% beta -sheet), and the unfolding trajectory started from an alpha -helical conformation (100% alpha -helical and 0% beta -sheet). The transition time of the folding reaction, beta right-arrowalpha , was defined as the first passage time (tau ) from the beta -hairpin to an alpha -helical structure with 100% alpha -helical content. Similarly, the transition time of the unfolding reaction, alpha right-arrowbeta , was defined as the time required to unwind the helix and to obtain a structure with 0% alpha -helical content and 83% beta -sheet. The alpha -helical and beta -sheet contents of any conformation of the conformation samples and of the folding/unfolding trajectories were calculated by using the DSSP program (35). In addition, to estimate the effect of solvation on the stability of the two states of Ala12 in vacuum, we collected 200 alpha -helical conformations and 200 beta -hairpin conformations of Ala12, each sampled along 0.4-ns trajectories at 300 K by using vacuum, the EEF1 implicit solvation model, and the TIP3P explicit water models (36).

The construction of the energy landscapes and the projections of the folding and unfolding trajectories onto a low-dimensional subspace (37-39) were obtained by applying the principal component analysis method (PCA) (38, 39). In this study, we use the root-mean-square distance of the backbone heavy atoms as the distance measure between conformations composed of the conformation samples in vacuum and in solvent. In this case, the principal two-dimensional subspace was found to represent the multidimensional data to an accuracy greater than 50%.

Solvent Effect on the Flexibility of Ala12. The flexibility of polypeptides in different solvents is expected to differ because of the different type of intermolecular interactions between the solvent and peptide molecules. For example, whereas the alpha -helical state of polyalanine is stable in organic media (4, 19, 23), its stability decreases in aqueous solutions, because water molecules interact with the peptide polar groups, which form the CO(i)right-arrowNH(i + 4) hydrogen bonds. Inevitably, the destabilization of the alpha -helix hydrogen bonds may result in destabilization and even unwinding of the entire polyalanine alpha -helix. Accordingly, the helical conformation of polyalanine may be more dominant in an organic medium than in water, whereas the coil structures, in which the polar groups are exposed to the solvent, are likely to be more accessible in water than in organic medium. The larger the number of conformations that a peptide can adopt, the larger are its entropy and flexibility.

The flexibility of Ala12 in each environment can be represented by the volume of conformation space that it occupies in that specific environment. To span the entire conformation space accessible to Ala12 in each environment, we sampled conformations by simulating the peptide at high temperatures (400 K and 500 K) starting from both alpha -helical and beta -sheet conformations. This sampling procedure was selected to overcome energetic and entropic barriers. The sampled conformations were annealed back to 300 K to ensure that the resulting conformations are accessible at physiological temperatures. The multidimensional conformation spaces of Ala12 in vacuum and in aqueous solvent were jointly projected onto the same three-dimensional subspace as shown in Fig. 1a. A comparison between the relative volume of each conformation sample and their overlap clearly outlines the effect of solvent on the flexibility of the peptide. For Ala12, a larger volume of conformation space is accessible in the aqueous solvent than in vacuum, indicating that in water the peptide can adopt conformations that are unlikely in vacuum. Fig. 2a shows three of the conformations sampled in an aqueous environment selected from the region that does not overlap with the conformation space of Ala12 in vacuum. These structures illustrate that the larger flexibility of Ala12 in aqueous solution is mainly because of unwinding of the helix and a solvation of all backbone polar groups that were involved in the stabilization of the alpha -helix. In water, random-coil structures are reasonable because the loss of intramolecular interaction is compensated by intermolecular interactions with water molecules.



View larger version (52K):
[in this window]
[in a new window]
 
Fig. 1.   A joint two-dimensional projection of the conformation space of Ala12. (a) Conformations sampled in vacuum (full circles) and with implicit water molecules (empty triangles). (b) Conformations in both samples with alpha -helical content larger than 50%. (c) Conformations in both samples with beta -sheet content larger than 50%. Ellipses are drawn to emphasize the Ala12 conformations in vacuum (solid line) and in solvent (dashed line). Conformations 1-9 are shown in Fig. 2.



View larger version (29K):
[in this window]
[in a new window]
 
Fig. 2.   Conformations of Ala12. (a) Conformations selected from the region in the conformation space, which is accessible only in an aqueous solvent (Fig. 1a). (b) Conformations selected from the region of alpha -helical conformations (Fig. 1b). (c) Conformations selected from the region of beta -sheet conformations (Fig. 1c).

A partial projection of only those conformations (in both conformation samples) that have an alpha -helical content larger than 50% is shown in Fig. 1b. A partial projection of conformations (in both conformation samples) that are characterized by a beta -sheet content larger than 50% is shown in Fig. 1c. The relatively smaller region of the helical conformations in both conformation spaces (Fig. 1b) in comparison to the region of beta -sheet conformations (Fig. 1c) indicates that helical structures are much more rigid than beta -hairpin structures. Moreover, whereas the regions corresponding to alpha -helical conformations sampled in vacuum and in aqueous solvent significantly overlap, as expected, less overlap is observed between the regions corresponding to beta -sheet structures. Three alpha -helical conformations of Ala12 placed in the alpha -helical region (Fig. 1b) are shown in Fig. 2b. Similarly, Fig. 2c shows three beta -structures selected from the corresponding region in conformation spaces (Fig. 1c). These structures illustrate the larger flexibility of the beta -sheet in comparison with the alpha -helical structures.

The accessibility and equilibrium population of the alpha -helical and the beta -sheet structures of both conformational samples, in vacuum and in aqueous solution, is expected to be different because their relative stability and the barrier heights that separate them are affected by the environment. Only the energy landscape can explain the relative stability of the conformations of the sampled conformation space and the dynamics between them. For Ala12 in vacuum as well as in water, the energy landscape can be constructed, being based on the conformational samples, by plotting the energy of each conformation versus the geometrical axes obtained by applying the PCA procedure (38).

Energy Landscape of Ala12 in Vacuum. The energy landscape of Ala12 in vacuum (Fig. 3a) is composed of two main basins. One basin is shallow and broad, and the other is narrow and deep. In fact, these two basins, which are different in energy, correspond to conformations with different secondary structures. The broad and shallow basin mainly corresponds to coil and beta -structures, whereas the narrow and deep basin corresponds to alpha -helical conformations (Fig. 3b). A comparison between the sizes of the two basins reflects the relative flexibility of the two corresponding structures. The broad basin of the coil and the beta -structures, in comparison to the narrow basin of the alpha -helical structures, illustrates that the coil and beta -structures are much more flexible and occupy most of the peptide conformation space. Namely, the basin of the beta -structures includes various structures such as beta -hairpin, beta -strand, and random coil, whereas the basin corresponding to the alpha -helical structures includes structures with different alpha -helical content (examples for alpha -helical and beta -structures are shown in Fig. 2 b and c, respectively).



View larger version (26K):
[in this window]
[in a new window]
 
Fig. 3.   (a) The two-state energy landscape of Ala12 in vacuum. (b) The secondary structure characteristic of the two states that compose the Ala12 energy landscape.

The Ala12 energy landscape in vacuum indicates that the peptide has two structural states: an alpha -helical and a beta -state (which includes coil and beta -structures). The alpha -helical conformation in vacuum is more stable by 11.6 kcal/mol than the beta -hairpin conformation (Table 1). The greater stability of the alpha -helical state in vacuum is in accord with common helical propensity scales, which rank alanine as the residue with the highest helical propensity. However, whereas Ala12 adopts an alpha -helix at physiological temperatures, the energy landscape suggests that the energetic barrier crossing is possible at higher temperatures, and a transition from the alpha -helical to the beta -basin can occur. In the alpha right-arrowbeta (beta right-arrowalpha ) processes in vacuum, we expect that Delta E > 0 (<0) and Delta S > 0 (<0). When a barrier crossing is possible, the reaction alpha  left-right-arrow beta  is in favor of the beta -state because of its large configurational entropy that decreases the free energy change of the reaction as the temperature increases.


                              
View this table:
[in this window]
[in a new window]
 
Table 1.   Properties of Ala12 in vacuum and in aqueous solvent based on its energy landscapes in both environments

Fig. 4a shows a 350-K folding trajectory, from a beta -hairpin conformation to an alpha -helical conformation, projected on the two-state vacuum energy landscape of Ala12. The projection of the trajectory onto the energy landscape demonstrates the two competing factors of a protein-folding process: energy versus entropy. In the case of Ala12, the peptide spends 10.7 ns in the highly entropic beta -basin (nonnative state) before it undergoes a transition to the energetically favorable alpha -helical state (native state) and adopts an ideal alpha -helix conformation after 11.2 ns (Table 1). In the opposite case of the 350-K unfolding simulation shown in Fig. 4b, starting with the ideal alpha -helix of Ala12, it took 8.8 ns to escape from the alpha -helix basin into the nonnative beta -basin adopting a beta -hairpin structure after 9.2 ns. The heavy lines in the projected folding/unfolding trajectories of Fig. 4 a and b emphasize the successive structures before entering/escaping the alpha -helical basin. These two trajectories show that there are several, perhaps many, pathways for these reactions, and that this system has a transition region with more than one transition state, rather than a single transition state, as classical chemical kinetics would suggest (35, 40). The two trajectories are presented here as an indication of the effect of entropy on the folding and unfolding processes. However, to infer from the transition times of the two processes, more trajectories should be collected, starting with different initial conditions.



View larger version (55K):
[in this window]
[in a new window]
 
Fig. 4.   Folding (a) and unfolding (b) trajectories of Ala12 in vacuum at 350 K projected onto the underlying energy landscape. The heavy lines indicate the 20 successive conformations before entering/escaping the native basin in the folding/unfolding trajectories, respectively.

The folding and unfolding pathways of Ala12 once projected on the two-state energy landscape illustrate the significance of the energy and entropy components of the "folded"/"unfolded" states to the folding/unfolding kinetics. The inherent flexibility of proteins results in broad basins, which dominate the folding and unfolding thermodynamics and kinetics. The large difference in the entropy of the two states may suggest that with increasing temperature, the free energy change for the folding process will be more negative and consequently may result in non-Arrhenius kinetic behavior (41).

Energy Landscape of Ala12 in an Aqueous Solvent. The gross features of the energy landscape of Ala12 in aqueous solution (Fig. 5a) are qualitatively similar to the corresponding energy landscape in vacuum, which is composed of two basins. One basin is associated with alpha -helical structures, and the other corresponds to coil and beta -structures. However, the quantitative features of the two basins are drastically different in aqueous solution and in vacuum. Although in vacuum the alpha -helical state constitutes the stable state (Fig. 3a), this state becomes destabilized in an aqueous environment, and the stable state is the one related to the beta -structure. Fig. 5a depicts the two-state energy landscape of Ala12 in implicit water projected onto a two-dimensional subspace. The secondary structure characteristic of the two basins is shown in Fig. 5b. In this case, unlike in vacuum, the two first principal coordinates obtained by the PCA procedure are not sufficient to distinguish between the two states, and additional coordinates are required to better characterize the energy landscape. To overcome this difficulty, we have separately plotted (Fig. 5a) the landscape for all structures with an alpha -helical content of 50% or more (the higher and narrower surface) and the landscape corresponding to all other conformations of this polypeptide. The lesser success of the PCA method in capturing the topography of the energy landscape in water originates from the fact that in water, Ala12 is more flexible than in vacuum and can populate a broader range of conformations with the root-mean-square distance being comparable to that between alpha -helical and beta -sheet conformations.



View larger version (30K):
[in this window]
[in a new window]
 
Fig. 5.   (a) The two-state energy landscape of Ala12 in implicit water. The broad basin (white surface) corresponds to coil and beta -structures, and the narrower and higher basin (pink surface) corresponds to alpha -helical structures. (b) The secondary structure characteristic of the two states composed the Ala12 energy landscape.

As is seen in Fig. 5a, the alpha -helical state in an aqueous medium is less stable by 9.4 kcal/mol relative to the beta -hairpin conformation, which is the global minimum of Ala12 in water (Table 1). For the beta right-arrowalpha (alpha right-arrowbeta ) reaction in solution, we have Delta E > 0 (<0) and expect Delta S < 0 (>0), because of the large configurational entropy of the beta -structures. The lower stability of the helical state of the peptide in water is because of solvation of the polar group, which otherwise participates in the CO(i)right-arrowNH(i + 4) hydrogen-bond network that defines the alpha -helix structure. The energy gap between the beta -hairpin and the other beta -structures is about 4 kcal/mol. In water there is only a relatively small thermodynamic preference for the beta -hairpin structure in comparison to the higher preference of the alpha -helix in vacuum. This hints to the fact that in aqueous solution, polyalanine may adopt a variety of open structures. The high flexibility of Ala12 in water is reflected by the broadness of the native basin that corresponds to the coil and beta -structures. To estimate the effect of solvation on the stability of the two states of Ala12 in vacuum, we calculate the average energies of Ala12's alpha -helical and beta -hairpin conformations collected from 0.4-ns trajectories at 300 K by using vacuum, implicit solvent, and explicit solvent models (see Table 2). The results of both implicit and explicit solvent calculations indicate that solvation destabilizes the helical structure of Ala12 and that in a hydrophilic environment, the beta -hairpin structure is the more stable state of Ala12. However, whereas the energy gap between the two states is about 10 kcal/mol in vacuum, the gap in aqueous solvent is only 4-5 kcal/mol.


                              
View this table:
[in this window]
[in a new window]
 
Table 2.   The stability of alpha -helical and beta -hairpin conformations of Ala12 in vacuum, implicit and explicit solvent, kcal/mol*

Clearly the difference between the energy landscapes of Ala12 in water and in vacuum indicates that the peptide exhibits different thermodynamic and kinetic behaviors in the two environments. Fig. 6a shows a 350-K trajectory of helixright-arrowcoil transition in implicit water projected onto the energy landscape of the solvated Ala12. The time required for the helix to unwind is 3.9 ns, being faster by more than a factor of two than for the same process in vacuum. In explicit water at 300 K, Ala12 unfolds after 1 ns, indicating that although the unfolding process in implicit solvent is faster than in vacuum, it does not entirely represent the solvent effect on the folding/unfolding reactions. Recall that in vacuum, the helix-to-coil process is driven only by increasing entropy, whereas in an aqueous medium, the process is also driven by decreasing energy, which, together with the increase in entropy, makes the helixright-arrowcoil transition much faster. On the other hand, the reverse process, coilright-arrowhelix, is predicted to be much slower in an aqueous environment relative to its rate in vacuum. Although in vacuum this process proceeds because of decreasing energy, in solvent the process is characterized by positive free energy change because of an increase in energy and a decrease in entropy (see Table 1). Fig. 6b shows a projection of a 100-ns 350-K trajectory on the energy landscape of Ala12 in aqueous solution. Because during these 100 ns a transition to the nonnative state, which corresponds to the alpha -helical structures, was not observed, the unfolding trajectory is presented only with respect to the native basin. The coilright-arrowhelix transition time in aqueous solution is longer than 100 ns, being slower by more than an order of magnitude than the corresponding transition time in vacuum (approx 11 ns, according to Fig. 4a). The different kinetics behavior can be traced to the environmental effect on the energy landscapes (Figs. 3a and 5a) that makes the free energy change of the coilright-arrowhelix transition positive over the entire temperature range.



View larger version (51K):
[in this window]
[in a new window]
 
Fig. 6.   Folding (a) and unfolding (b) trajectories of Ala12 in an aqueous solvent at 350 K projected onto the underlying energy landscape. The heavy line indicates the 20 successive conformations before entering the native basin in the folding trajectory. Because a transition to the unfolded state was not observed during the 100-ns unfolding trajectory, only the native basin is shown (Fig. 5b).


    Conclusion
Top
Abstract
Introduction
Methods
Conclusion
References

Our study demonstrates that the environment has a significant effect on the energy landscape of Ala12 polypeptide. Although the energy landscapes of Ala12 in hydrophilic and hydrophobic media are both composed of two basins, one associated with alpha -helical and the other with coil and beta -structures, their relative properties are markedly affected by the environment. In both media, the large flexibility of the beta -structures implies that most of the peptide configurational entropy originates from this state, whereas the contribution of the alpha -helical state to the total entropy is smaller. The main effect of the environment is exerted on the relative stability of the two states that make up the energy landscape. Although the alpha -helical state is the stable state in a hydrophobic environment, it is destabilized under hydrophilic conditions, and then the beta -state becomes the more stable one. The inversion of the stability of the two states because of polar solvation results in a native state in aqueous solution that is characterized by larger entropy and, consequently, by a very efficient folding process.

Accordingly, polyalanine is expected to adopt an alpha -helical conformation in a nonpolar organic solvent and beta -structures with coil conformations in a polar aqueous solution. Solid-state (i.e., hydrophobic environment) (42-45) and aqueous solution (12-14) measurements of polyalanine support its conformational dependence on the chemical environments with the ubiquity of the alpha -helix and beta -structures prevailing in hydrophobic and in polar environments, respectively. The dependence of the conformation of alanine on the solvent characteristics (i.e., the polarity and dielectric constant of the solvent) is supported by experimental evidence that the helical propensity of alanine in water shows a dramatic increase on addition of certain alcohols (e.g., trifluoroethanol) (46, 47). Recently, it was proposed that alcohol may act on the exposed CO and NH groups by diminishing their exposure to the solvent, i.e., shifting the conformational equilibrium toward more compact structures, such as alpha -helical conformation (23). An alternative proposal for shifting the helix-coil equilibrium of polyalanine toward an alpha -helical conformation is to introduce charged residues into the sequence (48, 49); however, in this case, the tendency to form the alpha -helix should not be associated with the alpha -helical propensity of alanine discussed herein.

Additional information emerges concerning temperature-dependent conformations. The energy landscape of Ala12 in a nonpolar medium implies that the alpha right-arrowbeta reaction is thermodynamically characterized by Delta E > 0 and Delta S > 0, whereon Delta G for this reaction will decrease with increasing temperature. Consequently, the equilibrium of the alpha right-arrowbeta reaction in a hydrophobic environment is predicted to shift toward the beta -state with increasing temperature. Thus, the significant preference of the alpha -helical conformation of polyalanine in a hydrophobic environment at 300 K can be reduced at higher temperatures. Experimentally, solid-state measurements of polyalanine indicate that the conformation of Ala200 is shifted from alpha -helix at lower temperatures to beta -sheet at higher temperatures (44). The energy landscape of Ala12 in an aqueous solvent suggests that because the beta right-arrowalpha reaction is characterized by Delta E > 0 and Delta S < 0 (i.e., Delta G > 0 for all of the temperature range), the temperature does not affect the ratio of the alpha - and beta -conformations of polyalanine in water, and thus Ala12 adopts beta -structures in aqueous solutions.


    Acknowledgements

We are grateful to Professor R. S. Berry for stimulating comments.


    Abbreviation

PCA, principal component analysis method.


    Footnotes

* To whom reprint requests should be addressed. E-mail: jortner{at}chemsg1.tau.ac.il.

dagger Present address: Bio Information Technologies (Bio-I.T.) Ltd., 1 Betzalel St., Ramat Gan 52521, Israel.


    References
Top
Abstract
Introduction
Methods
Conclusion
References

1. Creighton, T. E. (1993) Proteins: Structures and Molecular Properties (Freeman, New York).
2. Brooks, C. L. I. , Karplus, M. & Pettit, B. M. (1988) Proteins: A Theoretical Perspective of Dynamics, Structure, and Thermodynamics (Wiley, New York).
3. Fersht, A. (1999) Structure and Mechanism in Protein Science: A Guide to Enzyme Catalysis and Protein Folding (Freeman, New York).
4. Doruker, P. & Bahar, I. (1997) Biophys. J. 72, 2445-2456.
5. Daggett, V. & Levitt, M. (1992) J. Mol. Biol. 223, 1121-1138.
6. Soman, K. V. , Karimi, A. & Case, D. A. (1991) Biopolymers 31, 1351-1361.
7. Zimm, B. H. & Bragg, J. K. (1959) J. Chem. Phys. 31, 526-535.
8. Buuren, A. R. V. & Berendsen, H. J. C. (1993) Biopolymers 33, 1159-1166.
9. Chou, P. Y. & Fasman, G. D. (1978) Annu. Rev. Biochem. 47, 251-276.
10. O'Neil, K. T. & DeGrado, W. F. (1990) Science 250, 646-651.
11. Chakrabartty, A. , Kortemme, T. & Baldwin, R. L. (1994) Protein Sci. 3, 843-852.
12. Ingwall, R. T. , Scheraga, H. A. , Lotan, N. , Berger, A. & Katchalski, E. (1968) Biopolymers 6, 331-368.
13. Platzer, K. E. B. , Ananthanarayanan, V. S. , Andreatta, R. H. & Scheraga, H. A. (1972) Macromolecules 5, 177-187.
14. Blondelle, S. E. , Forood, B. , Houghten, R. A. & Perez-Paya, E. (1997) Biochemistry 36, 8393-8400.
15. Daggett, V. , Kollman, P. A. & Kuntz, I. D. (1991) Biopolymers 31, 1115-1134.
16. Sung, S. S. (1994) Biophys. J. 66, 1796-1803.
17. Yang, A. S. & Honig, B. (1995) J. Mol. Biol. 252, 351-365.
18. Yang, A. S. & Honig, B. (1995) J. Mol. Biol. 252, 366-376.
19. Shen, L. , Bassolino, D. & Stouch, T. (1997) Biophys. J. 73, 3-20.
20. Bertsch, R. A. , Vaidehi, N. , Chan, S. I. & Goddard, W. A. I. (1998) Proteins Struct. Funct. Genet. 33, 343-357.
21. Huo, S. & Straub, J. E. (1999) Proteins Struct. Funct. Genet. 36, 249-261.
22. Takano, M. , Yamato, T. , Higo, J. , Suyama, A. & Nagayama, K. (1999) J. Am. Chem. Soc. 121, 605-612.
23. Vila, J. A. , Ripoll, D. R. & Scheraga, H. A. (2000) Proc. Natl. Acad. Sci. USA 97, 13075-13079. (First Published November 14, 2000; 10.1073/pnas.240455797)
24. Dobson, C. M. , Sali, A. & Karplus, M. (1998) Angew. Chem. 37, 868-893.
25. Frauenfelder, H. & McMahon, B. (2000) Ann. Phys. (Leipzig) 9, 655-667.
26. Ball, K. D. & Berry, R. S. (1998) J. Chem. Phys. 109, 8557-8572.
27. Miller, M. A. , Doye, J. P. K. & Wales, D. J. (1999) Phys. Rev. E 60, 3701-3718.
28. Becker, O. M. & Karplus, M. (1997) J. Chem. Phys. 106, 1495-1517.
29. Levy, Y. & Becker, O. M. (1998) Phys. Rev. Lett. 81, 1126-1129.
30. Levy, Y. & Becker, O. M. (2001) J. Chem. Phys. 114, 993-1010.
31. Levy, Y. & Becker, O. M. (2001) in Conformational Diseases, eds. Katzir, E., Solomon, B. & Taraboulos, A. (Center for the Study of Emerging Diseases, Jerusalem).
32. Brooks, B. R. , Bruccoleri, R. E. , Olafson, B. D. , States, D. J. , Swaminathan, S. & Karplus, M. (1983) J. Comput. Chem. 4, 187-217.
33. Lazaridis, T. & Karplus, M. (1997) Science 278, 1928-1931.
34. Lazaridis, T. & Karplus, M. (1999) Proteins Struct. Funct. Genet. 35, 133-152.
35. Martinez, J. C. , Pisabarro, M. T. & Serrano, L. (1998) Nat. Struct. Biol. 5, 721-729.
36. Jorgensen, W. L. , Chandrasekhar, J. , Madura, J. D. , Impey, R. W. & Klien, M. L. (1983) J. Chem. Phys. 79, 926-935.
37. Elmaci, N. & Berry, R. S. (1999) J. Chem. Phys. 110, 10606-10622.
38. Becker, O. M. (1998) J. Comput. Chem. 19, 1255-1267.
39. Becker, O. M. , Levy, Y. & Ravitz, O. (2000) J. Phys. Chem. 104, 2123-2135.
40. Onuchic, J. N. , Socci, N. D. , Luthey-Schulten, Z. & Wolynes, P. G. (1996) Folding Des. 1, 441-450.
41. Karplus, M. (2000) J. Phys. Chem. B 104, 11-27.
42. Shoji, A. , Ozaki, T. , Fujito, T. , Deguchi, K. , Ando, S. & Ando, I. (1990) J. Am. Chem. Soc. 112, 4693-4697.
43. Kimura, H. , Ozaki, T. , Sugisawa, H. , Deguchi, K. & Shoji, A. (1998) Macromolecules 31, 7398-7403.
44. Lee, D.-K. & Ramamoorthy, A. (1999) J. Phys. Chem. B 103, 271-275.
45. Warras, R. , Wieruszeski, J.-M. , Boutillon, C. & Lippens, G. (2000) J. Am. Chem. Soc. 122, 1789-1795.
46. Rohl, C. A. , Chakrabartty, A. & Baldwin, R. L. (1996) Protein Sci. 5, 2623-2637.
47. Cammers-Goodwin, A. , Allen, T. J. , Oslick, S. L. , McClure, K. F. , Lee, J. K. & Kemp, D. S. (1996) J. Am. Chem. Soc. 118, 3082-3090.
48. Marqusee, S. , Robbins, V. H. & Baldwin, R. L. (1989) Proc. Natl. Acad. Sci. USA 86, 5286-5290.
49. Williams, L. , Kather, K. & Kemp, D. S. (1998) J. Am. Chem. Soc. 120, 11033-11043.
www.pnas.org/cgi/doi/10.1073/pnas.041611998
Add to CiteULike CiteULike   Add to Complore Complore   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us   Add to Digg Digg    What's this?


This article has been cited by other articles in HighWire Press-hosted journals:


Home page
Proc. Natl. Acad. Sci. USAHome page
L. Dougan, G. Feng, H. Lu, and J. M. Fernandez
Solvent molecules bridge the mechanical unfolding transition state of a protein
PNAS, March 4, 2008; 105(9): 3185 - 3190.
[Abstract] [Full Text] [PDF]


Home page
Biophys. JHome page
S. Tanizaki, J. Clifford, B. D. Connelly, and M. Feig
Conformational Sampling of Peptides in Cellular Environments
Biophys. J., February 1, 2008; 94(3): 747 - 759.
[Abstract] [Full Text] [PDF]


Home page
Proc. Natl. Acad. Sci. USAHome page
U. Heugen, G. Schwaab, E. Brundermann, M. Heyden, X. Yu, D. M. Leitner, and M. Havenith
Solute-induced retardation of water dynamics probed directly by terahertz spectroscopy
PNAS, August 15, 2006; 103(33): 12301 - 12306.
[Abstract] [Full Text] [PDF]


Home page
Biophys. JHome page
J. Gumbart and K. Schulten
Molecular Dynamics Studies of the Archaeal Translocon
Biophys. J., April 1, 2006; 90(7): 2356 - 2367.
[Abstract] [Full Text] [PDF]


Home page
Biophys. JHome page
V. M. Dadarlat
Potentials of Mean Force for the Interaction of Blocked Alanine Dipeptide Molecules in Water and Gas Phase from MD Simulations
Biophys. J., September 1, 2005; 89(3): 1433 - 1445.
[Abstract] [Full Text] [PDF]


Home page
Biophys. JHome page
D. Shental-Bechor, S. Kirca, N. Ben-Tal, and T. Haliloglu
Monte Carlo Studies of Folding, Dynamics, and Stability in {alpha}-Helices
Biophys. J., April 1, 2005; 88(4): 2391 - 2402.
[Abstract] [Full Text] [PDF]


Home page
Protein Sci.Home page
H. D. Nguyen, A. J. Marchut, and C. K. Hall
Solvent effects on the conformational transition of a model polyalanine peptide
Protein Sci., November 1, 2004; 13(11): 2909 - 2924.
[Abstract] [Full Text] [PDF]


Home page
Proc. Natl. Acad. Sci. USAHome page
S. V. Krivov and M. Karplus
Hidden complexity of free energy surfaces for peptide (protein) folding
PNAS, October 12, 2004; 101(41): 14766 - 14770.
[Abstract] [Full Text] [PDF]


Home page
Protein Sci.Home page
I. Chang, M. Cieplak, J. R. Banavar, and A. Maritan
What can one learn from experiments about the elusive transition state?
Protein Sci., September 1, 2004; 13(9): 2446 - 2457.
[Abstract] [Full Text] [PDF]


Home page
Hum Mol GenetHome page
H. Lavoie, F. Debeane, Q.-D. Trinh, J.-F. Turcotte, L.-P. Corbeil-Girard, M.-J. Dicaire, A. Saint-Denis, M. Page, G. A. Rouleau, and B. Brais
Polymorphism, shared functions and convergent evolution of genes with sequences coding for polyalanine domains
Hum. Mol. Genet., November 15, 2003; 12(22): 2967 - 2979.
[Abstract] [Full Text] [PDF]


Home page
Proc. Natl. Acad. Sci. USAHome page
C.-Y. Huang, Z. Getahun, Y. Zhu, J. W. Klemke, W. F. DeGrado, and F. Gai
Helix formation via conformation diffusion search
PNAS, February 20, 2002; (2002) 52700099.
[Abstract] [Full Text] [PDF]


Home page
Proc. Natl. Acad. Sci. USAHome page
C.-Y. Huang, Z. Getahun, Y. Zhu, J. W. Klemke, W. F. DeGrado, and F. Gai
Helix formation via conformation diffusion search
PNAS, March 5, 2002; 99(5): 2788 - 2793.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a colleague
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science