Attracting cavities for docking. Replacing the rough energy landscape of the protein by a smooth attracting landscape.

Molecular docking is a computational approach for predicting the most probable position of ligands in the binding sites of macromolecules and constitutes the cornerstone of structure-based computer-aided drug design. Here, we present a new algorithm called Attracting Cavities that allows molecular docking to be performed by simple energy minimizations only. The approach consists in transiently replacing the rough potential energy hypersurface of the protein by a smooth attracting potential driving the ligands into protein cavities. The actual protein energy landscape is reintroduced in a second step to refine the ligand position. The scoring function of Attracting Cavities is based on the CHARMM force field and the FACTS solvation model. The approach was tested on the 85 experimental ligand-protein structures included in the Astex diverse set and achieved a success rate of 80% in reproducing the experimental binding mode starting from a completely randomized ligand conformer. The algorithm thus compares favorably with current state-of-the-art docking programs

Structure and Plasticity of Indoleamine 2,3-Dioxygenase 1 (IDO1).

Since the discovery of the implication of indoleamine 2,3-dioxygenase 1 (IDO1) in tumoral immune resistance in 2003, the search for inhibitors has been intensely pursued both in academia and in pharmaceutical companies, supported by the publication of the first crystal structure of IDO1 in 2006. More recently, it has become clear that IDO1 is an important player in various biological pathways and diseases ranging from neurodegenerative diseases to infection and autoimmunity. Its inhibition may lead to clinical benefit in different therapeutic settings. At present, over 50 experimental structures of IDO1 in complex with different ligands are available in the Protein Data Bank. Our analysis of this wealth of structural data sheds new light on several open issues regarding IDO1's structure and function

On-the-Fly QM/MM Docking with Attracting Cavities.

We developed a hybrid quantum mechanical/molecular mechanical (QM/MM) on-the-fly docking algorithm to address the challenges of treating polarization and selected metal interactions in docking. The algorithm is based on our classical docking algorithm Attracting Cavities and relies on the semiempirical self-consistent charge density functional tight-binding (SCC-DFTB) method and the CHARMM force field. We benchmarked the performance of this approach on three very diverse data sets: (1) the Astex Diverse set of 85 common noncovalent drug/target complexes formed both by hydrophobic and electrostatic interactions; (2) a zinc metalloprotein data set of 281 complexes, where polarization is strong and ligand/protein interactions are dominated by electrostatic interactions; and (3) a heme protein data set of 72 complexes, where ligand/protein interactions are dominated by covalent ligand/iron binding. Redocking performance of the on-the-fly QM/MM docking algorithm was compared to the performance of classical Attracting Cavities, AutoDock, AutoDock Vina, and GOLD. The results demonstrate that the QM/MM code preserves the high accuracy of most classical scores on the Astex Diverse set, while it yields significant improvements on both sets of metalloproteins at moderate computational cost

Toward on-the-fly quantum mechanical/molecular mechanical (QM/MM) docking: development and benchmark of a scoring function.

We address the challenges of treating polarization and covalent interactions in docking by developing a hybrid quantum mechanical/molecular mechanical (QM/MM) scoring function based on the semiempirical self-consistent charge density functional tight-binding (SCC-DFTB) method and the CHARMM force field. To benchmark this scoring function within the EADock DSS docking algorithm, we created a publicly available dataset of high-quality X-ray structures of zinc metalloproteins ( http://www.molecular-modelling.ch/resources.php ). For zinc-bound ligands (226 complexes), the QM/MM scoring yielded a substantially improved success rate compared to the classical scoring function (77.0% vs 61.5%), while, for allosteric ligands (55 complexes), the success rate remained constant (49.1%). The QM/MM scoring significantly improved the detection of correct zinc-binding geometries and improved the docking success rate by more than 20% for several important drug targets. The performance of both the classical and the QM/MM scoring functions compare favorably to the performance of AutoDock4, AutoDock4Zn, and AutoDock Vina

Two-Step Covalent Docking with Attracting Cavities.

Due to their various advantages, interest in the development of covalent drugs has been renewed in the past few years. It is therefore important to accurately describe and predict their interactions with biological targets by computer-aided drug design tools such as docking algorithms. Here, we report a covalent docking procedure for our in-house docking code Attracting Cavities (AC), which mimics the two-step mechanism of covalent ligand binding. Ligand binding to the protein cavity is driven by nonbonded interactions, followed by the formation of a covalent bond between the ligand and the protein through a chemical reaction. To test the performance of this method, we developed a diverse, high-quality, openly accessible re-docking benchmark set of 95 covalent complexes bound by 8 chemical reactions to 5 different reactive amino acids. Combination with structures from previous studies resulted in a set of 304 complexes, on which AC obtained a success rate (rmsd ≤ 2 Å) of 78%, outperforming two state-of-the-art covalent docking codes, genetic optimization for ligand docking (GOLD (66%)) and AutoDock (AD (35%)). Using a more stringent success criterion (rmsd ≤ 1.5 Å), AC reached a success rate of 71 vs 55% for GOLD and 26% for AD. We additionally assessed the cross-docking performance of AC on a set of 76 covalent complexes of the SARS-CoV-2 main protease. On this challenging test set of mainly small and highly solvent-exposed ligands, AC yielded success rates of 58 and 28% for re-docking and cross-docking, respectively, compared to 45 and 17% for GOLD

The Binding Mode of N-Hydroxyamidines to Indoleamine 2,3-Dioxygenase 1 (IDO1).

Indoleamine 2,3-dioxygenase 1 (IDO1) is an important target in cancer immunotherapy. The most advanced clinical compound, epacadostat (INCB024360), binds to the heme cofactor of IDO1 through an N-hydroxyamidine function. Conflicting binding modes have recently been proposed, reporting iron binding either through the hydroxyamidine oxygen or through the hydroxyamidine nitrogen atom. Here, we use quantum chemical calculations, docking, and quantum mechanics/molecular mechanics calculations based on available X-ray data to resolve this issue and to propose a physically meaningful binding mode. Our findings will aid the design of novel IDO1 ligands based on this pharmacophore

Attracting Cavities 2.0: Improving the Flexibility and Robustness for Small-Molecule Docking.

Molecular docking is a computational approach for predicting the most probable position of a ligand in the binding site of a target macromolecule. Our docking algorithm Attracting Cavities (AC) has been shown to compare favorably to other widely used docking algorithms [Zoete, V.; et al. J. Comput. Chem. 2016, 37, 437]. Here we describe several improvements of AC, making the sampling more robust and providing more flexibility for either fast or high-accuracy docking. We benchmark the performance of AC 2.0 using the 285 complexes of the PDBbind Core set, version 2016. For redocking from randomized ligand conformations, AC 2.0 reaches a success rate of 73.3%, compared to 63.9% for GOLD and 58.0% for AutoDock Vina. Due to its force-field-based scoring function and its thorough sampling procedure, AC 2.0 also performs well for blind docking on the entire receptor surface. The accuracy of its scoring function allows for the detection of problematic experimental structures in the benchmark set. For cross-docking, the AC 2.0 success rate is about 30% lower than for redocking (42.5%), similar to GOLD (42.8%) and better than AutoDock Vina (33.1%), and it can be improved by an informed choice of flexible protein residues. For selected targets with a high success rate in cross-docking, AC 2.0 also achieves good enrichment factors in virtual screening

Challenges in the Discovery of Indoleamine 2,3-Dioxygenase 1 (IDO1) Inhibitors.

Since the discovery of indoleamine 2,3-dioxygenase 1 (IDO1) as an attractive target for anticancer therapy in 2003, the search for inhibitors has been intensely pursued both in academia and in pharmaceutical companies. Many novel IDO1 inhibitor scaffolds have been described, and a few potent compounds have entered clinical trials. However, a significant number of the reported compounds contain problematic functional groups, suggesting that enzyme inhibition could be the result of undesirable side reactions instead of selective binding to IDO1. Here, we describe issues in the employed experimental protocols, review and classify reported IDO1 inhibitors, and suggest different approaches for confirming viable inhibitor scaffolds

Structural Prediction of Peptide-MHC Binding Modes.

The immune system is constantly protecting its host from the invasion of pathogens and the development of cancer cells. The specific CD8 <sup>+</sup> T-cell immune response against virus-infected cells and tumor cells is based on the T-cell receptor recognition of antigenic peptides bound to class I major histocompatibility complexes (MHC) at the surface of antigen presenting cells. Consequently, the peptide binding specificities of the highly polymorphic MHC have important implications for the design of vaccines, for the treatment of autoimmune diseases, and for personalized cancer immunotherapy. Evidence-based machine-learning approaches have been successfully used for the prediction of peptide binders and are currently being developed for the prediction of peptide immunogenicity. However, understanding and modeling the structural details of peptide/MHC binding is crucial for a better understanding of the molecular mechanisms triggering the immunological processes, estimating peptide/MHC affinity using universal physics-based approaches, and driving the design of novel peptide ligands. Unfortunately, due to the large diversity of MHC allotypes and possible peptides, the growing number of 3D structures of peptide/MHC (pMHC) complexes in the Protein Data Bank only covers a small fraction of the possibilities. Consequently, there is a growing need for rapid and efficient approaches to predict 3D structures of pMHC complexes. Here, we review the key characteristics of the 3D structure of pMHC complexes before listing databases and other sources of information on pMHC structures and MHC specificities. Finally, we discuss some of the most prominent pMHC docking software