Premium
Large‐scale molecular dynamics simulations of HLA‐A*0201 complexed with a tumor‐specific antigenic peptide: Can the α3 and β 2 m domains be neglected?
Author(s) -
Wan Shunzhou,
Coveney Peter,
Flower Darren R.
Publication year - 2004
Publication title -
journal of computational chemistry
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.907
H-Index - 188
eISSN - 1096-987X
pISSN - 0192-8651
DOI - 10.1002/jcc.20100
Subject(s) - molecular dynamics , peptide , massively parallel , domain (mathematical analysis) , scalability , major histocompatibility complex , human leukocyte antigen , computer science , scale (ratio) , chemistry , computational science , physics , parallel computing , antigen , computational chemistry , biology , database , biochemistry , mathematics , genetics , gene , mathematical analysis , quantum mechanics
Large‐scale massively parallel molecular dynamics (MD) simulations of the human class I major histocompatibility complex (MHC) protein HLA‐A*0201 bound to a decameric tumor‐specific antigenic peptide GVYDGREHTV were performed using a scalable MD code on high‐performance computing platforms. Such computational capabilities put us in reach of simulations of various scales and complexities. The supercomputing resources available for this study allow us to compare directly differences in the behavior of very large molecular models; in this case, the entire extracellular portion of the peptide–MHC complex vs. the isolated peptide binding domain. Comparison of the results from the partial and the whole system simulations indicates that the peptide is less tightly bound in the partial system than in the whole system. From a detailed study of conformations, solvent‐accessible surface area, the nature of the water network structure, and the binding energies, we conclude that, when considering the conformation of the α1–α2 domain, the α3 and β 2 m domains cannot be neglected. © 2004 Wiley Periodicals, Inc. J Comput Chem 25: 1803–1813, 2004