[Frontiers in Bioscience 14, 1292-1303, January 1, 2009]

Free-energy landscapes of proteins in solution by generalized-ensemble simulations

Yuji Sugita1,2

1Advanced Science Institute, RIKEN, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan, 2CREST & BIRD, Japan Science and Technology Agency, 4-1-8 Honcho Kawaguchi, Saitama 332-0012, Japan


1. Abstract
2. Introduction
3. Generalized-ensemble algorithms
3.1. Multicanonical algorithm (MUCA)
3.2. Replica-exchange method (REM)
3.3. Further extensions of MUCA and REM
4. Folding simulations of the C-peptide in explicit solvent
5. Structural changes in the cytoplasmic domain of phospholamban by phosphorylation at Ser16
6. Perspective
7. Acknowledgements
8. References


Free-energy landscapes of proteins in solution are essential for understanding molecular mechanism of protein folding, stability, and dynamics. Because of the multiple-minima problem (or quasi-ergodicity problem), the conventional molecular dynamics or Monte Carlo methods cannot provide the landscapes accurately at low temperatures. By contrast, the simulations based on the generalized-ensemble algorithms can sample wider conformational spaces than the conventional approaches, thereby providing better free-energy landscapes of proteins at low temperatures. In this article, we review two well-known generalized-ensemble algorithms, namely, multicanonical algorithm and replica-exchange method, and then introduce further extensions of the above two methods, which are applicable to larger systems with rugged energy landscapes. These simulation methods have been applied to the protein folding simulations of the C-peptide in ribonuclease A with explicit solvent. We also demonstrate how the methods and the free-energy landscapes of proteins are useful for the biological research, by showing the simulation results on the phospholamban, a reversible regulator of sarco(end)plasmic reticulum Ca2+-pump.