Adsorption of Peptides onto Carbon Nanotubes Grafted with Poly(ethylene Oxide) Chains: A Molecular Dynamics Simulation Study

Carbon nanotubes (CNTs) display exceptional properties that predispose them to wide use in technological or biomedical applications. To remove the toxicity of CNTs and to protect them against undesired protein adsorption, coverage of the CNT sidewall with poly(ethylene oxide) (PEO) is often considered. However, controversial results on the antifouling effectiveness of PEO layers have been reported so far. In this work, the interactions of pristine CNT and CNT covered with the PEO chains at different grafting densities with polyglycine, polyserine, and polyvaline are studied using molecular dynamics simulations in vacuum, water, and saline environments. The peptides are adsorbed on CNT in all investigated systems; however, the adsorption strength is reduced in aqueous environments. Save for one case, addition of NaCl at a physiological concentration to water does not appreciably influence the adsorption and structure of the peptides or the grafted PEO layer. It turns out that the flexibility of the peptide backbone allows the peptide to adopt more asymmetric conformations which may be inserted deeper into the grafted PEO layer. Water molecules disrupt the internal hydrogen bonds in the peptides, as well as the hydrogen bonds formed between the peptides and the PEO chains.

Toxicity Assessment of the Binary Mixtures of Aquatic Organisms Based on Different Hypothetical Descriptors

Modern industrialization has led to the creation of a wide range of organic chemicals, especially in the form of multicomponent mixtures, thus making the evaluation of environmental pollution more difficult by normal methods. In this paper, we attempt to use forward stepwise multiple linear regression (MLR) and nonlinear radial basis function neural networks (RBFNN) to establish quantitative structure–activity relationship models (QSARs) to predict the toxicity of 79 binary mixtures of aquatic organisms using different hypothetical descriptors. To search for the proper mixture descriptors, 11 mixture rules were performed and tested based on preliminary modeling results. The statistical parameters of the best derived MLR model were Ntrain = 62, R2 = 0.727, RMS = 0.494, F = 159.537, Q2LOO = 0.727, and Q2pred = 0.725 for the training set; and Ntest = 17, R2 = 0.721, RMS = 0.508, F = 38.773, and q2ext = 0.720 for the external test set. The RBFNN model gave the following statistical results: Ntrain = 62, R2 = 0.956, RMS = 0.199, F = 1279.919, Q2LOO = 0.955, and Q2pred = 0.855 for the training set; and Ntest = 17, R2 = 0.880, RMS = 0.367, F = 110.980, and q2ext = 0.853 for the external test set. The quality of the models was assessed by validating the relevant parameters, and the final results showed that the developed models are predictive and can be used for the toxicity prediction of binary mixtures within their applicability domain.

Computational Modelling and Sustainable Synthesis of a Highly Selective Electrochemical MIP-Based Sensor for Citalopram Detection

A novel molecularly imprinted polymer (MIP) has been developed based on a simple and sustainable strategy for the selective determination of citalopram (CTL) using screen-printed carbon electrodes (SPCEs). The MIP layer was prepared by electrochemical in situ polymerization of the 3-amino-4 hydroxybenzoic acid (AHBA) functional monomer and CTL as a template molecule. To simulate the polymerization mixture and predict the most suitable ratio between the template and functional monomer, computational studies, namely molecular dynamics (MD) simulations, were carried out. During the experimental preparation process, essential parameters controlling the performance of the MIP sensor, including CTL:AHBA concentration, number of polymerization cycles, and square wave voltammetry (SWV) frequency were investigated and optimized. The electrochemical characteristics of the prepared MIP sensor were evaluated by both cyclic voltammetry (CV) and electrochemical impedance spectroscopy (EIS) techniques. Based on the optimal conditions, a linear electrochemical response of the sensor was obtained by SWV measurements from 0.1 to 1.25 µmol L-1 with a limit of detection (LOD) of 0.162 µmol L-1 (S/N = 3). Moreover, the MIP sensor revealed excellent CTL selectivity against very close analogues, as well as high imprinting factor of 22. Its applicability in spiked river water samples demonstrated its potential for adequate monitoring of CTL. This sensor offers a facile strategy to achieve portability while expressing a willingness to care for the environment.

Moving Average-Based Multitasking In Silico Classification Modeling: Where Do We Stand and What Is Next?

Conventional in silico modeling is often viewed as 'one-target' or 'single-task' computer-aided modeling since it mainly relies on forecasting an endpoint of interest from similar input data. Multitasking or multitarget in silico modeling, in contrast, embraces a set of computational techniques that efficiently integrate multiple types of input data for setting up unique in silico models able to predict the outcome(s) relating to various experimental and/or theoretical conditions. The latter, specifically, based upon the Box-Jenkins moving average approach, has been applied in the last decade to several research fields including drug and materials design, environmental sciences, and nanotechnology. The present review discusses the current status of multitasking computer-aided modeling efforts, meanwhile describing both the existing challenges and future opportunities of its underlying techniques. Some important applications are also discussed to exemplify the ability of multitasking modeling in deriving holistic and reliable in silico classification-based models as well as in designing new chemical entities, either through fragment-based design or virtual screening. Focus will also be given to some software recently developed to automate and accelerate such types of modeling. Overall, this review may serve as a guideline for researchers to grasp the scope of multitasking computer-aided modeling as a promising in silico tool.

N2O Hydrogenation on Silver Doped Gold Catalysts, A DFT Study

In this study, the full reaction mechanism for N₂O hydrogenation on silver doped Au(210) surfaces was investigated in order to clarify the experimental observations. Density functional theory (DFT) calculations were used to state the most favorable reaction paths for individual steps involved in the N₂O hydrogenation. From the DFT results, the activation energy barriers, rate constants and reaction energies for the individual steps were determined, which made it possible to elucidate the most favorable reaction mechanism for the global catalytic process. It was found that the N₂O dissociation occurs in surface regions where silver atoms are present, while hydrogen dissociation occurs in pure gold regions of the catalyst or in regions with a low silver content. Likewise, N₂O dissociation is the rate determining step of the global process, while water formation from O adatoms double hydrogenation and N₂ and H₂O desorptions are reaction steps limited by low activation energy barriers, and therefore, the latter are easily carried out. Moreover, water formation occurs in the edges between the regions where hydrogen and N₂O are dissociated. Interestingly, a good dispersion of the silver atoms in the surface is necessary to avoid catalyst poison by O adatoms accumulation, which are strongly adsorbed on the surface.

Structural behavior of monomer of SARS-CoV-2 spike protein during initial stage of adsorption on graphene

Spike glycoprotein of the SARS-CoV-2 virus and its structure play a crucial role in the infections of cells containing angiotensin-converting enzyme 2 (ACE2) as well as in the interactions of this virus with surfaces. Protection against viruses and often even their deactivation is one of the great varieties of graphene applications. The structural changes of the non-glycosylated monomer of the spike glycoprotein trimer (denoted as S-protein in this work) triggered by its adsorption onto graphene at the initial stage are investigated by means of atomistic molecular dynamics simulations. The adsorption of the S-protein happens readily during the first 10 ns. The shape of the S-protein becomes more prolate during the adsorption, but this trend, albeit less pronounced, is observed also for the freely relaxing S-protein in water. The receptor-binding domain (RBD) of the free and adsorbed S-protein manifests itself as the most rigid fragment of the whole S-protein. The adsorption even enhances the rigidity of the whole S-protein as well as its subunits. Only one residue of the RBD involved in the specific interactions with ACE2 during the cell infection is involved in the direct contact of the adsorbed S-protein with the graphene. The new intramolecular hydrogen bonds formed during the S-protein adsorption replace the S-protein-water hydrogen bonds; this trend, although less apparent, is observed also during the relaxation of the free S-protein in water. In the initial phase, the secondary structure of the RBD fragment specifically interacting with ACE2 receptor is not affected during the S-protein adsorption onto the graphene.

Developing a Multi-target Model to Predict the Activity of Monoamine Oxidase A and B Drugs

INTRODUCTION

Monoamine oxidase inhibitors (MAOIs) are compounds largely used in the treatment of Parkinson's disease (PD), Alzheimer's disease and other neuropsychiatric disorders since they are closely related to the MAO enzymes activity. The two isoforms of the MAO enzymes, MAO-A and MAO-B, are responsible for the degradation of monoamine neurotransmitters and due to this, relevant efforts have been devoted to finding new compounds with more selectivity and less side effects. One of the most used approaches is based on the use of computational approaches since they are time and money-saving and may allow us to find a more relevant structure-activity relationship.

OBJECTIVE

In this manuscript, we will review the most relevant computational approaches aimed at the prediction and development of new MAO inhibitors. Subsequently, we will also introduce a new multitask model aimed at predicting MAO-A and MAO-B inhibitors.

METHODS

The QSAR multi-task model herein developed was based on the use of the linear discriminant analysis. This model was developed gathering 5,759 compounds from the public dataset Chembl. The molecular descriptors used was calculated using the Dragon software. Classical statistical tests were performed to check the validity and robustness of the model.

RESULTS

The herein proposed model is able to correctly classify all the 5,759 compounds. All the statistical performed tests indicated that this model is robust and reproducible.

CONCLUSION

MAOIs are compounds of large interest since they are largely used in the treatment of very serious illness. These inhibitors may lose efficacy and produce severe side effects. Due to this, the development of selective MAO-A or MAO-B inhibitors is crucial for the treatment of these diseases and their effects. The herein proposed multi-target QSAR model may be a relevant tool in the development of new and more selective MAO inhibitors.

Janus Gold Nanoparticles from Nanodroplets of Alkyl Thiols: A Molecular Dynamics Study

Janus particles provide an asymmetry in structure, which can impart diverse functionalities leading to immense importance in various applications, ranging from targeted delivery to interfacial phenomena, including catalysis, electronics, and optics. In this work, we present results of a molecular dynamics study of the growth mechanism of coating on gold nanoparticles (AuNPs) from droplets of n-alkyl thiols with different chain lengths (C5 and C11) and terminal groups (CH₃ and COOH). The effect of chain lengths and functional groups on the formation of a monolayer of alkyl thiols on AuNPs is investigated. A two-step mechanism, initiated by the binding of the droplet to the nanoparticle surface with a time constant on the order of ∼1 ns, followed by the diffusion-driven growth with a larger time constant (on the order of 100 ns), is shown to capture the growth dynamics of the monolayer. It is observed that the time required for complete wetting increases with an increase in the chain length. Moreover, the monolayer formation is slowed down in the presence of carboxyl groups because of strong hydrogen bonding. The kinetics of the n-alkyl thiols coating on the nanoparticles is found to be independent of the droplet size but carboxyl-terminated thiols spread more with increasing droplet size. Furthermore, different time constants for different chains and functional groups yield Janus coating when two droplets of alkyl thiols with different terminal groups are allowed to form monolayers on the nanoparticle. The Janus balance (β) for different combinations of alkyl thiols and nanoparticle sizes varies in the range of 0.42-0.71.

A Machine Learning Approach for Hot-Spot Detection at Protein-Protein Interfaces

Understanding protein-protein interactions is a key challenge in biochemistry. In this work, we describe a more accurate methodology to predict Hot-Spots (HS) in protein-protein interfaces from their native complex structure compared to previous published Machine Learning (ML) techniques. Our model is trained on a large number of complexes and on a significantly larger number of different structural- and evolutionary sequence-based features. In particular, we added interface size, type of interaction between residues at the interface of the complex, number of different types of residues at the interface and the Position-Specific Scoring Matrix (PSSM), for a total of 79 features. We used twenty-seven algorithms from a simple linear-based function to support-vector machine models with different cost functions. The best model was achieved by the use of the conditional inference random forest (c-forest) algorithm with a dataset pre-processed by the normalization of features and with up-sampling of the minor class. The method has an overall accuracy of 0.80, an F1-score of 0.73, a sensitivity of 0.76 and a specificity of 0.82 for the independent test set.

Multi-Target QSAR Approaches for Modeling Protein Inhibitors. Simultaneous Prediction of Activities Against Biomacromolecules Present in Gram-Negative Bacteria

Drug discovery is aimed at finding therapeutic agents for the treatment of many diverse diseases and infections. However, this is a very slow an expensive process, and for this reason, in silico approaches are needed to rationalize the search for new molecular entities with desired biological profiles. Models focused on quantitative structure-activity relationships (QSAR) have constituted useful complementary tools in medicinal chemistry, allowing the virtual predictions of dissimilar pharmacological activities of compounds. In the last 10 years, multi-target (mt) QSAR models have been reported, representing great advances with respect to those models generated from classical approaches. Thus, mt- QSAR models can simultaneously predict activities against different biological targets (proteins, microorganisms, cell lines, etc.) by using large and heterogeneous datasets of chemicals. The present review is devoted to discuss the most promising mt-QSAR models, particularly those developed for the prediction of protein inhibitors. We also report the first multi-tasking QSAR (mtk-QSAR) model for simultaneous prediction of inhibitors against biomacromolecules (specifically proteins) present in Gram-negative bacteria. This model allowed us to consider both different proteins and multiple experimental conditions under which the inhibitory activities of the chemicals were determined. The mtk-QSAR model exhibited accuracies higher than 98% in both training and prediction sets, also displaying a very good performance in the classification of active and inactive cases that depended on the specific elements of the experimental conditions. The physicochemical interpretations of the molecular descriptors were also analyzed, providing important insights regarding the molecular patterns associated with the appearance/enhancement of the inhibitory potency.

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.

A general ANN-based multitasking model for the discovery of potent and safer antibacterial agents

Bacteria have been one of the world's most dangerous and deadliest pathogens for mankind, nowadays giving rise to significant public health concerns. Given the prevalence of these microbial pathogens and their increasing resistance to existing antibiotics, there is a pressing need for new antibacterial drugs. However, development of a successful drug is a complex, costly, and time-consuming process. Quantitative Structure-Activity Relationships (QSAR)-based approaches are valuable tools for shortening the time of lead compound identification but also for focusing and limiting time-costly synthetic activities and in vitro/vivo evaluations. QSAR-based approaches, supported by powerful statistical techniques such as artificial neural networks (ANNs), have evolved to the point of integrating dissimilar types of chemical and biological data. This chapter reports an overview of the current research and potential applications of QSAR modeling tools toward the rational design of more efficient antibacterial agents. Particular emphasis is given to the setup of multitasking models along with ANNs aimed at jointly predicting different antibacterial activities and safety profiles of drugs/chemicals under diverse experimental conditions.

Role of ligand-based drug design methodologies toward the discovery of new anti- Alzheimer agents: futures perspectives in Fragment-Based Ligand Design

Alzheimer's disease (AD), a degenerative disease affecting the brain, is the single most common source of dementia in adults. The cause and the progression of AD still remains a mystery among medical experts. As a result, a cure has not yet been discovered, even after decade's worth of research that started since 1906, when the disease was first identified. Despite the efforts of the scientific community, several of the biological receptors associated with AD have not been sufficiently studied to date, limiting in turn the design of new and more potent anti-AD agents. Thus, the search for new drug candidates as inhibitors of different targets associated with AD constitutes an essential part towards the discovery of new and more efficient anti-AD therapies. The present work is focused on the role of the Ligand-Based Drug Design (LBDD) methodologies which have been applied for the elucidation of new molecular entities with high inhibitory activity against targets related with AD. Particular emphasis is given also to the current state of fragment-based ligand approaches as alternatives of the Fragment-Based Drug Discovery (FBDD) methodologies. Finally, several guidelines are offered to show how the use of fragment-based descriptors can be determinant for the design of multi-target inhibitors of proteins associated with AD.

Quantitative structure-activity relationship modelling of the carcinogenic risk of nitroso compounds using regression analysis and the TOPS-MODE approach

Worldwide, legislative and governmental efforts are focusing on establishing simple screening tools for identifying those chemicals most likely to cause adverse effects without experimentally testing all chemicals of regulatory concern. This is because even the most basic biological testing of compounds of concern, apart from requiring a huge number of test animals, would be neither resource nor time effective. Thus, alternative approaches such as the one proposed here, quantitative structure-activity relationship (QSAR) modelling, are increasingly being used for identifying the potential health hazards and subsequent regulation of new industrial chemicals. This paper follows up on our earlier work that demonstrated the use of the TOPological Substructural MOlecular DEsign (TOPS-MODE) approach to QSAR modelling for predictions of the carcinogenic potency of nitroso compounds. The data set comprises 56 nitroso compounds which have been bio-assayed in female rats and administered by the oral water route. The QSAR model was able to account for about 81% of the variance in the experimental activity and exhibited good cross-validation statistics. A reasonable interpretation of the TOPS-MODE descriptors was achieved by means of bond contributions, which in turn afforded the recognition of structural alerts (SAs) regarding carcinogenicity. A comparison of the SAs obtained from different data sets showed that experimental factors, such as the sex and the oral administration route, exert a major influence on the carcinogenicity of nitroso compounds. The present and previous QSAR models combined together provide a reliable tool for estimating the carcinogenic potency of yet untested nitroso compounds and they should allow the identification of SAs, which can be used as the basis of prediction systems for the rodent carcinogenicity of these compounds.

Quantitative structure carcinogenicity relationship for detecting structural alerts in nitroso-compounds

Prevention of environmentally induced cancers is a major health problem of which solutions depend on the rapid and accurate screening of potential chemical hazards. Lately, theoretical approaches such as the one proposed here - Quantitative Structure-Activity Relationship (QSAR) - are increasingly used for assessing the risks of environmental chemicals, since they can markedly reduce costs, avoid animal testing, and speed up policy decisions. This paper reports a QSAR study based on the Topological Substructural Molecular Design (TOPS-MODE) approach, aiming at predicting the rodent carcinogenicity of a set of nitroso-compounds selected from the Carcinogenic Potency Data Base (CPDB). The set comprises nitrosoureas (14 chemicals), N-nitrosamines (18 chemicals) C-nitroso-compounds (1 chemical), nitrosourethane (1 chemical) and nitrosoguanidine (1 chemical), which have been bioassayed in male rat using gavage as the route of administration. Here we are especially concerned in gathering the role of both parameters on the carcinogenic activity of this family of compounds. First, the regression model was derived, upon removal of one identified nitrosamine outlier, and was able to account for more than 84% of the variance in the experimental activity. Second, the TOPS-MODE approach afforded the bond contributions -- expressed as fragment contributions to the carcinogenic activity -- that can be interpreted and provide tools for better understanding the mechanisms of carcinogenesis. Finally, and most importantly, we demonstrate the potentialities of this approach towards the recognition of structural alerts for carcinogenicity predictions.

Phenolic acid derivatives with potential anticancer properties--a structure-activity relationship study. Part 1: methyl, propyl and octyl esters of caffeic and gallic acids

The antiproliferative and cytotoxic properties of polyphenolic acid derivatives, structurally related with the natural models caffeic and gallic acids, have been tested in human cervix adenocarcinoma cells (HeLa). Simultaneous structural information was obtained for these compounds through theoretical ab initio methods. This study was conducted for the following esters: methyl caffeate (MC, 1), propyl caffeate (PC, 2), octyl caffeate (OC, 3), methyl gallate (MG, 4), propyl gallate (PG, 5) and octyl gallate (OG, 6). A significant growth-inhibition effect was assessed for some of these compounds, clearly dependent on their structural characteristics. Marked structure-activity relationships (SARs)--namely the number of hydroxyl ring substituents--were found to rule the biological effect of such systems.

Synthesis and QSAR study of the anticancer activity of some novel indane carbocyclic nucleosides

A set of 14 indane carbocyclic nucleosides were synthesized and experimentally assayed for their inhibitory effects in the proliferation of murine leukemia (L1210/0) and human T-lymphocyte (Molt4/C8, CEM/0) cells. The compounds have promising inhibitory activity judging from the IC(50) values obtained for all these cellular lines. Multiple linear regression analysis was then applied to build up consistent QSAR models based on quantum mechanics-derived molecular descriptors. The derived models reproduce well the experimental data of both three cells (r(2) >/=0.90), display a good predictive power and are, above all, easily interpretable. They show that frontier-orbital energies and hydrophobicity are mainly responsible for the activity of the synthesized compounds and also, suggest similar mechanisms of action. The final QSAR-models involve only two descriptors: the lowest unoccupied molecular orbital energy and the solvent accessible-hydrophobic surface area, but describe a sound correlation between predicted and experimental activity data (r(2)=0.931, r(2)=0.936 and r(2)=0.931 for the cells L1210/0, Molt4/C8 and CEM/0, respectively).