A toolkit for covalent docking with GOLD: from automated ligand preparation with KNIME to bound protein-ligand complexes

Motivation

Current covalent docking tools have limitations that make them difficult to use for performing large-scale structure-based covalent virtual screening (VS). They require time-consuming tasks for the preparation of proteins and compounds (standardization, filtering according to the type of warheads), as well as for setting up covalent reactions. We have developed a toolkit to help accelerate drug discovery projects in the phases of hit identification by VS of ultra-large covalent libraries and hit expansion by exploration of the binding of known covalent compounds. With this application note, we offer the community a toolkit for performing automated covalent docking in a fast and efficient way.

Results

The toolkit comprises a KNIME workflow for ligand preparation and a Python program to perform the covalent docking of ligands with the GOLD docking engine running in a parallelized fashion.

Availability and implementation

The KNIME workflow entitled 'Evotec_Covalent_Processing_forGOLD.knwf' for the preparation of the ligands is available in the KNIME Hub https://hub.knime.com/emilie_pihan/spaces.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Fine tuning for success in structure-based virtual screening

Structure-based virtual screening plays a significant role in drug-discovery. The method virtually docks millions of compounds from corporate or public libraries into a binding site of a disease-related protein structure, allowing for the selection of a small list of potential ligands for experimental testing. Many algorithms are available for docking and assessing the affinity of compounds for a targeted protein site. The performance of affinity estimation calculations is highly dependent on the size and nature of the site, therefore a rationale for selecting the best protocol is required. To address this issue, we have developed an automated calibration process, implemented in a Knime workflow. It consists of four steps: preparation of a protein test set with structures and models of the target, preparation of a compound test set with target-related ligands and decoys, automatic test of 24 scoring/rescoring protocols for each target structure and model, and graphical display of results. The automation of the process combined with execution on high performance computing resources greatly reduces the duration of the calibration phase, and the test of many combinations of algorithms on various target conformations results in a rational and optimal choice of the best protocol. Here, we present this tool and exemplify its application in setting-up an optimal protocol for SBVS against Retinoid X Receptor alpha.

Modeling and Deorphanization of Orphan GPCRs

Despite tremendous efforts, approximately 120 GPCRs remain orphan. Their physiological functions and their potential roles in diseases are poorly understood. Orphan GPCRs are extremely important because they may provide novel therapeutic targets for unmet medical needs. As a complement to experimental approaches, molecular modeling and virtual screening are efficient techniques to discover synthetic surrogate ligands which can help to elucidate the role of oGPCRs. Constitutively activated mutants and recently published active structures of GPCRs provide stimulating opportunities for building active molecular models for oGPCRs and identifying activators using virtual screening of compound libraries. We describe the molecular modeling and virtual screening process we have applied in the discovery of surrogate ligands, and provide examples for CCKA, a simulated oGPCR, and for two oGPCRs, GPR52 and GPR34.

Computational and biophysical approaches to protein-protein interaction inhibition of Plasmodium falciparum AMA1/RON2 complex

Invasion of the red blood cell by Plasmodium falciparum parasites requires formation of an electron dense circumferential ring called the Moving Junction (MJ). The MJ is anchored by a high affinity complex of two parasite proteins: Apical Membrane Antigen 1 (PfAMA1) displayed on the surface of the parasite and Rhoptry Neck Protein 2 that is discharged from the parasite and imbedded in the membrane of the host cell. Structural studies of PfAMA1 revealed a conserved hydrophobic groove localized to the apical surface that coordinates RON2 and invasion inhibitory peptides. In the present work, we employed computational and biophysical methods to identify competitive P. falciparum AMA1-RON2 inhibitors with the goal of exploring the 'druggability' of this attractive antimalarial target. A virtual screen followed by molecular docking with the PfAMA1 crystal structure was performed using an eight million compound collection that included commercial molecules, the ChEMBL malaria library and approved drugs. The consensus approach resulted in the selection of inhibitor candidates. We also developed a fluorescence anisotropy assay using a modified inhibitory peptide to experimentally validate the ability of the selected compounds to inhibit the AMA1-RON2 interaction. Among those, we identified one compound that displayed significant inhibition. This study offers interesting clues to improve the throughput and reliability of screening for new drug leads.

Plasmodium falciparum CTP:phosphocholine cytidylyltransferase possesses two functional catalytic domains and is inhibited by a CDP-choline analog selected from a virtual screening

Phosphatidylcholine is the major lipid component of the malaria parasite membranes and is required for parasite multiplication in human erythrocytes. Plasmodium falciparum CTP:phosphocholine cytidylyltransferase (PfCCT) is the rate-limiting enzyme of the phosphatidylcholine biosynthesis pathway and thus considered as a potential antimalarial target. In contrast to its mammalian orthologs, PfCCT contains a duplicated catalytic domain. Here, we show that both domains are catalytically active with similar kinetic parameters. A virtual screening strategy allowed the identification of a drug-size molecule competitively inhibiting the enzyme. This compound also prevented phosphatidylcholine biosynthesis in parasites and exerted an antimalarial effect. This study constitutes the first step towards a rationalized design of future new antimalarial agents targeting PfCCT.

e-Drug3D: 3D structure collections dedicated to drug repurposing and fragment-based drug design

MOTIVATION

In the drug discovery field, new uses for old drugs, selective optimization of side activities and fragment-based drug design (FBDD) have proved to be successful alternatives to high-throughput screening. e-Drug3D is a database of 3D chemical structures of drugs that provides several collections of ready-to-screen SD files of drugs and commercial drug fragments. They are natural inputs in studies dedicated to drug repurposing and FBDD.

AVAILABILITY

e-Drug3D collections are freely available at http://chemoinfo.ipmc.cnrs.fr/e-drug3d.html either for download or for direct in silico web-based screenings.