Search

Open-source workflows have become more and more an integral part of computer-aided drug design (CADD) projects since they allow reproducible and shareable research that can be easily transferred to other projects. Setting up, understanding, and applying such workflows involves either coding or using workflow managers that offer a graphical user interface. We previously reported the TeachOpenCADD teaching platform that provides interactive Jupyter Notebooks (talktorials) on central CADD topics using open-source data and Python packages. Here we present the conversion of these talktorials to KNIME workflows that allow users to explore our teaching material without any line of code. TeachOpenCADD KNIME workflows are freely available on the KNIME Hub: https://hub.knime.com/volkamerlab/space/TeachOpenCADD .

Alchemical free-energy calculations are now widely used to drive or maintain potency in small-molecule lead optimization with a roughly 1 kcal/mol accuracy. Despite this, the potential to use free-energy calculations to drive optimization of compound selectivity among two similar targets has been relatively unexplored in published studies. In the most optimistic scenario, the similarity of binding sites might lead to a fortuitous cancellation of errors and allow selectivity to be predicted more accurately than affinity. Here, we assess the accuracy with which selectivity can be predicted in the context of small-molecule kinase inhibitors, considering the very similar binding sites of human kinases CDK2 and CDK9 as well as another series of ligands attempting to achieve selectivity between the more distantly related kinases CDK2 and ERK2. Using a Bayesian analysis approach, we separate systematic from statistical errors and quantify the correlation in systematic errors between selectivity targets. We find that, in the CDK2/CDK9 case, a high correlation in systematic errors suggests that free-energy calculations can have significant impact in aiding chemists in achieving selectivity, while in more distantly related kinases (CDK2/ERK2), the correlation in systematic error suggests that fortuitous cancellation may even occur between systems that are not as closely related. In both cases, the correlation in systematic error suggests that longer simulations are beneficial to properly balance statistical error with systematic error to take full advantage of the increase in apparent free-energy calculation accuracy in selectivity prediction.

Many doctoral students and postdoctoral fellows face at some point in their career the decision between continuing in academia or pursuing a job in industry. Both career paths come with their advantages and disadvantages as well as associated clichés. Our scientific journeys have led us from an university Ph.D. degree to an industrial postdoctoral stay and back to a young faculty position in academia. In this perspective, we share our experiences while changing perspectives. We will discuss the insights we gained through the phase as industrial postdoctoral fellows, the motivation to return and take up a young faculty position in academia, and the freedom and the burden of starting out as a principal investigator (PI). We end with our thoughts on "quo vadis" computational chemistry.

The identification and design of selective compounds is important for the reduction of unwanted side effects as well as for the development of tool compounds for target validation studies. This is, in particular, true for therapeutically important protein families that possess conserved folds and have numerous members such as kinases. To support the design of selective kinase inhibitors, we developed a novel approach that allows identification of specificity determining subpockets between closely related kinases solely based on their three-dimensional structures. To account for the intrinsic flexibility of the proteins, multiple X-ray structures of the target protein of interest as well as of unwanted off-target(s) are taken into account. The binding pockets of these protein structures are calculated and fused to a combined target and off-target pocket, respectively. Subsequently, shape differences between these two combined pockets are identified via fusion rules. The approach provides a user-friendly visualization of target-specific areas in a binding pocket which should be explored when designing selective compounds. Furthermore, the approach can be easily combined with in silico alanine mutation studies to identify selectivity determining residues. The potential impact of the approach is demonstrated in four retrospective experiments on closely related kinases, i.e., p38α vs Erk2, PAK1 vs PAK4, ITK vs AurA, and BRAF vs VEGFR2. Overall, the presented approach does not require any profiling data for training purposes, provides an intuitive visualization of a large number of protein structures at once, and could also be applied to other target classes.

Protein kinases are involved in a variety of diseases including cancer, inflammation, and autoimmune disorders. Although the development of new kinase inhibitors is a major focus in pharmaceutical research, a large number of kinases remained so far unexplored in drug discovery projects. The selection and assessment of targets is an essential but challenging area. Today, a few thousands of experimentally determined kinase structures are available, covering about half of the human kinome. This large structural source allows guiding the target selection via structure-based druggability prediction approaches such as DoGSiteScorer. Here, a thorough analysis of the ATP pockets of the entire human kinome in the DFG-in state is presented in order to prioritize novel kinase structures for drug discovery projects. For this, all human kinase X-ray structures available in the PDB were collected, and homology models were generated for the missing part of the kinome. DoGSiteScorer was used to calculate geometrical and physicochemical properties of the ATP pockets and to predict the potential of each kinase to be druggable. The results indicate that about 75% of the kinome are in principle druggable. Top ranking structures comprise kinases that are primary targets of known approved drugs but additionally point to so far less explored kinases. The presented analysis provides new insights into the druggability of ATP binding pockets of the entire kinome. We anticipate this comprehensive druggability assessment of protein kinases to be helpful for the community to prioritize so far untapped kinases for drug discovery efforts.

A common task in drug development is the selection of compounds fulfilling specific structural features from a large data pool. While several methods that iteratively search through such data sets exist, their application is limited compared to the infinite character of molecular space. The introduction of the concept of fragment spaces (FSs), which are composed of molecular fragments and their connection rules, made the representation of large combinatorial data sets feasible. At the same time, search algorithms face the problem of structural features spanning over multiple fragments. Due to the combinatorial nature of FSs, an enumeration of all products is impossible. In order to overcome these time and storage issues, we present a method that is able to find substructures in FSs without explicit product enumeration. This is accomplished by splitting substructures into subsubstructures and mapping them onto fragments with respect to fragment connectivity rules. The method has been evaluated on three different drug discovery scenarios considering the exploration of a molecule class, the elaboration of decoration patterns for a molecular core, and the exhaustive query for peptides in FSs. FSs can be searched in seconds, and found products contain novel compounds not present in the PubChem database which may serve as hints for new lead structures.

Automated prediction of protein active sites is essential for large-scale protein function prediction, classification, and druggability estimates. In this work, we present DoGSite, a new structure-based method to predict active sites in proteins based on a Difference of Gaussian (DoG) approach which originates from image processing. In contrast to existing methods, DoGSite splits predicted pockets into subpockets, revealing a refined description of the topology of active sites. DoGSite correctly predicts binding pockets for over 92% of the PDBBind and the scPDB data set, being in line with the best-performing methods available. In 63% of the PDBBind data set the detected pockets can be subdivided into smaller subpockets. The cocrystallized ligand is contained in exactly one subpocket in 87% of the predictions. Furthermore, we introduce a more precise prediction performance measure by taking the pairwise ligand and pocket coverage into account. In 90% of the cases DoGSite predicts a pocket that contains at least half of the ligand. In 70% of the cases additionally more than a quarter of the respective pocket itself is covered by the cocrystallized ligand. Consideration of subpockets produces an increase in coverage yielding a success rate of 83% for the latter measure.

We present TrixP, a new index-based method for fast protein binding site comparison and function prediction. TrixP determines binding site similarities based on the comparison of descriptors that encode pharmacophoric and spatial features. Therefore, it adopts the efficient core components of TrixX, a structure-based virtual screening technology for large compound libraries. TrixP expands this technology by new components in order to allow a screening of protein libraries. TrixP accounts for the inherent flexibility of proteins employing a partial shape matching routine. After the identification of structures with matching pharmacophoric features and geometric shape, TrixP superimposes the binding sites and, finally, assesses their similarity according to the fit of pharmacophoric properties. TrixP is able to find analogies between closely and distantly related binding sites. Recovery rates of 81.8% for similar binding site pairs, assisted by rejecting rates of 99.5% for dissimilar pairs on a test data set containing 1331 pairs, confirm this ability. TrixP exclusively identifies members of the same protein family on top ranking positions out of a library consisting of 9802 binding sites. Furthermore, 30 predicted kinase binding sites can almost perfectly be classified into their known subfamilies.

Predicting druggability and prioritizing certain disease modifying targets for the drug development process is of high practical relevance in pharmaceutical research. DoGSiteScorer is a fully automatic algorithm for pocket and druggability prediction. Besides consideration of global properties of the pocket, also local similarities shared between pockets are reflected. Druggability scores are predicted by means of a support vector machine (SVM), trained, and tested on the druggability data set (DD) and its nonredundant version (NRDD). The DD consists of 1069 targets with assigned druggable, difficult, and undruggable classes. In 90% of the NRDD, the SVM model based on global descriptors correctly classifies a target as either druggable or undruggable. Nevertheless, global properties suffer from binding site changes due to ligand binding and from the pocket boundary definition. Therefore, local pocket properties are additionally investigated in terms of a nearest neighbor search. Local similarities are described by distance dependent histograms between atom pairs. In 88% of the DD pocket set, the nearest neighbor and the structure itself conform with their druggability type. A discriminant feature between druggable and undruggable pockets is having less short-range hydrophilic-hydrophilic pairs and more short-range lipophilic-lipophilic pairs. Our findings for global pocket descriptors coincide with previously published methods affirming that size, shape, and hydrophobicity are important global pocket descriptors for automatic druggability prediction. Nevertheless, the variety of pocket shapes and their flexibility upon ligand binding limit the automatic projection of druggable features onto descriptors. Incorporating local pocket properties is another step toward a reliable descriptor-based druggability prediction.