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×.
We present a simulation study of how properties of symmetric diblock copolymers depend on the invariant degree of polymerization N̅, focusing on the vicinity of the order–disorder transition (ODT). ...Results from several coarse-grained simulation models are combined to cover a range of N̅ ≃ 200–8000 that includes most of the experimentally relevant range. Results are presented for the free energy per chain, the value of χe N at the ODT, the latent heat of transition, the layer spacing, the composition profile, and compression modulus in the ordered phase. Universality (i.e., model independence) is demonstrated by showing that equivalent results for all properties are obtained from corresponding thermodynamic states of different simulation models. Corresponding states of symmetric copolymers are states with equal values of the parameters χe N and N̅, where χe is an effective Flory–Huggins interaction parameter and N is a degree of polymerization. The underlying universality becomes apparent, however, only if data are analyzed using an adequate estimate of χe, which we obtain by fitting the structure factor in the disordered state to recent theoretical predictions. The results show that behavior near the ODT exhibits a different character at moderate and high values of N̅, with a crossover near N̅ ≃ 104. Within the range N̅ ≲ 104 studied here, the ordered and disordered phases near the ODT both contain strongly segregated domains of nearly pure A and B, in contrast to the assumption of weak segregation underlying the Fredrickson–Helfand (FH) theory. In this regime, the FH theory is inaccurate and substantially underestimates the value of χe N at the ODT. Results for the highest values of N̅ studied here agree reasonably well with FH predictions, suggesting that the theory may be accurate for N̅ ≳ 104. Self-consistent field theory (SCFT) grossly underestimates (χe N)ODT for modest N̅ because it cannot describe strong correlations in the disordered phase. SCFT is found, however, to yield accurate predictions for several properties of the ordered lamellar phase.
Depletion interactions arise from entropic forces, and their ability to induce aggregation and even ordering of colloidal particles through self-assembly is well established, especially for spherical ...colloids. We vary the size and concentration of penetrable hard sphere depletants in a system of cuboctahedra, and we show how depletion changes the preferential facet alignment of the colloids and thereby selects different crystal structures. Moreover, we explain the cuboctahedra phase behavior using perturbative free energy calculations. We find that cuboctahedra can form a stable simple cubic phase, and, remarkably, that the stability of this phase can be rationalized only by considering the effects of both the colloid and depletant entropy. We corroborate our results by analyzing how the depletant concentration and size affect the emergent directional entropic forces and hence the effective particle shape. We propose the use of depletants as a means of easily changing the effective shape of self-assembling anisotropic colloids.
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, Letnik:
10, Številka:
1
Journal Article
Recenzirano
Odprti dostop
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.
Shape guides colloidal nanoparticles to form complex assemblies, but its role in defining interfaces in biomolecular complexes is less clear. In this work, we isolate the role of shape in protein ...complexes by studying the reversible binding processes of 46 protein dimer pairs, and investigate when entropic effects from shape complementarity alone are sufficient to predict the native protein binding interface. We employ depletants using a generic, implicit depletion model to amplify the magnitude of the entropic forces arising from lock-and-key binding and isolate the effect of shape complementarity in protein dimerization. For 13% of the complexes studied here, protein shape is sufficient to predict native complexes as equilibrium assemblies. We elucidate the results by analyzing the importance of competing binding configurations and how it affects the assembly. Furthermore, a machine learning classifier, with a precision of 89.14% and a recall of 77.11%, is able to identify the cases where shape alone predicts the native protein interface.
Shape guides colloidal nanoparticles to form complex assemblies, but its role in defining interfaces in biomolecular complexes is less clear. In this work, we isolate the role of shape in protein ...complexes by studying the reversible binding processes of 46 protein dimer pairs, and investigate when entropic effects from shape complementarity alone are sufficient to predict the native protein binding interface. We employ depletants using a generic, implicit depletion model to amplify the magnitude of the entropic forces arising from lock-and-key binding and isolate the effect of shape complementarity in protein dimerization. For 13% of the complexes studied here, protein shape is sufficient to predict native complexes as equilibrium assemblies. We elucidate the results by analyzing the importance of competing binding configurations and how it affects the assembly. A machine learning classifier, with a precision of 89.14% and a recall of 77.11%, is able to identify the cases where shape alone predicts the native protein interface.
Shape alone guides colloidal nanoparticles to form complex assemblies, and in some cases can define interfaces in biomolecular complexes.