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.

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.

The cis-trans isomerization of the peptide bond preceding a proline plays important roles in protein folding and biological function. Although many experimental and theoretical studies have been done, the mechanism has not yet been clearly elucidated. We studied the cis-trans isomerization of the proline dipeptide (Ace-Pro-NMe) in explicit water by molecular dynamics simulations using a combined potential derived from ab initio quantum mechanics and empirical molecular mechanics. We obtained the free energy landscape during the isomerization by using the umbrella sampling method. The free energy landscape is in good accordance with previous experimental and theoretical values. We observed that in the middle of the isomerization, the prolyl nitrogen transiently takes pyramidal conformations in two polarized directions and that, simultaneously, the prolyl C-N bond extends. We show that these geometrical changes cooperatively transform the prolyl nitrogen from a sp(2)-hybridized electronic state into a sp(3)-hybridized one, and thus realize a transition state that reduces the rotational barriers separating the cis- and trans-states. We also found that the hydration of the prolyl nitrogen stabilizes the negative pyramidal conformation, while an intramolecular interaction mainly stabilizes the positive one. Fluctuations in the polarity and magnitude of the pyramidal conformation during the isomerization are interpreted as a competition between the hydrogen-bonding partners for the prolyl nitrogen between different sides of the pyrrolidine ring.