D3R Grand Challenge 4: prospective pose prediction of BACE1 ligands with AutoDock-GPU

Reactive Docking: A Computational Method for High-Throughput Virtual Screenings of Reactive Species

We describe the formalization of the reactive docking protocol, a method developed to model and predict reactions between small molecules and biological macromolecules. The method has been successfully used in a number of applications already, including recapitulating large proteomics data sets, performing structure-reactivity target optimizations, and prospective virtual screenings. By modeling a near-attack conformation-like state, no QM calculations are required to model the ligand and receptor geometries. Here, we present its generalization using a large data set containing more than 400 ligand-target complexes, 8 nucleophilic modifiable residue types, and more than 30 warheads. The method correctly predicts the modified residue in ∼85% of complexes and shows enrichments comparable to standard focused virtual screenings in ranking ligands. This performance supports this approach for the docking and screening of reactive ligands in virtual chemoproteomics and drug design campaigns.

Brief literature review and comprehensive bioinformatics analytics unravel the potential mechanism of curcumin in the treatment of periodontitis

OBJECTIVE

Periodontitis is a chronic oral disease prevalent worldwide, and natural products are recommended as adjunctive therapy due to their minor side effects. Curcumin, a widely used ancient compound, has been reported to possess therapeutic effects in periodontitis. However, the exact mechanism underlying its activity remains unclear. In this context, the present study aimed to conduct computational simulations to uncover the potential mechanism of action of Curcumin in the treatment of periodontitis.

MATERIALS AND METHODS

Single-cell analysis was conducted using a dataset (i.e., GSE164241) curated from the Gene Expression Omnibus (GEO) database through an R package "Seurat package." Bulk RNA sequencing data were curated from GSE10334 and GSE16134 and processed by R package "Limma." Then, the marker genes in the single-cell transcriptome and differentially expressed genes (DEGs) in the bulk transcriptome were integrated. Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analyses were also carried out to reveal their functionalities. Key targets were mined from their protein-protein interaction (PPI) network topologically. Afterward, molecular docking was performed. The top-ranked pose was subjected to molecular dynamics simulations to investigate the stability of the docking result.

RESULTS

FOS, CXCL1, CXCL8, and IL1B, were filtered after a series of selected processes. The results of molecular modeling suggested that except for IL1B, the Vena Scores of the rest exceeded -5 kcal/mol. Furthermore, the molecular dynamic simulation indicated that the binding of the CXCL8-Curcumin complex was stable over the entire 100 ns simulation.

CONCLUSION

The present study unlocked the binding modes of CXCL1, FOS, and CXCL8 with the Curcumin molecule, which were relatively stable, especially for CXCL8, hindering its promising potential to serve as the critical targets of Curcumin in periodontitis treatment.

iORandLigandDB: A Website for Three-Dimensional Structure Prediction of Insect Odorant Receptors and Docking with Odorants

The use of insect-specific odorants to control the behavior of insects has always been a hot spot in research on "green" control strategies of insects. However, it is generally time-consuming and laborious to explore insect-specific odorants with traditional reverse chemical ecology methods. Here, an insect odorant receptor (OR) and ligand database website (iORandLigandDB) was developed for the specific exploration of insect-specific odorants by using deep learning algorithms. The website provides a range of specific odorants before molecular biology experiments as well as the properties of ORs in closely related insects. At present, the existing three-dimensional structures of ORs in insects and the docking data with related odorants can be retrieved from the database and further analyzed.

Recent breakthroughs in computational structural biology harnessing the power of sequences and structures

Recent advances in computational approaches and their integration into structural biology enable tackling increasingly complex questions. Here, we discuss several key areas, highlighting breakthroughs and remaining challenges. Theoretical modeling has provided tools to accurately predict and design protein structures on a scale currently difficult to achieve using experimental approaches. Molecular Dynamics simulations have become faster and more precise, delivering actionable information inaccessible by current experimental methods. Virtual screening workflows allow a high-throughput approach to discover ligands that bind and modulate protein function, while Machine Learning methods enable the design of proteins with new functionalities. Integrative structural biology combines several of these approaches, pushing the frontiers of structural and functional characterization to ever larger systems, advancing towards a complete understanding of the living cell. These breakthroughs will accelerate and significantly impact diverse areas of science.

Accelerating AutoDock Vina with GPUs

AutoDock Vina is one of the most popular molecular docking tools. In the latest benchmark CASF-2016 for comparative assessment of scoring functions, AutoDock Vina won the best docking power among all the docking tools. Modern drug discovery is facing a common scenario of large virtual screening of drug hits from huge compound databases. Due to the seriality characteristic of the AutoDock Vina algorithm, there is no successful report on its parallel acceleration with GPUs. Current acceleration of AutoDock Vina typically relies on the stack of computing power as well as the allocation of resource and tasks, such as the VirtualFlow platform. The vast resource expenditure and the high access threshold of users will greatly limit the popularity of AutoDock Vina and the flexibility of its usage in modern drug discovery. In this work, we proposed a new method, Vina-GPU, for accelerating AutoDock Vina with GPUs, which is greatly needed for reducing the investment for large virtual screens and also for wider application in large-scale virtual screening on personal computers, station servers or cloud computing, etc. Our proposed method is based on a modified Monte Carlo using simulating annealing AI algorithm. It greatly raises the number of initial random conformations and reduces the search depth of each thread. Moreover, a classic optimizer named BFGS is adopted to optimize the ligand conformations during the docking progress, before a heterogeneous OpenCL implementation was developed to realize its parallel acceleration leveraging thousands of GPU cores. Large benchmark tests show that Vina-GPU reaches an average of 21-fold and a maximum of 50-fold docking acceleration against the original AutoDock Vina while ensuring their comparable docking accuracy, indicating its potential for pushing the popularization of AutoDock Vina in large virtual screens.

Benchmarking the Performance of Irregular Computations in AutoDock-GPU Molecular Docking

Irregular applications can be found in different scientific fields. In computer-aided drug design, molecular docking simulations play an important role in finding promising drug candidates. AutoDock is a software application widely used for predicting molecular interactions at close distances. It is characterized by irregular computations and long execution runtimes. In recent years, a hardware-accelerated version of AutoDock, called AutoDock-GPU, has been under active development. This work benchmarks the recent code and algorithmic enhancements incorporated into AutoDock-GPU. Particularly, we analyze the impact on execution runtime of techniques based on early termination. These enable AutoDock-GPU to explore the molecular space as necessary, while safely avoiding redundant computations. Our results indicate that it is possible to achieve average runtime reductions of 50% by using these techniques. Furthermore, a comprehensive literature review is also provided, where our work is compared to relevant approaches leveraging hardware acceleration for molecular docking.

Performance evaluation of flexible macrocycle docking in AutoDock

Macrocycles represent an important class of ligands, both in natural products and designed drugs. In drug design, macrocyclizations can impart specific ligand conformations and contribute to passive permeation by encouraging intramolecular H-bonds. AutoDock-GPU and Vina can model macrocyclic ligands flexibly, without requiring the enumeration of macrocyclic conformers before docking. Here, we characterize the performance of the method for handling macrocyclic compounds, which is implemented and the default behaviour for ligand preparation with our ligand preparation pipeline, Meeko. A pseudoatom is used to encode bond geometry and produce an anisotropic closure force for macrocyclic rings. This method is evaluated on a diverse set of small molecule and peptide macrocycles, ranging from 7- to 33-membered rings, showing little accuracy loss compared to rigid redocking of the X-ray macrocycle conformers. This suggests that for conformationally flexible macrocycles with unknown binding modes, this method can be effectively used to predict the macrocycle conformation.

AutoDock Vina 1.2.0: New Docking Methods, Expanded Force Field, and Python Bindings

AutoDock Vina is arguably one of the fastest and most widely used open-source programs for molecular docking. However, compared to other programs in the AutoDock Suite, it lacks support for modeling specific features such as macrocycles or explicit water molecules. Here, we describe the implementation of this functionality in AutoDock Vina 1.2.0. Additionally, AutoDock Vina 1.2.0 supports the AutoDock4.2 scoring function, simultaneous docking of multiple ligands, and a batch mode for docking a large number of ligands. Furthermore, we implemented Python bindings to facilitate scripting and the development of docking workflows. This work is an effort toward the unification of the features of the AutoDock4 and AutoDock Vina programs. The source code is available at https://github.com/ccsb-scripps/AutoDock-Vina.

Accelerating AutoDock4 with GPUs and Gradient-Based Local Search

AutoDock4 is a widely used program for docking small molecules to macromolecular targets. It describes ligand-receptor interactions using a physics-inspired scoring function that has been proven useful in a variety of drug discovery projects. However, compared to more modern and recent software, AutoDock4 has longer execution times, limiting its applicability to large scale dockings. To address this problem, we describe an OpenCL implementation of AutoDock4, called AutoDock-GPU, that leverages the highly parallel architecture of GPU hardware to reduce docking runtime by up to 350-fold with respect to a single-threaded process. Moreover, we introduce the gradient-based local search method ADADELTA, as well as an improved version of the Solis-Wets random optimizer from AutoDock4. These efficient local search algorithms significantly reduce the number of calls to the scoring function that are needed to produce good results. The improvements reported here, both in terms of docking throughput and search efficiency, facilitate the use of the AutoDock4 scoring function in large scale virtual screening.

Supercomputing Pipelines Search for Therapeutics Against COVID-19

The urgent search for drugs to combat SARS-CoV-2 has included the use of supercomputers. The use of general-purpose graphical processing units (GPUs), massive parallelism, and new software for high-performance computing (HPC) has allowed researchers to search the vast chemical space of potential drugs faster than ever before. We developed a new drug discovery pipeline using the Summit supercomputer at Oak Ridge National Laboratory to help pioneer this effort, with new platforms that incorporate GPU-accelerated simulation and allow for the virtual screening of billions of potential drug compounds in days compared to weeks or months for their ability to inhibit SARS-COV-2 proteins. This effort will accelerate the process of developing drugs to combat the current COVID-19 pandemic and other diseases.

The AutoDock suite at 30

The AutoDock suite provides a comprehensive toolset for computational ligand docking and drug design and development. The suite builds on 30 years of methods development, including empirical free energy force fields, docking engines, methods for site prediction, and interactive tools for visualization and analysis. Specialized tools are available for challenging systems, including covalent inhibitors, peptides, compounds with macrocycles, systems where ordered hydration plays a key role, and systems with substantial receptor flexibility. All methods in the AutoDock suite are freely available for use and reuse, which has engendered the continued growth of a diverse community of primary users and third-party developers.