Display omitted
HOOMD-blue is a particle simulation engine designed for nano- and colloidal-scale molecular dynamics and hard particle Monte Carlo simulations. It has been actively developed since ...March 2007 and available open source since August 2008. HOOMD-blue is a Python package with a high performance C++/CUDA backend that we built from the ground up for GPU acceleration. The Python interface allows users to combine HOOMD-blue with other packages in the Python ecosystem to create simulation and analysis workflows. We employ software engineering practices to develop, test, maintain, and expand the code.
We describe a highly optimized implementation of MPI domain decomposition in a GPU-enabled, general-purpose molecular dynamics code, HOOMD-blue (Anderson and Glotzer, 2013). Our approach is inspired ...by a traditional CPU-based code, LAMMPS (Plimpton, 1995), but is implemented within a code that was designed for execution on GPUs from the start (Anderson et al., 2008). The software supports short-ranged pair force and bond force fields and achieves optimal GPU performance using an autotuning algorithm. We are able to demonstrate equivalent or superior scaling on up to 3375 GPUs in Lennard-Jones and dissipative particle dynamics (DPD) simulations of up to 108 million particles. GPUDirect RDMA capabilities in recent GPU generations provide better performance in full double precision calculations. For a representative polymer physics application, HOOMD-blue 1.0 provides an effective GPU vs. CPU node speed-up of 12.5×.
The vast size of chemical space necessitates computational approaches to automate and accelerate the design of molecular sequences to guide experimental efforts for drug discovery. Genetic algorithms ...provide a useful framework to incrementally generate molecules by applying mutations to known chemical structures. Recently, masked language models have been applied to automate the mutation process by leveraging large compound libraries to learn commonly occurring chemical sequences (i.e., using tokenization) and predict rearrangements (i.e., using mask prediction). Here, we consider how language models can be adapted to improve molecule generation for different optimization tasks. We use two different generation strategies for comparison, fixed and adaptive. The fixed strategy uses a pre-trained model to generate mutations; the adaptive strategy trains the language model on each new generation of molecules selected for target properties during optimization. Our results show that the adaptive strategy allows the language model to more closely fit the distribution of molecules in the population. Therefore, for enhanced fitness optimization, we suggest the use of the fixed strategy during an initial phase followed by the use of the adaptive strategy. We demonstrate the impact of adaptive training by searching for molecules that optimize both heuristic metrics, drug-likeness and synthesizability, as well as predicted protein binding affinity from a surrogate model. Our results show that the adaptive strategy provides a significant improvement in fitness optimization compared to the fixed pre-trained model, empowering the application of language models to molecular design tasks.
The unique mechanical performance of animal cells and tissues is attributed mostly to their internal biopolymer meshworks. Its perplexing universality and robustness against structural modifications ...by drugs and mutations is an enigma in cell biology and provides formidable challenges to materials science. Recent investigations could pinpoint highly universal patterns in the soft glassy rheology and nonlinear elasticity of cells and reconstituted networks. Here, we report observations of a glass transition in semidilute F-actin solutions, which could hold the key to a unified explanation of these phenomena. Combining suitable rheological protocols with high-precision dynamic light scattering, we can establish a remarkable rheological redundancy and trace it back to a highly universal exponential stretching of the single-polymer relaxation spectrum of a "glassy wormlike chain." By exploiting the ensuing generalized time-temperature superposition principle, the time domain accessible to microrheometry can be extended by several orders of magnitude, thus opening promising new metrological opportunities.
SARS-CoV2 billion-compound docking Rogers, David M; Agarwal, Rupesh; Vermaas, Josh V ...
Scientific data,
03/2023, Volume:
10, Issue:
1
Journal Article
Peer reviewed
Open access
This dataset contains ligand conformations and docking scores for 1.4 billion molecules docked against 6 structural targets from SARS-CoV2, representing 5 unique proteins: MPro, NSP15, PLPro, RDRP, ...and the Spike protein. Docking was carried out using the AutoDock-GPU platform on the Summit supercomputer and Google Cloud. The docking procedure employed the Solis Wets search method to generate 20 independent ligand binding poses per compound. Each compound geometry was scored using the AutoDock free energy estimate, and rescored using RFScore v3 and DUD-E machine-learned rescoring models. Input protein structures are included, suitable for use by AutoDock-GPU and other docking programs. As the result of an exceptionally large docking campaign, this dataset represents a valuable resource for discovering trends across small molecule and protein binding sites, training AI models, and comparing to inhibitor compounds targeting SARS-CoV-2. The work also gives an example of how to organize and process data from ultra-large docking screens.
Display omitted
We present a detailed derivation of the expression for the pressure in MD simulations that contain rigid bodies, where two equivalent formulations have been developed. One of these ...formulations was used in HOOMD-blue v1.x, but implemented incorrectly. We point out the precise reason for this implementation issue, the difference with the current and correct implementation in HOOMD-blue v2.x, and lessons learned. We perform numerical validation tests using dumbbell models, a mixture of cubic and spherical particles, and the SPC/E water model.