Development of massive multilevel molecular dynamics simulation program, platypus (PLATform for dYnamic protein unified simulation), for the elucidation of protein functions

A massively parallel program for quantum mechanical‐molecular mechanical (QM/MM) molecular dynamics simulation, called Platypus (PLATform for dYnamic Protein Unified Simulation), was developed to elucidate protein functions. The speedup and the parallelization ratio of Platypus in the QM and QM/MM calculations were assessed for a bacteriochlorophyll dimer in the photosynthetic reaction center (DIMER) on the K computer, a massively parallel computer achieving 10 PetaFLOPs with 705,024 cores. Platypus exhibited the increase in speedup up to 20,000 core processors at the HF/cc‐pVDZ and B3LYP/cc‐pVDZ, and up to 10,000 core processors by the CASCI(16,16)/6‐31G** calculations. We also performed excited QM/MM‐MD simulations on the chromophore of Sirius (SIRIUS) in water. Sirius is a pH‐insensitive and photo‐stable ultramarine fluorescent protein. Platypus accelerated on‐the‐fly excited‐state QM/MM‐MD simulations for SIRIUS in water, using over 4000 core processors. In addition, it also succeeded in 50‐ps (200,000‐step) on‐the‐fly excited‐state QM/MM‐MD simulations for the SIRIUS in water. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

[1]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[2]  Victor Eijkhout,et al.  Introduction to High Performance Scientific Computing , 2015 .

[3]  Takeharu Nagai,et al.  An ultramarine fluorescent protein with increased photostability and pH insensitivity , 2009, Nature Methods.

[4]  T. H. Dunning.  Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[5]  Michael J. Frisch,et al.  MP2 energy evaluation by direct methods , 1988 .

[6]  Toshiyuki Imamura,et al.  High-performance computing for exact numerical approaches to quantum many-body problems on the earth simulator , 2006, SC.

[7]  M. W. van der Kamp,et al.  Combined quantum mechanics/molecular mechanics (QM/MM) methods in computational enzymology. , 2013, Biochemistry.

[8]  Haruki Nakamura,et al.  Intra- and intermolecular interaction inducing pyramidalization on both sides of a proline dipeptide during isomerization: an ab initio QM/MM molecular dynamics simulation study in explicit water. , 2009, Journal of the American Chemical Society.

[9]  J. Pople,et al.  Self‐Consistent Molecular‐Orbital Methods. IX. An Extended Gaussian‐Type Basis for Molecular‐Orbital Studies of Organic Molecules , 1971 .

[10]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

[11]  Michael J. Frisch,et al.  Toward a systematic molecular orbital theory for excited states , 1992 .

[12]  F. Young.  Biochemistry , 1955, The Indian Medical Gazette.

[13]  T. Takada,et al.  A parallel algorithm for generating molecular integrals over MO basis sets , 2002 .

[14]  Timothy G. Mattson.  Parallel Computing in Computational Chemistry , 1995 .

[15]  O. Dym,et al.  Thermal stabilization of the protozoan Entamoeba histolytica alcohol dehydrogenase by a single proline substitution , 2008 .

[16]  B. Roos,et al.  A new method for large-scale Cl calculations , 1972 .

[17]  P. Kollman,et al.  How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? , 2000 .

[18]  Alan H. Karp,et al.  Measuring parallel processor performance , 1990, CACM.

[19]  Walter Thiel,et al.  QM/MM methods for biomolecular systems. , 2009, Angewandte Chemie.

[20]  S. Brereton.  Life , 1876, The Indian medical gazette.

[21]  P. Kollman.  Advances and Continuing Challenges in Achieving Realistic and Predictive Simulations of the Properties of Organic and Biological Molecules , 1996 .

[22]  Kazuya Ishimura,et al.  Massively parallel MP2‐F12 calculations on the K computer , 2015 .

[23]  M. Ratner.  Molecular electronic-structure theory , 2000 .

[24]  Trygve Helgaker,et al.  Molecular Electronic-Structure Theory: Helgaker/Molecular Electronic-Structure Theory , 2000 .

[25]  Naehyuck Chang,et al.  Proceedings of the 17th IEEE/ACM international symposium on Low-power electronics and design , 2011, ISLPED 2011.

[26]  L. Christophorou.  Science , 2018, Emerging Dynamics: Science, Energy, Society and Values.

[27]  Daniel V. Pryor,et al.  Proceedings of the 1990 ACM/IEEE conference on Supercomputing , 1990 .

[28]  A. Becke.  Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[29]  Y. Fukunishi,et al.  The Filling Potential Method: A Method for Estimating the Free Energy Surface for Protein−Ligand Docking , 2003 .

[31]  D. Truhlar,et al.  QM/MM: what have we learned, where are we, and where do we go from here? , 2007 .

[32]  S. Huzinaga,et al.  A systematic preparation of new contracted Gaussian‐type orbital sets. III. Second‐row atoms from Li through ne , 1980 .

[33]  M. Levitt,et al.  Theoretical studies of enzymic reactions: dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme. , 1976, Journal of molecular biology.

[34]  Shamoon Jamshed,et al.  Using HPC for Computational Fluid Dynamics: A Guide to High Performance Computing for CFD Engineers , 2015 .

[35]  Manabu Oumi,et al.  A doubles correction to electronic excited states from configuration interaction in the space of single substitutions , 1994 .

[36]  D C Rees,et al.  Light-induced structural changes in photosynthetic reaction center: implications for mechanism of electron-proton transfer. , 1997, Science.