Phenolic Acid-β-Cyclodextrin Complexation Study to Mask Bitterness in Wheat Bran: A Machine Learning-Based QSAR Study

The need to solvate and encapsulate hydro-sensitive molecules drives noticeable trends in the applications of cyclodextrins in the pharmaceutical industry, in foods, polymers, materials, and in agricultural science. Among them, β-cyclodextrin is one of the most used for the entrapment of phenolic acid compounds to mask the bitterness of wheat bran. In this regard, there is still a need for good data and especially for a robust predictive model that assesses the bitterness masking capabilities of β-cyclodextrin for various phenolic compounds. This study uses a dataset of 20 phenolic acids docked into the β-cyclodextrin cavity to generate three different binding constants. The data from the docking study were combined with topological, topographical, and quantum-chemical features from the ligands in a machine learning-based structure-activity relationship study. Three different models for each binding constant were computed using a combination of the genetic algorithm (GA) and multiple linear regression (MLR) approaches. The developed ML/QSAR models showed a very good performance, with high predictive ability and correlation coefficients of 0.969 and 0.984 for the training and test sets, respectively. The models revealed several factors responsible for binding with cyclodextrin, showing positive contributions toward the binding affinity values, including such features as the presence of six-membered rings in the molecule, branching, electronegativity values, and polar surface area.

Development of a Selective and Stable Antimicrobial Peptide

Antimicrobial peptides (AMPs) are presented as potential scaffolds for antibiotic development due to their desirable qualities including broad-spectrum activity, rapid action, and general lack of susceptibility to current resistance mechanisms. However, they often lose antibacterial activity under physiological conditions and/or display mammalian cell toxicity, which limits their potential use. Identification of AMPs that overcome these barriers will help develop rules for how this antibacterial class can be developed to treat infection. Here we describe the development of our novel synthetic AMP, from discovery through in vivo application. Our evolved AMP, DTr18-dab, has broad-spectrum antibacterial activity and is nonhemolytic. It is active against planktonic bacteria and biofilm, is unaffected by colistin resistance, and importantly is active in both human serum and a Galleria mellonella infection model. Several modifications, including the incorporation of noncanonical amino acids, were used to arrive at this robust sequence. We observed that the impact on antibacterial activity with noncanonical amino acids was dependent on assay conditions and therefore not entirely predictable. Overall, our results demonstrate how a relatively weak lead can be developed into a robust AMP with qualities important for potential therapeutic translation.

Quantitative Structure-Permittivity Relationship Study of a Series of Polymers

Dielectric constant is an important property which is widely utilized in many scientific fields and characterizes the degree of polarization of substances under the external electric field. In this work, a structure-property relationship of the dielectric constants (ε) for a diverse set of polymers was investigated. A transparent mechanistic model was developed with the application of a machine learning approach that combines genetic algorithm and multiple linear regression analysis, to obtain a mechanistically explainable and transparent model. Based on the evaluation conducted using various validation criteria, four- and eight-variable models were proposed. The best model showed a high predictive performance for training and test sets, with R2 values of 0.905 and 0.812, respectively. Obtained statistical performance results and selected descriptors in the best models were analyzed and discussed. With the validation procedures applied, the models were proven to have a good predictive ability and robustness for further applications in polymer permittivity prediction.

Stability Constant and Potentiometric Sensitivity of Heavy Metal-Organic Fluorescent Compound Complexes: QSPR Models for Prediction and Design of Novel Coumarin-like Ligands

Industrial wastewater often consists of toxic chemicals and pollutants, which are extremely harmful to the environment. Heavy metals are toxic chemicals and considered one of the major hazards to the aquatic ecosystem. Analytical techniques, such as potentiometric methods, are some of the methods to detect heavy metals in wastewaters. In this work, the quantitative structure-property relationship (QSPR) was applied using a range of machine learning techniques to predict the stability constant (logβ_ML) and potentiometric sensitivity (PS_ML) of 200 ligands in complexes with the heavy metal ions Cu²⁺, Cd²⁺, and Pb²⁺. In result, the logβML models developed for four ions showed good performance with square correlation coefficients (R²) ranging from 0.80 to 1.00 for the training and 0.72 to 0.85 for the test sets. Likewise, the PSML displayed acceptable performance with an R² of 0.87 to 1.00 for the training and 0.73 to 0.95 for the test sets. By screening a virtual database of coumarin-like structures, several new ligands bearing the coumarin moiety were identified. Three of them, namely NEW02, NEW03, and NEW07, showed very good sensitivity and stability in the metal complexes. Subsequent quantum-chemical calculations, as well as physicochemical/toxicological profiling were performed to investigate their metal-binding ability and developability of the designed sensors. Finally, synthesis schemes are proposed to obtain these three ligands with major efficiency from simple resources. The three coumarins designed clearly demonstrated capability to be suitable as good florescent chemosensors towards heavy metals. Overall, the computational methods applied in this study showed a very good performance as useful tools for designing novel fluorescent probes and assessing their sensing abilities.

Development of quantitative structure-retention relationship models to improve the identification of leachables in food packaging using non-targeted analysis

Quantitative structure-retention relationship (QSRR) models can be used to predict the chromatographic retention time of chemicals and facilitate the identification of unknown compounds, notably with non-targeted analysis. In this study, QSRR models were developed from the data obtained for 178 pure chemical standards and four types of analytical columns (C18, phenylhexyl, pentafluorophenyl, cyano) in liquid chromatography quadrupole time-of-flight mass spectrometry (LC-Q-TOF-MS). First, different data partitioning ratios and feature selection methods [random forest (RF) and support vector machine (SVM)] were tested to build models to predict chromatographic retention times based on 2D molecular descriptors. The internal and external performances of the non-linear (RF) and corresponding linear predictive models were systematically compared, and RF models resulted in better predictive capacities [p < 0.05, with an average PVE (proportion of variance explained) value of 0.89 ± 0.02] than linear models (0.79 ± 0.03). For each column, the resulting model was applied to identify leachables from actual plastic packaging samples. An in-depth investigation of the top 20 most intense molecular features revealed that all false-positives could be identified as outliers in the QSRR models (outside of the 95% prediction bands). Furthermore, analyzing a sample on multiple chromatographic columns and applying the associated QSRR models increased the capacity to filter false positives. Such an approach will contribute to a more effective identification of unknown or unexpected leachables in plastics (e.g. non-intended added substances), therefore refining our understanding of the chemical risks associated with food contact materials.

Principles and Procedures for Assessment of Acute Toxicity Incorporating In Silico Methods

Acute toxicity in silico models are being used to support an increasing number of application areas including (1) product research and development, (2) product approval and registration as well as (3) the transport, storage and handling of chemicals. The adoption of such models is being hindered, in part, because of a lack of guidance describing how to perform and document an in silico analysis. To address this issue, a framework for an acute toxicity hazard assessment is proposed. This framework combines results from different sources including in silico methods and in vitro or in vivo experiments. In silico methods that can assist the prediction of in vivo outcomes (i.e., LD50) are analyzed concluding that predictions obtained using in silico approaches are now well-suited for reliably supporting assessment of LD50-based acute toxicity for the purpose of GHS classification. A general overview is provided of the endpoints from in vitro studies commonly evaluated for predicting acute toxicity (e.g., cytotoxicity/cytolethality as well as assays targeting specific mechanisms). The increased understanding of pathways and key triggering mechanisms underlying toxicity and the increased availability of in vitro data allow for a shift away from assessments solely based on endpoints such as LD50, to mechanism-based endpoints that can be accurately assessed in vitro or by using in silico prediction models. This paper also highlights the importance of an expert review of all available information using weight-of-evidence considerations and illustrates, using a series of diverse practical use cases, how in silico approaches support the assessment of acute toxicity.

The Hydrolysis Rate of Paraoxonase-1 Q and R Isoenzymes: An In Silico Study Based on In Vitro Data

Human serum paraoxonase-1 (PON1) is an important hydrolase-type enzyme found in numerous tissues. Notably, it can exist in two isozyme-forms, Q and R, that exhibit different activities. This study presents an in silico (QSAR, Docking, MD and QM/MM) study of a set of compounds on the activity towards the PON1 isoenzymes (QPON1 and RPON1). Different rates of reaction for the Q and R isoenzymes were analyzed by modelling the effect of Q192R mutation on active sites. It was concluded that the Q192R mutation is not even close to the active site, while it is still changing the geometry of it. Using the combined genetic algorithm with multiple linear regression (GA-MLR) technique, several QSAR models were developed and relative activity rates of the isozymes of PON1 explained. From these, two QSAR models were selected, one each for the QPON1 and RPON1. Best selected models are four-variable MLR models for both Q and R isozymes with squared correlation coefficient R² values of 0.87 and 0.83, respectively. In addition, the applicability domain of the models was analyzed based on the Williams plot. The results were discussed in the light of the main factors that influence the hydrolysis activity of the PON1 isozymes.

Identification of concealed structural alerts using QSTR modeling for Pseudokirchneriella subcapitata

In the present work, QSTR modeling was conducted for microalga Pseudokirchneriella subcapitata using a data set of 271 molecules belonging to different types of chemical classes for the prediction of EC₅₀ for 72 hr based assays. The balanced QSTR model encompasses seven easily interpretable molecular descriptors and possesses statistical robustness with high predictive ability. This Genetic Algorithm Multi-linear regression (GA-MLR) model was subjected to internal validation, Y-randomization test, applicability domain analysis, and external validation as per the recommended OECD guidelines. The newly developed model fulfilled the threshold values for more than 20 recommended validation parameters including R² = 0.72, Q²_LOO = 0.70, etc. The developed QSTR model was successful in identifying the type of hybridization or specific type of atoms of previously reported and newer structural alerts. Thus, the model could be useful for data gap filling and expanding mechanistic interpretation of toxicity for different chemicals.

In vivo toxicity of nitroaromatics: A comprehensive quantitative structure-activity relationship study

The toxicity data of 90 nitroaromatic compounds related to their 50% lethal dose concentration for rats (LD50) were analyzed to develop quantitative structure-activity relationship (QSAR) models. Quantum-chemically calculated descriptors together with molecular descriptors generated by DRAGON, PaDEL, and HiT-QSAR software were utilized to build QSAR models. Quality and validity of the models were determined by internal and external validation techniques. The results show that the toxicity of nitroaromatic compounds depends on various factors, such as the number of nitro-groups, the topological state, and the presence of certain structural fragments. The developed models based on the largest (to date) dataset of nitroaromatics in vivo toxicity showed a good predictive ability. The results provide important input that could be applied in a preliminary assessment of nitroaromatic compounds' toxicity to mammals. Environ Toxicol Chem 2017;36:2227-2233. © 2017 SETAC.

Recent Developments in 3D QSAR and Molecular Docking Studies of Organic and Nanostructures

The development of quantitative structure–activity relationship (QSAR) methods is going very fast for the last decades. OSAR approach already plays an important role in lead structure optimization, and nowadays, with development of big data approaches and computer power, it can even handle a huge amount of data associated with combinatorial chemistry. One of the recent developments is a three-dimensional QSAR, i.e., 3D QSAR. For the last two decades, 3D-OSAR has already been successfully applied to many datasets, especially of enzyme and receptor ligands. Moreover, quite often 3D QSAR investigations are going together with protein–ligand docking studies and this combination works synergistically. In this review, we outline recent advances in development and applications of 3D QSAR and protein–ligand docking approaches, as well as combined approaches for conventional organic compounds and for nanostructured materials, such as fullerenes and carbon nanotubes.

Polymer Coating Materials and Their Fouling Release Activity: A Cheminformatics Approach to Predict Properties

A novel cheminformatics-based approach has been employed to investigate a set of polymer coating materials designed to mitigate the accumulation of marine biofouling on surfaces immersed in the sea. Specifically, a set of 27 nontoxic, amphiphilic polysiloxane-based polymer coatings was synthesized using a combinatorial, high-throughput approach and characterized for fouling-release (FR) activity toward a number of relevant marine fouling organisms, including bacteria, microalgae, and adult barnacles. In order to model these complex systems adequately, a new computational technique was used in which all investigated polymer-based coating materials were considered as mixture systems comprising several compositional variables at a range of concentrations. By applying a combination of methodologies for mixture systems and a quantitative structure-activity relationship approach (QSAR), seven unique QSAR models were developed that were able to successfully predict the desired FR properties. Furthermore, the developed models identified several significant descriptors responsible for FR activity of investigated polymer-based coating materials, with correlation coefficients ranging from r_test² = 0.63 to 0.94. The computational models derived from this study may serve as a powerful set of tools to predict optimal combinations of source components to produce amphiphilic polysiloxane-based coating systems with effective, broad-spectrum FR properties.

A Quantitative Structure Activity Relationship for acute oral toxicity of pesticides on rats: Validation, domain of application and prediction

Quantitative Structure Activity Relationship (QSAR) models are expected to play an important role in the risk assessment of chemicals on humans and the environment. In this study, we developed a validated QSAR model to predict acute oral toxicity of 329 pesticides to rats because a few QSAR models have been devoted to predict the Lethal Dose 50 (LD50) of pesticides on rats. This QSAR model is based on 17 molecular descriptors, and is robust, externally predictive and characterized by a good applicability domain. The best results were obtained with a 17/9/1 Artificial Neural Network model trained with the Quasi Newton back propagation (BFGS) algorithm. The prediction accuracy for the external validation set was estimated by the Q(2)ext and the root mean square error (RMS) which are equal to 0.948 and 0.201, respectively. 98.6% of external validation set is correctly predicted and the present model proved to be superior to models previously published. Accordingly, the model developed in this study provides excellent predictions and can be used to predict the acute oral toxicity of pesticides, particularly for those that have not been tested as well as new pesticides.

Amino substituted nitrogen heterocycle ureas as kinase insert domain containing receptor (KDR) inhibitors: Performance of structure–activity relationship approaches

A quantitative structure–activity relationship (QSAR) study was performed on a set of amino-substituted nitrogen heterocyclic urea derivatives. Two novel approaches were applied: (1) the simplified molecular input-line entry systems (SMILES) based optimal descriptors approach; and (2) the fragment-based simplex representation of molecular structure (SiRMS) approach. Comparison with the classic scheme of building up the model and balance of correlation (BC) for optimal descriptors approach shows that the BC scheme provides more robust predictions than the classic scheme for the considered pIC₅₀ of the heterocyclic urea derivatives. Comparison of the SMILES-based optimal descriptors and SiRMS approaches has confirmed good performance of both techniques in prediction of kinase insert domain containing receptor (KDR) inhibitory activity, expressed as a logarithm of inhibitory concentration (pIC50) of studied compounds.

Aconitum and delphinium diterpenoid alkaloids of local anesthetic activity: comparative QSAR analysis based on GA-MLRA/PLS and optimal descriptors approach

The duration of anesthesia (related to protein binding of a drug) and the onset time (determined by the pKa) are important characteristics in assessment of local anesthetic agents. They are known to be affected by a number of factors. Early studies of antiarrhythmic diterpenoid alkaloids from plants Aconitum and Delphinium suggested that they possess local anesthetic activity due to their ability to suppress sodium currents of excited membranes. In this study we utilized toxicity, duration, and onset of action as endpoints to construct Quantitative Structure-Activity Relationship (QSAR) models for the series of 34 diterpenoid alkaloids characterized by local anesthetic activity using genetic algorithm-based multiple linear regression analysis/partial least squares and simplified molecular input line entry system (SMILES)-based optimal descriptors approach. The developed QSAR models correctly reflected factors that determine three endpoints of interest. Toxicity correlates with descriptors describing partition and reactivity of compounds. The duration of anesthesia was encoded by the parameters defining the ability of a compound to bind at the receptor site. The size and number of H-bond acceptor atoms were found not to favor the speed of onset, while topographic electronic descriptor demonstrated strong positive effect on it. SMILES-based optimal descriptors approach resulted in overall improvement of models. This approach was shown to be more sensitive to structural peculiarities of molecules than regression methods. The results clearly indicate that obtained QSARs are able to provide distinct rationales for compounds optimization with respect to particular endpoint.