Optimization of parallel implementation of UNRES package for coarse‐grained simulations to treat large proteins
Adam K. Sieradzan,Jordi Sans‐Duñó,Emilia A. Lubecka,Cezary Czaplewski,Agnieszka G. Lipska,Henryk Leszczyński,Krzysztof M. Ocetkiewicz,Jerzy Proficz,Paweł Czarnul,Henryk Krawczyk,Adam Liwo
DOI: https://doi.org/10.1002/jcc.27026
2022-11-17
Journal of Computational Chemistry
Abstract:Molecular dynamics with the optimized and efficiently parallelized implementation of the highly‐reduced physics‐based UNRES model of polypeptide chains enables us to reach, with moderate computer resources, several nanosecond/day of MD time, for protein systems with size over 100,000 residues, which translates to microsecond/day real time given a 1000 fold faster occurrence of events in the simulations with the UNRES model compared to the all‐atom representation. We report major algorithmic improvements of the UNRES package for physics‐based coarse‐grained simulations of proteins. These include (i) introduction of interaction lists to optimize computations, (ii) transforming the inertia matrix to a pentadiagonal form to reduce computing and memory requirements, (iii) removing explicit angles and dihedral angles from energy expressions and recoding the most time‐consuming energy/force terms to minimize the number of operations and to improve numerical stability, (iv) using OpenMP to parallelize those sections of the code for which distributed‐memory parallelization involves unfavorable computing/communication time ratio, and (v) careful memory management to minimize simultaneous access of distant memory sections. The new code enables us to run molecular dynamics simulations of protein systems with size exceeding 100,000 amino‐acid residues, reaching over 1 ns/day (1 μs/day in all‐atom timescale) with 24 cores for proteins of this size. Parallel performance of the code and comparison of its performance with that of AMBER, GROMACS and MARTINI 3 is presented.
chemistry, multidisciplinary