PgpRules: a decision tree based prediction server for P-glycoprotein substrates and inhibitors.

SUMMARY: P-glycoprotein (P-gp) is a member of ABC transporter family that actively pumps xenobiotics out of cells to protect organisms from toxic compounds. P-gp substrates can be easily pumped out of the cells to reduce their absorption; conversely P-gp inhibitors can reduce such pumping activity. Hence, it is crucial to know if a drug is a P-gp substrate or inhibitor in view of pharmacokinetics. Here we present PgpRules, an online P-gp substrate and P-gp inhibitor prediction server with ruled-sets. The two models were built using classification and regression tree algorithm. For each compound uploaded, PgpRules not only predicts whether the compound is a P-gp substrate or a P-gp inhibitor, but also provides the rules containing chemical structural features for further structural optimization. AVAILABILITY AND IMPLEMENTATION: PgpRules is freely accessible at https://pgprules.cmdm.tw/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

G.A.M.E.: GPU-accelerated mixture elucidator.

GPU acceleration is useful in solving complex chemical information problems. Identifying unknown structures from the mass spectra of natural product mixtures has been a desirable yet unresolved issue in metabolomics. However, this elucidation process has been hampered by complex experimental data and the inability of instruments to completely separate different compounds. Fortunately, with current high-resolution mass spectrometry, one feasible strategy is to define this problem as extending a scaffold database with sidechains of different probabilities to match the high-resolution mass obtained from a high-resolution mass spectrum. By introducing a dynamic programming (DP) algorithm, it is possible to solve this NP-complete problem in pseudo-polynomial time. However, the running time of the DP algorithm grows by orders of magnitude as the number of mass decimal digits increases, thus limiting the boost in structural prediction capabilities. By harnessing the heavily parallel architecture of modern GPUs, we designed a "compute unified device architecture" (CUDA)-based GPU-accelerated mixture elucidator (G.A.M.E.) that considerably improves the performance of the DP, allowing up to five decimal digits for input mass data. As exemplified by four testing datasets with verified constitutions from natural products, G.A.M.E. allows for efficient and automatic structural elucidation of unknown mixtures for practical procedures. Graphical abstract .

Trowaglerix Venom Polypeptides As a Novel Antithrombotic Agent by Targeting Immunoglobulin-Like Domains of Glycoprotein VI in Platelet.

OBJECTIVE: Currently prescribed antiplatelet drugs have 1 common side effect-an increased risk of hemorrhage and thrombocytopenia. On the contrary, bleeding defects associated with glycoprotein VI (GPVI) expression deficiency are usually slightly prolonged bleeding times. However, GPVI antagonists are lacking in clinic. APPROACH AND RESULTS: Using reverse-phase high-performance liquid chromatography and sequencing, we revealed the partial sequence of trowaglerix α subunit, a potent specific GPVI-targeting snaclec (snake venom C-type lectin protein). Hexapeptide (Troα6 [trowaglerix a chain hexapeptide, CKWMNV]) and decapeptide (Troα10) derived from trowaglerix specifically inhibited collagen-induced platelet aggregation through blocking platelet GPVI receptor. Computational peptide design helped to design a series of Troα6/Troα10 peptides. Protein docking studies on these decapeptides and GPVI suggest that Troα10 was bound at the lower surface of D1 domain and outer surface of D2 domain, which was at the different place of the collagen-binding site and the scFv (single-chain variable fragment) D2-binding site. The newly discovered site was confirmed by inhibitory effects of polyclonal antibodies on collagen-induced platelet aggregation. This indicates that D2 domain of GPVI is a novel and important binding epitope on GPVI-mediated platelet aggregation. Troα6/Troα10 displayed prominent inhibitory effect of thrombus formation in fluorescein sodium-induced platelet thrombus formation of mesenteric venules and ferric chloride-induced carotid artery injury thrombosis model without prolonging the in vivo bleeding time. CONCLUSIONS: We develop a novel antithrombotic peptides derived from trowaglerix that acts through GPVI antagonism with greater safety-no severe bleeding. The binding epitope of polypeptides on GPVI is novel and important. These hexa/decapeptides have therapeutic potential for developing ideal small-mass GPVI antagonists for arterial thrombogenic diseases.

Effects of sarcosine and N, N-dimethylglycine on NMDA receptor-mediated excitatory field potentials.

BACKGROUND: Sarcosine, a glycine transporter type 1 inhibitor and an N-methyl-D-aspartate (NMDA) receptor co-agonist at the glycine binding site, potentiates NMDA receptor function. Structurally similar to sarcosine, N,N-dimethylglycine (DMG) is also N-methyl glycine-derivative amino acid and commonly used as a dietary supplement. The present study compared the effects of sarcosine and DMG on NMDA receptor-mediated excitatory field potentials (EFPs) in mouse medial prefrontal cortex brain slices using a multi-electrode array system. RESULTS: Glycine, sarcosine and DMG alone did not alter the NMDA receptor-mediated EFPs, but in combination with glutamate, glycine and its N-methyl derivatives significantly increased the frequency and amplitude of EFPs. The enhancing effects of glycine analogs in combination with glutamate on EFPs were remarkably reduced by the glycine binding site antagonist 7-chlorokynurenate (7-CK). However, DMG, but not sarcosine, reduced the frequency and amplitude of EFPs elicited by co-application of glutamate plus glycine. D-cycloserine, a partial agonist at the glycine binding site on NMDA receptors, affected EFPs in a similar manner to DMG. Furthermore, DMG, but not sarcosine, reduced the frequencies and amplitudes of EFPs elicited by glutamate plus D-serine, another endogenous ligand for glycine binding site. CONCLUSIONS: These findings suggest that sarcosine acts as a full agonist, yet DMG is a partial agonist at glycine binding site of NMDA receptors. The molecular docking analysis indicated that the interactions of glycine, sarcosine, and DMG to NMDA receptors are highly similar, supporting that the glycine binding site of NMDA receptors is a critical target site for sarcosine and DMG.

Mutagenicity in a Molecule: Identification of Core Structural Features of Mutagenicity Using a Scaffold Analysis.

With advances in the development and application of Ames mutagenicity in silico prediction tools, the International Conference on Harmonisation (ICH) has amended its M7 guideline to reflect the use of such prediction models for the detection of mutagenic activity in early drug safety evaluation processes. Since current Ames mutagenicity prediction tools only focus on functional group alerts or side chain modifications of an analog series, these tools are unable to identify mutagenicity derived from core structures or specific scaffolds of a compound. In this study, a large collection of 6512 compounds are used to perform scaffold tree analysis. By relating different scaffolds on constructed scaffold trees with Ames mutagenicity, four major and one minor novel mutagenic groups of scaffold are identified. The recognized mutagenic groups of scaffold can serve as a guide for medicinal chemists to prevent the development of potentially mutagenic therapeutic agents in early drug design or development phases, by modifying the core structures of mutagenic compounds to form non-mutagenic compounds. In addition, five series of substructures are provided as recommendations, for direct modification of potentially mutagenic scaffolds to decrease associated mutagenic activities.

Rule-Based Prediction Models of Cytochrome P450 Inhibition.

Hepatotoxicity, drug-induced liver injury, and competitive Cytochrome P-450 (CYP) isozyme binding are serious problems associated with drug use. It would be favorable to avoid or to understand potential CYP inhibition at the developmental stages. However, current in silico CYP prediction models or available public prediction servers can provide only yes/no classification results for just one or a few CYP enzymes. In this study, we utilized a rule-based C5.0 algorithm with different descriptors, including PaDEL, Mold(2), and PubChem fingerprints, to construct rule-based inhibition prediction models for five major CYP enzymes-CYP1A2, CYP2C9, CYP2C19, CYP2D6 and CYP3A4-that account for 90% of drug oxidation or hydrolysis. We also developed a rational sampling algorithm for the selection of compounds in the training data set, to enhance the performance of these CYP prediction models. The optimized models include several improved features. First, the final models significantly outperformed all of the currently available models. Second, the final models can also be used for rapid virtual screening of a large set of compounds due to their ruleset-based nature. Moreover, such rule-based prediction models can provide rulesets for structural features related to the five major CYP enzymes. The five most significant rules for CYP inhibition were identified for each CYP enzymes and discussed. An example was chosen for each of the five CYP enzymes to demonstrate how rule-based models can be used to gain insights into structural features that correspond with CYP inhibitions. A newer version of the freely accessible CYP prediction server, CypRules, is presented here as a result of the aforementioned improvements.

CypRules: a rule-based P450 inhibition prediction server.

UNLABELLED: Cytochrome P450 (CYPs) are the major enzymes involved in drug metabolism and bioactivation. Inhibition models were constructed for five of the most popular enzymes from the CYP superfamily in human liver. The five enzymes chosen for this study, namely CYP1A2, CYP2D6, CYP2C19, CYP2C9 and CYP3A4, account for 90% of the xenobiotic and drug metabolism in human body. CYP enzymes can be inhibited or induced by various drugs or chemical compounds. In this work, a rule-based CYP inhibition prediction online server, CypRules, was created based on predictive models generated by the rule-based C5.0 algorithm. CypRules can predict and provide structural rulesets for CYP inhibition for each compound uploaded to the server. Capable of fast execution performance, it can be used for virtual high-throughput screening (VHTS) of a large set of testing compounds. AVAILABILITY AND IMPLEMENTATION: CypRules is freely accessible at http://cyprules.cmdm.tw/ and models, descriptor and program files for all compounds are publically available at http://cyprules.cmdm.tw/sources/sources.rar.

Rule-based classification models of molecular autofluorescence.

Fluorescence-based detection has been commonly used in high-throughput screening (HTS) assays. Autofluorescent compounds, which can emit light in the absence of artificial fluorescent markers, often interfere with the detection of fluorophores and result in false positive signals in these assays. This interference presents a major issue in fluorescence-based screening techniques. In an effort to reduce the time and cost that will be spent on prescreening of autofluorescent compounds, in silico autofluorescence prediction models were developed for selected fluorescence-based assays in this study. Five prediction models were developed based on the respective fluorophores used in these HTS assays, which absorb and emit light at specific wavelengths (excitation/emission): Alexa Fluor 350 (A350) (340 nm/450 nm), 7-amino-4-trifluoromethyl-coumarin (AFC) (405 nm/520 nm), Alexa Fluor 488 (A488) (480 nm/540 nm), Rhodamine (547 nm/598 nm), and Texas Red (547 nm/618 nm). The C5.0 rule-based classification algorithm and PubChem 2D chemical structure fingerprints were used to develop prediction models. To optimize the accuracies of these prediction models despite the highly imbalanced ratio of fluorescent versus nonfluorescent compounds presented in the collected data sets, oversampling and undersampling strategies were applied. The average final accuracy achieved for the training set was 97%, and that for the testing set was 92%. In addition, five external data sets were used to further validate the models. Ultimately, 14 representative structural features (or rules) were determined to efficiently predict autofluorescence in data sets containing both fluorescent and nonfluorescent compounds. Several cases were illustrated in this study to demonstrate the applicability of these rules.

Template-based de novo design for type II kinase inhibitors and its extented application to acetylcholinesterase inhibitors.

There is a compelling need to discover type II inhibitors targeting the unique DFG-out inactive kinase conformation since they are likely to possess greater potency and selectivity relative to traditional type I inhibitors. Using a known inhibitor, such as a currently available and approved drug or inhibitor, as a template to design new drugs via computational de novo design is helpful when working with known ligand-receptor interactions. This study proposes a new template-based de novo design protocol to discover new inhibitors that preserve and also optimize the binding interactions of the type II kinase template. First, sorafenib (Nexavar) and nilotinib (Tasigna), two type II inhibitors with different ligand-receptor interactions, were selected as the template compounds. The five-step protocol can reassemble each drug from a large fragment library. Our procedure demonstrates that the selected template compounds can be successfully reassembled while the key ligand-receptor interactions are preserved. Furthermore, to demonstrate that the algorithm is able to construct more potent compounds, we considered kinase inhibitors and other protein dataset, acetylcholinesterase (AChE) inhibitors. The de novo optimization was initiated using a template compound possessing a less than optimal activity from a series of aminoisoquinoline and TAK-285 inhibiting type II kinases, and E2020 derivatives inhibiting AChE respectively. Three compounds with greater potency than the template compound were discovered that were also included in the original congeneric series. This template-based lead optimization protocol with the fragment library can help to design compounds with preferred binding interactions of known inhibitors automatically and further optimize the compounds in the binding pockets.

Predictive toxicology modeling: protocols for exploring hERG classification and Tetrahymena pyriformis end point predictions.

The inclusion and accessibility of different methodologies to explore chemical data sets has been beneficial to the field of predictive modeling, specifically in the chemical sciences in the field of Quantitative Structure-Activity Relationship (QSAR) modeling. This study discusses using contemporary protocols and QSAR modeling methods to properly model two biomolecular systems that have historically not performed well using traditional and three-dimensional QSAR methodologies. Herein, we explore, analyze, and discuss the creation of a classification human Ether-a-go-go Related Gene (hERG) potassium channel model and a continuous Tetrahymena pyriformis (T. pyriformis) model using Support Vector Machine (SVM) and Support Vector Regression (SVR), respectively. The models are constructed with three types of molecular descriptors that capture the gross physicochemical features of the compounds: (i) 2D, 2 1/2D, and 3D physical features, (ii) VolSurf-like molecular interaction fields, and (iii) 4D-Fingerprints. The best hERG SVM model achieved 89% accuracy and the three-best SVM models were able to screen a Pubchem data set with an accuracy of 97%. The best T. pyriformis model had an R(2) value of 0.924 for the training set and was able to predict the continuous end points for two test sets with R(2) values of 0.832 and 0.620, respectively. The studies presented within demonstrate the predictive ability (classification and continuous end points) of QSAR models constructed from curated data sets, biologically relevant molecular descriptors, and Support Vector Machines and Support Vector Regression. The ability of these protocols and methodologies to accommodate large data sets (several thousands compounds) that are chemically diverse - and in the case of classification modeling unbalanced (one experimental outcome dominates the data set) - allows scientists to further explore a remarkable amount of biological and chemical information.