Replica Exchange and Multicanonical Algorithms with the coarse-grained UNRES force field.
Three algorithms, namely a Replica Exchange method (REM), a Replica Exchange Multicanonical method (REMUCA), and Replica Exchange Multicanonical with Replica Exchange (REMUCAREM), were implemented with the coarse-grained united-residue force field (UNRES) in both Monte Carlo and Molecular Dynamics versions. The MD algorithms use the constant-temperature Berendsen thermostat, with the velocity Verlet algorithm and variable time step. The algorithms were applied to one peptide (20 residues of Alanine with free ends; ala(20)) and two small proteins, namely an α-helical protein of 46 residues (the B-domain of the staphylococal protein A; 1BDD), and an α+β-protein of 48 residues (the E. Coli Mltd Lysm Domain; 1E0G). Calculated thermodynamic averages, such as canonical average energy and heat capacity, are in good agreement among all simulations for poly-L-alanine, showing that the algorithms were implemented correctly, and that all three algorithms are equally effective for small systems. For protein A, all algorithms performed reasonably well, although some variability in the calculated results was observed whereas, for a more complicated α+β-protein (1E0G), only Replica Exchange was capable of producing reliable statistics for calculating thermodynamic quantities. Finally, from the Replica Exchange molecular dynamics results, we calculated free energy maps as functions of RMSD and radius of gyration for different temperatures. The free energy calculations show correct folding behavior for poly-L-alanine and protein A while, for 1E0G, the native structure had the lowest free energy only at very low temperatures. Hence, the entropy contribution for 1E0G is larger than that for protein A at the same temperature. A larger contribution from entropy means that there are more accessible conformations at a given temperature, making it more difficult to obtain an efficient coverage of conformational space to obtain reliable thermodynamic properties. At the same temperature, ala(20) has the smallest entropy contribution, followed by protein A, and then by 1E0G.
Nanias
Marian
M
Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, NY 14853-1301, U.S.A.
Czaplewski
Cezary
C
Scheraga
Harold A
HA
J Phys Chem B. 2005 Jul 21;109(28):13798-810
16852728
Proc Natl Acad Sci U S A. 2001 Dec 18;98(26):14931-6
11752441
Proteins. 2003 Aug 15;52(3):436-45
12866054
Phys Rev Lett. 1993 Jul 12;71(2):211-214
10054892
J Mol Biol. 1997 Mar 14;266(5):859-65
9086265
J Protein Chem. 1989 Apr;8(2):263-87
2736043
Proc Natl Acad Sci U S A. 1987 Oct;84(19):6611-5
3477791
Phys Rev Lett. 1989 Sep 18;63(12):1195-1198
10040500
Phys Rev Lett. 1992 Jan 6;68(1):9-12
10045099
Proc Natl Acad Sci U S A. 2005 Sep 27;102(39):13749-54
16172406
Phys Rev Lett. 1986 Nov 24;57(21):2607-2609
10033814
Proteins. 2002 Feb 1;46(2):225-34
11807951
J Mol Biol. 2000 Jun 16;299(4):1113-9
10843862
Biochemistry. 1992 Oct 13;31(40):9665-72
1390743
Protein Sci. 1993 Oct;2(10):1715-31
8251944
J Comput Chem. 2005 Nov 15;26(14):1472-86
16088925
Phys Rev Lett. 2003 Aug 1;91(5):058305
12906640
J Am Stat Assoc. 1949 Sep;44(247):335-41
18139350
Proc Natl Acad Sci U S A. 2005 Feb 15;102(7):2362-7
15677316
Phys Rev Lett. 2005 Sep 23;95(13):138102
16197184
Proc Natl Acad Sci U S A. 2003 Nov 25;100(24):13898-903
14623983
Nat Struct Biol. 2000 May;7(5):375-9
10802733
Protein Sci. 1993 Oct;2(10):1697-714
7504550
Proc Natl Acad Sci U S A. 2002 Feb 19;99(4):1937-42
11854494
J Mol Biol. 2006 Jan 20;355(3):536-47
16324712
J Phys Chem B. 2005 Jul 21;109(28):13785-97
16852727
Biopolymers. 2003 Jul;69(3):399-405
12833266
Macromolecules. 1974 Nov-Dec;7(6):797-806
4437206
Biopolymers. 1988 Aug;27(8):1283-303
3219397
Biophys J. 2003 Feb;84(2 Pt 1):775-86
12547762
J Chem Phys. 2004 Aug 8;121(6):2491-504
15281846
J Mol Biol. 1994 Nov 11;243(5):906-18
7966308
Proc Natl Acad Sci U S A. 2005 May 24;102(21):7547-52
15894609