PoSSuM v.3: A Major Expansion of the PoSSuM Database for Finding Similar Binding Sites of Proteins.

Information on structures of protein-ligand complexes, including comparisons of known and putative protein-ligand-binding pockets, is valuable for protein annotation and drug discovery and development. To facilitate biomedical and pharmaceutical research, we developed PoSSuM (https://possum.cbrc.pj.aist.go.jp/PoSSuM/), a database for identifying similar binding pockets in proteins. The current PoSSuM database includes 191 million similar pairs among almost 10 million identified pockets. PoSSuM drug search (PoSSuMds) is a resource for investigating ligand and receptor diversity among a set of pockets that can bind to an approved drug compound. The enhanced PoSSuMds covers pockets associated with both approved drugs and drug candidates in clinical trials from the latest release of ChEMBL. Additionally, we developed two new databases: PoSSuMAg for investigating antibody-antigen interactions and PoSSuMAF to simplify exploring putative pockets in AlphaFold human protein models.

PoSSuM v.2.0: data update and a new function for investigating ligand analogs and target proteins of small-molecule drugs.

PoSSuM (http://possum.cbrc.jp/PoSSuM/) is a database for detecting similar small-molecule binding sites on proteins. Since its initial release in 2011, PoSSuM has grown to provide information related to 49 million pairs of similar binding sites discovered among 5.5 million known and putative binding sites. This enlargement of the database is expected to enhance opportunities for biological and pharmaceutical applications, such as predictions of new functions and drug discovery. In this release, we have provided a new service named PoSSuM drug search (PoSSuMds) at http://possum.cbrc.jp/PoSSuM/drug_search/, in which we selected 194 approved drug compounds retrieved from ChEMBL, and detected their known binding pockets and pockets that are similar to them. Users can access and download all of the search results via a new web interface, which is useful for finding ligand analogs as well as potential target proteins. Furthermore, PoSSuMds enables users to explore the binding pocket universe within PoSSuM. Additionally, we have improved the web interface with new functions, including sortable tables and a viewer for visualizing and downloading superimposed pockets.

EzCatDB: the enzyme reaction database, 2015 update.

The EzCatDB database (http://ezcatdb.cbrc.jp/EzCatDB/) has emphasized manual classification of enzyme reactions from the viewpoints of enzyme active-site structures and their catalytic mechanisms based on literature information, amino acid sequences of enzymes (UniProtKB) and the corresponding tertiary structures from the Protein Data Bank (PDB). Reaction types such as hydrolysis, transfer, addition, elimination, isomerization, hydride transfer and electron transfer have been included in the reaction classification, RLCP. This database includes information related to ligand molecules on the enzyme structures in the PDB data, classified in terms of cofactors, substrates, products and intermediates, which are also necessary to elucidate the catalytic mechanisms. Recently, the database system was updated. The 3D structures of active sites for each PDB entry can be viewed using Jmol or Rasmol software. Moreover, sequence search systems of two types were developed for the EzCatDB database: EzCat-BLAST and EzCat-FORTE. EzCat-BLAST is suitable for quick searches, adopting the BLAST algorithm, whereas EzCat-FORTE is more suitable for detecting remote homologues, adopting the algorithm for FORTE protein structure prediction software. Another system, EzMetAct, is also available to searching for major active-site structures in EzCatDB, for which PDB-formatted queries can be searched.

Antibody informatics for drug discovery.

More and more antibody therapeutics are being approved every year, mainly due to their high efficacy and antigen selectivity. However, it is still difficult to identify the antigen, and thereby the function, of an antibody if no other information is available. There are obstacles inherent to the antibody science in every project in antibody drug discovery. Recent experimental technologies allow for the rapid generation of large-scale data on antibody sequences, affinity, potency, structures, and biological functions; this should accelerate drug discovery research. Therefore, a robust bioinformatic infrastructure for these large data sets has become necessary. In this article, we first identify and discuss the typical obstacles faced during the antibody drug discovery process. We then summarize the current status of three sub-fields of antibody informatics as follows: (i) recent progress in technologies for antibody rational design using computational approaches to affinity and stability improvement, as well as ab-initio and homology-based antibody modeling; (ii) resources for antibody sequences, structures, and immune epitopes and open drug discovery resources for development of antibody drugs; and (iii) antibody numbering and IMGT. Here, we review "antibody informatics," which may integrate the above three fields so that bridging the gaps between industrial needs and academic solutions can be accelerated. This article is part of a Special Issue entitled: Recent advances in molecular engineering of antibody.

Kotai Antibody Builder: automated high-resolution structural modeling of antibodies.

MOTIVATION: Kotai Antibody Builder is a Web service for tertiary structural modeling of antibody variable regions. It consists of three main steps: hybrid template selection by sequence alignment and canonical rules, 3D rendering of alignments and CDR-H3 loop modeling. For the last step, in addition to rule-based heuristics used to build the initial model, a refinement option is available that uses fragment assembly followed by knowledge-based scoring. Using targets from the Second Antibody Modeling Assessment, we demonstrate that Kotai Antibody Builder generates models with an overall accuracy equal to that of the best-performing semi-automated predictors using expert knowledge. AVAILABILITY AND IMPLEMENTATION: Kotai Antibody Builder is available at http://kotaiab.org CONTACT: standley@ifrec.osaka-u.ac.jp.

High-resolution modeling of antibody structures by a combination of bioinformatics, expert knowledge, and molecular simulations.

In the second antibody modeling assessment, we used a semiautomated template-based structure modeling approach for 11 blinded antibody variable region (Fv) targets. The structural modeling method involved several steps, including template selection for framework and canonical structures of complementary determining regions (CDRs), homology modeling, energy minimization, and expert inspection. The submitted models for Fv modeling in Stage 1 had the lowest average backbone root mean square deviation (RMSD) (1.06 Å). Comparison to crystal structures showed the most accurate Fv models were generated for 4 out of 11 targets. We found that the successful modeling in Stage 1 mainly was due to expert-guided template selection for CDRs, especially for CDR-H3, based on our previously proposed empirical method (H3-rules) and the use of position specific scoring matrix-based scoring. Loop refinement using fragment assembly and multicanonical molecular dynamics (McMD) was applied to CDR-H3 loop modeling in Stage 2. Fragment assembly and McMD produced putative structural ensembles with low free energy values that were scored based on the OSCAR all-atom force field and conformation density in principal component analysis space, respectively, as well as the degree of consensus between the two sampling methods. The quality of 8 out of 10 targets improved as compared with Stage 1. For 4 out of 10 Stage-2 targets, our method generated top-scoring models with RMSD values of less than 1 Å. In this article, we discuss the strengths and weaknesses of our approach as well as possible directions for improvement to generate better predictions in the future.

A new workflow for the effective curation of membrane permeability data from open ADME information.

Membrane permeability is an in vitro parameter that represents the apparent permeability (Papp) of a compound, and is a key absorption, distribution, metabolism, and excretion parameter in drug development. Although the Caco-2 cell lines are the most used cell lines to measure Papp, other cell lines, such as the Madin-Darby Canine Kidney (MDCK), LLC-Pig Kidney 1 (LLC-PK1), and Ralph Russ Canine Kidney (RRCK) cell lines, can also be used to estimate Papp. Therefore, constructing in silico models for Papp estimation using the MDCK, LLC-PK1, and RRCK cell lines requires collecting extensive amounts of in vitro Papp data. An open database offers extensive measurements of various compounds covering a vast chemical space; however, concerns were reported on the use of data published in open databases without the appropriate accuracy and quality checks. Ensuring the quality of datasets for training in silico models is critical because artificial intelligence (AI, including deep learning) was used to develop models to predict various pharmacokinetic properties, and data quality affects the performance of these models. Hence, careful curation of the collected data is imperative. Herein, we developed a new workflow that supports automatic curation of Papp data measured in the MDCK, LLC-PK1, and RRCK cell lines collected from ChEMBL using KNIME. The workflow consisted of four main phases. Data were extracted from ChEMBL and filtered to identify the target protocols. A total of 1661 high-quality entries were retained after checking 436 articles. The workflow is freely available, can be updated, and has high reusability. Our study provides a novel approach for data quality analysis and accelerates the development of helpful in silico models for effective drug discovery. Scientific Contribution: The cost of building highly accurate predictive models can be significantly reduced by automating the collection of reliable measurement data. Our tool reduces the time and effort required for data collection and will enable researchers to focus on constructing high-performance in silico models for other types of analysis. To the best of our knowledge, no such tool is available in the literature.

THOUSAND-GRAIN WEIGHT 6, which is an IAA-glucose hydrolase, preferentially recognizes the structure of the indole ring.

An indole-3-acetic acid (IAA)-glucose hydrolase, THOUSAND-GRAIN WEIGHT 6 (TGW6), negatively regulates the grain weight in rice. TGW6 has been used as a target for breeding increased rice yield. Moreover, the activity of TGW6 has been thought to involve auxin homeostasis, yet the details of this putative TGW6 activity remain unclear. Here, we show the three-dimensional structure and substrate preference of TGW6 using X-ray crystallography, thermal shift assays and fluorine nuclear magnetic resonance (19F NMR). The crystal structure of TGW6 was determined at 2.6 Å resolution and exhibited a six-bladed ß-propeller structure. Thermal shift assays revealed that TGW6 preferably interacted with indole compounds among the tested substrates, enzyme products and their analogs. Further analysis using 19F NMR with 1,134 fluorinated fragments emphasized the importance of indole fragments in recognition by TGW6. Finally, docking simulation analyses of the substrate and related fragments in the presence of TGW6 supported the interaction specificity for indole compounds. Herein, we describe the structure and substrate preference of TGW6 for interacting with indole fragments during substrate recognition. Uncovering the molecular details of TGW6 activity will stimulate the use of this enzyme for increasing crop yields and contributes to functional studies of IAA glycoconjugate hydrolases in auxin homeostasis.

Applying deep learning to iterative screening of medium-sized molecules for protein-protein interaction-targeted drug discovery.

We combined a library of medium-sized molecules with iterative screening using multiple machine learning algorithms that were ligand-based, which resulted in a large increase of the hit rate against a protein-protein interaction target. This was demonstrated by inhibition assays using a PPI target, Kelch-like ECH-associated protein 1/nuclear factor erythroid 2-related factor 2 (Keap1/Nrf2), and a deep neural network model based on the first-round assay data showed a highest hit rate of 27.3%. Using the models, we identified novel active and non-flat compounds far from public datasets, expanding the chemical space.

AI-driven molecular generation of not-patented pharmaceutical compounds using world open patent data.

Developing compounds with novel structures is important for the production of new drugs. From an intellectual perspective, confirming the patent status of newly developed compounds is essential, particularly for pharmaceutical companies. The generation of a large number of compounds has been made possible because of the recent advances in artificial intelligence (AI). However, confirming the patent status of these generated molecules has been a challenge because there are no free and easy-to-use tools that can be used to determine the novelty of the generated compounds in terms of patents in a timely manner; additionally, there are no appropriate reference databases for pharmaceutical patents in the world. In this study, two public databases, SureChEMBL and Google Patents Public Datasets, were used to create a reference database of drug-related patented compounds using international patent classification. An exact structure search system was constructed using InChIKey and a relational database system to rapidly search for compounds in the reference database. Because drug-related patented compounds are a good source for generative AI to learn useful chemical structures, they were used as the training data. Furthermore, molecule generation was successfully directed by increasing and decreasing the number of generated patented compounds through incorporation of patent status (i.e., patented or not) into learning. The use of patent status enabled generation of novel molecules with high drug-likeness. The generation using generative AI with patent information would help efficiently propose novel compounds in terms of pharmaceutical patents. Scientific contribution: In this study, a new molecule-generation method that takes into account the patent status of molecules, which has rarely been considered but is an important feature in drug discovery, was developed. The method enables the generation of novel molecules based on pharmaceutical patents with high drug-likeness and will help in the efficient development of effective drug compounds.

DLiP-PPI library: An integrated chemical database of small-to-medium-sized molecules targeting protein-protein interactions.

Protein-protein interactions (PPIs) are recognized as important targets in drug discovery. The characteristics of molecules that inhibit PPIs differ from those of small-molecule compounds. We developed a novel chemical library database system (DLiP) to design PPI inhibitors. A total of 32,647 PPI-related compounds are registered in the DLiP. It contains 15,214 newly synthesized compounds, with molecular weight ranging from 450 to 650, and 17,433 active and inactive compounds registered by extracting and integrating known compound data related to 105 PPI targets from public databases and published literature. Our analysis revealed that the compounds in this database contain unique chemical structures and have physicochemical properties suitable for binding to the protein-protein interface. In addition, advanced functions have been integrated with the web interface, which allows users to search for potential PPI inhibitor compounds based on types of protein-protein interfaces, filter results by drug-likeness indicators important for PPI targeting such as rule-of-4, and display known active and inactive compounds for each PPI target. The DLiP aids the search for new candidate molecules for PPI drug discovery and is available online (https://skb-insilico.com/dlip).

Collation and data-mining of literature bioactivity data for drug discovery.

The challenge of translating the huge amount of genomic and biochemical data into new drugs is a costly and challenging task. Historically, there has been comparatively little focus on linking the biochemical and chemical worlds. To address this need, we have developed ChEMBL, an online resource of small-molecule SAR (structure-activity relationship) data, which can be used to support chemical biology, lead discovery and target selection in drug discovery. The database contains the abstracted structures, properties and biological activities for over 700000 distinct compounds and in excess of more than 3 million bioactivity records abstracted from over 40000 publications. Additional public domain resources can be readily integrated into the same data model (e.g. PubChem BioAssay data). The compounds in ChEMBL are largely extracted from the primary medicinal chemistry literature, and are therefore usually 'drug-like' or 'lead-like' small molecules with full experimental context. The data cover a significant fraction of the discovery of modern drugs, and are useful in a wide range of drug design and discovery tasks. In addition to the compound data, ChEMBL also contains information for over 8000 protein, cell line and whole-organism 'targets', with over 4000 of those being proteins linked to their underlying genes. The database is searchable both chemically, using an interactive compound sketch tool, protein sequences, family hierarchies, SMILES strings, compound research codes and key words, and biologically, using a variety of gene identifiers, protein sequence similarity and protein families. The information retrieved can then be readily filtered and downloaded into various formats. ChEMBL can be accessed online at https://www.ebi.ac.uk/chembldb.

PreBINDS: An Interactive Web Tool to Create Appropriate Datasets for Predicting Compound-Protein Interactions.

Given the abundant computational resources and the huge amount of data of compound-protein interactions (CPIs), constructing appropriate datasets for learning and evaluating prediction models for CPIs is not always easy. For this study, we have developed a web server to facilitate the development and evaluation of prediction models by providing an appropriate dataset according to the task. Our web server provides an environment and dataset that aid model developers and evaluators in obtaining a suitable dataset for both proteins and compounds, in addition to attributes necessary for deep learning. With the web server interface, users can customize the CPI dataset derived from ChEMBL by setting positive and negative thresholds to be adjusted according to the user's definitions. We have also implemented a function for graphic display of the distribution of activity values in the dataset as a histogram to set appropriate thresholds for positive and negative examples. These functions enable effective development and evaluation of models. Furthermore, users can prepare their task-specific datasets by selecting a set of target proteins based on various criteria such as Pfam families, ChEMBL's classification, and sequence similarities. The accuracy and efficiency of in silico screening and drug design using machine learning including deep learning can therefore be improved by facilitating access to an appropriate dataset prepared using our web server (https://binds.lifematics.work/).

Identification of novel inhibitors of Keap1/Nrf2 by a promising method combining protein-protein interaction-oriented library and machine learning.

Protein-protein interactions (PPIs) are prospective but challenging targets for drug discovery, because screening using traditional small-molecule libraries often fails to identify hits. Recently, we developed a PPI-oriented library comprising 12,593 small-to-medium-sized newly synthesized molecules. This study validates a promising combined method using PPI-oriented library and ligand-based virtual screening (LBVS) to discover novel PPI inhibitory compounds for Kelch-like ECH-associated protein 1 (Keap1) and nuclear factor erythroid 2-related factor 2 (Nrf2). We performed LBVS with two random forest models against our PPI library and the following time-resolved fluorescence resonance energy transfer (TR-FRET) assays of 620 compounds identified 15 specific hit compounds. The high hit rates for the entire PPI library (estimated 0.56-1.3%) and the LBVS (maximum 5.4%) compared to a conventional screening library showed the utility of the library and the efficiency of LBVS. All the hit compounds possessed novel structures with Tanimoto similarity ≤ 0.26 to known Keap1/Nrf2 inhibitors and aqueous solubility (AlogP < 5). Reasonable binding modes were predicted using 3D alignment of five hit compounds and a Keap1/Nrf2 peptide crystal structure. Our results represent a new, efficient method combining the PPI library and LBVS to identify novel PPI inhibitory ligands with expanded chemical space.

Comprehensive Approach of 19F Nuclear Magnetic Resonance, Enzymatic, and In Silico Methods for Site-Specific Hit Selection and Validation of Fragment Molecules that Inhibit Methionine γ-Lyase Activity.

Fragment-based screening using 19F NMR (19F-FS) is an efficient method for exploring seed and lead compounds for drug discovery. Here, we demonstrate the utility and merits of using 19F-FS for methionine γ-lyase-binding fragments, together with a 19F NMR-based competition and mutation assay, as well as enzymatic and in silico methods. 19F NMR-based assays provided useful information on binding between 19F-FS hit fragments and target proteins. Although the 19F-FS and enzymatic assay were weakly correlated, they show that the 19F-FS hit fragments contained compounds with inhibitory activity. Furthermore, we found that in silico calculations partially account for the differences in activity levels between the 19F-FS hits as per NMR analysis. A comprehensive approach combining the 19F-FS and other methods not only identified fragment hits but also distinguished structural differences in chemical groups with diverse activity levels.

Mechanism of hERG inhibition by gating-modifier toxin, APETx1, deduced by functional characterization.

BACKGROUND: Human ether-à-go-go-related gene potassium channel 1 (hERG) is a voltage-gated potassium channel, the voltage-sensing domain (VSD) of which is targeted by a gating-modifier toxin, APETx1. APETx1 is a 42-residue peptide toxin of sea anemone Anthopleura elegantissima and inhibits hERG by stabilizing the resting state. A previous study that conducted cysteine-scanning analysis of hERG identified two residues in the S3-S4 region of the VSD that play important roles in hERG inhibition by APETx1. However, mutational analysis of APETx1 could not be conducted as only natural resources have been available until now. Therefore, it remains unclear where and how APETx1 interacts with the VSD in the resting state. RESULTS: We established a method for preparing recombinant APETx1 and determined the NMR structure of the recombinant APETx1, which is structurally equivalent to the natural product. Electrophysiological analyses using wild type and mutants of APETx1 and hERG revealed that their hydrophobic residues, F15, Y32, F33, and L34, in APETx1, and F508 and I521 in hERG, in addition to a previously reported acidic hERG residue, E518, play key roles in the inhibition of hERG by APETx1. Our hypothetical docking models of the APETx1-VSD complex satisfied the results of mutational analysis. CONCLUSIONS: The present study identified the key residues of APETx1 and hERG that are involved in hERG inhibition by APETx1. These results would help advance understanding of the inhibitory mechanism of APETx1, which could provide a structural basis for designing novel ligands targeting the VSDs of KV channels.

Assessing the druggability of protein-protein interactions by a supervised machine-learning method.

BACKGROUND: Protein-protein interactions (PPIs) are challenging but attractive targets of small molecule drugs for therapeutic interventions of human diseases. In this era of rapid accumulation of PPI data, there is great need for a methodology that can efficiently select drug target PPIs by holistically assessing the druggability of PPIs. To address this need, we propose here a novel approach based on a supervised machine-learning method, support vector machine (SVM). RESULTS: To assess the druggability of the PPIs, 69 attributes were selected to cover a wide range of structural, drug and chemical, and functional information on the PPIs. These attributes were used as feature vectors in the SVM-based method. Thirty PPIs known to be druggable were carefully selected from previous studies; these were used as positive instances. Our approach was applied to 1,295 human PPIs with tertiary structures of their protein complexes already solved. The best SVM model constructed discriminated the already-known target PPIs from others at an accuracy of 81% (sensitivity, 82%; specificity, 79%) in cross-validation. Among the attributes, the two with the greatest discriminative power in the best SVM model were the number of interacting proteins and the number of pathways. CONCLUSION: Using the model, we predicted several promising candidates for druggable PPIs, such as SMAD4/SKI. As more PPI data are accumulated in the near future, our method will have increased ability to accelerate the discovery of druggable PPIs.

Universal partitioning of the hierarchical fold network of 50-residue segments in proteins.

BACKGROUND: Several studies have demonstrated that protein fold space is structured hierarchically and that power-law statistics are satisfied in relation between the numbers of protein families and protein folds (or superfamilies). We examined the internal structure and statistics in the fold space of 50 amino-acid residue segments taken from various protein folds. We used inter-residue contact patterns to measure the tertiary structural similarity among segments. Using this similarity measure, the segments were classified into a number (Kc) of clusters. We examined various Kc values for the clustering. The special resolution to differentiate the segment tertiary structures increases with increasing Kc. Furthermore, we constructed networks by linking structurally similar clusters. RESULTS: The network was partitioned persistently into four regions for Kc >or= 1000. This main partitioning is consistent with results of earlier studies, where similar partitioning was reported in classifying protein domain structures. Furthermore, the network was partitioned naturally into several dozens of sub-networks (i.e., communities). Therefore, intra-sub-network clusters were mutually connected with numerous links, although inter-sub-network ones were rarely done with few links. For Kc >or= 1000, the major sub-networks were about 40; the contents of the major sub-networks were conserved. This sub-partitioning is a novel finding, suggesting that the network is structured hierarchically: Segments construct a cluster, clusters form a sub-network, and sub-networks constitute a region. Additionally, the network was characterized by non-power-law statistics, which is also a novel finding. CONCLUSION: Main findings are: (1) The universe of 50 residue segments found here was characterized by non-power-law statistics. Therefore, the universe differs from those ever reported for the protein domains. (2) The 50-residue segments were partitioned persistently and universally into some dozens (ca. 40) of major sub-networks, irrespective of the number of clusters. (3) These major sub-networks encompassed 90% of all segments. Consequently, the protein tertiary structure is constructed using the dozens of elements (sub-networks).

Structural bases of GM1 gangliosidosis and Morquio B disease.

Allelic mutations of the lysosomal beta-galactosidase gene cause heterogeneous clinical phenotypes, such as GM1 gangliosidosis and Morquio B disease, the former being further classified into three variants, namely infantile, juvenile and adult forms; and heterogeneous biochemical phenotypes were shown in these forms. We tried to elucidate the bases of these diseases from a structural viewpoint. We first constructed a three-dimensional structural model of human beta-galactosidase by means of homology modeling. The human beta-galactosidase consists of three domains, such as, a TIM barrel fold domain, which functions as a catalytic domain, and two galactose-binding domain-like fold domains. We then constructed structural models of representative mutant beta-galactosidase proteins (G123R, R201C, I51T and Y83H) and predicted the structural change associated with each phenotype by calculating the number of affected atoms, determining the root-mean-square deviation and the solvent-accessible surface area, and by color imaging. The results show that there is a good correlation between the structural changes caused by amino-acid substitutions in the beta-galactosidase molecule, as well as biochemical and clinical phenotypes in these representative cases. Protein structural study is useful for elucidating the bases of these diseases.

A prospective compound screening contest identified broader inhibitors for Sirtuin 1.

Potential inhibitors of a target biomolecule, NAD-dependent deacetylase Sirtuin 1, were identified by a contest-based approach, in which participants were asked to propose a prioritized list of 400 compounds from a designated compound library containing 2.5 million compounds using in silico methods and scoring. Our aim was to identify target enzyme inhibitors and to benchmark computer-aided drug discovery methods under the same experimental conditions. Collecting compound lists derived from various methods is advantageous for aggregating compounds with structurally diversified properties compared with the use of a single method. The inhibitory action on Sirtuin 1 of approximately half of the proposed compounds was experimentally accessed. Ultimately, seven structurally diverse compounds were identified.