QSAR modeling of endpoints for peptides which is based on representation of the molecular structure by a sequence of amino acids

QSAR: Using the Past to Study the Present

Quantitative structure-activity relationships (QSAR) is a method for predicting the physical and biological properties of small molecules; it is in use in industry and public services. However, as any scientific method, it is challenged by more and more requests, especially considering its possible role in assessing the safety of new chemicals. To answer the question whether QSAR, by exploiting available knowledge, can build new knowledge, the chapter reviews QSAR methods in search of a QSAR epistemology. QSAR stands on tree pillars, i.e., biological data, chemical knowledge, and modeling algorithms. Usually the biological data, resulting from good experimental practice, are taken as a true picture of the world; chemical knowledge has scientific bases; so if a QSAR model is not working, blame modeling. The role of modeling in developing scientific theories, and in producing knowledge, is so analyzed. QSAR is a mature technology and is part of a large body of in silico methods and other computational methods. The active debate about the acceptability of the QSAR models, about the way to communicate them, and the explanation to provide accompanies the development of today QSAR models. An example about predicting possible endocrine-disrupting chemicals (EDC) shows the many faces of modern QSAR methods.

Using the local symmetry in amino acids sequences of polypeptides to improve the predictive potential of models of their inhibitor activity

The minimal inhibitory concentrations (pMIC) are a valuable measure of the biological activity of polypeptides. Numerical data on the pMIC are necessary to systematize knowledge on polypeptides' biochemical behaviour. The model of negative decimal logarithm of pMIC of polypeptides in the form of a mathematical function of a sequence of amino acids is suggested. The suggested model is based on the so-called correlation weights of amino acids together with the correlation weights of fragments of local symmetry (FLS). Three kinds of the FLS are considered: (i) three-symbol fragments '…xyx…', (ii) four-symbol fragments '…xyyx…', and (iii) five-symbol fragments '…xyzyx…'. The models built using the Monte Carlo technique improved by applying the index of ideality of correlation (IIC) and correlation intensity index (CII).

First structure-activity relationship analysis of SARS-CoV-2 virus main protease (Mpro) inhibitors: an endeavor on COVID-19 drug discovery

Main protease (Mpro) of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) intervenes in the replication and transcription processes of the virus. Hence, it is a lucrative target for anti-viral drug development. In this study, molecular modeling analyses were performed on the structure activity data of recently reported diverse SARS-CoV-2 Mpro inhibitors to understand the structural requirements for higher inhibitory activity. The classification-based quantitative structure-activity relationship (QSAR) models were generated between SARS-CoV-2 Mpro inhibitory activities and different descriptors. Identification of structural fingerprints to increase or decrease in the inhibitory activity was mapped for possible inclusion/exclusion of these fingerprints in the lead optimization process. Challenges in ADME properties of protease inhibitors were also discussed to overcome the problems of oral bioavailability. Further, depending on the modeling results, we have proposed novel as well as potent SARS-CoV-2 Mpro inhibitors.

The index of ideality of correlation improves the predictive potential of models of the antioxidant activity of tripeptides from frog skin (Litoria rubella)

It is usually held that good-quality models for the biological activity of peptides must take into account their 3D architecture and descriptors of quantum mechanics. However, the present study shows that it is possible to build up models without these complex calculations. The structure of tripeptides represented by sequences of one-symbol abbreviations of the corresponding amino acids serves to build up quantitative structure-activity relationships for the antioxidant activity of tripeptides from frog skin. The statistical quality of the best model for the validation set is n = 27, r² = 0.93, RMSE = 0.15.

The sequence of amino acids as the basis for the model of biological activity of peptides

The algorithm of building up a model for the biological activity of peptides as a mathematical function of a sequence of amino acids is suggested. The general scheme is the following: The total set of available data is distributed into the active training set, passive training set, calibration set, and validation set. The training (both active and passive) and calibration sets are a system of generation of a model of biological activity where each amino acid obtains special correlation weight. The numerical data on the correlation weights calculated by the Monte Carlo method using the CORAL software (http://www.insilico.eu/coral). The target function aimed to give the best result for the calibration set (not for the training set). The final checkup of the model is carried out with data on the validation set (peptides, which are not visible during the creation of the model). Described computational experiments confirm the ability of the approach to be a tool for the design of predictive models for the biological activity of peptides (expressed by pIC50).

The Monte Carlo Method as a Tool to Build up Predictive QSPR/QSAR

Background:

The Monte Carlo method has a wide application in various scientific researches. For the development of predictive models in a form of the quantitative structure-property / activity relationships (QSPRs/QSARs), the Monte Carlo approach also can be useful. The CORAL software provides the Monte Carlo calculations aimed to build up QSPR/QSAR models for different endpoints.

Methods:

Molecular descriptors are a mathematical function of so-called correlation weights of various molecular features. The numerical values of the correlation weights give the maximal value of a target function. The target function leads to a correlation between endpoint and optimal descriptor for the visible training set. The predictive potential of the model is estimated with the validation set, i.e. compounds that are not involved in the process of building up the model.

Results:

The approach gave quite good models for a large number of various physicochemical, biochemical, ecological, and medicinal endpoints. Bibliography and basic statistical characteristics of several CORAL models are collected in the present review. In addition, the extended version of the approach for more complex systems (nanomaterials and peptides), where behaviour of systems is defined by a group of conditions besides the molecular structure is demonstrated.

Conclusion:

The Monte Carlo technique available via the CORAL software can be a useful and convenient tool for the QSPR/QSAR analysis.

Pesticides, cosmetics, drugs: identical and opposite influences of various molecular features as measures of endpoints similarity and dissimilarity

The similarity is an important category in natural sciences. A measure of similarity for a group of various biochemical endpoints is suggested. The list of examined endpoints contains (1) toxicity of pesticides towards rainbow trout; (2) human skin sensitization; (3) mutagenicity; (4) toxicity of psychotropic drugs; and (5) anti HIV activity. Further applying and evolution of the suggested approach is discussed. In particular, the conception of the similarity (dissimilarity) of endpoints can play the role of a "useful bridge" between quantitative structure property/activity relationships (QSPRs/QSARs) and read-across technique.

De Novo Design and In Vitro Testing of Antimicrobial Peptides against Gram-Negative Bacteria

Antimicrobial peptides (AMPs) have been identified as a potentially new class of antibiotics to combat bacterial resistance to conventional drugs. The design of de novo AMPs with high therapeutic indexes, low cost of synthesis, high resistance to proteases and high bioavailability remains a challenge. Such design requires computational modeling of antimicrobial properties. Currently, most computational methods cannot accurately calculate antimicrobial potency against particular strains of bacterial pathogens. We developed a tool for AMP prediction (Special Prediction (SP) tool) and made it available on our Web site (https://dbaasp.org/prediction). Based on this tool, a simple algorithm for the design of de novo AMPs (DSP) was created. We used DSP to design short peptides with high therapeutic indexes against gram-negative bacteria. The predicted peptides have been synthesized and tested in vitro against a panel of gram-negative bacteria, including drug resistant ones. Predicted activity against Escherichia coli ATCC 25922 was experimentally confirmed for 14 out of 15 peptides. Further improvements for designed peptides included the synthesis of D-enantiomers, which are traditionally used to increase resistance against proteases. One synthetic D-peptide (SP15D) possesses one of the lowest values of minimum inhibitory concentration (MIC) among all DBAASP database short peptides at the time of the submission of this article, while being highly stable against proteases and having a high therapeutic index. The mode of anti-bacterial action, assessed by fluorescence microscopy, shows that SP15D acts similarly to cell penetrating peptides. SP15D can be considered a promising candidate for the development of peptide antibiotics. We plan further exploratory studies with the SP tool, aiming at finding peptides which are active against other pathogenic organisms.

Recent trends and analytical challenges in plant bioactive peptide separation, identification and validation

Interest in research into bioactive peptides (BPs) is growing because of their health-promoting ability. Several bioactivities have been ascribed to peptides, including antioxidant, antihypertensive and antimicrobial properties. As they can be produced from precursor proteins, the investigation of BPs in foods is becoming increasingly popular. For the same reason, production of BPs from by-products has also emerged as a possible means of reducing waste and recovering value-added compounds suitable for functional food production and supplements. Milk, meat and fish are the most investigated sources of BPs, but vegetable-derived peptides are also of interest. Vegetables are commonly consumed, and agro-industrial wastes constitute a cheap, large and lower environmental impact source of proteins. The use of advanced analytical techniques for separation and identification of peptides would greatly benefit the discovery of new BPs. In this context, this review provides an overview of the most recent applications in BP investigations for vegetable food and by-products. The most important issues regarding peptide isolation and separation, by single or multiple chromatographic techniques, are discussed. Additionally, problems connected with peptide identification in plants and non-model plants are discussed regarding the particular case of BP identification. Finally, the issue of peptide validation to confirm sequence and bioactivity is presented. Graphical representation of the analytical workflow needed for investigation of bioactive peptides and applied to vegetables and vegetable wastes Graphical Abstract.

QSPR/QSAR Analyses by Means of the CORAL Software

In this chapter, the methodology of building up quantitative structure—property/activity relationships (QSPRs/QSARs)—by means of the CORAL software is described. The Monte Carlo method is the basis of this approach. Simplified Molecular Input-Line Entry System (SMILES) is used as the representation of the molecular structure. The conversion of SMILES into the molecular graph is available for QSPR/QSAR analysis using the CORAL software. The model for an endpoint is a mathematical function of the correlation weights for various features of the molecular structure. Hybrid models that are based on features extracted from both SMILES and a graph also can be built up by the CORAL software. The conceptually new ideas collected and revealed through the CORAL software are: (1) any QSPR/QSAR model is a random event; and (2) optimal descriptor can be a translator of eclectic information into an endpoint prediction.

Enabling the Discovery and Virtual Screening of Potent and Safe Antimicrobial Peptides. Simultaneous Prediction of Antibacterial Activity and Cytotoxicity

Antimicrobial peptides (AMPs) represent promising alternatives to fight against bacterial pathogens. However, cellular toxicity remains one of the main concerns in the early development of peptide-based drugs. This work introduces the first multitasking (mtk) computational model focused on performing simultaneous predictions of antibacterial activities, and cytotoxicities of peptides. The model was created from a data set containing 3592 cases, and it displayed accuracy higher than 96% for classifying/predicting peptides in both training and prediction (test) sets. The technique known as alanine scanning was computationally applied to illustrate the calculation of the quantitative contributions of the amino acids (in their respective positions of the sequence) to the biological effects of a defined peptide. A small library formed by 10 peptides was generated, where peptides were designed by considering the interpretations of the different descriptors in the mtk-computational model. All the peptides were predicted to exhibit high antibacterial activities against multiple bacterial strains, and low cytotoxicity against various cell types. The present mtk-computational model can be considered a very useful tool to support high throughput research for the discovery of potent and safe AMPs.

First Multitarget Chemo-Bioinformatic Model To Enable the Discovery of Antibacterial Peptides against Multiple Gram-Positive Pathogens

Antimicrobial peptides (AMPs) have emerged as promising therapeutic alternatives to fight against the diverse infections caused by different pathogenic microorganisms. In this context, theoretical approaches in bioinformatics have paved the way toward the creation of several in silico models capable of predicting antimicrobial activities of peptides. All current models have several significant handicaps, which prevent the efficient search for highly active AMPs. Here, we introduce the first multitarget (mt) chemo-bioinformatic model devoted to performing alignment-free prediction of antibacterial activity of peptides against multiple Gram-positive bacterial strains. The model was constructed from a data set containing 2488 cases of AMPs sequences assayed against at least 1 out of 50 Gram-positive bacterial strains. This mt-chemo-bioinformatic model displayed percentages of correct classification higher than 90.00% in both training and prediction (test) sets. For the first time, two computational approaches derived from basic concepts in genetics and molecular biology were applied, allowing the calculations of the relative contributions of any amino acid (in a defined position) to the antibacterial activity of an AMP and depending on the bacterial strain used in the biological assay. The present mt-chemo-bioinformatic model constitutes a powerful tool to enable the discovery of potent and versatile AMPs.

CORAL: model for no observed adverse effect level (NOAEL)

The in vivo repeated dose toxicity (RDT) test is intended to provide information on the possible risk caused by repeated exposure to a substance over a limited period of time. The measure of the RDT is the no observed adverse effect level (NOAEL) that is the dose at which no effects are observed, i.e., this endpoint indicates the safety level for a substance. The need to replace in vivo tests, as required by some European Regulations (registration, evaluation authorization and restriction of chemicals) is leading to the searching for reliable alternative methods such as quantitative structure-activity relationships (QSAR). Considering the complexity of the RDT endpoint, for which data quality is limited and depends anyway on the study design, the development of QSAR for this endpoint is an attractive task. Starting from a dataset of 140 organic compounds with NOAEL values related to oral short term toxicity in rats, we developed a QSAR model based on optimal descriptors calculated with simplified molecular input-line entry systems and the graph of atomic orbitals by the Monte Carlo method, using CORAL software. Three different splits into the training, calibration, and validation sets are studied. The mechanistic interpretation of these models in terms of molecular fragment with positive or negative contributions to the endpoint is discussed. The probabilistic definition for the domain of applicability is suggested.

QSPR/QSAR Analyses by Means of the CORAL Software

In this chapter, the methodology of building up quantitative structure—property/activity relationships (QSPRs/QSARs)—by means of the CORAL software is described. The Monte Carlo method is the basis of this approach. Simplified Molecular Input-Line Entry System (SMILES) is used as the representation of the molecular structure. The conversion of SMILES into the molecular graph is available for QSPR/QSAR analysis using the CORAL software. The model for an endpoint is a mathematical function of the correlation weights for various features of the molecular structure. Hybrid models that are based on features extracted from both SMILES and a graph also can be built up by the CORAL software. The conceptually new ideas collected and revealed through the CORAL software are: (1) any QSPR/QSAR model is a random event; and (2) optimal descriptor can be a translator of eclectic information into an endpoint prediction.

Comprehension of drug toxicity: software and databases

Quantitative structure-property/activity relationships (QSPRs/QSARs) are a tool (in silico) to rapidly predict various endpoints in general, and drug toxicity in particular. However, this dynamic evolution of experimental data (expansion of existing experimental data on drugs toxicity) leads to the problem of critical estimation of the data. The carcinogenicity, mutagenicity, liver effects and cardiac toxicity should be evaluated as the most important aspects of the drug toxicity. The toxicity is a multidimensional phenomenon. It is apparent that the main reasons for the increase in applications of in silico prediction of toxicity include the following: (i) the need to reduce animal testing; (ii) computational models provide reliable toxicity prediction; (iii) development of legislation that is related to use of new substances; (iv) filling data gaps; (v) reduction of cost and time; (vi) designing of new compounds; (vii) advancement of understanding of biology and chemistry. This mini-review provides analysis of existing databases and software which are necessary for use of robust computational assessments and robust prediction of potential drug toxicities by means of in silico methods.

In silico methods to predict drug toxicity

This review describes in silico methods to characterize the toxicity of pharmaceuticals, including tools which predict toxicity endpoints such as genotoxicity or organ-specific models, tools addressing ADME processes, and methods focusing on protein-ligand docking binding. These in silico tools are rapidly evolving. Nowadays, the interest has shifted from classical studies to support toxicity screening of candidates, toward the use of in silico methods to support the expert. These methods, previously considered useful only to provide a rough, initial estimation, currently have attracted interest as they can assist the expert in investigating toxic potential. They provide the expert with safety perspectives and insights within a weight-of-evidence strategy. This represents a shift of the general philosophy of in silico methodology, and it is likely to further evolve especially exploiting links with system biology.