GWOVina: A grey wolf optimization approach to rigid and flexible receptor docking

Ergothioneine boosts mitochondrial respiration and exercise performance via direct activation of MPST

Ergothioneine (EGT) is a diet-derived, atypical amino acid that accumulates to high levels in human tissues. Reduced EGT levels have been linked to age-related disorders, including neurodegenerative and cardiovascular diseases, while EGT supplementation is protective in a broad range of disease and aging models in mice. Despite these promising data, the direct and physiologically relevant molecular target of EGT has remained elusive. Here we use a systematic approach to identify how mitochondria remodel their metabolome in response to exercise training. From this data, we find that EGT accumulates in muscle mitochondria upon exercise training. Proteome-wide thermal stability studies identify 3-mercaptopyruvate sulfurtransferase (MPST) as a direct molecular target of EGT; EGT binds to and activates MPST, thereby boosting mitochondrial respiration and exercise training performance in mice. Together, these data identify the first physiologically relevant EGT target and establish the EGT-MPST axis as a molecular mechanism for regulating mitochondrial function and exercise performance.

SCORCH: Improving structure-based virtual screening with machine learning classifiers, data augmentation, and uncertainty estimation

INTRODUCTION

The discovery of a new drug is a costly and lengthy endeavour. The computational prediction of which small molecules can bind to a protein target can accelerate this process if the predictions are fast and accurate enough. Recent machine-learning scoring functions re-evaluate the output of molecular docking to achieve more accurate predictions. However, previous scoring functions were trained on crystalised protein-ligand complexes and datasets of decoys. The limited availability of crystal structures and biases in the decoy datasets can lower the performance of scoring functions.

OBJECTIVES

To address key limitations of previous scoring functions and thus improve the predictive performance of structure-based virtual screening.

METHODS

A novel machine-learning scoring function was created, named SCORCH (Scoring COnsensus for RMSD-based Classification of Hits). To develop SCORCH, training data is augmented by considering multiple ligand poses and labelling poses based on their RMSD from the native pose. Decoy bias is addressed by generating property-matched decoys for each ligand and using the same methodology for preparing and docking decoys and ligands. A consensus of 3 different machine learning approaches is also used to improve performance.

RESULTS

We find that multi-pose augmentation in SCORCH improves its docking power and screening power on independent benchmark datasets. SCORCH outperforms an equivalent scoring function trained on single poses, with a 1% enrichment factor (EF) of 13.78 vs. 10.86 on 18 DEKOIS 2.0 targets and a mean native pose rank of 5.9 vs 30.4 on CSAR 2014. Additionally, SCORCH outperforms widely used scoring functions in virtual screening and pose prediction on independent benchmark datasets.

CONCLUSION

By rationally addressing key limitations of previous scoring functions, SCORCH improves the performance of virtual screening. SCORCH also provides an estimate of its uncertainty, which can help reduce the cost and time required for drug discovery.

Structure Prediction, Evaluation, and Validation of GPR18 Lipid Receptor Using Free Programs

The GPR18 receptor, often referred to as the N-arachidonylglycine receptor, although assigned (along with GPR55 and GPR119) to the new class A GPCR subfamily-lipid receptors, officially still has the status of a class A GPCR orphan. While its signaling pathways and biological significance have not yet been fully elucidated, increasing evidence points to the therapeutic potential of GPR18 in relation to immune, neurodegenerative, and cancer processes to name a few. Therefore, it is necessary to understand the interactions of potential ligands with the receptor and the influence of particular structural elements on their activity. Thus, given the lack of an experimentally solved structure, the goal of the present study was to obtain a homology model of the GPR18 receptor in the inactive state, meeting all requirements in terms of protein structure quality and recognition of active ligands. To increase the reliability and precision of the predictions, different contemporary protein structure prediction methods and software were used and compared herein. To test the usability of the resulting models, we optimized and compared the selected structures followed by the assessment of the ability to recognize known, active ligands. The stability of the predicted poses was then evaluated by means of molecular dynamics simulations. On the other hand, most of the best-ranking contemporary CADD software/platforms for its full usability require rather expensive licenses. To overcome this down-to-earth obstacle, the overarching goal of these studies was to test whether it is possible to perform the thorough CADD experiments with high scientific confidence while using only license-free/academic software and online platforms. The obtained results indicate that a wide range of freely available software and/or academic licenses allow us to carry out meaningful molecular modelling/docking studies.

Parallel multi-swarm cooperative particle swarm optimization for protein-ligand docking and virtual screening

BACKGROUND

A high-quality docking method tends to yield multifold gains with half pains for the new drug development. Over the past few decades, great efforts have been made for the development of novel docking programs with great efficiency and intriguing accuracy. AutoDock Vina (Vina) is one of these achievements with improved speed and accuracy compared to AutoDock4. Since it was proposed, some of its variants, such as PSOVina and GWOVina, have also been developed. However, for all these docking programs, there is still large room for performance improvement.

RESULTS

In this work, we propose a parallel multi-swarm cooperative particle swarm model, in which one master swarm and several slave swarms mutually cooperate and co-evolve. Our experiments show that multi-swarm programs possess better docking robustness than PSOVina. Moreover, the multi-swarm program based on random drift PSO can achieve the best highest accuracy of protein-ligand docking, an outstanding enrichment effect for drug-like activate compounds, and the second best AUC screening accuracy among all the compared docking programs, but with less computation consumption than most of the other docking programs.

CONCLUSION

The proposed multi-swarm cooperative model is a novel algorithmic modeling suitable for protein-ligand docking and virtual screening. Owing to the existing coevolution between the master and the slave swarms, this model in parallel generates remarkable docking performance. The source code can be freely downloaded from https://github.com/li-jin-xing/MPSOVina .

Evaluation of Potential Anti-Hepatitis A Virus 3C Protease Inhibitors Using Molecular Docking

Hepatitis A virus (HAV) infection is a major cause of acute hepatitis worldwide and occasionally causes acute liver failure and can lead to death in the absence of liver transplantation. Although HAV vaccination is available, the prevalence of HAV vaccination is not adequate in some countries. Additionally, the improvements in public health reduced our immunity to HAV infection. These situations motivated us to develop potentially new anti-HAV therapeutic options. We carried out the in silico screening of anti-HAV compounds targeting the 3C protease enzyme using the Schrodinger Modeling software from the antiviral library of 25,000 compounds to evaluate anti-HAV 3C protease inhibitors. Additionally, in vitro studies were introduced to examine the inhibitory effects of HAV subgenomic replicon replication and HAV HA11-1299 genotype IIIA replication in hepatoma cell lines using luciferase assays and real-time RT-PCR. In silico studies enabled us to identify five lead candidates with optimal binding interactions in the active site of the target HAV 3C protease using the Schrodinger Glide program. In vitro studies substantiated our hypothesis from in silico findings. One of our lead compounds, Z10325150, showed 47% inhibitory effects on HAV genotype IB subgenomic replicon replication and 36% inhibitory effects on HAV genotype IIIA HA11-1299 replication in human hepatoma cell lines, with no cytotoxic effects at concentrations of 100 μg/mL. The effects of the combination therapy of Z10325150 and RNA-dependent RNA polymerase inhibitor, favipiravir on HAV genotype IB HM175 subgenomic replicon replication and HAV genotype IIIA HA11-1299 replication showed 64% and 48% inhibitory effects of HAV subgenomic replicon and HAV replication, respectively. We identified the HAV 3C protease inhibitor Z10325150 through in silico screening and confirmed the HAV replication inhibitory activity in human hepatocytes. Z10325150 may offer the potential for a useful HAV inhibitor in severe hepatitis A.

Computational methods to assist in the discovery of pharmacological chaperones for rare diseases

Pharmacological chaperones are chemical compounds able to bind proteins and stabilize them against denaturation and following degradation. Some pharmacological chaperones have been approved, or are under investigation, for the treatment of rare inborn errors of metabolism, caused by genetic mutations that often can destabilize the structure of the wild-type proteins expressed by that gene. Given that, for rare diseases, there is a general lack of pharmacological treatments, many expectations are poured out on this type of compounds. However, their discovery is not straightforward. In this review, we would like to focus on the computational methods that can assist and accelerate the search for these compounds, showing also examples in which these methods were successfully applied for the discovery of promising molecules belonging to this new category of pharmacologically active compounds.

MSLDOCK: Multi-Swarm Optimization for Flexible Ligand Docking and Virtual Screening

Autodock and its various variants are widely utilized docking approaches, which adopt optimization methods as search algorithms for flexible ligand docking and virtual screening. However, many of them have their limitations, such as poor accuracy for dockings with highly flexible ligands and low docking efficiency. In this paper, a multi-swarm optimization algorithm integrated with Autodock environment is proposed to design a high-performance and high-efficiency docking program, namely, MSLDOCK. The search algorithm is a combination of the random drift particle swarm optimization with a novel multi-swarm strategy and the Solis and Wets local search method with a modified implementation. Due to the algorithm's structure, MSLDOCK also has a multithread mode. The experimental results reveal that MSLDOCK outperforms other two Autodock-based approaches in many aspects, such as self-docking, cross-docking, and virtual screening accuracies as well as docking efficiency. Moreover, compared with three non-Autodock-based docking programs, MSLDOCK can be a reliable choice for self-docking and virtual screening, especially for dealing with highly flexible ligand docking problems. The source code of MSLDOCK can be downloaded for free from https://github.com/lcmeteor/MSLDOCK.

A multi-pronged approach targeting SARS-CoV-2 proteins using ultra-large virtual screening

The unparalleled global effort to combat the continuing severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic over the last year has resulted in promising prophylactic measures. However, a need still exists for cheap, effective therapeutics, and targeting multiple points in the viral life cycle could help tackle the current, as well as future, coronaviruses. Here, we leverage our recently developed, ultra-large-scale in silico screening platform, VirtualFlow, to search for inhibitors that target SARS-CoV-2. In this unprecedented structure-based virtual campaign, we screened roughly 1 billion molecules against each of 40 different target sites on 17 different potential viral and host targets. In addition to targeting the active sites of viral enzymes, we also targeted critical auxiliary sites such as functionally important protein-protein interactions.

Accounting of Receptor Flexibility in Ultra-Large Virtual Screens with VirtualFlow Using a Grey Wolf Optimization Method

Structure-based virtual screening approaches have the ability to dramatically reduce the time and costs associated to the discovery of new drug candidates. Studies have shown that the true hit rate of virtual screenings improves with the scale of the screened ligand libraries. Therefore, we have recently developed an open source drug discovery platform (VirtualFlow), which is able to routinely carry out ultra-large virtual screenings. One of the primary challenges of molecular docking is the circumstance when the protein is highly dynamic or when the structure of the protein cannot be captured by a static pose. To accommodate protein dynamics, we report the extension of VirtualFlow to allow the docking of ligands using a grey wolf optimization algorithm using the docking program GWOVina, which substantially improves the quality and efficiency of flexible receptor docking compared to AutoDock Vina. We demonstrate the linear scaling behavior of VirtualFlow utilizing GWOVina up to 128 000 CPUs. The newly supported docking method will be valuable for drug discovery projects in which protein dynamics and flexibility play a significant role.

GWOVina: A grey wolf optimization approach to rigid and flexible receptor docking

Protein–ligand docking programs are indispensable tools for predicting the binding pose of a ligand to the receptor protein. In this paper, we introduce an efficient flexible docking method, gwovina, which is a variant of the Vina implementation using the grey wolf optimizer (GWO) and random walk for the global search, and the Dunbrack rotamer library for side‐chain sampling. The new method was validated for rigid and flexible‐receptor docking using four independent datasets. In rigid docking, gwovina showed comparable docking performance to Vina in terms of ligand pose RMSD, success rate, and affinity prediction. In flexible‐receptor docking, gwovina has improved success rate compared to Vina and AutoDockFR. It ran 2 to 7 times faster than Vina and 40 to 100 times faster than AutoDockFR. Therefore, gwovina can play a role in solving the complex flexible‐receptor docking cases and is suitable for virtual screening of compound libraries. gwovina is freely available at https://cbbio.cis.um.edu.mo/software/gwovina for testing.