Premium
Quantum mechanics simulation of protein dynamics on long timescale
Author(s) -
Liu Haiyan,
Elstner Marcus,
Kaxiras Efthimios,
Frauenheim Thomas,
Hermans Jan,
Yang Weitao
Publication year - 2001
Publication title -
proteins: structure, function, and bioinformatics
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.699
H-Index - 191
eISSN - 1097-0134
pISSN - 0887-3585
DOI - 10.1002/prot.1114
Subject(s) - molecular dynamics , molecular mechanics , force field (fiction) , quantum , statistical physics , physics , quantum dynamics , quantum mechanics , chemistry , chemical physics , classical mechanics
Protein structure and dynamics are the keys to a wide range of problems in biology. In principle, both can be fully understood by using quantum mechanics as the ultimate tool to unveil the molecular interactions involved. Indeed, quantum mechanics of atoms and molecules have come to play a central role in chemistry and physics. In practice, however, direct application of quantum mechanics to protein systems has been prohibited by the large molecular size of proteins. As a consequence, there is no general quantum mechanical treatment that not only exceeds the accuracy of state‐of‐the‐art empirical models for proteins but also maintains the efficiency needed for extensive sampling in the conformational space, a requirement mandated by the complexity of protein systems. Here we show that, given recent developments in methods, a general quantum mechanical‐based treatment can be constructed. We report a molecular dynamics simulation of a protein, crambin, in solution for 350 ps in which we combine a semiempirical quantum‐mechanical description of the entire protein with a description of the surrounding solvent, and solvent‐protein interactions based on a molecular mechanics force field. Comparison with a recent very high‐resolution crystal structure of crambin (Jelsch et al., Proc Natl Acad Sci USA 2000;102:2246–2251) shows that geometrical detail is better reproduced in this simulation than when several alternate molecular mechanics force fields are used to describe the entire system of protein and solvent, even though the structure is no less flexible. Individual atomic charges deviate in both directions from “canonical” values, and some charge transfer is found between the N and C‐termini. The capability of simulating protein dynamics on and beyond the few hundred ps timescale with a demonstrably accurate quantum mechanical model will bring new opportunities to extend our understanding of a range of basic processes in biology such as molecular recognition and enzyme catalysis. Proteins 2001;44:484–489. © 2001 Wiley‐Liss, Inc.