Comprehensive structural and functional characterization of the human kinome by protein structure modeling and ligand virtual screening

FRAGSITE2: A structure and fragment-based approach for virtual ligand screening

Protein function annotation and drug discovery often involve finding small molecule binders. In the early stages of drug discovery, virtual ligand screening (VLS) is frequently applied to identify possible hits before experimental testing. While our recent ligand homology modeling (LHM)-machine learning VLS method FRAGSITE outperformed approaches that combined traditional docking to generate protein-ligand poses and deep learning scoring functions to rank ligands, a more robust approach that could identify a more diverse set of binding ligands is needed. Here, we describe FRAGSITE2 that shows significant improvement on protein targets lacking known small molecule binders and no confident LHM identified template ligands when benchmarked on two commonly used VLS datasets: For both the DUD-E set and DEKOIS2.0 set and ligands having a Tanimoto coefficient (TC) < 0.7 to the template ligands, the 1% enrichment factor (EF1% ) of FRAGSITE2 is significantly better than those for FINDSITEcomb2.0 , an earlier LHM algorithm. For the DUD-E set, FRAGSITE2 also shows better ROC enrichment factor and AUPR (area under the precision-recall curve) than the deep learning DenseFS scoring function. Comparison with the RF-score-VS on the 76 target subset of DEKOIS2.0 and a TC < 0.99 to training DUD-E ligands, FRAGSITE2 has double the EF1% . Its boosted tree regression method provides for more robust performance than a deep learning multiple layer perceptron method. When compared with the pretrained language model for protein target features, FRAGSITE2 also shows much better performance. Thus, FRAGSITE2 is a promising approach that can discover novel hits for protein targets. FRAGSITE2's web service is freely available to academic users at http://sites.gatech.edu/cssb/FRAGSITE2.

Structure prediction of eukaryotic elongation factor-2 kinase and identification of the binding mechanisms of its inhibitors: homology modeling, molecular docking, and molecular dynamics simulation

Protein kinases emerged as one of the most successful families of drug targets due to their increased activity and involvement in mediating critical signal transduction pathways in cancer cells. Recent evidence suggests that eukaryotic elongation factor 2 kinase (eEF-2K) is a potential therapeutic target for treating some highly aggressive solid cancers, including lung, pancreatic and triple-negative breast cancers. Thus, several compounds have been developed for the inhibition of the enzyme activity, but they are not sufficiently specific and potent. Besides, the crystal structure of this kinase remains unknown. Hence, the functional organization and regulation of eEF-2K remain poorly characterized. For this purpose, we constructed a homology model of eEF-2K and then used docking methodology to better understanding the binding mechanism of eEF-2K with 58 compounds that have been proposed as existing inhibitors. The results of this analysis were compared with the experimental results and the compounds effective against eEF-2K were determined against eEF-2K as a result of both studies. And finally, molecular dynamics (MD) simulations were performed for the stability of eEF-2K with these compounds. According to these study defined that the binding mechanism of eEF-2K with inhibitors at the molecular level and elucidated the residues of eEF-2K that play an important role in enzyme selectivity and ligand affinity. This information may lead to new selective and potential drug molecules to be for inhibition of eEF-2K.Communicated by Ramaswamy H. Sarma.

Artificial intelligence to guide precision anticancer therapy with multitargeted kinase inhibitors

BACKGROUND

Vast amounts of rapidly accumulating biological data related to cancer and a remarkable progress in the field of artificial intelligence (AI) have paved the way for precision oncology. Our recent contribution to this area of research is CancerOmicsNet, an AI-based system to predict the therapeutic effects of multitargeted kinase inhibitors across various cancers. This approach was previously demonstrated to outperform other deep learning methods, graph kernel models, molecular docking, and drug binding pocket matching.

METHODS

CancerOmicsNet integrates multiple heterogeneous data by utilizing a deep graph learning model with sophisticated attention propagation mechanisms to extract highly predictive features from cancer-specific networks. The AI-based system was devised to provide more accurate and robust predictions than data-driven therapeutic discovery using gene signature reversion.

RESULTS

Selected CancerOmicsNet predictions obtained for "unseen" data are positively validated against the biomedical literature and by live-cell time course inhibition assays performed against breast, pancreatic, and prostate cancer cell lines. Encouragingly, six molecules exhibited dose-dependent antiproliferative activities, with pan-CDK inhibitor JNJ-7706621 and Src inhibitor PP1 being the most potent against the pancreatic cancer cell line Panc 04.03.

CONCLUSIONS

CancerOmicsNet is a promising AI-based platform to help guide the development of new approaches in precision oncology involving a variety of tumor types and therapeutics.

Structure-based virtual screening and biological evaluation of novel small-molecule BTK inhibitors

Bruton tyrosine kinase (BTK) is linked to multiple signalling pathways that regulate cellular survival, activation, and proliferation. A covalent BTK inhibitor has shown favourable outcomes for treating B cell malignant leukaemia. However, covalent inhibitors require a high reactive warhead that may contribute to unexpected toxicity, poor selectivity, or reduced effectiveness in solid tumours. Herein, we report the identification of a novel noncovalent BTK inhibitor. The binding interactions (i.e. interactions from known BTK inhibitors) for the BTK binding site were identified and incorporated into a structure-based virtual screening (SBVS). Top-rank compounds were selected and testing revealed a BTK inhibitor with >50% inhibition at 10 µM concentration. Examining analogues revealed further BTK inhibitors. When tested across solid tumour cell lines, one inhibitor showed favourable inhibitory activity, suggesting its potential for targeting BTK malignant tumours. This inhibitor could serve as a basis for developing an effective BTK inhibitor targeting solid cancers.

GraphGR: A graph neural network to predict the effect of pharmacotherapy on the cancer cell growth.. [DOI: 10.1101/2020.05.20.107458]

Genomic profiles of cancer cells provide valuable information on genetic alterations in cancer. Several recent studies employed these data to predict the response of cancer cell lines to treatment with drugs. Nonetheless, due to the multifactorial phenotypes and intricate mechanisms of cancer, the accurate prediction of the effect of pharmacotherapy on a specific cell line based on the genetic information alone is problematic. High prediction accuracies reported in the literature likely result from significant overlaps among training, validation, and testing sets, making many predictors inapplicable to new data. To address these issues, we developed GraphGR, a graph neural network with sophisticated attention propagation mechanisms to predict the therapeutic effects of kinase inhibitors across various tumors. Emphasizing on the system-level complexity of cancer, GraphGR integrates multiple heterogeneous data, such as biological networks, genomics, inhibitor profiling, and genedisease associations, into a unified graph structure. In order to construct diverse and information-rich cancer-specific networks, we devised a novel graph reduction protocol based on not only the topological information, but also the biological knowledge. The performance of GraphGR, properly cross-validated at the tissue level, is 0.83 in terms of the area under the receiver operating characteristics, which is notably higher than those measured for other approaches on the same data. Finally, several new predictions are validated against the biomedical literature demonstrating that GraphGR generalizes well to unseen data, i.e. it can predict therapeutic effects across a variety of cancer cell lines and inhibitors. GraphGR is freely available to the academic community at https://github.com/pulimeng/GraphGR.

Kinase Atlas: Druggability Analysis of Potential Allosteric Sites in Kinases

The inhibition of kinases has been pursued by the pharmaceutical industry for over 20 years. While the locations of the sites that bind type II and III inhibitors at or near the adenosine 5'-triphosphate binding sites are well defined, the literature describes 10 different regions that were reported as regulatory hot spots in some kinases and thus are potential target sites for type IV inhibitors. Kinase Atlas is a systematic collection of binding hot spots located at the above ten sites in 4910 structures of 376 distinct kinases available in the Protein Data Bank. The hot spots are identified by FTMap, a computational analogue of experimental fragment screening. Users of Kinase Atlas ( https://kinase-atlas.bu.edu ) may view summarized results for all structures of a particular kinase, such as which binding sites are present and how druggable they are, or they may view hot spot information for a particular kinase structure of interest.

FINDSITEcomb2.0: A New Approach for Virtual Ligand Screening of Proteins and Virtual Target Screening of Biomolecules

Computational approaches for predicting protein-ligand interactions can facilitate drug lead discovery and drug target determination. We have previously developed a threading/structural-based approach, FINDSITE^comb, for the virtual ligand screening of proteins that has been extensively experimentally validated. Even when low resolution predicted protein structures are employed, FINDSITE^comb has the advantage of being faster and more accurate than traditional high-resolution structure-based docking methods. It also overcomes the limitations of traditional QSAR methods that require a known set of seed ligands that bind to the given protein target. Here, we further improve FINDSITE^comb by enhancing its template ligand selection from the PDB/DrugBank/ChEMBL libraries of known protein-ligand interactions by (1) parsing the template proteins and their corresponding binding ligands in the DrugBank and ChEMBL libraries into domains so that the ligands with falsely matched domains to the targets will not be selected as template ligands; (2) applying various thresholds to filter out falsely matched template structures in the structure comparison process and thus their corresponding ligands for template ligand selection. With a sequence identity cutoff of 30% of target to templates and modeled target structures, FINDSITE^comb2.0 is shown to significantly improve upon FINDSITE^comb on the DUD-E benchmark set by increasing the 1% enrichment factor from 16.7 to 22.1, with a p-value of 4.3 × 10^-3 by the Student t-test. With an 80% sequence identity cutoff of target to templates for the DUD-E set and modeled target structures, FINDSITE^comb2.0, having a 1% ROC enrichment factor of 52.39, also outperforms state-of-the-art methods that employ machine learning such as a deep convolutional neural network, CNN, with an enrichment of 29.65. Thus, FINDSITE^comb2.0 represents a significant improvement in the state-of-the-art. The FINDSITE^comb2.0 web service is freely available for academic users at http://pwp.gatech.edu/cssb/FINDSITE-COMB-2 .

Docking Methodologies and Recent Advances

Docking, a molecular modelling method, has wide applications in identification and optimization in modern drug discovery. This chapter addresses the recent advances in the docking methodologies like fragment docking, covalent docking, inverse docking, post processing, hybrid techniques, homology modeling etc. and its protocol like searching and scoring functions. Advances in scoring functions for e.g. consensus scoring, quantum mechanics methods, clustering and entropy based methods, fingerprinting, etc. are used to overcome the limitations of the commonly used force-field, empirical and knowledge based scoring functions. It will cover crucial necessities and different algorithms of docking and scoring. Further different aspects like protein flexibility, ligand sampling and flexibility, and the performance of scoring function will be discussed.

Predicting Ligand Binding Sites on Protein Surfaces by 3-Dimensional Probability Density Distributions of Interacting Atoms

Predicting ligand binding sites (LBSs) on protein structures, which are obtained either from experimental or computational methods, is a useful first step in functional annotation or structure-based drug design for the protein structures. In this work, the structure-based machine learning algorithm ISMBLab-LIG was developed to predict LBSs on protein surfaces with input attributes derived from the three-dimensional probability density maps of interacting atoms, which were reconstructed on the query protein surfaces and were relatively insensitive to local conformational variations of the tentative ligand binding sites. The prediction accuracy of the ISMBLab-LIG predictors is comparable to that of the best LBS predictors benchmarked on several well-established testing datasets. More importantly, the ISMBLab-LIG algorithm has substantial tolerance to the prediction uncertainties of computationally derived protein structure models. As such, the method is particularly useful for predicting LBSs not only on experimental protein structures without known LBS templates in the database but also on computationally predicted model protein structures with structural uncertainties in the tentative ligand binding sites.

GeauxDock: A novel approach for mixed-resolution ligand docking using a descriptor-based force field

Molecular docking is an important component of computer-aided drug discovery. In this communication, we describe GeauxDock, a new docking approach that builds on the ideas of ligand homology modeling. GeauxDock features a descriptor-based scoring function integrating evolutionary constraints with physics-based energy terms, a mixed-resolution molecular representation of protein-ligand complexes, and an efficient Monte Carlo sampling protocol. To drive docking simulations toward experimental conformations, the scoring function was carefully optimized to produce a correlation between the total pseudoenergy and the native-likeness of binding poses. Indeed, benchmarking calculations demonstrate that GeauxDock has a strong capacity to identify near-native conformations across docking trajectories with the area under receiver operating characteristics of 0.85. By excluding closely related templates, we show that GeauxDock maintains its accuracy at lower levels of homology through the increased contribution from physics-based energy terms compensating for weak evolutionary constraints. GeauxDock is available at http://www.institute.loni.org/lasigma/package/dock/.

Pocketome of human kinases: prioritizing the ATP binding sites of (yet) untapped protein kinases for drug discovery

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.

Virtual screening: an in silico tool for interlacing the chemical universe with the proteome

In silico screening both in the forward (traditional virtual screening) and reverse sense (inverse virtual screening (IVS)) are helpful techniques for interlacing the chemical universe of small molecules with the proteome. The former, which is using a protein structure and a large chemical database, is well-known by the scientific community. We have chosen here to provide an overview on the latter, focusing on validation and target prioritization strategies. By comparing it to complementary or alternative wet-lab approaches, we put IVS in the broader context of chemical genomics, target discovery and drug design. By giving examples from the literature and an own example on how to validate the approach, we provide guidance on the issues related to IVS.

AFAL: a web service for profiling amino acids surrounding ligands in proteins

With advancements in crystallographic technology and the increasing wealth of information populating structural databases, there is an increasing need for prediction tools based on spatial information that will support the characterization of proteins and protein-ligand interactions. Herein, a new web service is presented termed amino acid frequency around ligand (AFAL) for determining amino acids type and frequencies surrounding ligands within proteins deposited in the Protein Data Bank and for assessing the atoms and atom-ligand distances involved in each interaction (availability: http://structuralbio.utalca.cl/AFAL/index.html ). AFAL allows the user to define a wide variety of filtering criteria (protein family, source organism, resolution, sequence redundancy and distance) in order to uncover trends and evolutionary differences in amino acid preferences that define interactions with particular ligands. Results obtained from AFAL provide valuable statistical information about amino acids that may be responsible for establishing particular ligand-protein interactions. The analysis will enable investigators to compare ligand-binding sites of different proteins and to uncover general as well as specific interaction patterns from existing data. Such patterns can be used subsequently to predict ligand binding in proteins that currently have no structural information and to refine the interpretation of existing protein models. The application of AFAL is illustrated by the analysis of proteins interacting with adenosine-5'-triphosphate.

Modelling three-dimensional protein structures for applications in drug design

A structural perspective of drug target and anti-target proteins, and their molecular interactions with biologically active molecules, largely advances many areas of drug discovery, including target validation, hit and lead finding and lead optimisation. In the absence of experimental 3D structures, protein structure prediction often offers a suitable alternative to facilitate structure-based studies. This review outlines recent methodical advances in homology modelling, with a focus on those techniques that necessitate consideration of ligand binding. In this context, model quality estimation deserves special attention because the accuracy and reliability of different structure prediction techniques vary considerably, and the quality of a model ultimately determines its usefulness for structure-based drug discovery. Examples of G-protein-coupled receptors (GPCRs) and ADMET-related proteins were selected to illustrate recent progress and current limitations of protein structure prediction. Basic guidelines for good modelling practice are also provided.

Proteome-scale docking: myth and reality

Docking is the computational method of choice to quickly predict how a low molecular-weight ligand binds to its macromolecular target. Despite persistent problems in predicting binding free energies, docking has undergone significant advances in numerous topics (throughput, target flexibility). The ever increasing availability of high-resolution X-ray structures and the development of more reliable comparative models for proteins of pharmacological interest paved the way to apply protein–ligand docking to multiple targets to predict main and off-targets for bioactive compounds and even to repurpose existing drugs. Applying docking to multiple targets brings an additional level of complexity in scoring numerous and heterogeneous docking poses. Despite undeniable successes, proteomewide docking should, however, be considered with caution with regard to recall and precision of the predictions.

A new target for an old drug: identifying mitoxantrone as a nanomolar inhibitor of PIM1 kinase via kinome-wide selectivity modeling

The rational design of selective kinase inhibitors remains a great challenge. Here we describe a physics-based approach to computationally modeling the kinase inhibitor selectivity profile. We retrospectively assessed this protocol by computing the binding profiles of 17 well-known kinase inhibitors against 143 kinases. Next, we predicted the binding profile of the chemotherapy drug mitoxantrone, and chose the predicted top five kinase targets for in vitro kinase assays. Remarkably, mitoxantrone was shown to possess low nanomolar inhibitory activity against PIM1 kinase and to inhibit the PIM1-mediated phosphorylation in cancer cells. We further determined the crystal complex structure of PIM1 bound with mitoxantrone, which reveals the structural and mechanistic basis for a novel mode of PIM1 inhibition. Although mitoxantrone's mechanism of action had been originally thought to act through DNA intercalation and type II topoisomerase inhibition, we hypothesize that PIM1 kinase inhibition might also contribute to mitoxantrone's therapeutic efficacy and specificity.

Are predicted protein structures of any value for binding site prediction and virtual ligand screening?

The recently developed field of ligand homology modeling (LHM) that extends the ideas of protein homology modeling to the prediction of ligand binding sites and for use in virtual ligand screening has emerged as a powerful new approach. Unlike traditional docking methodologies, LHM can be applied to low-to-moderate resolution predicted as well as experimental structures with little if any diminution in performance; thereby enabling ≈ 75% of an average proteome to have potentially significant virtual screening predictions. In large scale benchmarking, LHM is able to predict off-target ligand binding. Thus, despite the widespread belief to the contrary, low-to-moderate resolution predicted structures have considerable utility for biochemical function prediction.

Finding protein targets for small biologically relevant ligands across fold space using inverse ligand binding predictions

Inverse ligand binding prediction utilizes a few protein-ligand (drug) complexes to predict other secondary therapeutic and off-targets of a given drug molecule on a proteomic scale. We adapt two binding site predictors, FINDSITE and SMAP, to perform the inverse predictions and evaluate them on over 30 representative ligands. Use of just one complex allows the identification of other protein targets; the availability of additional complexes improves the results. Both methods offer comparable quality when using three complexes with diverse proteins. SMAP is better when fewer complexes are available, while FINDSITE provides stronger predictions for smaller ligands. We propose a consensus that combines (and outperforms) the two complementary approaches implemented by FINDSITE and SMAP. Most importantly, we demonstrate that these methods successfully find distant targets that belong to structurally different folds compared to the proteins in the input complexes.

Kinome-wide activity modeling from diverse public high-quality data sets

Large corpora of kinase small molecule inhibitor data are accessible to public sector research from thousands of journal article and patent publications. These data have been generated employing a wide variety of assay methodologies and experimental procedures by numerous laboratories. Here we ask the question how applicable these heterogeneous data sets are to predict kinase activities and which characteristics of the data sets contribute to their utility. We accessed almost 500,000 molecules from the Kinase Knowledge Base (KKB) and after rigorous aggregation and standardization generated over 180 distinct data sets covering all major groups of the human kinome. To assess the value of the data sets, we generated hundreds of classification and regression models. Their rigorous cross-validation and characterization demonstrated highly predictive classification and quantitative models for the majority of kinase targets if a minimum required number of active compounds or structure-activity data points were available. We then applied the best classifiers to compounds most recently profiled in the NIH Library of Integrated Network-based Cellular Signatures (LINCS) program and found good agreement of profiling results with predicted activities. Our results indicate that, although heterogeneous in nature, the publically accessible data sets are exceedingly valuable and well suited to develop highly accurate predictors for practical Kinome-wide virtual screening applications and to complement experimental kinase profiling.

FINDSITE(comb): a threading/structure-based, proteomic-scale virtual ligand screening approach

Virtual ligand screening is an integral part of the modern drug discovery process. Traditional ligand-based, virtual screening approaches are fast but require a set of structurally diverse ligands known to bind to the target. Traditional structure-based approaches require high-resolution target protein structures and are computationally demanding. In contrast, the recently developed threading/structure-based FINDSITE-based approaches have the advantage that they are as fast as traditional ligand-based approaches and yet overcome the limitations of traditional ligand- or structure-based approaches. These new methods can use predicted low-resolution structures and infer the likelihood of a ligand binding to a target by utilizing ligand information excised from the target's remote or close homologous proteins and/or libraries of ligand binding databases. Here, we develop an improved version of FINDSITE, FINDSITE(filt), that filters out false positive ligands in threading identified templates by a better binding site detection procedure that includes information about the binding site amino acid similarity. We then combine FINDSITE(filt) with FINDSITE(X) that uses publicly available binding databases ChEMBL and DrugBank for virtual ligand screening. The combined approach, FINDSITE(comb), is compared to two traditional docking methods, AUTODOCK Vina and DOCK 6, on the DUD benchmark set. It is shown to be significantly better in terms of enrichment factor, dependence on target structure quality, and speed. FINDSITE(comb) is then tested for virtual ligand screening on a large set of 3576 generic targets from the DrugBank database as well as a set of 168 Human GPCRs. Excluding close homologues, FINDSITE(comb) gives an average enrichment factor of 52.1 for generic targets and 22.3 for GPCRs within the top 1% of the screened compound library. Around 65% of the targets have better than random enrichment factors. The performance is insensitive to target structure quality, as long as it has a TM-score ≥ 0.4 to native. Thus, FINDSITE(comb) makes the screening of millions of compounds across entire proteomes feasible. The FINDSITE(comb) web service is freely available for academic users at http://cssb.biology.gatech.edu/skolnick/webservice/FINDSITE-COMB/index.html.

Molecular simulations of drug–receptor complexes in anticancer research

Molecular modeling and computer simulation techniques have matured significantly in recent years and proved their value in the study of drug–DNA, drug–DNA–protein, drug–protein and protein–protein interactions. Evolution in this area has gone hand-in-hand with an increased availability of structural data on biological macromolecules, major advances in molecular mechanics force fields and considerable improvements in computer technologies, most significantly processing speeds, multiprocessor programming and data-storage capacity. The information derived from molecular simulations of drug–receptor complexes can be used to extract structural and energetic information that is usually beyond current experimental possibilities, provide independent accounts of experimentally observed behavior, help in the interpretation of biochemical or pharmacological results, and open new avenues for research by posing novel relevant questions that can guide the design of new experiments. As drug-screening tools, ligand- and fragment-docking platforms stand out as powerful techniques that can provide candidate molecules for hit and lead development. This review provides an overall perspective of the main methods and focuses on some selected applications to both classical and novel anticancer targets.

FINDSITE(X): a structure-based, small molecule virtual screening approach with application to all identified human GPCRs

We have developed FINDSITE(X), an extension of FINDSITE, a protein threading based algorithm for the inference of protein binding sites, biochemical function and virtual ligand screening, that removes the limitation that holo protein structures (those containing bound ligands) of a sufficiently large set of distant evolutionarily related proteins to the target be solved; rather, predicted protein structures and experimental ligand binding information are employed. To provide the predicted protein structures, a fast and accurate version of our recently developed TASSER(VMT), TASSER(VMT)-lite, for template-based protein structural modeling applicable up to 1000 residues is developed and tested, with comparable performance to the top CASP9 servers. Then, a hybrid approach that combines structure alignments with an evolutionary similarity score for identifying functional relationships between target and proteins with binding data has been developed. By way of illustration, FINDSITE(X) is applied to 998 identified human G-protein coupled receptors (GPCRs). First, TASSER(VMT)-lite provides updates of all human GPCR structures previously modeled in our lab. We then use these structures and the new function similarity detection algorithm to screen all human GPCRs against the ZINC8 nonredundant (TC < 0.7) ligand set combined with ligands from the GLIDA database (a total of 88,949 compounds). Testing (excluding GPCRs whose sequence identity > 30% to the target from the binding data library) on a 168 human GPCR set with known binding data, the average enrichment factor in the top 1% of the compound library (EF(0.01)) is 22.7, whereas EF(0.01) by FINDSITE is 7.1. For virtual screening when just the target and its native ligands are excluded, the average EF(0.01) reaches 41.4. We also analyze off-target interactions for the 168 protein test set. All predicted structures, virtual screening data and off-target interactions for the 998 human GPCRs are available at http://cssb.biology.gatech.edu/skolnick/webservice/gpcr/index.html .

Structural modelling and dynamics of proteins for insights into drug interactions

Proteins are the workhorses of biomolecules and their function is affected by their structure and their structural rearrangements during ligand entry, ligand binding and protein-protein interactions. Hence, the knowledge of protein structure and, importantly, the dynamic behaviour of the structure are critical for understanding how the protein performs its function. The predictions of the structure and the dynamic behaviour can be performed by combinations of structure modelling and molecular dynamics simulations. The simulations also need to be sensitive to the constraints of the environment in which the protein resides. Standard computational methods now exist in this field to support the experimental effort of solving protein structures. This review presents a comprehensive overview of the basis of the calculations and the well-established computational methods used to generate and understand protein structure and function and the study of their dynamic behaviour with the reference to lung-related targets.

Virtual screening filters for the design of type II p38 MAP kinase inhibitors: a fragment based library generation approach

In this work, we introduce the development and application of a three-step scoring and filtering procedure for the design of type II p38 MAP kinase leads using allosteric fragments extracted from virtual screening hits. The design of the virtual screening filters is based on a thorough evaluation of docking methods, DFG-loop conformation, binding interactions and chemotype specificity of the 138 p38 MAP kinase inhibitors from Protein Data Bank bound to DFG-in and DFG-out conformations using Glide, GOLD and CDOCKER. A 40 ns molecular dynamics simulation with the apo, type I with DFG-in and type II with DFG-out forms was carried out to delineate the effects of structural variations on inhibitor binding. The designed docking-score and sub-structure filters were first tested on a dataset of 249 potent p38 MAP kinase inhibitors from seven diverse series and 18,842 kinase inhibitors from PDB, to gauge their capacity to discriminate between kinase and non-kinase inhibitors and likewise to selectively filter-in target-specific inhibitors. The designed filters were then applied in the virtual screening of a database of ten million (10⁷) compounds resulting in the identification of 100 hits. Based on their binding modes, 98 allosteric fragments were extracted from the hits and a fragment library was generated. New type II p38 MAP kinase leads were designed by tailoring the existing type I ATP site binders with allosteric fragments using a common urea linker. Target specific virtual screening filters can thus be easily developed for other kinases based on this strategy to retrieve target selective compounds.

Evaluation of model quality predictions in CASP9

CASP has been assessing the state of the art in the a priori estimation of accuracy of protein structure prediction since 2006. The inclusion of model quality assessment category in CASP contributed to a rapid development of methods in this area. In the last experiment, 46 quality assessment groups tested their approaches to estimate the accuracy of protein models as a whole and/or on a per-residue basis. We assessed the performance of these methods predominantly on the basis of the correlation between the predicted and observed quality of the models on both global and local scales. The ability of the methods to identify the models closest to the best one, to differentiate between good and bad models, and to identify well modeled regions was also analyzed. Our evaluations demonstrate that even though global quality assessment methods seem to approach perfection point (weighted average per-target Pearson's correlation coefficients are as high as 0.97 for the best groups), there is still room for improvement. First, all top-performing methods use consensus approaches to generate quality estimates, and this strategy has its own limitations. Second, the methods that are based on the analysis of individual models lag far behind clustering techniques and need a boost in performance. The methods for estimating per-residue accuracy of models are less accurate than global quality assessment methods, with an average weighted per-model correlation coefficient in the range of 0.63-0.72 for the best 10 groups.

Prediction of ligand binding using an approach designed to accommodate diversity in protein-ligand interactions

Computational determination of protein-ligand interaction potential is important for many biological applications including virtual screening for therapeutic drugs. The novel internal consensus scoring strategy is an empirical approach with an extended set of 9 binding terms combined with a neural network capable of analysis of diverse complexes. Like conventional consensus methods, internal consensus is capable of maintaining multiple distinct representations of protein-ligand interactions. In a typical use the method was trained using ligand classification data (binding/no binding) for a single receptor. The internal consensus analyses successfully distinguished protein-ligand complexes from decoys (r², 0.895 for a series of typical proteins). Results are superior to other tested empirical methods. In virtual screening experiments, internal consensus analyses provide consistent enrichment as determined by ROC-AUC and pROC metrics.

Synthesis, structural elucidation, DNA-PK inhibition, homology modelling and anti-platelet activity of morpholino-substituted-1,3-naphth-oxazines

A number of new angular 2-morpholino-(substituted)-naphth-1,3-oxazines (compound 10b), linear 2-morpholino-(substituted)-naphth-1,3-oxazines (compounds 13b-c), linear 6, 7 and 9-O-substituted-2-morpholino-(substituted)-naphth-1,3-oxazines (compounds 17-22, 24, and 25) and angular compounds 14-16 and 23 were synthesised. The O-substituent was pyridin-2yl-methyl (15, 18, and 21) pyridin-3yl-methyl (16, 19, and 22) and 4-methylpipreazin-1-yl-ethoxy (23-25). Twelve compounds were tested for their inhibitory effect on collagen induced platelet aggregation and it was found that the most active compounds were compounds 19 and 22 with IC(50)=55±4 and 85±4 μM, respectively. Furthermore, the compounds were also assayed for their ability to inhibit DNA-dependent protein kinase (DNA-PK) activity. The most active compounds were 18 IC(50)=0.091 μM, 24 IC(50)=0.191 μM, and 22 IC(50)=0.331 μM. Homology modelling was used to build a 3D model of DNA-PK based on the X-ray structure of phosphatidylinositol 3-kinases (PI3Ks). Docking of synthesised compounds within the binding pocket and structure-activity relationships (SAR) analyses of the poses were performed and results agreed well with observed activity.

Regulation of the water channel aquaporin-2 by posttranslational modification

The cellular functions of many eukaryotic membrane proteins, including the vasopressin-regulated water channel aquaporin-2 (AQP2), are regulated by posttranslational modifications. In this article, we discuss the experimental discoveries that have advanced our understanding of how posttranslational modifications affect AQP2 function, especially as they relate to the role of AQP2 in the kidney. We review the most recent data demonstrating that glycosylation and, in particular, phosphorylation and ubiquitination are mechanisms that regulate AQP2 activity, subcellular sorting and distribution, degradation, and protein interactions. From a clinical perspective, posttranslational modification resulting in protein misrouting or degradation may explain certain forms of nephrogenic diabetes insipidus. In addition to providing major insight into the function and dynamics of renal AQP2 regulation, the analysis of AQP2 posttranslational modification may provide general clues as to the role of posttranslational modification for regulation of other membrane proteins.

Cross-reactivity virtual profiling of the human kinome by X-react(KIN): a chemical systems biology approach

Many drug candidates fail in clinical development due to their insufficient selectivity that may cause undesired side effects. Therefore, modern drug discovery is routinely supported by computational techniques, which can identify alternate molecular targets with a significant potential for cross-reactivity. In particular, the development of highly selective kinase inhibitors is complicated by the strong conservation of the ATP-binding site across the kinase family. In this paper, we describe X-React(KIN), a new machine learning approach that extends the modeling and virtual screening of individual protein kinases to a system level in order to construct a cross-reactivity virtual profile for the human kinome. To maximize the coverage of the kinome, X-React(KIN) relies solely on the predicted target structures and employs state-of-the-art modeling techniques. Benchmark tests carried out against available selectivity data from high-throughput kinase profiling experiments demonstrate that, for almost 70% of the inhibitors, their alternate molecular targets can be effectively identified in the human kinome with a high (>0.5) sensitivity at the expense of a relatively low false positive rate (<0.5). Furthermore, in a case study, we demonstrate how X-React(KIN) can support the development of selective inhibitors by optimizing the selection of kinase targets for small-scale counter-screen experiments. The constructed cross-reactivity profiles for the human kinome are freely available to the academic community at http://cssb.biology.gatech.edu/kinomelhm/ .