Desirability-based methods of multiobjective optimization and ranking for global QSAR studies. Filtering safe and potent drug candidates from combinatorial libraries

Immune-Mediated Drug-Induced Liver Injury: Immunogenetics and Experimental Models

Drug-induced liver injury (DILI) is a challenging clinical event in medicine, particularly because of its ability to present with a variety of phenotypes including that of autoimmune hepatitis or other immune mediated liver injuries. Limited diagnostic and therapeutic tools are available, mostly because its pathogenesis has remained poorly understood for decades. The recent scientific and technological advancements in genomics and immunology are paving the way for a better understanding of the molecular aspects of DILI. This review provides an updated overview of the genetic predisposition and immunological mechanisms behind the pathogenesis of DILI and presents the state-of-the-art experimental models to study DILI at the pre-clinical level.

Multi-objective optimization methods in novel drug design

Introduction: In multi-objective drug design, optimization gains importance, being upgraded to a discipline that attracts its own research. Current strategies are broadly classified into single - objective optimization (SOO) and multi-objective optimization (MOO).Areas covered: Starting with SOO and the ways used to incorporate multiple criteria into it, the present review focuses on MOO techniques, their comparison, advantages, and restrictions. Pareto analysis and the concept of dominance stand in the core of MOO. The Pareto front, Pareto ranking, and limitations of Pareto-based methods, due to high dimensions and data uncertainty, are outlined. Desirability functions and the weighted sum approaches are described as stand-alone techniques to transform the MOO problem to SOO or in combination with pareto analysis and evolutionary algorithms. Representative applications in different drug research areas are also discussed.Expert opinion: Despite their limitations, the use of combined MOO techniques, as well as being complementary to SOO or in conjunction with artificial intelligence, contributes dramatically to efficient drug design, assisting decisions and increasing success probabilities. For multi-target drug design, optimization is supported by network approaches, while applicability of MOO to other fields like drug technology or biological complexity opens new perspectives in the interrelated fields of medicinal chemistry and molecular biology.

Uncertainty quantification in drug design

Machine learning and artificial intelligence are increasingly being applied to the drug-design process as a result of the development of novel algorithms, growing access, the falling cost of computation and the development of novel technologies for generating chemically and biologically relevant data. There has been recent progress in fields such as molecular de novo generation, synthetic route prediction and, to some extent, property predictions. Despite this, most research in these fields has focused on improving the accuracy of the technologies, rather than on quantifying the uncertainty in the predictions. Uncertainty quantification will become a key component in autonomous decision making and will be crucial for integrating machine learning and chemistry automation to create an autonomous design-make-test-analyse cycle. This review covers the empirical, frequentist and Bayesian approaches to uncertainty quantification, and outlines how they can be used for drug design. We also outline the impact of uncertainty quantification on decision making.

Comparative study between deep learning and QSAR classifications for TNBC inhibitors and novel GPCR agonist discovery

Machine learning is a well-known approach for virtual screening. Recently, deep learning, a machine learning algorithm in artificial neural networks, has been applied to the advancement of precision medicine and drug discovery. In this study, we performed comparative studies between deep neural networks (DNN) and other ligand-based virtual screening (LBVS) methods to demonstrate that DNN and random forest (RF) were superior in hit prediction efficiency. By using DNN, several triple-negative breast cancer (TNBC) inhibitors were identified as potent hits from a screening of an in-house database of 165,000 compounds. In broadening the application of this method, we harnessed the predictive properties of trained model in the discovery of G protein-coupled receptor (GPCR) agonist, by which computational structure-based design of molecules could be greatly hindered by lack of structural information. Notably, a potent (~ 500 nM) mu-opioid receptor (MOR) agonist was identified as a hit from a small-size training set of 63 compounds. Our results show that DNN could be an efficient module in hit prediction and provide experimental evidence that machine learning could identify potent hits in silico from a limited training set.

Computer-Aided Design of Antimicrobial Peptides: Are We Generating Effective Drug Candidates?

Antimicrobial peptides (AMPs), especially antibacterial peptides, have been widely investigated as potential alternatives to antibiotic-based therapies. Indeed, naturally occurring and synthetic AMPs have shown promising results against a series of clinically relevant bacteria. Even so, this class of antimicrobials has continuously failed clinical trials at some point, highlighting the importance of AMP optimization. In this context, the computer-aided design of AMPs has put together crucial information on chemical parameters and bioactivities in AMP sequences, thus providing modes of prediction to evaluate the antibacterial potential of a candidate sequence before synthesis. Quantitative structure-activity relationship (QSAR) computational models, for instance, have greatly contributed to AMP sequence optimization aimed at improved biological activities. In addition to machine-learning methods, the de novo design, linguistic model, pattern insertion methods, and genetic algorithms, have shown the potential to boost the automated design of AMPs. However, how successful have these approaches been in generating effective antibacterial drug candidates? Bearing this in mind, this review will focus on the main computational strategies that have generated AMPs with promising activities against pathogenic bacteria, as well as anti-infective potential in different animal models, including sepsis and cutaneous infections. Moreover, we will point out recent studies on the computer-aided design of antibiofilm peptides. As expected from automated design strategies, diverse candidate sequences with different structural arrangements have been generated and deposited in databases. We will, therefore, also discuss the structural diversity that has been engendered.

Does Size Really Matter? Probing the Efficacy of Structural Reduction in the Optimization of Bioderived Compounds - A Computational "Proof-of-Concept"

Over the years, numerous synthetic approaches have been utilized in drug design to improve the pharmacological properties of naturally derived compounds and most importantly, minimize toxic effects associated with their transition to drugs. The reduction of complex bioderived compounds to simpler bioactive fragments has been identified as a viable strategy to develop lead compounds with improved activities and minimal toxicities. Although this ‘reductive’ strategy has been widely exemplified, underlying biological events remain unresolved, hence the unanswered question remains how does the fragmentation of a natural compound improve its bioactivity and reduce toxicities? Herein, using a combinatorial approach, we initialize a computational “proof-of- concept” to expound the differential pharmacological and antagonistic activities of a natural compound, Anguinomycin D, and its synthetic fragment, SB640 towards Exportin Chromosome Region Maintenance 1 (CRM1). Interestingly, our findings revealed that in comparison with the parent compound, SB640 exhibited improved pharmacological attributes, while toxicities and off-target activities were relatively minimal. Moreover, we observed that the reduced size of SB640 allowed ‘deep access’ at the Nuclear Export Signals (NES) binding groove of CRM1, which favored optimal and proximal positioning towards crucial residues while the presence of the long polyketide tail in Anguinomycin D constrained its burial at the hydrophobic groove. Furthermore, with regards to their antagonistic functions, structural inactivation (rigidity) was more pronounced in CRM1 when bound by SB640 as compared to Anguinomycin D. These findings provide essential insights that portray synthetic fragmentation of natural compounds as a feasible approach towards the discovery of potential leads in disease treatment.

Influence of Descriptor Implementation on Compound Ranking Based on Multiparameter Assessment

Most of the common molecular descriptors have numerous different implementations. This can influence the results of compound prioritization based on the multiparameter assessment (MPA) approach that allows a medicinal chemist to simultaneously analyze and achieve the desired balance of the diverse and often conflicting molecular and pharmacological properties. In this study, we analyzed the feasibility of using different implementations of common descriptors (logP, logS, TPSA, logBB, hERG, nHBA) interchangeably in predesigned sets of requirements in the course of multiparameter compound optimization. The influence of methods of descriptor calculation, continuity or discreteness of their values, their applicability domains, as well as of the nature of desirability functions in an MPA profile were examined in terms of the stability of MPA compound ranking. It was shown that the interchangeable use of different methods of descriptor calculation is reliably acceptable only for continuously distributed parameters transformed by a smooth desirability function. If a descriptor in an MPA scheme is discretely distributed, only the implementation that was used for building the scoring profile may be used for assessment. An inconsistency of assessment due to different applicability domains of descriptors was also demonstrated.

From flamingo dance to (desirable) drug discovery: a nature-inspired approach

The therapeutic effects of drugs are well known to result from their interaction with multiple intracellular targets. Accordingly, the pharma industry is currently moving from a reductionist approach based on a 'one-target fixation' to a holistic multitarget approach. However, many drug discovery practices are still procedural abstractions resulting from the attempt to understand and address the action of biologically active compounds while preventing adverse effects. Here, we discuss how drug discovery can benefit from the principles of evolutionary biology and report two real-life case studies. We do so by focusing on the desirability principle, and its many features and applications, such as machine learning-based multicriteria virtual screening.

Integrated in Silico Methods for the Design and Optimization of Novel Drug Candidates

Although almost fully automated, the discovery of novel, effective, and safe drugs is still a long-term and highly expensive process. Consequently, the need for fleet, rational, and cost-efficient development of novel drugs is crucial, and nowadays the advanced in silico drug design methodologies seem to effectively meet these issues. The aim of this chapter is to provide a comprehensive overview of some of the current trends and advances in the in silico design of novel drug candidates with a special emphasis on 6-fluoroquinolone (6-FQ) antibacterials as potential novel Mycobacterium tuberculosis DNA gyrase inhibitors. In particular, the chapter covers some of the recent aspects of a wide range of in silico drug discovery approaches including multidimensional machine-learning methods, ligand-based and structure-based methodologies, as well as their proficient combination and integration into an intelligent virtual screening protocol for design and optimization of novel 6-FQ analogs.

Multi-objective optimization methods in drug design

Drug discovery is a challenging multi-objective problem where numerous pharmaceutically important objectives need to be adequately satisfied for a solution to be found. The problem is characterized by vast, complex solution spaces further perplexed by the presence of conflicting objectives. Multi-objective optimization methods, designed specifically to address such problems, have been introduced to the drug discovery field over a decade ago and have steadily gained in acceptance ever since. This paper reviews the latest multi-objective methods and applications reported in the literature, specifically in quantitative structure–activity modeling, docking, de novo design and library design. Further, the paper reports on related developments in drug discovery research and advances in the multi-objective optimization field.

Quantifying the chemical beauty of drugs

Drug-likeness is a key consideration when selecting compounds during the early stages of drug discovery. However, evaluation of drug-likeness in absolute terms does not reflect adequately the whole spectrum of compound quality. More worryingly, widely used rules may inadvertently foster undesirable molecular property inflation as they permit the encroachment of rule-compliant compounds towards their boundaries. We propose a measure of drug-likeness based on the concept of desirability called the quantitative estimate of drug-likeness (QED). The empirical rationale of QED reflects the underlying distribution of molecular properties. QED is intuitive, transparent, straightforward to implement in many practical settings and allows compounds to be ranked by their relative merit. We extended the utility of QED by applying it to the problem of molecular target druggability assessment by prioritizing a large set of published bioactive compounds. The measure may also capture the abstract notion of aesthetics in medicinal chemistry.

Functional polymer thin films designed for antifouling materials and biosensors

Polymer thin films offer a versatile and ubiquitous platform for a wide variety of real-world applications in biomedicine, nanotechnology, catalysis, photovoltaic devices, and energy conversion and storage. Depending on the chemical composition of the polymers and the associated microenvironment, the physicochemical properties (biocompatibility, stability, wettability, adhesion, morphology, surface free energy, and others) of polymer films can be tuned for a specific application through precisely controlled surface synthesis and the incorporation of desirable and responsive functional groups. In this short review, we first summarise the methods most commonly used for the fabrication of polymer thin films. Then we discuss how these polymer thin films can be used in a selection of biomedical applications in antifouling materials and biosensors. Some directions for the rational design of polymer thin films to achieve a specific function or application are also provided.

How desirable are your IC50s? A way to enhance screening-based decision making

Dose-response curves, resulting in estimates of endpoints such as the IC(50), are fundamental to drug discovery. However, some estimates are more reliable than others. It is important to know just how reliable an estimate is if we want to base decisions on it or use it in further modeling. In this study, the authors propose a new measure of endpoint reliability, based on the concept of desirability first introduced by Harrington. The solution is not dependent on the application used to analyze the experimental data, provided a number of parameters to characterize the dose-response curve are available. The authors show how this score can be used as an objective and consistent measure to rank screening results, combine information from groups of experiments, and determine optimal levels of characterization of a compound's biological activity.

Prioritizing Hits with Appropriate Trade-Offs Between HIV-1 Reverse Transcriptase Inhibitory Efficacy and MT4 Blood Cells Toxicity Through Desirability-Based Multiobjective Optimization and Ranking

Nonnucleoside reverse transcriptase (RT) inhibitors (NNRTIs) constitute a promising therapeutic option for AIDS. However, the emergence of virus-NNRTIs resistance was found to be a major problem in the field. Toward that goal, a "knock-out" strategy stands out between the several options to circumvent the problem. However the high drug or drug-drug concentrations often used generate additional safety concerns. The need for approaches able to early integrate drug- or lead-likeness, toxicity and bioavailability criteria in the drug discovery phase is an emergent issue. Given that, we propose a combined strategy based on desirability-based multiobjective optimization (MOOP) and ranking for the prioritization of HIV-1 NNRTIs hits with appropriate trade-offs between inhibitory efficacy over the HIV-1 RT and toxic effects over MT4 blood cells. Through the MOOP process, the theoretical levels of the predictive variables required to reach a desirable RT inhibitor candidate with the best possible compromise between efficacy and safety were found. This information is used as a pattern to rank the library of compounds according to a similarity-based structural criterion, providing a ranking quality of 64 %/71 %/73 % in training/validation/test set. A comparative study between the sequential, parallel and multiobjective virtual screening revealed that the multiobjective approach can outperform the other approaches. These results suggest that the identification of NNRTIs hits with appropriate trade-offs between potency and safety, rather than fully optimized hits solely based on potency, can facilitate the hit to lead transition and increase the likelihood of the candidate to evolve into a successful antiretroviral drug.

Multi-target spectral moment QSAR versus ANN for antiparasitic drugs against different parasite species

There are many of pathogen parasite species with different susceptibility profile to antiparasitic drugs. Unfortunately, almost QSAR models predict the biological activity of drugs against only one parasite species. Consequently, predicting the probability with which a drug is active against different species with a single unify model is a goal of the major importance. In so doing, we use Markov Chains theory to calculate new multi-target spectral moments to fit a QSAR model that predict by the first time a mt-QSAR model for 500 drugs tested in the literature against 16 parasite species and other 207 drugs no tested in the literature using spectral moments. The data was processed by linear discriminant analysis (LDA) classifying drugs as active or non-active against the different tested parasite species. The model correctly classifies 311 out of 358 active compounds (86.9%) and 2328 out of 2577 non-active compounds (90.3%) in training series. Overall training performance was 89.9%. Validation of the model was carried out by means of external predicting series. In these series the model classified correctly 157 out 190, 82.6% of antiparasitic compounds and 1151 out of 1277 non-active compounds (90.1%). Overall predictability performance was 89.2%. In addition we developed four types of non Linear Artificial neural networks (ANN) and we compared with the mt-QSAR model. The improved ANN model had an overall training performance was 87%. The present work report the first attempts to calculate within a unify framework probabilities of antiparasitic action of drugs against different parasite species based on spectral moment analysis.

QSAR models to predict mutagenicity of acrylates, methacrylates and alpha,beta-unsaturated carbonyl compounds

OBJECTIVE

The purpose of this study is to develop a quantitative structure-activity relationship (QSAR) model that can distinguish mutagenic from non-mutagenic species with alpha,beta-unsaturated carbonyl moiety using two endpoints for this activity - Ames test and mammalian cell gene mutation test - and also to gather information about the molecular features that most contribute to eliminate the mutagenic effects of these chemicals.

METHODS

Two data sets were used for modeling the two mutagenicity endpoints: (1) Ames test and (2) mammalian cells mutagenesis. The first one comprised 220 molecules, while the second one 48 substances, ranging from acrylates, methacrylates to alpha,beta-unsaturated carbonyl compounds. The QSAR models were developed by applying linear discriminant analysis (LDA) along with different sets of descriptors computed using the DRAGON software.

RESULTS

For both endpoints, there was a concordance of 89% in the prediction and 97% confidentiality by combining the three models for the Ames test mutagenicity. We have also identified several structural alerts to assist the design of new monomers.

SIGNIFICANCE

These individual models and especially their combination are attractive from the point of view of molecular modeling and could be used for the prediction and design of new monomers that do not pose a human health risk.

Multi-target spectral moment: QSAR for antiviral drugs vs. different viral species

The antiviral QSAR models have an important limitation today. They predict the biological activity of drugs against only one viral species. This is determined by the fact that most of the current reported molecular descriptors encode only information about the molecular structure. As a result, predicting the probability with which a drug is active against different viral species with a single unifying model is a goal of major importance. In this work, we use Markov Chain theory to calculate new multi-target spectral moments to fit a QSAR model for drugs active against 40 viral species. The model is based on 500 drugs (including active and non-active compounds) tested as antiviral agents in the recent literature; not all drugs were predicted against all viruses, but only those with experimental values. The database also contains 207 well-known compounds (not as recent as the previous ones) reported in the Merck Index with other activities that do not include antiviral action against any virus species. We used Linear Discriminant Analysis (LDA) to classify all these drugs into two classes as active or non-active against the different viral species tested, whose data we processed. The model correctly classifies 5129 out of 5594 non-active compounds (91.69%) and 412 out of 422 active compounds (97.63%). Overall training predictability was 92.34%. The validation of the model was carried out by means of external predicting series, the model classifying, thus, 2568 out of 2779 non-active compounds and 224 out of 229 active compounds. Overall training predictability was 92.82%. The present work reports the first attempts to calculate within a unified framework the probabilities of antiviral drugs against different virus species based on a spectral moment analysis.

Alignment-free prediction of a drug-target complex network based on parameters of drug connectivity and protein sequence of receptors

There are many drugs described with very different affinity to a large number of receptors. In this work, we selected drug-receptor pairs (DRPs) of affinity/nonaffinity drugs to similar/dissimilar receptors and we represented them as a large network, which may be used to identify drugs that can act on a receptor. Computational chemistry prediction of the biological activity based on quantitative structure-activity relationships (QSAR) substantially increases the potentialities of this kind of networks avoiding time- and resource-consuming experiments. Unfortunately, most QSAR models are unspecific or predict activity against only one receptor. To solve this problem, we developed here a multitarget QSAR (mt-QSAR) classification model. Overall model classification accuracy was 72.25% (1390/1924 compounds) in training, 72.28% (459/635) in cross-validation. Outputs of this mt-QSAR model were used as inputs to construct a network. The observed network has 1735 nodes (DRPs), 1754 edges or pairs of DRPs with similar drug-target affinity (sPDRPs), and low coverage density d = 0.12%. The predicted network has 1735 DRPs, 1857 sPDRPs, and also low coverage density d = 0.12%. After an edge-to-edge comparison (chi-square = 9420.3; p < 0.005), we have demonstrated that the predicted network is significantly similar to the one observed and both have a distribution closer to exponential than to normal.

Global antifungal profile optimization of chlorophenyl derivatives against Botrytis cinerea and Colletotrichum gloeosporioides

Twenty-two aromatic derivatives bearing a chlorine atom and a different chain in the para or meta position were prepared and evaluated for their in vitro antifungal activity against the phytopathogenic fungi Botrytis cinerea and Colletotrichum gloeosporioides. The results showed that maximum inhibition of the growth of these fungi was exhibited for enantiomers S and R of 1-(4'-chlorophenyl)-2-phenylethanol (3 and 4). Furthermore, their antifungal activity showed a clear structure-activity relationship (SAR) trend confirming the importance of the benzyl hydroxyl group in the inhibitory mechanism of the compounds studied. Additionally, a multiobjective optimization study of the global antifungal profile of chlorophenyl derivatives was conducted in order to establish a rational strategy for the filtering of new fungicide candidates from combinatorial libraries. The MOOP-DESIRE methodology was used for this purpose providing reliable ranking models that can be used later.