Abstract:The Grand Canonical Monte Carlo (GCMC) ensemble defined by the excess chemical potential, μex , volume, and temperature, in the context of molecular simulations allows for variations in the number of particles in the system. In practice, GCMC simulations have been widely applied for the sampling of rare gasses and water, but limited in the context of larger molecules. To overcome this limitation, the oscillating μex GCMC method was introduced and shown to be of utility for sampling small solutes, such as formamide, propane, and benzene, as well as for ionic species such as monocations, acetate, and methylammonium. However, the acceptance of GCMC insertions is low, and the method is computationally demanding. In the present study, we improved the sampling efficiency of the GCMC method using known cavity-bias and configurational-bias algorithms in the context of GPU architecture. Specifically, for GCMC simulations of aqueous solution systems, the configurational-bias algorithm was extended by applying system partitioning in conjunction with a random interval extraction algorithm, thereby improving the efficiency in a highly parallel computing environment. The method is parallelized on the GPU using CUDA and OpenCL, allowing for the code to run on both Nvidia and AMD GPUs, respectively. Notably, the method is particularly well suited for GPU computing as the large number of threads allows for simultaneous sampling of a large number of configurations during insertion attempts without additional computational overhead. In addition, the partitioning scheme allows for simultaneous insertion attempts for large systems, offering considerable efficiency. Calculations on the BK Channel, a transporter, including a lipid bilayer with over 760,000 atoms, show a speed up of ~53-fold through the use of system partitioning. The improved algorithm is then combined with an enhanced μex oscillation protocol and shown to be of utility in the context of the site-identification by ligand competitive saturation (SILCS) co-solvent sampling approach as illustrated through application to the protein CDK2.

CUDASW++4.0: ultra-fast GPU-based Smith–Waterman protein sequence database search

GPU Accelerated Biological Sequence Alignment

CUDAMPF++: A Proactive Resource Exhaustion Scheme for Accelerating Homologous Sequence Search on CUDA-enabled GPU

Accelerating Biological Sequence Alignment Algorithm on GPU with CUDA

GPU-accelerated homology search with MMseqs2

Fast Alignment of Biological Sequences Based on General Purpose Graphic Processor Unit

Efficient Distributed Smith-Waterman Algorithm Based on Apache Spark

gpuPairHMM: High-speed Pair-HMM Forward Algorithm for DNA Variant Calling on GPUs

MetaCache-GPU: Ultra-Fast Metagenomic Classification

Fast and Exact Sequence Alignment with the Smith-Waterman Algorithm: The SwissAlign Webserver

An Efficient Parallel Algorithm For Accelerating Computational Protein Design

Sparksw: Scalable Distributed Computing System For Large-Scale Biological Sequence Alignment

Fast and accurate protein structure search with Foldseek

AGAThA: Fast and Efficient GPU Acceleration of Guided Sequence Alignment for Long Read Mapping

Hardware Acceleration for the Banded Smith-Waterman Algorithm with the Cycled Systolic Array

hLCS. A Hybrid GPGPU Approach for Solving Multiple Short and Unbalanced LCS Problems

Rapid multiple protein sequence search by parallel and heterogeneous computation

High-Throughput Sequence Translation Using Cuda

Implementation of the Smith-Waterman Algorithm on a Reconfigurable Supercomputing Platform

GPU-specific algorithms for improved solute sampling in grand canonical Monte Carlo simulations

Fast Calculation of Protein Molecular Field and Its Application on Protein Recog-nition Research??