SwissDock 2024: major enhancements for small-molecule docking with Attracting Cavities and AutoDock Vina.

Drug discovery aims to identify potential therapeutic compounds capable of modulating the activity of specific biological targets. Molecular docking can efficiently support this process by predicting binding interactions between small molecules and macromolecular targets and potentially accelerating screening campaigns. SwissDock is a computational tool released in 2011 as part of the SwissDrugDesign project, providing a free web-based service for small-molecule docking after automatized preparation of ligands and targets. Here, we present the latest version of SwissDock, in which EADock DSS has been replaced by two state-of-the-art docking programs, i.e. Attracting Cavities and AutoDock Vina. AutoDock Vina provides faster docking predictions, while Attracting Cavities offers more accurate results. Ligands can be imported in various ways, including as files, SMILES notation or molecular sketches. Targets can be imported as PDB files or identified by their PDB ID. In addition, advanced search options are available both for ligands and targets, giving users automatized access to widely-used databases. The web interface has been completely redesigned for interactive submission and analysis of docking results. Moreover, we developed a user-friendly command-line access which, in addition to all options of the web site, also enables covalent ligand docking with Attracting Cavities. The new version of SwissDock is freely available at https://www.swissdock.ch/.

The impact of structural bioinformatics tools and resources on SARS-CoV-2 research and therapeutic strategies.

SARS-CoV-2 is the causative agent of COVID-19, the ongoing global pandemic. It has posed a worldwide challenge to human health as no effective treatment is currently available to combat the disease. Its severity has led to unprecedented collaborative initiatives for therapeutic solutions against COVID-19. Studies resorting to structure-based drug design for COVID-19 are plethoric and show good promise. Structural biology provides key insights into 3D structures, critical residues/mutations in SARS-CoV-2 proteins, implicated in infectivity, molecular recognition and susceptibility to a broad range of host species. The detailed understanding of viral proteins and their complexes with host receptors and candidate epitope/lead compounds is the key to developing a structure-guided therapeutic design. Since the discovery of SARS-CoV-2, several structures of its proteins have been determined experimentally at an unprecedented speed and deposited in the Protein Data Bank. Further, specialized structural bioinformatics tools and resources have been developed for theoretical models, data on protein dynamics from computer simulations, impact of variants/mutations and molecular therapeutics. Here, we provide an overview of ongoing efforts on developing structural bioinformatics tools and resources for COVID-19 research. We also discuss the impact of these resources and structure-based studies, to understand various aspects of SARS-CoV-2 infection and therapeutic development. These include (i) understanding differences between SARS-CoV-2 and SARS-CoV, leading to increased infectivity of SARS-CoV-2, (ii) deciphering key residues in the SARS-CoV-2 involved in receptor-antibody recognition, (iii) analysis of variants in host proteins that affect host susceptibility to infection and (iv) analyses facilitating structure-based drug and vaccine design against SARS-CoV-2.

SwissParam 2023: A Modern Web-Based Tool for Efficient Small Molecule Parametrization.

Most steps of drug discovery are now routinely supported and accelerated by computer-aided drug design tools. Among them, structure-based approaches use the three-dimensional structure of the targeted biomacromolecule as a major source of information. When it comes to calculating the interactions of small molecules with proteins using the equations of molecular mechanics, topologies, atom typing, and force field parameters are required. However, generating parameters for small molecules remains challenging due to the large number of existing chemical groups. The SwissParam web tool was first released in 2011 with the aim of generating parameters and topologies for small molecules based on the Merck molecular force field (MMFF) while being compatible with the CHARMM22/27 force field. Here, we present an updated version of SwissParam, providing various new features, including the possibility to setup covalent ligands. Molecules can now be imported from different file formats or via a molecular sketcher. The MMFF-based approach has been updated to provide parameters and topologies compatible with the CHARMM36 force field. An option was added to generate small molecule parametrizations following the CHARMM General Force Field via the multipurpose atom-typer for CHARMM (MATCH) approach. Additionally, SwissParam now generates information on probable alternative tautomers and protonation states of the query molecule so that the user can consider all microspecies relevant to its compound. The new version of SwissParam is freely available at www.swissparam.ch and can also be accessed through a newly implemented command-line interface.

The SwissSimilarity 2021 Web Tool: Novel Chemical Libraries and Additional Methods for an Enhanced Ligand-Based Virtual Screening Experience.

Hit finding, scaffold hopping, and structure-activity relationship studies are important tasks in rational drug discovery. Implementation of these tasks strongly depends on the availability of compounds similar to a known bioactive molecule. SwissSimilarity is a web tool for low-to-high-throughput virtual screening of multiple chemical libraries to find molecules similar to a compound of interest. According to the similarity principle, the output list of molecules generated by SwissSimilarity is expected to be enriched in compounds that are likely to share common protein targets with the query molecule and that can, therefore, be acquired and tested experimentally in priority. Compound libraries available for screening using SwissSimilarity include approved drugs, clinical candidates, known bioactive molecules, commercially available and synthetically accessible compounds. The first version of SwissSimilarity launched in 2015 made use of various 2D and 3D molecular descriptors, including path-based FP2 fingerprints and ElectroShape vectors. However, during the last few years, new fingerprinting methods for molecular description have been developed or have become popular. Here we would like to announce the launch of the new version of the SwissSimilarity web tool, which features additional 2D and 3D methods for estimation of molecular similarity: extended-connectivity, MinHash, 2D pharmacophore, extended reduced graph, and extended 3D fingerprints. Moreover, it is now possible to screen for molecular structures having the same scaffold as the query compound. Additionally, all compound libraries available for screening in SwissSimilarity have been updated, and several new ones have been added to the list. Finally, the interface of the website has been comprehensively rebuilt to provide a better user experience. The new version of SwissSimilarity is freely available starting from December 2021.

SwissTargetPrediction: updated data and new features for efficient prediction of protein targets of small molecules.

SwissTargetPrediction is a web tool, on-line since 2014, that aims to predict the most probable protein targets of small molecules. Predictions are based on the similarity principle, through reverse screening. Here, we describe the 2019 version, which represents a major update in terms of underlying data, backend and web interface. The bioactivity data were updated, the model retrained and similarity thresholds redefined. In the new version, the predictions are performed by searching for similar molecules, in 2D and 3D, within a larger collection of 376 342 compounds known to be experimentally active on an extended set of 3068 macromolecular targets. An efficient backend implementation allows to speed up the process that returns results for a druglike molecule on human proteins in 15-20 s. The refreshed web interface enhances user experience with new features for easy input and improved analysis. Interoperability capacity enables straightforward submission of any input or output molecule to other on-line computer-aided drug design tools, developed by the SIB Swiss Institute of Bioinformatics. High levels of predictive performance were maintained despite more extended biological and chemical spaces to be explored, e.g. achieving at least one correct human target in the top 15 predictions for >70% of external compounds. The new SwissTargetPrediction is available free of charge (www.swisstargetprediction.ch).

Application of the SwissDrugDesign Online Resources in Virtual Screening.

SwissDrugDesign is an important initiative led by the Molecular Modeling Group of the SIB Swiss Institute of Bioinformatics. This project provides a collection of freely available online tools for computer-aided drug design. Some of these web-based methods, i.e., SwissSimilarity and SwissTargetPrediction, were especially developed to perform virtual screening, while others such as SwissADME, SwissDock, SwissParam and SwissBioisostere can find applications in related activities. The present review aims at providing a short description of these methods together with examples of their application in virtual screening, where SwissDrugDesign tools successfully supported the discovery of bioactive small molecules.

Educational Tools to Introduce Computer-Aided Drug Design to Students and to the Public at Large.

The Drug Design Workshop initiative was put in place in 2015 and consists of a collection of educational tools especially developed to introduce computer-aided drug design to the general public and students of various levels. These presentations, hands-on sessions, physical material and on-line educational tools (http://www.drug-design-workshop.ch) have been used in a variety of settings including classrooms, universities, teacher training sessions and science fairs. The main goal is to inform an audience as broad as possible regarding the processes and challenges for the design, discovery and development of drugs. The present article describes what is presently available and the future direction for new activities currently under development.

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.

SwissSimilarity: A Web Tool for Low to Ultra High Throughput Ligand-Based Virtual Screening.

SwissSimilarity is a new web tool for rapid ligand-based virtual screening of small to unprecedented ultralarge libraries of small molecules. Screenable compounds include drugs, bioactive and commercial molecules, as well as 205 million of virtual compounds readily synthesizable from commercially available synthetic reagents. Predictions can be carried out on-the-fly using six different screening approaches, including 2D molecular fingerprints as well as superpositional and fast nonsuperpositional 3D similarity methodologies. SwissSimilarity is part of a large initiative of the SIB Swiss Institute of Bioinformatics to provide online tools for computer-aided drug design, such as SwissDock, SwissBioisostere or SwissTargetPrediction with which it can interoperate, and is linked to other well-established online tools and databases. User interface and backend have been designed for simplicity and ease of use, to provide proficient virtual screening capabilities to specialists and nonexperts in the field. SwissSimilarity is accessible free of charge or login at http://www.swisssimilarity.ch .

SwissTargetPrediction: a web server for target prediction of bioactive small molecules.

Bioactive small molecules, such as drugs or metabolites, bind to proteins or other macro-molecular targets to modulate their activity, which in turn results in the observed phenotypic effects. For this reason, mapping the targets of bioactive small molecules is a key step toward unraveling the molecular mechanisms underlying their bioactivity and predicting potential side effects or cross-reactivity. Recently, large datasets of protein-small molecule interactions have become available, providing a unique source of information for the development of knowledge-based approaches to computationally identify new targets for uncharacterized molecules or secondary targets for known molecules. Here, we introduce SwissTargetPrediction, a web server to accurately predict the targets of bioactive molecules based on a combination of 2D and 3D similarity measures with known ligands. Predictions can be carried out in five different organisms, and mapping predictions by homology within and between different species is enabled for close paralogs and orthologs. SwissTargetPrediction is accessible free of charge and without login requirement at http://www.swisstargetprediction.ch.

iLOGP: a simple, robust, and efficient description of n-octanol/water partition coefficient for drug design using the GB/SA approach.

The n-octanol/water partition coefficient (log Po/w) is a key physicochemical parameter for drug discovery, design, and development. Here, we present a physics-based approach that shows a strong linear correlation between the computed solvation free energy in implicit solvents and the experimental log Po/w on a cleansed data set of more than 17,500 molecules. After internal validation by five-fold cross-validation and data randomization, the predictive power of the most interesting multiple linear model, based on two GB/SA parameters solely, was tested on two different external sets of molecules. On the Martel druglike test set, the predictive power of the best model (N = 706, r = 0.64, MAE = 1.18, and RMSE = 1.40) is similar to six well-established empirical methods. On the 17-drug test set, our model outperformed all compared empirical methodologies (N = 17, r = 0.94, MAE = 0.38, and RMSE = 0.52). The physical basis of our original GB/SA approach together with its predictive capacity, computational efficiency (1 to 2 s per molecule), and tridimensional molecular graphics capability lay the foundations for a promising predictor, the implicit log P method (iLOGP), to complement the portfolio of drug design tools developed and provided by the SIB Swiss Institute of Bioinformatics.

Testing the predictive power of reverse screening to infer drug targets, with the help of machine learning.

Estimating protein targets of compounds based on the similarity principle-similar molecules are likely to show comparable bioactivity-is a long-standing strategy in drug research. Having previously quantified this principle, we present here a large-scale evaluation of its predictive power for inferring macromolecular targets by reverse screening an unprecedented vast external test set of more than 300,000 active small molecules against another bioactivity set of more than 500,000 compounds. We show that machine-learning can predict the correct targets, with the highest probability among 2069 proteins, for more than 51% of the external molecules. The strong enrichment thus obtained demonstrates its usefulness in supporting phenotypic screens, polypharmacology, or repurposing. Moreover, we quantified the impact of the bioactivity knowledge available for proteins in terms of number and diversity of actives. Finally, we advise that developers of such approaches follow an application-oriented benchmarking strategy and use large, high-quality, non-overlapping datasets as provided here.

Molecular docking using the molecular lipophilicity potential as hydrophobic descriptor: impact on GOLD docking performance.

GOLD is a molecular docking software widely used in drug design. In the initial steps of docking, it creates a list of hydrophobic fitting points inside protein cavities that steer the positioning of ligand hydrophobic moieties. These points are generated based on the Lennard-Jones potential between a carbon probe and each atom of the residues delimitating the binding site. To thoroughly describe hydrophobic regions in protein pockets and properly guide ligand hydrophobic moieties toward favorable areas, an in-house tool, the MLP filter, was developed and herein applied. This strategy only retains GOLD hydrophobic fitting points that match the rigorous definition of hydrophobicity given by the molecular lipophilicity potential (MLP), a molecular interaction field that relies on an atomic fragmental system based on 1-octanol/water experimental partition coefficients (log P(oct)). MLP computations in the binding sites of crystallographic protein structures revealed that a significant number of points considered hydrophobic by GOLD were actually polar according to the MLP definition of hydrophobicity. To examine the impact of this new tool, ligand-protein complexes from the Astex Diverse Set and the PDB bind core database were redocked with and without the use of the MLP filter. Reliable docking results were obtained by using the MLP filter that increased the quality of docking in nonpolar cavities and outperformed the standard GOLD docking approach.

Secondary IDH1 resistance mutations and oncogenic IDH2 mutations cause acquired resistance to ivosidenib in cholangiocarcinoma.

The mutant IDH1 inhibitor ivosidenib improves outcomes for patients with IDH1-mutated cholangiocarcinoma, but resistance inevitably develops. Mechanisms of resistance and strategies to overcome resistance are poorly understood. Here we describe two patients with IDH1 R132C-mutated metastatic cholangiocarcinoma who developed acquired resistance to ivosidenib. After disease progression, one patient developed an oncogenic IDH2 mutation, and the second patient acquired a secondary IDH1 D279N mutation. To characterize the putative IDH1 resistance mutation, cells expressing the double-mutant were generated. In vitro, IDH1 R132H/D279N produces (R)-2HG less efficiently than IDH1 R132H. However, its binding to ivosidenib is impaired and it retains the ability to produce (R)-2HG and promote cellular transformation in the presence of ivosidenib. The irreversible mutant IDH1 inhibitor LY3410738 binds and blocks (R)-2HG production and cellular transformation by IDH1 R132H/D279N. These resistance mechanisms suggest that IDH1-mutated cholangiocarcinomas remain dependent on (R)-2HG even after prolonged ivosidenib treatment. Sequential mutant IDH inhibitor therapy should be explored as a strategy to overcome acquired resistance to mutant IDH inhibitors.

FGFR2 Extracellular Domain In-Frame Deletions Are Therapeutically Targetable Genomic Alterations That Function as Oncogenic Drivers in Cholangiocarcinoma.

We conducted next-generation DNA sequencing on 335 biliary tract cancers and characterized the genomic landscape by anatomic site within the biliary tree. In addition to frequent FGFR2 fusions among patients with intrahepatic cholangiocarcinoma (IHCC), we identified FGFR2 extracellular domain in-frame deletions (EID) in 5 of 178 (2.8%) patients with IHCC, including two patients with FGFR2 p.H167_N173del. Expression of this FGFR2 EID in NIH3T3 cells resulted in constitutive FGFR2 activation, oncogenic transformation, and sensitivity to FGFR inhibitors. Three patients with FGFR2 EIDs were treated with Debio 1347, an oral FGFR1/2/3 inhibitor, and all showed partial responses. One patient developed an acquired L618F FGFR2 kinase domain mutation at disease progression and experienced a further partial response for 17 months to an irreversible FGFR2 inhibitor, futibatinib. Together, these findings reveal FGFR2 EIDs as an alternative mechanism of FGFR2 activation in IHCC that predicts sensitivity to FGFR inhibitors in the clinic. SIGNIFICANCE: FGFR2 EIDs are transforming genomic alterations that occur predominantly in patients with IHCC. These FGFR2 EIDs are sensitive to FGFR inhibition in vitro, and patients with these alterations benefited from treatment with FGFR inhibitors in the clinic.This article is highlighted in the In This Issue feature, p. 2355.

Molecular electrocatalysis for oxygen reduction by cobalt porphyrins adsorbed at liquid/liquid interfaces.

Molecular electrocatalysis for oxygen reduction at a polarized water/1,2-dichloroethane (DCE) interface was studied, involving aqueous protons, ferrocene (Fc) in DCE and amphiphilic cobalt porphyrin catalysts adsorbed at the interface. The catalyst, (2,8,13,17-tetraethyl-3,7,12,18-tetramethyl-5-p-amino-phenylporphyrin) cobalt(II) (CoAP), functions like conventional cobalt porphyrins, activating O(2) via coordination by the formation of a superoxide structure. Furthermore, due to the hydrophilic nature of the aminophenyl group, CoAP has a strong affinity for the water/DCE interface as evidenced by lipophilicity mapping calculations and surface tension measurements, facilitating the protonation of the CoAP-O(2) complex and its reduction by ferrocene. The reaction is electrocatalytic as its rate depends on the applied Galvani potential difference between the two phases.

Stereoselective block of hERG channel by bupivacaine scrutinized at molecular level.

In the heart, the hERG voltage-gated potassium channel mediates the I(Kr) current, which is crucial for the duration of cardiac action potential. Undesired block of the channel may prolong the QT interval with increased risk of malignant ventricular arrhythmia called torsades de pointes. Although the molecular determinants of hERG block are intensively studied, stereoselectivity has been poorly investigated. Levo-(S)-bupivacaine was the first drug reported to have higher affinity for hERG than its enantiomer. This study aims at understanding the principles underlying the stereoselectivity of bupivacaine block with the help of molecular modeling. Putative binding modes of levo-(S)- and dextro-(R)-bupivacaine inside an open form model of hERG channel were predicted by docking simulations, allowing a clear depiction of ligand-protein interactions. Estimated binding energies for both enantiomers to wild-type channel are in line with previously published electrophysiology measurements. These results may be considered as a confirmation at the molecular level of bupivacaine stereoselective binding towards hERG. Moreover this information lays the foundations for a structural guideline to filter out potentially cardiotoxic drug candidates in silico.

Rational Design, Synthesis, and Pharmacological Characterization of Novel Ghrelin Receptor Inverse Agonists as Potential Treatment against Obesity-Related Metabolic Diseases.

A new chemotype of ghrelin inverse agonists was discovered through chimeric design based on molecular scaffolds known as growth-hormone secretagogue receptor (GHSR) modulators but with divergent pharmacodynamic and pharmacokinetic properties. The structure-activities/properties exploration led to compound 47, which displayed potent human GHSR antagonism and inverse agonism in cellular assays (IC50 = 68 nM, EC50 = 29 nM), moderate oral bioavailability, and notable brain penetration in rat ( F = 27%, B/ P ratio = 1.9). First in vivo studies demonstrated effective reduction of food intake after oral or parenteral administration to mouse (78% at 1 h and 38% at 8 h, respectively). Further preclinical studies are needed to evaluate the most suited mode of administration with the aim of promoting a first central-acting ghrelin inverse agonist molecule to development, which would represent a significant step toward therapeutic agents to treat metabolic disorders related to obesity, such as type 2 diabetes mellitus.

SwissADME: a free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules.

To be effective as a drug, a potent molecule must reach its target in the body in sufficient concentration, and stay there in a bioactive form long enough for the expected biologic events to occur. Drug development involves assessment of absorption, distribution, metabolism and excretion (ADME) increasingly earlier in the discovery process, at a stage when considered compounds are numerous but access to the physical samples is limited. In that context, computer models constitute valid alternatives to experiments. Here, we present the new SwissADME web tool that gives free access to a pool of fast yet robust predictive models for physicochemical properties, pharmacokinetics, drug-likeness and medicinal chemistry friendliness, among which in-house proficient methods such as the BOILED-Egg, iLOGP and Bioavailability Radar. Easy efficient input and interpretation are ensured thanks to a user-friendly interface through the login-free website http://www.swissadme.ch. Specialists, but also nonexpert in cheminformatics or computational chemistry can predict rapidly key parameters for a collection of molecules to support their drug discovery endeavours.

Impact of species-dependent differences on screening, design, and development of MAO B inhibitors.

The impact of species-dependent differences between human and rat MAO B on inhibitor screening was evidenced for two classes of compounds, coumarin and 5H-indeno[1,2-c]pyridazin-5-one derivatives. All examined compounds have shown a greater inhibitor potency toward human MAO B than toward rat MAO B. Moreover, no correlation was found between human and rat pIC(50) values. These divergences have important implications for the design and development of drugs involved in the MAO B metabolic pathway, suggesting that results obtained using rat enzyme cannot be extrapolated to human CNS, a priori. Indeed, the selection of a hit compound for lead generation could be different using human rather than rat enzyme. Moreover, the influence of substituents on the in vitro inhibition of human MAO B was markedly different between homogeneous series of coumarin and 5H-indeno[1,2-c]pyridazin-5-one derivatives, suggesting different binding modes, a hypothesis clearly supported by molecular docking simulations of inhibitors into the active site of human MAO B.