Exploring the epigenetic drug discovery landscape.

INTRODUCTION: Epigenetic modification has been implicated in a wide range of diseases and the ability to modulate such systems is a lucrative therapeutic strategy in drug discovery. Areas covered: This article focuses on the concepts and drug discovery aspects of epigenomics. This is achieved by providing a survey of the following concepts: (i) factors influencing epigenetics, (ii) diseases arising from epigenetics, (iii) epigenetic enzymes as druggable targets along with coverage of existing FDA-approved drugs and pharmacological agents, and (iv) drug repurposing/repositioning as a means for rapid discovery of pharmacological agents targeting epigenetics. Expert opinion: Despite significant interests in targeting epigenetic modifiers as a therapeutic route, certain classes of target proteins are heavily studied while some are less characterized. Thus, such orphan target proteins are not yet druggable with limited report of active modulators. Current research points towards a great future with novel drugs directed to the many complex multifactorial diseases of humans, which are still often poorly understood and difficult to treat.

Towards Predicting the Cytochrome P450 Modulation: From QSAR to Proteochemometric Modeling.

Drug metabolism determines the fate of a drug when it enters the human body and is a critical factor in defining their absorption, distribution, metabolism, excretion and toxicity (ADMET) characteristics. Among the various drug metabolizing enzymes, cytochrome P450s (CYP450) constitute an important protein family that aside from functioning in xenobiotic metabolism, is also responsible for a diverse array of other roles encompassing steroid and cholesterol biosynthesis, fatty acid metabolism, calcium homeostasis, neuroendocrine functions and growth regulation. Although CYP450 typically converts xenobiotics into safe metabolites, there are some situations whereby the metabolite is more toxic than its parent molecule. Computational modeling has been instrumental in CYP450 research by rationalizing the nature of the binding event (i.e. inhibit or induce CYP450s) or metabolic stability of query compounds of interest. A plethora of computational approaches encompassing ligand, structure and systems based approaches have been utilized to model CYP450-ligand interactions. This review provides a brief background on the CYP450 family (i.e. its roles, advantages and disadvantages as well as its modulators) and then discusses the various computational approaches that have been used to model CYP450-ligand interaction. Particular focus was given to the use of quantitative structure-activity relationship (QSAR) and more recent proteochemometric modeling studies. Finally, a perspective on the current state of the art and future trends of the field is also provided.

Probing the origins of human acetylcholinesterase inhibition via QSAR modeling and molecular docking.

Alzheimer's disease (AD) is a chronic neurodegenerative disease which leads to the gradual loss of neuronal cells. Several hypotheses for AD exists (e.g., cholinergic, amyloid, tau hypotheses, etc.). As per the cholinergic hypothesis, the deficiency of choline is responsible for AD; therefore, the inhibition of AChE is a lucrative therapeutic strategy for the treatment of AD. Acetylcholinesterase (AChE) is an enzyme that catalyzes the breakdown of the neurotransmitter acetylcholine that is essential for cognition and memory. A large non-redundant data set of 2,570 compounds with reported IC50 values against AChE was obtained from ChEMBL and employed in quantitative structure-activity relationship (QSAR) study so as to gain insights on their origin of bioactivity. AChE inhibitors were described by a set of 12 fingerprint descriptors and predictive models were constructed from 100 different data splits using random forest. Generated models afforded R (2), [Formula: see text] and [Formula: see text] values in ranges of 0.66-0.93, 0.55-0.79 and 0.56-0.81 for the training set, 10-fold cross-validated set and external set, respectively. The best model built using the substructure count was selected according to the OECD guidelines and it afforded R (2), [Formula: see text] and [Formula: see text] values of 0.92 ± 0.01, 0.78 ± 0.06 and 0.78 ± 0.05, respectively. Furthermore, Y-scrambling was applied to evaluate the possibility of chance correlation of the predictive model. Subsequently, a thorough analysis of the substructure fingerprint count was conducted to provide informative insights on the inhibitory activity of AChE inhibitors. Moreover, Kennard-Stone sampling of the actives were applied to select 30 diverse compounds for further molecular docking studies in order to gain structural insights on the origin of AChE inhibition. Site-moiety mapping of compounds from the diversity set revealed three binding anchors encompassing both hydrogen bonding and van der Waals interaction. Molecular docking revealed that compounds 13, 5 and 28 exhibited the lowest binding energies of -12.2, -12.0 and -12.0 kcal/mol, respectively, against human AChE, which is modulated by hydrogen bonding, π-π stacking and hydrophobic interaction inside the binding pocket. These information may be used as guidelines for the design of novel and robust AChE inhibitors.

Large-scale ligand-based predictive modelling using support vector machines.

The increasing size of datasets in drug discovery makes it challenging to build robust and accurate predictive models within a reasonable amount of time. In order to investigate the effect of dataset sizes on predictive performance and modelling time, ligand-based regression models were trained on open datasets of varying sizes of up to 1.2 million chemical structures. For modelling, two implementations of support vector machines (SVM) were used. Chemical structures were described by the signatures molecular descriptor. Results showed that for the larger datasets, the LIBLINEAR SVM implementation performed on par with the well-established libsvm with a radial basis function kernel, but with dramatically less time for model building even on modest computer resources. Using a non-linear kernel proved to be infeasible for large data sizes, even with substantial computational resources on a computer cluster. To deploy the resulting models, we extended the Bioclipse decision support framework to support models from LIBLINEAR and made our models of logD and solubility available from within Bioclipse.

Roles of d-Amino Acids on the Bioactivity of Host Defense Peptides.

Host defense peptides (HDPs) are positively-charged and amphipathic components of the innate immune system that have demonstrated great potential to become the next generation of broad spectrum therapeutic agents effective against a vast array of pathogens and tumor. As such, many approaches have been taken to improve the therapeutic efficacy of HDPs. Amongst these methods, the incorporation of d-amino acids (d-AA) is an approach that has demonstrated consistent success in improving HDPs. Although, virtually all HDP review articles briefly mentioned about the role of d-AA, however it is rather surprising that no systematic review specifically dedicated to this topic exists. Given the impact that d-AA incorporation has on HDPs, this review aims to fill that void with a systematic discussion of the impact of d-AA on HDPs.

New Polyfunctional Phragmalin Limonoids from Neobeguea mahafalensis.

Neobeguea mahafalensis is used as a medicinal plant in Madagascar. A decoction of the stem bark of this species is reported to treat back pain. Recently, it was reported that a decoction of the root bark, containing two novel phragmalin limonoids as identified active constituents, exhibited an extraordinarily high potency and remarkably long duration in augmenting sexual activity in male rodents.From the dichloromethane extract of the root barks of N. mahafalensis, nine phragmalin limonoids were isolated, of which eight were novel compounds. The structures were established mainly by extensive use of 2D NMR spectroscopic techniques and high-resolution mass spectrometry. One of the new compounds named dodoguin displayed sleep-inducing activity in Swiss albino mice. The amount of 3-10 mg/kg of this compound induced sleep 18-22 min after its administration with a duration of 16-18 min.

Origin of aromatase inhibitory activity via proteochemometric modeling.

Aromatase, the rate-limiting enzyme that catalyzes the conversion of androgen to estrogen, plays an essential role in the development of estrogen-dependent breast cancer. Side effects due to aromatase inhibitors (AIs) necessitate the pursuit of novel inhibitor candidates with high selectivity, lower toxicity and increased potency. Designing a novel therapeutic agent against aromatase could be achieved computationally by means of ligand-based and structure-based methods. For over a decade, we have utilized both approaches to design potential AIs for which quantitative structure-activity relationships and molecular docking were used to explore inhibitory mechanisms of AIs towards aromatase. However, such approaches do not consider the effects that aromatase variants have on different AIs. In this study, proteochemometrics modeling was applied to analyze the interaction space between AIs and aromatase variants as a function of their substructural and amino acid features. Good predictive performance was achieved, as rigorously verified by 10-fold cross-validation, external validation, leave-one-compound-out cross-validation, leave-one-protein-out cross-validation and Y-scrambling tests. The investigations presented herein provide important insights into the mechanisms of aromatase inhibitory activity that could aid in the design of novel potent AIs as breast cancer therapeutic agents.

osFP: a web server for predicting the oligomeric states of fluorescent proteins.

BACKGROUND: Currently, monomeric fluorescent proteins (FP) are ideal markers for protein tagging. The prediction of oligomeric states is helpful for enhancing live biomedical imaging. Computational prediction of FP oligomeric states can accelerate the effort of protein engineering efforts of creating monomeric FPs. To the best of our knowledge, this study represents the first computational model for predicting and analyzing FP oligomerization directly from the amino acid sequence. RESULTS: After data curation, an exhaustive data set consisting of 397 non-redundant FP oligomeric states was compiled from the literature. Results from benchmarking of the protein descriptors revealed that the model built with amino acid composition descriptors was the top performing model with accuracy, sensitivity and specificity in excess of 80% and MCC greater than 0.6 for all three data subsets (e.g. training, tenfold cross-validation and external sets). The model provided insights on the important residues governing the oligomerization of FP. To maximize the benefit of the generated predictive model, it was implemented as a web server under the R programming environment. CONCLUSION: osFP affords a user-friendly interface that can be used to predict the oligomeric state of FP using the protein sequence. The advantage of osFP is that it is platform-independent meaning that it can be accessed via a web browser on any operating system and device. osFP is freely accessible at http://codes.bio/osfp/ while the source code and data set is provided on GitHub at https://github.com/chaninn/osFP/.Graphical Abstract.

Proteochemometric model for predicting the inhibition of penicillin-binding proteins.

Neisseria gonorrhoeae infection threatens to become an untreatable sexually transmitted disease in the near future owing to the increasing emergence of N. gonorrhoeae strains with reduced susceptibility and resistance to the extended-spectrum cephalosporins (ESCs), i.e. ceftriaxone and cefixime, which are the last remaining option for first-line treatment of gonorrhea. Alteration of the penA gene, encoding penicillin-binding protein 2 (PBP2), is the main mechanism conferring penicillin resistance including reduced susceptibility and resistance to ESCs. To predict and investigate putative amino acid mutations causing ß-lactam resistance particularly for ESCs, we applied proteochemometric modeling to generalize N. gonorrhoeae susceptibility data for predicting the interaction of PBP2 with therapeutic ß-lactam antibiotics. This was afforded by correlating publicly available data on antimicrobial susceptibility of wild-type and mutant N. gonorrhoeae strains for penicillin-G, cefixime and ceftriaxone with 50 PBP2 protein sequence data using partial least-squares projections to latent structures. The generated model revealed excellent predictability (R2=0.91, Q2=0.77, QExt2=0.78). Moreover, our model identified amino acid mutations in PBP2 with the highest impact on antimicrobial susceptibility and provided information on physicochemical properties of amino acid mutations affecting antimicrobial susceptibility. Our model thus provided insight into the physicochemical basis for resistance development in PBP2 suggesting its use for predicting and monitoring novel PBP2 mutations that may emerge in the future.

Benchmarking study of parameter variation when using signature fingerprints together with support vector machines.

QSAR modeling using molecular signatures and support vector machines with a radial basis function is increasingly used for virtual screening in the drug discovery field. This method has three free parameters: C, γ, and signature height. C is a penalty parameter that limits overfitting, γ controls the width of the radial basis function kernel, and the signature height determines how much of the molecule is described by each atom signature. Determination of optimal values for these parameters is time-consuming. Good default values could therefore save considerable computational cost. The goal of this project was to investigate whether such default values could be found by using seven public QSAR data sets spanning a wide range of end points and using both a bit version and a count version of the molecular signatures. On the basis of the experiments performed, we recommend a parameter set of heights 0 to 2 for the count version of the signature fingerprints and heights 0 to 3 for the bit version. These are in combination with a support vector machine using C in the range of 1 to 100 and γ in the range of 0.001 to 0.1. When data sets are small or longer run times are not a problem, then there is reason to consider the addition of height 3 to the count fingerprint and a wider grid search. However, marked improvements should not be expected.

Ligand-based target prediction with signature fingerprints.

When evaluating a potential drug candidate it is desirable to predict target interactions in silico prior to synthesis in order to assess, e.g., secondary pharmacology. This can be done by looking at known target binding profiles of similar compounds using chemical similarity searching. The purpose of this study was to construct and evaluate the performance of chemical fingerprints based on the molecular signature descriptor for performing target binding predictions. For the comparison we used the area under the receiver operating characteristics curve (AUC) complemented with net reclassification improvement (NRI). We created two open source signature fingerprints, a bit and a count version, and evaluated their performance compared to a set of established fingerprints with regards to predictions of binding targets using Tanimoto-based similarity searching on publicly available data sets extracted from ChEMBL. The results showed that the count version of the signature fingerprint performed on par with well-established fingerprints such as ECFP. The count version outperformed the bit version slightly; however, the count version is more complex and takes more computing time and memory to run so its usage should probably be evaluated on a case-by-case basis. The NRI based tests complemented the AUC based ones and showed signs of higher power.

Illuminating the origins of spectral properties of green fluorescent proteins via proteochemometric and molecular modeling.

Green fluorescent protein (GFP) has immense utility in biomedical imaging owing to its autofluorescent nature. In efforts to broaden the spectral diversity of GFP, there have been several reports of engineered mutants via rational design and random mutagenesis. Understanding the origins of spectral properties of GFP could be achieved by means of investigating its structure-activity relationship. The first quantitative structure-property relationship study for modeling the spectral properties, particularly the excitation and emission maximas, of GFP was previously proposed by us some years ago in which quantum chemical descriptors were used for model development. However, such simplified model does not consider possible effects that neighboring amino acids have on the conjugated π-system of GFP chromophore. This study describes the development of a unified proteochemometric model in which the GFP chromophore and amino acids in its vicinity are both considered in the same model. The predictive performance of the model was verified by internal and external validation as well as Y-scrambling. Our strategy provides a general solution for elucidating the contribution that specific ligand and protein descriptors have on the investigated spectral property, which may be useful in engineering novel GFP variants with desired characteristics.

Semisynthesis of libiguin A and its analogues by trans-lactonization of phragmalin.

Libiguins are limonoids with highly potent sexual activity enhancing effects, originally isolated from the Madagascarian Meliaceae species Neobeguea mahafalensis, where they exist in only minute quantities. Their low natural abundance has hampered mapping of their biological effects. Here we describe an approach to the semisynthesis of libiguin A and its close analogues 1-3 starting from phragmalin, which is a limonoid present in high amounts in a commercially cultivated Meliaceae species, Chukrasia tabularis, allowing the preparation of libiguins in appreciable quantities.

Libiguins A and B: novel phragmalin limonoids isolated from Neobeguea mahafalensis causing profound enhancement of sexual activity.

In a screening programme directed towards the discovery of drugs that could enhance sexual activity, we found that a decoction of the root bark of Neobeguea mahafalensis displayed an extraordinarily high potency and remarkably long duration in augmenting sexual activity in male rodents. Bioassay-guided fractionation led to the isolation of two pharmacoactive constituents, which turned out to be novel 1,8,9-orthoacetate phragmalin limonoids that we named libiguins A and B, each with a C-16/30 δ-lactone ring. Chemical structures were established by the interpretation of their 1D and 2D NMR data. In vivo pharmacological tests showed that starting with a treatment from 0.004-0.4 mg/kg/day for three consecutive days, over a 3-h sampling period, these limonoids induced a long-lasting augmentation of frequency and sustainment of mounting behaviour in male rodents, with an effect lasting for up to 11 days post-treatment. Libiguin A proved to be markedly more potent than libiguin B. This report is the first of limonoids having such an effect, and the findings could lead to novel therapies for the treatment of sexual dysfunction. Moreover, the results can serve as an opening to elucidate the central physiological control of mating behaviour, which is still not well mapped out.

A unified proteochemometric model for prediction of inhibition of cytochrome p450 isoforms.

A unified proteochemometric (PCM) model for the prediction of the ability of drug-like chemicals to inhibit five major drug metabolizing CYP isoforms (i.e. CYP1A2, CYP2C9, CYP2C19, CYP2D6 and CYP3A4) was created and made publicly available under the Bioclipse Decision Support open source system at www.cyp450model.org. In regards to the proteochemometric modeling we represented the chemical compounds by molecular signature descriptors and the CYP-isoforms by alignment-independent description of composition and transition of amino acid properties of their protein primary sequences. The entire training dataset contained 63 391 interactions and the best PCM model was obtained using signature descriptors of height 1, 2 and 3 and inducing the model with a support vector machine. The model showed excellent predictive ability with internal AUC = 0.923 and an external AUC = 0.940, as evaluated on a large external dataset. The advantage of PCM models is their extensibility making it possible to extend our model for new CYP isoforms and polymorphic CYP forms. A key benefit of PCM is that all proteins are confined in one single model, which makes it generally more stable and predictive as compared with single target models. The inclusion of the model in Bioclipse Decision Support makes it possible to make virtual instantaneous predictions (∼100 ms per prediction) while interactively drawing or modifying chemical structures in the Bioclipse chemical structure editor.

Design and evaluation of substrate-based octapeptide and non substrate-based tetrapeptide inhibitors of dengue virus NS2B-NS3 proteases.

A series of 45 peptide inhibitors was designed, synthesized, and evaluated against the NS2B-NS3 proteases of the four subtypes of dengue virus, DEN-1-4. The design was based on proteochemometric models for Michaelis (Km) and cleavage rate constants (kcat) of protease substrates. This led first to octapeptides showing submicromolar or low micromolar inhibitory activities on the four proteases. Stepwise removal of cationic substrate non-prime side residues and variations in the prime side sequence resulted finally in an uncharged tetrapeptide, WYCW-NH2, with inhibitory Ki values of 4.2, 4.8, 24.4, and 11.2 µM for the DEN-1-4 proteases, respectively. Analysis of the inhibition data by proteochemometric modeling suggested the possibility for different binding poses of the shortened peptides compared to the octapeptides, which was supported by results of docking of WYCW-NH2 into the X-ray structure of DEN-3 protease.

Bioclipse-R: integrating management and visualization of life science data with statistical analysis.

SUMMARY: Bioclipse, a graphical workbench for the life sciences, provides functionality for managing and visualizing life science data. We introduce Bioclipse-R, which integrates Bioclipse and the statistical programming language R. The synergy between Bioclipse and R is demonstrated by the construction of a decision support system for anticancer drug screening and mutagenicity prediction, which shows how Bioclipse-R can be used to perform complex tasks from within a single software system. AVAILABILITY AND IMPLEMENTATION: Bioclipse-R is implemented as a set of Java plug-ins for Bioclipse based on the R-package rj. Source code and binary packages are available from https://github.com/bioclipse and http://www.bioclipse.net/bioclipse-r, respectively. CONTACT: martin.eklund@farmbio.uu.se SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Enzymatic analysis of recombinant Japanese encephalitis virus NS2B(H)-NS3pro protease with fluorogenic model peptide substrates.

BACKGROUND: Japanese encephalitis virus (JEV), a member of the Flaviviridae family, causes around 68,000 encephalitis cases annually, of which 20-30% are fatal, while 30-50% of the recovered cases develop severe neurological sequelae. Specific antivirals for JEV would be of great importance, particularly in those cases where the infection has become persistent. Being indispensable for flaviviral replication, the NS2B-NS3 protease is a promising target for design of anti-flaviviral inhibitors. Contrary to related flaviviral proteases, the JEV NS2B-NS3 protease is structurally and mechanistically much less characterized. Here we aimed at establishing a straightforward procedure for cloning, expression, purification and biochemical characterization of JEV NS2B(H)-NS3pro protease. METHODOLOGY/PRINCIPAL FINDINGS: The full-length sequence of JEV NS2B-NS3 genotype III strain JaOArS 982 was obtained as a synthetic gene. The sequence of NS2B(H)-NS3pro was generated by splicing by overlap extension PCR (SOE-PCR) and cloned into the pTrcHisA vector. Hexahistidine-tagged NS2B(H)-NS3pro, expressed in E. coli as soluble protein, was purified to >95% purity by a single-step immobilized metal affinity chromatography. SDS-PAGE and immunoblotting of the purified enzyme demonstrated NS2B(H)-NS3pro precursor and its autocleavage products, NS3pro and NS2B(H), as 36, 21, and 10 kDa bands, respectively. Kinetic parameters, K(m) and k(cat), for fluorogenic protease model substrates, Boc-GRR-amc, Boc-LRR-amc, Ac-nKRR-amc, Bz-nKRR-amc, Pyr-RTKR-amc and Abz-(R)(4)SAG-nY-amide, were obtained using inner filter effect correction. The highest catalytic efficiency k(cat)/K(m) was found for Pyr-RTKR-amc (k(cat)/K(m): 1962.96 ± 85.0 M(-1) s(-1)) and the lowest for Boc-LRR-amc (k(cat)/K(m): 3.74±0.3 M(-1) s(-1)). JEV NS3pro is inhibited by aprotinin but to a lesser extent than DEN and WNV NS3pro. CONCLUSIONS/SIGNIFICANCE: A simplified procedure for the cloning, overexpression and purification of the NS2B(H)-NS3pro was established which is generally applicable to other flaviviral proteases. Kinetic parameters obtained for a number of model substrates and inhibitors, are useful for the characterization of substrate specificity and eventually for the design of high-throughput assays aimed at antiviral inhibitor discovery.

Brunn: an open source laboratory information system for microplates with a graphical plate layout design process.

BACKGROUND: Compound profiling and drug screening generates large amounts of data and is generally based on microplate assays. Current information systems used for handling this are mainly commercial, closed source, expensive, and heavyweight and there is a need for a flexible lightweight open system for handling plate design, and validation and preparation of data. RESULTS: A Bioclipse plugin consisting of a client part and a relational database was constructed. A multiple-step plate layout point-and-click interface was implemented inside Bioclipse. The system contains a data validation step, where outliers can be removed, and finally a plate report with all relevant calculated data, including dose-response curves. CONCLUSIONS: Brunn is capable of handling the data from microplate assays. It can create dose-response curves and calculate IC50 values. Using a system of this sort facilitates work in the laboratory. Being able to reuse already constructed plates and plate layouts by starting out from an earlier step in the plate layout design process saves time and cuts down on error sources.

Services for prediction of drug susceptibility for HIV proteases and reverse transcriptases at the HIV drug research centre.

SUMMARY: The HIV Drug Research Centre (HIVDRC) has established Web services for prediction of drug susceptibility for HIV proteases and reverse transcriptases. The services are based on two proteochemometric models which accepts a protease or reverse transcriptase sequence in amino acid form, and outputs the predicted drug susceptibility values. The predictions are based on a comprehensive analysis where all the relevant inhibitors are included, resulting in models with excellent predictive capabilities. AVAILABILITY AND IMPLEMENTATION: The services are implemented as interoperable Web services (REST and XMPP), with supporting web pages to allow for individual analyses. A set of plugins were also developed which make the services available from the Bioclipse workbench for life science. Services are available at http://www.hivdrc.org/services.