Abstract
There has been an increasing demand for longer-timescale molecular dynamics (MD) simulations of larger biomolecular systems. To meet these demands, using the C++ API of OpenMM, we developed a fast and flexible MD software, OpenCafeMol, for residue-resolution protein and lipid models that shows high performance on graphics processing unit (GPU) machines. We validated OpenCafeMol for folding small proteins, lipid membrane dynamics, and membrane protein structures. Benchmark tests of the computation times showed that OpenCafeMol with one GPU for proteins and lipid membranes is approximately 100 and 240 times faster than the corresponding simulations on a typical CPU machine (eight cores), respectively. Taking advantage of the high speed of OpenCafeMol, we applied it to two sets of vesicle fusion simulations; one driven by force and the other coupled with conformational dynamics of a soluble N-ethylmaleimide-sensitive factor attachment protein receptor (SNARE) complex. In the latter MD simulation at a high temperature resulted in vesicle docking and pore formation, followed by fusion, which are coupled with local folding of linkers in the SNARE complex. This opens up a new avenue to study membrane-fusion mechanisms via MD simulations. The source code for OpenCafeMol is fully available.