Homology Modeling and Molecular Dynamics-Driven Search for Natural Inhibitors That Universally Target Receptor-Binding Domain of Spike Glycoprotein in SARS-CoV-2 Variants.

The rapid spread of SARS-CoV-2 required immediate actions to control the transmission of the virus and minimize its impact on humanity. An extensive mutation rate of this viral genome contributes to the virus' ability to quickly adapt to environmental changes, impacts transmissibility and antigenicity, and may facilitate immune escape. Therefore, it is of great interest for researchers working in vaccine development and drug design to consider the impact of mutations on virus-drug interactions. Here, we propose a multitarget drug discovery pipeline for identifying potential drug candidates which can efficiently inhibit the Receptor Binding Domain (RBD) of spike glycoproteins from different variants of SARS-CoV-2. Eight homology models of RBDs for selected variants were created and validated using reference crystal structures. We then investigated interactions between host receptor ACE2 and RBDs from nine variants of SARS-CoV-2. It led us to conclude that efficient multi-variant targeting drugs should be capable of blocking residues Q(R)493 and N487 in RBDs. Using methods of molecular docking, molecular mechanics, and molecular dynamics, we identified three lead compounds (hesperidin, narirutin, and neohesperidin) suitable for multitarget SARS-CoV-2 inhibition. These compounds are flavanone glycosides found in citrus fruits - an active ingredient of Traditional Chinese Medicines. The developed pipeline can be further used to (1) model mutants for which crystal structures are not yet available and (2) scan a more extensive library of compounds against other mutated viral proteins.

Using quasi-SMILES for the predictive modeling of the safety of 574 metal oxide nanoparticles measured in different experimental conditions.

The production of nanomaterials continues its rapid growth; however, newly manufactured nanomaterials' environmental and health safety are among the most significant concerns. A safety assessment is usually a lengthy and costly process, so computational studies are often used to complement experimental testing. One of the most time-efficient techniques is structure-activity relationships (SAR) modeling. In this project, we analyzed the Sustainable Nanotechnology (S2NANO) dataset that contains 574 experimental cell viability and toxicity datapoints for Al2O3, CuO, Fe2O3, Fe3O4, SiO2, TiO2, and ZnO measured in different conditions. We aimed to develop classification- and regression-based structure-activity relationship models using quasi-SMILES molecular representation. Introduced quasi-SMILES took into consideration all available information, including structural features of nanoparticles (molecular structure, core size, etc.) and related experimental parameters (cell line, dose, exposure time, assay, hydrodynamic size, surface charge, etc.). Resultant regression models demonstrated sufficient predictive power, while classification models demonstrated higher accuracy.

Zeta potentials (ζ) of metal oxide nanoparticles: A meta-analysis of experimental data and a predictive neural networks modeling.

Zeta potential is usually measured to estimate the surface charge and the stability of nanomaterials, as changes in these characteristics directly influence the biological activity of a given nanoparticle. Nowadays, theoretical methods are commonly used for a pre-screening safety assessments of nanomaterials. At the same time, the consistency of data on zeta potential measurements in the context of environmental impact is an important challenge. The inconsistency of data measurements leads to inaccuracies in predictive modeling. In this article, we report a new curated dataset of zeta potentials measured for 208 silica- and metal oxide nanoparticles in different media. We discuss the data curation framework for zeta potentials designed to assess the quality and usefulness of the literature data for further computational modeling. We also provide an analysis of specific trends for the datapoints harvested from different literature sources. In addition to that, we present for the first time a structure-property relationship model for nanoparticles (nano-SPR) that predicts values of zeta potential values measured in different environmental conditions (i.e., biological media and pH).

N-arylnaphthylamines as inhibitors of human immunodeficiency virus integrase - lens epithelium-derived growth factor interactions: theoretical studies.

Presented work reports a comprehensive theoretical study on the inhibitory nature of N-arylnaphthylamines in Human Immunodeficiency Virus Integrase (HIV IN) - Lens Epithelium-Derived Growth Factor (LEDGF/p75) complexes. Factors influencing the inhibition efficiency in AlphaScreen% assay are evaluated and explained through the structure- and ligand-based studies; including molecular docking, molecular dynamics calculations, and quantitative structure-activity relationship (QSAR) approach. It has been shown that N-arylnaphthylamines possess a wide variety of binding poses. Three QSAR models have been developed using structural descriptors and descriptors derived from docking calculations. The activity of untested N-arylnaphthylamines have been predicted using the most successful model. Proposed here technique could become a useful tool for ligand selection, accelerating the development of a new generation of anti-HIV medications. [Formula: see text] Communicated by Ramaswamy H. Sarma.

Predictive Modeling of Critical Temperatures in Superconducting Materials.

In this study, we have investigated quantitative relationships between critical temperatures of superconductive inorganic materials and the basic physicochemical attributes of these materials (also called quantitative structure-property relationships). We demonstrated that one of the most recent studies (titled "A data-driven statistical model for predicting the critical temperature of a superconductor" and published in Computational Materials Science by K. Hamidieh in 2018) reports on models that were based on the dataset that contains 27% of duplicate entries. We aimed to deliver stable models for a properly cleaned dataset using the same modeling techniques (multiple linear regression, MLR, and gradient boosting decision trees, XGBoost). The predictive ability of our best XGBoost model (R2 = 0.924, RMSE = 9.336 using 10-fold cross-validation) is comparable to the XGBoost model by the author of the initial dataset (R2 = 0.920 and RMSE = 9.5 K in ten-fold cross-validation). At the same time, our best model is based on less sophisticated parameters, which allows one to make more accurate interpretations while maintaining a generalizable model. In particular, we found that the highest relative influence is attributed to variables that represent the thermal conductivity of materials. In addition to MLR and XGBoost, we explored the potential of other machine learning techniques (NN, neural networks and RF, random forests).

Evaluating genotoxicity of metal oxide nanoparticles: Application of advanced supervised and unsupervised machine learning techniques.

Presence of missing data points in datasets is among main challenges in handling the toxicological data for nanomaterials. As the processing of missing data is an important part of data analysis, we have introduced a read-across approach that uses a combination of supervised and unsupervised machine learning techniques to fill the missing values. A series of classification models (supervised learning) was developed to predict class label, and self-organizing map approach (unsupervised learning) was used to estimate relative distances between nanoparticles and refine results obtained during supervised learning. In this study, genotoxicity of 49 silicon and metal oxide nanoparticles in Ames and Comet tests. Collected literature data did not demonstrate significant variations related to the change of size including selected bulk materials. Genotoxicity-related features of nanomaterials were represented by ionic characteristics. General tendencies found in the current study were convincingly linked to known theories of genotoxic action at nano-level. Mechanisms of primary and secondary genotoxic effects were discussed in the context of developed models.

A chemoinformatics approach for the characterization of hybrid nanomaterials: safer and efficient design perspective.

In this study, photocatalytic properties and in vitro cytotoxicity of 29 TiO2-based multi-component nanomaterials (i.e., hybrids of more than two composition types of nanoparticles) were evaluated using a combination of the experimental testing and supervised machine learning modeling. TiO2-based multi-component nanomaterials with metal clusters of silver, and their mixtures with gold, palladium, and platinum were successfully synthesized. Two activities, photocatalytic activity and cytotoxicity, were studied. A novel cheminformatic approach was developed and applied for the computational representation of the photocatalytic activity and cytotoxicity effect. In this approach, features of investigated TiO2-based hybrid nanomaterials were reflected by a series of novel additive descriptors for hybrid and hybrid nanostructures (denoted as "hybrid nanosctructure descriptors"). These descriptors are based on quantum chemical calculations and the Smoluchowski equation. The obtained experimental data and calculated hybrid-nanostructure descriptors were used to develop novel predictive Quantitative Structure-Activity Relationship computational models (called "nano-QSARmix"). The proposed modeling approach is an initial step in the understanding of the relationships between physicochemical properties of hybrid nanoparticles, their toxicity, and photochemical activity under UV-vis irradiation. Acquired knowledge supports the safe-by-design approaches relevant to the development of efficient hybrid nanomaterials with reduced hazardous effects.

Modeling of Glass Transition Temperatures for Polymeric Coating Materials: Application of QSPR Mixture-based Approach.

Cross-linked block copolymers are structurally complex, and utilization of traditional methods of molecular representation in chemoinformatics is only of limited applicability. Therefore, we introduced new techniques of structural representation for block copolymers. We developed additive and combinatorial approaches that treat a copolymer as a mixture system. In this approach, DRAGON descriptors are concentration-weighted for all chemicals in the reaction mixture. As a proof of concept, we have studied glass transition temperatures of block copolymers of hydroxyalkyl- and dihydroxyalkyl carbamate terminated poly(dimethylsiloxane) oligomers with poly(-caprolactone) and developed four quantitative structure-property relationships (QSPR) models. The correlation coefficient (R2 ) for mentioned QSPR models ranges from 0.851 to 0.911 for the training set. In addition to the newly introduced technique we found that the octanol-water partition coefficient and 3D-MoRSE unweighted descriptors were the most important descriptors for the studied property. The results of the study demonstrated that all chemicals in reaction mixture influenced the glass transition temperatures.

Multiple e-Pharmacophore modeling to identify a single molecule that could target both streptomycin and paromomycin binding sites for 30S ribosomal subunit inhibition.

The bacterial ribosome is an established target for anti-bacterial therapy since decades. Several inhibitors have already been developed targeting both defined subunits (50S and 30S) of the ribosome. Aminoglycosides and tetracyclines are two classes of antibiotics that bind to the 30S ribosomal subunit. These inhibitors can target multiple active sites on ribosome that have a complex structure. To screen putative inhibitors against 30S subunit of the ribosome, the crystal structures in complex with various known inhibitors were analyzed using pharmacophore modeling approach. Multiple active sites were considered for building energy-based three-dimensional (3D) pharmacophore models. The generated models were validated using enrichment factor on decoy data-set. Virtual screening was performed using the developed 3D pharmacophore models and molecular interaction towards the 30S ribosomal unit was analyzed using the hits obtained for each pharmacophore model. The hits that were common to both streptomycin and paromomycin binding sites were identified. Further, to predict the activity of these hits a robust 2D-QSAR model with good predictive ability was developed using 16 streptomycin analogs. Hence, the developed models were able to identify novel inhibitors that are capable of binding to multiple active sites present on 30S ribosomal subunit.

Genotoxicity induced by metal oxide nanoparticles: a weight of evidence study and effect of particle surface and electronic properties.

The genetic toxicology of nanomaterials is a crucial toxicology issue and one of the least investigated topics. Substantially, the genotoxicity of metal oxide nanomaterials' data is resulting from in vitro comet assay. Current contributions to the genotoxicity data assessed by the comet assay provide a case-by-case evaluation of different types of metal oxides. The existing inconsistency in the literature regarding the genotoxicity testing data requires intelligent assessment strategies, such as weight of evidence evaluation. Two main tasks were performed in the present study. First, the genotoxicity data from comet assay for 16 noncoated metal oxide nanomaterials with different core composition were collected. An evaluation criterion was applied to establish which of these individual lines of evidence were of sufficient quality and what weight could have been given to them in inferring genotoxic results. The collected data were surveyed on (1) minimum necessary characterization points for nanomaterials and (2) principals of correct comet assay testing for nanomaterials. Second, in this study the genotoxicity effect of metal oxide nanomaterials was investigated by quantitative nanostructure-activity relationship approach. A set of quantum-chemical descriptors was developed for all investigated metal oxide nanomaterials. A classification model based on decision tree was developed for the investigated dataset. Thus, three descriptors were identified as the most responsible factors for genotoxicity effect: heat of formation, molecular weight, and surface area of the oxide cluster based on the conductor-like screening model. Conclusively, the proposed genotoxicity assessment strategy is useful to prioritize the study of the nanomaterials for further risk assessment evaluations.

Towards the Development of Global Nano-Quantitative Structure-Property Relationship Models: Zeta Potentials of Metal Oxide Nanoparticles.

Zeta potential indirectly reflects a charge of the surface of nanoparticles in solutions and could be used to represent the stability of the colloidal solution. As processes of synthesis, testing and evaluation of new nanomaterials are expensive and time-consuming, so it would be helpful to estimate an approximate range of properties for untested nanomaterials using computational modeling. We collected the largest dataset of zeta potential measurements of bare metal oxide nanoparticles in water (87 data points). The dataset was used to develop quantitative structure-property relationship (QSPR) models. Essential features of nanoparticles were represented using a modified simplified molecular input line entry system (SMILES). SMILES strings reflected the size-dependent behavior of zeta potentials, as the considered quasi-SMILES modification included information about both chemical composition and the size of the nanoparticles. Three mathematical models were generated using the Monte Carlo method, and their statistical quality was evaluated (R² for the training set varied from 0.71 to 0.87; for the validation set, from 0.67 to 0.82; root mean square errors for both training and validation sets ranged from 11.3 to 17.2 mV). The developed models were analyzed and linked to aggregation effects in aqueous solutions.

QSPR modeling of optical rotation of amino acids using specific quantum chemical descriptors.

Many chemical phenomena occur in solution. Different solvents can change the optical activity of chiral molecules. The optical rotation angles of solutes of 75 amino acids in dimethylformamide, water and methanol were analyzed using the quantitative structure-activity relationships approach. For an accurate description of chirality, we used specific quantum chemical descriptors, which reflect the properties of a chiral center, and continuous symmetry measures. The set of specific quantum chemical descriptors for atoms located near the chiral center, such as Mulliken charges, the sum of Mulliken charges on an atom (with the hydrogen charges summed up with the adjacent non-hydrogen atoms), and nuclear magnetic resoncance tensors was applied. To represent solvent effects, we used mixture-like structural simplex descriptors and quantum chemical descriptors obtained for structures optimized for specified solvent using PBE1PBE/6-31G** level of theory with the polarizable continuum model. Multiple linear regression, M5P, and locally weighted learning techniques were used to achieve accurate predictions. The specific quantum chemical descriptors proposed here demonstrated high specificity in the majority of the developed models and established direct quantitative structure-property relationships.

How the toxicity of nanomaterials towards different species could be simultaneously evaluated: a novel multi-nano-read-across approach.

Application of predictive modeling approaches can solve the problem of missing data. Numerous studies have investigated the effects of missing values on qualitative or quantitative modeling, but only a few studies have discussed it for the case of applications in nanotechnology-related data. The present study is aimed at the development of a multi-nano-read-across modeling technique that helps in predicting the toxicity of different species such as bacteria, algae, protozoa, and mammalian cell lines. Herein, the experimental toxicity of 184 metal and silica oxide (30 unique chemical types) nanoparticles from 15 datasets is analyzed. A hybrid quantitative multi-nano-read-across approach that combines interspecies correlation analysis and self-organizing map analysis is developed. In the first step, hidden patterns of toxicity among nanoparticles are identified using a combination of methods. Subsequently, the developed model based on categorization of the toxicity of the metal oxide nanoparticle outcomes is evaluated via the combination of supervised and unsupervised machine learning techniques to determine the underlying factors responsible for the toxicity.

Evaluating the toxicity of TiO2-based nanoparticles to Chinese hamster ovary cells and Escherichia coli: a complementary experimental and computational approach.

Titania-supported palladium, gold and bimetallic nanoparticles (second-generation nanoparticles) demonstrate promising photocatalytic properties. However, due to unusual reactivity, second-generation nanoparticles can be hazardous for living organisms. Considering the ever-growing number of new types of nanoparticles that can potentially contaminate the environment, a determination of their toxicity is extremely important. The main aim of presented study was to investigate the cytotoxic effect of surface modified TiO2-based nanoparticles, to model their quantitative nanostructure-toxicity relationships and to reveal the toxicity mechanism. In this context, toxicity tests for surface-modified TiO2-based nanoparticles were performed in vitro, using Gram-negative bacteria Escherichia coli and Chinese hamster ovary (CHO-K1) cells. The obtained cytotoxicity data were analyzed by means of computational methods (quantitative structure-activity relationships, QSAR approach). Based on a combined experimental and computational approach, predictive models were developed, and relationships between cytotoxicity, size, and specific surface area (Brunauer-Emmett-Teller surface, BET) of nanoparticles were discussed.

Binding of Alkali Metal Ions with 1,3,5-Tri(phenyl)benzene and 1,3,5-Tri(naphthyl)benzene: The Effect of Phenyl and Naphthyl Ring Substitution on Cation-π Interactions Revealed by DFT Study.

The effect of substitution of phenyl and naphthyl rings to benzene was examined to elucidate the cation-π interactions involving alkali metal ions with 1,3,5-tri(phenyl)benzene (TPB) and 1,3,5-tri(naphthyl)benzene (TNB). Benzene, TPB, and four TNB isomers (with ααα, ααß, αßß, and ßßß types of fusion) and their complexes with Li+, Na+, K+, Rb+, and Cs+ were optimized using DFT approach with B3LYP and M06-2X functionals in conjunction with the def2-QZVP basis set. Higher relative stability of ß,ß,ß-TNB over α,α,α-TNB can be attributed to peri repulsion, which is defined as the nonbonding repulsive interaction between substituents in the 1- and the 8-positions on the naphthalene core. Binding energies, distances between ring centroid and the metal ions, and the distance to metal ions from the center of other six-membered rings were compared for all complexes. Our computational study reveals that the binding affinity of alkali metal cations increases significantly with the 1,3,5-trisubstitution of phenyl and naphthyl rings to benzene. The detailed computational analyses of geometries, partial charges, binding energies, and ligand organization energies reveal the possibility of favorable C-H···M+ interactions when a α-naphthyl group exists in complexes of TNB structures. Like benzene-alkali metal ion complexes, the binding affinity of metal ions follows the order: Li+ > Na+ > K+ > Rb+ > Cs+ for any considered 1,3,5-trisubstituted benzene systems. In case of TNB, we found that the strength of interactions increases as the fusion point changes from α to ß position of naphthalene.

Modeling of Interactions between the Zebrafish Hatching Enzyme ZHE1 and A Series of Metal Oxide Nanoparticles: Nano-QSAR and Causal Analysis of Inactivation Mechanisms.

The quantitative relationships between the activity of zebrafish ZHE1 enzyme and a series of experimental and physicochemical features of 24 metal oxide nanoparticles were revealed. Vital characteristics of the nanoparticles' structure were reflected using both experimental and theoretical descriptors. The developed quantitative structure-activity relationship model for nanoparticles (nano-QSAR) was capable of predicting the enzyme inactivation based on four descriptors: the hydrodynamic radius, mass density, the Wigner-Seitz radius, and the covalent index. The nano-QSAR model was calculated using the non-linear regression tree M5P algorithm. The developed model is characterized by high robustness R²bagging = 0.90 and external predictivity Q²EXT = 0.93. This model is in agreement with modern theories of aquatic toxicity. Dissolution and size-dependent characteristics are among the key driving forces for enzyme inactivation. It was proven that ZnO, CuO, Cr2O3, and NiO nanoparticles demonstrated strong inhibitory effects because of their solubility. The proposed approach could be used as a non-experimental alternative to animal testing. Additionally, methods of causal discovery were applied to shed light on the mechanisms and modes of action.

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.

Chlorophenol sorption on multi-walled carbon nanotubes: DFT modeling and structure-property relationship analysis.

The presence of chlorophenols in drinking water can be hazardous to human health. Understanding the mechanisms of adsorption under specific experimental conditions would be beneficial when developing methods to remove toxic substances from drinking water during water treatment in order to limit human exposure to these contaminants. In this study, we investigated the sorption of chlorophenols on multi-walled carbon nanotubes using a density functional theory (DFT) approach. This was applied to study selected interactions between six solvents, five types of nanotubes, and six chlorophenols. Experimental data were used to construct structure-adsorption relationship (SAR) models that describe the recovery process. Specific interactions between solvents and chlorophenols were taken into account in the calculations by using novel specific mixture descriptors.

A DFT-based QSAR study on inhibition of human dihydrofolate reductase.

Diaminopyrimidine derivatives are frequently used as inhibitors of human dihydrofolate reductase, for example in treatment of patients whose immune system are affected by human immunodeficiency virus. Forty-seven dicyclic and tricyclic potential inhibitors of human dihydrofolate reductase were analyzed using the quantitative structure-activity analysis supported by DFT-based and DRAGON-based descriptors. The developed model yielded an RMSE deviation of 1.1 a correlation coefficient of 0.81. The prediction set was characterized by R2=0.60 and RMSE=3.59. Factors responsible for inhibition process were identified and discussed. The resulting model was validated via cross validation and Y-scrambling procedure. From the best model, we found several mass-related descriptors and Sanderson electronegativity-related descriptors that have the best correlations with the investigated inhibitory concentration. These descriptors reflect results from QSAR studies based on characteristics of human dihydrofolate reductase inhibitors.