期刊名称：Proceedings of the National Academy of Sciences
印刷版ISSN：0027-8424
电子版ISSN：1091-6490
出版年度：2003
卷号：100
期号：24
页码：13934-13939
DOI：10.1073/pnas.2232868100
语种：English
出版社：The National Academy of Sciences of the United States of America
摘要：We compare simulations using the generalized Born/surface area (GB/SA) implicit solvent model with simulations using explicit solvent (transferable intermolecular potential 3 point, TIP3P) to test the GB/SA algorithm. We use the replica exchange molecular dynamics method to sample the conformational phase space of two -helical peptides, A21 and the Fs, by using two different classical potentials and both water models. We find that when using GB/SA: (i) A21 is predicted to be more helical than the Fs peptide at all temperatures; (ii) the native structure of the Fs peptide is predicted to be a helical bundle instead of a single helix; and (iii) the persistence length and most probable end-to-end distance are too large in the unfolded state when compared against the explicit solvent simulations. We find that the potential of mean force in the {phi}{psi} plane is markedly different in the two solvents, making the two simulated peptides respond differently when the backbone torsions are perturbed. A fit of the temperature melting curves obtained in these simulations to a Lifson-Roig model finds that the GB/SA model has an unphysically large nucleation parameter, whereas the explicit solvent model produces values similar to experiment.
